Genome-wide analysis reveals regulatory role of G4 DNA in gene transcription

Investigation of mRNA quadruplex formation in Escherichia coli

The protocol presented here allows for the investigation of the formation of unusual nucleic acid structures in the 5'-untranslated region (UTR) of bacteria by correlating gene expression levels to the in vitro stability of the respective structure. In particular, we describe the introduction of G-quadruplex forming sequences close to the ribosome-binding site (RBS) on the mRNA of a reporter gene and the subsequent read-out of the expression levels. Insertion of a stable secondary structure results in the cloaking of RBS and eventually reduced gene expression levels. The structures and stability of the introduced sequences are further characterized by circular dichroism (CD) spectroscopy and thermal melting experiments. The extent of inhibition is then correlated to the stability of the respective quadruplex structure, allowing judgement of whether factors other than thermodynamic stability affect the formation of a given quadruplex sequence in vivo. Measuring gene expression levels takes 2 d including cloning; CD experiments take 5 hours per experiment.

Molecular models for intrastrand DNA G-quadruplexes

BACKGROUND

Independent surveys of human gene promoter regions have demonstrated an overrepresentation of G(3)X(n1)G3X(n2)G(3)X(n3)G(3) motifs which are known to be capable of forming intrastrand quadruple helix structures. In spite of the widely recognized importance of G-quadruplex structures in gene regulation and growing interest around this unusual DNA structure, there are at present only few such structures available in the Nucleic Acid Database. In the present work we generate by molecular modeling feasible G-quadruplex structures which may be useful for interpretation of experimental data.

RESULTS

We have used all quadruplex DNA structures deposited in the Nucleic Acid Database in order to select a list of fragments entailing a strand of three adjacent G's paired with another strand of three adjacent G's separated by a loop of one to four residues. These fragments were further clustered and representative fragments were finally selected. Further fragments were generated by assemblying the two strands of each fragment with loops from different fragments whenever the anchor G's were superimposable. The fragments were used to assemble G quadruplex based on a superimposability criterion.

CONCLUSION

Molecular models have been generated for a large number of G(3)X(n1)G(3)X(n2)G3X(n3)G(3) sequences. For a given sequence not all topologies are possible with the available repertoire of fragments due to steric hindrance and low superimposability. Since all molecular models are generated by fragments coming from observed quadruplex structures, molecular models are in principle reliable and may be used for interpretation of experimental data. Some examples of applications are given.

Genome-wide colonization of gene regulatory elements by G4 DNA motifs

G-quadruplex (or G4 DNA), a stable four-stranded structure found in guanine-rich regions, is implicated in the transcriptional regulation of genes involved in growth and development. Previous studies on the role of G4 DNA in gene regulation mostly focused on genomic regions proximal to transcription start sites (TSSs). To gain a more comprehensive understanding of the regulatory role of G4 DNA, we examined the landscape of potential G4 DNA (PG4Ms) motifs in the human genome and found that G4 motifs, not restricted to those found in the TSS-proximal regions, are bias toward gene-associated regions. Significantly, analyses of G4 motifs in seven types of well-known gene regulatory elements revealed a constitutive enrichment pattern and the clusters of G4 motifs tend to be colocalized with regulatory elements. Considering our analysis from a genome evolutionary perspective, we found evidence that the occurrence and accumulation of certain progenitors and canonical G4 DNA motifs within regulatory regions were progressively favored by natural selection. Our results suggest that G4 DNA motifs are ‘colonized’ in regulatory regions, supporting a likely genome-wide role of G4 DNA in gene regulation. We hypothesize that G4 DNA is a regulatory apparatus situated in regulatory elements, acting as a molecular switch that can modulate the role of the host functional regions, by transition in DNA structure.

Genome-wide analysis of a G-quadruplex-specific single-chain antibody that regulates gene expression

