Epigenetic Reprogramming and Inheritance of the Cellular Differentiation Status Following Transient Expression of a Nonfunctional DN Rb Mutant in Murine MSCs

1. Introduction

Rb1 is a transcriptional corepressor 1, and one of its primary functions is to regulate the cell cycle, preventing uncontrolled cell proliferation [1,2,3]. Rb1 function is modulated by phosphorylation mediated by cyclin-dependent kinases and can be achieved through the use of a combination of phosphoacceptor sites, cyclin/cdk docking sites, and different cyclin/cdk complexes [4]. pRb1 can exist in only three states: unphosphorylated, monophosphorylated, which is considered “hypophosphorylated” active pRb, and hyperphosphorylated inactive pRb, and each phosphorylation-dependent state has a unique cellular function [5]. When pRb is in its active hypophosphorylated form, it inhibits the progression of the cell cycle from the G1 phase to the S phase by binding to and inhibiting E2F1, 2, and 3a transcription factors that serve as activators of many proteins involved in the S phase [6]. This arrest of the cell cycle gives the cells time to begin their differentiation processes. In addition to its role in controlling the cell cycle, pRb is recognized as a relevant player in the differentiation of a wide range of cells. [7,8]. Calo et al. reported that pRb plays a key role in the fate choice of differentiating common mesenchymal precursors between osteoblasts and adipocytes [9]. Although adipogenesis is orchestrated by the interaction of multiple transcription factors and epigenetic modifiers, encompassing multiple signaling pathways, two main transcription factors, peroxisome proliferator-activated receptor γ2 (PPARγ2) and CCAAT/enhancer-binding protein alpha (C/EBPα), are essential to complete adipocyte differentiation process. It has been shown that the physical interaction between pRb1 and C/EBPα influences the transcriptional activity of C/EBPα [10]. At the molecular level, pRb has been shown to contain a pocket motif containing two functionally essential domains, A and B, which are sufficient for transcriptional repressor activity [11]. It has been shown that even a single amino acid substitution, Cys706Phe, in domain B of the Rb1 pocket can block Rb1 function, inhibiting repressor activity [12]. Interestingly, overexpression of PPARγΔ5, a naturally occurring truncated isoform of Pparγ, has been shown to modify the PPARγ-induced transcriptional network, significantly impairing the differentiation capacity of adipocyte progenitor cells by acting as a dominant-negative splice isoform that reduces PPARγ activity [13].

The ΔS/N dominant-negative mutant (DN) of Rb1 used in the generation of the stably transfected 10T1/2 cell line was designed to eliminate 6 amino acids adjacent to the C-terminus of Rb1 (Figure 1a). This mutant form loses the ability to bind the protein that interacts with Rb1 and is predominantly in an unphosphorylated state, while retaining the ability to interact nonspecifically [14]. In the presence of an exogenously expressed ΔS/N DN Rb1 mutant in the nucleoplasm, this mutant interferes predominantly with the function of endogenous Rb1 such that it is able to abolish adipocyte differentiation through epigenetic silencing of PPARγ2. This silencing status, as we have previously shown, becomes heritable, representing the epigenetic reprogramming event of a key regulator of adipogenesis whose promoter is resistant to CpG demethylation in differentiated ΔS/N cells and is unable to express Pparγ2 [15]. Given that Cebpa and Pparγ2 are considered key regulators of terminal adipocyte differentiation and that Pparγ2 is constitutively suppressed by DNA methylation, we examined the epigenetic status of the Cebpa 5′ flanking region and mRNA expression to investigate whether Cebpa gene silencing also occurs in ΔS/N cells and whether gene-silent status is inherited. The methylation status of the promoters of the Rb1, Enhancer of zeste homolog 2 (Ezh2), Ubiquitously Transcribed Tetratricopeptide Repeat on chromosome X (Utx), Mixed-lineage leukemia 4 (Mll4) and Ten-eleven translocation 2 (Tet2) DNA dioxygenase genes associated with adipogenesis was also assessed to clarify the directional changes in the Pparγ2 and Cebpa genes in ΔS/N cells.

