Skip to main content

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

Abstract

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.

Background

Nucleic acids containing guanine tracts can form quadruplex structures via non-Watson-Crick base pairing [1, 2]. These structures, known as G-quadruplexes, are composed of stacked G-tetrads - square co-planar arrays of four guanine bases each (Figure 1). Each tetrad in the G-quadruplex is stabilized by cyclic Hoogsteen hydrogen bonding between the four guanines [3–6]. G-quadruplexes can be formed by unimolecular interactions via repeated folding of a single nucleic acid molecule or through multimolecular interactions between two or four strands. Our work focuses on unimolecular G-quadruplex formation which has been the target of most genomic studies and is more likely to be encountered in physiological conditions [6, 7].

Figure 1
figure 1

G-quadruplex structure. Left: an intramolecular G-quadruplex formed by a G3AG3CUG3CG3 RNA motif. Right: a G-tetrad structure.

Formation of G-quadruplexes in DNA and RNA is associated with regulation of important biological functions [8]. There have been several reports of G-quadruplex involvement in a variety of human diseases [9–11]. Consequently, G-quadruplexes are being considered as targets for therapeutic purposes [9, 12, 13]. DNA G-quadruplexes can influence transcription, telomere biology, genomic rearrangements and instability, DNA recombination, DNA repair, DNA replication, and epigenetic mechanisms [14–18]. RNA G-quadruplexes readily form in vivo[19] and are considered more stable than DNA G-quadruplexes [20, 21].

G-quadruplexes have been reported in untranslated as well as translated regions of mRNAs [22, 23] and have been shown to function as cis-regulatory elements of translation and alternative splicing [24–29]. For example, G-quadruplex structures located in the 5′-UTR (untranslated region) of human fibroblast growth factor 2 (FGF2) [30] and human vascular endothelial growth factor (VEGF) [25] mRNAs act as internal ribosomal entry sites (IRES) for cap-independent translation initiation. Deletion of a conserved G-quadruplex motif in the 5′-UTR of MECP2 gene has been associated with reduced protein product and autism, suggesting that a G-quadruplex plays a role in normal protein synthesis [11]. Formation of G-quadruplexes in this region can also inhibit translation of neuroblastoma Ras (NRAS) oncogene [24] and ADAM10 - responsible for anti-amyloidogenic processing of the amyloid precursor protein (APP) [31]. G-rich sequences capable of forming G-quadruplexes in the 3′-UTR have been shown to enhance cleavage and polyadenylation of mammalian pre-mRNAs [32–35], influence RNA turnover [29], and to play a role in subcellular mRNA sorting [36].

G-quadruplexes can be detected in the lab using UV spectroscopy and fluorescent tags [37, 38]. Structure-specific antibodies can be used to quantitatively visualize DNA G-quadruplexes in human cells [39, 40]. However, for large-scale analysis, researchers turn to computational approaches for identifying regions likely to form G-quadruplexes. Bioinformatics studies have been conducted to map potential G-quadruplex motifs in genomic regions [41–46], splice regions [47], and untranslated regions [48] of mammals, and genomes of yeast [49] and prokaryotes [50]. Most computational algorithms rely on identifying sequence motifs capable of forming G-quadruplex structures. These algorithms primarily analyze G-tract lengths and loop lengths in the nucleotide sequences.

In the cell, many factors - such as the presence of a cation within or between co-planar guanines - may contribute to G-quadruplex formation. Furthermore, G-quadruplex formation can be transient - they may form and dissolve under physiological conditions intermittently. These complexities can make computational mappings subject to false-positives, where motifs occurring by chance will be identified along with motifs representing real G-quadruplexes with significant biological impact.

One way to reduce the likelihood of false-positives when predicting G-quadruplex formation is by examining the motif's conservation across species. Nucleotides which are integral to G-quadruplex formation are known to be conserved and less likely to be polymorphic [51]. This increased level of conservation has been specifically observed in Saccharomyces[52]. Capra and colleagues measured this conservation in two ways: at the motif level and at the nucleotide level. Motif-level conservation was defined as a binary result based on whether an overlapping motif was found in two aligned sequences. Nucleotide conservation was calculated by comparing the number of conserved bases within a motif against conserved nucleotides within the neighboring 100 nt.

In this paper we present quadruplex forming G-rich sequences (QGRS)-Conserve, a strategy for identifying highly conserved G-quadruplex motifs in homologous nucleotide sequences. Our approach differs from Capra and co-workers in that it evaluates not only location conservation but also the conservation of structural features of the G-quadruplex motif - such as loop lengths, number of tetrads, and the total length of the structure. We also present strategies for dealing with specific challenges relating to overlapping G-quadruplex motifs and the impact they have on conservation analysis. Our technique allows for the filtering of motifs based on qualitative conservation and promotes accurate wide-scale analysis of G-quadruplexes within exomes, transcriptomes, and genomes. Since conservation suggests biological relevance, our tools can help researchers focus on motifs likely to be of high scientific interest. We present the summative results from a database of G-quadruplexes created using the QGRS-Conserve.

Methods

QGRS-Conserve is a computational method for calculating motif conservation across exomes that supports filtering to provide researchers with more precise methods of studying distributions patterns. The work presented in this paper focuses on G-quadruplexes found in mRNA; however, the algorithms and principles could be similarly applied when working with DNA.

Identifying G-quadruplexes

QGRS-Conserve identifies QGRS (predicted G-quadruplexes) in nucleotide sequences based on algorithms described previously [53, 54]. Putative G-quadruplexes are identified using the motif G x Ny1G x Ny2G x Ny3G x . The motif consists of four guanine tracts of equal size interspersed by three loops. In this expression, x represents the number of guanine tetrads, and y1–y3 represent the loops (gaps). The size of each G-tract corresponds to the number of stacked G-tetrads forming the quadruplex structure (Figure 1).

This motif expression is consistent with those in other published methods, although some methods place differing restrictions on the number of tetrads and/or loop lengths applied [55, 56].

The numeric quantification of predicted stability is useful in facilitating filtering and searching across large datasets. To quantify stability, we calculate ‘G-scores’ based on a previously published strategy [53, 54] for each identified QGRS. While this stability is not used in our conservation measurements, it can be used to filter which QGRS that are ultimately used for conservation calculations.

Measuring conservation

