Stem Cells International
Volume 2016 (2016), Article ID 2019498, 13 pages
http://dx.doi.org/10.1155/2016/2019498

Research Article

Senescence-Associated Molecular and Epigenetic Alterations in Mesenchymal Stem Cell Cultures from Amniotic Fluid of Normal and Fetus-Affected Pregnancy

Jūratė Savickienė, Sandra Baronaitė, Aistė Zentelytė, Gražina Treigytė, and Rūta Navakauskienė

Department of Molecular Cell Biology, Institute of Biochemistry, Vilnius University, LT-10257 Vilnius, Lithuania

Received 27 May 2016; Accepted 25 August 2016

Academic Editor: Gary E. Lyons

Copyright © 2016 Jūratė Savickienė et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Human amniotic-fluid-derived mesenchymal stem cells (AF-MSCs) are interesting for their multilineage differentiation potential and wide range of therapeutic applications due to the ease of culture expansion. However, MSCs undergo replicative senescence. So far, the molecular mechanisms that underlie fetal diseases and cell senescence are still poorly understood. Here, we analyzed senescence-associated morphologic, molecular, and epigenetic characteristics during propagation of MSCs derived from AF of normal and fetus-affected pregnancy. AF-MSCs cultures from both cell sources displayed quite similar morphology and expression of specific cell surface (CD44, CD90, and CD105) and stemness (Oct4, Nanog, Sox2, and Rex1) markers but had interindividual variability in proliferation capability and time to reach senescence. Within passages 4 and 8, senescent cultures exhibited typical morphological features, senescence-associated β-galactosidase activity, increased levels of p16, and decreased levels of miR-17 and miR-21 but showed differential expression of p21, p53, and ATM dependently on the onset of cell senescence. These differences correlated with changes in the level of chromatin modifiers (DNMT1 and HDAC1) and polycomb group proteins (EZH2, SUZ12, and BMI1) paralleling with changes in the expression of repressive histone marks (H3K9me3 and H3K27me3) and stemness markers (Oct4, Nanog, Sox2, and Rex1). Therefore epigenetic factors are important for AF-MSCs senescence process that may be related with individuality of donor or a fetus malignancy status.

1. Introduction

Human amniotic fluid (AF) contains multiple cell types, including about 1% of cells with mesenchymal stem cell characteristics [1]. These amniotic-fluid-derived stem cells (AF-MSCs) can be obtained from a small amount of AF of second trimester pregnancy through procedure of amniocentesis used for prenatal diagnosis to determine the fetus status and genetic diseases [2, 3]. AF possesses several advantages such as safe and easy sampling, little maternal and child trauma, and no ethical restrictions. To date, AF is considered as an attractive source of MSCs having great potential application in regenerative medicine because of their extensive self-renewal, multilineage differentiation capacities [4–6], low or negligible immunogenicity, and inability to form tumors after implantation [7, 8]. However, the extended stem cell cultivation through sequential passaging in vitro may lead to a state of replicative senescence. It has been suggested that human MSCs derived from different sources become senescent as indicated by their morphological changes, attenuated expression of specific surface markers, decreased proliferation and differentiation potential, telomere length shortening, and alterations in global gene expression, DNA methylation patterns, and miRNA profiles [9–11]. Clearly, cellular senescence is a complex process and molecular events associated with AF-MSCs in vitro expansion are thus far unknown.

Primary somatic cells can replicate in culture about 50 cumulative population doublings, after which the cultures stop dividing. This phenomenon is termed Hayflick’s limit and is known as replicative senescence [12]. MSCs have a limited lifespan in vitro as any normal, somatic cells and are also subject to the Hayflick limit. It was suggested that replicative senescence of MSCs is a continuous process starting from the first passage [11], leading to seriously affected proliferative and clonogenic potential [11, 13] and related to morphological and molecular characteristics. Proliferation arrest within 10–20 passages was observed in senescent MSCs undergoing morphological and phenotype alterations [11, 14, 15] together with changes in the global gene expression pattern at different passages and significant downregulation of genes involved in cell cycle, DNA replication, and DNA repair in late passages [11]. By the use of gel-based proteomic method [16], the variety of fluctuations was observed for chaperones, signaling, antioxidant, proteasome, cytoskeleton, and connective tissue proteins in human AF stem cell line CD117/2 culture along passages 5–7 and 25. Two-dimensional gel electrophoresis-based quantitative proteomics in bone marrow MSCs cultures at passages 3–7 demonstrated also differential protein expression profiles associated with the impairment of cytoskeleton remodeling and/or organization and the repair of damaged proteins [17]. There is also evidence that senescence involves DNA damage, oxidative stress, a decrease in multilineage differentiation potential, and telomerase activity [15, 18, 19].

Several biomarkers have been used for quantitative assessment of cell senescence, including senescence-associated β-galactosidase (SA-β gal) and the cyclin-dependent kinase inhibitors, p16INK4A and p21WAF1, involved in the control of growth arrest by two major tumor suppressor pathways, p16INK4A/pRb and p53/p21WAF1 [20, 21]. A potential marker of the senescent state of stem cells may serve decreased expression of pluripotency transcription factors, such as Nanog and Oct4 [22, 23].

Replicative senescence seems to be epigenetically controlled [24]. Recent studies indicate that senescence-associated DNA methylation changes are associated with repressive histone marks and with targets of the histone methyltransferase EZH2, a component of polycomb complex PRC2 [25–27]. However, epigenetic regulation mechanisms underlying morphologic and phenotypic changes in MSCs during culture expansion are still unclear.

It is known that there exists different proliferative capability of MSCs cultures isolated from different donors [11, 13], undergoing senescence and considerable property changes. In order to further explore the observation about growth disadvantages, the present study aimed to examine how donor individuality and fetus abnormalities may influence amniotic-fluid-derived MSCs state during culture expansion. Here we provide evidence that MSCs cultures from AF of normal gestation and with fetus abnormalities exhibited diversity in their proliferation and senescence. We describe senescence-associated molecular and epigenetic changes during MSCs cultivation related with the extent of MSCs characteristics.

2. Material and Methods

2.1. Amniotic-Fluid-Derived MSCs Isolation, Culture, and Expansion