We found that the transient presence of a ΔS/N DN Rb1 mutant is able to elicit permanent silencing, in addition to Pparγ2, the proximal Cebpa promoter, thereby suppressing two major transcription factors responsible for the terminal phase of adipocyte differentiation and not affecting the methylation status of the Rb1, Ezh2, Utx, Mll4 and Tet2 genes involved in adipogenesis. Pparγ2 and Cebpa gene silencing status is heritable, as confirmed by recultivating cells over an extended period of time after they have lost ΔS/N DN Rb expression and rechecking their promoter methylation status. This is the first experimental evidence of apparently “targeted” epigenetic reprogramming through transient expression of a dominant-negative Rb1 mutant.

2. Results

2.1. Alterations in Rb1 Expression in Cells That Temporarily Expressed ΔS/N Rb1

2.1.1. ΔS/N Cells Overexpress Rb1 in the Undifferentiated State

To further study the consequences of transient expression of a nonfunctional mutant, ΔS/N Rb1, which has a deletion of 6 amino acids, structurally representing an alpha helix of the B domain Rb1 pocket motif (Figure 1a,b), in 10T1/2 cells in the presence of functional endogenous Rb1, we compared the expression levels of endogenous Rb1 and the member Rb family p130 protein in differentiated and undifferentiated (D, UD) ΔS/N cells using real-time RT‒qPCR analysis. The primer sequences are shown in Table 1. As shown in Figure 1d, moderate induction of Rb1 was observed in UD ΔS/N cells compared to wildt-type one. It has been shown that. Rb1 promoter activity is directly stimulated by its own gene product through the ATF-2 binding site [16]. Analysis of the ATF-2 binding site in the Rb1 promoter revealed that it contains a CpG, but this dinucleotide is not methylated in ΔS/N cells (Figure 1e,f); therefore Rb1 transcription is not blocked. Apparently, the ability of Rb1 to autoinduce its own expression was used to avoid the insufficient background of endogenous Rb1 due to the existence of the chimeric Rb1:Rb1 mutant in coassembled complexes during transient ΔS/N Rb1 expression. Although the Rb1 level was moderately elevated in UD ΔS/N cells, in response to adipogenic stimuli, it was almost twofold lower than the Rb1 level in wild-type D cells (Figure 1d). How the overexpression status of endogenous Rb1 is maintained in ΔS/N Rb1 postexpression culture remains an open question. A thorough investigation of the altered epigenetic status of the Rb1 gene is necessary.

2.1.2. Rb1 and p130 Show Weak Induction during Adipogenic Differentiation in ΔS/N Cells

Considering that RB1/Rb1 and RB2/p130 are indispensable in the second step of adipocyte differentiation and that the absence of RB1/Rb1 or RB2/P130 results in dysregulated adipose cells [17], we assessed the level of p130 expression in DN Rb1 mutant post-expressing cells. The p130 expression profiles in ΔS/N UD and D cells detected by RT‒qPCR are similar to those of Rb1 (Figure 1c,d) in terms of the response to adipogenic differentiation, highlighting the failure of both to be induced in response to adipogenesis stimulation, as occurs in wild-type cells. It appears that the changes, if any, affected the expression of both genes similarly.

2.2. Dnmt3a Is Upregulated in Undifferentiated ΔS/N Cells

DNA methylation/hydroxylation regulatory enzymes are key modulators of gene transcription regulation, acting as repressors for genes unsuitable for a particular type of differentiation and activating those required by the demands of cellular differentiation, as evidenced by numerous studies of MSC differentiation into specific cell lineages. To evaluate changes in Dnmt1, 3a, and 3b expression at the transcriptional level, we quantified mRNA transcripts in ΔS/N cells and their normal counterparts under differentiated and undifferentiated conditions (Figure 2a). We observed significant induction of Dnmt3a under adipogenic differentiation (AD) in wild-type and ΔS/N cells, which is consistent with a role for Dnmt3a in adipogenesis, but we also observed some upregulation of Dnmt3a in UD ΔS/N cells compared with wild-type cells (Figure 2a). It has been suggested that Dnmt3a-mediated DNA methylation is required for adipose tissue development because cells with heterozygous null mutations in Dnmt3a (lacking one allele) begin differentiating with lower DNA methylation and fail to experience the dynamic methylation alterations usually observed in WT cells [18]. However, ΔS/N cells, in which CpG methylation is still retained near the TSS of Cebpa (Figure 3b), fail to differentiate and do not express Cebpa. In contrast to Dnmt3a, which is strongly induced by AD stimuli in these cells (Figure 2a), Dnmt1 is moderately activated by AD in wild-type and ΔS/N cells, confirming the important role of Dnmt3a in adipocyte differentiation. Notably, Dnmt3b, a de novo DNA methyltransferase, is more inert to induction in response to adipogenic stimuli in ΔS/N cells (Figure 2a), highlighting its less prominent role in adipogenesis than Dnmt3a is. The induction of Dnmt1, 3a, 3b in wild-type D cells is consistent with the function of these enzymes in adipogenic differentiation, but misregulation of Dnmt3a in UD ΔS/N cells and Dnmt3b in D ΔS/N cells appears to be caused by altered epigenetics of these genes or their signaling pathways through transient expression of the ΔS/N DN Rb1 mutant.