Finding a QGRS motif - especially one that yields a high G-score - is of interest. However, finding QGRS that appear to be under evolutionary constraint - and hence likely to play important biological roles - is our larger goal. Measuring a QGRS's conservation across homologous genes is a way of further measuring, indirectly, the likelihood that the QGRS does indeed form a G-quadruplex. The combination of evolutionary conservation and high G-scores can be a strong indicator that the nucleotide region in question requires more direct study.

QGRS-Conserve measures both location and structural conservation of QGRS, using four numerical components: (1) their location (overlap) within aligned sequences, (2) tetrad similarity, (3) loop length similarity, and (4) total QGRS length similarity. For each pairing of QGRS between two homologs, a weighted score between 0 and 1 is computed - where 1 indicates high conservation. We will refer to this score as the pair s ‘conservation score.’

Location conservation

Aligned homologous nucleotide sequences are expected to contain significant conserved stretches, representing evolutionarily conserved function. It is therefore reasonable to assume that QGRS motifs in the same or similar location across aligned nucleotide sequence homologs have undergone similar evolutionary pressure and can be assigned with homologous status. To evaluate how well the two QGRS motifs align, a semi-global alignment between homologs is performed; the principal nucleotide sequence is typically from Homo sapiens and the comparison sequence is from another species. A semi-global alignment is preferred over the full global alignment because of the variation between the lengths of untranslated regions across orthologous mRNAs.

Sixty-five percent of our conservation score is dedicated to the overall location similarity between the two QGRS. Once a set of homologous sequences has been aligned, each pair of QGRS is assigned an overlap score. The overlap score, like the overall conservation score, ranges from 0.0 to 1.0. The premise is that pairs that do not occupy the same relative region should have low scores approaching 0, while pairs highly proximate (while not necessarily precisely aligned) receive scores nearing 1.0.

To compute the overlap score for a given QGRS pair, a start and ending nucleotide position is determined for the QGRS in the principal and comparison sequences. Recognizing that up/downstream sequence variations can easily cause QGRS pairs to be misaligned by several nucleotides, we apply padding to each start and end point of the QGRS of 50% the total length of the smaller QGRS. The percentage overlap between the modified endpoints is computed as shown in Figure 2. Note that padding prevents the outright elimination of proximate but non-overlapping QGRS.

Figure 2
figure 2

Overlap percentage computation. The extension of sequence length using 50% of the total sequence length as padding inflates overlap scores for proximate sequences.

Any pairing with an overlap score greater than 0 is transformed through a step function to eliminate the penalty of slightly misaligned motifs. Overlap percentages of 85% or greater is given a score of 1.0, and all other scores are assigned an interpolation between 0 and 1.

Motif structural conservation

An additional 20% of the conservation score is based on similarity in the number of tetrads (length of G-tracts). The number of tetrads can determine stability, as well as regulatory characteristics of G-quadruplexes [57]. Each QGRS has a number of tetrads ranging from 2 to N, and the percentage similarity among the two is determined by dividing the smaller over the larger. This initial ratio is transformed by a step function to assign a numeric score, where any ratio at or below 50% is assigned a 0 and between 50% and 100% is mapped between 0.5 and 1.0. For example, two QGRS of three and four tetrads would have a similarity score of 0.75, where two with two and four tetrads would be assigned a score of 0.

Since the number of tetrads tends to be quite small (most QGRS have tetrad counts less than 5), we have also accounted for fractional differences. Rather than considering QGRS pairs with two and three tetrads to have a difference of a full tetrad, we consider the individual tracts within the motif with shorter G-tracts to examine how close it is to becoming a full three-tetrad QGRS. For example, the QGRS GGGACGTTCGGATTGGGTTACCAGGG would be considered (only for the purposes of the tetrad similarity computation) to have a tetrad length of 2.75 - as three of the four tracts have an adjacent extra guanine, meaning that a single mutation in the other tract would have yielded a three-tetrad QGRS.

This method creates a more favorable tetrad similarity score, and therefore a higher conservation score, and more accurately measures the similarity of G-quadruplex structures putatively formed by QGRS motifs in question.

Similarity among loop lengths is also factored into the conservation score calculation, accounting for 10% of the overall score. The loop length similarity is taken as the average percentage similarity of each of the three loops. Each loop's similarity is simply the smaller loop length divided by the larger, and a similar step function is applied to the score to transform any percentage less than 50 to 0.

The last 5% of the conservation score is the overall length score. If there is 60% overall length similarity or lower, a length score of 0 is assigned. Larger percentages are interpolated to scale along a range from 0 to 1. The overall length score only accounts for 5% because the parameter of the overall length is already somewhat incorporated in the tetrad and loop length similarity scores but is still important to consider and record as it can be used to filter data as needed.

Three of the four numerical components are based on structural determinants of G-quadruplex formation. For these components, a relatively higher score (20%) was chosen for tetrad similarity among structural components because the number of tetrads can determine stability as well as the regulatory characteristics of G-quadruplexes [57]. A fair evaluation of conservation between a QGRS pair should, therefore, involve a close comparison of the number of tetrads.

Since the location of G-quadruplex motifs in the mRNAs is of high biological relevance, a 65% weighting was assigned to location/overlap numerical component of conservation score.

We have developed our quantitative scoring methods and component weightings in an effort to generally reflect the current state of awareness in the G-quadruplex field. It is an active area of research and to the best of our knowledge, ours is one of the first attempts to quantify G-quadruplex motif conservation.

It should be noted that other factors observed in literature have not been incorporated into our conservation score, including loop composition similarity and ion binding preference similarity of G-quadruplexes [58–60]. These characteristics of G-quadruplex structures are still being actively researched. Until a clear picture emerges, we felt that it would not be possible to incorporate these factors into the conservation score at present time.

Implementation

QGRS-Conserve is implemented as a series of three stages (Figure 3). Once the full sequence data for the principal and comparison genes have been obtained, the first stage is to identify all QGRS within each sequence independently.

Figure 3
figure 3

Algorithm stages for SNCA mRNA. Overview of the algorithm stages for QGRS identification (independently) on each mRNA, performing a semi-global alignment on the homologs, and finally evaluating the QGRS pair for conservation.

Q p = set of all QGRS found within principal mRNA
(1)
Q c = set of all QGRS found within comparison mRNA
(2)

