Reciprocal Regulation of the Cardiac Epigenome by Chromatin Structural Proteins Hmgb and Ctcf

Transcriptome remodeling in heart disease occurs through the coordinated actions of transcription factors, histone modifications, and other chromatin features at pathology-associated genes. The extent to which genome-wide chromatin reorganization also contributes to the resultant changes in gene expression remains unknown. We examined the roles of two chromatin structural proteins, Ctcf (CCCTC-binding factor) and Hmgb2 (high mobility group protein B2), in regulating pathologic transcription and chromatin remodeling. Our data demonstrate a reciprocal relationship between Hmgb2 and Ctcf in controlling aspects of chromatin structure and gene expression. Both proteins regulate each others' expression as well as transcription in cardiac myocytes; however, only Hmgb2 does so in a manner that involves global reprogramming of chromatin accessibility. We demonstrate that the actions of Hmgb2 on local chromatin accessibility are conserved across genomic loci, whereas the effects on transcription are loci-dependent and emerge in concert with histone modification and other chromatin features. Finally, although both proteins share gene targets, Hmgb2 and Ctcf, neither binds these genes simultaneously nor do they physically colocalize in myocyte nuclei. Our study uncovers a previously unknown relationship between these two ubiquitous chromatin proteins and provides a mechanistic explanation for how Hmgb2 regulates gene expression and cellular phenotype. Furthermore, we provide direct evidence for structural remodeling of chromatin on a genome-wide scale in the setting of cardiac disease.

Transcriptome remodeling in heart disease occurs through the coordinated actions of transcription factors, histone modifications, and other chromatin features at pathology-associated genes. The extent to which genome-wide chromatin reorganization also contributes to the resultant changes in gene expression remains unknown. We examined the roles of two chromatin structural proteins, Ctcf (CCCTC-binding factor) and Hmgb2 (high mobility group protein B2), in regulating pathologic transcription and chromatin remodeling. Our data demonstrate a reciprocal relationship between Hmgb2 and Ctcf in controlling aspects of chromatin structure and gene expression. Both proteins regulate each others' expression as well as transcription in cardiac myocytes; however, only Hmgb2 does so in a manner that involves global reprogramming of chromatin accessibility. We demonstrate that the actions of Hmgb2 on local chromatin accessibility are conserved across genomic loci, whereas the effects on transcription are loci-dependent and emerge in concert with histone modification and other chromatin features. Finally, although both proteins share gene targets, Hmgb2 and Ctcf, neither binds these genes simultaneously nor do they physically colocalize in myocyte nuclei. Our study uncovers a previously unknown relationship between these two ubiquitous chromatin proteins and provides a mechanistic explanation for how Hmgb2 regulates gene expression and cellular phenotype. Furthermore, we provide direct evidence for structural remodeling of chromatin on a genome-wide scale in the setting of cardiac disease.
The genome is organized hierarchically with the functional unit being the nucleosome, a heteromultimeric complex of two copies each of four core histones around which ϳ147 base pairs of DNA wrap (1). The histone tails can be modified according to many combinations, resulting in recruitment of chromatin remodeling complexes that in turn determine the spacing of nucleosomes along the genome (2). Additional structural regulation of transcription is conferred by DNA methylation and nucleosome packaging (3) at a larger scale by the formation of chromatin domains (4) and chromosomal territories (5), which controls the interaction of genes with distal regulatory DNA sequences and chromatin-modifying proteins. Chromatin structural proteins in turn regulate these multiple tiers of genomic structure to influence cell type-specific transcription, although the mechanisms for this phenomenon are poorly understood.
Transcriptome remodeling during pathologic stress to the heart has been well documented, as has the role of histonemodifying proteins in this process (6 -8). However, chromatin structure remodeling in disease requires coordination with other chromatin features, including DNA methylation (9) and transcription factors (10). How chromatin structure is reorganized in a genome-wide manner to carry out disease-associated gene expression remains poorly understood. In cardiac development, the patterning of histone modifications changes as the cell commits to a lineage, molding the transcriptome for the appropriate phenotype (11,12). Similar observations have been made in disease; in a transverse aortic constriction model of heart failure in mice, alterations to cis-acting histone posttranslational modifications coordinate expression changes of 325 genes (13).
However, chromatin regulation in the setting of cardiac disease involves features more complex than the effects of a local histone modification on gene expression, notably, higher order chromatin structure, which remains unexplored in the heart. Different cell types ought to have different chromatin structure underpinning their different transcriptomes, although most of the knowledge on endogenous chromatin structure (4,14,15) and the non-nucleosome proteins that regulate it (16,17) comes from non-cardiac cells, with some exceptions (18). At the level of the whole nucleus, we have reported a decrease in trimethylation of lysine 9 on histone H3 (H3K9me3), 8 a marker of constitutively silenced DNA (19), and an increase in H3K4me3 abundance, a marker of active expression (20), in failing hearts after transverse aortic constriction (21). Similarly, decreased H3K9me2 and increased H3K4me2 in the heart were observed in a mouse model of diabetes with glomerulosclerosis, a condition that can lead to heart disease in humans and induces hypertrophy of the cardiomyocytes in mice (22). Nuclear organization in a larger context is critical for cardiomyocyte function, as evinced by cardiomyopathy resulting from mutations in the nuclear envelope protein lamin (23) and, more recently, by studies of high mobility group nucleosomebinding domain-containing protein 5 (Hmgn5), indicating that chromatin decompaction drives disease by upsetting the normal role of heterochromatin to withstand the forces of myocyte contraction (24). These observations together suggest a more plastic chromatin environment, on a genome-wide scale, underlies gene expression remodeling during heart failure.
CCCTC-binding factor (Ctcf) is an 11 zinc finger protein that organizes higher order chromatin structure by one or more of the following actions: insulating genes from their enhancers (25); orchestrating DNA looping to bring together genes and their regulatory elements (26); and/or localizing to the boundaries between heterochromatin (compact and silenced DNA) and euchromatin (loosely packed and accessible DNA) to prevent heterochromatin spreading (27). Despite the well established role for Ctcf in genome organization, virtually nothing is known about its role in the normal or diseased cardiomyocyte.
High mobility group protein B2 (Hmgb2) is a non-nucleosomal chromatin structural protein, which, by binding to and bending DNA, can alter gene expression (28). We previously found that Hmgb2 abundance is altered in heart disease (21), and in this study we sought to test the molecular mechanisms for its actions at the level of both the chromatin fiber and the entire genome. Herein, we uncover a previously unknown relationship between Hmgb2 and Ctcf, and we use this relationship to explore the role of chromatin architectural proteins in regulating cardiac gene expression in disease.

