CoChat
Sign Up
Cited by 269DNA and Nucleic Acid ChemistryGenomics and Chromatin DynamicsEpigenetics and DNA Methylation

The DOI matches its Crossref record; 94 of the 97 references it cites were found in a registry or an index.

Promoter G-quadruplexes and transcription factors cooperate to shape the cell type-specific transcriptome

Open access
Lago, Sara, Nadai, Matteo, Cernilogar, Filippo M., Kazerani, Maryam, Domíniguez Moreno, Helena, Schotta, Gunnar, Richter, Sara N.·2021·Nature Communications
DOI 10.1038/s41467-021-24198-2PMC PMC8222265

Published version of a preprint

Open this paper in CoChat — free accountOpen PDF

Abstract

Cell identity is maintained by activation of cell-specific gene programs, regulated by epigenetic marks, transcription factors and chromatin organization. DNA G-quadruplex (G4)-folded regions in cells were reported to be associated with either increased or decreased transcriptional activity. By G4-ChIP-seq/RNA-seq analysis on liposarcoma cells we confirmed that G4s in promoters are invariably associated with high transcription levels in open chromatin. Comparing G4 presence, location and transcript levels in liposarcoma cells to available data on keratinocytes, we showed that the same promoter sequences of the same genes in the two cell lines had different G4-folding state: high transcript levels consistently associated with G4-folding. Transcription factors AP-1 and SP1, whose binding sites were the most significantly represented in G4-folded sequences, coimmunoprecipitated with their G4-folded promoters. Thus, G4s and their associated transcription factors cooperate to determine cell-specific transcriptional programs, making G4s to strongly emerge as new epigenetic regulators of the transcription machinery.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-021-24198-2

Cites