2.3. Ezh2 Is Highly Overexpressed in Undifferentiated ΔS/N Cells

During adipocyte differentiation, Ezh2 dynamically regulates the expression of key genes involved in this process by controlling the chromatin structure and accessibility of their promoters. Ezh2-mediated silencing of negative regulators of adipogenesis and activation of proadipogenic genes are crucial steps in promoting adipocyte differentiation. As shown in Figure 2b, Ezh2 is highly upregulated in ΔS/N cells compared with wild-type cells, and propagation of the repressive H3K27me3 mark appears to occur among target genes in UD ΔS/N cells, whereas Ezh2 expression is dramatically reduced upon stimulation of adipogenic differentiation in ΔS/N cells and does not change in wild-type cells upon AD (Figure 2b).

2.4. Moderate UTX Induction Occurs in Differentiated ΔS/N Cells

The histone demethylase UTX is known to regulate brown adipocyte-to-myocyte remodeling in mature brown adipocytes. Demethylation of the repressive H3K27me3 mark at the Prdm16 promoter by UTX leads to high Prdm16 expression. PRDM16 then recruits the DNA methyltransferase DNMT1 to the Myod1 promoter, causing Myod1 promoter hypermethylation and suppressing its expression [19]. In addition to its H3K27me2/3 demethylase-independent function, UTX has been shown to promote chromatin remodeling through the BRG1-containing SWI/SNF remodeling complex [20]. UTX binds directly to multiple genes that encode proteins that interact physically with pRB, including RB1 itself. Since we observed the greatest induction of Dnmt1 and moderate induction of Utx during the AD of ΔS/N cells (Figure 2b), we hypothesize that Dnmt1 may be recruited to the Cebpa promoter by Utx to silence the gene state, similar to the hypermethylation of the Myod1 promoter mentioned above.

2.5. Heavily Methylated Cebpa in Wild-Type and ΔS/N Cells Loses CpG Methylation under AD, Maintaining an Unerased CpG Pattern near the TSS in ΔS/N Cells

Previously, we reported that cells that overexpress a nonfunctional exogenous ΔS/N variant of Rb1 that is lost during cell culture exhibit significant demethylation of CpG in the distal part of the Pparγ2 promoter, but the methylation of CpG at position 60 of its proximal part is still retained [15]. We assessed the methylation status of the 5’-flanking region of Cebpa, which contains 39 CpG sites located within a 510 bp core promoter. The fragment encompassing the proximal promoter was amplified, cloned, and sequenced from either wild-type or ΔS/N cells under D and UD conditions (Figure 3a). We found that wild-type or ΔS/N cells had a densely methylated Cebpa promoter, showing a trend toward greater methylation in ΔS/N cells in the UD state (Figure 3b). Interestingly, as with Pparγ2, the Cebpa promoter in ΔS/N cells was generally unmethylated in D cells but retained a group of methylated CpGs downstream of the TSS (Figure 3b). Analysis of Cebpa expression revealed complete abrogation of its expression at the mRNA level under AD in ΔS/N cells, but strong induction of Cebpa expression was observed in wild-type cells (Figure 3c).

Notably, Cebpb, another TF involved in adipogenesis, could not be induced in ΔS/N cells. A transient increase in Cebpb expression is one of the first events observed in the process of adipocyte differentiation. Cebpb acts directly on the Pparγ2 promoter to activate its transcription, but the failure of Pparγ2 to initiate transcription may result in the abrogation of Cebpb induction, or other mechanisms are responsible for the downregulation of Cebpb induction in ΔS/N cells.

2.5.1. Unerased CpG Methylation near the TSS May Interfere with the Formation of the Overall Transcription Complex