AF samples (about five milliliters) were obtained by biopsy (amniocentesis) from mid-second-trimester (16–24 weeks) or third-trimester (28–34 weeks) pregnant woman who needed prenatal diagnosis, and no abnormalities were revealed by genetic analysis using protocols approved by the Ethics Committee of Biomedical Researches of Vilnius District, number 158200-123-428-122. Samples were maintained at room temperature for about 4 hours prior to isolation of amniotic cells using two-stage protocol [23]. The sample was centrifuged at 1,800 rpm for 20 min, the supernatant was removed, and the cell pellet was washed once in DMEM medium (Sigma-Aldrich Ltd.) without serum to remove blood and cell debris. After centrifugation, the cell pellet was resuspended in 5 mL of growth medium (AmnioMAX™-C100 basal medium with AmnioMAX™-C100 supplement (Gibco, Life Technologies, Grand Island, NY, USA)), containing 100 U/mL penicillin, and 100 μg/mL streptomycin (Gibco, Grand Island, NY, USA) and plated in a 25 cm² culture flask (TPP, Switzerland). Amniocytes were incubated for 10–15 days at 37°C in 5% CO₂, when first colonies appeared (first stage). For culturing MSC (second stage), nonadhering AF cells were collected from primary culture and further expanded in a new 25 cm² culture flask at 37°C in 5% CO₂. After the appearance of the cell colonies, the growing medium was changed every 3 days. When cells reached confluence at 80%, subculturing into higher passages was performed by trypsinization with 0.05% trypsin-EDTA (Gibco, Life Technologies, Grand Island, NY, USA) for 3 min. A morphologically homogeneous population of fibroblast-like cells was obtained after two rounds of subculture. Cell morphology was observed by phase-contrast microscope (Nicon Eclipse TS100). Cell population doubling levels are presented as the average values calculated as follows: cell culture time (days)/number of passages.

2.2. Flow Cytometry Analysis

For identification of the phenotype of AF-MSCs from passages 4-5, cells were collected by centrifugation at 1,200 rpm for 6 min, washed once in phosphate buffered saline (PBS) with 0.2% fetal calf serum (FCS), and centrifuged again. A total of 5 × 10⁵ cells were resuspended in 50 μL of PBS with 1% BSA and incubated with fluorescein isothiocyanate- (FITC-) conjugated mouse anti-human antibodies against CD44 (Invitrogen), CD34 (Miltenyi Biotech), CD90 (Molecular Probes, Life Technologies), or phycoerythrin- (PE-) labeled CD105 (Invitrogen) and appropriate isotype control—mouse IgG2A-FITC (Miltenyi Biotech) or IgG1-PE (Molecular Probes, Life Technologies). Samples were incubated in the dark at 4°C for 30 min and finally analyzed with the Millipore Guava® easyCyte 8HT flow cytometer, using the InCyte 2.2.2 software. Ten thousand events were collected for each sample.

2.3. Senescence-Associated β-Galactosidase Assay

Cellular senescence was assessed using the β-galactosidase staining kit (Cell Signaling) following the manufacturer’s instruction. Cells were seeded in 48-well plates and cultivated for 48 hours. After washing with PBS, cells were fixed with 4% formaldehyde in PBS for 15 min and stained for β-galactosidase activity using staining solution overnight at 37°C in a thermostat without CO₂. Stained cells were viewed under a phase-contrast microscope (Nicon Eclipse TS100) and the data was expressed as the percentage of β-galactosidase positive cells.

2.4. RNA Isolation and RT-qPCR

MSCs from AF samples of normal gestation and with fetus abnormalities were cultured for 3–8 passages and analyzed for stem cell specific markers. Total RNA was extracted by TRIzol (Applied Biosystems, USA), as recommended by the manufacturer, and then reverse-transcribed into cDNA using Maxima First Strand cDNA Synthesis Kit (Thermo Scientific). qPCR was performed with Maxima® SYBR Green qPCR Master Mix (Thermo Scientific) on the Rotor-Gene 6000 system (Corbett Life Science). The amount of mRNA was normalized to GAPDH. The relative gene expression was calculated by a comparative threshold cycle delta-delta Ct method. Statistical analysis was performed using Student’s -test.

Forward (F) and reverse (R) primers (5′-3′) used in RT-qPCR are as follows:

 OCT4: F: CGAGAAGGATGTGGTCCGAG; R: CAGAGGAAAGGACACTGGTC Nanog: F: AGATGCCTCACACGGAGACT; R: GTTTGCCTTTGGGACTGGTG Sox2: F: TGGACAGTTACGCGCACAT; R: CGAGTAGGACATGCTGTAGGT Rex1: F: GCCTTATGTGATGGCTATGTGT; R: ACCCCTTATGACGCATTCTATGT p16: F: GCTGCCCAACGCACCGAATA; R: ACCACCAGCGTGTCCAGGAA p21: F: GGCAGACCAGCATGACAGATT; R: GCGGATTAGGGCTTCCTCT p53: F: TAACAGTTCCTGCATGGGCGGC; R: AGGACAGGCACAAACACGCACC ATM: F: CTCTGAGTGGCAGCTGGAAGA; R: TTTAGGCTGGGATTGTTCGCT GAPDH: F: TCCATGACAACTTTGGTATCG; R: TGTAGCCAAATTCGTTGTCA

The amount of mRNA was normalized to GAPDH. Normalization and fold changes were calculated using a comparative threshold cycle delta-delta Ct method. Student’s test was used for statistical analysis. All values are expressed as mean ± SD. value ≤ 0.05 was considered significant.

2.5. MicroRNA Expression

The expression of target miRNAs in AF-MSC samples was evaluated by RT-qPCR analysis using probe-based TaqMan MicroRNA Assays. Reverse transcription was performed using TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, USA); amplification was performed using TaqMan Universal PCR Master Mix (Applied Biosystems, USA) according to the manufacturer’s recommendations. Reactions were conducted in Finnzymes Thermocycler with a program as follows: 16°C for 30 min, 42°C for 30 min, and 85°C for 5 min. The product cDNA of each reaction was used as template for qPCR. RT-qPCR was carried out in a Rotor-Gene 6000 system (Corbett Life Science). The reaction was performed at 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 60 s. RNU 48 gene was used as a reference gene to normalize all experimental data. All reactions were performed in triplicate. The threshold cycle (Ct) was determined using the default threshold settings, and relative quantification of miRNA was calculated with the method. Student’s -test was used for statistical analysis. All values are expressed as mean ± SD. values ≤ 0.05 were considered significant.

2.6. Preparation of Proteins and Western Blotting