Results
Hmgb2 and Ctcf Are Inversely Regulated in the Heart-To uncover the contribution of Ctcf and Hmgb2 to the cardiac phenotype, we analyzed microarray data (29) from the hearts of 84 classical inbred and recombinant strains of mice in the basal state and mice that were treated for 3 weeks with isoproterenol, a ␤-adrenergic agonist that increases the inotropy of the heart and is used to model hypertrophy and failure in animal models (30). In the basal condition, hmgb2 and ctcf mRNA levels are inversely correlated with each other, and this relationship is weakened, but still significant, after isoproterenol treatment (Fig. 1A). This correlation is striking given that 66 of the 84 strains down-regulated Ctcf, whereas the response of Hmgb2 is genotype-dependent (a finding validated in NRVMs, in which Ctcf was down-regulated at the protein level in response to isoproterenol, phenylephrine, or endothelin-1 (data not shown), and the response of Hmgb2 is agonist-dependent (21)). hmgb2 mRNA abundance correlated with that of neither hmgb1 nor hmgb3 in the basal setting or after isoproterenol (Fig. 1A). Strains that respond to isoproterenol by going into heart failure exhibit a significantly direct correlation between the ratio of Ctcf to Hmgb2 and heart size, whereas the mice that are resistant or develop hypertrophy show the opposite trend (Fig. 1B). In support of these observations having functional significance in the heart, microarray data from the same panel of mice taken from other tissues (73 strains analyzed for liver (31) and 98 strains analyzed for macrophages (32)) showed that the relationship between hmgb2 and ctcf levels is organ-dependent (Fig. 1C). Immunohistochemistry to label for Hmgb2 and Ctcf in mouse cardiac tissue sections confirmed that both proteins are expressed in the nuclei of adult myocytes (Fig. 1D).
We next sought to functionally validate the inverse regulation between Hmgb2 and Ctcf. We used adenoviruses to overexpress or siRNAs to knock down hmgb2 and ctcf in neonatal rat ventricular myocytes (NRVMs) (Fig. 2). Increased Hmgb2 expression resulted in a dose-dependent reduction in Ctcf at the protein and mRNA level, whereas hmgb2 knockdown caused up-regulation of ctcf at the mRNA level with no change in protein by 72 h (Fig. 2, A and B). hmgb2 knockdown did not affect levels of histone H1, another chromatin structural protein (Fig. 2B). Likewise, ctcf knockdown up-regulated and ctcf overexpression down-regulated Hmgb2 (Fig. 2C). By microscopy, we observed an increase in the overall abundance of Hmgb2 in nuclei depleted of Ctcf (Fig. 2D). Finally, in cardiacspecific ctcf knock-out mice, we observe a doubling (204.9%) in hmgb2 abundance by RNA-seq from isolated adult cardiomyocytes. Together, these findings extend our observation of an endogenous inverse relationship between Ctcf and Hmgb2 levels in the mouse heart by demonstrating that this relationship is dynamic and responsive to experimental perturbation.
Hmgb2 and Ctcf Target the Same Loci at Distinct Times-We next examined available Ctcf chromatin immunoprecipitation and DNA sequencing (ChIP-seq) data from human CD4 ϩ cells (gene expression omnibus accession GSM325895), adult mouse heart (UCSC accession wgEncodeEM001684), mouse embryonic stem cells (gene expression omnibus accession GSM69916), and rat liver (33). We performed our own ChIPseq for Hmgb2 in NRVMs and used liftOver to compare Hmgb2 binding to Ctcf peaks from the other four samples. In all four comparisons, Hmgb2 reads were strongly enriched around Ctcf binding peaks as compared with a randomized set of reads of similar length and number (Fig. 3A). We also compared Ctcf ChIP-seq data in the heart to Hmgb2 binding peaks separated by whether they fell in genes, promoters, or intergenic regions, and we found the greatest enrichment in intergenic regions using basal Hmgb2 ChIP-seq and promoter regions using Hmgb1 ChIP-seq from hypertrophic (phenylephrine-treated) NRVMs (Fig. 3B). For comparison, we mapped Hmgb2 enrichment around peaks for Nkx2.5, a cardiac transcription factor, in HL1 cells ((34) an atrial myocyte cell line) and did not see enrichment (Fig. 3A). Together, these suggest that Hmgb2 and Ctcf bind the same regions of the genome and that a portion of Hmgb2 peaks may be cell type-independent.
We validated the Hmgb2 ChIP-seq by selecting 20 peaks that were within 2 kb upstream of the TSS and performing ChIP- Basal CTCF 5 FIGURE 1. Ctcf and Hmgb2 are coregulated in the mouse heart. A, hmgb2 and ctcf abundances exhibit an inverse relationship in the basal state that is maintained after isoproterenol treatment. Plotted are microarray data for hmgb2 (x axis) and ctcf (y axis) across 84 mouse strains in the basal state or after treatment with isoproterenol. As controls, we found that hmgb2 abundance had relationship to abundance of neither hmgb1 nor hmgb3 (yellow indicates p value Ͻ0.01 (converted from R 2 value), red line represents linear regression). B, strains were grouped by their response to isoproterenol: hypertrophy (n ϭ 22), failure (n ϭ 13), minimal change or resistant (n ϭ 9), unassigned (showed traits of different disease states, n ϭ 40, not shown), and their cardiac phenotype compared with hmgb2 and ctcf abundance (x axis ratio of ctcf/hmgb2, y axis heart/heart chamber mass normalized to body weight). Total heart mass and left ventricular mass normalized to body weight showed significant positive correlation with the ratio of ctcf to hmgb2 expression in the isoproterenol-treated hearts for failing mice, but not for hypertrophic or resistant mice (p value Ͻ0.05 indicated above each plot and color-coded by strain subset; line represents linear regression). C, unlike the heart, the liver (73 mouse strains) showed no correlation between hmgb2 and ctcf abundance, and the bone marrow (98 mouse strains) had a direct correlation. D, immunohistochemistry demonstrates abundant nuclear expression of Ctcf (left) and Hmgb2 (right) in myocytes in tissue sections from adult mouse heart. Bar, 25 m.
PCR in NRVMs. In our analysis, 19 of the 20 peaks showed enriched pulldown over an IgG control. ChIP-PCR for Ctcf showed Ctcf binding five of the 19 promoter regions (Fig. 3C). These five promoters also showed Ctcf binding in the adult mouse heart (UCSC accession wgEncodeEM001684). We then examined whether Hmgb2 and Ctcf co-occupied these promot-ers in cardiac cells using ChIP-reChIP (Fig. 3C). We immunoprecipitated Ctcf with one antibody, eluted the protein and DNA complex from the beads, and then reimmunoprecipitated for Ctcf using a different antibody. There was some loss of DNA in this process, but in all cases the ChIP-reChIP successfully pulled down the five promoter sequences, serving as a positive control for the assay. We then immunoprecipitated Hmgb2 with one antibody, followed by immunoprecipitation with a second Hmgb2 antibody. In this case, the loss in sample was greater, due to the poor utility of the second Hmgb2 antibody for immunoprecipitation. Despite these limitations, we still achieved enrichment of one of the five Hmgb2 targets. Finally, we immunoprecipitated for Hmgb2 followed by immunoprecipitation for Ctcf (and vice versa, with the first immunopre-  . Ctcf and Hmgb2 can occupy the same loci but not coincidently. A, ChIP-seq data for Hmgb2 in NRVMs was compared with published ChIP-seq data for Ctcf in multiple tissues. Hmgb2 reads (blue) are enriched around Ctcf binding peaks (point 0 on x axis) but not peaks for the cardiac transcription factor Nkx2.5. Randomly generated reads of the same number and size as the Hmgb2 dataset show no enrichment (red). B, Hmgb2 ChIP-seq peaks were subset by whether they localized to gene bodies, promoters (Ϫ2 kb, ϩ500 bp of TSS), or intergenic regions. Ctcf ChIP-seq reads from the adult mouse heart were aligned across Hmgb2 ChIP-seq peaks from basal or phenylephrine-treated NRVMs. C, ChIP-reChIP experiments in NRVMs for Ctcf and Hmgb2 revealed both proteins could bind to the promoter of the five indicated genes; however, there was a loss of DNA recovery at these five loci when immunoprecipitation for Ctcf was followed by immunoprecipitation for Hmgb2 (orange, white arrowhead) but not when followed by a second Ctcf immunoprecipitation (purple, black arrowhead) as compared with Ctcf (green) or Hmgb2 (red) immunoprecipitation alone. This suggests Ctcf does not bind these regions at the same time as Hmgb2. No binding was found in a negative control region, chosen based on absence of reads in the Hmgb2 ChIP-seq data. Data correspond to one representative ChIP assay from a total of two independent assays, each of them with different NRVM isolates. Data are shown as average Ϯ S.D. Statistical significance was assessed by two-tailed Student's t test: *, p Ͻ 0.05. D, immunolabeling for Hmgb2 and Ctcf in NRVMs was detected by STED microscopy (zoom in right panel) and confirmed a lack of colocalization of these two proteins (colocalization would appear as yellow). Bar, 5 m. E, for each Hmgb2 cluster we calculated the distance to the closest Ctcf cluster. Histogram represents distribution of all clusters across six nuclei. The same calculation was repeated for each Ctcf cluster. Only ϳ9% of clusters showed colocalization (less than 50 nm, the resolution of the STED microscope). cipitation for Ctcf followed by immunoprecipitation for Hmgb2). In these experiments, the better performing Hmgb2 antibody was used. However, unlike the control experiments using different Ctcf antibodies against the same protein, here we saw loss of enrichment when immunoprecipitating for both Ctcf and Hmgb2 on the same sample. Together, these experiments indicate that both Ctcf and Hmgb2 bind these five regions, but not at the same time, hence the ability of Ctcf and Hmgb2 to pull down non-overlapping pools of these DNA fragments. Super-resolution imaging of immunolabeled Hmgb2 and Ctcf in NRVM nuclei confirmed the lack of colocalization of Hmgb2 and Ctcf in cardiomyocytes (Fig. 3D); only ϳ9% of the Hmgb2 or Ctcf puncta were within 50 nm of each other (we define colocalization as within 50 nm based on the resolution of the microscope, Fig. 3E). HMGB2 and CTCF also did not colocalize in 293T cells (data not shown). By contrast, the ChIP-seq data predicts 20 -24% overlap of these two proteins (ϳ41,000 total Hmgb2 peaks, ϳ33,000 total cardiac Ctcf peaks, and ϳ8000 overlapping Hmgb2 and Ctcf peaks indicates ϳ20% of Hmgb2 peaks and ϳ24% of Ctcf peaks should overlap).
Hmgb2 Regulates Ribosomal RNA Transcription-To characterize the phenotypic implications of disrupting the balance of Hmgb2 or Ctcf in myocytes, we analyzed the global effect of Hmgb2 and Ctcf on cardiac gene expression. Previous studies in mouse embryonic fibroblasts have shown that loss of Ctcf increases the nucleolar area; however, in that cell type, the nucleoli number decreased (35). Furthermore, Ctcf depletion in ES cells causes a modest decrease in ribosomal RNA transcripts (35).
In mammals, the ribosome is made up of ϳ85 proteins (30 -50 for 40S subunit and 40 -50 for 60S subunit) and four ribosomal RNAs (rRNAs: 18S, 5S, 5.8S, and 25S) (36). The sequences of DNA encoding these rRNAs colocalize within the nucleus, forming so-called "nucleolar organizing regions" surrounding the nucleolus (37). Mammalian genomes contain Ͼ100 repeats of ribosomal DNA units (38). In the rat, the 45S rDNA gene (the precursor for 18S, 5.8S, and 28S) is clustered in repeats on chromosomes 3, 11, and 12, with ϳ35 copies per chromosome (NCBI and RefSeq 2012). Like mRNA, rRNA expression is regulated by histone post-translational modifications and DNA methylation (39). The majority of gene copies are silenced in mammalian cells (40).
5Ј-Fluorouridine is a uracil analogue that incorporates into newly transcribed RNA when added to the media of living cells. We overexpressed Gfp alone or Hmgb2 or Ctcf tagged with GFP in 293T cells, incubated them in 5Ј-fluorouridine for 30 min and then stopped the reaction, fixed the cells, and used immunocytochemistry to detect the localization and intensity of 5Ј-fluorouridine labeling. We found no overt change to 5Ј-fluorouridine signal in either the ctcf or Gfp-only overexpressing cells. Similarly, ctcf knockdown in NRVMs had only a modest effect (11% decrease in 5Ј-fluorouridine intensity, p Ͻ 0.001, n ϭ 264 control, n ϭ 287 ctcf knockdown). By contrast, we found a stark absence of transcription in the nucleoli of hmgb2-overexpressing cells (Fig. 4A, open arrowheads). Cells in the same plate that did not actively express the hmgb2 plasmid, and which therefore were not green, also did not show transcriptional inhibition (Fig. 4A, solid arrowheads). We quanti-fied the top 250 cells with the most Gfp expression. There was Ͼ50% reduction in the mean 5Ј-fluorouridine intensity of hmgb2-overexpressing cells (Fig. 4B) and an inverse relationship between 5Ј-fluorouridine intensity and Gfp intensity in the overexpressing cells (Fig. 4B). Both of these observations were also true when analyzing all (n ϭ 329) hmgb2-overexpressing cells, rather than focusing only on the cells with the greatest overexpression (reduction in median level of 5ЈFU by 49% and a significant correlation between GFP and 5ЈFU in overexpressing cells, p value Ͻ0.001). Analyses of hmgb2 overexpression in NRVMs showed the same effect (Fig. 4, C and D).
This global reduction in transcription agreed with our previous finding (21) that hmgb2 knockdown in NRVMs increased the abundance of H3K4 methylation, a modification associated with active promoters and enhancers (20,41,42), whereas loss of Hmgb2 decreased the abundance of H3K9me3, a marker of constitutively silenced DNA (19). In this study, we now show that hmgb2 knockdown also decreases the abundance of H3K27me3 (Fig. 4E), a marker of facultative heterochromatin (43), i.e. heterochromatin more likely to be dynamically regulated over the lifetime of the cell.
We repeated the global transcription analyses in NRVMs after hmgb2 knockdown. hmgb2 knockdown decreased nucleolar transcription with no change in total transcription (Fig. 5A, panels i-ii), increasing the ratio of nucleoplasmic transcription to nucleolar transcription (Fig. 5A, panel iv). Localization of transcription was determined by costaining for nucleolin to label nucleoli. We tested whether the loss of nucleolar transcription could be due to alterations in the nucleolar structure, and we found no significant difference in nuclear, or nucleolar, size (Fig. 5B, panels i-iii). However, we did find fewer nucleoli on average per nuclei in the knockdown cells (Fig. 5B, panel iv). Additionally, we observed a reduction in the total levels of nucleolin ( Fig. 5B, panel v). The ratio of nucleolin localized within the nucleoli versus within the nucleoplasm did not change with knockdown. This suggests that hmgb2 knockdown is disrupting rRNA transcription without disrupting the gross morphology of the nucleoli. Whether the decrease in nucleolin levels with hmgb2 knockdown is a cause or a consequence of alterations in rRNA transcription remains unknown.
During cardiac hypertrophy, both translational efficiency (translational rate of ribosomes) and translational capacity (number of ribosomes) increase to support the elevated protein synthesis necessary for cellular growth (44). rRNA transcription and subsequent synthesis of new ribosomes are necessary for cardiomyocyte hypertrophy induced by phenylephrine (45), and treatment with endothelin-1 and angiotensin II also increases rRNA synthesis (44). Thus, the regulation of rRNA synthesis by Hmgb2 could be one mechanism by which Hmgb2 regulates hypertrophy.
We next sought to understand the observation that both hmgb2 knockdown and overexpression decreased ribosomal transcription. We examined Hmgb2 binding to ribosomal genes and found enrichment across the rDNA gene bodies in NRVMs as compared with the entire genome (Fig. 5C). We hypothesized that these differences could be explained by a concentration-dependent functionality of Hmgb2, such that baseline levels of Hmgb2 are necessary to promote rRNA tran-scription (in concert with other associated factors); however, an exorbitant amount of Hmgb2 overloads rDNA genes and promotes nonspecific chromatin condensation.
To test this, we partially digested chromatin from NRVMs with micrococcal nuclease, isolated heterochromatic, euchromatic, and intermediately compacted DNA based on the level of digestion, and used qPCR to amplify a region of the rDNA gene (designated H42.1 (35)) ( Fig. 5D; see also Fig. 7A for a schematic representation of this assay). We normalized the distribution to the level of the DNA sequence in the most hetero-chromatic fraction, which accounted for the majority of the cardiac genome. We then compared how this ratio changed with hmgb2 knockdown or overexpression and found that hmgb2 overexpression had minimal effect, whereas phenylephrine treatment increased the ratio of euchromatic to heterochromatic DNA (Fig. 5D). hmgb2 knockdown caused an increase in the ratio of intermediately compacted DNA as compared with heterochromatin ( Fig. 5D), illustrating that hmgb2 knockdown and overexpression have opposing effects on rDNA genes.  It was not previously known that Hmgb2 regulates nucleolar transcription. We thus sought to further explore the specific phenotypes regulated by Hmgb2 acting on rRNA transcription. We hypothesized that the dramatic decrease in rRNA synthesis upon hmgb2 overexpression would disrupt cell growth. We labeled cells with crystal violet and then gently washed the plate to remove dead cells. In 3T3 cells, hmgb2 overexpression was lethal; as predicted, however, this was not the case in HeLa cells, where hmgb2 knockdown but not overexpression resulted in a loss of cells (Fig. 5E). ctcf knockdown was also lethal in HeLa cells (Fig. 5E). In NRVMs, hmgb2 knockdown, hmgb2 overexpression, ctcf knockdown, or ctcf overexpression all had no effect on cell death, possibly because the cells are not dividing, and can therefore better withstand disrupted rRNA synthesis (Fig. 5E). This observation is in agreement with data from COS-1 cells, in which Hmgb2 knockdown suppresses cell division (46). Furthermore, the cell type-specific lethality of Ctcf levels is also in agreement with data showing that Ctcf knockout is embryonic lethal in mice (47,48). Ctcf depletion in isolated cells often affects cell division or cell death processes, but in a cell type-dependent manner, potentially due to the cell type-specific localization of Ctcf and arrangement of a higher order structure coordinated by the protein (49).
Interaction of Ctcf and Hmgb2 with Local Chromatin Features to Influence Gene Expression-We previously found that hmgb2-regulated genes are enriched in pathways important to cardiac hypertrophy (21). We next asked how Hmgb2 and Ctcf regulate mRNA expression to explain why some hmgb2-regulated genes are up-regulated and others are down-regulated in response to hmgb2 knockdown in cardiomyocytes (21). We used ENCODE date from the adult mouse heart and our Hmgb2 ChIP-seq data. Hmgb2 occupancy was grossly similar between genes up-regulated or down-regulated by hmgb2 knockdown, as was Ctcf enrichment at ctcf-regulated genes (genes up-regulated or down-regulated by ctcf knock-out in the mouse heart; Fig. 6A). Compared with all genes and down-regulated genes, genes that were up-regulated by hmgb2 knockdown had greater levels of the activating marks H3K4me3 and RNA pol II in their promoters in the basal setting (1 kb upstream of TSS; Fig. 6B), suggesting that removing hmgb2 potentiates the transcriptional effect of the local epigenomic code specified through histone post-translational modifications. Remarkably, ctcf-regulated genes showed the opposite pattern, with up-regulated genes depleted in activating marks (Fig. 6B). This suggested that the interaction of these two pro-  Nucleoli area represents area of individual nucleoli, and nucleolar area represents total area of all nucleoli in a nucleus.) However, we found a decrease in the abundance of nucleolin levels (* indicates p Ͻ 0.05 (Mann-Whitney); one representative experiment of three). C, Hmgb2 ChIP-seq reads were aligned across ribosomal RNA genes showing Hmgb2 is enriched at these loci in the basal state and 48 h after phenylephrine treatment. D, partial chromatin digestion by micrococcal nuclease was used to isolate euchromatic and heterochromatic DNA in NRVMs followed by qPCR to determine the relative distribution of H42.1, a region of rDNA (see Fig. 7A for schematic). The ratio of intermediately packed chromatin to heterochromatin (center of plot) or euchromatin to heterochromatin (right end of plot) was plotted as a change in the ratio after hmgb2 knockdown or overexpression as compared with basal cells. In the control setting, the majority of H42.1 sequences was in the most heterochromatic fraction. hmgb2 overexpression had little effect on the ratios of heterochromatic and euchromatic rDNA, whereas hmgb2 knockdown increased the ratio of intermediately packed to tightly packed DNA. 48 h after phenylephrine, rDNA was shifted to a more euchromatic environment (average of three separate experiments). E, cell viability was assayed by labeling with crystal violet, which stains living cells. Although the effect of overexpression of hmgb2 was minimal in NRVMs, in dividing 3T3 cells there was a dramatic cell killing effect. In contrast, ctcf or hmgb2 knockdown induced cell death in HeLa cells (representative example of three separate experiments).
teins with local chromatin environment could have opposing effects on gene expression.
We next asked if the local chromatin environment at the shared HMGB2 and CTCF loci correlated with which of the two proteins bound. Ctcf binding has been shown to be sensitive to DNA methylation (50 -52) specifically in the context of cell type-dependent binding events (53). We used our bisulfite sequencing data from adult mouse hearts (54) to assess the average DNA methylation status of all CpGs across Hmgb2and Ctcf-binding sites. Bisulfite sequencing data exhibit a bimodal distribution, with the majority of CpGs either lowly methylated or highly methylated, reflective, we posit, of a relative entrainment of CpG methylation status across the population of sampled cells (i.e. methylated or not). As compared with Hmgb2-only or Ctcf-only intergenic peaks, the shared intergenic please were less likely to be highly methylated and were more likely to have intermediate levels of methylation (Fig. 6C). This is also true when examining methylation data from isoproterenol-treated mouse hearts. We thus speculate that DNA methylation at these peaks might be more dynamic and thus capable of modulating Hmgb2 and Ctcf differential binding.
Others have shown that depleting DNA methylation promotes new Ctcf binding at less than 1.5% of Ctcf's consensus motifs (55), suggesting DNA methylation may be insufficient to regulate differential binding of Ctcf to Hmgb2 and the Ctcfshared peaks. This previous study also found that the consensus motifs, which did recruit Ctcf as a result of depleted methylation, were sites that would otherwise have bound Ctcf in other cell types and were at newly formed DNase HS sites (55). We thus asked whether the Hmgb2-and Ctcf-shared peaks also had similar features, specifically the ability to bind Ctcf in other cells and an open chromatin phenotype. We combined ENCODE A, Hmgb2 reads or Ctcf reads were plotted across genes up-regulated or down-regulated by hmgb2 knockdown or ctcf knock-out (left and right panels, respectively) or across all genes in the genome, revealing that the distribution of Hmgb2/Ctcf is not the major determinant of whether a gene's expression will change with depletion of these proteins (Hmgb2-and Ctcf-regulated genes use expression data averaged from three experiments). B, distribution of ChIP-seq (ENCODE, basal adult mouse heart) reads across hmgb2-regulated genes suggests genes that are up-regulated by knockdown are enriched in the activating marks H3K4me3 and RNA pol II in the basal state as compared with genes that are down-regulated by knockdown (top panel). Strikingly, the opposite is true for genes that are regulated by ctcf (bottom panel; ENCODE data represents one biological replicate). C, we next examined DNA methylation from the mouse heart. Plotted is the bimodal distribution of average DNA methylation across Hmgb2-only and Ctcf-only peaks in intergenic regions. Shared peaks are less likely to have high methylation but rather show low or intermediate level methylation. This is true using DNA methylation data from basal adult mice and mice treated with isoproterenol for 3 weeks (DNA methylation is average of six mice (54)).
ChIP-seq data for Ctcf from 10 adult mouse organs and We next asked whether the local chromatin environment at the shared Hmgb2 and Ctcf binding loci correlated with either of the two proteins. s to define a set of conserved Ctcf-binding sites. Only 15% of intergenic Ctcf-only peaks in the heart are also conserved Ctcf-binding sites in other organs. By contrast, 64% of Ctcf peaks shared with Hmgb2 overlapped with conserved Ctcf-binding sites (data not shown). On average, Ctcfonly and Hmgb2-and Ctcf-shared peaks both overlap with one DNase HS site in the mouse heart (ENCODE). These findings suggest that other chromatin features may be responsible for the differential binding of Hmgb2 and Ctcf and that DNA methylation could be one such feature. A complementary interpretation is that Hmgb2's ability to promote heterochromatin and decrease accessibility (as evidenced by our MNase assay, Figs. 5 and 7) could prevent Ctcf binding. However, future analyses are necessary to determine the temporal regulation of other chromatin features that influence the differential binding of Hmgb2 or Ctcf to their shared sites in the heart. Hmgb2 and Ctcf Exert Opposing Effects on Promoter Accessibility-We next sought to investigate direct effects of Hmgb2 and Ctcf to structurally modify local chromatin environment. We isolated nuclei from NRVMs and treated them with 0.001 units of MNase, an enzyme that preferentially digests DNA that is accessible, i.e. not bound by nucleosomes. The digested genome of control NRVMs, when run on an agarose gel, gives bands of multiple sizes as follows: the smallest bands migrate around 150 -200 bp (the size of a mono-nucleosome); bands migrating at increasing molecular masses correspond to sections of chromatin which, endogenously, reside in states of greater compaction (Fig. 7, A, concept, and B, data). Electrophoretically less mobile regions at the top of the gel represent DNA that is more compact and heterochromatic; those at the bottom are more open and euchromatic. hmgb2 knockdown shifted the distribution of the genome toward more euchromatic DNA (Fig. 7B), whereas ctcf knockdown or phenylephrine treatment had minimal effect on the global pattern of DNA compaction (Fig. 7B).
To investigate the behavior of individual genes, we repeated this experiment, and we cut the agarose gel containing the digested genome into three sections representing DNA that had been in heterochromatic environments (compact), euchromatic environments (open), or that came from an intermediate region of the genome, and we analyzed by qPCR the distribution of specific promoter sequences for genes of interest. Each plot represents the change in the ratio of intermediate regions or open regions divided by compact regions between basal and treated cells. For all promoters tested (with the exception of that for the gene dhrs7c), the alteration of local accessibility conferred by phenylephrine was very similar to the effect of hmgb2 knockdown (Fig. 7, C and D), i.e. phenylephrine and hmgb2 knockdown made the chromatin more open (although in some cases the effect of hmgb2 knockdown was less pronounced), whereas hmgb2 overexpression had the antithetic effect. This was true even for genes that had distinct transcriptional responses to phenylephrine treatment (56) and hmgb2 knockdown (Fig. 7D) (21). This observation suggests a partial transcription-independent effect on chromatin structure that is shared in both hmgb2 knockdown-induced and phenylephrine-induced hypertrophy, and which we hypothesize can account for the fact that hmgb2 knockdown can prime the genome for hypertrophy in NRVMs by facilitating the actions of other local chromatin features (Fig. 6). Similar to hmgb2 overexpression, ctcf knockdown generally resulted in more heterochromatic packing of the promoter (Fig. 7E), providing chromatin fiber level support for the opposing relationship between these two proteins described thus far in this study. We hypothesize that Ctcf forms chromatin boundaries, with Hmgb2 promoting compaction within these boundaries (Fig. 7F).
Hmgb2 Binding Overlaps Ctcf at Three-dimensional Domain Boundaries and Cardiac Enhancers to Oppositely Regulate mRNA Expression-We next investigated the role of Ctcf at the sites that both Ctcf and Hmgb2 can bind. Previous work has shown Ctcf is enriched within 100 kb of lamina-associated domains (LADs, 0.1-10 Mb domains of DNA interacting with the nuclear membrane and associated with repressed transcription (57)). Furthermore, analysis of LADs in four different mouse cell types revealed 33% are shared between cell types (58,59). We determined the overlap for these four datasets and compared their location to the proximity of Hmgb2 peaks that do not overlap with Ctcf (Hmgb2-only), Ctcf peaks that do not overlap with Hmgb2 (Ctcf-only), and sites than can be bound by either Hmgb2 or Ctcf (Hmgb2-and Ctcf-shared peaks). This analysis confirmed enrichment of Ctcf-only peaks at LAD boundaries and further revealed Hmgb2-only peaks are enriched at boundaries, suggesting possible new functionality for Hmgb2. However, Hmgb2-and Ctcf-shared peaks were not enriched to the same extent, suggesting that the overlapping peaks correspond to a different functional role of Ctcf (Fig. 8, A  and B).
Topologically associating domains (TADs, megabase-scale domains (median in mESC ϭ 880 kb) of DNA enriched for self-interaction) have Ctcf enriched at their boundaries (15). Furthermore, TAD boundaries are conserved between cell types (ϳ69% shared between cell types of the same species and ϳ65% shared between the same cell type in mouse and human) (15). We compared our peaks to TAD boundaries measured in the mouse cortex and surprisingly found both Hmgb2-only and Hmgb2-and Ctcf-shared peaks were enriched at the boundary (Fig. 8C). We further compared Rad21 (a known component of the cohesin complex) mouse ChIP-seq data with our Ctcf peaks to assess the level of enrichment of Ctcf-and cohesin-shared peaks. As expected, these sites were enriched at TAD boundaries, and the Hmgb2-and Ctcf-shared peaks were enriched to almost the same extent (Fig. 8C). These findings highlight novel functionalities for Hmgb2 (Fig. 8A) to regulate different types of heterochromatin as well as endogenous genomic architecture.
Ctcf has also been shown to act as an enhancer blocker, although recent DNA conformation capture data suggest that whether Ctcf promotes or blocks enhancer interaction with its target gene can depend on whether the interaction occurs within or between TADs (60). We explored whether Ctcf-binding sites regulating enhancer function overlapped with Hmgb2 binding, using cardiac enhancers defined by H3K27ac in adult mouse cardiomyocytes (13). Both Ctcf-only and Hmgb2-only cells; as shown in the 1st panel, these ratios are represented in the ensuing graphs as upward inflections for a shift to more accessible chromatin by a given gene and downward deflections for a gene that shifts to more compact DNA. Gene expressions for hmgb2 knockdown and phenylephrine from microarray data in NRVMs (21) were used to distinguish between genes with similar expression changes induced by hmgb2 and phenylephrine (C) or genes regulated differently by hmgb2 and phenylephrine (D). Interestingly, hmgb2 knockdown and phenylephrine showed similar trends for shifting promoter sequences between these categories even when these stimuli (hmgb2 knockdown or phenylephrine) had different effects on the transcription of the gene. E, promoter sequences for the genes bound by Ctcf and Hmgb2 by ChIP were also examined. In four of the five cases, ctcf knockdown shifted these sequences to more compact regions of chromatin. (All MNase data are the average of at least three experiments.) F, model for relationship between Hmgb2 and Ctcf. Ctcf serves as a boundary preventing the spread of heterochromatin, whereas Hmgb2 promotes heterochromatin formation.