We used the UCSC Genome Browser transcription factor binding search tool in mice (GRCm39/mm39) to analyze a 35-bp region containing an unerased methylated CpG located downstream of the TSS (Figure 3b). The results revealed that the Nrf1 transcription factor (TF), which is sensitive to CpG methylation, has a high binding score to the locus where unerased methylation is present (Figure 4a, b). The Nrf2-ECH homology 5-like subdomain of the acidic domain 1 of Nrf1 is known to be a major component that interacts with the basal transcription machinery, p300/CBP, to activate transcription [21]. In addition, the Cap’n’Collar (CNC) domain of Nrf1, a common domain shared by all CNC proteins and proteins in the CNC family, functions as a critical regulator of terminal differentiation in the blood [22], bone [23], and adipose tissue [24]. Thus, it cannot be ruled out that Cebpa transcription initiation is abolished due to an inability to assemble a general transcription factor complex because of the absence of Nrf1, which is unable to bind to the methylated site. Methylation-sensitive NRF1 was shown to bind only in the absence of DNA methylation, and binding was abolished by direct inhibition by CG methylation but not by indirect inhibition by MBD mproteins [25]. Thus, we conclude that double CG methylation abolished Nrf1 binding to this site. The lack of Cebpa expression under AD in ΔS/N cells (Figure 2c) led us to speculate that Cebpa silencing in ΔS/N cells is highly selective.

2.5.2. Unerased CpG Methylation Is Observed in the Cebpa Promoter Away from the TSS

Notably, some CpGs in the central region of the Cebpa promoter in wild-type cells presented a patterned, unerased distribution of methylated CpGs (Figure 3 b). Since the discovery of the phenomenon of mCpG-dependent TF binding, which was observed via large-scale analyses of gene expression profiles and DNA methylomes, it has been suggested that some TFs bind to methylated regulatory elements and activate gene expression [26,27]. CEBPα and CEBPβ were found to specifically bind to a methylated sequence [28]. Again, we analyzed the 35-bp region containing the methylated CpG using the UCSC Genome Browser transcription factor binding search tool. Having discovered a putative TF binding site of CEBPG that recognizes the same motif as Cebpb/a (Figure 5a, b), we propose that the unerased mCpG may represent a potential methylation-dependent motif for self-activation of Cebpa. In addition, other studies have suggested a role for members of the KLF family, including Klf5, Klf6 and Klf15, in promoting adipogenesis [29]. Ectopic expression of KLF2 in preadipocytes inhibits Pparγ2 transcription, indicating its suppressive role in adipogenesis. [30]. It remains to be explored whether these patterned CpGs are specific for CpG methylation readers that promote increased Cebpa expression in response to adipogenic stimuli or whether they represent an intermediate hydroxymethylation state toward full demethylation.

2.6. CpG Methylation at Position 60 of the Pparγ2 Promoter Still Occurs in Postexpressing DN Rb1 Mutant ΔS/N Cells

To confirm the existence of Pparγ2 promoter methylation, we assessed the methylation status of the Pparγ2 5’ flanking region, which contains 7 CpG sites located near the transcription start site. The Pparγ2 promoter was slightly methylated in UD wild-type cells but maintained a densely methylated promoter in ΔS/N cells postexpressing ΔS/N nonfunctional Rb1 (Figure 6). Although the Pparγ2 promoter showed pronounced CpG demethylation, constitutive CpG methylation was detected in almost all the plasmid clones analyzed at position 60 of the proximal Pparγ2 promoter in the presence of differentiation stimuli (Figure 6).

2.7. The Promoters of the Adipogenesis-Related Genes Ezh2, Utx, Mll4, Rb1, and Tet2 Are Unmethylated Regardless of the Differentiation State

To assess the effects of transient expression of the ΔS/N DN Rb1 mutant on promoter methylation and mRNA expression of the Ezh2, Utx, Mll4, Rb1, and Tet2 genes involved in adipogenesis, their promoter methylation and mRNA expression were assessed, with the exception of Mll4 expression. All selected genes are CGIs and are required for adipogenic differentiation. The corresponding promoter regions of the genes were amplified (Figure 3a) with bisulfite-specific primers. The primer sequences, number of CpGs analyzed, and amplified fragment sizes are given in Table 2, and a representative bisulfite sequencing chromatogram is shown in Supplementary Figs. S1, 2, 3, and 4.

