Folding Stability of Pax9 Intronic G-Quadruplex Correlates with Relative Molar Size in Eutherians.

Eutherian dentition has been the focus of a great deal of studies in the areas of evolution, development, and genomics. The development of molar teeth is regulated by an antero-to-posterior cascade mechanism of activators and inhibitors molecules, where the relative sizes of the second (M2) and third (M3) molars are dependent of the inhibitory influence of the first molar (M1). Higher activator/inhibitor ratios will result in higher M2/M1 or M3/M1. Pax9 has been shown to play a key role in tooth development. We have previously shown that a G-quadruplex in the first intron of Pax9 can modulate the splicing efficiency. Using a sliding window approach with we analyzed the association of the folding energy (Mfe) of the Pax9 first intron with the relative molar sizes in 42 mammalian species, representing 9 orders. The Mfe of two regions located in the first intron of Pax9 were shown to be significantly associated with the M2/M1 and M3/M1 areas and mesiodistal lengths. The first region is located at the intron beginning and can fold into a stable G4 structure, whereas the second is downstream the G4 and 265 bp from intron start. Across species, the first intron of Pax9 varied in G-quadruplex structural stability. The correlations were further increased when the Mfe of the two sequences were added. Our results indicate that this region has a role in the evolution of the mammalian dental pattern by influencing the relative size of the molars.

Cruciform DNA Sequences in Gene Promoters Can Impact Transcription upon Oxidative Modification of 2'-Deoxyguanosine.

Sequences of DNA typically adopt B-form duplexes in genomes, although noncanonical structures such as G-quadruplexes, i-motifs, Z-DNA, and cruciform structures can occur. A challenge is to determine the functions of these various structures in cellular processes. We and others have hypothesized that G-rich G-quadruplex-forming sequences in human genome promoters serve to sense oxidative damage generated during oxidative stress impacting gene regulation. Herein, chemical tools and a cell-based assay were used to study the oxidation of guanine to 8-oxo-7,8-dihydroguanine (OG) in the context of a cruciform-forming sequence in a gene promoter to determine the impact on transcription. We found that OG in the nontemplate strand in the loop of a cruciform-forming sequence could induce gene expression; conversely when OG was in the same sequence on the template strand, gene expression was inhibited. A model for the transcriptional changes observed is proposed in which OG focuses the DNA repair process on the promoter to impact expression. Our cellular and biophysical studies and literature sources support the idea that removal of OG from duplex DNA by OGG1 yields an abasic site (AP) that triggers a structural shift to the cruciform fold. The AP-bearing cruciform structure is presented to APE1, which functions as a conduit between DNA repair and gene regulation. The significance is enhanced by a bioinformatic study of all human gene promoters and transcription termination sites for inverted repeats (IRs). Comparison of the two regions showed that promoters have stable and G-rich IRs at a low frequency and termination sites have many AT-rich IRs with low stability.

Potential G-Quadruplex Forming Sequences and N6-Methyladenosine Colocalize at Human Pre-mRNA Intron Splice Sites.

Maturation of mRNA in humans involves modifying the 5' and 3' ends, splicing introns, and installing epitranscriptomic modifications that are essential for mRNA biogenesis. With respect to epitranscriptomic modifications, they are usually installed in specific consensus motifs, although not all sequences are modified suggesting a secondary structural component to site selection. Using bioinformatic analysis of published data, we identify in human mature-mRNA that potential RNA G-quadruplex (rG4) sequences colocalize with the epitranscriptomic modifications N6-methyladenosine (m6A), pseudouridine (Ψ), and inosine (I). Using the only available pre-mRNA data sets from the literature, we demonstrate colocalization of potential rG4s and m6A was greatest overall and occurred in introns near 5' and 3' splice sites. The loop lengths and sequence context of the m6A-bearing potential rG4s exhibited short loops most commonly comprised of single A nucleotides. This observation is consistent with a literature report of intronic m6A found in SAG (S = C or G) consensus motifs that are also recognized by splicing factors. The localization of m6A and potential rG4s in pre-mRNA at intron splice junctions suggests that these features could function together in alternative splicing. A similar analysis for potential rG4s around sites of Ψ installation or A-to-I editing in mRNA also found a colocalization; however, the frequency was less than that observed with m6A. These bioinformatic analyses guide a discussion of future experiments to understand how noncanonical rG4 structures may collaborate with epitranscriptomic modifications in the human cellular context to impact cellular phenotype.

DNA G-quadruplex stability, position and chromatin accessibility are associated with CpG island methylation.

CpG islands (CGI) are genomic regions associated with gene promoters and involved in gene expression regulation. Despite their high CpG content and unlike bulk genomic DNA methylation pattern, these regions are usually hypomethylated. So far, the mechanisms controlling the CGI methylation patterning remain unclear. G-quadruplex (G4) structures can inhibit DNA methyltransferases 1 enzymatic activity, leading to CGI hypomethylation. Our aim was to analyse the association of G4 forming sequences (G4FS) and CGI methylation as well as to determine the intrinsic and extrinsic characteristics of G4FS that may modulate this phenomenon. Using methylation data from human embryonic stem cells (hESCs) and three hESC-derived populations, we showed that hypomethylated CpGs located inside CGI (CGI/CpG) tend to be associated with highly stable G4FS (Minimum free energy ≤ -30 kcal·mol-1 ). The association of highly stable G4FS and hypomethylation tend to be stronger when these structures are located at shorter distances (~ < 150 bp) from GCI/CpGs, when G4FS and CpGs are located within open chromatin and G4FS are inside CGI. Moreover, this association is not strongly influenced by the CpG content of CGI. Conversely, highly methylated CGI/CpG tend to be associated with low stability G4FS. Although CpGs inside CGIs without a G4FS tend to be more methylated, high stability G4FS within CGI neighbourhood were associated with decreased methylation. In summary, our data indicate that G4FS may act as protective cis elements against CGI methylation, and this effect seems to be influenced by the G4FS folding potential, its presence within CGI, CpG distance from G4FS and chromatin accessibility.