FIGURE 8. Hmgb2 and Ctcf target shared loci near three-dimensional domain boundaries and cardiac enhancers.
A, Ctcf is known to be enriched at LAD and TAD boundaries and to serve as a blocker of enhancer promoter interactions. We show Hmgb2-and Ctcf-shared binding sites are also enriched at TAD boundaries and near cardiac enhancers, with opposing effects on the transcription of genes near these enhancers. B, Hmgb2 peaks that do not overlap Ctcf (Hmgb2-only, red), Ctcf peaks that do not overlap Hmgb2 (Ctcf-only, green), and Hmgb2-and Ctcf-shared peaks (blue) were examined to determine their distance to the nearest conserved mouse lamina-associated domain (LAD, conserved across four mouse cell types). Density plots display the distribution of distances, revealing enrichment of both Ctcf and Hmgb2 at LAD boundaries, with overlapping peaks less enriched (highest peak is further from site of LAD, and peak is shorter and more broad). Signal at 0 bp represents peaks that fall within LADs. C, density plots display the distribution of distances between Hmgb2-only, Ctcf-only, and shared peaks to topologically associating domains (TAD, domain boundaries from mouse cortex but we see similar results with boundaries from mESC). All three sets of peaks are enriched at TAD boundaries (left panel). Cohesin (a complex including Rad21) is known to physically interact with Ctcf at TAD boundaries. Hmgb2-and Ctcf-shared peaks show similar levels of enrichment at TAD boundaries as the levels exhibited by Rad21-and Ctcf-shared peaks (right panel). D, Hmgb2-only and Ctcf-only peaks are enriched near cardiac enhancers (defined by H3K27ac from mouse adult cardiomyocytes), with Hmgb2-and Ctcf-shared peaks enriched to a greater degree (left panel). The nearest genes to cardiac enhancers were found and grouped by the distance of the enhancer to a Ctcf peak (green) or shared peak (blue). Close (C) indicates peaks within 30 bp of enhancer, medium (M) indicates 31-381 bp distance, and far (F) indicates 381 bp to 1 kb. Plotted is the percentage of these nearest genes up-or down-regulated by ctcf KO (middle panel) or hmgb2 KD (left panel) at each distance. When Hmgb2-and Ctcf-shared peaks are within 380 bp to 1 kb of an enhancer, loss of ctcf is three times more likely to cause the nearest gene to be up-regulated than down-regulated, whereas the loss of hmgb2 is four times more likely to cause down-regulation. E, for Ctcf-only, Hmgb2-only, and shared peaks that occur within 1 kb of a gene, but not in a gene, we plotted the percentages that were near active genes (P, presence of RNA pol II peak in promoter) or inactive genes (NP, no pol II). Only the shared peaks showed bias toward being near active genes (left panel). When Hmgb2-only, Ctcf-only, or shared peaks are within 1 kb of an active gene, loss of ctcf is two times more likely to cause up-regulation (dark color) than down-regulation (light color), but no bias is seen for inactive genes or peaks within genes. Rad21 and RNA pol II ChIP-seq represent one biological replication from ENCODE. Location of cardiac enhancers (13), LAD boundaries (59), and TAD boundaries (15) come from published work which cite one, two, and multiple (number not provided) biological replicates, respectively.
peaks are enriched at these enhancers, with Hmgb2-and Ctcfshared peaks enriched the most (Fig. 8D). This supports a role for these ubiquitous proteins to regulate cell type-specific gene expression. We then asked for Hmgb2 or Ctcf peaks close to an enhancer (within 30 bp), at a medium distance (31-380 bp) or at a far distance (381 bp to 1 kb), what the effect was on the expression of the nearest gene to the enhancer using RNA-seq data from ctcf knock-out mouse hearts and microarray data from hmgb2 knockdown NRVMs. Loss of ctcf from regions a medium or far distance from an enhancer is biased toward upregulation of the target gene (although the majority of genes do not change in expression, Fig. 8D), mimicking the expected behavior of an enhancer blocker. The same is seen when ctcf is lost from sites that can also bind Hmgb2. At Hmgb2-and Ctcfshared peaks, loss of ctcf is three times more likely to result in up-regulation than down-regulation of the nearest gene to the enhancer, whereas loss of hmgb2 is four times more likely to cause down-regulation than up-regulation (Fig. 8D). One possible explanation is that loss of hmgb2 allows Ctcf binding, which promotes enhancer blocking.
We then performed a similar analysis, this time asking how Hmgb2 or Ctcf peaks near genes that regulate the expression of the nearest gene. First, we found that Hmgb2 and Ctcf are just as likely to bind near (designated as within 1 kb of a gene but not within a gene) an active gene as an inactive one, when using the presence of an RNA pol II promoter peak to define active genes (ENCODE ChIP-seq dataset). However, Hmgb2-and Ctcfshared peaks are biased to binding near the active genes (Fig.  8E, left panel). Second, we found that for Hmgb2-only, Ctcfonly, and Hmgb2-and Ctcf-shared peaks that occur near a gene, loss of ctcf was biased toward up-regulation of the gene if the gene was also active in the basal state (Fig. 8E, right panel). However, there was no bias for genes that did not have RNA pol II nor for genes where Hmgb2 or Ctcf were binding within the gene as opposed to near the gene (Fig. 8A). This finding is similar to our observation that hmgb2 target genes are likely to be up-regulated if bound by RNA pol II in their promoter in the basal state and also is in line with our observations from enhancers, wherein Ctcf binding near, but not in, enhancers or genes is biased to being repressive, a function that is conserved at sites that can also bind Hmgb2.
It remains unknown whether and how Hmgb2 binding can modulate the repressive function of Ctcf on nearby genes. We hypothesize this action to involve mutually exclusive binding of these two proteins, wherein Hmgb2's binding prevents that of Ctcf by changing the chromatin landscape and/or specifically promoting heterochromatin (Fig. 7F).
Additionally, gene ontology analysis of genes with a nearby (within 1 kb) shared Hmgb2 and Ctcf peak are enriched for annotation relating to cytoskeleton, ribosome, and nucleus (enrichment score 4.14, 3.51, and 3.25 respectively, DAVID), whereas gene ontology analysis for the nearest genes to enhancers with Hmgb2-and Ctcf-shared peaks were enriched for transcription (enrichment score 3.62, DAVID), with no enrichment for any disease processes after Bonferroni correction (hypertrophic cardiomyopathy and insulin signaling were the two lowest p values for KEGG pathways, although not significant after correction). These data suggest that although Hmgb2 and Ctcf target cardiac-specific genes and can regulate pathologic pathways (21), their shared binding sites may be more important for regulating gross, rather than stress-responsive, cardiac genome organization.