The demethylation of the Cebpa promoter observed during the differentiation of 10T1/2 cells is mediated by the active oxygenase Tet2, which is known to be recruited to the CpG islands of Cebpa by CREB [31]. Extensive demethylation of Cebpa was detected in wild-type and ΔS/N cells, confirming that Tet 2 is expressed and functional. MLL4, a major H3K4me1/2 methyltransferase on mammalian enhancers, is required for the activation of cell type-specific enhancers during differentiation. During adipogenesis, MLL4 exhibits cell type- and differentiation stage-specific genomic binding and colocalizes with lineage-determining TFs and H3K4me1/2 on active enhancers [32]. Methylation analysis of 84 CpGs in the 696 bp promoter revealed a lack of CGI methylation of Mll4 in ΔS/N and wild-type cells regardless of their differentiation state (Supplementary Figure 4). Therefore, it can be assumed that at the epigenetic level, Mll4 regulation is not subject to changes in DNA methylation, and Mll4 may function in adipogenic regulation. Since no tendency for methylation was observed among the Ezh2, Utx, Mll4, Rb1, and Tet2 genes either before or after differentiation, the existing effect of the transiently expressed ΔS/N DN Rb1 mutant may not affect the expression of these genes, at least through the mechanism of promoter methylation.

3. Discussion

In this study, we describe the phenomenon of “targeted” epigenetic reprogramming resulting from transient expression of the nonfunctional ΔS/N Rb1 mutant and confirm the maintenance of the reprogrammed epigenetic state in postexpression culture. Unlike mutations that prevent protein translation or disrupt protein structure, thereby causing loss of function, other mutations have alternative mechanisms, such as a dominant negative effect when the mutant protein directly or indirectly blocks the normal biological function of the wild-type protein [33]. They can thus cause (>50%) loss of function, even though only half of the protein is mutated [34]. This phenomenon is often observed for proteins that are able to coassemble into a complex with wild-type subjunints, in which mutant subunits can interfere with assembly/function [35]. During the expression of the ΔS/N DN Rb1 mutant, until it is lost, we hypothesize that the following scenario may occur. Given that the deletion of 6 amino acids results in the dysfunction of exogenously expressed Rb1 as a transcriptional corepressor but renders Rb1 unphosphorylatable so that it is not inactivated by a natural mechanism and assuming that the binding activity of the protein is retained, we hypothesize that mutant Rb1 interferes with the function of wild-type Rb1 by assembling Rb1 multiprotein complexes that have a normal, mutant Rb1 chimera in which the mutant protein can disrapt the activity of this complex and cause a disproportionate loss of function. Some proteins, such as BRCA1, have been shown to have RB-binding sites in both the C- and N-termini [36,37]. According to the mechanism of competitive binding, if one of the Rb molecules is a DN mutant and the other is functional, this leads to the formation of a dysfunctional molecular complex and does not allow the hybrid complex to fulfill its natural role, reducing the efficiency of the processes carried out by Rb1 and the endogenous Rb1 context, resulting in an Rb1-null phenotype. Surprisingly, in the case of the ΔS/N Rb1 mutant, we observed the induction of Rb1 instead of a decrease in endogenous Rb1 (Figure 1d).

Considering that Rb interacts with more than 300 proteins associated with multiple metabolic/signaling pathways, the existence of many Rb-nonproducing complexes may be stressful for cells. The ability of Rb to autoinduce its own expression is used to induce Rb1 overproduction in ΔS/N cells, resulting in a lack of loss of function because of the existence of faulty Rb1-containing complexes. It has been shown that Rb1 promoter activity is directly stimulated by its own gene product through the ATF-2 binding site [17]. Analysis of the Rb1 promoter, as well as the Rb1 ATF-2 binding site, revealed that it does not contain CpG methylation in ΔS/N cells and that Rb1 transcription is not blocked (Figure 1e,f). Since Rb1 is involved in the differentiation and fate of MSCs [9], the alteration of endogenous Rb1pull caused by the formation of nonproductive complexes in which Rb1 participates dramatically affects the fate of the cell. Our results revealed that in the context of an insufficient background of functional Rb1, the expression of two major adipogenic factors, Cebpa and Pparγ2, is suppressed by CpG methylation (Figure 3b, Figure 6), whereas the methylation of the Rb1, Ezh2, Mll4, Utx and Tet2 genes, which are required for adipocyte differentiation, does not occur (Supplementary Figures 3,4 ). The mRNA expression of Ezh2, Utx and Rb1 was detected to a certain degree in wild-type and ΔS/N cells (Figure 2b, Figure 3d), indicating that these epigenetic regulators are not targets for epigenetic silencing in DN Rb1 ΔS/N cells. This finding is consistent with data from other researchers on the role of Rb1 in cell differentiation and fate [9].