G-quadruplex nucleic acids have been proposed to play a role in a number of fundamental biological processes that include transcription and translation. We have developed a single-chain antibody that is selective for G-quadruplex DNA over double-stranded DNA, and here show that when it is expressed in human cells, it significantly affects the expression of a wide variety of genes, in a manner that correlates with the presence of predicted G-quadruplexes. We observe cases where gene expression is increased or decreased, and that there are apparent interactions with G-quadruplex motifs at the beginning and end of the genes, and on either strand. The outcomes of this genome-wide study demonstrate that G-quadruplex recognition by the antibody has physiological consequences, and provides insights into some of the complexity associated with G-quadruplex-based regulation.

Predictable suppression of gene expression by 5'-UTR-based RNA quadruplexes

Four-stranded DNA and RNA quadruplexes or G4 motifs are non-B DNA conformations that are presumed to form in vivo, although only few explicit evidence has been reported. Using bioinformatics the presence of putative DNA G-quadruplexes within critical promoter regions has been demonstrated and a regulatory role in transcription has been suspected. However, in genomic DNA the presence of the complementary strand interferes with the potential to form a quadruplex motif. Contrarily RNA G4 motifs have no such limitation and consequently strong interference with gene expression is suspected. Nevertheless, experimental evidence is scarce. Here we show a well-defined structure–function relationship of synthetic quadruplex sequences in 5′-UTRs in multiple mammalian cell-lines. We establish a universal ‘translational suppressor’ effect of these motifs on gene expression at the translational level and show for the first time that specific features such as loop-length and the number of ‘GGG’-repeats further determine the suppressive impact. Moreover, a consistent and predictable repression of gene expression is observed for naturally occurring RNA G4 motifs, augmenting the functional relevance of these unusual nucleic acid structures.

Comparative genomics of chondrichthyan Hoxa clusters

BACKGROUND

The chondrichthyan or cartilaginous fish (chimeras, sharks, skates and rays) occupy an important phylogenetic position as the sister group to all other jawed vertebrates and as an early lineage to diverge from the vertebrate lineage following two whole genome duplication events in vertebrate evolution. There have been few comparative genomic analyses incorporating data from chondrichthyan fish and none comparing genomic information from within the group. We have sequenced the complete Hoxa cluster of the Little Skate (Leucoraja erinacea) and compared to the published Hoxa cluster of the Horn Shark (Heterodontus francisci) and to available data from the Elephant Shark (Callorhinchus milii) genome project.

RESULTS

A BAC clone containing the full Little Skate Hoxa cluster was fully sequenced and assembled. Analyses of coding sequences and conserved non-coding elements reveal a strikingly high level of conservation across the cartilaginous fish, with twenty ultraconserved elements (100%,100 bp) found between Skate and Horn Shark, compared to three between human and marsupials. We have also identified novel potential non-coding RNAs in the Skate BAC clone, some of which are conserved to other species.

CONCLUSION

We find that the Little Skate Hoxa cluster is remarkably similar to the previously published Horn Shark Hoxa cluster with respect to sequence identity, gene size and intergenic distance despite over 180 million years of separation between the two lineages. We suggest that the genomes of cartilaginous fish are more highly conserved than those of tetrapods or teleost fish and so are more likely to have retained ancestral non-coding elements. While useful for isolating homologous DNA, this complicates bioinformatic approaches to identify chondrichthyan-specific non-coding DNA elements.

The disruptive positions in human G-quadruplex motifs are less polymorphic and more conserved than their neutral counterparts

Specific guanine-rich sequence motifs in the human genome have considerable potential to form four-stranded structures known as G-quadruplexes or G4 DNA. The enrichment of these motifs in key chromosomal regions has suggested a functional role for the G-quadruplex structure in genomic regulation. In this work, we have examined the spectrum of nucleotide substitutions in G4 motifs, and related this spectrum to G4 prevalence. Data collected from the large repository of human SNPs indicates that the core feature of G-quadruplex motifs, 5′-GGG-3′, exhibits specific mutational patterns that preserve the potential for G4 formation. In particular, we find a genome-wide pattern in which sites that disrupt the guanine triplets are more conserved and less polymorphic than their neutral counterparts. This also holds when considering non-CpG sites only. However, the low level of polymorphisms in guanine tracts is not only confined to G4 motifs. A complete mapping of DNA three-mers at guanine polymorphisms indicated that short guanine tracts are the most under-represented sequence context at polymorphic sites. Furthermore, we provide evidence for a strand bias upstream of human genes. Here, a significantly lower rate of G4-disruptive SNPs on the non-template strand supports a higher relative influence of G4 formation on this strand during transcription.

