- Research article
- Open access
- Published:
Hmga2 protein loss alters nuclear envelope and 3D chromatin structure
BMC Biology volume 20, Article number: 171 (2022)
Abstract
Background
The high-mobility group Hmga family of proteins are non-histone chromatin-interacting proteins which have been associated with a number of nuclear functions, including heterochromatin formation, replication, recombination, DNA repair, transcription, and formation of enhanceosomes. Due to its role based on dynamic interaction with chromatin, Hmga2 has a pathogenic role in diverse tumors and has been mainly studied in a cancer context; however, whether Hmga2 has similar physiological functions in normal cells remains less explored. Hmga2 was additionally shown to be required during the exit of embryonic stem cells (ESCs) from the ground state of pluripotency, to allow their transition into epiblast-like cells (EpiLCs), and here, we use that system to gain further understanding of normal Hmga2 function.
Results
We demonstrated that Hmga2 KO pluripotent stem cells fail to develop into EpiLCs. By using this experimental system, we studied the chromatin changes that take place upon the induction of EpiLCs and we observed that the loss of Hmga2 affects the histone mark H3K27me3, whose levels are higher in Hmga2 KO cells. Accordingly, a sustained expression of polycomb repressive complex 2 (PRC2), responsible for H3K27me3 deposition, was observed in KO cells. However, gene expression differences between differentiating wt vs Hmga2 KO cells did not show any significant enrichments of PRC2 targets. Similarly, endogenous Hmga2 association to chromatin in epiblast stem cells did not show any clear relationships with gene expression modification observed in Hmga2 KO. Hmga2 ChIP-seq confirmed that this protein preferentially binds to the chromatin regions associated with nuclear lamina. Starting from this observation, we demonstrated that nuclear lamina underwent severe alterations when Hmga2 KO or KD cells were induced to exit from the naïve state and this phenomenon is accompanied by a mislocalization of the heterochromatin mark H3K9me3 within the nucleus. As nuclear lamina (NL) is involved in the organization of 3D chromatin structure, we explored the possible effects of Hmga2 loss on this phenomenon. The analysis of Hi-C data in wt and Hmga2 KO cells allowed us to observe that inter-TAD (topologically associated domains) interactions in Hmga2 KO cells are different from those observed in wt cells. These differences clearly show a peculiar compartmentalization of inter-TAD interactions in chromatin regions associated or not to nuclear lamina.
Conclusions
Overall, our results indicate that Hmga2 interacts with heterochromatic lamin-associated domains, and highlight a role for Hmga2 in the crosstalk between chromatin and nuclear lamina, affecting the establishment of inter-TAD interactions.
Background
Hmga are non-histone chromatin proteins belonging to the high mobility group protein family [1]. Hmga proteins contain three different DNA binding domains, named AT-hooks domains, by which these proteins bind short AT-rich regions on DNA, and an acidic C-terminal tail [2]. The relevance of Hmga2 functions is demonstrated by the phenotypes of mouse loss of function mutants. Indeed, Hmga2 gene KO causes the “pigmy” phenotype, mainly consisting in reduced body weight, about 40% of the normal weight [3]. This severe phenotype is due to the reduction of the adipose tissue mass and of that of skeletal muscles [4]. Hmga1, a protein similar to Hmga2, probably partially compensates for the absence of Hmga2: indeed, the Hmga1/Hmga2 double KO mice have a dramatic phenotype with a 70% reduction of the body weight at birth, not compatible with pups survival [5]. In addition, Hmga1/Hmga2 double KO also causes a certain degree of embryonic lethality, thus confirming the crucial role of these proteins during development [5].
Hmga proteins appear to dynamically interact with chromatin [6] and a vast literature indicates that they could be involved in various nuclear functions, among these heterochromatin formation, replication, recombination, and DNA repair [7,8,9]. Hmga proteins compete with histone H1 for the binding to internucleosomal linker DNA [10] and they also influence transcription, contributing to the assembly of complexes at initiator sequences [11] and to the formation of enhanceosomes [12,13,14].
Hmga2 is expressed during embryogenesis; its levels are undetectable in most adult cells, while it is strongly upregulated in cancer tissues [15]. Many studies aimed at addressing the functions of Hmga2 have been done in cancer cells, mostly because of its possible pathogenic role in various types of tumors. Hmga2 has a physiologic role in ESCs where it is induced upon differentiation and whose suppression alters the exit of ESCs from the ground state of pluripotency, preventing their transition into the primed state of pluripotency, i.e., into EpiLCs [16]. Thus, the transition from ESCs into EpiLCs represent an ideal experimental system where to study the functions of Hmga2 in a normal environment.
In this study, we observed that, in epiblast stem cells (EpiSCs), endogenous Hmga2 is mostly associated with heterochromatin regions interacting with Lamins and we demonstrated that, in pluripotent stem cells (PSCs) induced to exit from the naïve state, the loss of Hmga2 impacts the deposition and distribution of histone modifications and the integrity of the nuclear envelope. These results led to demonstrate that Hmga2 is crucial to guarantee the proper 3D structure of the chromatin, maintaining the right interactions among TADs inside the nucleus and their relationship with Lamin-associated chromatin domains.
Results
The absence of Hmga2 alters H3K27me3 deposition upon the induction of the transition from undifferentiated pluripotent stem cells into EpiLCs
Considering the general thought that Hmga2 is a chromatin protein, we explored the effects of Hmga2 suppression (Additional file 1: Fig. S1A) on histone modification dynamics taking place upon the transition of PSCs into EpiLCs. As shown in Fig. 1A (and Additional file 1: Fig. S1B), the trimethylation of K27 of histone H3 (H3K27me3) is affected in Hmga2 KO cells [16] induced to develop into EpiLCs.
Indeed, while in basal conditions (undifferentiated cells) the levels of H3K27me3 were comparable in Hmga2 wt and KO cells, they were maintained higher in Hmga2 KO cells throughout the three days of EpiLC induction. This difference between wt and KO clones is clearly appreciable in time course experiments, where the expected decrease of H3K27me3 is less pronounced in Hmga2 KO cells than in the wt cells, while the amount of H3K27me2 and H3K27ac showed similar dynamics in both cells (Fig. 1B).
To understand whether the observed effect on the variations of the epigenetic state due to Hmga2 absence is a general mechanism to allow pluripotent stem cells to leave the naïve state, we induced the differentiation of wt and KO Hmga2 PSC clones through the formation serum-free embryoid bodies (SFEBs). This 3D-model mimics different aspects of cell differentiation during early mammalian embryogenesis and it is particularly useful for the interpretation of KO phenotypes, allowing to by-pass the limitations related to 2D-systems [17, 18]. In this experimental setting, the different behavior of H3K27me3 between wt and KO cells was clearly confirmed, with a sustained presence of this histone mark at 2, 3, and 4 days after the induction of differentiation (Fig. 1C). Considering that H3K27me3 modification is catalyzed by polycomb repressive complex 2 (PRC2), we analyzed the expression of two subunits of the core PRC2, Suz12, and Ezh2, and we observed a slight but significant change of Suz12 levels that remained higher in KO cells throughout the induction of the transition (Fig. 1D). A similar behavior, although less pronounced, was also observed for Ezh2 (Fig. 1D).
Gene expression changes upon the induction of EpiLCs in the absence of Hmga2
The above-described results revealed that the absence of Hmga2 influences the abundance of the repressive histone modification H3K27me3 upon the exit of PSCs from the undifferentiated naïve state. To address whether this phenomenon is accompanied by a Hmga2-dependent change in the gene expression profile, we performed an RNA sequencing (RNA-seq) experiment to compare Hmga2 wt and KO cells at 24h after the induction of EpiLCs. The expression profiles of wt cells and two independent clones of Hmga2 KO cells were similar, with some exceptions (Fig. 2A). Indeed, we found 8341 genes whose expression was significantly modified 24 h after the exit from the naïve state of wt cells. Most of these genes were similarly modified in Hmga2 KO cells (Additional file 2: Table S1). There were 2311 genes whose expression was significantly modified only in wt cells and in none of the two Hmga2 KO clones, while the expression of 343 genes was significantly modified only in both Hmga2 KO clones and not in the wt cells. Then, we compared the list of genes, whose expression changes upon EpiLC induction appeared to be Hmga2-dependent, with the list of PRC2 putative target genes, available in GEO GSE74330 [21]. This comparison showed no significant enrichments.
Hmga2 association with chromatin in EpiSCs
There are numerous results indicating that Hmga2 interacts with chromatin at specific sites, thus contributing to gene regulation [11,12,13,14]. A genome-wide analysis of Hmga2 interaction with chromatin was already done in undifferentiated mouse ESCs overexpressing the recombinant protein [22] that demonstrated a diffuse interaction of Hmga2 with chromatin. We sought to explore endogenous Hmga2 association with chromatin and this was possible in EpiSCs (see the “Methods” section) where the levels of Hmga2 are enough high.
Thus, we performed a ChIP-seq experiment using a primary antibody against Hmga2 and this experiment allowed us to detect a diffuse interaction of endogenous Hmga2 with chromatin, confirming the previously published results obtained in undifferentiated ESCs [22]. In addition, we observed 2673 sites where Hmga2 accumulated as discrete peaks in our non-replicated experiment (2-fold, p< 0.001, see Additional file 3: Table S2). The analysis of the sequences of these Hmga2 peaks, by MEME-ChIP software, allowed us to identify three significant motifs (Fig. 2B).
These peaks were preferentially located at the inter-nucleosomal spacers, as demonstrated by the alignment of Hmga2 ChIP-seq with the data deriving from the MNase-seq experiments (ref. [23]) indicating that the number of the overlapping peaks is negligible.
A small number of Hmga2 peaks overlapped or were in close spatial relationship with putative enhancer elements marked by H3K4me1, H3K27ac, or both (Fig. 2C), while no significant overlapping between Hmga2 peaks and H3K27me3 peaks was observed. In all the cases analyzed, the histone marks of the putative enhancers were not significantly modified upon the silencing of Hmga2 in EpiSCs (Additional file 4: Fig. S2). Furthermore, no significant associations were found between the genes whose transcription start site (TTS) is located within 100 kb from these putative enhancers and the genes whose expression was different in wt vs Hmga2 KO cells.
The above-reported results do not support a direct relationship between Hmga2 association with chromatin, the observed changes in histone mark deposition and/or with gene expression modification in Hmga2 KO cells. However, it should be considered that the ChIP-seq experiment was done with chromatin from EpiSCs; thus, a different interaction pattern of Hmga2 during the early steps of this transition cannot be excluded.
Hmga2 loss alters the structure of the nuclear lamina in differentiating cells
The analysis of Hmga2 ChIP-seq data showed a significant association (~51% Spearman correlation) of this protein with non-transcribed regions and with Lamin-associated domains (LADs) [24] (Fig. 2D), confirming the results already observed in ESCs overexpressing Hmga2 [22]. Lamins guarantee the proper structure of the nuclear periphery and are essential to maintain an adequate 3D genome organization [25]. We asked whether the cells lacking Hmga2 suffer from alterations of the NL. Therefore, we performed Lamin b1 (Lmnb1) immunofluorescence on Hmga2 wild-type and knock-out cells during the transition into EpiLCs. At day 1, Hmga2 KO cells showed a remarkable distortion of the NL compared to the wt cells, with many nuclear protrusions and a clear nuclear blebbing (Fig. 3A, B).
Tridimensional analysis of the z-slices maximum projections showed that Hmga2 KO cells suffered from nuclear weakness, and they were also characterized by a reduced porosity compared to the wt cells that, instead, had an intact NL, which appeared as circular, thick, and definite (Fig. 3C).
The observed altered immunostaining of Lmnb1 was only dependent on the structure of the NL, because Lmnb1 expression was unchanged in Hmga2 KO cells compared to wt cells (Additional file 5: Fig. S3A).
Noteworthy, no differences were observed between Hmga2 wild-type and KO cells in the undifferentiated state (Additional file 5: Fig. S3B). This evidence indicated that the distortion of the NL in Hmga2 KO cells occurred only when the cells were induced into EpiLCs, thus coinciding with the block of transition caused by Hmga2 loss [16].
To confirm the data obtained in Hmga2 KO cells, we examined the effects of Hmga2 silencing in wt mouse ESCs, upon the induction of EpiLCs. Hmga2 downregulation (see Additional file 5: Fig. S3C and D) strongly altered the structure of the NL, causing nuclear abnormalities similar to those observed in Hmga2 KO cells (Fig. 3D, E).
The nuclear lamina interposes between condensed chromatin and the inner nuclear membrane [26]. To examine the status of the inner nuclear membrane, we carried out the immunostaining with the anti-Emerin antibody in Hmga2 wt and KO cells as well as in Hmga2-silenced EpiLCs, as this protein of the inner nuclear membrane is molecularly connected to the nuclear lamina. As reported in Fig. 4, the structural abnormalities detected using Lmnb1 antibody were also observed with Emerin antibody: the immunostaining revealed that the inner nuclear membrane of the Hmga2 KO and Hmga2-silenced cells suffered from the distortion of the nuclear lamina. As for Lmnb1, no nuclear envelope defects were observed in Hmga2 KO cells compared to the wild-type cells in the undifferentiated conditions, confirming that the distortion occurred during the exit from the undifferentiated state (Additional file 6: Fig. S4).
Altogether, these data demonstrate that the absence of Hmga2 inside the nucleus prevents maintaining the proper structure of the nuclear lamina.
Hmga2 loss affects the intranuclear distribution of H3K9me3
Chromatin associated with nuclear lamina is marked by heterochromatic histone modifications; thus, we investigated whether changes regarding the distribution and the abundance of these histone modifications occurred in Hmga2 KO cells. First, we conducted immunofluorescence experiments to investigate H3K9me2/3 and H3K27me3 localization in Hmga2 KO cells at day 1 of transition into EpiLCs. These experiments allowed us to detect a widespread mis-localization of H3K9me3 in differentiating Hmga2 KO cells (Fig. 5A).
In detail, in the absence of Hmga2, H3K9me3 signal was diffuse inside the nucleus and not organized in intense dots at chromocenters (pericentric constitutive heterochromatin represented by DAPI-dense foci), as observed in the wild-type cells (Fig. 5A). H3K9me2 signal was mostly located close to nuclear envelope so its distribution in KO cells confirmed the distortion of nuclear lamina (Additional file 7: Fig. S5A). These results suggested that Hmga2 loss affected H3K9me3 distribution, with effects on peripheral and internal heterochromatin. The distribution of H3K27me3 histone modification was not affected: Hmga2 KO cells were characterized by small intense dots marked by H3K27me3 whose distribution was similar to that observed in wt cells (Additional file 7: Fig. S5B). A similar behavior was also observed in case of H3K4me3 and H3K27ac, which were similar in wt and KO cells (Additional file 7: Fig. S5C-D). In addition, we explored whether the increased levels of H3K27me3 and the slight accumulation of PRC2 subunits, Suz12 and Ezh2 (see Fig. 1), could have any role in the lamina-associated phenotype observed in Hmga2 KO cells. To this aim, we inhibited PRC2 by treating the cells with GSK126 for 72 h and inducing the transition towards EpiLCs for 24h. Although the treatment led to a clear decrease of H3K27me3, we observed the same alterations of the nuclear lamina as in mock treated cells (Additional file 8: Fig. S6).
Hmga2 silencing in ESCs induced toward EpiLCs also confirmed the altered localization of H3K9me3 (Fig. 5B). It is worth noting that H3K9me3 and H3K27me3 accumulate within or close to nuclear blebs observed in Hmga2 KO cells thus confirming the heterochromatin enrichment in nuclear abnormalities (HENA) observed in other experimental settings [27].
Inter-TAD interactions and their association with nuclear lamina are altered in Hmga2 KO cells
The crosstalk between LADs with whole genome organization [25] and with the dynamics of TAD cliques during differentiation [20, 28] was clearly demonstrated in several conditions. On the other hand, in some cases, LAD borders are marked by H3K27me3 [24, 29], and the suppression of Ezh2, the catalytic subunits responsible for H3K27me3 deposition, prevents NL interaction with chromatin [27]. Thus, we hypothesized that defects in nuclear shape and heterochromatin mis-localization accompanying Hmga2 suppression could be related to defects in 3D genome assembly.
To explore the effects of Hmga2 on genome 3D structure, we performed the Hi-C mapping in two independent clones of wt cells and in two clones of Hmga2 KO cells, at day 3 after the induction of the transition into EpiLCs. After the filtering and mapping of Hi-C data of the four clones, we obtained 38,483,563; 55,451,307; 45,019,690; and 53,506,473 valid interactions, respectively (Additional file 9: Table S3).
The contact maps of all the chromosome and of an example (chromosome 11) are shown in Fig. 6A. Visual inspection of these maps indicate that the contact pattern is very similar in wt cells and in KO cells. To quantify this observation, we performed PCA analysis and identified A/B compartments genome-wide in wt and KO (Additional file 10: Fig. S7A) [30]. Comparison between compartments in different conditions reveals a general high correlation (genome-wide average Pearson r=0.96, Additional file 10: Fig. S7B), with few, short regions, typically located at the compartment periphery, exhibiting changes in compartment identity (Additional file 10: Fig. S7C, Additional file 11: Table S4). Next, we quantified the similarity of TAD boundaries [31, 32] in the different systems, called by their Insulation Score [33]. Such an analysis confirmed that TAD boundaries are similar in the four samples, with an overlap, in the case of chromosome 11, that was always close to 90% at 25 kb (see Fig. 6B and C and the “Methods” section), i.e., comparable to replicate similarity. Genome-wide analysis made with data at 25 kb indicated a global overlap close to 90%. Thus, the results indicated that the suppression of Hmga2 does not alter TAD structure. The correlation of boundaries with epigenetic marks and with the binding of Hmga2 to chromatin confirmed that H3K4me3 is tightly associated with the boundaries [31], while no significant associations were observed in the case of H3K27me3 as previously reported. Interestingly, though, Hmga2 abundance on chromatin appears to be negatively associated with the boundaries (see Fig. 6D), highlighting a novel potential involvement of Hmga2 in TAD architecture.
Next, we explored the effects of Hmga2 KO on inter-TAD interactions. To this aim, we considered the average interaction between two given TADs and calculated the log ratio (or log fold-change, indicated as Log FC) between their interaction in wt cells vs Hmga2 KO cells. Figure 7A summarizes the computational approach to address this point.
The analysis of the log fold changes allowed us to identify the TAD pairs with most prominent differences between wt and KO cells, defined as those where LogFC >1 or <−1 (Fig. 7B, chromosome 1). In a significant number of cases, the interaction present in wt cells is decreased or increased by a factor of 2, highlighting that Hmga2, although not affecting TAD boundaries, has a role in the formation of chromatin 3D structure. Considering the phenotype of Hmga2 KO EpiLCs, and the demonstrated role of Lamins in the organization of global 3D genome structure [34], we asked whether there is any relationship between the changes we observed in the interactions among TADs and the chromatin regions preferentially associated with the nuclear lamina (DamID data taken from ref. [34]). To this aim, we selected all 787 pairs of TADs whose interaction is significantly decreased in Hmga2 KO cells (Log FC >1) and all 645 pairs that interact more in Hmga2 KO cells than in wt cells (LogFC<-1). The results show that, in most cases (77%) of TAD pairs with decreased interaction in KO cells, only one of the TADs is located in a region associated with Lamin, while only 23% of the pairs are either both associated with Lamin or not associated with Lamin (Fig. 7C). On the contrary, the TAD pairs showing an increased interaction in Hmga2 KO cells compared to wt cells are almost uniformly distributed (54% with only one of the two TADs associated with Lamin and 46% both associated with Lamin or not associated with Lamin). To check the robustness of these results, we repeated the above-described analysis using a random control system where we permuted (bootstrap) the boundary positions along the genome (Additional file 12: Fig. S8). Although TAD pairs were selected with the same criterion, no specific DamID enrichment pattern is observed in TAD pairs with increased or decreased interaction in KO as compared to the wt condition.
Finally, we checked that both types of TAD pairs exhibit analogous size distributions (Kolmogoroff-Smirnov test p-val > 0.05), so to exclude size-related effects. Taken together, our results highlight that the known role of Hmga2 in defining chromatin interaction with the lamina affects the establishment of TAD borders and inter-TAD interactions.
Discussion
Hmga proteins have been studied for many years, in particular to address their role in cancer, where they are often overexpressed and were suggested to have a crucial role in tumorigenesis. However, the functions of these proteins are still not definitively addressed. In particular, their interaction with chromatin and the possible effects of this interaction have been explored in several studies. By using EpiSCs, where endogenous levels of Hmga2 are sufficiently high, we have studied the association of the protein to chromatin and we observed that it widely binds to chromatin. Only in a relatively small number of regions, an accumulation of Hmga2 was detectable, corresponding to very A/T-rich sequences. In any cases, we did not find any significant overlaps between the genes whose transcription is affected by Hmga2 suppression during differentiation and the genes whose transcription start sites are located within 100 kb from Hmga2 peaks. This suggests that the known effects of Hmga2 on the transcription of specific genes do not necessarily require a stable and continuous binding to chromatin, while Hmga2 interaction with chromatin appears to vary during the subsequent phases of the differentiation process. Thus, we cannot exclude that, in other cell types and in particular in cancer cells, Hmga2 could contribute to transcription regulation through the association with specific chromatin regions.
Our results are in good agreement with those obtained in ESCs overexpressing Hmga2 [22] that demonstrated the attitude of Hmga proteins to bind any chromatin region, with a particular preference for heterochromatic regions. In particular, the observation that Hmga2 is preferentially found in the chromatin regions associated with nuclear lamina allowed us to demonstrate that Hmga2 KO cells, upon the exit from naïve state, show evident distortions of the nuclear envelope. The latter is composed by nuclear membranes, crossed by numerous associated proteins [35], and the inner nuclear membrane is strictly connected to a filamentous meshwork of proteins, called nuclear lamina [26]. Nuclear lamina interacts with DNA, histones, transcription factors, and chromatin proteins [20], and, more in general, with large heterochromatin regions located at the nuclear periphery, defined lamina-associated domains (LADs) [36], marked by H3K9me2/me3 [29, 36]. We observed that the induction of EpiLCs from Hmga2 KO and KD PSCs is characterized by a distortion of the nuclear lamina that results in evident abnormalities of the nuclear shape. This phenomenon is accompanied by a mis-localization of H3K9me3 within the nucleus. These observations suggest that the absence of Hmga2 hampers the establishment of the proper interactions between nuclear lamina and heterochromatin.
The absence of Hmga2 interferes with the deposition and/or the removal of H3K27me3, which marks some LAD boundaries [24, 29]. This histone modification is rapidly suppressed in wt cells soon after the exit from the naïve state, while continuing to be present at relatively high levels in Hmga2 KO cells. This sustained presence of H3K27me3 could be the consequence of the persistence of the PRC2 activity, as suggested by the slightly increased levels of Suz12 and Ezh2 proteins in KO cells, likely due to post-transcriptional events. However, we did not find a significant relationship between genes regulated by PRC2 and genes whose expression was modified in Hmga2 KO cells, therefore it appears that the sustained deposition of H3K27me3 by PRC2 specifically concerns with nuclear lamina distortion and not with a general effect on bivalent promoters. The structure of Hmga2 suggests that it could be involved in the interaction of distant chromatin regions; indeed, the three AT-hooks domains present in this protein could bind distant DNA elements thus contributing to the 3D structure of the chromatin [10]. Accordingly, we observed that inter-TAD interactions are altered in the absence of Hmga2. This lack of interaction between cognate TADs and the wrong associations between TADs normally not interacting could result in a significant alteration of global 3D chromatin structure. The normal inter-TAD interactions that mostly suffer from the absence of Hmga2 are those involving one TAD endowed in Lamin-associated chromatin and the other one not associated with Lamins. This observation suggests the possibility that Hmga2 is involved in tethering euchromatic regions (TADs not directly interacting with Lamins) to the repressive environment where heterochromatic TADs are directly associated with the nuclear lamina. The increased deposition of H3K27me3, and the sustained expression of PRC2, could be interpreted as an attempt to keep repressed genes that, as a consequence of the absence of Hmga2, are localized in TADs which have lost their interaction with the heterochromatic TADs directly associated to the nuclear lamina. The compensatory effect of H3K27me3 was already observed in different experimental contexts, aimed at maintaining the silent state of altered heterochromatin [37, 38]. On the other hand, the TAD pairs showing an increased interaction in Hmga2 KO cells could result from inappropriate interactions between TADs, so explaining their uniform distribution among various chromatin regions, irrespective of their association with nuclear lamina.
The possible relationship between the loss of Hmga2 and the deposition of H3K27me3 by PRC2 is also suggested by a certain similarity of the phenotypes of the mice carrying mutations of Hmga2 or of PRC2. Indeed, Hmga2 KO mice display a pygmy phenotype, with reduced body weight and length [3] and this phenotype is exacerbated by the contemporary inactivation of the Hmga1 gene [5]. An opposite phenotype was observed in transgenic mice where a mutant form of Hmga2 is constitutively expressed that show a giant phenotype [39]. Some rare human genetic diseases, associated with prenatal and postnatal overgrowth, are known as PRC2-related overgrowth syndromes (Weaver syndrome, Cohen–Gibson syndrome, and SUZ12-related overgrowth) [40, 41]. Accordingly, mice carrying a Weaver syndrome mutation, which show a reduction of H2K27me2-3, show overgrowth and organomegaly [42]. Although the similarity of these phenotypes could be a mere coincidence, it probably deserves further studies.
The findings reported herein could also contribute to the understanding of the role of Hmga2 in cancer. It is evident that Hmga2 expression is increased in many different types of tumors and in some cases these tumors appear to be Hmga2 addicted, as its suppression results in a dramatic change in cancer behavior and aggressiveness [15]. Most results available suggest that these phenomena are the consequence of Hmga2-induced changes in gene expression. To address whether these changes are due to the direct regulation of transcription by Hmga2 and also to a general role of this protein in the 3D genome organization, as it emerges from our results, deserves further studies.
Conclusions
Here, we have demonstrated that the loss of Hmga2 leads to an inappropriate persistence of H3K27me3 histone mark after the exit from the naïve state of PSCs. Moreover, the absence of Hmga2 has a dramatic effect on the structure of NL and on the proper localization of the heterochromatin mark H3K9me3 within the nucleus. These nuclear alterations found in Hmga2 KO cells influence the inter-TAD interactions indicating that Hmga2 is required for the establishment of the proper 3D structure of the chromatin.
Methods
Cell culture, transfection, EpiLC and EpiSC generation, and SFEB differentiation
Hmga2 wt and KO induced PSCs (iPSCs) were generated as previously described [16]. The cells were grown on mitotically inactivated fibroblast feeder layers and expanded in the following iPSC medium: GMEM supplemented with 2 mM glutamine, 100 U/mL penicillin/streptomycin, 1 mM sodium pyruvate, 1X non-essential amino acids, 0.1 mM β-mercaptoethanol (all from Invitrogen), 15% FBS (Hyclone), 103 U/ml leukaemia inhibitory factor (EMD Millipore), and 2i (MEK and GSK3 inhibitors; 1 μΜ PD0325901 and 3 μΜ CHIR-99021, both from Selleckchem). Two clones of Hmga2 wt iPSCs and two clones of Hmga2 KO iPSCs were used throughout this study.
Mouse ESCs (E14Tg2a, Bay Genomics) were expanded on gelatin-coated plates in the following ESC medium: GMEM supplemented with 2 mM glutamine, 100 U/mL penicillin/streptomycin, 1 mM sodium pyruvate, 1X non-essential amino acids, 0.1 mM β-mercaptoethanol, 10% FBS, and 103 U/ml leukaemia inhibitory factor.
iPSCs and ESCs were induced into EpiLCs according to the following protocol: 2×104 cells/cm2 were plated on plates coated with Geltrex™ LDEV-Free Reduced Growth Factor Basement Membrane Matrix (Gibco) in EpiLC medium (DMEMF12/Neurobasal medium 1:1 (Gibco), 0.5% N2 supplement (Thermo Fisher), 1% B27 supplement (Thermo Fisher), BSA (Sigma Aldrich), 2 mM glutamine (Invitrogen), 20 ng/mL Activin A (R&D Systems), and 12 ng/mL bFGF (Invitrogen). For RNA and protein extraction and immunofluorescence, the cells were harvested at different time points as specified in the figures after 1, 2, and 3 days of EpiLC induction. For the treatment with PRC2 inhibitor, iPSCs were plated on gelatin-coated plates and grown in the presence or in the absence of 10 μM GSK126 (S7061, Selleckchem) for 72 h. DMSO was used as vehicle control. Then, the cells were differentiated into EpiLCs according to the protocol described above.
Hmga2 wt and KO iPSCs were differentiated into SFEBs using the protocol in ref. [43]. Briefly, 1×106 cells were plated in 100-mm Petri dishes in the following differentiation medium: Glasgow minimum essential medium supplemented with 2 mM glutamine, 1 mM sodium pyruvate, 1X non-essential amino acids, 0.1 mM β-mercaptoethanol, and 10% KO serum replacement (ThermoFisher). The aggregates were harvested for protein extraction at 2, 3, and 4 days of differentiation.
To induce Hmga2 knockdown, 2×104 cells/cm2 ESCs were plated in EpiLC medium. Twenty-four hours after plating (at 1 day of EpiLC induction), the cells were transfected with a pool of three different small interference RNAs (siRNAs, Invitrogen) directed against Hmga2 transcript (siHmga2 MSS274854, siHmga2 MSS274855, siHmga2 MSS274856) using Lipofectamine 2000 (Invitrogen), according to the manufacturer’s protocol. The cells were harvested at day 3 for protein and RNA extraction as well as for immunofluorescence experiments.
For ChIP-seq experiments Hmga2 wt and KO iPSCs were induced into EpiSCs by using the protocol in ref. [44]. In detail, 3×103/cm2 cells were grown in a medium containing N2 and B27 supplements (Thermo Scientific) complemented with 20 ng/mL Activin A (R&D Systems) and 12 ng/mL bFGF (Invitrogen). The cells were cultured in this medium for 3 days dissociated with dispase (Gibco) and cultured on ™ LDEV-Free Reduced Growth Factor Basement Membrane Matrix (Gibco)-coated plates. The cells were expanded repeating this procedure every 2–3 days for 2 weeks and then collected for ChIP-seq experiments.
Protein extraction and western blotting
For protein extraction, the cells were washed twice with ice-cold 1X phosphate saline buffer (1X PBS) and lysed with RIPA Lysis and Extraction Buffer (Thermo Scientific), complemented with protease (Sigma-Aldrich) and phosphatase inhibitors, according to the manufacturer’s protocol. After lysis, the samples were sonicated with Minisonix XL 2000 and centrifugated at 14000×g for 15 min to isolate total nuclear proteins. Bradford Protein assay (Bradford Reagent, Bio-Rad) was used to determine the concentration of total proteins. Fifteen micrograms of total protein was loaded in each lane. For western blotting experiments, the following antibodies were used: anti-Hmga2 #5269 (Cell Signaling) anti-H3K27me3 #3915S, anti-H3 (clone 0301) #39064, anti-H3K4me3 #61379, anti-H3K9me2 #39239, anti-H3K9me3 #39161 (all from Active Motif), anti-H3K27ac ab4729 (Abcam), anti-H3K27me2 ab26684, anti-Lmnb1 D4Q4Z (Cell Signalling), and anti-Gapdh sc-32233 (Santa Cruz Biotechnology). The secondary antibodies used in this study were goat anti-rabbit IgG (HeL), HRP conjugate GtxRb-003-DHRPX (ImmunoReagents) or goat anti-mouse IgG (HeL), and HRP conjugate GtxMu-003-DHRPX. Clarity TM western ECL substrate (170-5060, Bio-Rad) was used for signals detection by enhanced chemiluminescence. Western blotting images were acquired with ChemiDoc MP Imaging System (Bio-Rad). Band intensities were quantified using Image J software. The experiments were carried out in triplicate.
RNA extraction, reverse transcription, real-time qPCR, and RNA-sequencing
Total RNA was extracted using Tri-Sure (Bioline), according to the manufacturer’s protocol. Full-length first cDNA strand was synthetized by RevertAid Reverse Transcriptase (Thermo Scientific). RNA concentration and purity was determined using Nanodrop 2000c (Thermo Scientific). qRT-PCR was run on QuantStudio 7 Flex Real Time PCR System using Fast SYBR Green PCR Master Mix (Thermo Fisher Scientific). The sequence of the primers used to determine the expression of Hmga2 transcript was as follow: 5′-AAAACGGCCAAGAGGCAGAC-3′, 5′-ATGTCTCTTCAGTCTCCTGAGCA-3′. The expression of glyceraldehyde-3-phosphate dehydrogenase (Gapdh, 5′-GTATGACTCCACTCACGGCAAA-3′, 5′-TTCCCATTCTCGGCCTTG-3′) was used as internal control. Gene expression was evaluated by ΔΔCt method. The experiments were carried out in triplicate. For RNA-seq and analysis, total RNA was extracted from Hmga2 wt and KO cells induced into EpiLCs as described above and sequenced using the pipeline reported in our previous study [45].
Next-generation sequencing experiments were performed by Genomix4life S.R.L. (Baronissi, Salerno, Italy). RNA concentrations were determined by using NanoDropOne spectrophotometer (Thermo Fisher) and its quality assessed with the TapeStation 4200 (Agilent Technologies). Indexed libraries were prepared from 500 ng of purified RNA with TruSeq Stranded totalRNA Sample Prep Kit (Illumina) according to the manufacturer’s instructions. Libraries were quantified using the TapeStation 4200 (Agilent Technologies) and Qubit fluorometer (Invitrogen Co.), then pooled such that each index-tagged sample was present in equimolar amounts. The pooled samples were subject to cluster generation and sequencing using an Illumina NovaSeq6000 System (Illumina) in a 2×75 paired-end format.
Sequencing reads were processed using “nf-core’s” [46] “rnaseq” pipeline (ver.: 1.2.2) run within Dockers [47] container through Nextflow process scheduler [48]. In short, raw reads were evaluated with FastQC, the 5′ and 3′ adapters were trimmed using Trim Galore, and genome/ribosomal RNA sequences were discarded with BBSplit/SortMeRN. Filtered reads were aligned to Ensembl’s GRCm38 reference genome from Amazon’s iGenomes database using STAR sort, indexed with SAMtools, duplicate marked using Picard and gene/transcript quantified applying Salmon/StringTie. Finally, gene expression levels were compared between cell lines using DESeq2. Statistical analysis performed using the computing environment R [49].
Total RNA isolated from three independent biological replicates was used for RNA-seq experiments.
Chromatin-immunoprecipitation (ChIP)-sequencing
ChIP experiments were performed on EpiSCs. Briefly, the cells were crosslinked with 1% formaldehyde for 10 min at room temperature and next quenched by addition of glycine (125 mM). The samples were sonicated to obtain fragments with a length of 200–500bp. Chromatin extracts were incubated with anti-Hmga2 antibody (Cell Signaling, #5269, 1:100). Appropriate IgG (Abcam) were used as a negative control, while the supernatant without antibody was used as input control. Oligonucleotide pairs used are the following:
-
Peak 7040: for-TAGATCTGTGGGGGTCCAGG, rev-CTGGCCCTGTTGTAACCAGT
-
Peak 4169: for-GGGGAGAGCATGCCTTTGAT, rev- TGGAGAGATGGCTGGACTCA
-
Peak 10319: for-ACATTTACTGTTCCTCACAGCCT, rev-ACGTATGAATCCCTTCGCCA
-
Peak 10994: for-GTCTTGTCTAAACCCAGGTGGA, rev-GTCTGGAATGGACACGGTGA
-
Peak 15807: for-GAGCCCACCTACGTCAATCA, rev-ACCGAGCAAATTCCTCCCTG
-
miR-23a: for-CTTGGCTCCCTGTACCTGTC, rev-TTACCTTCTCAGGCCACCCT
ChIP-seq library was prepared using ChIP-Seq Sample Prep Kit (Illumina) and sequenced on the Illumina HiScanSQ Platform, as in ref. [50]. The obtained reads were mapped to the mouse genome (mm9 assembly). Sequencing reads were processed using “nf-core’s” [46] “chipseq” pipeline (ver.: 1.2.2) run within Docker [47] container through Nextflow process scheduler. In short, raw reads were evaluated with FastQC [48], and the 5′ and 3′ adapters were trimmed using Trim Galore. Then, trimmed reads were aligned against Ensembl’s GRCm38 reference genome from Amazon’s iGenomes database with BWA and filtered out from ambiguous reads using Picard, SAMtools, BAMtools, and Pysam. Afterwards, the peaks were called with MACS2, providing input samples from the chip-seq experiment as a control, and Hmga2 enrichments were annotated using HOMER. Finally, Hmga2 peaks were compared between cell lines using DESeq2. Statistical analysis was performed using the computing environment R [49]. H3K4me1, H3K27ac, and H3K27me3 ChIP-seq data are from ref. [19].
Immunofluorescence
Immunofluorescence experiments on undifferentiated iPSCs were conducted plating the cells on gelatin-coated Nunc TM Lab-Tek TM chamber slides (Thermo Scientific). For immunofluorescence experiments on EpiLCs, the cells were plated on coverslip coated with Geltrex™ LDEV-Free Reduced Growth Factor Basement Membrane Matrix (Gibco) and induced into EpiLCs. The cells were harvested at day 1 or 3 of EpiLC transition, washed with 1X PBS and fixed with 4% paraformaldehyde for 10 min at room temperature. Next, they were permeabilized with 0.2% Triton-X100 (PanReac Applichem) for ten minutes and incubated with the blocking solution (10% FBS, 1%BSA, 1X PBS) for 30 minutes at room temperature. For antigen detection, the cells were incubated with anti-H3K27me3 #3915S, anti-H3K9me2 #39239, anti-H3K9me3 #39161 anti-H3K27me2 ab26684 antibodies (all from Active Motif), anti-H3K27ac ab4729 (Abcam), anti-Lmnb1 (B-10) sc-374015 (Santa Cruz Biotechnology), and anti-Emerin (D3B9G) XP ® Rabbit mAB#30853 (Cell Signalling) overnight at 4°C. The day after, the coverslips were washed with 1X PBS and incubated with the secondary antibody donkey anti-Rabbit IgG (H + L), Alexa Fluor™ 488, or goat anti-Rabbit IgG (H+L) Alexa Fluor™ 594 (Invitrogen, USA) for 1 h at room temperature. After washes, the cells were stained with DAPI (268298, Calbiochem) for 5 min at room temperature. Images were captured using Leica Thunder Imaging System (Leica Microsystems) equipped with a LEICA DFC9000 GTC camera, lumencor fluorescence LED light source, and 100× oil immersion objective was used to acquire Z-slice images. Small volume computational clearing was used to remove the background signal derived from out-of-focus blur.
Hi-C sequencing
Hmga2 wt and KO cells were induced into EpiLCs and at day 3 were washed with 1X PBS and trypsinized for 5 min at 37°C. The cells were centrifuged for 10 min at 450×g at room temperature and the pellet gentled resuspended in EpiLC fresh medium. Next, the cells were cross-linked with 1% formaldehyde for 10 min at room temperature and incubated with 3M Tris pH 7.5 for 5 min at room temperature and then on ice for 15 min to quench the formaldehyde. The cells were centrifugated for 10 min at 800×g at 4° C, the pellet was resuspended in cold 1X PBS, and the cell density was determined. Hi-C sequencing was performed by Genomix4life S.R.L. (Baronissi, Salerno, Italy) [51]. Indexed libraries were prepared from 2×106 cells, with EpiTect® Hi-C (Qiagen) according to the manufacturer’s instructions. Libraries were quantified using the TapeStation 4200 (Agilent Technologies) and Invitrogen Qubit fluorometer (Thermo Fisher Scientific), then pooled such that each index-tagged sample was present in equimolar amounts, with final concentration of the pooled samples of 2nM. The pooled samples were subject to cluster generation and sequencing using an Illumina NextSeq 550 System (Illumina) in a 2×150 paired-end format. Paired-end reads were aligned to the NCBI reference sequence (mm10) using the read-analysis for Hi-C reads published [52] pipeline named Hi-C-Pro. Two different clones of Hmga2 wt cells and two different clones of Hmga2 KO cells were used for Hi-C experiments.
Hi-C data analysis
Hi-C data were analyzed using the platform Juicer [53]. For most analyses, we used Hi-C maps binned at 25kb resolution and KR normalized [54]. TAD boundaries were called using the insulation score approach [33], with a standard insulation window having a 500Kb × 500Kb size. We checked the robustness of the results by applying the procedure on Hi-C data binned at 50kb and returned, as expected, similar boundaries. Overlapping boundaries in wt and Hmga2 KO are defined as boundaries which have a genomic distance less than a tolerance threshold set to 3 bins (75kb at 25kb resolution and 150kb at 50kb resolution), i.e., if two boundaries are within this threshold they are considered overlapping, otherwise not [55]. TADs are defined as genomic regions between two consecutive boundaries. Interaction between different TADs (say TAD1 and TAD2, Fig. 7a) is defined as the average Hi-C score: \({X}_{1,2}=\frac{\sum_{i,j}{x}_{i,j}}{a1\ast a2}\) where a1 (a2) is the size of TAD1 (TAD2), xi, j is the Hi-C score between bin i and bin j, i (j) runs over the entries of the Hi-C matrix corresponding to TAD1 (TAD2) coordinates. Taking as reference the TAD boundaries called in wt conditions, fold change is the ratio between the average Hi-C scores \(F{C}_{1,2}={X}_{1,2}^{wt}/{X}_{1,2}^{Hgma2 KO}\). In order to have Log FC values symmetrically distributed around 0, the distribution is shifted so to have a null median (Fig. 7B). In general, TAD pairs with low (high) Log FC are classified as more (less) interacting in Hmga2 KO than wt condition. In addition, TAD pairs are classified according to their “TAD distance”, i.e., the number of TADs separating two TADs in the pair. For each distance, we consider the two most/least interacting pairs in KO, so to equally consider different genomic distances. Finally, we consider TADs having a size higher than 500kb. For each TAD, we evaluated lamina association level by measuring the average DamID signal (binned at 100kb, from ref. [34]) contained in the TAD. As DamID was available for mm9 assembly, we employed the liftover tool to remap our TAD boundaries from mm10 to mm9 assembly in this analysis. Average DamID levels for TAD pairs are then represented on a 2D plane (Fig. 7c). As control case, we generated control boundaries configurations by randomly permuting real boundaries on each chromosome, so to conserve the same TAD size distribution. Then, we calculated on this new arrangement of randomly located TADs the average DamID level, as described above. Enrichment profiles of H3K4me3, H3K27me3 and Hmga2 marks at wt TAD boundaries are computed by aggregating the average signal (binned at 25kb) of a 300-kb interval centered around each boundary. For each point, average and standard error is estimated over the genome-wide boundary distribution.
PCA analysis was performed on Hi-C data at 100 Kb resolution using the function eigenvector in Juicer, as described in [52]. Then, the first eigenvector was correlated with GC content (https://genome.ucsc.edu/cgi-bin/hgTables) to call A/B compartments [54,55,56,57]. Finally, genomic regions exhibiting changes in A/B compartment membership from wt to KO condition were simply defined according to the change in the eigenvector signal, i.e., regions marked as A in WT and B in KO (A→B) and, vice-versa, as B in WT and A in KO (B→A), as shown in Additional file 10: Fig. S7C.
Availability of data and materials
All data generated or analyzed during this study are included in this published article, its supplementary information files, and publicly available repositories. RNAseq, CHIPseq, and Hi-C data are available at GEO (GSE200673). Uncropped blots and supporting data values are in the Additional files 13 and 14, respectively. Publicly available datasets analyzed in this study are accessible at GEO with the reference ID GSE17051 for DamID in mouse ESCs [20], GSM1382217 and GSM1382219 for chromatin markers H3K4me1 and H3K27ac [19], GSE74330 for PRC2 ChIP-seq [21], and GSE58101 for MNase-seq [23].
Abbreviations
- ChIP:
-
Chromatin immunoprecipitation
- ESCs:
-
Embryonic stem cells
- EpiLCs:
-
Epiblast-like stem cells
- EpiSCs:
-
Epiblast stem cells
- NL:
-
Nuclear lamina
- LADs:
-
Lamin-associated domains
- PSCs:
-
Luripotent stem cells
- qPCR:
-
Quantitative real time PCR
- TADs:
-
Topologically associated domains
References
Ozturk N, Singh I, Mehta A, Braun T, Barreto G. HMGA proteins as modulators of chromatin structure during transcriptional activation. Front Cell Dev Biol. 2014;2:5.
Reeves R. Nuclear functions of the HMG proteins. Biochim Biophys Acta Gene Regul Mech. 1799;2010:3–14.
Zhou X, Benson KF, Ashar HR, Chada K. Mutation responsible for the mouse pygmy phenotype in the developmentally regulated factor HMGI-C. Nature. 1995;376:771–4.
Li Z, Gilbert JA, Zhang Y, Zhang M, Qiu Q, Ramanujan K, et al. An HMGA2-IGF2BP2 axis regulates myoblast proliferation and myogenesis. Dev Cell. 2012;23:1176–88.
Federico A, Forzati F, Esposito F, Arra C, Palma G, Barbieri A, et al. Hmga1/Hmga2 double knock-out mice display a “superpygmy” phenotype. Biol Open. 2014;3:372–8.
Phair RD, Misteli T. High mobility of proteins in the mammalian cell nucleus. Nature. 2000;404:604–9.
Narita M, Narita M, Krizhanovsky V, Nuñez S, Chicas A, Hearn SA, et al. A novel role for high-mobility group a proteins in cellular senescence and heterochromatin formation. Cell. 2006;126:503–14.
Mansoori B, Mohammadi A, Ditzel HJ, Duijf PHG, Khaze V, Gjerstorff MF, et al. HMGA2 as a critical regulator in cancer development. Genes (Basel). 2021;12:269.
Yu H, Lim HH, Tjokro NO, Sathiyanathan P, Natarajan S, Chew TW, et al. Chaperoning HMGA2 protein protects stalled replication forks in stem and cancer cells. Cell Rep. 2014;6:684–97.
Vogel B, Löschberger A, Sauer M, Hock R. Cross-linking of DNA through HMGA1 suggests a DNA scaffold. Nucleic Acids Res. 2011;39:7124–33.
Xu MY, Sharma P, Pan SQ, Malik S, Roeder RG, Martinez E. Core promoter-selective function of HMGA1 and mediator in initiator-dependent transcription. Genes Dev. 2011;25:2513–24.
Munshi N, Agalioti T, Lomvardas S, Merika M, Chen GY, Thanos D. Coordination of a transcriptional switch by HMGI(Y) acetylation. Science. 2001;293:1133–6.
Munshi N, Merika M, Yie J, Senger K, Chen G, Thanos D. Acetylation of HMG I(Y) by CBP turns off IFN beta expression by disrupting the enhanceosome. Mol Cell. 1998;2:457–67.
Yie J, Merika M, Munshi N, Chen G, Thanos D. The role of HMG I(Y) in the assembly and function of the IFN-beta enhanceosome. EMBO J. 1999;18:3074–89.
Fusco A, Fedele M. Roles of HMGA proteins in cancer. Nat Rev Cancer. 2007;7:899–910.
Navarra A, Musto A, Gargiulo A, Petrosino G, Pierantoni GM, Fusco A, et al. Hmga2 is necessary for Otx2-dependent exit of embryonic stem cells from the pluripotent ground state. BMC Biol. 2016;14:24.
Höpfl G, Gassmann M, Desbaillets I. Differentiating embryonic stem cells into embryoid bodies. Methods Mol Biol. 2004;254:79–98.
Phillips AW, Nestor JE, Nestor MW. Developing HiPSC derived serum free embryoid bodies for the interrogation of 3-D stem cell cultures using physiologically relevant assays. J Vis Exp. 2017;125:55799.
Factor DC, Corradin O, Zentner GE, Saiakhova A, Song L, Chenoweth JG, et al. Epigenomic comparison reveals activation of “seed” enhancers during transition from naive to primed pluripotency. Cell Stem Cell. 2014;14(6):854–63.
Peric-Hupkes D, Meuleman W, Pagie L, Bruggeman SW, Solovei I, Brugman W, et al. Molecular maps of the reorganization of genome-nuclear lamina interactions during differentiation. Mol Cell. 2010;38:603–13.
Kloet SL, Makowski MM, Baymaz HI, van Voorthuijsen L, Karemaker ID, Santanach A, et al. The dynamic interactome and genomic targets of Polycomb complexes during stem-cell differentiation. Nat Struct Mol Biol. 2016;23:682–90.
Colombo DF, Burger L, Baubec T, Schübeler D. Binding of high mobility group A proteins to the mammalian genome occurs as a function of AT-content. PLoS Genet. 2017;13:e1007102.
Carone BR, Hung JH, Hainer SJ, Chou MT, Carone DM, Weng Z, et al. High-resolution mapping of chromatin packaging in mouse embryonic stem cells and sperm. Dev Cell. 2014;30:11–22.
Guelen L, Pagie L, Brasset E, Meuleman W, Faza MB, Talhout W, et al. Domain organization of human chromosomes revealed by mapping of nuclear lamina interactions. Nature. 2008;453:948–51.
van Steensel B, Belmont AS. Lamina-associated domains: links with chromosome architecture, heterochromatin, and gene repression. Cell. 2017;169:780–91.
Karoutas A, Akhtar A. Functional mechanisms and abnormalities of the nuclear lamina. Nat Cell Biol. 2021;23:116–26.
Karoutas A, Szymanski W, Rausch T, Guhathakurta S, Rog-Zielinska EA, Peyronnet R, et al. The NSL complex maintains nuclear architecture stability via LaminA/C acetylation. Nat Cell Biol. 2019;21:1248–60.
Paulsen J, Liyakat Ali TM, Nekrasov M, Delbarre E, Baudement MO, Kurscheid S, et al. Long-range interactions between topologically associating domains shape the four dimensional genome during differentiation. Nat Genet. 2019;51:835–43.
Harr JC, Luperchio TR, Wong X, Cohen E, Wheelan SJ, Reddy KL. Directed targeting of chromatin to the nuclear lamina is mediated by chromatin state and A-type lamins. J Cell Biol. 2015;208:33–52.
Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93. https://doi.org/10.1126/science.1181369.
Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80.
Nora EP, Lajoie BR, Schulz EG, Giorgetti L, Okamoto I, Servant N, et al. Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature. 2012;485:381–5.
Crane E, Bian Q, McCord RP, Lajoie BR, Wheeler BS, Ralston EJ, et al. Condensin-driven remodelling of X chromosome topology during dosage compensation. Nature. 2015;523:240–4.
Zheng X, Hu J, Yue S, Kristiani L, Kim M, Sauria M, et al. Lamins organize the global three-dimensional genome from the nuclear periphery. Mol Cell. 2018;71:802–815.e7.
Dickmanns A, Kehlenbach RH, Fahrenkrog B. Nuclear pore complexes and nucleocytoplasmic transport: from structure to function to disease. Int Rev Cell Mol Biol. 2015;320:171–233.
Kind J, Pagie L, Ortabozkoyun H, Boyle S, de Vries SS, Janssen H, et al. Single-cell dynamics of genome-nuclear lamina interactions. Cell. 2013;153:178–92.
Saksouk N, Barth TK, Ziegler-Birling C, Olova N, Nowak A, Rey E, et al. Redundant mechanisms to form silent chromatin at pericentromeric regions rely on BEND3 and DNA methylation. Mol Cell. 2014;56:580–94.
Peters AHFM, Kubicek S, Mechtler K, O'Sullivan RJ, Derijck AAHA, Perez-Burgos L, et al. Partitioning and plasticity of repressive histone methylation states in mammalian chromatin. Mol Cell. 2003;12:1577–89.
Battista S, Fidanza V, Fedele M, Klein-Szanto AJ, Outwater E, Brunner H, et al. The expression of a truncated HMGI-C gene induces gigantism associated with lipomatosis. Cancer Res. 1999;59:4793–7.
Cyrus S, Burkardt D, Weaver DD, Gibson WT. PRC2-complex related dysfunction in overgrowth syndromes: A review of EZH2, EED, and SUZ12 and their syndromic phenotypes. Am J Med Genet C: Semin Med Genet. 2019;181:519–31.
Cyrus SS, Cohen ASA, Agbahovbe R, Avela K, Yeung KS, Chung BHY, et al. Rare SUZ12 variants commonly cause an overgrowth phenotype. Am J Med Genet C: Semin Med Genet. 2019;181:532–47.
Lui JC, Barnes KM, Dong L, Yue S, Graber E, Rapaport R, et al. Ezh2 mutations found in the weaver overgrowth syndrome cause a partial loss of H3K27 histone methyltransferase activity. J Clin Endocrinol Metab. 2018;103:1470–8.
Parisi S, Passaro F, Russo L, Musto A, Navarra A, Romano S, et al. Lin28 is induced in primed embryonic stem cells and regulates let-7-independent events. FASEB J. 2017;31:1046–58.
Acampora D, Di Giovannantonio LG, Simeone A. Otx2 is an intrinsic determinant of the embryonic stem cell state and is required for transition to a stable epiblast stem cell condition. Development. 2013;140:43–55.
Passaro F, De Martino I, Zambelli F, Di Benedetto G, Barbato M, D'Erchia AM, et al. YAP contributes to DNA methylation remodeling upon mouse embryonic stem cell differentiation. J Biol Chem. 2021;296:100138.
Ewels PA, Peltzer A, Fillinger S, Patel H, Alneberg J, Wilm A, et al. The nf-core framework for community-curated bioinformatics pipelines. Nat Biotechnol. 2020;38:276–8.
Merkel D. Docker: lightweight linux containers for consistent development and deployment. Linux J. 2014;239:2.
Di Tommaso P, Chatzou M, Floden EW, Prieto Barja P, Palumbo E, Notredame C, et al. Nextflow enables reproducible computational workflows. Nat Biotechnol. 2017;35:316–9.
Core Team R. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2022.
Krepelova A, Neri F, Maldotti M, Rapelli S, Oliviero S. Myc and max genome-wide binding sites analysis links the Myc regulatory network with the polycomb and the core pluripotency networks in mouse embryonic stem cells. PLoS One. 2014;9:e88933.
Belaghzal H, Dekker J, Gibcus JH. Hi-C 2.0: An optimized Hi-C procedure for high-resolution genome-wide mapping of chromosome conformation. Methods. 2017;123:56–65.
Servant N, Varoquaux N, Lajoie BR, Viara E, Chen CJ, Vert JP, et al. Hi-C-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015;16:259.
Durand NC, Shamim MS, Machol I, Rao SS, Huntley MH, Lander ES, et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 2016;3:95–8.
Knight PA, Ruiz D. A fast algorithm for matrix balancing. IMA J Numer Anal. 2012;33:1029–47.
Fraser J, Ferrai C, Chiariello AM, Schueler M, Rito T, Laudanno G, et al. Hierarchical folding and reorganization of chromosomes are linked to transcriptional changes in cellular differentiation. Mol Syst Biol. 2015;11:852.
Beagrie RA, Scialdone A, Schueler M, Kraemer DC, Chotalia M, Xie SQ, et al. Complex multi-enhancer contacts captured by genome architecture mapping. Nature. 2017;543:519–24. https://doi.org/10.1038/nature21411.
Imakaev M, Fudenberg G, McCord RP, Naumova N, Goloborodko A, Lajoie BR, et al. Iterative correction of Hi-C data reveals hallmarks of chromosome organization. Nat Methods. 2012;9:999–1003. https://doi.org/10.1038/nmeth.2148.
Acknowledgements
We thank V.E. Avvedimento for helpful suggestions.
Funding
This work was supported by PRIN 2017CH4RNP from Italian Ministry of Research to SP, PRIN 20157JF8P5 from Italian Ministry of Research to TR, POR-Campania SATIN to TR, Associazione Italiana Ricerca sul Cancro (AIRC) IG 2017-20240 to SO, and PRIN 217 P352Z4 from Italian Ministry of Research to SO.
Author information
Authors and Affiliations
Contributions
Conceptualization, supervision, and writing the manuscript: SP2 and TR. Performing experiments and analyzing results: GD, AMC, AE, SP1, and FP. Analyzing results: PZ, FZ, GP, DI, and SO. Critical reviewing of the manuscript: SO, MN, and AMC. Funding acquisition: SO, SP2, and TR. The authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Hmga2 western blot in wt and KO cells; quantification of western blot of Fig. 1A. (A) western blot analysis of wt and KO cells with Hmga2 ab; Gapdh was used as loading control. (B) The data in the graphs represent the mean ± SD (n=3 biological replicates) of H3K27me3, H3K4me3, H3K9me2, H3K9me3 and H3K27ac against total H3. Student’s t-test, two tailed (ns: not significant, *p<0.05).
Additional file 2: Table S1.
Expression changes at 24h after the induction of differentiation in Hmga2 wt and KO cells.
Additional file 3: Table S2.
Hmga2 ChIP-seq data.
Additional file 4: Figure S2.
Histone marks at Hmga2 peaks in wt and Hmga2 KD EpiSCs. Chromatin from EpiSCs transfected with control (WT) and Hmga2 targeting siRNAs (KD) was immunoprecipitated with IgG or anti-H3K4me1 or anti-H3K27ac antibodies. Immunoprecipitated DNA was amplified by using oligo pairs annealing to the indicated Hmga2 peaks, which were associated by ChIP-seq with H3K4me1 and/or H3K27ac peaks. A region overlapping the miR-23a gene enhancer was used as a control; Hmga2 peak 10319, not associated with histone modifications, was used as a negative control.
Additional file 5: Figure S3.
Lmnb1 abnormalities occur during the transition from the naïve state to the primed state of pluripotency. (A) Western blotting image showing Lmnb1 expression in Hmga2 wt and KO cells at day 3 after the induction of EpiLCs. Two different clones of PSCs were induced into EpiLCs for each condition. Western blotting quantification revealed no differences in terms of Lmnb1 expression between Hmga2 wild type and knock-out cells. The data in the graphs represent the mean ± SD of the Lmnb1 signal against Gapdh (n=3 biological replicates). ns: not significant (Student’s t-test, two tailed). (B) Lmnb1 immunofluorescence (red) on undifferentiated Hmga2 wt and KO cells showing no structural abnormalities of the nuclei. DAPI (blue) was used to counterstain the nuclei. A single plane of z-stack projection is shown. Scale bar = 50 μm. The data in the graphs represent the mean ± SD of the Lmnb1 signal against Gapdh (n=3 biological replicates). ns: not significant (Student’s t-test, two tailed). (C) q-PCR analysis to evaluate the levels of Hmga2 mRNA in EpiLCs upon transfection of siRNA. Non-silencing siRNA was used as control. The data in the graphs are reported as mean ± SD of relative values (fold changes) obtained normalizing the data of Hmga2 silencing with those of non-silencing control siRNA. **p<0.01, Student’s t-test, two tailed. (D) Western blotting image showing Hmga2 protein expression levels in EpiLCs upon transfection of non-silencing control siRNA and Hmga2 siRNA, respectively.
Additional file 6: Figure S4.
Absence of nuclear alterations in Hmga2 KO undifferentiated cells. Lmnb1 (red) and Emerin (green) immunofluorescence on undifferentiated Hmga2 wt and KO cells showing the absence of nuclear abnormalities. DAPI (blue) was used to counterstain the nuclei. Maximum projection of z-slices (ROI 1024x1024) is shown. Scale bar = 50 μm.
Additional file 7: Figure S5.
Hmga2 wt and KO cells showed a similar localization of several histone H3 modifications upon the induction of EpiLCs. Immunoflurescence experiments on Hmga2 wt and KO cells at day 1 after the induction of EpiLCs. Lmnb1 (red) and H3K4me3, H3K9me2, H3K27me3 and H3K27ac (green) as indicated. Scale bars=50 μm. Quantification graphs obtained by counting >200 nuclei/condition. Error bars represent standard deviation. Statistical significance on three biological replicate experiments was determined using the student’s t-test, ns: not significant).
Additional file 8: Figure S6.
Treatment of wt and KO cells with PRC2 inhibitor and its effect on nuclear lamina. (A) Western blotting image showing the level of H3K27me3 histone modifications in Hmga2 wt and KO PSCs treated or not with PRC2 inhibitor at day 0 and at day 1 after the induction of the transition into EpiLCs. The levels of the histone H3 were used as control. (B) Immunoflurescence experiments on Hmga2 wt and KO cells treated or not with PRC2 inhibitor at day 1 after the induction of EpiLCs showing the nuclear lamina phenotype. All pictures are shown as single z-plane, ROI 1024x1024, scale bar=50μm. Quantification graph showing the percentage of cells with nuclear lamina distortion in Hmga2 wt and KO cells treated or not with PRC2 inhibitor (n≤100 nuclei per condition). The count was performed at day 1 of EpiLC transition on three biological replicates. Error bars represent standard deviation. Statistical significance was determined using the student’s t-test (ns: not significant).
Additional file 9: Table S3.
Hi-C data summary.
Additional file 10: Figure S7.
Identification of A/B compartment. (A) 1st eigenvector from the PCA analysis is shown in wt and KO cells. GC content (reported below) is used to call A/B compartment. All signals are binned at 100kb. (B) Pearson correlation between the 1st eigenvector in wt and KO for each chromosome. Dashed line represents the average value. (C) Examples of genomic regions exhibiting changes in A/B compartment membership from wt to KO condition.
Additional file 11: Table S4.
Genomic coordinates of A->B and B->A compartment changes observed in Hmga2 KO cells.
Additional file 12: Figure S8.
Simulated comparison of inter-TAD interactions. The analysis reported in panel C of Fig. 7 was repeated using a random control system where the boundary positions along the genome were permuted.
Additional file 13.
Uncropped blots.
Additional file 14.
Data.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Divisato, G., Chiariello, A.M., Esposito, A. et al. Hmga2 protein loss alters nuclear envelope and 3D chromatin structure. BMC Biol 20, 171 (2022). https://doi.org/10.1186/s12915-022-01375-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12915-022-01375-3