AF-MSCs (about 5 × 10⁵) were harvested by centrifugation (500 ×g, 6 min) after trypsinization with 0.05% trypsin-EDTA, washed twice in ice cold PBS, and resuspended in 10 volumes of lysis solution (62.5 mM Tris, pH 6.8, 100 mM DTT and 2% SDS, 10% glycerol). Benzonase (Pure Grade, Merck) was added to give a final concentration of 2.5 units/mL. Cell lysate was prepared by homogenization through the needle Nr 21 on ice and then centrifuged at 20,000 ×g for 10 min, 4°C. The supernatants were immediately subjected to electrophoresis or frozen at −76°C. Protein concentrations were measured using commercial RCDC protein assay (Bio Rad). The lysates were separated on a 7–15% polyacrylamide gradient SDS-PAGE gel and then transferred to a PVDF membrane. The filters were incubated with the primary antibody according to the manufacturer’s recommendations and then with horseradish peroxidase-conjugated (HPR) secondary antibody (Dako Cytomation, Glostrup) at room temperature for 1 h. The bands were developed using enhanced chemiluminescence detection (Amersham Pharmacia) according to manufacturer’s instruction. For Western blotting, the following antibodies were used: DNMT1 and EZH2 were from Cell Signaling; H3K9me3 was from Upstate; H3K27me3 and BMI1 were from Millipore; ATM, ATM (phospho S1981), and GAPDH were from Abcam; goat anti-rabbit or rabbit anti-goat HPR (horseradish peroxidase) linked secondary antibodies were from Dako Cytomation A/S.

3. Results

3.1. Characterization of MSCs Cultures from AF of Normal and Fetus-Affected Pregnancy during Passaging

MSCs were independently isolated from AF of healthy (N) donors (D1) and those carrying fetus abnormalities (P), such as D2, fetal central nervous system pathology and dilated brain ventricles; D3, circulatory disorders and increased heart; D4, trisomy 21 (Down’s syndrome). MSCs were cultured for 3–9 consecutive passages using AmnioMAX™-C100 complete medium. At each passage, MSCs were replated at 80–90% confluence. At early passages (p3-4), cell cultures from both AF sources formed and maintained homogeneous populations of typical elongated mesenchymal-type and spindle-shaped morphology (Figure 1(a)). The proliferation of MSCs achieved maximum on passage 3 (Figure 1(b)). After passages 5-6, there was a decline in the efficiency of proliferation in all cell populations. Time to reach senescence differed between cultures. As shown in Figure 1(b), some donors samples, including AF of normal and fetus-pathological gestation (D1, D2), stopped proliferation earlier (at p5-6), while others, slower proliferating cultures of fetus-pathological gestation (D3, D4), took longer time to achieve senescence (at p7-8). Morphological changes of cells were observed during passaging, where AF-MSCs cultures from different donors presented a typical MSCs morphology with the appearance of a varied proportion of flattened cells (Figure 1(a)). The comparison of images of cell cultures during passaging showed enlarged and flattened morphology cells with increased cytoplasmic granularity and frequency for positive staining with senescence-associated β-galactosidase (SA-β-gal) at late passages (Figure 1(c)). These senescence-associated changes occurred later in AF-MSCs samples (D3, D4) that showed relatively lower level of proliferation.

(a)

(b)

(c)

Figure 1: Characteristics of MSCs derived from AF of normal and pathological gestation during cultivation. (a) Representative images of initially spindle-shaped MSCs morphology in early (p3) and senescent passages (p7-8) with a 40x magnification are presented. (b) Differential proliferation of MSCs from AF of individual donors (D) during cultivation to passage 8. (c) Senescence-associated β-galactosidase staining in the late passage (p7-8); NS: nonstained cells. Representative results of four independent MSC preparations are demonstrated.

Concomitantly, the analysis of flow cytometry was performed to determine immunophenotype representing one of the major parameters for the characterization of MSCs cultures. As shown in Figure 2(a), MSCs from AF samples of normal gestation at early passages were characteristic in their immunophenotype, showing over 90% of cells positive for CD44 and CD90 and over 75% for CD105 but were negative for the hematopoietic marker CD34. The expression of those markers decreased at late senescent passages (by about 27–15%, resp.). Fetus-pathological samples (Figure 2(b)) at passage 3 maintained quite similar but variable level of cell surface markers CD44 and CD105 with a lower extent of mesenchymal marker CD90. At senescent passage 8, the expression of all phenotypic markers decreased about 25%, while the expression of CD105, which is involved in cell proliferation, dropped by about 40%.

(a)

(b)

(c)

Figure 2: Immunophenotype and stemness markers of MSCs derived from AF of normal and pathological gestation during cultivation. (a) Representative flow cytometry histograms of cell surface antigen markers of MSCs obtained from AF of normal gestation (N) at passage 3 and the late passage 8. Cells were positive for staining with CD44, CD90, and CD105 antibodies and negative for CD34 antibody. The appropriate IgG isotype was used as a control. (b) The percentage mean of CD44, CD90, and CD105 positive cells in MSCs cultures derived from AF of pathological gestation (P) at passages 3 and 8. (c) The expression of Nanog, Oct4, Sox2, and Rex1 from AF samples of normal (N) and pathological gestation (P) at passages 3, 6, and 8 was analyzed by RT-qPCR and normalized to GAPDH expression levels. Results are presented as the mean ± SD (, from each group). and were considered as significant changes.

To confirm stem cell origin of AF-MSCs, we performed RT-qPCR analysis of the main transcription factors from AF samples of healthy donors () and with fetal abnormalities (). Oct-4 and Sox2, involved in pluripotency or self-renewal, and Rex1, critically important in maintaining proliferative state in MSCs, were expressed at quite similar level in all AF samples at early and late passages (Figure 2(c)) while the expression of Nanog was relatively lower in fetus-pathological samples at passage 3. mRNA expression of stemness markers revealed alterations in the levels during passaging with variability in Oct4 expression in fetus-pathological samples. At late passage (p8), the increased expressional levels of Nanog and Rex1 were found in those samples as compared with passage 6, although the differences were not statistically significant.

3.2. Senescence-Associated Molecular Changes in MSCs Cultures

