G-quadruplexes (G4) are noncanonical secondary genome structures. Aberrant formation of G4s can impair genome integrity. Investigation of the relationship between G4s and somatic structural variants (SV) in cancers could provide a better understanding of the role of G4 formation in cancer development and progression. In this study, we combined bioinformatic approaches and multiomics data to investigate the connection between G4s and the somatic SVs. Somatic SV breakpoints were significantly enriched in G4 regions, regardless of SV subtypes. This enrichment was only observed in regions demonstrated to form G4s in cells (“active quadruplexes”), rather than in regions with a sequence compatible with G4 formation but without confirmed G4 formation (“potential quadruplexes”). Several genomic features affected the connection between G4s and SVs, with the enrichment being notably strengthened at the boundary of topologically associated domains. Somatic breakpoints were also preferentially associated with G4 regions with earlier replication timing and open chromatin status. In patients with cancer with homologous recombination repair defects, G4s and somatic breakpoints were substantially more strongly associated. Machine learning models were constructed that showed that G4 propensity is a potent feature for predicting the density of SV breakpoints. Altogether, these findings suggest that the G4 structures play a critical role in modulating the production of somatic SVs in cancers. Significance: G-quadruplex structure formation constitutes a critical step in the production of somatic structural variants in cancers, suggesting G-quadruplex structures as potential targets for future cancer prevention and treatment strategies.

Cis-regulatory elements (CREs) are crucial for regulating gene expression, and G-quadruplexes (G4s), as prototypal non-canonical DNA structures, may play a role in this regulation. However, the relationship between G4s and CREs, especially with non-promoter-like functional elements, requires further systematic investigation. We aimed to investigate the associations between G4s and human cCREs (candidate CREs) inferred from the Encyclopedia of DNA Elements (ENCODE) data. We found that G4s are prominently enriched in most types of cCREs, especially those with promoter-like signatures (PLS). The co-occurrence of CTCF signals with H3K4me3 or H3K27ac signals strengthens the association between cCREs and G4s. Genetic variants in G4s, particularly within their G-runs, exhibit higher regulatory potential and deleterious effects compared to cCREs. The G-runs within G4s near transcriptional start sites (TSSs) are more evolutionarily constrained compared to G-runs in cCREs, while those far from the TSS are relatively less conserved. The presence of G4s is often linked to a more favorable local chromatin environment for the activation and execution of regulatory function of cCREs, potentially attributable to the formation of G4 secondary structures. Finally, we discovered that G4-associated cCREs exhibit widespread activation in a variety of cancers. Our study suggests that G4s are integral components of human cis-regulatory elements, extending beyond their potential role in promoters. The G4 primary sequences are associated with the localization of CREs, while the G4 structures are linked to the activation of these elements. Therefore, we propose defining G4s as pivotal regulatory elements in the human genome.

Cis -regulatory elements have the ability to precisely regulate gene expression in cells, and G-quadruplexes (G4s), as non-canonical secondary structures, may potentially play a role in this regulation. However, a lack of systematic studies exists to uncover the connection between these two elements. Here, we comprehensively characterized the associations between G4s and human candidate cis -regulatory elements (cCREs) inferred from the Encyclopedia of DNA Elements (ENCODE) data. Our findings demonstrate that G4s are prominently enriched in most types of cCREs, particularly in elements with promoter-like signatures (PLS). Furthermore, we observed that the co-occurrence of CTCF signals with H3K4me3 or H3K27ac signals significantly strengthens the association between cCREs and G4s. This association becomes more pronounced when examining cell type-specific cCREs and G4s. Interestingly, compared to cCREs, genetic variants occurring in G4s, especially within their G-runs, often exhibit higher regulatory potential and deleterious effects. Runs of several consecutive guanines (G-runs) in the vicinity of transcriptional start sites tend to be more constrained in mammalian evolution than cCREs’s G-runs. Besides, the presence of G4s, is often linked to a more favorable local chromatin environment for the activation and execution of regulatory functions of cCREs, potentially attributable to the formation of G4 secondary structures. Finally, we discovered that G4-associated cCREs exhibit widespread activation in a variety of cancers. Altogether, our study suggests that G4s are integral components of human cis -regulatory elements, and the G4 primary sequences are associated with the localization of the cCREs, while the G4 structures are linked to the activation of the cCREs. Therefore, we propose to define G4s as pivotal regulatory elements in the human genome.

DNA G-quadruplexes (G4) are non-canonical DNA structures that play key roles in various biological processes. Antibody-dependent sequencing is an important tool for identifying intracellularly formed DNA G4s, and peak calling is a crucial step in processing the sequencing data. As the applicability of existing peak calling algorithms to intracellular G4 data has not been previously assessed, we systematically compared and evaluated these algorithms to determine those best suited for G4 detection. We selected seven representative candidates from 43 published peak calling algorithms for detailed evaluation. The performance of each candidate on six published intracellular G4 sequencing datasets (GSE107690, GSE145090, GSE133379, GSE178668ChIP-seq, GSE178668CUT&Tag, GSE221437) were assessed by precision and recall against customized benchmarks integrating results from multiple algorithms, as well as consistency with known G4 information (pG4 predicted by pqsfinder, oG4 from GSE63874, and multi-cell-line conserved G4s) and epigenetic signals. We identified MACS2, PeakRanger, and GoPeaks as the most effective algorithms for analyzing intracellular G4 sequencing data, and attributed their superior performance partially to the distribution model of sequencing reads/fragments used in the hypothesis testing step of the peak calling procedures. These findings provide guidance and rationale for selecting peak callers appropriate for intracellular G4 data.

G-quadruplexes (G4) are noncanonical secondary genome structures. Aberrant formation of G4s can impair genome integrity. Investigation of the relationship between G4s and somatic structural variants (SV) in cancers could provide a better understanding of the role of G4 formation in cancer development and progression. In this study, we combined bioinformatic approaches and multiomics data to investigate the connection between G4s and the somatic SVs. Somatic SV breakpoints were significantly enriched in G4 regions, regardless of SV subtypes. This enrichment was only observed in regions demonstrated to form G4s in cells (“active quadruplexes”), rather than in regions with a sequence compatible with G4 formation but without confirmed G4 formation (“potential quadruplexes”). Several genomic features affected the connection between G4s and SVs, with the enrichment being notably strengthened at the boundary of topologically associated domains. Somatic breakpoints were also preferentially associated with G4 regions with earlier replication timing and open chromatin status. In patients with cancer with homologous recombination repair defects, G4s and somatic breakpoints were substantially more strongly associated. Machine learning models were constructed that showed that G4 propensity is a potent feature for predicting the density of SV breakpoints. Altogether, these findings suggest that the G4 structures play a critical role in modulating the production of somatic SVs in cancers.Significance:G-quadruplex structure formation constitutes a critical step in the production of somatic structural variants in cancers, suggesting G-quadruplex structures as potential targets for future cancer prevention and treatment strategies.

In vivo G-quadruplex data used in this study