The second stage is a semi-global alignment of the two sequences. A semi-global alignment is performed without consideration of the location of the QGRS found in the first stage and produces an optimal alignment over the entire sequence - avoiding penalizing the gaps at the ends of the strands. The semi-global alignment typically injects gaps into each sequence representing differences with the aligned sequence. Care is taken within the program to preserve a record of the location of each QGRS found in the first stage both before and after these gaps are introduced. Our program uses the Needleman-Wunsch Semi-Global alignment web service provided by EMBOSS, which can be accessed at (http://www.ebi.ac.uk/Tools/webservices/services/psa/emboss_needle_soap).

The final stage is the calculation of the conservation score. Here, each QGRS in Qp is paired with each QGRS in Qc. Many of these pairings are extraneous, they represent QGRS that maybe in very different locations within their respective genes and are not necessarily homologous - the conservation score computation will automatically dispose of low scoring pairs, as their overlap score will be 0. For each pair, the composite conservation score consisting of scores for overlap (using nucleotide positions after the gaps from alignment are added), loop length similarity, tetrad similarity, and overall length similarity is computed. Loop length and overall length similarities are based on nucleotide counts without the gaps resulting from alignment.

Efficient large-scale analysis

The QGRS-Conserve algorithm has been made available through a web application previously called QGRS Predictor [61], the functionality of which has been published elsewhere [62]. The interactive web interface allows a user to examine a single pair of mRNA sequences at a time, which is useful when researchers already know which set of homologs they are interested in studying. The broader goal of our research is to study the distribution patterns of conserved QGRS across entire exomes/transcriptomes for better understanding the biological roles of G-quadruplex motifs in post-transcriptional regulation of gene expression.

We have created a database of conserved QGRS using QGRS-Conserve in conjunction with listings of homologous mRNA pairings available through NCBI's Homologene database [63]. At the time of writing, our master database contains conservation computations of the entire human exome compared against homologs within the exomes of Mus musculus, Pan troglodytes, Canis lupus familiaris, Danio rerio, Caenorhabditis elegans, and Kluyveromyces lactis. Large-scale analysis presents many challenges beyond the interactive use-case, and we now describe additional solutions we have developed specifically to allow efficient storage and analyses of QGRS across entire exomes.

Overlapping QGRS

QGRS within a sequence are said to overlap when their positions (start to end) in the nucleotide sequence overlap (Figure 4). Overlapping QGRS may differ in total length, loop size, or loop sequence; any of them have the potential to form intramolecular quadruplexes in the cell, many of which may be stable. Several strategies for dealing with overlapping QGRS can be found in the literature. One approach is to focus on regions where G-quadruplexes may be formed, rather than specific motifs. These regions, which contain overlapping and proximate G-quadruplexes, are called regions with ‘G4 potential,’ or G4P [43]. If the presentation of specific motifs is desired (such as when evaluating structural conservation), overlapping motifs may be preserved [41] or they may be filtered based on the selection of the motif with the highest number of G-tracts, the higher predicted stability [54], or other characteristics. Typically, this filtering is acceptable, as the tools merely highlight regions in a sequence where G-quadruplexes are likely to be found; they do not claim the specific motif presented will be the one and only G-quadruplex formed among the overlapping motifs.

Figure 4
figure 4

Overlapping QGRS yielding different G-scores. A small sample of hundreds of overlapping QGRS found in the ‘GGATGAGCCTGTGGGTGGAGGGCAGTGGAGCGGGCTGGG’ sequence of human GREB1 gene.

When presenting the analysis between two mRNA sequences in the interactive interface, all pairings of QGRS across the species are considered, and all homologous pairs meeting the user's specified filter criteria (i.e., minimum conservation score) are presented. This listing includes the pairing of overlapping QGRS. We chose this approach because when performing analysis to measure conservation across species, the removal of overlapping QGRS cannot be done without impact on results. When highlighting conservation, we not only consider the characteristics of an individual QGRS but also how similar it is with QGRS in the comparison species sequence. Due to this complication, it is not optimal to remove overlaps until conservation computations across all QGRS pairs can be computed.

Overlapping QGRS and conservation calculations

Unlike when performing an interactive analysis, in large-scale analysis, the consideration of all overlapping pairings of QGRS is impractical and ineffective. To illustrate the challenge of dealing with overlapping QGRS when computing conservation scores, consider a situation where multiple overlapping QGRS are found in similar areas of two aligned mRNA sequences [principal and comparison (homologous) sequences] as shown in Figure 5. Each QGRS found in the principal sequence have their own characteristics (such as number of G-tracts and loop lengths), yielding a different G-score (integer to the left of the nucleotide sequence).

Figure 5
figure 5

Choosing QGRS considered for conservation calculations. The number next to each QGRS motif represents its calculated G-score. This figure demonstrates how it is possible that QGRS with highest G-score do not necessarily yield the highest conservation score.

We have adopted the term ‘QGRS family’ to represent a set of overlapping QGRS within a sequence. When choosing which QGRS within a family to compare with QGRS in the homologous sequence, we could choose the one with the highest predicted stability (G-score); however, optimizing on predicted stability first may not yield the highest conservation score present among QGRS families. It is quite possible that comparing QGRS instances with slightly suboptimal predicted stability (i.e., relatively lower G-scores) may yield better similarity. As an example, in Figure 5, selecting the QGRS instance within each family with the highest G-score might yield a conservation score of 0.81, where the optimal score across the two families might be far larger (0.92 in this case).

Comparing only QGRS with the highest G-Score can potentially underestimate conservation. However, computing and storing all combinations of overlapping motif present serious problems when conducting any large-scale analysis, not only for performance reason but because the presence of overlapping pairs can skew aggregate statistics concerning the numbers of QGRS found in particular regions of the sequences.

Filtering overlapping motifs

To address overlapping motifs, we used a compromise that allows us to report the highest conservation measurement across families of overlapping motifs while maintaining only a single specific representative motif from each family. The solution has the characteristic of identifying representative QGRS within specific mRNA independently from homologs it is compared against - allowing for the cross-referencing of the same QGRS in one gene (e.g., from human) against many homologs from other species.

Step 1: locating QGRS instances in each mRNA

For a given pair of mRNA (principal and comparison), the first stage of processing simply identifies all QGRS - including overlapped motifs, giving us two sets of QGRS.

Q p = set of all QGRS found within principal mRNA , including overlaps
(3)
Q c = set of all QGRS found within comparison mRNA , including overlaps
(4)

Step 2: group QGRS into overlapping families

For Qp and Qc, the QGRS are grouped into a collection of QGRS families, where each QGRS within a given family differ in nucleotide position by less than five bases with all other QGRS in the family. Note that it is possible (while not necessarily common) for QGRS to fall into multiple families. This provides us with two sets of families:

F p = subsets ofQ p where each QGRS in a given family overlap
(5)
F c = subsets ofQ c where each QGRS in a given family overlap
(6)

Step 3: persisting QGRS instances to the database

A QGRS record is created for each family in Fp and Fc. The QGRS record's identifier in the database will correspond to the family number (ordered from 5′ to 3′ by starting nucleotide position). The structural characteristics (G-tracts, loop lengths, etc.) and G-score will be saved corresponding to the QGRS within the family with the highest G-score. Note while only one QGRS per family is stored in the database as a QGRS record, the set of all QGRS found in the family are not discarded immediately - they are utilized in step 5.

Step 4: pair families in Fp and Fcfor conservation computation

Each QGRS family in Fp is paired with each family in Fc, yielding another set X f , where each element in Xf is a pair of QGRS families (one from principal and one from comparison mRNA).

X f = cross product of F p and F c
(7)

For each pair of QGRS families in Xf, a conservation record will be created in the database. A conservation record will contain the identifier of the QGRS from the principal and comparison mRNA. In this case, the conservation record will reference the representatives selected in step 3.

Step 5: conservation computations

In order to compute the conservation between two QGRS families, a new set of pairings, Xq, is computed. Xq represents each QGRS instance from the principal family paired with each from the comparison family.

X q = cross product of QGRS instances in family from principal with family from comparison
(8)

The highest conservation score (calculated using the procedure described in previous sections) in Xq for each pair of families in Xf is then persisted to the database.

Our filtering design has several implications. Firstly, all possible QGRS pairs are evaluated across a set of homologous mRNA, not only the ones with the highest G-score among them. This prevents the analysis from missing highly conserved motifs that may warrant special investigation. Secondly, individual QGRS recorded for a particular mRNA correspond to the QGRS within each family with the highest G-score. This means the set of actual QGRS stored for a particular mRNA is independent from conservation scores across homologs. Since the same QGRS represent a single family of overlapping QGRS regardless of which homolog it is being compared with, it is possible to cross-reference QGRS in a single mRNA with potentially dozens of homologs. This feature allows for queries such as ‘find all QGRS in the human SNCA gene which have highly homologous representations in three or more other homologs from other species.’

Results and discussion

The algorithms and methods behind QGRS-Conserve allow for efficient identification of conserved QGRS across large sets of nucleotide sequences and include strategies to handle overlapping QGRS without sacrificing the fidelity of the conservation scoring. These strategies are critical to make large-scale analysis efficient and tractable.

We have used QGRS-Conserve to create a database of QGRS records and corresponding conservation scores focused on human mRNA compared against mRNA homologs from seven species. At the time of writing, over 68,000 mRNA homolog pairs have been analyzed, yielding over 433,000 computationally mapped QGRS instances in the human exome alone and over 1 million across all seven species. This dataset includes information about 343,153 human QGRS (in 16,684 human mRNAs) conserved in at least one ortholog. Conservation (conservation score > 0.95) of human QGRS was calculated against homologs in seven species: M musculus, P troglodytes, C. lupus familiaris, D rerio, C. elegans, and K. lactis.

We have collected information for the categorization of QGRS by tetrads, loop lengths, G-score, conservation, absolute location (both nucleotide positions, and regional - 5′-UTR, CDS, 3′-UTR), relative location (e.g. distance from polyA signals/sites), mRNA similarity (alignment scores), and gene ontology terms.

We have created a web interface, open to the public, which allows researchers to browse the database. This site can be found at http://qdb.ramapo.edu. The data is fully searchable and allows the user to obtain listings of QGRS filtered by mRNA sequence, associated gene ontology, stability (G-score), number of tetrads, and conservation score. As we continue to incorporate more species and search functionality, such features will be added to the web site.

Table 1 shows an overview of our database and highlights the ability to filter out potential false-positives at the motif level using conservation. When considering conservation between human and other species such as dog and mouse, conservation scores filtered at 0.95 or above eliminate over 50% of the mapped human QGRS. When coupled with region considerations (for example, only considering 5′ or 3′-UTR), this computation allows researchers to locate QGRS of high potential interest quite quickly. The following sections discuss the categorizations and trends in Table 1 in more detail.

Table 1 An overview of human QGRS conservation across seven orthologous exomes

Categorizing conservation

Table 1 shows the percentage of QGRS found in H. sapiens mRNAs where homologous QGRS were found in another organism. Human QGRS were found to be conserved in untranslated as well as translated regions of orthologous mRNAs. Predicted homologous QGRS with two tetrads were expectedly more prevalent than the corresponding motifs with three tetrads, which in turn were identified with higher frequency than QGRS with four or more tetrads. The rate of conservation is highest within the CDS region, across all species examined, in line with the fact that the CDS tends to be more highly conserved as compared to the 5′-UTR and 3′-UTR regions. The untranslated mRNA regions are known to be more variable. Conserved regions in the UTRs are likely to represent regulatory motifs. Therefore, while fewer in number, the identified homologous QGRS in these otherwise less conserved regions are of particular interest for future exploration.

Conservation among phylogenetic groups

We observed that the rate of QGRS conservation generally increases as the phylogenetic distance from H. sapiens decreases. Our analysis noted significant conservation of human QGRS among mammalian species (Table 1), including conservation in the 5′- and 3′-UTR regions. In contrast, almost all of the observed QGRS conservation in non-mammalian organisms could be attributed to the CDS. Non-homeothermic organisms had significantly lower level of conservation for three-tetrad QGRS among the only two species (D. rerio and Drosophila melangoaster) in which they were detected. In general, there is a distinct lack of conserved human G-quadruplex motifs in the untranslated exomes of the lower organisms. Based on overall sequence variation between such species and H. sapiens, these results are expected.

Conservation of established G-quadruplexes

Conservation is, of course, not a prerequisite for QGRS to actually form a G-quadruplex in a living cell - it is simply a strong clue. To evaluate the potential for our algorithm to provide false negatives (where low conservation scores are assigned to actual G-quadruplexes, such that they may be filtered out during a search), we investigated how our procedure worked with known G-quadruplexes reported in the literature.

We targeted ten previously published human G-quadruplexes reported for their post-transcriptional regulatory roles in different genes. These genes are involved in a wide variety of important biological processes and diseases (Table 2). G-quadruplexes mapped to 5′- or 3′-UTRs in these genes have been reported for their involvement in regulating 3′-end RNA processing [64, 65], translation [11, 24, 31, 66–69], or subcellular localization [36]. Table 2 shows human QGRS mapped to orthologous M. musculus genes along with their corresponding QGRS-conservation scores. While similar mappings can be made with other species, we chose M. musculus due to its reasonable phylogenetic distance from humans (i.e., data from chimpanzee has extremely high nucleotide conservation overall, while very little conservation of any kind can be found in distant organisms such as fruit flies). Of these ten instances, six were found to have highly homologous QGRS in mouse orthologous genes - meaning using only computational means our tool accurately highlights QGRS already known to form G-quadruplexes. Three of the G-quadruplexes in the human genes could not immediately be matched to QGRS in mouse because the latter were found only in the untranscribed region immediately upstream of the mouse transcription start site (TSS). 5′-UTRs of the currently available mouse NRAS, ADAM10, and YY1 mRNA sequences are shorter than their corresponding established human orthologs (Table 2). To test our method for these G-quadruplexes, we manually entered the upstream data taken from corresponding genomic sequence sources. Using our method with this adjusted sequence data, we found highly homologous QGRS. The homology scores marked with double asterisk in the table are each derived from calculations using sequence data mapped in the upstream genomic sequence prepended to the mRNA sequence of the mouse homolog.

Table 2 Mapping QGRS orthologs of previously published human G-quadruplexes

Among the ten selected G-quadruplexes, the QGRS found in the KISS1 gene offer an example of where our method does not locate a highly conserved (conservation score ≥0.95) homologous QGRS. The G-quadruplex in the 3′-UTR of KISS1 was matched with a QGRS yielding a relatively lower conservation score of 0.706 (Table 2). Close examination of the KISS1 QGRS pair shows that the nearest mouse QGRS differs significantly by number of tetrads (2 vs 4), loop lengths (7, 3, 16 vs 8, 4, 3 nt) and overall length (42 vs 23 nt) when compared to the published human G-quadruplex [67]. The alignment between KISS1 human and mouse mRNA sequences is relatively low (49% overall similarity), particularly in the 3′-UTR; and it is difficult to determine if this particular human G-quadruplex is conserved. In this case, a known human G-quadruplex would likely be filtered out of search results if a user of our system was examining QGRS with only high conservation scores (ex., 0.95+). This is simply a characteristic of our approach - while conservation strongly suggests biological relevance, it is not necessarily a prerequisite for a QGRS to form a G-quadruplex in living cells - and users of our tools must recognize this important point.

Conclusions

QGRS-Conserve can identify G-quadruplex-forming sequence motifs conserved among homologous nucleic acids. Identifying evolutionarily conserved motifs helps validate computations and provides evidence for their biological relevance. Ours is one of the first strategies for quantitatively evaluating conservation of G-quadruplex motifs (QGRS) using both the location as well as structural characteristics of motif pairs. Coupled with filtering approaches related to predicted stability [53, 54], conservation scoring is very helpful in limiting false positives when identifying QGRS, especially for analyzing large datasets.

We have applied QGRS-Conserve to identify multiple homologous motifs of several previously published G-quadruplexes involved in regulating gene expression and human disease; however, we have only begun the analysis of exomic data across multiple species.

Our analysis has focused on a variety of information for categorization of QGRS including their structural characteristics, relative location, and mRNA characteristics. We have identified many homologous G-quadruplexes in the regulatory 5′- and 3′-untranslated regions and discovered significant differences in their distribution patterns among phylogenetic divisions. Our data is searchable through an interactive web application at http://qdb.ramapo.edu.

Looking at aggregate distribution patterns can be helpful for gaining insights into overall conservation of QGRS and their conserved biological roles. However, one of the most important benefits of QGRS-Conserve is the ability to quickly and efficiently identify conserved QGRS in regions where there is relatively little overall sequence conservation. QGRS regions that exhibit higher conservation than their general location in an mRNA offer exciting areas of study, as their conservation suggests regulatory roles. It will be important, as this research continues, to apply new computational methods to compare the QGRS conservation with the surrounding areas of the nucleic acid sequence to isolate regions with different evolutionary constraints and correlate it to their biological relevance.

General observations can be made when looking at entire transcriptomes/exomes; however, when nucleic acids are grouped by more specific characteristics (such as gene function), these analyses may yield more actionable results and help reveal information concerning the biological roles of the G-quadruplex in the genomes. Furthermore, studying the quantitative distribution of G-quadruplex motifs relative to functional hot spots on nucleic acids - such as polyA sites, translation initiation and termination sites, and splice sites - within a variety of mRNA clusters, can provide valuable insights into both the role of the G-quadruplexes and their relationships with the mRNAs themselves. The development of computational tools such as those presented allows us to build datasets that can be filtered and categorized in a variety of dimensions and allows these critical research questions to be better addressed.

Abbreviations

ADAM10:

A Disintegrin and Metallopeptidase Domain 10

APP:

amyloid precursor protein

CDS:

coding sequence

DNA:

deoxyribose nucleic acid

EMBOSS:

European Molecular Biology Open Software Suite

FGF2:

fibroblast growth factor 2

IGF2:

insulin-like growth factor 2

KISS1:

Kisspeptin1

IRES:

Internal Ribosomal Entry Site

LRP5:

low-density lipoprotein receptor-related protein 5

MeCP2:

methyl CpG binding protein-2

mRNA:

messenger RNA

NCBI:

National Center for Biotechnology Information

NRAS:

neuroblastoma Ras

QGRS:

quadruplex forming G-rich sequence(s)

PIM1:

Pim1 proto-oncogene

PSD95:

postsynaptic density protein 95

RNA:

ribonucleic Acid

SNCA:

synuclein, alpha (non-A4 component of amyloid precursor)

TERF2:

telomeric repeat binding factor 2

TSS:

transcription start site

UTR:

untranslated region

VEGF:

vascular endothelial growth factor

YY1:

Ying Yang 1.

References

  1. Gellert M, Lipsett MN, Davies DR: Helix formation by guanylic acid. Proc Natl Acad Sci U S A. 1962, 48: 2013-2018. 10.1073/pnas.48.12.2013.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Simonsson T: G-quadruplex DNA structures–variations on a theme. Biol Chem. 2001, 382: 621-628.

    Article  CAS  PubMed  Google Scholar 

  3. Davis JT: G-quartets 40 years later: from 5′-GMP to molecular biology and supramolecular chemistry. Angew Chem Int Ed Engl. 2004, 43: 668-698. 10.1002/anie.200300589.

    Article  CAS  PubMed  Google Scholar 

  4. Burge S, Parkinson GN, Hazel P, Todd AK, Neidle S: Quadruplex DNA: sequence, topology and structure. Nucleic Acids Res. 2006, 34: 5402-5415. 10.1093/nar/gkl655.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Ji X, Sun H, Zhou H, Xiang J, Tang Y, Zhao C: Research progress of RNA quadruplex. Nucleic Acid Ther. 2011, 21: 185-200. 10.1089/nat.2010.0272.

    Article  CAS  PubMed  Google Scholar 

  6. Halder K, Hartig JS: RNA quadruplexes. Met Ions Life Sci. 2011, 9: 125-139.

    Article  CAS  PubMed  Google Scholar 

  7. Halder K, Chowdhury S: Kinetic resolution of bimolecular hybridization versus intramolecular folding in nucleic acids by surface plasmon resonance: application to G-quadruplex/duplex competition in human c-myc promoter. Nucleic Acids Res. 2005, 33: 4466-4474. 10.1093/nar/gki750.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Maizels N, Gray LT: The G4 genome. PLoS Genet. 2013, 9: e1003468-10.1371/journal.pgen.1003468.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Balasubramanian S, Neidle S: G-quadruplex nucleic acids as therapeutic targets. Curr Opin Chem Biol. 2009, 13: 345-353. 10.1016/j.cbpa.2009.04.637.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Wu Y, Brosh RM: G-quadruplex nucleic acids and human disease. Febs J. 2010, 277: 3470-3488. 10.1111/j.1742-4658.2010.07760.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Bagga JS, D’Antonio LA: Role of conserved cis-regulatory elements in the post-transcriptional regulation of the human MECP2 gene involved in autism. Hum Genomics. 2013, 7: 19-10.1186/1479-7364-7-19.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Collie GW, Parkinson GN: The application of DNA and RNA G-quadruplexes to therapeutic medicines. Chem Soc Rev. 2011, 40: 5867-5892. 10.1039/c1cs15067g.

    Article  CAS  PubMed  Google Scholar 

  13. Faudale M, Cogoi S, Xodo LE: Photoactivated cationic alkyl-substituted porphyrin binding to g4-RNA in the 5'-UTR of KRAS oncogene represses translation. Chem Commun (Camb). 2012, 48: 874-876. 10.1039/c1cc15850c.

    Article  CAS  Google Scholar 

  14. Kumar P, Yadav VK, Baral A, Kumar P, Saha D, Chowdhury S: Zinc-finger transcription factors are associated with guanine quadruplex motifs in human, chimpanzee, mouse and rat promoters genome-wide. Nucleic Acids Res. 2011, 39: 8005-8016. 10.1093/nar/gkr536.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Zaug AJ, Podell ER, Cech TR: Human POT1 disrupts telomeric G-quadruplexes allowing telomerase extension in vitro. Proc Natl Acad Sci U S A. 2005, 102: 10864-10869. 10.1073/pnas.0504744102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Baral A, Kumar P, Pathak R, Chowdhury S: Emerging trends in G-quadruplex biology - role in epigenetic and evolutionary events. Mol Biosyst. 2013, 9 (7): 1568-1575. 10.1039/c3mb25492e.

    Article  CAS  PubMed  Google Scholar 

  17. Bena F, Gimelli S, Migliavacca E, Brun-Druc N, Buiting K, Antonarakis SE, Sharp AJ: A recurrent 14q32.2 microdeletion mediated by expanded TGG repeats. Hum Mol Genet. 2010, 19: 1967-1973. 10.1093/hmg/ddq075.

    Article  CAS  PubMed  Google Scholar 

  18. Lopes J, Piazza A, Bermejo R, Kriegsman B, Colosio A, Teulade-Fichou MP, Foiani M, Nicolas A: G-quadruplex-induced instability during leading-strand replication. Embo J. 2011, 30: 4033-4046. 10.1038/emboj.2011.316.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Wieland M, Hartig JS: RNA quadruplex-based modulation of gene expression. Chem Biol. 2007, 14: 757-763. 10.1016/j.chembiol.2007.06.005.

    Article  CAS  PubMed  Google Scholar 

  20. Joachimi A, Benz A, Hartig JS: A comparison of DNA and RNA quadruplex structures and stabilities. Bioorg Med Chem. 2009, 17: 6811-6815. 10.1016/j.bmc.2009.08.043.

    Article  CAS  PubMed  Google Scholar 

  21. Mergny JL, De Cian A, Ghelab A, Sacca B, Lacroix L: Kinetics of tetramolecular quadruplexes. Nucleic Acids Res. 2005, 33: 81-94. 10.1093/nar/gki148.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Abdelmohsen K, Tominaga K, Lee EK, Srikantan S, Kang MJ, Kim MM, Selimyan R, Martindale JL, Yang X, Carrier F, Zhan M, Becker KG, Gorospe M: Enhanced translation by Nucleolin via G-rich elements in coding and non-coding regions of target mRNAs. Nucleic Acids Res. 2011, 39: 8513-8530. 10.1093/nar/gkr488.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Endoh T, Kawasaki Y, Sugimoto N: Stability of RNA quadruplex in open reading frame determines proteolysis of human estrogen receptor alpha. Nucleic Acids Res. 2013, doi:10.1093/nar/gkt286

    Google Scholar 

  24. Kumari S, Bugaut A, Huppert JL, Balasubramanian S: An RNA G-quadruplex in the 5' UTR of the NRAS proto-oncogene modulates translation. Nat Chem Biol. 2007, 3: 218-221. 10.1038/nchembio864.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Morris MJ, Negishi Y, Pazsint C, Schonhoft JD, Basu S: An RNA G-quadruplex is essential for cap-independent translation initiation in human VEGF IRES. J Am Chem Soc. 2010, 132: 17831-17839. 10.1021/ja106287x.

    Article  CAS  PubMed  Google Scholar 

  26. Didiot MC, Tian Z, Schaeffer C, Subramanian M, Mandel JL, Moine H: The G-quartet containing FMRP binding site in FMR1 mRNA is a potent exonic splicing enhancer. Nucleic Acids Res. 2008, 36: 4902-4912. 10.1093/nar/gkn472.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Gomez D, Lemarteleur T, Lacroix L, Mailliet P, Mergny JL, Riou JF: Telomerase downregulation induced by the G-quadruplex ligand 12459 in A549 cells is mediated by hTERT RNA alternative splicing. Nucleic Acids Res. 2004, 32: 371-379. 10.1093/nar/gkh181.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Marcel V, Tran PL, Sagne C, Martel-Planche G, Vaslin L, Teulade-Fichou MP, Hall J, Mergny JL, Hainaut P, Van Dyck E: G-quadruplex structures in TP53 intron 3: role in alternative splicing and in production of p53 mRNA isoforms. Carcinogenesis. 2010, 32: 271-278.

    Article  PubMed  Google Scholar 

  29. Millevoi S, Moine H, Vagner S: G-quadruplexes in RNA biology. Wiley Interdiscip Rev RNA. 2012, 3: 495-507. 10.1002/wrna.1113.

    Article  CAS  PubMed  Google Scholar 

  30. Bonnal S, Boutonnet C, Prado-Lourenco L, Vagner S: IRESdb: the Internal Ribosome Entry Site database. Nucleic Acids Res. 2003, 31: 427-428. 10.1093/nar/gkg003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Lammich S, Kamp F, Wagner J, Nuscher B, Zilow S, Ludwig AK, Willem M, Haass C: Translational repression of the disintegrin and metalloprotease ADAM10 by a stable G-quadruplex secondary structure in its 5'-untranslated region. J Biol Chem. 2011, 286: 45063-45072. 10.1074/jbc.M111.296921.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Bagga PS, Ford LP, Chen F, Wilusz J: The G-rich auxiliary downstream element has distinct sequence and position requirements and mediates efficient 3' end pre-mRNA processing through a trans-acting factor. Nucleic Acids Res. 1995, 23: 1625-1631. 10.1093/nar/23.9.1625.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Bagga PS, Arhin GK, Wilusz J: DSEF-1 is a member of the hnRNP H family of RNA-binding proteins and stimulates pre-mRNA cleavage and polyadenylation in vitro. Nucleic Acids Res. 1998, 26: 5343-5350. 10.1093/nar/26.23.5343.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Newnham CM, Hall-Pogar T, Liang S, Wu J, Tian B, Hu J, Lutz CS: Alternative polyadenylation of MeCP2: influence of cis-acting elements and trans-acting factors. RNA Biol. 2010, 7: 361-372. 10.4161/rna.7.3.11564.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Decorsière A, Cayrel A, Vagner S, Millevoi S: Essential role for the interaction between hnRNP H/F and a G quadruplex in maintaining p53 pre-mRNA 3'-end processing and function during DNA damage. Genes Dev. 2011, 25: 220-225. 10.1101/gad.607011.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Subramanian M, Rage F, Tabet R, Flatter E, Mandel JL, Moine H: G-quadruplex RNA structure as a signal for neurite mRNA targeting. EMBO Rep. 2011, 12: 697-704. 10.1038/embor.2011.76.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Mergny JL, Phan AT, Lacroix L: Following G-quartet formation by UV-spectroscopy. FEBS Lett. 1998, 435: 74-78. 10.1016/S0014-5793(98)01043-6.

    Article  CAS  PubMed  Google Scholar 

  38. Ying L, Green JJ, Li H, Klenerman D, Balasubramanian S: Studies on the structure and dynamics of the human telomeric G quadruplex by single-molecule fluorescence resonance energy transfer. Proc Natl Acad Sci U S A. 2003, 100: 14629-14634. 10.1073/pnas.2433350100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Biffi G, Tannahill D, McCafferty J, Balasubramanian S: Quantitative visualization of DNA G-quadruplex structures in human cells. Nat Chem. 2013, 5: 182-186. 10.1038/nchem.1548.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Henderson A, Wu Y, Huang YC, Chavez EA, Platt J, Johnson FB, Brosh RM, Sen D, Lansdorp PM: Detection of G-quadruplex DNA in mammalian cells. Nucleic Acids Res:. 2013, 42: 860-869.

    Article  Google Scholar 

  41. Huppert JL, Balasubramanian S: Prevalence of quadruplexes in the human genome. Nucleic Acids Res. 2005, 33: 2908-2916. 10.1093/nar/gki609.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Todd AK, Johnston M, Neidle S: Highly prevalent putative quadruplex sequence motifs in human DNA. Nucleic Acids Res. 2005, 33: 2901-2907. 10.1093/nar/gki553.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Eddy J, Maizels N: Gene function correlates with potential for G4 DNA formation in the human genome. Nucleic Acids Res. 2006, 34: 3887-3896. 10.1093/nar/gkl529.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Verma A, Halder K, Halder R, Yadav VK, Rawal P, Thakur RK, Mohd F, Sharma A, Chowdhury S: Genome-wide computational and expression analyses reveal G-quadruplex DNA motifs as conserved cis-regulatory elements in human and related species. J Med Chem. 2008, 51: 5641-5649. 10.1021/jm800448a.

    Article  CAS  PubMed  Google Scholar 

  45. Yadav VK, Abraham JK, Mani P, Kulshrestha R, Chowdhury S: QuadBase: genome-wide database of G4 DNA–occurrence and conservation in human, chimpanzee, mouse and rat promoters and 146 microbes. Nucleic Acids Res. 2008, 36: D381-385.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Zhang R, Lin Y, Zhang CT: Greglist: a database listing potential G-quadruplex regulated genes. Nucleic Acids Res. 2008, 36: D372-376. 10.1093/nar/gkn252.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Kikin O, Zappala Z, D’Antonio L, Bagga PS: GRSDB2 and GRS_UTRdb: databases of quadruplex forming G-rich sequences in pre-mRNAs and mRNAs. Nucleic Acids Res. 2008, 36: D141-148. 10.1093/nar/gkn705.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Huppert JL, Bugaut A, Kumari S, Balasubramanian S: G-quadruplexes: the beginning and end of UTRs. Nucleic Acids Res. 2008, 36: 6260-6268. 10.1093/nar/gkn511.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Hershman SG, Chen Q, Lee JY, Kozak ML, Yue P, Wang LS, Johnson FB: Genomic distribution and functional analyses of potential G-quadruplex-forming sequences in Saccharomyces cerevisiae. Nucleic Acids Res. 2008, 36: 144-156. 10.1093/nar/gkn735.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Rawal P, Kummarasetti VB, Ravindran J, Kumar N, Halder K, Sharma R, Mukerji M, Das SK, Chowdhury S: Genome-wide prediction of G4 DNA as regulatory motifs: role in Escherichia coli global regulation. Genome Res. 2006, 16: 644-655. 10.1101/gr.4508806.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Nakken S, Rognes T, Hovig E: The disruptive positions in human G-quadruplex motifs are less polymorphic and more conserved than their neutral counterparts. Nucleic Acids Res. 2009, 37: 5749-5756. 10.1093/nar/gkp590.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Capra JA, Paeschke K, Singh M, Zakian VA: G-quadruplex DNA sequences are evolutionarily conserved and associated with distinct genomic features in Saccharomyces cerevisiae. PLoS Comput Biol. 2010, 6: e1000861-10.1371/journal.pcbi.1000861.

    Article  PubMed Central  PubMed  Google Scholar 

  53. D’Antonio L, Bagga PS: Computational methods for predicting intramolecular G-quadruplexes in nucleotide sequences. Proceedings of the IEEE Conference on Computational Systems Bioinformatics (CSB 2004). 2004, 561-562.

    Google Scholar 

  54. Kikin O, D’Antonio L, Bagga PS: QGRS Mapper: a web-based server for predicting G-quadruplexes in nucleotide sequences. Nucleic Acids Res. 2006, 34: W676-682. 10.1093/nar/gkl253.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Huppert JL: Hunting G-quadruplexes. Biochimie. 2008, 90: 1140-1148. 10.1016/j.biochi.2008.01.014.

    Article  CAS  PubMed  Google Scholar 

  56. Huppert JL: Structure, location and interactions of G-quadruplexes. Febs J. 2010, 277: 3452-3458. 10.1111/j.1742-4658.2010.07758.x.

    Article  CAS  PubMed  Google Scholar 

  57. Halder K, Wieland M, Hartig JS: Predictable suppression of gene expression by 5'-UTR-based RNA quadruplexes. Nucleic Acids Res. 2009, 37: 6811-6817. 10.1093/nar/gkp696.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Marathias VM, Bolton PH: Determinants of DNA quadruplex structural type: sequence and potassium binding. Biochemistry. 1999, 38: 4355-4364. 10.1021/bi982604+.

    Article  CAS  PubMed  Google Scholar 

  59. Olsen CM, Lee HT, Marky LA: Unfolding thermodynamics of intramolecular G-quadruplexes: base sequence contributions of the loops. J Phys Chem B. 2009, 113: 2587-2595. 10.1021/jp806853n.

    Article  CAS  PubMed  Google Scholar 

  60. Sacca B, Lacroix L, Mergny JL: The effect of chemical modifications on the thermal stability of different G-quadruplex-forming oligonucleotides. Nucleic Acids Res. 2005, 33: 1182-1192. 10.1093/nar/gki257.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Menendez C, Frees S, Bagga PS: QGRS-H Predictor. http://quadruplex.ramapo.edu/qgrs/app/start. Accessed 28 April 2014

  62. Menendez C, Frees S, Bagga PS: QGRS-H Predictor: a web server for predicting homologous quadruplex forming G-rich sequence motifs in nucleotide sequences. Nucleic Acids Res. 2012, 40: W96-W103. 10.1093/nar/gks422.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Acland AAR, Barrett T, Beck J, Benson DA, Bollin C, Bolton E, Bryant SH, Canese K, Church DM, Clark K, DiCuccio M, Dondoshansky I, Federhen S, Feolo M, Geer LY, Gorelenkov V, Hoeppner M, Johnson M, Kelly C, Khotomlianski V, Kimchi A, Kimelman M, Kitts P, Krasnov S, Kuznetsov A, Landsman D, Lipman DJ, Lu Z, Madden TL, Madej T, Maglott DR, Marchler-Bauer A, Karsch-Mizrachi I, Murphy T, Ostell J, O'Sullivan C, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2013, 41: D8-D20.

    Article  Google Scholar 

  64. Beaudoin JD, Perreault JP: Exploring mRNA 3'-UTR G-quadruplexes: evidence of roles in both alternative polyadenylation and mRNA shortening. Nucleic Acids Res. 2013, 41: 5898-5911. 10.1093/nar/gkt265.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Christiansen J, Kofod M, Nielsen FC: A guanosine quadruplex and two stable hairpins flank a major cleavage site in insulin-like growth factor II mRNA. Nucleic Acids Res. 1994, 22: 5709-5716. 10.1093/nar/22.25.5709.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Arora A, Suess B: An RNA G-quadruplex in the 3' UTR of the proto-oncogene PIM1 represses translation. RNA Biol. 2011, 8: 802-805. 10.4161/rna.8.5.16038.

    Article  CAS  PubMed  Google Scholar 

  67. Huijbregts L, Roze C, Bonafe G, Houang M, Le Bouc Y, Carel JC, Leger J, Alberti P, Roux N: DNA polymorphisms of the KiSS1 3' untranslated region interfere with the folding of a G-rich sequence into G-quadruplex. Mol Cell Endocrinol. 2012, 351: 239-248. 10.1016/j.mce.2011.12.014.

    Article  CAS  PubMed  Google Scholar 

  68. Huang W, Smaldino PJ, Zhang Q, Miller LD, Cao P, Stadelman K, Wan M, Giri B, Lei M, Nagamine Y, Vaugn JP, Akman SA, Sui G: Yin Yang 1 contains G-quadruplex structures in its promoter and 5'-UTR and its expression is modulated by G4 resolvase 1. Nucleic Acids Res. 2011, 2011: 2011-

    Google Scholar 

  69. Gomez D, Guedin A, Mergny JL, Salles B, Riou JF, Teulade-Fichou MP, Calsou P: A G-quadruplex structure within the 5'-UTR of TRF2 mRNA represses translation in human cells. Nucleic Acids Res. 2010, 38: 7187-7198. 10.1093/nar/gkq563.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported in part by a grant from the Provost Office of Ramapo College of New Jersey.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Scott Frees.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CM and MC participated in the implementation of application code and performed extensive testing and validation on the software. SF designed the application architecture, including the design and implementation of all algorithms and the methods used to manage and present the collected data. PB conceived of this study and provided expertise and leadership in the development of the algorithms and methods from a biological perspective. All authors made significant contributions to the creation of this manuscript and have read and approved its final version.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Frees, S., Menendez, C., Crum, M. et al. QGRS-Conserve: a computational method for discovering evolutionarily conserved G-quadruplex motifs. Hum Genomics 8, 8 (2014). https://doi.org/10.1186/1479-7364-8-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1479-7364-8-8

Keywords