The structures of quadruplex nucleic acids and their drug complexes

Quadruplex nucleic acids are four-stranded structures formed from short tracts of G-rich sequence associating together. Their formation from eukaryotic telomeric DNA sequences is well established, at least in vitro, and has more recently been the focus of attention as novel anticancer targets since their formation inhibits the telomerase complex from maintaining telomere length in cancer cells. Structural studies have revealed a diversity of topologies for telomeric quadruplexes, which are sensitive to the nature of the cations present, to the flanking sequences, and probably also to concentration. The small number of structures determined for ligand-quadruplex complexes have to date shown a marked preference for the quadruplex component to have all strands in a parallel folding arrangement, which maximises the planar surface available for ligand binding. Informatics has revealed that there are a large number of quadruplex sequences present in human and other genomes, both eukaryotic and prokaryotic, with over-representation in upstream promoter regions. Structures for a few of these are now available, which emphasise the role played by loop sequences in determining fold.

The yeast Pif1 helicase prevents genomic instability caused by G-quadruplex-forming CEB1 sequences in vivo

In budding yeast, the Pif1 DNA helicase is involved in the maintenance of both nuclear and mitochondrial genomes, but its role in these processes is still poorly understood. Here, we provide evidence for a new Pif1 function by demonstrating that its absence promotes genetic instability of alleles of the G-rich human minisatellite CEB1 inserted in the Saccharomyces cerevisiae genome, but not of other tandem repeats. Inactivation of other DNA helicases, including Sgs1, had no effect on CEB1 stability. In vitro, we show that CEB1 repeats formed stable G-quadruplex (G4) secondary structures and the Pif1 protein unwinds these structures more efficiently than regular B-DNA. Finally, synthetic CEB1 arrays in which we mutated the potential G4-forming sequences were no longer destabilized in pif1Δ cells. Hence, we conclude that CEB1 instability in pif1Δ cells depends on the potential to form G-quadruplex structures, suggesting that Pif1 could play a role in the metabolism of G4-forming sequences.

Changes in the primary DNA sequence are a major source of pathologies and cancers. The hereditary information also resides in secondary DNA structures, a layer of genetic information that remains poorly understood. Biophysical and structural studies have long established that, in vitro, the DNA molecule can adopt diverse structures different from the canonical Watson-Crick conformations. However, for a long time their existence in vivo has been regarded with a certain skepticism and their functional role elusive. One example is the G-quadruplex structure, which involves G-quartets that form between four DNA strands. Here, using in vitro and in vivo assays in the yeast S. cerevisiae, we reveal the unexpected role of the Pif1 helicase in maintaining the stability of the human CEB1 G-rich tandem repeat array. By site-directed mutagenesis, we show that the genomic instability of CEB1 repeats in absence of Pif1 and is directly dependent on the ability of CEB1 to form G-quadruplex structures. We show that Pif1 is very efficient in vitro in processing G-quadruplex structures formed by CEB1. We propose that Pif1 maintains CEB1 repeats by its ability to resolve G-quadruplex structures, thus providing circumstantial evidence of their formation in vivo.

Identification and prevention of genotyping errors caused by G-quadruplex- and i-motif-like sequences

BACKGROUND

Reliable PCR amplification of DNA fragments is the prerequisite for most genetic assays. We investigated the impact of G-quadruplex- or i-motif-like sequences on the reliability of PCR-based genetic analyses.

METHODS