Discussion
We propose a model (Figs. 7F and 8A) whereby Ctcf, acting in an insulator capacity, serves as a boundary for heterochromatin spreading. In the absence of Ctcf, heterochromatin can spread, silencing nearby regions, a phenomenon that is accompanied by the increased presence of Hmgb2, which maintains the compact environment. Inversely, increased abundance of Hmgb2 can promote heterochromatin spreading and thereby evict Ctcf from specific loci. By disrupting the boundaries of heterochromatin, Hmgb2 and Ctcf can affect multiple genes in a given region, depending on the type and extent of other modifications.
We favor this model, as opposed to one in which Hmgb2 preferentially targets (and differentially regulates) individual genes, because Hmgb2 lacks DNA sequence specificity and is not a cardiac-specific protein. These observations, coupled with our new understanding of the finely regulated balance between Hmgb2 and Ctcf, indicate that although the overall chromatin structure of the myocyte may not be regulated with single gene resolution, this structure is critical for regulating myocyte physiology in health and disease. It is also important to note that the conclusions in this study were made from a combination of in vivo adult mouse models and isolated neonatal ventricular myocytes. Although this approach provides some experimental advantages, there are likely important differences in how adult and neonatal cardiomyocytes package chromatin (related to differences in regenerative, proliferative, and stress response capacities in these cell populations), which will have to be resolved by further experimentation.
An open question when we began these studies was the molecular basis for how Hmgb2 ostensibly promotes the transcription of some genes while inhibiting the expression of others. Previous studies have implicated Hmgb2 in transcriptional activation or repression (61)(62)(63), attributing these actions to cooperation with distinct proteins. The predominant hypothesis for the role of Hmgbs in gene expression is that they bend DNA to promote binding by other proteins. Hmgbs may only transiently interact with these client proteins (if at all), thereby acting as promiscuous chaperones at different loci (28). Our data support this model, indicating that Hmgb2 has locus-specific effects on gene expression notwithstanding a conserved effect to compact chromatin in a locus-independent manner. Future studies using locus-specific proteomic analyses will be required to determine which Hmgb2 binding partners encode, in a combinatorial manner, different transcriptional logic. Although informative, this approach may obfuscate the issue of trans effects, particularly through non-coding regions of the genome, by arbitrarily restricting examination of Hmgb2's functions within the physical unit of a gene.
The increased DNA flexibility conferred by Hmgb2 binding (64) that facilitates formation of locus-specific complexes can also more generally facilitate tighter packaging of DNA. Previous reports have shown Hmgb1 is enriched in euchromatic regions of photoreceptor nuclei (65). However, Hmgb1 and Hmgb2 are also bound to highly heterochromatic DNA formed during mitosis (66). Compared with linker histone H1, Hmgb1 also compacts DNA, although to a lesser degree (67), and can directly compete with histone H1 for binding to linker DNA (68). We show that hmgb2 knockdown disrupts global measures of heterochromatin (H3K27me3 abundance) and the chromatin environment at specific genes without altering global levels of histone H1. One potential explanation is that hmgb2 knockdown decreases nucleosome abundance. Others have shown that Hmgb1-deficient mouse embryonic fibroblast cells have reduced nucleosome number, and yeast deficient in the hmgb2 homologue display a shift to euchromatin when assayed by MNase (69), similar to what we observed following hmgb2 knockdown. Indeed, in vitro, Hmgb1 can facilitate nucleosome deposition (70). However, we see no difference in the abundance of histone H3 with Hmgb2 knockdown, although we did not directly measure nucleosome assembly. Hmgb2 can also alter nucleosome distribution by facilitating nucleosome sliding via interaction with SWI/SNF ATP-dependent chromatin remodeling complexes, as shown in vitro (71).
Here, we propose a model where Hmgb2 targeting is partially regulated by the distribution of heterochromatin, such that Ctcf mediates the boundaries between hetero-and euchromatin, and Hmgb2 maintains the integrity of facultative heterochromatin, i.e. genes that are silenced in a given cell type. We reason that the overlap between Hmgb2 ChIP-seq reads from rat cardiomyocytes with Ctcf ChIP-seq peaks in other species and tissues indicates cell type independent functions of the proteins. Indeed, topological domains are largely conserved between cell types and species (15), and Ctcf is both enriched at the topological domains (15) and is critical for maintaining them (72). Unlike transcription factors, Ctcf does not preferentially localize to genes that belong to a similar class (73). However, we see preferential binding of Hmgb2-and Ctcf-shared peaks to genes actively expressed in the heart, suggesting that a subset of the Hmgb2 and Ctcf binding is cell type-specific. Furthermore, the changes to chromatin accessibility at individual promoters induced by hmgb2 knockdown largely mimic the effects of phenylephrine treatment. This is especially interesting given that phenylephrine treatment, unlike hmgb2 knockdown, does not cause global changes in genome accessibility, suggesting that alterations in chromatin packaging, although important for phenotype, may be decoupled from transcriptional changes.
We also find both Ctcf and Hmgb2 regulate nucleolar transcription in multiple cell types, having particular implications for cardiac cells. Nucleolar disruption occurs with cardiac stress (34), and rRNA synthesis is up-regulated in hypertrophy (74). Previous findings identified Hmgb3 as a component of the T-cell nucleolar proteome (75) and found evidence for interaction between Hmgb2 and nucleolin outside of the nucleolus (76).
The relationship between Hmgb2 and Ctcf suggests a mechanism by which both chromatin structural proteins are regulated (in abundance and localization) in part by the chromatin environment. Ctcf, unlike Hmgb2, has DNA binding consensus motifs, and perhaps their coregulation involves sequestering of Ctcf-binding sites into heterochromatin by Hmgb2 and/or is mediated through changes in DNA methylation, although this will take additional experiments to fully elucidate. Our data also indicate that these two proteins confer opposite regulation of chromatin accessibility when they target the same promoters and opposite regulation of gene expression when they target the same cardiac enhancers.
In this model, Ctcf organizes the framework of the genome within which the cell type-specific chromatin factors operate. Hmgb2 also acts within the boundaries of this model to maintain heterochromatic regions (with a high density of Hmgb2 to allow for tight packaging) and facilitate complex formation (with a low density of Hmgb2 priming DNA for binding by other proteins) whose specific functions are dependent on the cell type-specific proteome. In support of this model, we observe a promoter-specific effect of hmgb2 knockdown on transcription but a uniform effect of Hmgb2 to regulate chromatin accessibility at the genomic scale. Thus, in a cell type-dependent way, the nucleus regulates the regions established by Ctcf. However, alterations to balance between Hmgb2 and Ctcf disrupt the boundaries of heterochromatin, undoing the cell type-specific silencing. We hypothesize that the changes in the ratio of Hmgb2 to Ctcf that we observe with cardiac pathology and across genetic backgrounds allow for varied genomic plasticity, supporting a general theory in which global chromatin accessibility is an important component to transcriptome remodeling in disease.