The formation of H3K27me3 by EZH2 has been shown to be a powerful repressive mechanism to silence HIV-1 [38]. EZH2 was also shown to interact with DNMT1, DNMT3A, and DNMT3B in cancer cells and to result in hypermethylation of genes, leading to increased permanent silencing of target genes [39]. The recruitment of DNMT3A to the 5’LTR CpG of HIV-1 upon the inhibition of UTX/JMJD3 by GSK-J4 likely results in the induction of DNA methylation. Demethylation of the repressive mark H3K27me3 on the Prdm16 promoter by UTX has been shown to result in high Prdm16 expression. PRDM16 then recruits the DNA methyltransferase DNMT1 to the Myod1 promoter, causing Myod1 promoter hypermethylation and suppressing its expression [40]. Since we observed strong induction of Dnmt1 and a moderate increase in Utx during the adipogenic differentiation of ΔS/N cells (Figure 2a,b), we hypothesize that Dnmt1 may be recruited to the Cebpa promoter to prolong the gene silencing state, similar to the hypermethylation of the Myod1 promoter mentioned above. In addition to its H3K27me2/3 demethylase-independent function, UTX has been shown to promote chromatin remodeling through the BRG1-containing SWI/SNF remodeling complex. UTX binds directly to multiple genes that encode proteins that interact physically with pRB, including the RB gene network implicated in cell fate control [41]. Taken together, these findings suggest that altered expression of Rb1, Ezh2 and Dnmt3a represents a new epigenetically reprogrammed state of genes, the combined action of altered activities which results in the permanent silencing of Cebpa and Pparγ2 in ΔS/N cells. Interestingly, CEBPA is predicted to interact with the N-terminus of DNMT3a, reducing the accessibility of DNMT3a to DNA and thereby preventing promoter methylation at target genes, including those containing the CEBPA binding motif. This could be seen as a mechanism by which CEBPA helps avoid gene promoter methylation in addition to the enzymatic Tet DNA demethylation that can be assumed to occur in AD in wild-type 10T1/2 cells on the basis of successful cell differentiation and expression of Cebpa and Pparγ2.

Given that CpG methylation plays a repressive role through a variety of actions, from direct inhibition of TF binding and hromatin structure effects to repositioning of the nucleosome via the (CpG)3 methylation element [43], we assume that unerased methylation of the CpG + (CpG)2 element in close proximity to the TSS may also contribute to Cebpa silencing via nucleosome repositioning. Interestingly, based on nucleosome assembly experiments performed on various DNA substrates, researchers have shown that the long azidohexynyl group (ahyC, an analog of 5-mC) present at multiple CpG sites promotes nucleosome assembly. Furthermore, ahyC DNA nucleosomes can be efficiently repositioned by a Snf2H chromatin remodeler and display similar thermal stability [44]. RB1 uses the ATPase activity of Brm or BRG1 to change nucleosome structure. This occurs in cooperation with histone deacetylases and/or histone demethylases to produce tight nucleosome structures and facilitate the formation of closed chromatin structures and RB1-mediated repression [45].

Our results suggest that epigenetic reprogramming occurs through the resetting of epigenetic marks such as DNA methylation patterns and presumably histone modifications and nucleosome repositioning that control gene expression, allowing cells to change their differentiation status or potential fate following stress-related events owing to transient expression of the mutant form of Rb1. By nature, large-scale epigenetic reprogramming occurs early in embryonic development and during gametogenesis, as well as artificially when various in vitro reprogramming conditions are used [46,47]. The described findings may reveal new conditions of epigenetic reprogramming caused by the temporal exposure of MSCs to DN Rb1. The unique effect of DN Rb1 overexpression in ΔS/N cells is that it selectively downregulates two major adipogenic genes through promoter methylation, alters the expression of Rb1 and the chromatin modifiers Ezh2 and Dnmt3a (Figure 2a,b; Figure 3d), and to counteract the stress situation, the cells eliminate the inserted DNA fragment encoding the mutant form of Rb1 while maintaining the epigenetically reprogrammed state of ΔS/N cells. How the demand for cells to differentiate into adipocytes is disrupted under the influence of the ΔS/N DN Rb1 mutant in more detail and how new epigenetic memory is established in ΔS/N cells are questions that remain to be answered. Importantly, the exact consequences of a dominant-negative pRB mutant can vary depending on the specific mutations involved and the cell type in question. Furthermore, the context in which these mutations occur can also influence the overall effects on cell differentiation and proliferation.

