Submitted:
21 August 2024
Posted:
22 August 2024
You are already at the latest version
Abstract
Keywords:
1. Introduction
2. Methods
2.1. Experimental Design
2.2. Cell Culture
2.3. Epigenetic Reprogramming
2.3.1. CRISPR-Cas9 Gene Editing
- ▪
- Target Genes: The genes p16INK4a and p53, associated with cellular senescence and aging, were selected for editing. Guide RNAs (gRNAs) specific to these genes were designed using online CRISPR design tools (e.g., CRISPR-ERA) and synthesized by a commercial provider.
- ▪
- Transfection: Fibroblasts were transfected with CRISPR-Cas9 plasmids containing the gRNAs using Lipofectamine 3000 (Thermo Fisher Scientific) according to the manufacturer's protocol. Successfully edited cells were selected using puromycin resistance and validated by Sanger sequencing.
2.3.2. Small Molecule Inhibitors
- ▪
- DNA Methylation: 5-Azacytidine (5-AzaC), a DNA methyltransferase inhibitor, was used to reduce DNA methylation levels. Cells were treated with 5 μM 5-AzaC for 72 hours.
- ▪
- Histone Deacetylation: Trichostatin A (TSA), a histone deacetylase inhibitor, was used to promote histone acetylation. Cells were treated with 100 nM TSA for 24 hours.
2.4. DNA Methylation Analysis
2.5. Transcriptomic Analysis
2.6. Functional Assays
2.6.1. Cellular Senescence
2.6.2. Proliferation Assay
2.6.3. Differentiation Assay
2.7. Statistical Analysis
3. Data Analysis
3.1. DNA Methylation Analysis
3.1.1. Data Processing
- ▪
- Raw bisulfite sequencing data were quality-checked using FastQC and trimmed for adapter sequences and low-quality bases using Trim Galore.
- ▪
- Cleaned reads were aligned to the human reference genome (GRCh38) using Bismark, a specialized bisulfite aligner.
- ▪
- Methylation calls were extracted using the Bismark methylation extractor and further processed using the methylKit package in R to identify differentially methylated regions (DMRs).
3.1.2. Statistical Analysis
- ▪
- Differential methylation analysis was performed comparing control and treated fibroblasts. Methylation differences at individual CpG sites and regions were statistically assessed using a logistic regression model.
- ▪
- A cutoff of p < 0.01 and an absolute methylation difference of >25% were applied to identify significant DMRs.
3.2. Transcriptomic Analysis
3.2.1. Data Processing
- ▪
- Raw RNA sequencing reads were quality-checked with FastQC and trimmed using Trim Galore to remove adapters and low-quality bases.
- ▪
- Trimmed reads were aligned to the human reference genome (GRCh38) using STAR aligner, and gene counts were quantified using featureCounts.
3.2.2. Differential Expression Analysis
- ▪
- Gene expression levels were normalized, and differential expression analysis was conducted using DESeq2 in R.
- ▪
- Genes with a false discovery rate (FDR) < 0.05 and absolute log2 fold change > 1 were considered significantly differentially expressed.
3.2.3. Functional Enrichment
- ▪
- Gene Ontology (GO) and KEGG pathway enrichment analyses were performed using the DAVID tool to identify biological processes and pathways significantly affected by epigenetic reprogramming.
- ▪
- Enrichment results were visualized using bar charts and network diagrams in Cytoscape.
3.3. Functional Assays
3.3.1. Cellular Senescence
- ▪
- The percentage of SA-β-gal positive cells was calculated by counting stained cells in randomly selected fields under a microscope.
- ▪
- Statistical significance was assessed using a two-tailed Student's t-test comparing treated and control groups.
3.3.2. Proliferation Assay
- ▪
- Absorbance readings from the CCK-8 assay were analyzed to determine cell proliferation rates. Data were normalized to the control group and statistically evaluated using one-way ANOVA followed by Tukey's post-hoc test.
3.3.3. Differentiation Assay
- ▪
- The efficiency of differentiation into adipogenic, osteogenic, and chondrogenic lineages was quantified by staining intensity measurements using ImageJ software.
- ▪
- Statistical comparisons were made using one-way ANOVA with post-hoc tests to determine significant differences between treated and control groups.
3.4. Statistical Tools and Software
3.4.1. R
3.4.2. GraphPad Prism
3.4.3. Cytoscape
3.4.4. ImageJ
3.5. Summary of Findings
3.5.1. DNA Methylation
3.5.2. Gene Expression
3.5.3. Functional Improvements
4. Results
4.1. DNA Methylation Changes
4.1.1. Global Methylation Levels
- ▪
- Treatment with 5-Azacytidine (5-AzaC) resulted in a significant reduction in global DNA methylation levels compared to the control group (p < 0.001).
- ▪
- Bisulfite sequencing data revealed a global decrease in methylation across the genome, with the most pronounced demethylation observed at CpG islands and gene promoters.
4.1.2. Differentially Methylated Regions (DMRs)
- ▪
- MethylKit analysis identified 2,345 DMRs between control and treated fibroblasts.
- ▪
- Among these, 1,530 regions showed hypomethylation, and 815 regions exhibited hypermethylation in the treated cells.
- ▪
- Significant DMRs were found in genes related to cellular senescence (e.g., CDKN2A), DNA repair (e.g., BRCA1), and metabolic processes (e.g., IGF1R), as shown in Figure 1.
4.2. Gene Expression Changes
4.2.1. Differential Expression Analysis
- ▪
- RNA sequencing (RNA-seq) revealed significant changes in gene expression profiles following epigenetic reprogramming.
- ▪
- DESeq2 identified 1,812 differentially expressed genes (DEGs) with an adjusted p-value < 0.05 and absolute log2 fold change > 1.
- ▪
- Of these, 1,045 genes were upregulated, and 767 genes were downregulated in treated cells.
4.2.2. Key Gene Changes
- ▪
- Upregulated genes included those involved in cell cycle regulation (e.g., CCNA2, CDK1), DNA repair (e.g., RAD51, XRCC5), and stem cell maintenance (e.g., SOX2, NANOG).
- ▪
- Downregulated genes were primarily associated with cellular senescence (e.g., CDKN2A, CDKN1A), inflammation (e.g., IL6, TNF), and extracellular matrix organization (e.g., COL1A1, COL3A1), as shown in Figure 2.
4.3. Functional Enrichment
- ▪
- Gene Ontology (GO) analysis revealed significant enrichment of biological processes related to cell proliferation, DNA repair, and stem cell differentiation.
- ▪
4.4. Functional Assays
4.4.1. Cellular Senescence
- ▪
- SA-β-gal staining indicated a 50% reduction in senescence-associated β-galactosidase positive cells in the treated group compared to controls (p < 0.01).
- ▪
- Microscopic images confirmed a notable decrease in the number of senescent cells after treatment, as shown in Figure 5.
4.4.2. Proliferation Assay
- ▪
- The CCK-8 assay demonstrated a significant increase in cell proliferation rates in the treated cells compared to control cells (p < 0.01).
- ▪
- Treated cells exhibited a higher absorbance at 450 nm, indicating enhanced metabolic activity and growth, as shown in Figure 6.
4.4.3. Differentiation Assay
- ▪
- Treated fibroblasts showed improved differentiation potential into adipogenic, osteogenic, and chondrogenic lineages.
- ▪
- Quantitative analysis of staining intensity revealed significant increases in Oil Red O, Alizarin Red S, and Alcian Blue staining in treated cells, indicating successful differentiation, as shown in Figure 7.
4.5. Summary of Findings
4.5.1. Epigenetic Reprogramming
- ▪
- Epigenetic interventions successfully induced global and gene-specific DNA methylation changes.
- ▪
- Significant alterations in gene expression profiles were observed, with increased expression of genes related to cellular rejuvenation and decreased expression of aging-related genes.
4.5.2. Functional Improvements
- ▪
- Treated fibroblasts exhibited reduced cellular senescence, enhanced proliferation rates, and improved differentiation potential, demonstrating the efficacy of epigenetic reprogramming in reversing aging phenotypes.
5. Discussion
Interpretation of Findings
5.2. Comparison with Previous Studies
5.3. Implications of the Study
5.4. Limitations of the Study
6. Future Directions
6.1. Optimization of Epigenetic Reprogramming Protocols
6.2. Evaluation in Animal Models
6.3. Clinical Translation
6.4. Long-Term Studies
6.5. Exploring Combination Therapies
7. Conclusion
Acknowledgments
Conflicts of Interest
References
- Bird, A. (2002). DNA methylation patterns and epigenetic memory. Genes & Development, 16(1), 6-21. [CrossRef]
- Horvath, S. (2013). DNA methylation age of human tissues and cell types. Genome Biology, 14(10), R115. [CrossRef]
- Issa, J. P. (2014). Aging and epigenetic drift: A vicious cycle. The Journal of Clinical Investigation, 124(1), 24-29. [CrossRef]
- Jones, P. A., & Takai, D. (2001). The role of DNA methylation in mammalian epigenetics. Science, 293(5532), 1068-1070. [CrossRef]
- López-Otín, C., Blasco, M. A., Partridge, L., Serrano, M., & Kroemer, G. (2013). The hallmarks of aging. Cell, 153(6), 1194-1217. [CrossRef]
- Rando, T. A., & Chang, H. Y. (2012). Aging rejuvenation and epigenetic reprogramming: Resetting the aging clock. Cell, 148(1-2), 46-57. [CrossRef]
- Sen, P., Shah, P. P., Nativio, R., & Berger, S. L. (2016). Epigenetic mechanisms of longevity and aging. Cell, 166(4), 822-839. [CrossRef]
- Esteller, M. (2008). Epigenetics in cancer. The New England Journal of Medicine, 358(11), 1148-1159. [CrossRef]
- Hargrove, P. W., & Hackett, P. B. (2014). Gene therapy using DNA-modifying enzymes. Molecular Therapy, 22(4), 702-707. [CrossRef]
- Blackburn, E. H., Epel, E. S., & Lin, J. (2015). Human telomere biology: A contributory and interactive factor in aging disease risks and protection. Science, 350(6265), 1193-1198. [CrossRef]
- Bocklandt, S., Lin, W., Sehl, M. E., Sánchez, F. J., Sinsheimer, J. S., & Horvath, S. (2011). Epigenetic predictor of age. PLoS ONE, 6(6), e14821. [CrossRef]
- Bolger, A. M., Lohse, M., & Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics, 30(15), 2114-2120. [CrossRef]
- Campisi, J. (2013). Aging cellular senescence and cancer. Annual Review of Physiology, 75, 685-705. [CrossRef]
- Cedar, H., & Bergman, Y. (2009). Linking DNA methylation and histone modification: Patterns and paradigms. Nature Reviews Genetics, 10(5), 295-304. [CrossRef]
- Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., ... & Gingeras, T. R. (2013). STAR: Ultrafast universal RNA-seq aligner. Bioinformatics, 29(1), 15-21. [CrossRef]
- Krueger, F., & Andrews, S. R. (2011). Bismark: A flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics, 27(11), 1571-1572. [CrossRef]
- Liao, Y., Smyth, G. K., & Shi, W. (2014). featureCounts: An efficient general-purpose program for assigning sequence reads to genomic features. Bioinformatics, 30(7), 923-930. [CrossRef]
- Love, M. I., Huber, W., & Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology, 15(12), 550. [CrossRef]
- Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., ... & Ideker, T. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Research, 13(11), 2498-2504. [CrossRef]
- Tollefsbol, T. O. (Ed.). (2010). Epigenetics of aging. Springer. [CrossRef]
- Zhang, Y., Rau, R. E., & Liu, Q. (2014). Epigenetic regulation of stem cell differentiation. Cell Stem Cell, 15(3), 301-315. [CrossRef]







Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).