As is shown, during passaging of cell cultures derived from two AF sources, MSCs maintained their characteristic features such as immunophenotype and stemness factors, but the proliferative potential was seriously affected at late passages. Based on the observation that MSCs cultures differently undergo senescence in the course of cultivation, group I represents faster proliferating cell culture that earlier became senescent (at p5-6) and group II includes slower proliferating and undergoing senescence cultures (at p8-9) (Figure 3(a)). We observed that MSCs senescence was accompanied by typical morphological changes and cells became enlarged and flattened with enhanced SA-β-gal activity (Figure 3(b)). The same amount (about 30%) of SA-β-gal-positive cells was detected in groups I and II at p6 and p9 (denoted as late passage), respectively. To gain insight into the molecular characteristics of MSCs senescence, we analyzed changes in mRNA expression profiles of classical senescence-associated markers, p16INK4A (p16), p21WAF1 (p21), p53, and ATM (Ataxia telangiectasia mutated protein kinase), in both groups of cell cultures at early (p3) and senescent late passages paralleled with the monitoring of cellular morphology.

(a)

(b)

(c)

Figure 3: Characterization of senescence-associated markers in MSCs cultures during cultivation. (a) Differential proliferation of MSCs in two groups of cell cultures (I and II) during cultivation to passages 3, 5, and 8 (, for each group). (b) Nonstained (NS) and SA-β-gal positive cells at passages 6 and 9 in two groups of MSCs cultures. (c) Differential mRNA expression of senescence-associated markers in MSCs cultures (I gr and II gr) at passage 3 and the late passage was determined by RT-qPCR. Normalization to GAPDH and fold expression difference compared with passage 3 calculated using a comparative threshold cycle delta-delta Ct method. The data is presented as the mean ± SD (, from each group). , , and were considered as significant changes.

As demonstrated mRNA expression analysis determined RT-qPCR, senescent AF-MSCs cultures showed increased p16 expression (Figure 3(c)) with more pronounced changes (about 14-fold versus p3) in group I as compared with group II (about 5.5-fold versus p3). The level of p16 expression was inversely correlated with the proliferation capability. Surprisingly, we observed the increase in p21 and p53 expression in faster senescent cells from group I and the decrease in the level of both in slower senescent cultures from group II at late passages compared with early passage cells. The analysis of the expression of ATM that plays a role in cell cycle delay after DNA damage showed a positive correlation between increased expression levels of ATM and p53/p21 in senescent cultures from group I and a negative one from group II. The data demonstrate that MSCs cultures differently underwent cycle arrest together with the decline of the efficiency of proliferation and the appearance of senescent morphology.

Next, we focused on the expression analysis of miR-21 and miR-17 using RT-qPCR to evaluate the role of those miRNAs in the fate of MSCs during cultivation. The data presented in Figure 4(a) shows significant downregulation of both miR-17 and miR-21 in senescent cell cultures from group I and group II as well. Here we observed difference in fold changes of miR-17 expression ranging up to 2.7- and 2-fold downregulation in group I and group II, respectively, and quite similar decrease in miR-21 expression in samples of both groups (2.3–2.5-fold, resp.).

(a)

(b)

Figure 4: Differential expression of miRNAs and stem cell markers in MSCs cultures during passaging. RT-qPCR analysis of miR-17 and miR-21 (a) and Nanog, Oct4, Sox2, and Rex1 (b) from cell cultures of two groups (I and II) at passage 3 and the late passage. Normalization to calibrator sample (a) or GAPDH (b) and fold expression difference as compared with passage 3 calculated using a comparative threshold cycle delta-delta Ct method. The data is presented as the mean ± SD (, from each group). and were considered as significant changes.

In order to evaluate a possible link between miRNAs and self-renewal regulating factors in senescent MSCs, we performed the comparative RT-qPCR analysis of essential transcription factors Oct4, Nanog, Sox2, and Rex1 in samples of MSCs cultures from group I and group II during culturing from early (p3) to late senescent passages. As shown in Figure 4(b), senescence of cells from group I caused a decrease in the expression of Oct4 and Nanog and nonsignificant changes in Sox2 and Rex1 levels. MSCs from group II showed a statistically significant decrease in expression levels of Oct4, Nanog, Sox2, and Rex1 at senescent passage compared with those at early passage (p3). Likewise, miR-17 and miR-21 may regulate the proliferation and senescence of MSCs through the effects on components of cell cycle machinery and possibly by the interaction with the transcription factors involved in cell proliferation and self-renewal.

3.3. Molecular and Epigenetic Alterations Associated with Senescence of MSCs from AF of Normal Gestation and with Fetal Abnormalities

Next, we analyzed the functioning of the epigenetic regulatory factors HDAC and DNMT during the cellular senescence process. As shown in Figure 5(a), the level of DNMT1, which is constitutively expressed in proliferating cells for the maintenance of preexisting DNA methylation, decreased in MSCs cultures from groups I and II with the onset of cell senescence process and became almost undetectable at passage 6 (in group I) and at passage 8 (in group II). Additionally, the ATM analysis at protein level showed MSCs culture state before entering a growth arrest typical of senescence. In MSCs cultures of group I, phosphorylated ATM protein (ATM-P) accumulated at passages 6 and 8 concomitantly with the increase in the level of p53, while nonphosphorylated ATM was found at passages 4 and 6 with the apparent decrease at late passage (p8) (Figure 5(b)), showing cell cycle arrest after DNA damage-induced formation of double-strand breaks. In MSCs cultures II, less remarkable expression of ATM-P and ATM, gradually decreasing during passaging and in parallel with a marked decrease in p53 expression, was evident similarly as was demonstrated in Figure 3(c).

(a)

(b)

(c)

Figure 5: Epigenetic changes during passaging of MSCs cultures derived from AF of normal and pathological gestation. Lysates from cell samples of healthy (N) donors (D) and with fetal abnormalities (P) were subjected to Western blot analysis to monitor the expression of proteins using the indicated antibodies. (a, b) Representative blots of proteins from cell cultures of two groups (I and II) during culturing at passages 4, 6, and 8 and (c) from cell cultures derived from AF of normal gestation (D1) and individual donors (D2–D4) carrying fetus abnormalities. The data is representative of at least two gels showing similar results.

The analysis of HDAC1 and PRC2 polycomb group proteins, EZH2 and SUZ12, in samples of groups I and II demonstrated the expression profile similar to that of DNMT1 and HDAC1 (Figure 5(a)). In faster senescent cells from group I, a mark of constitutive heterochromatin, H3K9me3, showed decreased expression levels during passaging by a contrast to unchanged levels of this mark in cells from group II (Figure 5(a)). As shown in Figure 5(c), a decrease in the levels of DNMT1 and HDAC1 at late passages occurred in all cell cultures derived from individual AF samples from donors (D) of normal (N) and fetus-pathological (P) gestation and represented the extent of replicative ability of individual culture. For example, in sample from donor D4 (carrying Down’s syndrome), where MSCs senescence process was delayed, DNMT1 and HDAC1 expression remained higher than in other samples at the same passage (p8). Similar expression changes in the levels of PRC2 proteins, EZH2 and SUZ12, and PRC1 component BMI1 were noticed in AF samples of the same donors, demonstrating their implication in MSCs senescence process (Figures 5(a) and 5(c)).