94 of 94 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Genomic position of G4s and association to gene expression. A Percentage distribution of G4 peaks in functional genomic regions according to Homer gene annotation. Percentages are normalized over the genomic abundance of each functional region. B Percentage proportion of expressed genes among the G4-containing genes (yellow). G4-depleted genes (no G4, violet) are reported as reference. One transcript per gene was considered as threshold. C Gene expression distribution expressed in transcripts per million (TPM) of all the G4-containing genes (yellow). Genes were grouped according to the functional annotation of the immunoprecipitated G4 region. G4-depleted genes (no G4, violet) are reported as reference. The box plots central line represents the median, the bottom, and upper bounds of the box represent the 25th and 75th percentile, respectively, and the whiskers represent the lowest and highest score, excluding outliers. The significance level of each gene category was calculated by two-sided T test (CI 95%) with respect to the no G4 group. ***p value < 0.001, **p value < 0.01, the absence of asterisks indicates that the difference is not statistically significant. Exact p values are the following: 5′UTR p = 3.3e−14, exon p = 0.0032, intergenic p = 1.2e−9, intron p < 2.22e−16, promoter-TSS p < 2.22e−16. Numerosity of each category is: 3′UTR n = 29, 5′UTR n = 49, exon n = 44, intergenic n = 729, intron n = 808, noncoding n = 36, promoter-TSS n = 1434, TTS n = 47, no G4 n = 23662. D Upper panel: percentage of G4-containing genes, in genes grouped according to their expression level (no expression, low, medium, or high) and their distance from the TSS of the closest gene (<1000 bp, between 1000–15,000 and >15,000 bp). Lower panel: detailed view of gene expression level (TPM) and density distribution of genes with folded G4s within 1000 bp from TSS in function of the G4 distance from the TSS. E Genomic view of representative regions showing the G4-ChIP peak position with respect to the TSS: G4-ChIP peaks in two gene promoters with noncoding upstream regions are displayed in the upper panels (METTL13 and SUCO); G4-ChIP peaks embedded in the coding regions of two adjacent genes with opposite transcription direction are shown in the lower panels (COMMD6 and UCHL3—left; CRTC-AS1 and BLM—right). Source data for each panel are provided or referenced in the Source data file.
Fig. 1 Genomic position of G4s and association to gene expression. A Percentage distribution of G4 peaks in functional genomic regions according to Homer gene annotation. Percentages are normalized over the genomic abundance of each functional region. B Percentage proportion of expressed genes among the G4-containing genes (yellow). G4-depleted genes (no G4, violet) are reported as reference. One transcript per gene was considered as threshold. C Gene expression distribution expressed in transcripts per million (TPM) of all the G4-containing genes (yellow). Genes were grouped according to the functional annotation of the immunoprecipitated G4 region. G4-depleted genes (no G4, violet) are reported as reference. The box plots central line represents the median, the bottom, and upper bounds of the box represent the 25th and 75th percentile, respectively, and the whiskers represent the lowest and highest score, excluding outliers. The significance level of each gene category was calculated by two-sided T test (CI 95%) with respect to the no G4 group. ***p value < 0.001, **p value < 0.01, the absence of asterisks indicates that the difference is not statistically significant. Exact p values are the following: 5′UTR p = 3.3e−14, exon p = 0.0032, intergenic p = 1.2e−9, intron p < 2.22e−16, promoter-TSS p < 2.22e−16. Numerosity of each category is: 3′UTR n = 29, 5′UTR n = 49, exon n = 44, intergenic n = 729, intron n = 808, noncoding n = 36, promoter-TSS n = 1434, TTS n = 47, no G4 n = 23662. D Upper panel: percentage of G4-containing genes, in genes grouped according to their expression level (no expression, low, medium, or high) and their distance from the TSS of the closest gene (<1000 bp, between 1000–15,000 and >15,000 bp). Lower panel: detailed view of gene expression level (TPM) and density distribution of genes with folded G4s within 1000 bp from TSS in function of the G4 distance from the TSS. E Genomic view of representative regions showing the G4-ChIP peak position with respect to the TSS: G4-ChIP peaks in two gene promoters with noncoding upstream regions are displayed in the upper panels (METTL13 and SUCO); G4-ChIP peaks embedded in the coding regions of two adjacent genes with opposite transcription direction are shown in the lower panels (COMMD6 and UCHL3—left; CRTC-AS1 and BLM—right). Source data for each panel are provided or referenced in the Source data file.CC BY · PubMed Central
Fig. 2 Relationship between G4s and open chromatin. A Genomic view showing input and G4 IP samples (blue tracks), ATAC-seq (purple track), and RNA-seq (gray track) peaks in the promoter of representative genes: CDK4 (upper panel), SMU1 (mid panel), and ESR1 as negative control gene (lower panel). B Venn diagram displaying the intersection between peak regions corresponding to IP G4s (light blue) and open chromatin regions (violet) mapped by ATAC-seq in promoters. C Percentage proportion of expressed genes grouped according to the presence of G4s and open chromatin signal in their promoter region. One transcript per gene was considered as expression threshold. D Expression distribution of all genes grouped according to the presence of ChIP-seq G4s and ATAC-seq signals in their promoter region. Gene expression is reported as TPM (transcript per million). In C and D, the presence and absence of G4 and ATAC-seq signals are indicated below the graphs. The box plots central line represents the median, the bottom, and upper bounds of the box represent the 25th and 75th percentile, respectively, and the whiskers represent the lowest and highest score, excluding outliers. The significance level of each category was calculated by two-sided T test (CI 95%) with respect to the G4:ATAC −/− group. ***p value < 0.001. Exact p values are <2.22e−16 for both G4:ATAC +/+ and −/+ categories. Numerosity of each category is the following: G4:ATAC +/+ n = 1351, −/+ n = 8893, −/− n = 16204. Source data for each panel are provided or referenced in the Source data file.
Fig. 2 Relationship between G4s and open chromatin. A Genomic view showing input and G4 IP samples (blue tracks), ATAC-seq (purple track), and RNA-seq (gray track) peaks in the promoter of representative genes: CDK4 (upper panel), SMU1 (mid panel), and ESR1 as negative control gene (lower panel). B Venn diagram displaying the intersection between peak regions corresponding to IP G4s (light blue) and open chromatin regions (violet) mapped by ATAC-seq in promoters. C Percentage proportion of expressed genes grouped according to the presence of G4s and open chromatin signal in their promoter region. One transcript per gene was considered as expression threshold. D Expression distribution of all genes grouped according to the presence of ChIP-seq G4s and ATAC-seq signals in their promoter region. Gene expression is reported as TPM (transcript per million). In C and D, the presence and absence of G4 and ATAC-seq signals are indicated below the graphs. The box plots central line represents the median, the bottom, and upper bounds of the box represent the 25th and 75th percentile, respectively, and the whiskers represent the lowest and highest score, excluding outliers. The significance level of each category was calculated by two-sided T test (CI 95%) with respect to the G4:ATAC −/− group. ***p value < 0.001. Exact p values are <2.22e−16 for both G4:ATAC +/+ and −/+ categories. Numerosity of each category is the following: G4:ATAC +/+ n = 1351, −/+ n = 8893, −/− n = 16204. Source data for each panel are provided or referenced in the Source data file.CC BY · PubMed Central
Fig. 3 Detection of G4s foci upon transcription perturbation treatments. A Representative fields of view showing G4 foci formation detected by immunofluorescence in control non-treated (CTRL), entinostat (2 and 4 μM) and actinomycin D (0.5 and 1 μM) treated 93T449 cells. Nuclear staining (blue), BG4 (green), and the merged channels are reported. Scale bars = 20 μm. The shown fields belong to one of two independent biological replicates. B Quantification of BG4 nuclear staining detected by immunofluorescence in control non-treated and entinostat (2 and 4 μM)—upper panel—or actinomycin D (0.5 and 1 μM) treated—lower panel. BG4 integrated fluorescence intensity within nuclei normalized by the corresponding nuclear area (μm2) is reported. The central line for each condition represents the mean ± standard deviation. Statistical significance was calculated by unpaired two-sided T test (CI 95%) with: ****p value < 0.0001, ***p value < 0.001, *p value < 0.01. Exact p values are the following: actinomycine D 0.5 μM p = 0.0214, actinomycine D 1 μM p = 0.0002, entinostat 2 μM p < 0.0001, entinostat 4 μM p < 0.0001. The number of quantified cells for each condition are the following: ctrl n = 150, entinostat 2 μM n = 147, entinostat 4 μM n = 104, actinomycin D 0.5 μM n = 74, actinomycin D 1 μM n = 128. Source data for each panel are provided or referenced in the Source data file.
Fig. 3 Detection of G4s foci upon transcription perturbation treatments. A Representative fields of view showing G4 foci formation detected by immunofluorescence in control non-treated (CTRL), entinostat (2 and 4 μM) and actinomycin D (0.5 and 1 μM) treated 93T449 cells. Nuclear staining (blue), BG4 (green), and the merged channels are reported. Scale bars = 20 μm. The shown fields belong to one of two independent biological replicates. B Quantification of BG4 nuclear staining detected by immunofluorescence in control non-treated and entinostat (2 and 4 μM)—upper panel—or actinomycin D (0.5 and 1 μM) treated—lower panel. BG4 integrated fluorescence intensity within nuclei normalized by the corresponding nuclear area (μm2) is reported. The central line for each condition represents the mean ± standard deviation. Statistical significance was calculated by unpaired two-sided T test (CI 95%) with: ****p value < 0.0001, ***p value < 0.001, *p value < 0.01. Exact p values are the following: actinomycine D 0.5 μM p = 0.0214, actinomycine D 1 μM p = 0.0002, entinostat 2 μM p < 0.0001, entinostat 4 μM p < 0.0001. The number of quantified cells for each condition are the following: ctrl n = 150, entinostat 2 μM n = 147, entinostat 4 μM n = 104, actinomycin D 0.5 μM n = 74, actinomycin D 1 μM n = 128. Source data for each panel are provided or referenced in the Source data file.CC BY · PubMed Central
Fig. 4 Comparison of 93T449 and HaCaT cell lines. A Venn diagram showing the intersection between G4 peaks found in 93T449 (yellow) and HaCaT (salmon) cell lines. B Percentage of genes containing at least one G4 in their promoter in 93T449 (upper panel) or HaCaT (lower panel) cells, in function of their differential expression in the two cell lines expressed as log2 fold change (FC). Orange and blue-violet bars in both panels correspond to genes that have higher expression in 93T449 cells and HaCaT cells, respectively. C Venn diagram showing the intersection between the ATAC-seq peaks found in promoters of 93T449 (light blue) and HaCaT (pink) cell lines. D Differential gene expression comparison of the same genes in 93T449 and HaCaT cells, based on the presence of G4s and open chromatin combinations. Orange and blue symbols indicate data for 93T449 and HaCaT cells, respectively. The presence (+) or absence (−) of G4s or ATAC signals are reported. Bars indicate gene expression distribution of the differentially expressed genes in 93T449 vs HaCaT cells evaluated by two-sided T test in comparison to the G4:ATAC ++/++ condition (CI 95%, ***p value < 0.001). Exact p values are the following: +−/+− p = 2.3e−16, −+/−+ p < 2.22e−16, −−/−+ p = 3e−8, −−/+− p = 4e−7. The box plots central line represents the median, the bottom, and upper bounds of the box represent the 25th and 75th percentile, respectively, and the whiskers represent the lowest and highest score, excluding outliers. G4:ATAC ++/++ n = 1098, +−/+− n = 61, −+/−+ n = 435, −−/−− n = 14,717, −−/−+ n = 649, −−/+− n = 2315. E Genomic view showing representative regions of 93T449 and HaCaT cell lines gene expression (RNA-seq track, gray) with respect to the presence of G4 peaks (ChIP-seq, blue) and open chromatin (ATAC-seq, purple). In particular, TTC14 gene is displayed, which shows both G4 and ATAC signals in 93T449 cells, while it shows limited accessibility and no G4 in HaCaT cells. These differences are reflected in the corresponding RNA amount, which is much lower in HaCaT cells. Source data for each panel are provided or referenced in the Source data file.
Fig. 4 Comparison of 93T449 and HaCaT cell lines. A Venn diagram showing the intersection between G4 peaks found in 93T449 (yellow) and HaCaT (salmon) cell lines. B Percentage of genes containing at least one G4 in their promoter in 93T449 (upper panel) or HaCaT (lower panel) cells, in function of their differential expression in the two cell lines expressed as log2 fold change (FC). Orange and blue-violet bars in both panels correspond to genes that have higher expression in 93T449 cells and HaCaT cells, respectively. C Venn diagram showing the intersection between the ATAC-seq peaks found in promoters of 93T449 (light blue) and HaCaT (pink) cell lines. D Differential gene expression comparison of the same genes in 93T449 and HaCaT cells, based on the presence of G4s and open chromatin combinations. Orange and blue symbols indicate data for 93T449 and HaCaT cells, respectively. The presence (+) or absence (−) of G4s or ATAC signals are reported. Bars indicate gene expression distribution of the differentially expressed genes in 93T449 vs HaCaT cells evaluated by two-sided T test in comparison to the G4:ATAC ++/++ condition (CI 95%, ***p value < 0.001). Exact p values are the following: +−/+− p = 2.3e−16, −+/−+ p < 2.22e−16, −−/−+ p = 3e−8, −−/+− p = 4e−7. The box plots central line represents the median, the bottom, and upper bounds of the box represent the 25th and 75th percentile, respectively, and the whiskers represent the lowest and highest score, excluding outliers. G4:ATAC ++/++ n = 1098, +−/+− n = 61, −+/−+ n = 435, −−/−− n = 14,717, −−/−+ n = 649, −−/+− n = 2315. E Genomic view showing representative regions of 93T449 and HaCaT cell lines gene expression (RNA-seq track, gray) with respect to the presence of G4 peaks (ChIP-seq, blue) and open chromatin (ATAC-seq, purple). In particular, TTC14 gene is displayed, which shows both G4 and ATAC signals in 93T449 cells, while it shows limited accessibility and no G4 in HaCaT cells. These differences are reflected in the corresponding RNA amount, which is much lower in HaCaT cells. Source data for each panel are provided or referenced in the Source data file.CC BY · PubMed Central
Fig. 5 Identification of TFs binding to BG4-IP regions in 93T449 cells. A Consensus sequences of TFBSs that are significantly enriched in BG4 ChIP peaks, as calculated by Homer software. AP-1: p value 1e−490, SP1: p value 1e−172. B Position and frequency of TFBSs with respect to the BG4 ChIP peak center. Data points for each TFBS occurrence are reported and fitted according to a nonparametric spline regression curve, the gray area surrounding the curve represents the confidence interval as measure of the regression likelihood. C Percentage of genes with (yellow) and without (darkblue) ChIP-G4s, oG4s, and pG4s containing validated TFBS for AP-1 and SP1, according to ENCODE database. D Western blot showing co-immunoprecipitation of G4s and the two TFs SP1 and AP-1. The INPUT lane 1 corresponds to the total fraction of the sheared chromatin used as starting material. G4s were immunoprecipitated by BG4 antibody (IP BG4), and AP-1 and SP1 were detected by immunoblotting (lane 2). AP-1 (lanes 4 and 5) and SP1 (lanes 6 and 7) TFs were immunoprecipitated from the sheared chromatin with or without previous incubation in the presence of BG4 antibody. Mock G (lane 3) and A (lane 8) are the negative controls immunoprecipitated in the absence of primary antibody using protein-G- or protein-A-coated beads, respectively. The shown blots belong to one of at least two independent biological replicates performed for each sample. Source data for each panel are provided or referenced in the Source data file.
Fig. 5 Identification of TFs binding to BG4-IP regions in 93T449 cells. A Consensus sequences of TFBSs that are significantly enriched in BG4 ChIP peaks, as calculated by Homer software. AP-1: p value 1e−490, SP1: p value 1e−172. B Position and frequency of TFBSs with respect to the BG4 ChIP peak center. Data points for each TFBS occurrence are reported and fitted according to a nonparametric spline regression curve, the gray area surrounding the curve represents the confidence interval as measure of the regression likelihood. C Percentage of genes with (yellow) and without (darkblue) ChIP-G4s, oG4s, and pG4s containing validated TFBS for AP-1 and SP1, according to ENCODE database. D Western blot showing co-immunoprecipitation of G4s and the two TFs SP1 and AP-1. The INPUT lane 1 corresponds to the total fraction of the sheared chromatin used as starting material. G4s were immunoprecipitated by BG4 antibody (IP BG4), and AP-1 and SP1 were detected by immunoblotting (lane 2). AP-1 (lanes 4 and 5) and SP1 (lanes 6 and 7) TFs were immunoprecipitated from the sheared chromatin with or without previous incubation in the presence of BG4 antibody. Mock G (lane 3) and A (lane 8) are the negative controls immunoprecipitated in the absence of primary antibody using protein-G- or protein-A-coated beads, respectively. The shown blots belong to one of at least two independent biological replicates performed for each sample. Source data for each panel are provided or referenced in the Source data file.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 97 references

