PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of narLink to Publisher's site
 
Nucleic Acids Res. May 2011; 39(10): 4076–4087.
Published online Feb 2, 2011. doi:  10.1093/nar/gkq1343
PMCID: PMC3105402
Effects of HMGN variants on the cellular transcription profile
Mark Rochman,corresponding author1 Leila Taher,corresponding author2 Toshihiro Kurahashi,1 Srujana Cherukuri,1 Vladimir N. Uversky,3,4,5 David Landsman,2 Ivan Ovcharenko,2 and Michael Bustin1*
1Protein Section, Laboratory of Metabolism, Center for Cancer Research, National Cancer Institute, 2Computational Biology Branch, National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, 9000 Rockville Pike, Bethesda, MD 20894, 3Institute for Intrinsically Disordered Protein Research, Center for Computational Biology and Bioinformatics, Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, IN 46202, 4Department of Molecular Medicine, University of South Florida, Tampa, FL 33612, USA and 5Institute for Biological Instrumentation, Russian Academy of Sciences, 142290 Pushchino, Moscow Region, Russia
corresponding authorCorresponding author.
*To whom correspondence should be addressed. Tel: Phone: +1 301 496 5234; Fax: +1 301 496 8419; Email: bustin/at/helix.nih.gov
The authors wish it to be known that, in their opinion, the first two authors should be regarded as joint First Authors.
Present address: Srujana Cherukuri, University of Maryland School of Medicine, Department of Pediatrics, Center for Stem Cell and Regenerative Medicine 20 Penn Street, HSF II Room S103, Baltimore, MD 212014, USA.
Received December 6, 2011; Accepted December 22, 2010.
High mobility group N (HMGN) is a family of intrinsically disordered nuclear proteins that bind to nucleosomes, alters the structure of chromatin and affects transcription. A major unresolved question is the extent of functional specificity, or redundancy, between the various members of the HMGN protein family. Here, we analyze the transcriptional profile of cells in which the expression of various HMGN proteins has been either deleted or doubled. We find that both up- and downregulation of HMGN expression altered the cellular transcription profile. Most, but not all of the changes were variant specific, suggesting limited redundancy in transcriptional regulation. Analysis of point and swap HMGN mutants revealed that the transcriptional specificity is determined by a unique combination of a functional nucleosome-binding domain and C-terminal domain. Doubling the amount of HMGN had a significantly larger effect on the transcription profile than total deletion, suggesting that the intrinsically disordered structure of HMGN proteins plays an important role in their function. The results reveal an HMGN-variant-specific effect on the fidelity of the cellular transcription profile, indicating that functionally the various HMGN subtypes are not fully redundant.
The dynamic architecture of the chromatin fiber plays a key role in regulating transcriptional processes necessary for proper cell function and mounting adequate responses to various internal and external biological signals. Architectural nucleosome-binding proteins such as the linker histone H1 protein family and the high mobility group (HMG) protein superfamily are known to continuously and reversibly bind to chromatin, transiently altering its structure and affecting the cellular transcription output (1,2). Although extensively studied, the cellular function and mechanism of action of these chromatin-binding architectural proteins are still not fully understood. A major question in this field is the extent of the functional specificity of the structural variants of histone H1 or of the various HMG families (3–6). Experiments with genetically altered mice lacking one or several H1 variants revealed that loss of one variant leads to increase synthesis of the remaining variants, suggesting functional redundancy between H1 variants (7,8). Yet, analysis of cells in which the levels of specific H1 variants have been altered suggests a certain degree of variant-specific effects on transcriptional output (9–11)
The HMG superfamily is composed of three families named HMGA, HMGB and high mobility group N (HMGN), each containing several protein members (3,4). It is known that HMG proteins affect transcription and modulate the cellular phenotype (12); however, the transcriptional specificity of the various HMG variants has not yet been systematically studied. Here, we examine the role of the various HMGN variants in the regulation of the cellular transcription profile.
The HMGN family of chromatin architectural proteins consists of five members with a similar structure (13). All contain a bipartite nuclear localization signal (NLS), a highly conserved nucleosome-binding domain (NBD) and a negatively charged and highly disordered C-terminal domain. The HMGNs are the only nuclear proteins known to specifically recognize generic structural features of the 147-bp nucleosome core particle (CP), the building block of the chromatin fiber (3,4). HMGN binds to chromatin and CP without any known specificity for the sequence of the underlying DNA. In the nucleus, HMGNs are highly mobile moving among nucleosomes in a stop-and-go manner (2,14). The fraction of time that an HMGN resides on a nucleosome (stop period) is longer than the time it takes to ‘hop’ from one nucleosome to another; therefore, most of the time, most of the HMGNs are bound to chromatin. The amount of HMGN present in most nuclei is sufficient to bind only ~1% of the nucleosomes; however, the dynamic binding of HMGNs to chromatin ensures that potentially every nucleosome will temporarily interact with an HMGN molecule. Thus, potentially, HMGNs may affect the transcription of numerous genes.
HMGN variants share several functional properties, such as binding affinity to nucleosomes in vitro and in vivo, competition with linker histone H1 for the binding sites on nucleosomes, and effects on chromatin architecture. Likewise, both HMGN1 and HMGN2, the most abundant and ubiquitous members of this protein family, form multiple complexes with nuclear proteins (15). These findings, and the similarity of their domain structure, suggest that by enlarge, HMGN proteins could be functionally redundant. Yet, several studies indicate that HMGN proteins are not fully redundant. Both in vivo and in vitro studies indicate that the interaction of HMGN variants with CPs lead to the formation of complexes containing two molecules of a single type of variant; CPs containing two different HMGN variants are not formed under physiological conditions (16,17). In addition, while HMGN1 and HMGN2 seem to be ubiquitously expressed, HMGN3 and HMGN5 proteins show distinct developmental and tissue-specific expression (18–20). Most significantly, analysis of genetically altered mice and cells revealed variant-specific phenotypes and indication that the variants are not fully functionally redundant (12).
It has been repeatedly shown that interaction of HMGNs with chromatin affects transcription (21–24). However, the extent of specificity of HMGN variants in transcriptional regulation and the level of functional redundancy between them remain largely unknown, mainly because of the lack of systematic analysis of the effect of HMGNs on gene expression in a unified experimental system.
To gain insights into the extent of transcriptional specificity of the HMGN variants, we compared expression profiles of mouse embryonic fibroblasts (MEFs) in which various HMGN variants were either knocked out or stably overexpressed, to double their cellular content. We found that loss of proteins affected the expression of a limited number of genes, while doubling the cellular levels of an HMGN variant affected the expression of hundreds of genes. While some of the genes were affected by more than one variant, the great majority of the genes were affected in a variant-specific manner.
Intrinsically disordered proteins are predicted to affect transcription even at low dosage overexpression because they form weak interactions with multiple partners (25,26). Thus, the significant transcriptional effects resulting from doubling the amount of HMGNs is in agreement with the highly intrinsically disordered structure of HMGN proteins and with their tendency to form multiple metastable protein complexes (15). We also found that specific variants affect the transcription profile in a cell-specific manner. Analysis of domain swap mutants suggests that the specificity of each HMGN is determined by a unique combination of a functional NBD and a C-terminal domain.
The results reveal an HMGN-variant-specific effect on the global transcription profile suggesting that these proteins fine tune the fidelity of the cellular transcription. We speculate that part of their specificity is due to their intrinsic highly disordered structure that enables each variant to form multiple types of complexes with nuclear components.
Isolation of MEFs and generation of stable cell lines
Mouse embryonic fibroblasts SV-40 transformed (MEFs) were purchased from ATCC.
MIN6 cell line was a gift from A.L.Notkins, NIDCR, NIH. Primary MEFs from variant-specific knock out mice were isolated from two embryos as described (27) and analyzed separately. Cells were grown in Dulbecco's Modified Eagle Medium (DMEM) supplemented with 10% Fetal Calf Serum (FCS).
Retroviruses were produced in Phoenix helper cell line transfected with the pHAN vector bearing various HMGN proteins tagged with FLAG and HA at C-terminus. Stable cell lines were generated by retroviral infection in the presence of polybrene at the concentration of 5 µg/ml and subsequent selection with puromycin at the concentration of 1 µg/ml for 7 days. Cells were grown without antibiotics for 1 day prior to collecting samples for expression analysis and western blotting.
Antibodies and western blotting
All the antibodies used in the study were from our laboratory. Secondary Horseradish Peroxidase (HRP)-Conjugated antibodies were from Pierce.
Whole cell lysates were prepared in 2× Laemmli sample buffer (Bio-Rad) supplemented with protease inhibitors. Samples were separated on 15% pre-cast Criterion gels, transferred by semi-dry method to polyvinylidene difluoride (PVDF) membrane, blocked with non-fat milk in Phosphate Buffered Saline (PBS) and probed with indicated antibodies. Chemiluminiscent detection by enhanced chemiluminescence (ECL) Detection Reagent (Amersham) was done according to the manufacturer’s recommendations.
RNA preparation
RNA was prepared by TRIzol reagent according to the manufacturer’s protocol. Subsequently, RNA was cleaned up by Qiagen RNeasy kit with on-column DNaseI treatment.
Gene arrays
Microarray expression analysis was performed using Affymetrix Mouse GeneChips 430 2 (430v2). Hybridization of biotin-labeled cRNA fragment to Mouse Genome 430 2.0 array, washing, staining with streptavidin–phycoerythrin (Molecular Probes), and signal amplification were performed according to the manufacturer’s instructions at the Laboratory of Molecular Technology (LMT, Frederic, NCI).
Statistical analysis
We analyzed 51 array data sets (n = 3 for each particular experiment) to search for genes whose expression levels were significantly altered. All analyses were performed using R and BioConductor (28). R packages ‘affy’ (29) and ‘simpleaffy’ (30) and ‘affyQCReport’ (31) were employed to evaluate the quality of the arrays by means of images, histograms, box plots, degradation plots and scatter plots. Expression values were derived using the Robust Multichip Average protocol (32) with default settings. All analyses were done at the so-called sequence level, i.e. data from probes representing the same gene were combined. We did not apply any unspecific filter on the expression values.
Differentially expressed genes were identified using an empirical Bayes method implemented in the R package ‘Limma’ (33). P-values were corrected for multiple testing using a false discovery rate method (34). Genes for which the adjusted P-value was <0.001 (overexpression of different HMGN variants) or <0.05 (knock out of different HMGN variants) in at least one of the comparisons were considered differentially expressed. No fold-change cut-off was applied.
Mouse Genome 430 2.0 array has 45 101 probe sets associated with approximately 20 000 Mouse Genome Informatics (MGI) gene identifiers. Probe sets were mapped to MGI identifiers using information provided by the Jackson Laboratory (http://www.informatics.jax.org/).
Venn diagrams
Differentially expressed genes in all experiments were compared to controls and represented as Venn diagrams (R package ‘Vennerable’) (35).
Functional analysis
Functional analysis of microarray data was based on overrepresentation of GO terms (36). P-values were corrected for multiple comparisons using Bonferroni’s method.
Bioinformatics structural analysis
Composition profiling
Analysis of amino acid composition of HMGN proteins was performed using Composition Profiler online service (http://www.cprofiler.org) (37) with default settings. The following reference protein sets were used: DisProt 3.4 (38), PDB Select 25 (39) and mouse HMGN proteins. The set DisProt 3.4 comprises consensus sequences of experimentally determined disordered regions; PDB Select 25 contains PDB structures with <25% sequence identity, biased toward the composition of proteins amenable to crystallization studies. Amino acids are arranged in the order of increase of their disorder propensity, according to the scale by Radivojac et al. (40).
Intrinsic disorder prediction
Per-residue predictions of intrinsic disorder in HMGN proteins were performed using a PONDR® VLXT predictor, access to which was provided by Molecular Kinetics, Inc. (http://www.pondr.com). PONDR® (Predictor Of Natural Disordered Regions) is a set of neural network predictors of disordered regions on the basis of local amino acid composition, flexibility, hydropathy, coordination number and other factors. These predictors classify each residue within a sequence as either ordered or disordered. PONDR® VL-XT integrates three feed forward neural networks: the Variously characterized Long, version 1 (VL1) predictor (41), which predicts non-terminal residues, and the X-ray characterized N- and C-terminal predictors (XT) (42), which predicts terminal residues. Output for the VL1 predictor starts and ends 11 amino acids from the termini. The XT predictor output provides predictions up to 14 amino acids from their respective ends. A simple average is taken for the overlapping predictions; a sliding window of nine amino acids is used to smooth the prediction values along the length of the sequence. Unsmoothed prediction values from the XT predictors are used for the first and last four sequence positions.
Structural characterization of HMGN variants
Examination of the structure of genes coding for the various members of the HMGN family suggests that they originated from a common ancestor. All the genes contain relatively long 5′- and 3′-untranslated regions, six exons and the boundaries of the first four exons are highly conserved (Figure 1A). The gene coding for HMGN5 evolved recently because it is found only in mammals. All the proteins encoded by the genes contain a positively charged, highly conserved, NBD (Figure 1A) that serves as their main chromatin-binding site. Embedded in the NBD is the sequence RRSARLSA(K,M)P that has been shown to be the core sequence that specifically anchors HMGN proteins to the 147-bp nucleosome CP, the building block of the chromatin fiber (43). A NLS that is localized at the N-terminal part of the proteins is also highly conserved in all HMGN variants. The C-terminal region of the proteins, encoded by exons 5 and 6, differs significantly among the HMGN variants. The HMGN5 C-terminal domain is especially long and contains several repeats of a negatively charged sequence motif (18). The alignment shown in Figure 1A illustrates the major similarities and differences between the mouse HMGN variants. Mouse HMGN1, HMGN2 and HMGN3a are similar in size, ranging between 89 and 95 amino acids, and are more similar to each other than to HMGN5, which is 406 amino acid long. The alignment does not contain the splice variant HMGN3b, which lacks the 21 C-terminal residues of HMGN3a, nor the HMGN4 variant, which has not yet been investigated in detail.
Figure 1.
Figure 1.
HMGN proteins are intrinsically disordered. (A) Multiple sequence alignment of mouse HMGN1, HMGN2, HMGN3a and HMGN5 proteins by ClustalW. Only the first 94 amino acids of HMGN5 are aligned. The positively charged NBD, the hallmark of HMGN proteins, is (more ...)
Analysis of the amino acid composition of the HMGN proteins in comparison to ordered proteins listed in the Protein Data Bank (http://www.pdb.org/pdb/home/home.do) reveals that potentially all HMGNs are highly disordered proteins (Figure 1B). in fact, HMGNs are expected to be significantly more disordered than an ‘average’ disordered protein, because they are much more depleted in major order-promoting residues (compare the colored bars with negative values for various HMGNs with the black bars for Intrinsically Disordered Proteins at the left side of the plot) and are significantly enriched in major disorder-promoting residues (right side of the plot). Interestingly, the HMGN variants are different from each other and show significant variability in amino acid compositions, as exemplified by the large variations in R, T, D, G, A, S, E and P.
In agreement with intrinsic disorder prediction of HMGN proteins based on amino acid composition, PONDR analysis (41,44) predicts a high degree of structural disorder in all HMGN variants, with a few short regions with increased order propensity (Figure 1C, dips in the graph). These relatively ordered regions often correspond to potential binding sites that fold upon interaction with binding partners (45–47). Notably, all HMGNs contain several regions that according to the ANCHOR algorithm (48) are predicted to serve as binding sites for other interacting proteins (Figure 1D), an observation fully compatible with our previous findings that both HMGN1 and HMGN2 can be found in numerous metastable multiprotein complexes (15). The number and localization of the predicted protein binding sites although highly similar are not identical between HMGN variants.
In summary, although all HMGNs share several physical properties and are nuclear proteins that bind to nucleosome CPs through a highly conserved domain, each variant has a distinct structure and has several sites for interacting with other proteins. These characteristics raise the possibility of HMGN-variant-specific effects on the cellular transcription profile.
Transcriptional impact of HMGN proteins
To investigate the transcriptional specificity of HMGN variants we first analyzed the transcriptional profile of primary MEFs isolated form Hmgn1/, Hmgn3/ and Hmgn5/ mice using mouse 430.2 Affymetrix expression arrays. Hmgn2/ are not available because these mice are embryonic lethal (M.B. unpublished data). The results reveal variant-specific changes in gene expression profile; no overlap was observed for the genes affected by the knockout of different HMGN variants (Figure 2A). The changes involved both up- and downregulation of transcript levels, a finding that is fully compatible with the notion that HMGNs enhance transcriptional fidelity by affecting chromatin structure and optimizing the fidelity of transcription. Even though transcription of a relatively small number of genes was affected, Gene Ontology (GO) analysis revealed significant enrichment in a few non-overlapping pathways for Hmgn1/, Hmgn3/ MEFs (Figure 2B). These results are in agreement with our previous observations that the phenotypes of Hmgn1/ and Hmgn3/ mice are distinct but not severe (19,49).
Figure 2.
Figure 2.
Effects of HMGNs knock out on transcription in primary MEFs. (A) Venn diagrams of down- and upregulated genes in primary MEFs. (B) GO analysis of affected genes (P <0.05).
Because HMGN proteins are intrinsically disordered proteins (Figure 1) and because dosage changes in such proteins may lead to large changes in transcription (25), we reasoned that a mild increase in the cellular levels of HMGN variants, in a uniform system, may give a more sensitive indication of the potential transcriptional specificity of the HMGN variants. To compare the effect of the overexpression of HMGN variants on transcription in a uniform system, we used retroviral infection to generate MEFs cell lines stably expressing specific HMGN variants tagged with FLAG and HA at their C-terminus. Vectors expressing HMGN1, HMGN2, HMGN3a, HMGN5 and the HMGN5-S17,21E double-point mutant, which does not bind to chromatin (18), were generated and efficiently expressed in MEFs. Following infection, cells were subjected to the selection procedure and all the cells that passed the selection were analyzed as a pool.
Western blot analysis of the infected cells revealed that the level of expression of each exogenous protein was comparable to the level of its endogenous counterpart (Figure 3A). Thus, stably infected MEFs express ~2-fold higher levels of a specific HMGN variant. The HMGN5-S17,21E protein contains mutations in two serine residues in the NBD which abolish its binding to nucleosomes (18) and thus served as a control for transcriptional effects due to chromatin binding. As an additional control, we infected MEFs with virus carrying an empty vector. For each variant, we analyzed the transcription profile of three independently infected pools of MEFs using mouse 430.2 Affymetrix expression arrays. We compared the transcription profile of MEFs overexpressing specific HMGN variants to the control cell lines transfected with empty vectors or with the HMGN5-S17,21E double-point mutant.
Figure 3.
Figure 3.
Effects of elevated expression of HMGNs on transcription in MEFs. (A) Western blot analysis of stably infected MEFs. Shown are western analysis of MEFs stably expressing FLAG and HA tagged HMGN variants. Endogenous and exogenous HMGN proteins are indicated. (more ...)
Class comparison between cell lines indicated that overexpression of HMGNs altered the expression level of 5203 genes. Three-dimensional clustering of these transcripts, based on principal component analysis (PCA), revealed that the various cell lines formed four distinct expression clusters (Figure 3B). Three of the clusters were formed by the cell overexpressing either HMGN1, or HMGN2, or HMGN5. The fourth cluster was formed by cell lines overexpressing either the HMGN5-S17,21E double-point mutant, cell transfected with an empty vector, or by the cells overexpressing the HMGN3a variant. These results demonstrate that in MEFs, each HMGN variant has specific effects on transcription.
The cells differed not only in the specificity of genes affected but also in the number of genes affected. Doubling the levels of HMGN1, HMGN2 and HMGN5 affected the expression of 1268, 2753 and 3183 genes, whereas HMGN3a and HMGN5-S17,21E caused no significant changes in gene expression. For HMGN1, HMGN2 and HMGN5 proteins, the proportion of up- and down regulated genes was roughly equal, indicating that the proteins did not preferentially activate or inhibit transcription (Figure 3C).
More detailed comparison of the genes affected by each of the HMGN variants revealed that while each protein either up- or downregulated the expression of a unique set of genes, a fraction of the genes was affected by more than one HMGN protein, suggesting partial redundancy in transcriptional regulation (Figure 3D). Thus, of the 457 genes that were downregulated by overexpressing HMGN1, 40% were uniquely affected, 17% were also downregulated by HMGN5, 22% were also downregulated by HMGN2 and 21% were downregulated by all three HMGNs. Of the 811 genes that were upregulated in HMGN1-overexpressing cells, 44% were uniquely affected, 12% were also affected by HMGN5, 31% were also affected by HMGN2 and 13% were upregulated by all the HMGNs. For HMGN2, a total of 1263 genes were downregulated, of these 60% were specifically downregulated only by HMGN2, and 24% were also downregulated by HMGN5. Likewise, ~50% of the 1490 genes were specifically upregulated by HMGN2 and ~70% of the genes up- or downregulated by HMGN5 were specifically affected by HMGN5.
Only a small proportion of the 5203 genes whose transcription changed by overexpression of the HMGNs was commonly affected by all the three HMGN proteins. In all, 96 genes (2%) were downregulated and 103 (2%) were upregulated. The most extensive overlap was observed between HMGN1 and HMGN2 proteins; 44% of the genes upregulated by HMGN1 were also upregulated by HMGN2. We also found a large number of genes regulated by both HMGN5 and HMGN2; 478 and 403 genes were up- and downregulated by both proteins, respectively. While the total number of the genes whose transcription levels changes was statistically significant was relatively large, the transcription levels of most of the genes changed ~2-fold (Figure 3E). These findings are in agreement with previous studies indicating that while HMGN proteins affect the expression of many genes, the changes in transcription levels are relatively small (19,27,50).
Next, we performed a functional analysis on the sets of genes exclusively regulated by each individual protein, as well as on the sets of genes regulated by combinations of several proteins (Table 1) for significantly overrepresented GO terms. The results indicate that overexpression of individual HMGN variants affected gene from different categories. Whereas HMGN1 affected genes involved in cell division and mitosis, HMGN2 regulated genes involved in regulation of transcription, development and chromatin binding. Notably, genes involved in cell cycle regulation were also preferentially affected in Hmgn1/ MEFs (Figure 2B). HMGN5-induced transcriptional changes were mainly associated with metabolic processes, protein, metal ion and transcription factor binding. We note, however, that the GO analysis suggests a certain degree of redundancy among the HMGN variants. For instance, GO term ‘response to virus’ (GO:0009615) was enriched for genes commonly regulated by HMGN1 and HMGN2 proteins. In addition, several biosynthetic processes, such as sterol biosynthesis (GO:0016126), cholesterol (GO:0006695) and lipid biosynthesis (GO:0008610), and others were enriched for the genes regulated by all three HMGN variants. In fact, of 199 genes regulated by all HMGNs, 52 genes are involved in various biosynthetic processes.
Table 1.
Table 1.
GO analysis of gene expression in MEFs with elevated HMGN levels
Taken together, the gene expression profiles and the GO analyses reveal a surprising degree of specificity in the effects of the various HMGN variants on the cellular transcription profile. The functional redundancy among the variants is lower than what would be expected from a set of proteins with structural similarities, which bind to nucleosomes with similar affinities, have highly similar NBDs and use an identical sequence motif to bind specifically to nucleosome CPs.
Transcriptional effects of HMGN swap mutants
Because the N-terminal half of the HMGNs are highly similar, while their C-terminal domains are clearly distinct (Figure 1A), we assumed that the functional specificity of the proteins resides in their C-terminal region. To test this assumption, we generated retroviral vectors expressing tail swap mutants in which C-terminal region of either HMGN2 or HMGN3a protein was fused to N-terminal part of HMGN1 protein, immediately after the conserved NBD (see Figure 1A for exact location of the regions swapped). We named these mutants as N1–N2 swap and N1–N3 swap. The correct expression of the swap mutants in MEFs infected with the retroviral vectors was verified by western blot analysis (Figure 4A) using an antibody elicited against the conserved NBD of the HMGN protein family, which recognizes all the HMGN variants (51).
Figure 4.
Figure 4.
Comparison of the effect of HMGN tail swap mutants on transcription in MEFs. (A) Western blot analysis of MEFs stably expressed swap mutants using an antibody that recognizes the conserved NBD. N1, endogenous HMGN1; c, control expression of empty vector. (more ...)
The transcription profile of the MEFs expressing the swap mutants was determined by mouse 430.2 Affymetrix expression arrays and compared with that of cells infected with vectors expressing the native proteins. Three-dimensional clustering of the results using PCA (Figure 4B) revealed that the effect of the swap mutants on transcription was distinct from that of their ‘source’ proteins. Thus, while HMGN3a did not affect transcription (Figure 3) in MEFs, the N1–N3 swap mutant significantly affected the transcription of 1522 genes, most of which were distinct from the genes affected by HMGN1 (Figure 4C, 3–4). Further comparison of the genes affected by the swap mutants with the genes affected by either HMGN1 or HMGN2 (Figure 4C) supported the notion that transcriptional changes induced by the tail swap mutants differed from those observed for HMGN1, HMGN2 or HMGN3a proteins. The swap mutants specifically downregulated 621 genes (Figure 4C 1,3) and upregulated 402 genes (Figure 4C 2,4). Interestingly, the two swap mutants had very similar effects on the cellular transcription profile (Figure 4D). Of the 1252 and the 1177 genes, respectively, downregulated by the N1–N2 and the N1–N3 swap mutants, close to 80% genes overlapped. Likewise, most of the genes that were up-regulated by the N1–N2 swap mutant were also upregulated by the N1–N3 swap protein (Figure 4D). The similarity in the genes regulated by the swap mutants points out to the importance of their shared NBD region in determining the effect on the transcription profile. Yet, the effects were clearly distinct from HMGN1, the donor of their shared NBD, an indication that ultimately, the structure of the entire protein, the combination of individual N- and C-terminal domains, determines the functional specificity of the HMGN variants.
Cell type-specific effects on transcription
Surprisingly, our analysis revealed that overexpression of HMGN3a had no effect on transcription profile of the transfected MEFs. We previously reported that in MEFs the protein levels of HMGN3 are lower than those of HMGN1 and HMGN2 that are robustly expressed in most cells. However, HMGN3a is highly expressed in MIN6 cell, a mouse pancreatic cell line that secretes insulin (19). In these cells, small interfering RNA-mediated downregulation of HMGN3, but not that of HMGN1 or HMGN2, affects the transcription of genes involved in insulin secretion, suggesting that HMGN3a may affect transcription in a cell type-specific manner. To test this possibility, we first re-examined the relative amount of HMGN3 protein in MIN6 and MEF cells. Western blot analyses revealed that indeed the HMGN3 protein levels in MIN6 were significantly higher than in MEFs (Figure 5A). Next, we infected MIN6 cells with retroviral vectors expressing either HMGN1 or HMGN3a proteins and verified protein expression by western blot analysis (Figure 5B and C).Transcriptional array analysis revealed that in MIN6 cells, HMGN3a significantly changed the expression of 1429 genes; of these 471 were up-regulated and 958 genes were down-regulated (Figure 5D). In contrast, overexpression of HMGN1 in MIN6 cell line had no significant effect on the transcription profile. Because HMGN1 had significant effects on the transcription profile of MEFs (Figure 3), these results suggest cell type-specific transcription effects of HMGN variants.
Figure 5.
Figure 5.
Effect of HMGN1 and HMGN3a on transcription in MIN6 cells. (A) Comparison of the protein levels of HMGN3 and HMGN1 in MIN6 and MEFs by western blot. CBB, Coomassie Blue staining. Western blot analysis of MIN6 cell lines stably expressing HMGN1 (B) or (more ...)
The major goal of this study is to examine whether the various members of the HMGN protein family can affect the cellular transcription profile in an HMGN-variant-specific manner. Although previous studies indicated that the binding of HMGN protein to chromatin alters the cellular transcription profile, the degree to which these changes are HMGN-variant specific has not yet been investigated.
The dynamic nature of HMGN binding to chromatin and the lack of any DNA sequence specificity in their chromatin interactions, taken together with the conservation of their nuclear-binding domain and similarities in their overall organization and physical properties, raised the possibility that the individual HMGN variants would be functionally redundant and have similar effects on the cellular transcription profile. Conversely, the widespread expression of HMGN1 and HMGN2, but not HMGN3 and HMGN5, in most tissues and the sequence specificity of their C-terminal domains suggest that potentially the proteins may have variant-specific effects on the transcription profile. Indeed, in vitro experiments revealed variant-specific effects on histone modifications, and experiments with genetically altered mice also suggest that the HMGN variants are not fully functionally redundant.
Our experiments suggest that each HMGN variant can affect the expression of numerous genes, especially when overexpressed, and by enlarge in an HMGN-specific manner. The amplitude of transcriptional changes was moderate; for most of the affected genes being in the limits of 2-fold difference. Importantly, nearly equal amount of genes were either up- or downregulated by each HMGN, suggesting that HMGNs are neither transcriptional activators nor repressors. The GO analyses indicated that multiple cellular processes were affected by individual HMGNs or by combinations of several HMGNs, suggesting that the HMGNs are general modulators of the cellular transcriptional fidelity.
Two molecular mechanisms whereby HMGN affect the transcription profile could be envisioned. One possibility is that by binding to nucleosomes, HMGN induce structural changes that alter the ability of transcriptional regulators, either positive or negative, to interact with their chromatin targets. A second possibility is that the HMGN interact with specific regulators and affect their chromatin interactions. Both possibilities suggest that the ability of HNGN variants to bind to chromatin is a major effect on the transcriptional output. Indeed, our previous experiments (50), and our present analyses of the HMGN5-S19,23E mutant, indicate that HMGNs affect transcription by binding to nucleosomes.
While nucleosome binding seems to be an absolute requirement for any noticeable effects on transcription, the variant-specific effects on the transcription profile suggest that additional properties of these proteins play a role in determining their biological specificity. Because the C-terminal domain of HMGNs is highly variable in sequence between individual HMGN proteins, we tested the possibility that the specific transcriptional effects of HMGNs reside in this domain and expressed several HMGN swap mutants in MEF cells. Surprisingly, the transcriptional outcome following the expression of these swap mutants with a common NBD from HMGN1, and a C-terminal domain from HMGN1, HMGN2 or HMGN3 was different from either one of their ‘source’ proteins. Thus, the variant-specific effects of HMGNs on transcription are the consequence of coordinate effects of the various structural domains of each variant. In other words, neither the NBD nor the C-terminal domain alone defines the transcriptional effect of each HMGN protein, but rather the entire structure of the protein defines its specific role in transcription (Figure 5E).
In considering the molecular mechanisms leading to HMGN-variant-specific effects on transcription, we note that early structural studies indicated that HMGNs have little ordered structure (52), and our computational analysis (Figure 1) reveal that HMGNs are among the most intrinsically disordered proteins known. Intrinsically disordered proteins can interact with multiple protein partners with relatively low affinity and acquire more ordered structures (45–47,53–58). It has been recently reported that the harmful effect of elevated cellular levels of many proteins is correlated with the degree of their disorderness (25). At the same time, cells with decreased amount of these proteins function robustly and do not demonstrate significant changes in cellular functions. Our observation that knock out of HMGNs has significantly smaller effect on transcription supports this theory and strongly argues that disordered structure of HMGNs is one of the major functional properties of these proteins.
Variations in structure of HMGN proteins due to interaction with different protein partners can modulate the effects of HMGN variants on local nucleosome structure, global chromatin architecture and transcription (Figure 5E). Indeed, both HMGN1 and HMGN2 have been shown to form multiple metastable macromolecular complexes (15), and specific protein partners have been identified for several HMGN variants. Thus, HMGN3 interacts specifically with the thyroid hormone receptor (59) and with the transcription factor PDX1 (19), HMGN1 forms a complex with ERalpha and SRF (15), and HMGN2 was shown to interact with PITX2 (60). Our observation of cell-specific effects of HMGN3a protein on transcription in the pancreatic derived MIN6 cell line, but not in MEFs (19), supports the idea of existence of specific protein partners for individual HMGN proteins.
In conclusion, our results reveal both specific and redundant roles of HMGN variants in the global regulation of gene expression. Each HMGN preferentially affects a unique set of genes with little or no specificity for defined cellular processes. Thus, changes in the expression of an HMGN may disrupt the fidelity of the cellular transcription and render the organism more susceptible to further damage. Indeed, experiments with genetically altered mice and with cells derived from these mice indicate that loss of HMGN1 leads to an impaired DNA damage repair response and increased tumorigenicity (27,49,61). Likewise, loss of HMGN3, which is highly expressed in beta cells of the pancreatic islets, affects insulin secretion leading to a mild diabetic phenotype (19). The transcriptional specificity of the HMGN variants is similar to that of the H1 variants. It seems that the dynamic interaction of HMGN, H1 and other structural proteins with chromatin is part of the mechanism that ultimately fine tunes the transcription profile to optimize cellular function.
FUNDING
Intramural program of the CCR, NCI, NIH and Research Program of the National Library of Medicine, NIH. Funding for open access charge: Intramural, CCR, NCI, NIH.
Conflict of interest statement. None declared.
ACKNOWLEDGEMENT
We thank Dr Nina Bubunenko, LMT, Frederick for help with the Affymetrix arrays.
1. Bustin M, Catez F, Lim JH. The dynamics of histone H1 function in chromatin. Mol. Cell. 2005;17:617–620. [PubMed]
2. Catez F, Hock R. Binding and interplay of HMG proteins on chromatin: lessons from live cell imaging. Biochim. Biophys. Acta-Gene Regul. Mech. 2010;1799:15–27. [PubMed]
3. Bianchi ME, Agresti A. HMG proteins: dynamic players in gene regulation and differentiation. Curr. Opin. Genet. Dev. 2005;15:496–506. [PubMed]
4. Bustin M. Regulation of DNA-dependent activities by the functional motifs of the high-mobility-group chromosomal proteins. Mol. Cell. Biol. 1999;19:5237–5246. [PMC free article] [PubMed]
5. Woodcock CL, Skoultchi AI, Fan Y. Role of linker histone in chromatin structure and function: H1 stoichiometry and nucleosome repeat length. Chromosome Res. 2006;14:17–25. [PubMed]
6. Khochbin S. Histone H1 diversity: bridging regulatory signals to linker histone function. Gene. 2001;271:1–12. [PubMed]
7. Fan Y, Skoultchi AI. Genetic analysis of H1 linker histone subtypes and their functions in mice. Methods Enzymol. 2004;377:85–107. [PubMed]
8. Fan Y, Sirotkin A, Russell RG, Ayala J, Skoultchi AI. Individual somatic H1 subtypes are dispensable for mouse development even in mice lacking the H1(0) replacement subtype. Mol. Cell. Biol. 2001;21:7933–7943. [PMC free article] [PubMed]
9. Alami R, Fan Y, Pack S, Sonbuchner TM, Besse A, Lin Q, Greally JM, Skoultchi AI, Bouhassira EE. Mammalian linker-histone subtypes differentially affect gene expression in vivo. Proc. Natl Acad. Sci. USA. 2003;100:5920–5925. [PubMed]
10. Bhan S, May W, Warren SL, Sittman DB. Global gene expression analysis reveals specific and redundant roles for H1 variants, H1c and H1(0), in gene expression regulation. Gene. 2008;414:10–18. [PMC free article] [PubMed]
11. Sancho M, Diani E, Beato M, Jordan A. Depletion of human histone H1 variants uncovers specific roles in gene expression and cell growth. PLoS Genet. 2008;4:e1000227. [PMC free article] [PubMed]
12. Hock R, Furusawa T, Ueda T, Bustin M. HMG chromosomal proteins in development and disease. Trends Cell Biol. 2007;17:72–79. [PMC free article] [PubMed]
13. Postnikov Y, Bustin M. Regulation of chromatin structure and function by HMGN proteins. Biochim. Biophys. Acta. 1799:62–68. [PMC free article] [PubMed]
14. Catez F, Brown DT, Misteli T, Bustin M. Competition between histone H1 and HMGN proteins for chromatin binding sites. EMBO Rep. 2002;3:760–766. [PubMed]
15. Lim JH, Bustin M, Ogryzko VV, Postnikov YV. Metastable macromolecular complexes containing high mobility group nucleosome-binding chromosomal proteins in HeLa nuclei. J. Biol. Chem. 2002;277:20774–20782. [PubMed]
16. Postnikov YV, Trieschmann L, Rickers A, Bustin M. Homodimers of chromosomal proteins HMG-14 and HMG-17 in nucleosome cores. J. Mol. Biol. 1995;252:423–432. [PubMed]
17. Postnikov YV, Herrera JE, Hock R, Scheer U, Bustin M. Clusters of nucleosomes containing chromosomal protein HMG-17 in chromatin. J. Mol. Biol. 1997;274:454–465. [PubMed]
18. Shirakawa H, Landsman D, Postnikov YV, Bustin M. NBP-45, a novel nucleosomal binding protein with a tissue-specific and developmentally regulated expression. J. Biol. Chem. 2000;275:6368–6374. [PubMed]
19. Ueda T, Furusawa T, Kurahashi T, Tessarollo L, Bustin M. The nucleosome binding protein HMGN3 modulates the transcription profile of pancreatic beta cells and affects insulin secretion. Mol. Cell. Biol. 2009;29:5264–5276. [PMC free article] [PubMed]
20. West KL, Ito Y, Birger Y, Postnikov Y, Shirakawa H, Bustin M. HMGN3a and HMGN3b, two protein isoforms with a tissue-specific expression pattern, expand the cellular repertoire of nucleosome-binding proteins. J. Biol. Chem. 2001;276:25959–25969. [PubMed]
21. Ding HF, Rimsky S, Batson SC, Bustin M, Hansen U. Stimulation of RNA polymerase II elongation by chromosomal protein HMG-14. Science. 1994;265:796–799. [PubMed]
22. Paranjape SM, Krumm A, Kadonaga JT. HMG17 is a chromatin-specific transcriptional coactivator that increases the efficiency of transcription initiation. Genes Dev. 1995;9:1978–1991. [PubMed]
23. Bustin M. Chromatin unfolding and activation by HMGN(*) chromosomal proteins. Trends Biochem. Sci. 2001;26:431–437. [PubMed]
24. Rochman M, Malicet C, Bustin M. HMGN5/NSBP1: a new member of the HMGN protein family that affects chromatin structure and function. Biochim. Biophys. Acta. 1799:86–92. [PMC free article] [PubMed]
25. Vavouri T, Semple JI, Garcia-Verdugo R, Lehner B. Intrinsic protein disorder and interaction promiscuity are widely associated with dosage sensitivity. Cell. 2009;138:198–208. [PubMed]
26. Hansen JC, Lu X, Ross ED, Woody RW. Intrinsic protein disorder, amino acid composition, and histone terminal domains. J. Biol. Chem. 2006;281:1853–1856. [PubMed]
27. Birger Y, West KL, Postnikov YV, Lim JH, Furusawa T, Wagner JP, Laufer CS, Kraemer KH, Bustin M. Chromosomal protein HMGN1 enhances the rate of DNA repair in chromatin. EMBO J. 2003;22:1665–1675. [PubMed]
28. Gentleman R, Carey V, Bates D, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. [PMC free article] [PubMed]
29. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307–315. [PubMed]
30. Wilson C, Miller C. Simpleaffy: a BioConductor package for Affymetrix Quality Control and data analysis. Bioinformatics. 2005;21:3683–3685. [PubMed]
31. Parman C, Halling C, Gentleman R. affyQCReport: QC Report Generation for affyBatch objects. R package version 1.18.0.
32. Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–193. [PubMed]
33. Smyth GK. Limma: linear models for microarray data. In: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W, editors. Bioinformatics and Computational Biology Solutions using R and Bioconductor. New York: Springer; 2005. pp. 397–420.
34. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. Roy. Stat. Soc. Ser. B (Methodological) 1995;57:289–300.
35. Swinton J. Vennerable: Venn and Euler area-proportional diagrams. R package version 2.0/r75. 2009
36. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000;25:25–29. [PMC free article] [PubMed]
37. Vacic V, Uversky VN, Dunker AK, Lonardi S. Composition Profiler: a tool for discovery and visualization of amino acid composition differences. BMC Bioinformatics. 2007;8:211. [PMC free article] [PubMed]
38. Sickmeier M, Hamilton JA, LeGall T, Vacic V, Cortese MS, Tantos A, Szabo B, Tompa P, Chen J, Uversky VN, et al. DisProt: the Database of Disordered Proteins. Nucleic Acids Res. 2007;35:D786–D793. [PubMed]
39. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE. The Protein Data Bank. Nucleic Acids Res. 2000;28:235–242. [PMC free article] [PubMed]
40. Radivojac P, Iakoucheva LM, Oldfield CJ, Obradovic Z, Uversky VN, Dunker AK. Intrinsic disorder and functional proteomics. Biophys. J. 2007;92:1439–1456. [PubMed]
41. Romero P, Obradovic Z, Li X, Garner EC, Brown CJ, Dunker AK. Sequence complexity of disordered protein. Proteins. 2001;42:38–48. [PubMed]
42. Li X, Romero P, Rani M, Dunker AK, Obradovic Z. Predicting protein disorder for N-, C-, and internal regions. Genome Inform. Ser. Workshop Genome Inform. 1999;10:30–40. [PubMed]
43. Ueda T, Catez F, Gerlitz G, Bustin M. Delineation of the protein module that anchors HMGN proteins to nucleosomes in the chromatin of living cells. Mol. Cell. Biol. 2008;28:2872–2883. [PMC free article] [PubMed]
44. Romero P, Obradovic Z, Dunker AK. Natively disordered proteins: functions and predictions. Appl Bioinformatics. 2004;3:105–13. [PubMed]
45. Garner E, Romero P, Dunker AK, Brown C, Obradovic Z. Predicting binding regions within disordered proteins. Genome Inform. Ser. Workshop Genome Inform. 1999;10:41–50. [PubMed]
46. Oldfield CJ, Cheng Y, Cortese MS, Romero P, Uversky VN, Dunker AK. Coupled folding and binding with alpha-helix-forming molecular recognition elements. Biochemistry. 2005;44:12454–12470. [PubMed]
47. Cheng Y, Oldfield CJ, Meng J, Romero P, Uversky VN, Dunker AK. Mining alpha-helix-forming molecular recognition features with cross species sequence alignments. Biochemistry. 2007;46:13468–13477. [PMC free article] [PubMed]
48. Dosztanyi Z, Meszaros B, Simon I. ANCHOR: web server for predicting protein binding regions in disordered proteins. Bioinformatics. 2009;25:2745–2746. [PMC free article] [PubMed]
49. Birger Y, Catez F, Furusawa T, Lim JH, Prymakowska-Bosak M, West KL, Postnikov YV, Haines DC, Bustin M. Increased tumorigenicity and sensitivity to ionizing radiation upon loss of chromosomal protein HMGN1. Cancer Res. 2005;65:6711–6718. [PubMed]
50. Rochman M, Postnikov Y, Correll S, Malicet C, Wincovitch S, Karpova TS, McNally JG, Wu X, Bubunenko NA, Grigoryev S, et al. The interaction of NSBP1/HMGN5 with nucleosomes in euchromatin counteracts linker histone-mediated chromatin compaction and modulates transcription. Mol. Cell. 2009;35:642–656. [PMC free article] [PubMed]
51. Bustin M, Crippa MP, Pash JM. Immunochemical analysis of the exposure of high mobility group protein 14 and 17 surfaces in chromatin. J. Biol. Chem. 1990;265:20077–20080. [PubMed]
52. Bradbury EM. Conformations and flexibilities of histones and high mobility group (HMG) proteins in chromatin structure and function. Ciba Found. Symp. 1983;93:246–270. [PubMed]
53. Tompa P. Intrinsically unstructured proteins. Trends Biochem. Sci. 2002;27:527–533. [PubMed]
54. Dunker AK, Silman I, Uversky VN, Sussman JL. Function and structure of inherently disordered proteins. Curr. Opin. Struct. Biol. 2008;18:756–764. [PubMed]
55. Oldfield CJ, Meng J, Yang JY, Yang MQ, Uversky VN, Dunker AK. Flexible nets: disorder and induced fit in the associations of p53 and 14-3-3 with their partners. BMC Genomics. 2008;9(Suppl. 1):S1. [PMC free article] [PubMed]
56. Tompa P, Fuxreiter M. Fuzzy complexes: polymorphism and structural disorder in protein-protein interactions. Trends Biochem. Sci. 2008;33:2–8. [PubMed]
57. Uversky VN, Dunker AK. Understanding protein non-folding. Biochim. Biophys. Acta. 1804:1231–1264. [PMC free article] [PubMed]
58. Lu X, Hamkalo B, Parseghian MH, Hansen JC. Chromatin condensing functions of the linker histone C-terminal domain are mediated by specific amino acid composition and intrinsic protein disorder. Biochemistry. 2009;48:164–172. [PMC free article] [PubMed]
59. Lee JW, Choi HS, Gyuris J, Brent R, Moore DD. Two classes of proteins dependent on either the presence or absence of thyroid hormone for interaction with the thyroid hormone receptor. Mol. Endocrinol. 1995;9:243–254. [PubMed]
60. Amen M, Espinoza HM, Cox C, Liang X, Wang J, Link TM, Brennan RG, Martin JF, Amendt BA. Chromatin-associated HMG-17 is a major regulator of homeodomain transcription factor activity modulated by Wnt/beta-catenin signaling. Nucleic Acids Res. 2008;36:462–476. [PMC free article] [PubMed]
61. Kim YC, Gerlitz G, Furusawa T, Catez F, Nussenzweig A, Oh KS, Kraemer KH, Shiloh Y, Bustin M. Activation of ATM depends on chromatin interactions occurring before induction of DNA damage. Nat. Cell Biol. 2009;11:92–96. [PMC free article] [PubMed]
Articles from Nucleic Acids Research are provided here courtesy of
Oxford University Press