The inverse correlation between expression levels of histone H3K27me3 and EZH2, which specifically trimethylates histone H3 at lysine (K) 27, was found in senescent cells from fetus-affected AF samples, showing dynamical changes in chromatin structure during cell senescence. This mark accumulation occurred earlier and was maintained over passaging of faster senescent culture from sample D1. The results indicate that the epigenetic modifying enzymes, PRC2 and PRC1 complex proteins, and repressive histone modifications together with miRNAs cooperatively participate in senescence process of AF-MSCs cultures derived from healthy donors and individual donors with fetus malignancy.

4. Discussion

In our previous study, we demonstrated that AF-MSCs cultures from second and third trimester of normal pregnancy exhibited characteristic stem cell futures by the expression of specific cell surface markers and their ability of multipotential differentiation [28]. To date, little is known how AF-MSCs subpopulations derived from normal and fetus-affected donors change during culture expansion. As shown in this study, some MSCs cultures derived from AF of defected pregnancy during passaging displayed some variations in their phenotype (CD44+, CD90+, and CD105+) and proliferation potential. MSCs from different donors presented a typical elongated and spindle-shaped morphology at early passages to the appearance of various proportions of enlarged cells compatible with morphology of cells entering senescence. AF-MSCs samples showed interindividual variability in proliferation capability and time to reach senescence as was described in previous studies [29–31]. The data [24] support the notion of heterogeneity in MSCs cultures and also with regard to replicative senescence. There are evidences about exiting differences in the long-term proliferative capability of MSCs isolated from different donors and a negative correlation between donor age and MSCs proliferative capacity [13]. Additionally, proliferation rates in different AF samples highly varied from indefinitive exponential growth length to abruption by senescence without any explanation [13]. The observations about growth disadvantages due to fetus abnormalities are in line with several other studies, which indicated unusual growth failures related to fetal aneuploidy [32], aneuploid karyotypes, and gestation age as well [33] with lower growth rate in cultures of third trimester than those of first trimester [34].

Here we demonstrated that cultivation of different MSCs samples derived from AF of normal and pathological gestation leads to senescence and is related to morphological and molecular characteristics expressed differently in rather early passages. During cultivation of AF-MSC samples (from p3 to p8), the alterations in gene expression levels of senescence-associated markers (p16, p21, p53, and ATM) were correlated with the population proliferation and senescence abilities (according to the percentage of SA-β-gal-positive cells). In faster senescent MSCs cultures, p21, p53, and ATM expression notably increased but decreased in cultures with delayed senescence (Figure 3). However, both cultures contained elevated levels of senescence marker p16. Controversial data exists in senescent bone marrow MSCs regarding the expression of p21, p53, and p16 from the upregulation of all of them [35] or of only p16 [36] to a reduced expression of p53 [37]. As known, cellular senescence is mainly regulated by p53/p21 and p16/pRb pathways, where the p53/p21 pathway mediates the replicative senescence and plays important roles in DNA damage response, while the p16/pRb pathway mediates stress-induced and premature senescence [38]. Based on our results, we suggest the involvement of both pathways in faster senescent MSCs, when p53 is induced and stabilized by phosphorylation of upstream kinases, including ATM and Chk2 [39–41]. Subsequently p53 upregulates transcription of p21, which activates Rb through the inhibition of cyclin E/Cdk2 complex, while activated Rb inhibits the transcription of E2F target genes, such as cyclin A and PCNA, causing a long-term cell cycle arrest [42]. Another Cdk inhibitor p16, which is mediated through transcriptional activation by Ets transcription factors, activates Rb through inhibition of cyclin D/Cdk4 and 6 complexes and accumulates in senescent cells causing cell growth retardation and arrest at G1 [21, 43, 44]. In a case of AF-MSCs from group II, where senescence process was delayed, only p16 was upregulated, suggesting the involvement of p16/pRb pathway and supporting a notion that p16 and p21 play different roles in the initiation and maintenance of senescence cell cycle arrest [45].

It is not unexpected that DNMT controls cell senescence. The proliferative properties of MSCs are maintained by the expression of key pluripotency genes that downregulates cell cycle regulators, such as p16 and p21 [22], through a direct binding of Oct4 and Nanog to the promoter of DNMT1, which enhances its expression maintaining DNA methylation [23]. In our study, the reduction in the DNMT1 expression level during AF-MSCs senescence paralleled with downregulation of PRC2 complex proteins, EZH2 and SUZ12, and upregulation of p16 and p21 (Figures 3(c) and 5). The effective downregulation of EZH2 with greater effects in faster senescent AF-MSCs cultures (Figure 5) coincides with the notion that EZH2 is implicated in replicative senescence as well [46].

miRNAs are a class of small and noncoding RNAs of 18–25 nucleotides that generally degrade target gene expression at the posttranscriptional level [47]. Increasing evidence shows that miRNAs contribute to senescence-related changes in gene expression of many human cell types including MSCs [11, 48]. Recently, miR-17 (a member of miR-17-92 cluster) was demonstrated to be significantly downregulated in many aging model systems [48], suggesting its role as a novel biomarker of cellular aging. miR-17 has been shown to target genes involved in cell cycle control [49], including p21 [48, 50], through transcription activation by E2F and repression by p53, consequently activating the cyclinD1/Cdk4 complex [51, 52]. This is in agreement with our results, demonstrating a decrease in miR-17 expression (Figure 4), occurring together with increased levels of p21 and p53 at late passages in faster senescent cultures (Figure 3). In slower senescent cultures, miR-17 acts possibly through the regulation of E2F target genes involved in cell cycle control.