We found the sequence context of a common intronic polymorphism in the MEN1 gene (multiple endocrine neoplasia I) to be the cause of systematic genotyping errors by inducing preferential amplification of one allelic variant [allele dropout (ADO)]. Bioinformatic analyses and pyrosequencing-based allele quantification enabled the identification of the underlying DNA structures.

RESULTS

We showed that G-quadruplex- or i-motif-like sequences can reproducibly cause ADO. In these cases, amplification efficiency strongly depends on the PCR enzyme and buffer conditions, the magnesium concentration in particular. In a randomly chosen subset of candidate single-nucleotide polymorphisms (SNPs) defined by properties deduced from 2 originally identified ADO cases, we confirmed preferential PCR amplification in up to 50% of the SNPs. We subsequently identified G-quadruplex and i-motifs harboring a SNP that alters the typical motif as the cause of this phenomenon, and a genomewide search based on the respective motifs predicted 0.5% of all SNPs listed by dbSNP and Online Mendelian Inheritance in Man to be potentially affected.

CONCLUSIONS

Undetected, the described phenomenon produces systematic errors in genetic analyses that may lead to misdiagnoses in clinical settings. PCR products should be checked for G-quadruplex and i-motifs to avoid the formation of ADO-causing secondary structures. Truly affected assays can then be identified by a simple experimental procedure, which simultaneously provides the solution to the problem.

Non-B DNA conformations as determinants of mutagenesis and human disease

Repetitive DNA motifs may fold into non-B DNA structures, including cruciforms/hairpins, triplexes, slipped conformations, quadruplexes, and left-handed Z-DNA, thereby representing chromosomal targets for DNA repair, recombination, and aberrant DNA synthesis leading to repeat expansion or genomic rearrangements associated with neurodegenerative and genomic disorders. Hairpins and quadruplexes also determined the relative abundances of simple sequence repeats (SSR) in vertebrate genomes, whereas strong base stacking has permitted the expansion of purine.pyrimidine-rich SSR during evolutionary time. SSR are enriched in regulatory and cancer-related gene classes, where they have been actively recruited to participate in both gene and protein functions. SSR polymorphic alleles in the population are associated with cancer susceptibility, including within genes that appear to share regulatory circuits involving reactive oxygen species.

Elevated polyamines induce c-MYC overexpression by perturbing quadruplex-WC duplex equilibrium

The biological role of quadruplexes and polyamines has been independently associated with cancer. However, quadruplex-polyamine mediated transcriptional regulation remain unaddressed. Herein, using c-MYC quadruplex model, we have attempted to understand quadruplex–polyamine interaction and its role in transcriptional regulation. We initially employed biophysical approach involving CD, UV and FRET to understand the role of polyamines (spermidine and spermine) on conformation, stability, molecular recognition of quadruplex and to investigate the effect of polyamines on quadruplex–Watson Crick duplex transition. Our study demonstrates that polyamines affect the c-MYC quadruplex conformation, perturb its recognition properties and delays duplex formation. The relative free energy difference (ΔΔG°) between the duplex and quadruplex structure indicate that polyamines stabilize and favor c-MYC quadruplex over duplex. Further, we investigated the influence of polyamine mediated perturbation of this equilibrium on c-MYC expression. Our results suggest that polyamines induce structural transition of c-MYC quadruplex to a transcriptionally active motif with distinctive molecular recognition property, which drives c-MYC expression. These findings may allow exploiting quadruplex–polyamines interaction for developing antiproliferative strategies to combat aberrant gene expression.

Evidence of genome-wide G4 DNA-mediated gene expression in human cancer cells