Experimental Procedures
Analysis of the Hybrid Mouse Diversity Panel-Microarray (RNA isolated from the left ventricle) and phenotypic data from 84 classical inbred and recombinant strains of adult (aged 8 -10 weeks) female mice in the basal state or after treatment with isoproterenol (Alzet microosmotic pump releasing 30 mg/kg/ day) for 3 weeks (29) were analyzed. These data (29) were acquired from female mice for three reasons. First, choosing a single gender removed this as a variable and the scope of this study prohibited repeating all the analyses in both genders. Second, male mice are prone to establishing societal hierarchies when housed in the same cage, which could affect how the mice responded to stress signaling, such as isoproterenol. Third, pilot studies on both genders revealed a more reproducible and pronounced phenotype in female animals. Transcript abundances were correlated to determine an R 2 value, which was converted to a p value. Correlation of transcript abundance in liver (31) and bone marrow (32) was also assessed. Separately, the ratio of ctcf to hmgb2 (calculated by subtracting the logscaled hmgb2 value from the ctcf value) was plotted against cardiac phenotype, with different plots for different subsets of strains. Strains were grouped by whether they went into failure upon isoproterenol treatment (n ϭ 13 strains), developed cardiac hypertrophy (n ϭ 22), showed minimal change (n ϭ 9), or showed phenotypic traits inconsistent with a single disease state (n ϭ 40).
We used Genomic Ranges in R to determine the portion of Hmgb2 peaks that were overlapped by Ctcf peaks and designated them Hmgb2-and Ctcf-shared peaks. We then removed these from the other datasets to get Hmgb2-only and Ctcf-only peaks. We compared the proximity of these peaks to other genomic features (lamina-associated domains (57), topologically associating domains (15), and cardiac enhancers (13)) using BED Tools to determine the distance between two features and density plots to visualize the distribution of distances of the Hmgb2 or Ctcf peaks to a genomic feature (set to bp 0). We defined peaks as being close distance (in or within 30 bp), medium distance (31-380 bp), or far distance (381 bp to 1 kb) of a cardiac enhancer and then found the nearest gene (nearest function in Genomic Ranges) and asked how the nearest gene's expression changed in response to ctcf knock-out or hmgb2 knockdown. Bar graphs indicate the percent of genes up-or down-regulated at each distance (with the remaining genes unchanged; not plotted). We also found Hmgb2 or Ctcf peaks within 1 kb of, but not in, a gene and plotted the percent of nearest genes that were active or inactive (determined by the presence of a RNA pol II promoter peak, ENCODE, ENCFF001YAH). For peaks within a gene, or within 1 kb of a gene, we asked how the expression of the gene changed with ctcf knock-out depending on whether the gene was active or inactive in the basal state. Percentage indicates percent of genes up-or down-regulated out of the pool of all genes within that distance of a peak and with the same pol II binding status (promoter peak or not). Ctcf cell type-independent peaks were determined by finding overlapping peaks in 10 adult mouse tissues, including the heart, and mESC cells (ENCODE: ENCFF001YAF, ENCFF001YAC, ENCFF001-YAW, ENCFF001YAY, ENCFF001YBA, ENCFF001YAM, ENCFF001YAI, ENCFF001XZU, ENCFF001YBC, ENCFF001-YAO, and ENCFF001XZY). 9 R. Garrido and T. M. Vondriska, unpublished data.
Hmgb2 ChIP-seq and microarray data were compared with Hi-C chromatin conformation capture data (15) to identify genes regulated by Hmgb2 that fell in the same topological domain. Genes regulated by Hmgb2 in rat were found in the genomes of the Hi-C data sets (mm9, hg18). For domains with more than one gene, each gene was compared with every other gene to determine whether they had matching transcriptional responses to hmgb2 knockdown (either both up-regulated or both down-regulated). As control, all RefSeq genes were randomly sampled to create a mock list of genes of the same size as our Hmgb2-regulated list, which were then randomly designated as up or down-regulated, and used to determine the percentage of inter-domain comparisons that matched. This was repeated for a total of 70 times. In all cases, the average percentage that matched plateaued (ϳ50%) after 5-30 random samples.
Genes regulated by Hmgb2 (microarray) or Ctcf (RNA-seq) were found in the mouse genome and compared with ChIP-seq data in the adult mouse heart from ENCODE datasets as follows: RNA pol II, ENCFF001LKL; H3K4me3, ENCFF001KHV. Alignment of ChIP-seq data across Hmgb2-regulated genes was performed using SeqPlots with the following parameters: anchored features, 10-bp bins, extend targets 1 kb up and downstream. As control, alignment across all genes was determined using RefSeq genes and gene predictions for mm9 downloaded from the UCSC Genome Browser.
For ChIP-reChIP experiments (35) the manufacturer's guidelines for the ChIP kit were followed until elution. Elution was performed in elution buffer shaken at 65°C for 30 min. The second immunoprecipitation was performed on eluant diluted in IP dilution buffer. ChIP-reChIP experiments used anti-Ctcf (Abcam, ab70303) and anti-Hmgb2 (Abcam, ab67282) antibody except when the second immunoprecipitation was for the same protein as the first, in which case the second immunoprecipitation used anti-Ctcf (Active Motif, 61311) or anti-Hmgb2 (Abcam, ab55169). Primers target the promoter of the indicated genes; negative control was determined from Hmgb2 ChIP-seq. See under "Primers" for the primer sequences.
Microscopy-Images were acquired on a Nikon A1R confocal microscope and analyzed in ImageJ. For colocalization analysis, super-resolution was achieved using dual-color STED microscopy on a custom STED instrument developed at UCLA. Colocalization was measured using in-house software to measure distances between clusters.