EMBO/SourceData huggingface dataset:
Figure 2(A) Line histogram of DMR methylation levels in each sample. Kernel density estimate smoothed scatterplot comparison of DMR methylation levels in MCF7-control and (B) ZF-D3A +dox, or (C) ZF-D3A dox-wd.(D) Line histogram of UMR methylation levels in each sample. Scatterplot compari- son of UMR methylation levels in MCF7-control and (E) ZF-D3A +dox, or (F) ZF-D3A dox-wd.(G) Heat- maps of DMRs (ordered by decreasing H3K4me3 read density) showing: DNA methylation levels (±2.5 kb), differences (ΔmCG), classification as promoter or non-promoter located, CpG dinucleotide density (±2.5 kb), and histone modification ChIP-seq read density (±2.5 kb); (H) As in (G) but for the 6,000 DMRs with the lowest H3K4me3 levels, ordered by decreasing H3K4me1 read density, and DMR classification as enhancer or non-enhancer in MCF7 (ENCODE).
Figure 3(A) H3K4me3 ChIP-seq read density flanking (±1 kb, 50 bp bins) the centre of DMRs with ∆mCG >0.2 (top), and H3K4me3 ChIP-seq peaks with ∆mCG <0.2 (bottom).(B) H3K4me3 ChIP-seq normalized read counts in DMRs with ∆mCG >0.2, from H3K4me3 ChIP-seq in ZF-D3A no-dox and ZF-D3A +dox cells. Black line represents equal normalized read counts in ZF-D3A no-dox and ZF-D3A +dox cells; blue line represents a linear regression of the data.(C) Comparison of UMR ∆mCG and difference in H3K4me3 ChIP-seq normalized read counts between ZF-D3A no-dox and ZF-D3A +dox. Blue line represents a linear regression of the data.(D) DNA methylation levels of UMRs that intersect with H3K4me3 ChIP-seq peaks in ZF-D3A no-dox and ZF-D3A +dox cells, for genomic DNA or DNA immunoprecipitated with an anti-H3K4me3 antibody.(E) DNA methylation levels of UMRs that intersect with phos- pho-ser5 RNA polymerase II ChIP-seq peaks in ZF-D3A no-dox and ZF-D3A +dox cells, for genomic DNA or DNA immunoprecipitated with an anti-phospho-ser5 RNA polymerase II antibody.
Figure 4(A) Western blot measurement of ZF-D3A protein abundance upon doxycycline induction and withdrawal.(B) Box and whisker plots of the DNA methylation levels of the 121 CpG sites in the DACH1 promoter region through- out doxycycline induction and withdrawal. Whiskers indicate 1.5 times the interquartile range or the most extreme data point, whichever is lower.(C) Average DNA methylation level at the 121 CpG sites in the DACH1 promoter region throughout doxycycline induction followed by doxycycline withdrawal and cell cycle inhibition by growth in doxycycline-free media containing different cell cycle inhibitors (Lovastatin, RO-3306, thymidine).(D) DMR mCG and hmCG levels in MCF7-control and ZF-D3A +dox cells, ordered by the difference in DMR hmCG level (hmCG/CG) between MCF7-control and ZF-D3A +dox cells, and heatmap representation of H3K4me3 levels in ZF-D3A +dox cells flanking (±2.5 kb, 50 bp bins) the centre of the DMRs.(E) Heatmap of hmCG levels flanking (±2.5 kb, 50 bp bins) the centre of the DMRs, ordered by DMR hmCG level in MCF7-control cells.
Figure 5(A) Scatter plot of the difference in UMR DNA methylation levels, versus the fold change in mRNA abun- dance of UMR-associated expressed genes, between MCF7-control and ZF-D3A +dox. Point color indicates the gene differential expression significance: red indicates FDR <0.01, black indicates FDR >0.01. Trend line (green) was fitted using a generalized additive model, and contour lines represent the relative density of UMRs.(B) Distribution of the ratios of mRNA levels between ZF-D3A +dox and MCF7-control, for expressed genes that have an associated UMR with a ∆mCG >0.3. Column color indicates the class of mRNA abundance ratios between ZF-D3A +dox and MCF7-control: blue = no decrease; green = small decrease; red = decrease.(C) Representative genes with high induction of promoter DNA methylation in ZF-D3A +dox and no reduction in mRNA abundance.
Figure 6Fluorescence distribution of cells stained with RNA-FISH probes targeting (A) GAPDH or (B) PSME1 mRNA before and after FACS.mRNA quantitation after FACS sorting of (C) GAPDH by RT-qPCR, and (D) PSME1 by RNA-seq.DNA methylation levels of single CpGs at the promoter region of (E) GAPDH, and (F) PSME1, for the low (x axis) and high (y axis) expressing sorted cell populations.(G) Genome browser representation of the GAPDH (top) and PSME1 (bottom) promoter regions showing DNA methylation in MCF7-control and low and high expressing sorted cell populations.(H) Distribution of per-read methylation levels (blue line) for selected ZF-D3A +dox WGBS reads (≥10 cytosines in the CpG context) that align within robustly methylated DMRs (DMR mean methylation level [mCG/CG] of 0.4-0.6 in ZF-D3A +dox and <0.1 in MCF7-control), and average CG methylation levels from ZF-D3A +dox WGBS data of the DMR genomic intervals covered by the selected reads (red line).
Figure 1A. Left: Telomerase activity measured using ddTRAP assay. The horizontal pink lines represent the threshold used in these samples for background correction. Each dot on the ddPCR output represents a unique droplet that is either positive or negative for fluorescent signal. Right: Quantification of TRAP activity after normalization to SW39/TEL+ cells (n=3). Data information: Data are presented as mean ± SD. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise. PTM, post-translational modification.B. hTR and hTERT expression, normalized to ACTB mRNA and to SW39/TEL+. Data information: Data are presented as mean ± SD. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise. PTM, post-translational modification.C. RNA Pol II binding at hTR promoter in SI14/TEL+, SI24/ALT+ and 8G12/ALT+ cell lines, normalized to H3 ChIP. Control loci with, respectively, high or low transcriptional activity: c-fos promoter and subtelomeric 4qHox locus (n=3). Statistical analyses were performed using one-way ANOVA with post-hoc Tukey's test for multiple comparisons. Data information: Data are presented as mean ± SD. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise. PTM, post-translational modification.D. H3K9me3 and H3K27me3 (left) or H3K4me3 and H3Ac (right) densities at hTR promoter in SI14/TEL+ and SI24/ALT+, normalized to H3 (n=3). Control loci as in C. Data information: Data are presented as mean ± SD. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise. PTM, post-translational modification.E. NFYB binding at hTR promoter in SI14/TEL+, SI24/ALT+, SW39/TEL+ and IMRB/ALT+ cells normalized to H3. Control loci with, respectively, high and low NFYB binding: hnRNPA1 promoter and subtelomeric 4qHox (n=3). Data information: Data are presented as mean ± SD. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise. PTM, post-translational modification.
Figure 2A. Western blot for hTR chaperone proteins in TEL+ and ALT+ hybrids. β-Actin is used as loading control. Quantification of western blots is shown on the right, normalized to β-Actin and to SI14/TEL+ cells (n=2-4). Statistical analyses were performed using one-way ANOVA with post-hoc Tukey's test for multiple comparisons. Data information: Data are presented as mean ± SD.B. NHP2 mRNA expression levels, normalized to ACTB mRNA and to SI14/TEL+ cells (n=3). Data information: Data are presented as mean ± SD.C. Western blot for NHP2 in SI14/TEL+, SI24/ALT+ and 8G12/ALT+ cell lines treated with DMSO or 10 μM MG132 for 6 h. β-Actin served as loading control. Quantifications are shown on the right, normalized to β-Actin and to SI14/TEL+ cells (n=5). Statistical analyses were performed with Student's two-tailed unpaired t-test. Data information: Data are presented as mean ± SD.D. hTR RNA and NHP2 mRNA expression, normalized to ACTB mRNA and to SI14/TEL+. Data information: Data are presented as mean ± SD.E. Telomerase activity measured with ddTRAP assay in NHP2-overexpressing 8G12/ALT+cells, normalized to NHP2-overexpressing SI14/TEL+. Data information: Data are presented as mean ± SD.F. (Top) Western blot for NHP2 overexpression in 8G12/ALT+ clonal cell lines. β-Actin served as loading control. SI14/TEL+ cells are shown as reference for NHP2 expression in hybrids. (Bottom) Western blot for NHP2 overexpression in SI14/TEL+ and 8G12/ALT+ cells compared to empty vector control. β-Actin served as loading control. Data information: Data are presented as mean ± SD.G. hTR RNA, hTERT and NHP2 mRNA expression, normalized to ACTB mRNA and to SI14/TEL+ (n=3). Data information: Data are presented as mean ± SD.H. RNA Pol II binding at hTR promoter, normalized to H3 ChIP (n=3). See Fig. 1C legend for control loci. Statistical analyses were performed using one-way ANOVA with post-hoc Tukey's test for multiple comparisons. Data information: Data are presented as mean ± SD.
Figure 3A. hTR RNA, hTERT and NHP2 mRNA expression, normalized to ACTB mRNA and to 143B/TEL+ cells, in a series of TEL+ and ALT+ human cell lines as indicated (n=3). Data information: Data are presented as mean ± SD. Statistical analyses were performed using one-way ANOVA with post-hoc Tukey's test for multiple comparisons.B. Western blot for NHP2 expression. β-Actin served as loading control. Quantification is shown on the right after normalization to β-Actin and to 143B/TEL+ cells (n=3). Data information: Data are presented as mean ± SD. Statistical analyses were performed using one-way ANOVA with post-hoc Tukey's test for multiple comparisons.
Figure 4A. hTR RNA expression, normalized to ACTB mRNA and to U2OS cells (n=3). Data information: Data are shown as mean ± SD, with the red line marking the mean value.B. Western blot for NHP2 expression in the early passages of hTR-expressing IMRB or U2OS cells, with β-Actin as loading control.C. Quantification of doubling time for U2OS, U2OS-hTR, IMRB and IMRB-hTR cell lines (n=3). Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution Data information: Data are shown as mean ± SD, with the red line marking the mean value.D. Western blot in U2OS and U2OS-hTR cells treated with or without Hydroxyurea (2 mM for 24 h), for the indicated proteins. β-Actin and Vinculin served as loading controls.E. Representative images for co-localization between p-RPA32S33 (green) and telomeres (red FISH probe). White arrows indicate co-localization events. Scale bars: 5 μm. Quantification is shown on the right. At least 50 nuclei were scored in each of the three independent experiments performed. Each dot represents the number of p-RPA32S33-telomere co-localizations in one nucleus. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution. Data information: Data are shown as mean ± SD, with the red line marking the mean value.F. Quantification of live cell numbers with trypan blue staining and manual counting, 72 h post-transfection with siLuc, sihTR#1 or sihTR#2, measured as the percentage of siLuc-treated cells (n=3). On the right, representative images of U2OS cells treated with sihTR#1 for 72 h. Statistical analyses were performed using one-way ANOVA with post-hoc Tukey's test for multiple comparisons. Data information: Data are shown as mean ± SD, with the red line marking the mean value.G. Western blot analyses in the indicated cell lines cells transfected for 48 h with siLuc or sihTR, with the indicated antibodies. U2OS cells treated with Hydroxyurea (2 mM for 24 h) were used as positive controls.
Figure 6A. Representative images for co-localization between PML (green) and telomeres (red FISH probe). White arrows indicate co-localization events, marking APBs. Quantification of images is shown on the right. A total of 85-300 nuclei were scored from four independent experiments. Each dot represents the number of events per nucleus: PML foci or APBs. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution.C. Representative images for co-localization between EdU (green) and telomeres (red FISH probe) within PML bodies (cyan) in U2OS and U2OS-hTR cells. White arrows indicate co-localization events. D, E. Quantification of images shown in (C). Co-localization events between EdU foci and APBs were analyzed (D) and the percentage of APBs that are actively engaged in G2/M telomeric synthesis was calculated (E). A total of at least 50 nuclei were scored from three independent experiments. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution.
Figure 7A. Representative images for co-localization between RPA32 (green) and single stranded C-rich telomeric sequences (ss-Telo-C, red FISH probe), in G2/M-synchronized cells. White arrows indicate co-localization events. B. Quantification of images in (A). A total of at least 55 nuclei were scored from three independent experiments. Each dot represents the number of events per nucleus. C. Quantification of images in (A) as the percentage of single stranded C-rich telomeric sequences associated with RPA32 to total single stranded C-rich telomeric sequences A total of at least 55 nuclei were scored from three independent experiments. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise.D. Representative images for co-localization between RPA32 or p-RPA32S33 (green) and telomeres (red FISH probe), in the indicated cell lines treated with siLuc or sihnRNPA1 for 48 h. White arrows indicate co-localization events. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise.E. Quantification of images in (D). A total of 85-135 nuclei were scored from three independent experiments. Each dot represents the number of events per nucleus. Statistical analyses were performed using one-way ANOVA with post-hoc Games-Howell's test for multiple comparisons. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise.F. Representative images for co-localization between RPA32 (green) and telomeres (red FISH probe) in the indicated cell lines treated with either DMSO or 10 μM NU7026 DNA-PK-inhibitor (DNA-PK-i) for 48 h. White arrows indicate co-localization events. Quantifications are shown on the right. A total of 53-84 nuclei were scored from three independent experiments. Statistical analyses were performed using one-way ANOVA with post-hoc Games-Howell's test for multiple comparisons. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value. Statistical analyses were performed with Student's two-tailed unpaired t-test with Welch's correction for normality of distribution unless specified otherwise.
Figure 8A. hTR and NHP2 mRNA expression in U2OS-hTR cells 72 h post-transfection with siLuc or siNHP2, normalized to ACTB mRNA and to SI14/TEL+ cells (n=2). Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value.B. Western blot showing NHP2 protein levels in the indicated cells treated with siLuc or siNHP2, 72 h post-transfection. β-Actin served as loading control.C. Representative images for co-localization between p-RPA32S33 (green) and telomeres (red FISH probe) in U2OS-hTR cells treated with control siLuc or siNHP2, 72 h post-transfection. White arrows indicate co-localization events. Quantification of images is shown on the right. Each dot indicates the number of co-localization events per nucleus (n=2). Statistical analyses were performed with Student's two-tailed unpaired t-tests with Welch's correction for normality of distribution. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value.D. Representative images for TIF quantification: co-localization between 53BP1 (green) and telomeres (red FISH probe), measured 72 h post-transfection with siLuc or siNHP2 in cells further treated with either DMSO or Doxorubicin (DOXO, 1.5 μM) for 2 h. White arrows indicate TIF. Quantification of TIF is shown on the right. A total of at least 100 nuclei were scored from two independent experiments. Each dot represents the number of TIF per nucleus. Statistical analyses were performed using one-way ANOVA with post-hoc Games-Howell's test for multiple comparisons. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value.E. Western blot analyses in U2OS and U2OS-hTR cells treated with the indicated siRNAs, 72 h post-transfection to detect total protein levels using the indicated antibodies. β-Actin and Vinculin served as loading controls.F. Representative images for 53BP1 foci formation in U2OS and U2OS-hTR cells, 72 h post-transfection with siLuc or siNHP2, exposed to 1 Gy irradiation and fixed at various time points post-irradiation as indicated.G. Quantification of images shown in (F). Each dot represents the number of 53BP1 foci per nucleus. A total of 114-140 nuclei were scored from two independent experiments. Statistical analyses were performed using two-way ANOVA Mixed effect model. Data information: All scale bars: 5 μm. Data are shown as mean ± SD, with the red line marking the mean value.
Figure 1Left panel: Representative blots showing protein levels of YAP and phosphorylated YAP (Ser127) in control (si-Ctrl) and YAP knockdown (si-YAP) hOSEs. Protein levels were detected by Western blotting. β-Actin was used as a protein loading control. Right panel: Cell proliferation (cell number) in control (si-Ctrl) and YAP knockdown (si-YAP) hOSE cells.Left panel: Representative blots showing protein levels of YAP and phosphorylated YAP (Ser127) in control (si-Ctrl), YAP knockdown (si-YAP), TAZ knockdown (si-TAZ) and YAP/TAZ double knockdown (si-YAP/TAZ) hOSE Cells. β-Actin was used as a protein loading control. Right panel: Cell proliferation (cell number) of control (si-Ctrl), YAP knockdown (si-YAP), TAZ knockdown (si-TAZ) and YAP/TAZ double knockdown (si-YAP&si-TAZ) hOSE cells.Representative blots showing protein levels of YAP and phosphorylated YAP (Ser127) in control hOSEs (MXIV) and hOSEs expressing wild type YAP (YAP) or constitutively active YAP (YAPS127A). Cells were collected at passage 4. β-Actin was used as a protein loading control.Growth of hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells at passage 4. Each bar represents mean ± SEM of four independent samples. Growth curves of hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells at passage 7.Representative images showing the morphologic changes of hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells at passage 7. Scale bar: 50 µm.Representative images showing expression and location of YAP in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells at passage 7. YAP was visualized using an Alexa-488 (Green) conjugated secondary antibody. Nuclei were stained with DAPI. Scale bar: 50 µm.
Figure 2Representative images showing SA-β-gal staining in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells at passage 7. Scale bar: 50 µm.Quantitative data showing the ratio of SA-β-gal positive cells in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells.Representative blots showing expression of key proteins associated with cell cycle arrest in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells analyzed using Western blotting at the 7th passage. β-Actin was used as a protein loading control.Representative images showing expression and location of phosphorylated pRB (807/811) in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells at the 7th passage. Phosphorylated pRB (807/811) was visualized using an Alexa-488 (Green) conjugated secondary antibody. Nuclei were stained with DAPI (blue).
Figure 3Representative blots showing expression of TP53 and RB1 expression in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells in the presence or absence of HPV E6/E7 oncoproteins at the 7th passage. β-Actin was used as a protein loading control.Growth curves of hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells in the presence of HPV E6/E7 oncoproteins at the 9th passage. Each point represents mean ± SEM of four independent samples.Representative images showing SA-β-gal staining in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells in the presence or absence of HPV E6/E7 oncoproteins at the 7th passage.Quantitative data showing the ratio of SA-β-gal positive cells in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells in the presence or absence of HPV E6/E7 oncoproteins.Representative images showed SA-β-gal staining in cultured primary hOSE cells in the presence or absence of HPV E6/E7 oncoproteins. Cells were collected and stained at the 13th passage. Scale bar: 25 µm. Right graph showing the ratio of SA-β-gal positive cells in P13 hOSE cells transfected with empty vectors (CTRL-P13) or vectors expressing E6/E7 oncoproteins (E6/E7-P13).
Figure 4Representative blots showing the expression of YAP and RB1 in hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells with or without pRB knockdown. Cells were collected at the 7th passage and protein levels were analyzed using Western blotting. β-Actin was used as a protein loading control.Growth curves of hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells after RB1 knockdown.Representative images showing SA-β-gal staining in the hOSE-MXIV, hOSE-YAP and hOSE-YAPS127A cells with or without RB1 or TP53 knockdown. Cells were collected at the 7th passage.
Figure 5Representative blots showing LATS and YAP protein levels in hOSEs at the fourth, ninth and thirteenth passage. β-Actin was used as a protein loading control.Representative gel photographs showing mRNA levels of LATS1 and LATS2 in hOSEs with or without LATS2 knockout. CRISPR/Cas-9 system was used to knock out LATS2 in hOSEs. LATS1/2 mRNA levels were analyzed using qRT-PCR at the fourth, ninth and thirteenth passage.Representative images showing the expression and location of YAP in control and LATS2 knockout hOSE cells. Cells were collected at the 4th passage and YAP expression was examined by fluorescent immunohistochemistry. YAP protein was visualized using an Alexa-488 (Green) conjugated secondary antibody. Nuclei were stained with DAPI.Representative images showing YAP protein expression in the 13th passage hOSE Cells with or without LATS2 knockout. YAP protein was examined by fluorescent immunohistochemistry and visualized using an Alexa-488 (Green) conjugated secondary antibody. Nuclei were stained with DAPI (blue).Representative images showing morphological change (E) in the 13th passage (P13) hOSE cells with (LATS2 KO) or without (CTRL) LATS2 knockout.Representative images showing SA-β-gal staining (F) in the 13th passage (P13) hOSE cells with (LATS2 KO) or without (CTRL) LATS2 knockout. Scale bar: 50 µm.Quantitative data showing the ratio of SA-β-gal positive cells in the control (CTRL) and LATS2-knockout (LATS2KO) P13 hOSE cells.