Artificially inserted strong promoter containing multiple G-quadruplexes induces long-range chromatin modification
eLife assessment
This valuable study demonstrates that genomic insertion of a G4-containing sequence can be sufficient to induce chromosome loops and alter gene expression. The evidence supporting the conclusions is convincing. Effects were shown by Hi-C as well as qPCR for chromatin modifications and expression, and the specificity of the effects was controlled by mutating the G4-containing sequence or treating with LNA probes to abolish G4 structure formation. The work will be of interest to researchers working on chromatin organization and gene regulation.
https://doi.org/10.7554/eLife.96216.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Although the role of G-quadruplex (G4) DNA structures has been suggested in chromosomal looping this was not tested directly. Here, to test causal function, an array of G4s, or control sequence that does not form G4s, were inserted within chromatin in cells. In vivo G4 formation of the inserted G4 sequence array, and not the control sequence, was confirmed using G4-selective antibody. Compared to the control insert, we observed a remarkable increase in the number of 3D chromatin looping interactions from the inserted G4 array. This was evident within the immediate topologically associated domain (TAD) and throughout the genome. Locally, recruitment of enhancer histone marks and the transcriptional coactivator p300/Acetylated-p300 increased in the G4-array, but not in the control insertion. Resulting promoter-enhancer interactions and gene activation were clear up to 5 Mb away from the insertion site. Together, these show the causal role of G4s in enhancer function and long-range chromatin interactions. Mechanisms of 3D topology are primarily based on DNA-bound architectural proteins that induce/stabilize long-range interactions. Involvement of the underlying intrinsic DNA sequence/structure in 3D looping shown here therefore throws new light on how long-range chromosomal interactions might be induced or maintained.
Introduction
G-quadruplexes (G4s), non-canonical DNA secondary structures with quartets of Guanines bonded by Hoogsteen base pairing, are instrumental in regulating gene expression (Sengupta et al., 2021; Varshney et al., 2020). G4s were primarily observed to be able to regulate gene expression when present around transcription start sites (TSSs; Huppert and Balasubramanian, 2007; Rawal et al., 2006; Verma et al., 2008). G4s can regulate gene expression by directly regulating recruitment of transcription factors and RNA polymerase or via alteration of DNA accessibility by modulating the epigenetic state of the gene promoters (Hussain et al., 2017; Kumar et al., 2011; Lago et al., 2021; Mukherjee et al., 2019; Saha et al., 2017; Sharma et al., 2021; Varshney et al., 2020). Recent studies have implicated the role of G4s in long-distance gene regulation (Robinson et al., 2021).
High-throughput chromosome conformation capture techniques reveal that specific regions of the human genome interact in three dimensions (3D) via chromatin looping and formation of topologically associated domains (TADs; Bonev and Cavalli, 2016; Denker and de Laat, 2016; Roy et al., 2018). Interestingly, recent computational studies observed enrichment of G4s in TAD boundaries along with higher enrichment of architectural proteins like CTCF and cohesin (Hou et al., 2019). Further, multiple studies noted the presence of G4s correlated with enhancer histone marks like H3K27Ac and H3K4Me1, and predominantly open chromatin regions (Calo and Wysocka, 2013; Hou et al., 2021; Shlyueva et al., 2014).
Although these studies implicate the role of G4s in long-range interactions and/or enhancer function, this was not directly tested. Here we asked if G4 structures might directly alter 3D chromatin, and affect long-range interactions including the epigenetic state of chromatin. To address this, we inserted an array of G4s into an isolated locus devoid of G4-forming sequences using CRISPR-Cas9 genome editing. To evaluate the specific function of G4s, a similar sequence of identical length but devoid of G4-forming capability was introduced. Using these pair of cell lines, we observed insertion of G4s specifically led to the recruitment of enhancer histone marks and increased expression of genes in a 10 Mb window. 3 C and Hi-C results showed induced long-range interactions throughout the genome affecting topologically associated domains (TADs) that were specifically due to the incorporated G4s, and not found in case of the control insertion.
Results
Insertion of an array of G4s in an isolated locus
First, we sought to insert an array of G4s in a relatively isolated locus. We looked into Hi-C data from Rao et al., 2014 and identified a region that was markedly isolated with little or no interaction with its surrounding regions (as shown by snapshots of Hi-C interaction matrices obtained using the 3D genome browser Wang et al., 2018 in Figure 1—figure supplement 1). In addition, this region was devoid of any G4s in the vicinity (no G4 forming motifs in a±2.5 kb window). Thereafter, we artificially inserted an array of G4 forming sequences (275 bp long) at this region near the 79 millionth position of chromosome 12 (79M in following text, chr12:79872423–79872424, hg19 genome assembly) using CRISPR-Cas9 genome editing (Figure 1A, Figure 1—figure supplement 2). To study specific effects due to G4s, if any, a control sequence of identical length was inserted in HEK293T cells at the same locus where specific G/Cs necessary for G4 formation were substituted so that G4s are not formed by this sequence (G4-mutated control, Figure 1A, Figure 1—figure supplement 2); we also ensured that the GC content was minimally affected by the substitutions (72.4% from 76.73%). Homozygous insertion was confirmed by PCR using primers adjacent to the insertion site followed by Sanger sequencing (Figure 1B, Figure 1—figure supplement 3). The array of G4-forming sequences used for insertion was previously reported to form stable G4s in human cells (Lim et al., 2010; Monsen et al., 2020; Palumbo et al., 2009; Sharma et al., 2021).
Chromatin epigenetic landscape upon insertion of G4s
To understand how the formation of G4s altered the local chromatin, chromatin immunoprecipitation (ChIP) of different chromatin-modifying histone marks was done followed by qRT-PCR using primers spanning the inserted locus. PCR primers were designed such that none of the primers bind to any site of G/C alteration in the mutated control insert; either the forward/reverse primer is from the adjacent region for specificity; covers adjacent regions for studying any effects on chromatin; and, PCRs optimized keeping in mind the repeats within the inserted sequence. Given these, primer pairs R1-R4 were chosen for further work following optimizations (Figure 2, top panel). For G4 formation within cells by the G4-array insert sequence we used the reported G4 antibody BG4 (Hänsel-Hertsch et al., 2016). Using primer pairs R2, covering >100 bases of the inserted G4-array, or the G4-mutated control, BG4 ChIP followed by qPCR was performed. Significant BG4 binding was clear in the G4-array insert, and not in the G4-mutated insert, demonstrating formation of G4s by the inserted G4-array (Figure 2—figure supplement 1).
We observed significant increase in H3K4Me1 and H3K27Ac enhancer marks in the G4-array when compared to the G4-mutated control (Figure 2A and B). However, there was no G4-specific change in the presence of chromatin compaction marks, H3K27Me3 and H3K9Me3, or the promoter activation mark H3K4Me3 (Figure 2C, D and E). The G4-dependent recruitment of H3K4Me1 (associated with enhancers Heintzman et al., 2009; Heintzman et al., 2007) and H3K27Ac (associated with active enhancers and promoters Creyghton et al., 2010; Heintzman et al., 2009; Heintzman et al., 2007) indicated enhancer-like characteristics of the inserted G4s.
Enhancer-like features emerged upon insertion of G4s
We next asked how the insertion of the G4-array influenced the expression of surrounding genes. To understand the distance-dependent gene regulatory impacts of the inserted G4-array, the mRNA expression of the nearest three genes and then some arbitrarily chosen genes further away up to 5 megabases (Mb) both up and downstream from the insertion site was quantified. Notably, the expression of four of the tested genes (PAWR, PPP1R12A, NAV3, and SLC6A15) increased in the G4-array insert compared to the mutated insert control cells (Figure 3A). Based on this enhanced expression, we further tested and observed a somewhat concomitant increase in the recruitment of Ser5 phosphorylated RNA Pol II in the surrounding gene promoters (Figure 3B). Next, we tested if chromosomal looping interactions between the insertion site and the gene promoters were involved in these long-distance effects by using chromosome conformation capture (3C). The 3C assay between the insertion locus and the gene promoters could only be performed till the NAV3 promoter 1.6 Mb away. Beyond this distance, there was not any significantly detectable PCR amplification of 3C interaction products. The 3C assays revealed that there was a G4-dependent increase in chromosomal looping interactions of the insertion locus with the gene promoters (Figure 3C). These results suggested that the inserted G4-array sequence was acting like an enhancer element.
To understand the mechanism behind the enhancer-like property of the inserted G4-array we analyzed the recruitment of transcriptional coactivator p300 (Kalkhoven, 2004). There was a relatively modest increase in the recruitment of p300, and a more substantial increase in the recruitment of the more functionally active acetylated p300/CBP, was seen within the G4-array when compared against the mutated control (Figure 3D and E). Together, these results supported the enhancer-like function of the inserted G4-array.
LNA-mediated disruption of the inserted G4s reverses enhancer phenotype
To further establish that the enhancer effects upon the G4 array insertion are due to the formation of G4s, we wanted to see if some of the effects observed could be reversed upon disrupting the inserted structures. Specific Locked Nucleic Acid (LNA) probes were designed to target and disrupt the G4 using a similar approach as shown by others (Cadoni et al., 2021; Chowdhury et al., 2022; Kumar et al., 2008). Three probes were designed with stretches of mostly cytosines (Cs) as LNAs which would hybridize with stretches of guanines (Gs) in the G4-array insert important for the structure formation (Figure 4A; see Methods). We observed that there was a significant decrease in the expression of PPP1R12A and NAV3, two of the genes initially observed to have G4-dependent enhanced expression (Figure 3A), when the G4 array inserted cells were treated with the G4 targeting LNAs (Figure 4B). As expected, although modest, a decrease in the H3K4Me1 and H3K27Ac enhancer histone modifications was evident within the insert upon the LNAs treatment (Figure 4C and D). As a control experiment, we next tested whether the LNA probes affected surrounding gene expression in the G4-mutated insert cells. Changes in the expression of the genes were not significant across replicates in case of G4-mutated insert cells (Figure 4—figure supplement 1). Together these confirmed the decrease in the expression of PPP1R12A and NAV3 in the G4-array insert upon LNA treatment was likely specific to G4 disruption. These indicate that the disruption of the inserted G4s can reverse the enhancer functions observed upon G4 insertion, further supporting the role of the G4 structure in enhancer functions.
Domain-wide increase in looping interactions by G4s
For in-depth analysis of the long-range changes in chromatin architecture upon G4 insertion, we performed genome-wide interaction by Hi-C. First, we compared all the Hi-C contacts originating within a±10 kb window comprising the G4-array insert, or the G4-mutated control insert. Compared to the mutated control, the G4-inserted locus had more than twice as many genome-wide Hi-C interactions (6390 vs 3133; Figure 5A and B, Supplementary file 1). To rule out the possibility of artifacts due to the insertion we independently analyzed Hi-C data in HEK293T cells reported earlier (taken from GSE44267, Zuin and Dixon, 2014). After normalizing for sequencing depth, the number of Hi-C contacts from the same window in HEK293T was relatively similar to the G4-mutated insert control (3968 and 3133 respectively, Figure 5C, Supplementary file 1). Together, these showed that a significant number of new long-range interactions were induced throughout the genome due to the inserted G4s, but not from the inserted control sequence.
For closer analysis, we focused on intrachromosomal Hi-C interaction matrices of the G4-array insert, or the mutated control insert. This was centered on the insertion locus on chromosome 12 (chr12:7,80,72,423–8,16,72,423; insertion site marked with arrows in Figure 6A and B). The number of Hi-C interactions in the G4-array insert was clearly enriched compared to the G4-mutated insert control, as expected from the global Hi-C contacts noted above. We noted that while the interactions from the G4-array insert were significantly more, the insertion per se did not affect the overall domain architecture, which was largely similar between G4 or G4-mutated inserts as clear from Figure 6A and B. Further, we asked if the domain architecture was retained from that seen in HEK293T cells (with no insertion): Comparison using reported HiC data for the same region from HEK293T cells showed this to be the case confirming that the chromatin domain architecture remained relatively unchanged on introducing the G-array or G4-mutated regions (Figure 6—figure supplement 1).
-
Figure 6—source data 1
- https://cdn.elifesciences.org/articles/96216/elife-96216-fig6-data1-v1.xlsx
-
Figure 6—source data 2
- https://cdn.elifesciences.org/articles/96216/elife-96216-fig6-data2-v1.zip
-
Figure 6—source data 3
- https://cdn.elifesciences.org/articles/96216/elife-96216-fig6-data3-v1.zip
To evaluate the effect of G4s in more detail, we plotted a Hi-C heatmap to show the enhanced or reduced (differential) contacts in the G4-array insert compared to the G4-mutated insert control cells (Figure 6C; relatively enriched/reduced contacts in the G4-array insert w.r.t. the G4-mutated insert plotted in red or blue, respectively; using Juicebox for analysis). This clearly showed that the G4-array induced significantly more Hi-C interactions; interestingly this was particularly evident in the downstream regions. For a closer analysis, we mapped the interaction frequency in a±100 kb window centered on the insertion site. This clearly showed the difference in the number of interactions between the upstream region vis-a-vis the region downstream of the insertion (Figure 6D).
To further confirm we used an independent HiC analysis method, HOMER (Hypergeometric Optimization of Motif EnRichment, Heinz et al., 2018) to compute the enhanced/reduced long-range interactions in the G4-array insert, compared to the control G4-mutated insert. Differential analysis using HOMER showed that the inserted locus induced significantly higher number of interactions in the case of G4-array insert relative to the control G4-mutated case (Figure 6E). When we plotted the significantly different chromosomal interactions with minimum 20 interaction reads, it was again clear that the number of interactions with the G4-array insertion region was significantly enhanced in the downstream region relative to the upstream (Figure 6F).
Together these show a clear role of G4s in inducing long-range interactions. A similar sequence devoid of G4-forming capability did not induce such interactions. Furthermore, the overall nature of the TAD was not disturbed, and largely consistent with what is noted in cells with no insertion. Overall, these support that the insertion of G4s induced long-range interactions with minimal organizational changes in the 3D chromatin domain, underlining the molecular role of G4s in the arrangement of 3D chromatin.
A second significant feature was notable at the insertion locus. The number of induced long-range interactions was more significant downstream of the insertion site, compared to the upstream region (Figure 6C, D and F). A close look at the Hi-C contact matrices indicated that the site of insertion was very close and downstream to the TAD boundary (Figure 6A–C). We reasoned that the G4-dependent long-range interactions were largely within the TAD, and limited in the upstream region due to the TAD boundary. This is clearly seen in Figure 6C, akin to an ‘architectural stripe’ displaying that the inserted G4 array had enhanced Hi-C interactions across the domain, thus prominently featured in the downstream regions.
G4-array insertion at a second locus gives enhancer-like functions
Finally, we checked if enhancer-like effects were observed upon insertion of G4 array at another locus. Like the first site of insertion, we first identified an isolated locus devoid of G4s in the vicinity and with low interactions with surrounding regions near the 10 millionth position of chromosome 12 (10M hereafter, chr12:10588429–10588430, hg19; Figure 7—figure supplement 1). The G4-array, or its G4-mutated (control) sequences were inserted at the 10 M locus (Figure 7A and B).
As for the 79 M locus, to validate intracellular G4 formation and study the chromatin state at the inserted locus, PCR primers were designed keeping multiple points in mind (as described above). Here, for testing formation of G4 at the 10 M insertion, we used primer pairs R2 (scheme for 10 M shown in Figure 7 top panel), covering >100 bases of the inserted G4-array, or the G4-mutated control. BG4 ChIP-qPCR validated formation of intracellular G4s within the G4 array, and not the G4-mutated control sequence (Figure 7—figure supplement 2). Next, we checked for changes in chromatin and the surrounding gene expression due to G4 formation. A relative increase in the H3K4Me1 and H3K27Ac enhancer marks in the G4-array was evident compared to the G4-mutated control (Figure 7C and D), consistent with earlier observations following G4 insertion at the 79 M locus (Figure 2A and B). We noticed, however, that the enhanced levels of H3K27Ac were not as marked as the 79 M locus. On the other hand, interestingly, relative increase in the H3K27Me3 repressor mark compared to the control mutated-G4 insert, particularly at the downstream end of the insertion locus was seen (Figure 7E). There was no G4-specific change in the presence of the chromatin compaction mark H3K9Me3, or the promoter activation mark H3K4Me3 (Figure 7F and G). As expected from earlier observations and the enhancer histone marks, there was a G4-dependent increase in the expression of surrounding genes KLRC2, KLRC1 and NTF3; except for PTPRO, which had reduced expression (Figure 7H). Taken together, G4-specific chromatin changes were evident at the 10 M locus consistent with the 79 M locus. Notable variations however must be pointed out: like the presence of the H3K27Me3 repressor histone mark, along with H3K27Ac/H3K4Me1 enhancer histone marks, indicating a poised enhancer-like state as described earlier (Calo and Wysocka, 2013). These suggest the impact of G4 formation on chromatin is likely context-specific, that is, dependent on the chromatin state of the adjacent regions.
Discussion
To directly test if G4s affect long-range chromatin organization we artificially inserted an array of G4s in the chromatin. Hi-C experiments clearly showed an enhanced number of cis- and trans-chromosomal long-range interactions emanating from the introduced G4s. This was G4-specific because a similar sequence devoid of G4-forming capability introduced at the same site did not result in enhanced interactions. Furthermore, interestingly, most new long-range interactions following G4 incorporation were downstream from the site of insertion. This is likely because the G4 insertion locus was proximal to the upstream TAD boundary thereby restricting most new interactions to the downstream regions within the TAD (Figures 5 and 6).
The insertion of the G4 array led to enhanced expression of genes up to 5 Mb away compared to cells with the G4 mutated control insertion. Furthermore, there was enrichment in the H3K4Me1 and H3K27Ac enhancer histone marks, along with recruitment of transcriptional coactivator p300 and more prominently the functionally active acetylated p300/CBP. This was clearly due to the introduction of G4s and not found upon the introduction of the G4-mutated control sequence. The enhancer marks were relatively reduced, although not markedly, when the inserted G4s were specifically disrupted (Figure 4).
The enhancer-like effects observed upon G4 insertion at the 79 M locus were largely replicated when the same G4-array was inserted in another locus (10 M). Interestingly, the G4-dependent increase in H3K27Me3 at the 10 M locus, in addition to the enhancer mark H3K4Me1, supports formation of a poised enhacer-like state, as described earlier (Calo and Wysocka, 2013). Although both the inserted G4 regions induced enhancer-like chromatin, notable context-specific influence, likely due to the chromatin-state at the regions adjacent to the insertion locus, were evident. Taken together, findings here directly support the function of G4s as enhancer-like elements and as factors that enhance long-range chromatin interactions. It is possible that such interactions are also contextually dependent on the type of G4 structure, in addition to the adjacent sequence context, and further studies will be necessary to elucidate these.
To ensure that the observed effects were from intracellular G4 formation, we accounted for the following while designing the experiments. First, we introduced an array of G4s and confirmed in vivo G4 formation; the inserted sequence was from the hTERT promoter region with multiple arrayed G4s (Lim et al., 2010; Monsen et al., 2020; Palumbo et al., 2009). Second, we selected an insertion locus that was otherwise devoid of intrinsic G4s in a±2.5 kb window. Third, the selected insertion locus was relatively sparse in long-range interactions. Fourth, we independently inserted a sequence of identical length (and similar GC%) which does not form G4s at the same locus (G4-mutated control). All results were compared to the G4-mutated insertion. Although the introduced mutations for the control sequences were minimal, the impacts of such mutations on the binding of specific transcription factors associated with the sequence, particularly SP1, reported to bind G4s (Raiber et al., 2012), cannot be ruled out.
Existing literature shows promoter G4s are involved in regulating gene expression (Huppert and Balasubramanian, 2007; Rawal et al., 2006; Verma et al., 2008). Additionally, G4s have been reported to regulate chromatin epigenetics through both cytosine methylation and histone modifications (Halder et al., 2010; Mao et al., 2018; Sengupta et al., 2021). Previous studies by us further show that promoter G4s regulate gene expression by recruiting histone-modifying regulatory complexes (Hussain et al., 2017; Mukherjee et al., 2019; Saha et al., 2017; Sharma et al., 2021). Here, we aimed to study how G4s affect the expression of genes far from their location, and if this was through G4-induced modifications in long-range 3D chromatin interactions.
Multiple studies have correlated the presence of G4s with long-range associations. CTCF, an architectural protein primarily involved in TAD boundary formation, was observed to bind to G4s and G4 stabilization was noted to enhance CTCF occupancy (Tikhonova et al., 2021). In addition, G4s were noted to be enriched in TAD boundaries and associated with the formation of chromatin loops (Hou et al., 2019). G4s were also found to coincide with open chromatin regions and H3K27Ac and H3K4Me1 ChIP-Seq peaks, which are markers for transcriptional enhancers (Calo and Wysocka, 2013; Hou et al., 2021; Lyu et al., 2022; Shlyueva et al., 2014). Most of these regions were observed to overlap with annotated enhancers and promoters regulated by such enhancers were enriched in G4s (Williams et al., 2020). A recent G4 CUT&Tag study further noted G4 formation at both active promoters and active and poised enhancers (Lyu et al., 2022).
Further, it was proposed that inter-molecular G4 formation between distant stretches of Gs may lead to DNA looping (Hegyi, 2015; Liano et al., 2022). Consistent with this, we noted with interest promoters of the four genes (PAWR, PPP1R12A, NAV3 and SLC6A15; Figure 3A), activated on long-range interaction with the inserted loci, harbor potential G4-forming sequences (pG4) (Figure 3—figure supplement 1). Further, we analyzed the long-range contact regions shown in Figure 6F, along with the whole locus, for pG4s (Hon et al., 2017). Relative enrichment in pG4s was evident, particularly within the significantly enhanced contact points, at times spreading beyond the interacting region (Figure 6F, lower panel). Together, these support G4-induced long-range interactions.
The YY1 transcription factor was found to bind to G4s and dimerization of G4-bound YY1 led to chromatin looping interactions and consequent regulation of target gene expression (Li et al., 2021). A recent study also showed R-loops and possibly R-loop-associated G4 formation are enriched at CTCF binding sites, and stronger CTCF binding facilitated by G4s promotes chromatin looping (Wulfridge et al., 2023). In addition, it was shown that G4s assist in RNA polymerase II-associated chromatin looping (Yuan et al., 2023). In this context, further work will be required to understand whether and how formation of R-loops or RNA-DNA hybrid G4s (Fay et al., 2017), and/or association of factors like cohesion and CTCF, at the G4-array insertion sites impact chromatin looping.
In summary, our findings here demonstrate a causal role of G4s in inducing both long-range associations and enhancer function. Findings from a G4-forming stretch inserted at two independent loci illustrate the function of G4s in 3D gene regulation. Together these shed new mechanistic light on how DNA secondary structure motifs directly control the state of 3D chromatin and thereby biological function.
Materials and methods
Cell lines and cell culture conditions
Request a detailed protocolHEK293T cells were procured from the NCCS cell repository, the cell identity was authenticated using STR profiling, and the cells were tested negative for mycoplasma contamination. The cells were cultured in Dulbecco’s Modified Eagle’s Medium- High Glucose (DMEM-HG) supplemented with 10% FBS and 1XAnti-Anti (Gibco).
Primary antibodies
Request a detailed protocolHistone H3 rabbit polyclonal (Abcam ab1791), H3K4Me1 rabbit polyclonal (Abcam ab8895), H3K27Ac rabbit polyclonal (Abcam ab4729), H3K4Me3 mouse monoclonal (Abcam ab1012), H3K27Me3 mouse monoclonal (Abcam ab6002), H3K9Me3 rabbit polyclonal (Abcam ab8898), p300 rabbit monoclonal (CST 54062), Ac-p300/CBP rabbit polyclonal (CST 4771), BG4 antibody (Sigma-Aldrich MABE917).
Genomic insertions using CRISPR-Cas9 genome editing
Request a detailed protocolFor the genomic insertions CRISPR-Cas9 genome editing technique was used (Ran et al., 2013). For the G4 array insertion, 275 bp long hTERT promoter region was PCR amplified from HEK 293T genomic DNA. For the insertion of the mutated G4s, a synthetic DNA template was synthesized and cloned into pUC57 vector by Genscript Biotech Corp, where 12 Gs were substituted with Ts (see Supplementary file 2 for detailed sequences). Both the G4 array and the G4 mutated insertion templates were PCR amplified using longer primers where the short homology arms were introduced as overhangs of the primer for the accurate insertion at the 79 M locus via homologous recombination (see Supplementary file 2 for primer sequences) (Paix et al., 2017). For cleavage at the 79 M locus (chr12:79,872,423–79,872,424 (hg19)), the gRNA sequence, 5’-ACTATGTATGTACATCCAGG-3’, was cloned into the pX459 v2.0, a gift from Feng Zhang, that co-expresses cas9 protein and the gRNA. For cleavage at the 10 M locus (chr12:10,588,429–10,588,430 (hg19)), the gRNA sequence, 5’-ATCCTTCCCTGAATCATCAA-3’, was used. Guide RNAs (gRNAs) were designed using the CRISPOR tool (Haeussler et al., 2016). Once the gRNA cloned vector and the insertion donor templates were ready, they were transfected into HEK293T cells and the transfected cells were selected using puromycin, whose resistance gene was present in the pX459 vector. Then these selected cells were serially diluted to isolate clones originating from single cells. Many such clones were screened to detect cells with homozygous/heterozygous insertion of the G4 array or mutated G4 insert by performing locus-specific PCR. Either primers adjacent to the insertion site or cross primers, i.e., one primer within the insert and another from the adjacent region, were used to screen and identify insertions. While using adjacent primers, a shift in PCR product with an increase in amplicon size by 275 bp (size of the insert) indicated successful insertion (see Supplementary file 2 for primer sequences).
ChIP (chromatin immunoprecipitation)
Request a detailed protocolChIP assays were performed as per the protocol previously reported in Mukherjee et al., 2018. Immunoprecipitation was done using relevant primary antibodies. IgG was used for isotype control. Total histone H3 was used as a control for the histone modifications ChIP. Three million cells were harvested and crosslinked with ~1% formaldehyde for 10 min and lysed. Chromatin was sheared to an average size of ∼250–500 bp using Biorupter (Diagenode). Ten percent of sonicated fraction was processed as input using phenol–chloroform and ethanol precipitation. ChIP was performed using 3 μg of the respective antibody incubated overnight at 4 °C. Immune complexes were collected using salmon sperm DNA-saturated magnetic protein G Dynabeads (Anti-FLAG M2 magnetic beads for BG4 ChIP) and washed extensively using a series of low salt, high salt and LiCl Buffers. The Dynabeads were then resuspended in TE (Tris- EDTA pH 8.1) buffer and treated with proteinase K at 65 ° C for ~5 hrs. Then, phenol-chloroform-isoamyl alcohol was utilized to extract DNA. Extracted DNA was precipitated by centrifugation after incubating overnight at –20 ° C with isopropanol, 0.3 M sodium acetate and glycogen. The precipitated pellet was washed with freshly prepared 70% ethanol and resuspended in TE buffer. ChIP DNA was analyzed by qRT-PCR method. See Supplementary file 2 for primer sequences.
Real-time PCR for gene (mRNA) expression
Request a detailed protocolTotal RNA was isolated using TRIzol Reagent (Invitrogen, Life Technologies) according to the manufacturer’s instructions. RNA was quantified and cDNA was synthesized using iScript cDNA Synthesis Kits. A relative transcript expression level for genes was measured by quantitative real-time PCR using a SYBR Green based method (see Supplementary file 2 for primer sequences). Average fold change was calculated by the difference in threshold cycles (Ct) between test and control samples. GAPDH gene was used as internal control for normalizing the cDNA concentration of each sample.
Chromosome conformation capture (3C)
Request a detailed protocolChromosome Conformation Capture (3 C) assay was done as per the protocol reported in Cope and Fraser, 2009 with certain modifications. Briefly, about 5–6 million cells were crosslinked using 1% formaldehyde for 10 min and then lysed to isolate the nuclei. Nuclei were digested overnight by HindIII and then ligated in a diluted reaction so that intramolecular ligation is favored. After ligation, the reaction mixture was treated with proteinase K at 65 °C to de-crosslink the DNA, followed by RNase A treatment. Then, phenol-chloroform-isoamyl alcohol was utilized to extract DNA. Extracted DNA was precipitated by centrifugation after incubating overnight at –80 °C with 70% ethanol, 0.1 M sodium acetate and glycogen. The precipitated pellet was washed with freshly prepared 70% ethanol and resuspended in TE buffer. 3C looping interactions were analyzed by TaqMan qRT-PCR method. For comparison, each interaction frequency was normalized to the interaction between exons 2 and 8 of the human α-actin (ACTA2)(Hadjur et al., 2009). See Supplementary file 2 for primer sequences.
G4 disruption using LNA probes
Request a detailed protocolProbes were designed to specifically bind to regions of genomic DNA containing G repeats which would form the G stems of the G4 structure. The probes containing LNA nucleotides should hybridize with the target with higher stability than the stability of the G4 structure thus destabilizing the G4. The probes used to target the G4 array insert were: 5’-c*ccgacccctcc*c-3’, 5’-c*cagccccctcc*g-3’, 5’-c*ccctccccttc*c-3’. Stretches of three or more Cs are shown in bold, LNA nucleotides within the probes are underlined, the ends of the probes were protected using phosphorothioate bonds, shown as *. Approximately 0.8 μg of LNA probes (all three mixed in equimolar amounts) were transfected per million cells. Cells were treated with the LNA probes for 108 hr by transfecting thrice with a gap of 36 hr in between. Scheme 1 shows the LNA probes designed to disrupt the inserted G4 structures along with the inserted G4 array sequence to show the specific sites of hybridization by the LNA probes.
Hi-C
Request a detailed protocolHi-C was performed using the Arima-HiC Kit as per the manufacturer’s protocol. After the proximally-ligated Hi-C templates were generated, sequencing libraries were prepared using NEBNext Ultra II DNA Library Prep Kit as per the Arima-HiC Kit’s protocol. The quality of the sequencing libraries was cross-checked using TapeStation (Agilent Technologies) and the KAPA Library Quantification Kit (Roche) before proceeding with sequencing using NovaSeq 6000 (Illumina).
Hi-C data analysis
Request a detailed protocolHi-C reads were mapped to the hg19 human genome and processed using default parameters using Juicer (https://github.com/aidenlab/juicer, Aiden Lab, 2023a; Durand et al., 2016b). Hi-C count matrices were generated at 5 kb, 10 kb, 25 kb, 50 kb, 100 kb, and 250 kb using Juicer. Hi-C heatmap figures were rendered using Juicebox (https://github.com/aidenlab/Juicebox/wiki/Download, Aiden Lab, 2023b; Durand et al., 2016a). Hi-C contacts originating in the loci flanking the G4 insertion site were generated using bedtools (https://bedtools.readthedocs.io/en/latest/, Quinlan, 2023). The circos plots were rendered using Circos (https://circos.ca/, Krzywinski, 2009). To identify significant interaction the data was processed using HOMER (http://homer.ucsd.edu/homer/, Benner, 2024) using analyzeHiC function. The bins showing 2-fold enrichment in G4 WT over G4 Mut and vice-versa were retained for filtering contacts for representation on circos plots.
Materials availability
Request a detailed protocolHEK 293T cells with the hTERT promoter G4 array or the G4-mutated control insertions are available upon request. Such requests can be directed to the corresponding author.
Data availability
The sequencing data underlying this article are available in the NCBI Sequence Read Archive, accessible using the following link- https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1048044. The rest of the data are available within the manuscript and supporting files. Further inquiries can be directed to the corresponding author.
-
NCBI BioProjectID PRJNA1048044. Inserted G-quadruplexes enhance chromosomal looping interactions.
References
-
SoftwareHOMER: Hypergeometric optimization of motif enrichment, version v5.1University of California, San Diego.
-
Organization and function of the 3D genomeNature Reviews Genetics 17:661–678.https://doi.org/10.1038/nrg.2016.112
-
Beyond small molecules: targeting G-quadruplex structures with oligonucleotides and their analoguesNucleic Acids Research 49:6638–6659.https://doi.org/10.1093/nar/gkab334
-
Modification of enhancer chromatin: what, how, and why?Molecular Cell 49:825–837.https://doi.org/10.1016/j.molcel.2013.01.038
-
Short LNA-modified oligonucleotide probes as efficient disruptors of DNA G-quadruplexesNucleic Acids Research 50:7247–7259.https://doi.org/10.1093/nar/gkac569
-
The second decade of 3C technologies: detailed insights into nuclear organizationGenes & Development 30:1357–1382.https://doi.org/10.1101/gad.281964.116
-
RNA G-quadruplexes in biology: principles and molecular mechanismsJournal of Molecular Biology 429:2127–2147.https://doi.org/10.1016/j.jmb.2017.05.017
-
Guanine quadruplex DNA structure restricts methylation of CpG dinucleotides genome-wideMolecular bioSystems 6:2439–2447.https://doi.org/10.1039/c0mb00009d
-
G-quadruplex structures mark human regulatory chromatinNature Genetics 48:1267–1272.https://doi.org/10.1038/ng.3662
-
G-quadruplexes in promoters throughout the human genomeNucleic Acids Research 35:406–413.https://doi.org/10.1093/nar/gkl1057
-
CBP and p300: HATs for different occasionsBiochemical Pharmacology 68:1145–1155.https://doi.org/10.1016/j.bcp.2004.03.045
-
YY1 interacts with guanine quadruplexes to regulate DNA looping and gene expressionNature Chemical Biology 17:161–168.https://doi.org/10.1038/s41589-020-00695-1
-
Long-range DNA interactions: inter-molecular G-quadruplexes and their potential biological relevanceChemical Communications 58:12753–12762.https://doi.org/10.1039/d2cc04872h
-
Coexistence of two distinct G-quadruplex conformations in the hTERT promoterJournal of the American Chemical Society 132:12331–12342.https://doi.org/10.1021/ja101252n
-
Genome-wide mapping of G-quadruplex structures with CUT&TagNucleic Acids Research 50:e13.https://doi.org/10.1093/nar/gkab1073
-
DNA G-quadruplex structures mold the DNA methylomeNature Structural & Molecular Biology 25:951–957.https://doi.org/10.1038/s41594-018-0131-8
-
The hTERT core promoter forms three parallel G-quadruplexesNucleic Acids Research 48:5720–5734.https://doi.org/10.1093/nar/gkaa107
-
Formation of a unique end-to-end stacked pair of G-Quadruplexes in the hTERT core promoter with implications for inhibition of telomerase by G-quadruplex-interactive ligandsJournal of the American Chemical Society 131:10878–10891.https://doi.org/10.1021/ja902281d
-
A non-canonical DNA structure is A binding motif for the transcription factor SP1 in vitroNucleic Acids Research 40:1499–1508.https://doi.org/10.1093/nar/gkr882
-
Genome engineering using the CRISPR-Cas9 systemNature Protocols 8:2281–2308.https://doi.org/10.1038/nprot.2013.143
-
DNA G-quadruplex structures: more than simple roadblocks to transcription?Nucleic Acids Research 49:8419–8431.https://doi.org/10.1093/nar/gkab609
-
Epigenetic suppression of human telomerase (hTERT) is mediated by the metastasis suppressor NME2 in a G-quadruplex-dependent fashionThe Journal of Biological Chemistry 292:15205–15215.https://doi.org/10.1074/jbc.M117.792077
-
Transcriptional enhancers: from properties to genome-wide predictionsNature Reviews. Genetics 15:272–286.https://doi.org/10.1038/nrg3682
-
Dna g‐quadruplexes contribute to ctcf recruitmentInternational Journal of Molecular Sciences 22:7090.https://doi.org/10.3390/IJMS22137090/S1
-
The regulation and functions of DNA and RNA G-quadruplexesNature Reviews. Molecular Cell Biology 21:459–474.https://doi.org/10.1038/s41580-020-0236-x
-
G-quadruplexes associated with R-loops promote CTCF bindingMolecular Cell 83:3064–3079.https://doi.org/10.1016/j.molcel.2023.07.009
-
G-quadruplex DNA contributes to RNA polymerase II-mediated 3D chromatin architectureNucleic Acids Research 51:8434–8446.https://doi.org/10.1093/nar/gkad588
Article and author information
Author details
Funding
DBT/Wellcome Trust India Alliance (IA/S/18/2/504021)
- Shantanu Chowdhury
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the DBT/Wellcome Trust India Alliance (IA/S/18/2/504021).
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.96216. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Roy et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,063
- views
-
- 83
- downloads
-
- 2
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Chromosomes and Gene Expression
- Neuroscience
Pathogenic variants in subunits of RNA polymerase (Pol) III cause a spectrum of Polr3-related neurodegenerative diseases including 4H leukodystrophy. Disease onset occurs from infancy to early adulthood and is associated with a variable range and severity of neurological and non-neurological features. The molecular basis of Polr3-related disease pathogenesis is unknown. We developed a postnatal whole-body mouse model expressing pathogenic Polr3a mutations to examine the molecular mechanisms by which reduced Pol III transcription results primarily in central nervous system phenotypes. Polr3a mutant mice exhibit behavioral deficits, cerebral pathology and exocrine pancreatic atrophy. Transcriptome and immunohistochemistry analyses of cerebra during disease progression show a reduction in most Pol III transcripts, induction of innate immune and integrated stress responses and cell-type-specific gene expression changes reflecting neuron and oligodendrocyte loss and microglial activation. Earlier in the disease when integrated stress and innate immune responses are minimally induced, mature tRNA sequencing revealed a global reduction in tRNA levels and an altered tRNA profile but no changes in other Pol III transcripts. Thus, changes in the size and/or composition of the tRNA pool have a causal role in disease initiation. Our findings reveal different tissue- and brain region-specific sensitivities to a defect in Pol III transcription.
-
- Biochemistry and Chemical Biology
- Chromosomes and Gene Expression
The mRNA 5'-cap structure removal by the decapping enzyme DCP2 is a critical step in gene regulation. While DCP2 is the catalytic subunit in the decapping complex, its activity is strongly enhanced by multiple factors, particularly DCP1, which is the major activator in yeast. However, the precise role of DCP1 in metazoans has yet to be fully elucidated. Moreover, in humans, the specific biological functions of the two DCP1 paralogs, DCP1a and DCP1b, remain largely unknown. To investigate the role of human DCP1, we generated cell lines that were deficient in DCP1a, DCP1b, or both to evaluate the importance of DCP1 in the decapping machinery. Our results highlight the importance of human DCP1 in decapping process and show that the EVH1 domain of DCP1 enhances the mRNA-binding affinity of DCP2. Transcriptome and metabolome analyses outline the distinct functions of DCP1a and DCP1b in human cells, regulating specific endogenous mRNA targets and biological processes. Overall, our findings provide insights into the molecular mechanism of human DCP1 in mRNA decapping and shed light on the distinct functions of its paralogs.