Guanine-rich DNA of a particular sequence adopts four-stranded structural forms known as G-quadruplex or G4 DNA. Though in vitro formation of G4 DNA is known for several years, in vivo presence of G4 DNA was only recently noted in eukaryote telomeres. Recent bioinformatics analyses showing prevalence of G4 DNA within promoters of human and related species seems to implicate G4 DNA in a genome-wide cis-regulatory role. Herein we demonstrate that G4 DNA may present regulatory sites on a genome-wide scale by showing widespread effect on gene expression in response to the established intracellular G4 DNA-binding ligands. This is particularly relevant to genes that harbor conserved potential G4 DNA (PG4 DNA) forming sequence across human, mouse and rat promoters of orthologous genes. Genes with conserved PG4 DNA in promoters show co-regulated expression in 79 human and 61 mouse normal tissues (z-score > 3.5; P < 0.0001). Conservation of G4 DNA across related species also emphasizes the biological importance of G4 DNA and its role in transcriptional regulation of genes; shedding light on a relatively novel mechanism of regulation of gene expression in eukaryotes.

A novel chair-type G-quadruplex formed by a Bombyx mori telomeric sequence

Recently, the human telomeric d[TAGGG(TTAGGG)₃] sequence has been shown to form in K⁺ solution an intramolecular (3+1) G-quadruplex structure, whose G-tetrad core contains three strands oriented in one direction and the fourth in the opposite direction. Here we present a study on the structure of the Bombyx mori telomeric d[TAGG(TTAGG)₃] sequence, which differs from the human counterpart only by one G deletion in each repeat. We found that this sequence adopted multiple G-quadruplex structures in K⁺ solution. We have favored a major G-quadruplex form by a judicious U-for-T substitution in the sequence and determined the folding topology of this form. We showed by NMR that this was a new chair-type intramolecular G-quadruplex which involved a two-layer antiparallel G-tetrad core and three edgewise loops. Our result highlights the effect of G-tract length on the folding topology of G-quadruplexes, but also poses the question of whether a similar chair-type G-quadruplex fold exists in the human telomeric sequences.

Metastases suppressor NM23-H2 interaction with G-quadruplex DNA within c-MYC promoter nuclease hypersensitive element induces c-MYC expression

Regulatory influence of the G-quadruplex or G4 motif present within the nuclease hypersensitive element (NHE) in the promoter of c-MYC has been noted. On the other hand, association of NM23-H2 to the NHE leads to c-MYC activation. Therefore, NM23-H2 interaction with the G4 motif within the c-MYC NHE presents an interesting mechanistic possibility. Herein, using luciferase reporter assay and chromatin immunoprecipitation we show NM23-H2 mediated c-MYC activation involves NM23-H2-G4 motif binding within the c-MYC NHE. G4 motif complex formation with recombinant NM23-H2 was independently confirmed using fluorescence energy transfer, which also indicated that the G4 motif was resolved to an unfolded state within the protein-bound complex. Taken together, this supports transcriptional role of NM23-H2 via a G4 motif.

Recognition and discrimination of DNA quadruplexes by acridine-peptide conjugates

We have explored a series of trisubstituted acridine-peptide conjugates for their ability to recognize and discriminate between DNA quadruplexes derived from the human telomere, and the c-kit and N-ras proto-oncogenes. Quadruplex affinity was measured as the peptide sequences were varied, together with their substitution position on the acridine, and the identity of the C-terminus (acid or amide). Surface plasmon resonance measurements revealed that all compounds bound to the human telomeric quadruplex with sub-micromolar affinity. Docking calculations from molecular modelling studies were used to model the effects of substituent orientation and peptide sequence. Modelling and experiment were in agreement that placement of the peptide over the face of the acridine is detrimental to binding affinity. The highest degrees of selectivity were observed towards the N-ras quadruplex by compounds capable of forming simultaneous contacts with their acridine and peptide moieties. The ligands that bound best displayed quadruplex affinities in the 1-5 nM range and at least 10-fold discrimination between the quadruplexes studied.

Genome-wide computational and expression analyses reveal G-quadruplex DNA motifs as conserved cis-regulatory elements in human and related species