The role of miR-21 has been studied in many fields, including stem cell biology and aging [53]. In AF-MSCs, miR-21 has been shown to be expressed in high levels and is implicated in the regulation of proliferation potential and cell cycle arrest by direct targeting Sox2 and inhibiting its expression or reducing Oct4 and Nanog expression through an indirect mechanism [54]. Recent study [54] has shown that miR-21 induction in cycling, presenescence cells causes the decrease in the expression of Sox2, Nanog, and Oct4 and is associated with a decrease in cell proliferation rate and cell cycle arrest at G0/G1. In our study, miR-21 expression was significantly lowered in MSCs at senescent passages (p5–8) in coordination with decreased levels of Oct4, Nanog, and Sox2 (Figure 4). The suggestion that miR-21 is involved in promoting the senescence by targeting Sox2 is based on the data that Sox2 (but not Oct-4 and Nanog) was identified as a target of miR-21, being negatively regulated in spindle-shaped AF-MSCs [54]. In agreement with our data, the reduced expression of Oct-4 and Nanog with increased expression of p21 has previously been reported in MSCs at late passage. On the contrary, knock-down of p21 in late passages MSCs inhibited senescence and increased cell proliferation and expression of stemness markers [54]. Thus, downregulated expression of Oct4 and Nanog was pointed as a potential candidate marker of the senescent state in stem cells [22, 23]. We also observed decreased expression of Rex1 in slower senescent MSCs cultures in association with a markedly decreased level of Nanog, which is transcriptional activator of Rex1, sustaining its expression.

The causal factors that mediate senescence process of AF-MSCs might include successive changes in the epigenetic state. Polycomb group proteins, which form multimeric complexes PRC1 and PRC2, have been shown to be implicated in replicative senescence [25–27]. The involvement of miRNAs and active/inactive histone marks at the promoter regions of cell cycle regulating p21 and p16 genes by targeting histone methyltransferases EZH1 is known. In turn, the inhibition of EZH1 in senescent MSCs leads to the demethylation of H3K27 and the activation of p16 expression that is consistent with our results. In our study, we have demonstrated changes in repressive histone modifications, H3K9me3 (downregulated) and H3K27me3 (upregulated), during AF-MSCs culture expansion between passages 4 and 8. It is widely accepted that methylation of both marks is associated with gene silencing, but gene repression by H3K27me3 is related to facultative heterochromatin formation, whereas repression by H3K9me2/me3 is linked with constitutive (permanent) heterochromatin [55]. Methylation of H3K27 by EZH1 or EZH2 is associated with gene repression via the polycomb group proteins [46]. In proliferating cells, DNA is mostly found to be in a less condensed euchromatic state, allowing access by the transcription and DNA replication machinery. Conversely, senescent cells are characterized by the presence of densely packaged facultative heterochromatin, organized into structures named senescence-associated heterochromatin foci (SAHF) [56]. The data of recent study [24] indicate that DNA methylation changes during stem cells senescence are associated with repressive histone marks (including H3K9me3) rather than with bivalent modifications (H3K4me3 and H3K27me3). DNA methylation profiles revealed consistent senescence-induced hypermethylation in regions associated with H3K27me3 and H3K4me1/me3 marks, whereas hypomethylation is linked with chromatin containing H3K9me3. DNA hypermethylation was significantly enriched in genes that are either up- or downregulated at later passages [57]. A global loss of H3K9me3 and changes in heterochromatin architecture were detected in MSCs with futures of cell aging [58]. The causal role of HDAC1 in cellular senescence via changes in the structure of chromatin followed a decrease in HDAC1 levels and growth arrest was demonstrated also [59]. In this study we defined a significant reduction of HDAC1 linked with the extent of senescent process of AF-MSCs (Figure 5). Thus, within the context of chromatin, chromatin-modifying enzymes and repressive histone modifications are implicated in the defining of AF-MSCs fate via chromatin remodeling.

In conclusion, our study indicates that senescence process of MSCs during culture expansion and passaging includes alterations in the expression of cell cycle regulating genes, stemness transcription factors, miRNAs, and epigenetic modulators that have to be taken into account for AF-MSCs therapeutic application.

Abbreviations

AF-MSCs: | Amniotic-fluid-derived mesenchymal stem cells
BMI1: | B lymphoma Mo-MLV insertion region 1
EZH1 and 2: | Enhancer of Zeste histone H3K27 methyltransferase
GAPDH: | Glyceraldehyde-3-phosphate dehydrogenase
Oct-4: | Octamer-binding transcription factor 4
PRC: | Polycomb repressive complex
RT-qPCR: | Quantitative real-time polymerase chain reaction
Sox-2: | Sex determining region Y-box 2
Suz: | Suppressor of Zeste.

Ethical Approval

The study with the human amniotic fluid was approved by the Ethics Committee of Biomedical Researches of Vilnius District, no. 158200-123-428-122.

Competing Interests

No potential conflicts of interests were disclosed.

Authors’ Contributions

Jūratė Savickienė designed the study, analyzed experiments, and wrote the paper. Sandra Baronaitė did amniotic-fluid-derived stem cell isolation, characterization, immunophenotype analysis, and preparation for the further analysis. Aistė Zentelytė performed RT-qPCR and flow cytometry. Gražina Treigytė did Western blotting. Rūta Navakauskienė conceived experiments, contributed to the interpretation of data, and revised the manuscript.

Acknowledgments

This work was supported by the Research Council of Lithuania (Project MIP-57/2015).

References

L. Perin, S. Sedrakyan, S. Da Sacco, and R. De Filippo, “Characterization of human amniotic fluid stem cells and their pluripotential capability,” Methods in Cell Biology, vol. 86, pp. 85–99, 2008.

Q. You, L. Cai, J. Zheng, X. Tong, D. Zhang, and Y. Zhang, “Isolation of human mesenchymal stem cells from third-trimester amniotic fluid,” International Journal of Gynecology & Obstetrics, vol. 103, no. 2, pp. 149–152, 2008.

P. Dekoninck, J. Toelen, S. Zia et al., “Routine isolation and expansion late mid trimester amniotic fluid derived mesenchymal stem cells in a cohort of fetuses with congenital diaphragmatic hernia,” European Journal of Obstetrics Gynecology and Reproductive Biology, vol. 178, pp. 157–162, 2014.

N. Siegel, M. Rosner, M. Hanneder, A. Valli, and M. Hengstschläger, “Stem cells in amniotic fluid as new tools to study human genetic diseases,” Stem Cell Reviews, vol. 3, no. 4, pp. 256–264, 2007.

S. Da Sacco, R. E. De Filippo, and L. Perin, “Amniotic fluid as a source of pluripotent and multipotent stem cells for organ regeneration,” Current Opinion in Organ Transplantation, vol. 16, no. 1, pp. 101–105, 2011.