4. Materials and Methods

4.1. Cell Culture

Mouse embryonic polypotent fibroblasts (C3H10T1/2 (10T1/2)) were obtained from the American Type Culture Collection (ATCC). The cells were expanded in Eagle’s basal medium (Thermo Fisher Scientific Cat. 21010-046) supplemented with 10% heat-inactivated fetal bovine serum, 2 mM L-glutamine, and 50 μg/mL gentamicin in a CO2 incubator (5% CO2 and 100% humidity). The ΔS/N-cell line that has experienced transient expression of the DN ΔS/N Rb1 mutant was defrosted and cultured in Eagle’s basal medium (Thermo Fisher Scientific Cat. 21010-046) supplemented with 10% heat-inactivated fetal bovine serum and 2 mM L-glutamine in a CO2 incubator (5% CO2 and 100% humidity).

4.2. Preparation and Electrophoresis of RNA in Denaturing Gels and Its Quantitative Evaluation

RNA was extracted from cultured cells with acidic water-saturated phenol. The RNA integrity was monitored, and the RNA amount was normalized by electrophoresis in a denaturing agarose gel supplemented with formaldehyde. Upon electrophoresis, RNA was visualized in the gel by luminescence of the 18S and 28S bands stained with ethidium bromide. The RNA amounts in different samples were normalized to the 28S band intensity in the corresponding lanes via the TotalLab Quant program product.

4.3. Real-Time PCR

Total RNA was isolated from 2 × 10⁶ cells with a GeneJET RNA kit (Thermo Fisher Scientific, USA). Purified RNA was treated with DNase I (Thermo Fisher Scientific, USA) (3 units/35 μL at 37°C for 1 h). The reaction was stopped by the addition of 3.5 μL of 50 mM EDTA for 10 min at 65°C. The RNA was quantified via a spectrophotometer at 260 nm. Its integrity was verified by the presence of 28S and 18S bands after electrophoresis on denaturing agarose gels with formaldehyde. For cDNA synthesis, 2 μg of RNA was mixed with 0.5 μg of oligo-dT18 primer. The volume was adjusted to 12.5 μL and incubated for 5 min. The reaction mixture was cooled, 4 μL of fivefoldfold buffer, 1 μL (200 units) RevertAid reverse transcriptase (Thermo Fisher Scientific, USA), 2 μL of 10 mM dNTP mixture, and 0.5 μL of RNase inhibitor (Thermo Fisher Scientific, USA) were added, and the mixture was incubated at 42°C for 1 h. The reaction was stopped by heating to 70°C for 10 min. Real-time PCR was performed on an Applied Biosystems 7300 real PCR system. The parameters were as follows: denaturation at 95°C for 5 min, melting at 95°C for 15 s, annealing at 58°C for 30 s, and synthesis at 72°C for 20 s for 40 cycles. The reaction mixture was composed of 8 μL of a 2.5-fold mixture containing deoxynucleoside triphosphates (dNTPs), PCR buffer, MgCl2, Taq DNA polymerase SYBR Green I and ROX (Syntol, Russia), 0.2 μL of forward and reverse primers (Beagle, Russia) (Table 1), and 0.2 μL of cDNA and water to adjust the volume to 20 μL. The actin gene was used as a loading control. Relative gene expression was calculated via the formula R = 2–ΔΔC(T). These experiments were repeated 3 times.

4.4. DNA Extraction and Bisulfite Treatment