Using a combination of in silico and experimental approaches, we present evidence that the G-quadruplex (G4) motif (an alternative higher-order DNA conformation) has regulatory potential. Genome-wide analyses of 99980 human, chimpanzee, mouse, and rat promoters showed enrichment of sequence with potential to adopt G4 (potential G4 or PG4) motifs near transcription start sites (TSS; P < 0.0001), supporting earlier findings. Interestingly, we found >700 orthologously related promoters in human, mouse, and rat conserve PG4 motif(s). The corresponding genes have enriched (z score > 4.0) tissue-specific expression in 75 of 79 human tissues and are significantly overrepresented in signaling and regulation of cell-cycle (P < 10(-05)). This is supported by results from whole genome expression experiments in human HeLa S3 cells following treatment with TMPyP4 [5,10,15,20-tetra(N-methyl-4-pyridyl) porphine chloride], which is known to bind the G4 motif inside cells. Our results implicate G4-motif mediated regulation as a more general mode of transcription control than currently appreciated.

A thermodynamic overview of naturally occurring intramolecular DNA quadruplexes

Loop length and its composition are important for the structural and functional versatility of quadruplexes. To date studies on the loops have mainly concerned model sequences compared with naturally occurring quadruplex sequences which have diverse loop lengths and compositions. Herein, we have characterized 36 quadruplex-forming sequences from the promoter regions of various proto-oncogenes using CD, UV and native gel electrophoresis. We examined folding topologies and determined the thermodynamic profile for quadruplexes varying in total loop length (5–18 bases) and composition. We found that naturally occurring quadruplexes have variable thermodynamic stabilities (ΔG₃₇) ranging from −1.7 to −15.6 kcal/mol. Overall, our results suggest that both loop length and its composition affect quadruplex structure and thermodynamics, thus making it difficult to draw generalized correlations between loop length and thermodynamic stability. Additionally, we compared the thermodynamic stability of quadruplexes and their respective duplexes to understand quadruplex–duplex competition. Our findings invoke a discussion on whether biological function is associated with quadruplexes with lower thermodynamic stability which undergo facile formation and disruption, or by quadruplexes with high thermodynamic stability.

Mix and measure fluorescence screening for selective quadruplex binders

The human genome contains thousands of regions, including that of the telomere, that have the potential to form quadruplex structures. Many of these regions are potential targets for therapeutic intervention. There are many different folding patterns for quadruplex DNAs and the loops exhibit much more variation than do the quartets. The successful targeting of a particular quadruplex structure requires distinguishing that structure from all of the other quadruplex structures that may be present. A mix and measure fluorescent screening method has been developed, that utilizes multiple reporter molecules that bind to different features of quadruplex DNA. The reporter molecules are used in combination with DNAs that have a variety of quadruplex structures. The screening is based on observing the increase or decrease in the fluorescence of the reporter molecules. The selectivity of a set of test molecules has been determined by this approach.

Transcription induces strand-specific mutations at the 5' end of human genes

A regional analysis of nucleotide substitution rates along human genes and their flanking regions allows us to quantify the effect of mutational mechanisms associated with transcription in germ line cells. Our analysis reveals three distinct patterns of substitution rates. First, a sharp decline in the deamination rate of methylated CpG dinucleotides, which is observed in the vicinity of the 5' end of genes. Second, a strand asymmetry in complementary substitution rates, which extends from the 5' end to 1 kbp downstream from the 3' end, associated with transcription-coupled repair. Finally, a localized strand asymmetry, an excess of C-->T over G-->A substitution in the nontemplate strand confined to the first 1-2 kbp downstream of the 5' end of genes. We hypothesize that higher exposure of the nontemplate strand near the 5' end of genes leads to a higher cytosine deamination rate. Up to now, only the somatic hypermutation (SHM) pathway has been known to mediate localized and strand-specific mutagenic processes associated with transcription in mammalia. The mutational patterns in SHM are induced by cytosine deaminase, which just targets single-stranded DNA. This DNA conformation is induced by R-loops, which preferentially occur at the 5' ends of genes. We predict that R-loops are extensively formed in the beginning of transcribed regions in germ line cells.