O. Trohatou, N. P. Anagnou, and M. G. Roubelakis, “Human amniotic fluid stem cells as an attractive tool for clinical applications,” Current Stem Cell Research and Therapy, vol. 8, no. 2, pp. 125–132, 2013.

P. De Coppi, G. Bartsch Jr., M. M. Siddiqui et al., “Isolation of amniotic stem cell lines with potential for therapy,” Nature Biotechnology, vol. 25, no. 1, pp. 100–106, 2007.

M. Di Trapani, G. Bassi, E. Fontana et al., “Immune regulatory properties of amniotic fluid stem cells vary according to gestational age,” Stem Cells and Development, vol. 24, no. 1, pp. 132–143, 2015.

M. M. Bonab, K. Alimoghaddam, F. Talebian, S. H. Ghaffari, A. Ghavamzadeh, and B. Nikbin, “Aging of mesenchymal stem cell in vitro,” BMC Cell Biology, vol. 7, article 14, 2006.

J. Campisi and F. D. di Fagagna, “Cellular senescence: when bad things happen to good cells,” Nature Reviews Molecular Cell Biology, vol. 8, no. 9, pp. 729–740, 2007.

W. Wagner, P. Horn, M. Castoldi et al., “Replicative senescence of mesenchymal stem cells: a continuous and organized process,” PLoS ONE, vol. 3, no. 5, p. e2213, 2008.

L. Hayflick, “The limited in vitro lifetime of human diploid cell strains,” Experimental Cell Research, vol. 37, no. 3, pp. 614–636, 1965.

W. Wagner, S. Bork, P. Horn et al., “Aging and replicative senescence have related effects on human stem and progenitor cells,” PLoS ONE, vol. 4, no. 6, Article ID e5846, 2009.

Y. Jiang, B. N. Jahagirdar, R. L. Reinhardt et al., “Pluripotency of mesenchymal stem cells derived from adult marrow,” Nature, vol. 418, no. 6893, pp. 41–49, 2002.

R. Izadpanah, C. Trygg, B. Patel et al., “Biologic properties of mesenchymal stem cells derived from bone marrow and adipose tissue,” Journal of Cellular Biochemistry, vol. 99, no. 5, pp. 1285–1297, 2006.

W.-Q. Chen, N. Siegel, L. Li, A. Pollak, M. Hengstschläger, and G. Lubec, “Variations of protein levels in human amniotic fluid stem cells CD117/2 over passages 5–25,” Journal of Proteome Research, vol. 8, no. 11, pp. 5285–5295, 2009.

A. Madeira, C. L. da Silva, F. dos Santos, E. Camafeita, J. M. S. Cabral, and I. Sá-Correia, “Human mesenchymal stem cell expression program upon extended ex-vivo cultivation, as revealed by 2-de-based quantitative proteomics,” PLoS ONE, vol. 7, no. 8, Article ID e43523, 2012.

V. Janzen, R. Forkert, H. E. Fleming et al., “Stem-cell ageing modified by the cyclin-dependent kinase inhibitor ,” Nature, vol. 443, no. 7110, pp. 421–426, 2006.

T. Kiyono, S. A. Foster, J. I. Koop, J. K. McDougall, D. A. Galloway, and A. J. Klingelhutz, “Both Rb/ inactivation and telomerase activity are required to immortalize human epithelial cells,” Nature, vol. 396, no. 6706, pp. 84–88, 1998.

C. J. Sherr and F. McCormick, “The RB and p53 pathways in cancer,” Cancer Cell, vol. 2, no. 2, pp. 103–112, 2002.

J. Gil and G. Peters, “Regulation of the INK4b-ARF-INK4a tumour suppressor locus: all for one or one for all,” Nature Reviews Molecular Cell Biology, vol. 7, no. 9, pp. 667–677, 2006.

J. Nichols, B. Zevnik, K. Anastassiadis et al., “Formation of pluripotent stem cells in the mammalian embryo depends on the POU transcription factor Oct4,” Cell, vol. 95, no. 3, pp. 379–391, 1998.

C.-C. Tsai, P.-F. Su, Y.-F. Huang, T.-L. Yew, and S.-C. Hung, “Oct4 and nanog directly regulate Dnmt1 to maintain self-renewal and undifferentiated state in mesenchymal stem cells,” Molecular Cell, vol. 47, no. 2, pp. 169–182, 2012.

A. Schellenberg, Q. Lin, H. Schüler et al., “Replicative senescence of mesenchymal stem cells causes DNA-methylation changes which correlate with repressive histone marks,” Aging, vol. 3, no. 9, pp. 873–888, 2011.

A. Noer, L. C. Lindeman, and P. Collas, “Histone H3 modifications associated with differentiation and long-term culture of mesenchymal adipose stem cells,” Stem Cells and Development, vol. 18, no. 5, pp. 725–736, 2009.

J.-W. Jung, S. Lee, M.-S. Seo et al., “Histone deacetylase controls adult stem cell aging by balancing the expression of polycomb genes and jumonji domain containing 3,” Cellular and Molecular Life Sciences, vol. 67, no. 7, pp. 1165–1176, 2010.

A.-Y. So, J.-W. Jung, S. Lee, H.-S. Kim, and K.-S. Kang, “DNA methyltransferase controls stem cell aging by regulating BMI1 and EZH2 through microRNAs,” PLoS ONE, vol. 6, no. 5, Article ID e19503, 2011.

J. Savickiene, G. Treigyte, S. Baronaite et al., “Human amniotic fluid mesenchymal stem cells from second- and third-trimester amniocentesis: differentiation potential, molecular signature, and proteome analysis,” Stem Cells International, vol. 2015, Article ID 319238, 15 pages, 2015.

W. Wagner, A. D. Ho, and M. Zenke, “Different facets of aging in human mesenchymal stem cells,” Tissue Engineering Part B: Reviews, vol. 16, no. 4, pp. 445–453, 2010.

J. Ren, D. F. Stroncek, Y. Zhao et al., “Intra-subject variability in human bone marrow stromal cell (BMSC) replicative senescence: molecular changes associated with BMSC senescence,” Stem Cell Research, vol. 11, no. 3, pp. 1060–1073, 2013.

D. C. Colter, R. Class, C. M. DiGirolamo, and D. J. Prockop, “Rapid expansion of recycling stem cells in cultures of plastic-adherent cells from human bone marrow,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 7, pp. 3213–3218, 2000.