Genomic DNA was obtained by overnight cell incubation in TES buffer containing 0.1% SDS and 100 µg/ml proteinase K at 55°C with subsequent phenol/chloroform extraction and isopropanol precipitation. Two µg of DNA in 50 µl of TE buffer was denatured for 15 min in 0.3 M NaOH at 37°C. The denatured DNA was mixed with 550 µl of freshly prepared solutions of 10 mM hydroquinone and 3 M sodium bisulfite at pH 5.0 and incubated under mineral oil at 50°C for 12 h. Bisulfite-treated DNA was desalted by isopropanol precipitation, desulfonated with 0.3 M NaOH for 5 min at room temperature, and precipitated with ethanol. Converted DNA was dissolved in 100 µl of water and stored at -20°C. Three microliters of precipitated DNA was used for each PCR. DNA from peripheral white blood cells (WBCs) was used as a control for normal differentiated quiescent cells.

4.5. Bisulfite Sequencing Analysis

Bisulfite-treated DNA was used to amplify the promoter regions of genes involved in adipogenesis using primers specific for bisulfite-treated DNA (Table 2). The PCR was carried out via the use of 2.5 units of homemade Taq polymerase in a final volume of 50 µl and the following cycling conditions: 5 min at 95°C, followed by 35 cycles (30 sec denaturation at 95°C, annealing for 30 sec at 62°C and elongation at 72°C for 1 min). The PCR products were gel purified and cloned using TOPO TA cloning kit. To prevent clonal amplification of sequences, the competent transformed cells were plated immediately after heat shock; excluding shaking bacteria for 1 h. The recombinant plasmids were isolated via the DM miniprep method [48] and sequenced in both directions using M13 forward and reverse primers (Invitrogen) and an ABI BigDye Terminator Cycle Sequencing Kit v3.1 (Thermo Fisher, Waltham, MA, USA) with a Gene Amp 9700 PCR System (Thermo Fisher, Carlsbad, CA, USA). The sequences were detected with an ABI 3130XL Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Eighteen to 20 clones for each PCR product were sequenced and analyzed. Promoter methylation analysis was performed by aligning sequenced clones with the gene promoter region, where cytosine to uracil was converted in an in silico experiment. The efficiency of the conversion of cytosine to uracil in the bisulfite reaction was estimated as the ratio of cytosine in a non-CpG context to the total number of cytosines in the region. The clones with an efficiency of cytosine conversion less than 98% were omitted from the analysis.

4.6. Primer Design

The primer-BLAST software tool was used to design new target-specific primers for RT‒PCR experiments. This tool is available at http://www.ncbi.nlm.nih.gov/tools/primer-blast. All amplicon primers were designed to encompass exon‒exon boundaries to avoid genomic DNA amplification. The specificity of the amplicons and primer pairs was checked in silico using BLAST (National Center for Biotechnology Information) alignment tools.

To design primer-to-bisulfite region-specific matches for bisulfite-modified DNA, the following rule was applied: since the DNA strands are no longer complementary after bisulfite treatment, an individual primer set will amplify only one strand of the target sequence. The first primer was designed to anneal to the converted target sequence. The second primer was designed to anneal to the extension product of the first primer, not the opposite template strand. The PrimerSuite program (PrimerSuite, PrimerDimer and PrimerPlex) was used to generate robust primers for PCR bisulfite sequencing analysis. The PrimeSuite, PrimerDimer, and PrimerPlex modules are all available online at www.primersuite.com, www.primer-dimer.com, and www.primer-plex.com, respectively.

4.7. Bioinformatics

The Basic Local Alignment Search Tool (BLAST) commonly used in bioinformatics has been applied to search for similarities and identify homologous sequences.

The UCSC Genome Browser transcription factor binding search tool in mice (GRCm39/mm39) https://genome.ucsc.edu was used to analyze TF binding in the local region of the Cebpa promoter.

4.8. Statistical Methods and Associated Software

The CFX96™ real-time PCR detection system associated software and algorithm were used in this study. All experiments were performed in triplicate.

5. Conclusion

By identifying unerased, in response to adipogenic stimuli, CpG methylation near the TSS in the Pparγ2 and Cebpa promoters, we postulate that epigenetic reprogramming through resetting of DNA methylation patterns and, presumably, histone modifications and nucleosomal repositioning of these genes allows cells to change their differentiation status or potential fate to counteract a stressful situation unsuitable for cell differentiation. The absence of CpG methylation in other genes required for adipogenesis, regardless of their differentiation status, suggests a highly selective process to sort out the genes most important for adipogenic differentiation for their subsequent permanent repression. A precise understanding of the molecular pathways underlying stress-associated changes in MSCs may shed light on the development of stress-associated conditions for the “targeted” cell fate reprogramming that we observed in ΔS/N DN mutant postexpressing cells.