Research works printed but not linked: 1

  1. Hänsel-Hertsch, R. et al. G-quadruplex structures mark human regulatory chromatin. Nat Genet. 48, 1267–1272 (2016). · DOI · Paper page

Cited by

Sorted by citations

Items 201–204 of 204

Metadata

Authors
Lago, S., Nadai, M., Cernilogar, F. M., Kazerani, M., Domíniguez Moreno, H., Schotta, G., & Richter, S. N.
Affiliations
  • University of Padua (IT)
  • Ludwig-Maximilians-Universität München (DE)
Corresponding author
  • Lago, Sara
  • Schotta, Gunnar
  • Richter, Sara N.
ORCID
Venue
Nature Communications
Publisher
Springer Science and Business Media LLC
Published
2021-06-23
Volume
12
Issue
1
Pages
3885-3885
ISSN
2041-1723 (ISSN-L)
DOI
10.1038/s41467-021-24198-2
PMC
PMC8222265
OpenAlex
W3173676343
Cited by
269
References
97
Citations by year
2026: 26 · 2025: 53 · 2024: 72 · 2023: 65 · 2022: 41 · 2021: 10
FWCI
14.76
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Communications (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Padua Research Archive (University of Padova) (submittedVersion, cc-by)
Funders
  • Deutsche Forschungsgemeinschaft (213249687, SFB1064, SFB1321, 329628492)
  • Associazione Italiana per la Ricerca sul Cancro (21850)
  • Fondazione Italiana per la Ricerca sul Cancro (21850)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › DNA and Nucleic Acid Chemistry
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Epigenetics and DNA Methylation
Keywords
G-quadruplex, transcription factors, cell type-specific transcriptome, epigenetic regulation of transcription, AP-1, Sp1, keratinocytes, chromatin accessibility
MeSH terms
Base Sequence, Binding Sites, Cell Line, DNA, Humans, Nucleic Acid Conformation, Promoter Regions, Genetic, Protein Binding, Transcription Factors, Sp1 Transcription Factor, Transcription Factor AP-1, Gene Expression Profiling, Cell Line, Tumor, G-Quadruplexes, Transcriptome

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review