C. T. Maguire, B. L. Demarest, J. T. Hill et al., “Genome-wide analysis reveals the unique stem cell identity of human amniocytes,” PLoS ONE, vol. 8, no. 1, Article ID e53372, 2013.

W. H. Persutte and R. R. Lenke, “Failure of amniotic-fluid-cell growth: is it related to fetal aneuploidy?” The Lancet, vol. 345, no. 8942, pp. 96–97, 1995.

D. Moschidou, K. Drews, A. Eddaoudi, J. Adjaye, P. de Coppi, and P. V. Guillot, “Molecular signature of human amniotic fluid stem cells during fetal development,” Current Stem Cell Research and Therapy, vol. 8, no. 1, pp. 73–81, 2013.

H. Cheng, L. Qiu, J. Ma et al., “Replicative senescence of human bone marrow and umbilical cord derived mesenchymal stem cells and their differentiation to adipocytes and osteoblasts,” Molecular Biology Reports, vol. 38, no. 8, pp. 5161–5168, 2011.

K. R. Shibata, T. Aoyama, Y. Shima et al., “Expression of the p16INK4A gene is associated closely with senescence of human mesenchymal stem cells and is potentially silenced by DNA methylation during in vitro expansion,” Stem Cells, vol. 25, no. 9, pp. 2371–2382, 2007.

J. Kim, J. W. Kang, J. H. Park et al., “Biological characterization of long-term cultured human mesenchymal stem cells,” Archives of Pharmacal Research, vol. 32, no. 1, pp. 117–126, 2009.

I. Ben-Porath and R. A. Weinberg, “The signals and pathways activating cellular senescence,” The International Journal of Biochemistry & Cell Biology, vol. 37, no. 5, pp. 961–976, 2005.

M. Serrano, A. W. Lin, M. E. McCurrach, D. Beach, and S. W. Lowe, “Oncogenic ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a,” Cell, vol. 88, no. 5, pp. 593–602, 1997.

K. Webley, J. A. Bond, C. J. Jones et al., “Posttranslational modifications of p53 in replicative senescence overlapping but distinct from those induced by DNA damage,” Molecular and Cellular Biology, vol. 20, no. 8, pp. 2803–2808, 2000.

F. di Fagagna, P. M. Reaper, L. Clay-Farrace et al., “A DNA damage checkpoint response in telomere-initiated senescence,” Nature, vol. 426, no. 6963, pp. 194–198, 2003.

N. Dyson, “The regulation of E2F by pRB-family proteins,” Genes & Development, vol. 12, no. 15, pp. 2245–2262, 1998.

E. Hara, R. Smith, D. Parry, H. Tahara, S. Stone, and G. Peters, “Regulation of p16CDKN2 expression and its implications for cell immortalization and senescence,” Molecular and Cellular Biology, vol. 16, no. 3, pp. 859–867, 1996.

D. S. Peeper, T. M. Upton, M. H. Ladha et al., “Ras signalling linked to the cell-cycle machinery by the retinoblastoma protein,” Nature, vol. 386, no. 6621, pp. 177–181, 1997.

G. H. Stein, L. F. Drullinger, A. Soulard, and V. Dulić, “Differential roles for cyclin-dependent kinase inhibitors p21 and p16 in the mechanisms of senescence and differentiation in human fibroblasts,” Molecular and Cellular Biology, vol. 19, no. 3, pp. 2109–2117, 1999.

R. Margueron, G. Li, K. Sarma et al., “Ezh1 and Ezh2 maintain repressive chromatin through different mechanisms,” Molecular Cell, vol. 32, no. 4, pp. 503–518, 2008.

D. P. Bartel, “MicroRNAs: target recognition and regulatory functions,” Cell, vol. 136, no. 2, pp. 215–233, 2009.

M. Hackl, S. Brunner, K. Fortschegger et al., “miR-17, miR-19b, miR-20a, and miR-106a are down-regulated in human aging,” Aging Cell, vol. 9, no. 2, pp. 291–296, 2010.

L. He, J. M. Thomson, M. T. Hemann et al., “A microRNA polycistron as a potential human oncogene,” Nature, vol. 435, no. 7043, pp. 828–833, 2005.

J. Grillari, M. Hackl, and R. Grillari-Voglauer, “miR-17-92 cluster: ups and downs in cancer and aging,” Biogerontology, vol. 11, no. 4, pp. 501–506, 2010.

K. Woods, J. M. Thomson, and S. M. Hammond, “Direct regulation of an oncogenic micro-RNA cluster by E2F transcription factors,” The Journal of Biological Chemistry, vol. 282, no. 4, pp. 2130–2134, 2007.

H.-L. Yan, G. Xue, Q. Mei et al., “Repression of the miR-17-92 cluster by p53 has an important function in hypoxia-induced apoptosis,” The EMBO Journal, vol. 28, no. 18, pp. 2719–2732, 2009.

A. M. Krichevsky and G. Gabriely, “miR-21: a small multi-faceted RNA,” Journal of Cellular and Molecular Medicine, vol. 13, no. 1, pp. 39–53, 2009.

O. Trohatou, D. Zagoura, V. Bitsika et al., “Sox2 suppression by miR-21 governs human mesenchymal stem cell properties,” Stem Cells Translational Medicine, vol. 3, no. 1, pp. 54–68, 2014.

Q. Gan, T. Yoshida, O. G. McDonald, and G. K. Owens, “Concise review: epigenetic mechanisms contribute to pluripotency and cell lineage determination of embryonic stem cells,” Stem Cells, vol. 25, no. 1, pp. 2–9, 2007.

M. Narita, S. Nũnez, E. Heard et al., “Rb-mediated heterochromatin formation and silencing of E2F target genes during cellular senescence,” Cell, vol. 113, no. 6, pp. 703–716, 2003.

S. Hänzelmann, F. Beier, E. G. Gusmao et al., “Replicative senescence is associated with nuclear reorganization and with dna methylation at specific transcription factor binding sites,” Clinical Epigenetics, vol. 7, article 19, 2015.

W. Zhang, J. Li, K. Suzuki et al., “A Werner syndrome stem cell model unveils heterochromatin alterations as a driver of human aging,” Science, vol. 348, no. 6239, pp. 1160–1163, 2015.

D. Willis-Martinez, H. W. Richards, N. A. Timchenko, and E. E. Medrano, “Role of HDAC1 in senescence, aging, and cancer,” Experimental Gerontology, vol. 45, no. 4, pp. 279–285, 2010.