- Open Access
Heterogeneity in the kinetics of nuclear proteins and trajectories of substructures associated with heterochromatin
© Stixová et al; licensee BioMed Central Ltd. 2011
- Received: 1 December 2010
- Accepted: 18 March 2011
- Published: 18 March 2011
Protein exchange kinetics correlate with the level of chromatin condensation and, in many cases, with the level of transcription. We used fluorescence recovery after photobleaching (FRAP) to analyse the kinetics of 18 proteins and determine the relationships between nuclear arrangement, protein molecular weight, global transcription level, and recovery kinetics. In particular, we studied heterochromatin-specific heterochromatin protein 1β (HP1β) B lymphoma Mo-MLV insertion region 1 (BMI1), and telomeric-repeat binding factor 1 (TRF1) proteins, and nucleolus-related proteins, upstream binding factor (UBF) and RNA polymerase I large subunit (RPA194). We considered whether the trajectories and kinetics of particular proteins change in response to histone hyperacetylation by histone deacetylase (HDAC) inhibitors or after suppression of transcription by actinomycin D.
We show that protein dynamics are influenced by many factors and events, including nuclear pattern and transcription activity. A slower recovery after photobleaching was found when proteins, such as HP1β, BMI1, TRF1, and others accumulated at specific foci. In identical cells, proteins that were evenly dispersed throughout the nucleoplasm recovered more rapidly. Distinct trajectories for HP1β, BMI1, and TRF1 were observed after hyperacetylation or suppression of transcription. The relationship between protein trajectory and transcription level was confirmed for telomeric protein TRF1, but not for HP1β or BMI1 proteins. Moreover, heterogeneity of foci movement was especially observed when we made distinctions between centrally and peripherally positioned foci.
Based on our results, we propose that protein kinetics are likely influenced by several factors, including chromatin condensation, differentiation, local protein density, protein binding efficiency, and nuclear pattern. These factors and events likely cooperate to dictate the mobility of particular proteins.
- Mean Square Displacement
- Fluorescence Recovery After Photobleaching
- Nuclear Periphery
- Fluorescence Recovery
The eukaryotic nucleus is highly compartmentalised. The functional consequences of nuclear compartmentalisation have been described previously in both fixed and live cells [1–3]. Generally, chromatin consists of histones wrapped with DNA and the many proteins that are directly responsible for proper nuclear functions, such as replication, transcription, splicing, and DNA repair. Nuclear processes proceed in specific compartments, such as nuclear foci, transcription and replication factories, or nuclear speckles (summarised in [4, 5]). Moreover, interphase chromosomes are arranged into chromosome territories, which can intermingle to some degree [4, 6, 7]. Centromeres and telomeres are also essential structures for chromosome function. Centromeres are the sites of mitotic spindle attachment and are required for cell division . Telomeres, which include the shelterin proteins (telomeric-repeat binding factor 1 (TRF1), TRF2, protection of telomeres 1 (POT1), tripeptidyl peptidase 1 (TPP1), Ras-related protein 1 (RAP1), and TRF1-interacting nuclear factor 2 (TIN2)), protect the ends of chromosomes and are essential for chromosome stability [9, 10]. Another prominent nuclear structure is the nucleolus, which is the largest transcription factory and is necessary for the synthesis of ribosomal subunits [11–13]. Like the rest of the genome, nucleoli are remarkably compartmentalised. The nucleolar region consists of the fibrillar centre (FC), the dense fibrillar component (DFC), and the granular component (GC). The boundaries of nucleoli are surrounded by clusters of centromeric heterochromatin, called chromocentres . The nucleolar organiser regions (NORs) of specific acrocentric chromosomes are responsible for the structural and functional properties of nucleoli [15, 16]. Many specific proteins preferentially localise to the boundary between the FC and DFC, which is thought to be the site of transcription of ribosomal genes. The transcription machinery used for rRNA synthesis includes RNA polymerase I (RNA pol I) and upstream binding factor (UBF) [17, 18]. Gorski et al. measured the dynamic nature of these important nucleolar proteins using fluorescence recovery after photobleaching (FRAP). In that study, the fluorescence recovery kinetics of several RNA pol I subunits in the G1 and S phases of the cell cycle displayed biphasic behaviour, characterised by fast fluorescence increase followed by a gentle fluorescence recovery phase . These authors also showed that the RNA polymerase I recovery kinetics decrease as cell cycle-dependent transcription increases. The transcription of RNA polymerase I can be suppressed by exposure to a variety of stimuli. Actinomycin D is a widely used drug that intercalates into double stranded DNA and blocks the transcription elongation activity of all three polymerases (I, II, III) [20, 21]. Therefore, actinomycin D is a useful agent that enables the study of gene silencing and corresponding epigenetic events.
Green fluorescent protein (GFP) technologies combined with microscopy have opened a new avenue for studying the mobility of tagged proteins in living cells. However, most methods are limited by the confocal optical resolution of about 200 nm laterally and 600 nm axially . Thus, several microscopy techniques have been developed to improve the microscope resolution to 10 to 20 nm. The new methods include 4Pi microscopy , localisation microscopy approaches , three-dimensional structured illumination microscopy (3D-SIM) , and stimulated emission depletion microscopy (STED) . Photobleaching methods are used to analyse the kinetic properties of particular proteins [1, 26]. Many chromatin-related proteins have been studied using the FRAP technique to measure local protein dynamics, which are important functional characteristics of these proteins. For example, two distinct kinetic pools for Polycomb group-related protein, B lymphoma Mo-MLV insertion region 1 (BMI1), have been described . Similarly, Wang et al. showed the energy-dependent heterogeneous movement of telomeres. In other studies, it has been reported that heterochromatin protein 1 (HP1α, HP1β, HP1γ) accumulates more rapidly into euchromatin-rich nuclear regions than the HP1 subtypes of highly condensed heterochromatin foci [29, 30]. The recovery of HP1α protein is faster in pluripotent mouse embryonic stem cells (mESCs) than in differentiated cells . Furthermore, the mobile fraction of histones H2B and H3 is low compared to HP1 subtypes. However, an increased recovery after photobleaching of H2B-GFP and H3-yellow fluorescent protein (YFP) was measured in mESCs relative to differentiated cells. This implies that the kinetics of chromatin-related proteins are cell-type specific and dependent on the level of chromatin condensation .
Previous studies have provided interesting data for functionally important features of chromatin-related proteins. The observations described above have led to additional questions. Is the association of particular proteins with nuclear domains of specific function the most important criterion responsible for their kinetic properties or trajectories? Could molecular weight or protein binding efficiency also influence protein diffusion into a photobleached region? We analysed the trajectories and kinetics of several chromatin-related proteins and determined the relationships between protein nuclear arrangement, molecular weight, and recovery kinetics. In addition, we asked whether the trajectories and kinetics of particular proteins change in response to histone hyperacetylation by the histone deacetylase (HDAC) inhibitor trichostatin A (TSA) or after suppression of transcription by actinomycin D.
Cellular patterns of proteins studied
Correlation between molecular weight and kinetics of chromatin-related proteins
In contrast, we showed that there was no significant correlation between molecular weight and fluorescence recovery for proteins that accumulated at foci (Figure 2c) or nucleoli (Figure 2d). These results indicate that molecular weight is not the sole factor in determining protein kinetics, but that other factor can influence the dynamics of some biologically active molecules.
Trajectories of HP1β, BMI1, and TRF1 foci are influenced by histone hyperacetylation and suppression of transcription
Area of heterochromatin protein 1β (HP1β), B lymphoma Mo-MLV insertion region 1 (BMI1), and telomeric-repeat binding factor 1 (TRF1) foci
Area of minimal enclosing ellipse (μm2)
Peripheral foci: control
0.22 ± 0.12 (n = 10)
0.17 ± 0.09 (n = 4)
0.19 ± 0.06 (n = 8)
Peripheral foci: TSA
0.13 ± 0.10 (n = 19)
0.02 ± 0.02 (n = 12)
0.7 ± 0.54 (n = 15)
Peripheral foci: actinomycin D
0.59 ± 0.48 (n = 4)
0.09 ± 0.05 (n = 5)
0.08 ± 0.05 (n = 20)
Peripheral foci: vorinostat
0.17 ± 0.10 (n = 23)
0.17 ± 0.22 (n = 13)
0.28 ± 0.47 (n = 29)
Central foci: control
0.15 ± 0.07 (n = 7)
0.11 ± 0.05 (n = 10)
0.27 ± 0.21 (n = 8)
Central foci: TSA
0.16 ± 0.09 (n = 16)
NA (n = 0)
0.65 ± 0.42 (n = 8)
Central foci: actinomycin D
0.69 ± 0.46 (n = 6)
0.07 ± 0.06 (n = 5)
0.05 ± 0.02 (n = 11)
Central foci: vorinostat
0.21 ± 0.12 (n = 13)
0.09 ± 0.03 (n = 10)
0.15 ± 0.17 (n = 18)
Area mean: control
0.19 ± 0.11
0.13 ± 0.07
0.23 ± 0.16
Area mean: TSA
0.14 ± 0.09
0.68 ± 0.49
Area mean: actinomycin D
0.65 ± 0.44
0.08 ± 0.05
0.06 ± 0.04
Area mean: vorinostat
0.18 ± 0.11
0.14 ± 0.17
0.23 ± 0.39
Average velocity of heterochromatin protein 1β (HP1β), B lymphoma Mo-MLV insertion region 1 (BMI1), and telomeric-repeat binding factor 1 (TRF1) foci
Average velocity (μm/s)
Peripheral foci: control
0.034 ± 0.007
0.017 ± 0.003
0.024 ± 0.005
Peripheral foci: TSA
0.024 ± 0.011
0.012 ± 0.005
0.046 ± 0.012
Peripheral foci: actinomycin D
0.024 ± 0.003
0.016 ± 0.005
0.021 ± 0.015
Peripheral foci: vorinostat
0.020 ± 0.003
0.020 ± 0.007
0.026 ± 0.010
Central foci: control
0.036 ± 0.018
0.015 ± 0.003
0.021 ± 0.004
Central foci: TSA
0.028 ± 0.017
0.050 ± 0.020
Central foci: actinomycin D
0.022 ± 0.003
0.013 ± 0.003
0.018 ± 0.003
Central foci: vorinostat
0.020 ± 0.004
0.017 ± 0.003
0.026 ± 0.010
Velocity mean: control
0.035 ± 0.012
0.016 ± 0.003
0.023 ± 0.005
Velocity mean: TSA
0.026 ± 0.014
0.047 ± 0.015
Velocity mean: actinomycin D
0.023 ± 0.003
0.015 ± 0.005
0.020 ± 0.010
Velocity mean: vorinostat
0.020 ± 0.004
0.018 ± 0.005
0.026 ± 0.010
Diffusion coefficient of heterochromatin protein 1β (HP1β), B lymphoma Mo-MLV insertion region 1 (BMI1), and telomeric-repeat binding factor 1 (TRF1) foci
Diffusion coefficient D(10-4μm2/s)
The diffusion coefficient D for HP1β differed between control cells and those treated by HDAC inhibitors (TSA and vorinostat) and actinomycin D (Table 3). Conversely, changes in BMI1 diffusion were not so pronounced after selected treatments, whereas the diffusion coefficient for TRF1 significantly changed after all treatments (Table 3). Conclusions related to diffusion coefficient are closely related to trends expressed by the MSD (Figure 5).
Recovery of fluorescence after photobleaching of HP1β, BMI1, and TRF1 proteins after HDAC inhibition and suppression of transcription by actinomycin D
We confirmed that the fluorescence recovery kinetics of HP1β is more rapid in euchromatin than in heterochromatic chromocentres ( and Figure 6a). However, the recovery time after photobleaching was similar for HP1β in euchromatin away from foci (recovery after 10 s [R10] was 72%) and in nucleoli [R10] = 74% (Figure 6a). Interestingly, after TSA treatment, the curves were identical for HP1β in euchromatin ([R10] = 76%), in heterochromatic foci ([R10] = 75%), and in the nucleolus ([R10] = 79%) (Figure 6b). When transcription was suppressed by actinomycin D, the fluorescence recovery kinetics of HP1β was significantly slower only in chromocentres ([R10] = 47%) (foci) (Figure 6c). We did not observe statistically significant changes when we compared the kinetics of HP1β in euchromatin (Figure 6d) and within nucleoli (Figure 6e) after TSA and actinomycin D treatment. However, the recovery time after photobleaching for the HP1β in foci significantly decreased after actinomycin D and increased after TSA addition (Figure 6f). Thus, we confirmed the effect of TSA on heterochromatin regions, which has also been published by Cheutin et al.).
Fluorescence recovery of BMI1 was rapid when it was diffuse in the nucleoplasm ([R10] = 78%) (Figure 7a), and slower when it was accumulated in PcG bodies (foci) ([R10] = 36%) (Figure 7a). After TSA treatment, this trend was similar, but the differences between the BMI1 molecules outside ([R10] = 80%) and within PcG foci ([R10] = 28%) were more pronounced (Figure 7b). Moreover, actinomycin D significantly reduced the fluorescence recovery time after photobleaching of BMI1 protein accumulated in foci ([R10] = 20%) (Figure 7c). We measured a slower recovery for BMI1 away from PcG bodies when the cells were treated with actinomycin D (Figure 7d), but this was not statistically significant. Interestingly, the BMI1 in PcG bodies recovered more slowly after TSA treatment, but the difference was not as strong as after actinomycin D treatment (Figure 7e).
The telomeric TRF1 protein also showed significant differences depending on whether or not TRF1 was accumulated in foci. TRF1 foci are likely clusters consisting of several telomeres, consistent with the observation that cells expressing GFP-TRF1 had fewer TRF1-positive foci than telomeres. We measured [R10] = 81% for TRF1 outside of foci (likely the unbound form), and [R10] = 60% for TRF1 localised to telomere clusters (Figure 8a). Changes in TRF1 fluorescence recovery time were also observed in cells treated with TSA or actinomycin D (Figure 8b,c). After TSA treatment, [R10] = 72% for unbound TRF1 and [R10] = 52% for TRF1 within foci (Figure 8b). Similarly, when cells were treated with actinomycin D, [R10] = 64% for TRF1 outside of foci, but for TRF1 within foci, [R10] = 37% (Figure 8c). We also compared TRF1 fluorescence recovery after photobleaching in untreated control cells, and TSA-treated and actinomycin D-treated cells. After actinomycin D treatment, fluorescence recovery time after photobleaching was not significantly reduced for TRF1 outside of foci (Figure 8d). However, TRF1 recovery was significantly reduced within foci, particularly in cells treated with actinomycin D (Figure 8e).
Kinetic modelling of HP1β, RPA194, and UBF dynamics with modulated transcription levels
In the present work, we have uncovered a link between the dynamics and trajectories of chromatin-related proteins, and transcriptional regulation by TSA or an inhibitor of transcription elongation, actinomycin D. We have also addressed the effect of vorinostat (SAHA), another HDAC inhibitor that has promising therapeutic applications in oncology.
Cheutin et al. reported distinct kinetics for HP1 subtypes (HP1α, HP1β, HP1γ) depending on whether they are accumulated in clusters of heterochromatin or are in euchromatin. Meshorer et al. published that the recovery of heterochromatin-related HP1α is more rapid in mESCs than in terminally differentiated neurons derived from mESCs. Here, we report that in nucleoli, where HP1β appears in dispersed form , the kinetics of HP1β are similar to those described by Cheutin et al. at euchromatin (Figure 6d,e). According to these results, it is evident that protein density and protein affiliation with heterochromatin or euchromatin influence protein mobility. This is illustrated by the observations that proteins recovered more slowly when accumulated in foci (HP1β, BMI1, TRF1) than when evenly dispersed throughout the nucleoplasm (Figures 6, 7, 8, 9). Moreover, the values for fluorescence recovery correlated with protein molecular weight only for proteins that were homogeneously distributed throughout the nucleoplasm (Figure 2b). These data suggest that it is not molecular weight alone, but molecular weight together with nuclear pattern that influences protein kinetics. This is supported by the lack of correlation between molecular weight and recovery of 18 proteins that were analysed irrespective of nuclear pattern (Figure 2a). Our results are consistent with Jacobson and Wojcieszyn , who reported that diffusion coefficients do not depend on molecular weight. However, Braga et al. published that diffusion coefficients are lower in the nucleoplasm for higher molecular weight dextran, which is a temperature-dependent effect. Based on these results, it seems that biologically inactive molecules move differently than physiologically active proteins of similar size [38–40].
We also showed that, for all examples tested, suppression of transcription by actinomycin D had a strong impact on HP1β, BMI1, and TRF1 mobility. It is evident that actinomycin D, an agent that intercalates into double stranded DNA, robustly disturbed nuclear structure and induced changes in the mobility of proteins. Intriguingly, we observed more pronounced differences in the recovery time of HP1β and BMI1 in heterochromatic foci than in euchromatin regions after treatment with actinomycin D or HDAC inhibitor (Figures 6 and 7). For the most part, distinct global transcription levels were associated with decreased recovery kinetics of the studied proteins. Generally, protein residence time was more frequently reduced when core histones were hyperacetylated (; Figures 7b,e and 8e). Similarly to Cheutin et al., we observed that HP1β was an exception to this rule. In this case, histone hyperacetylation-induced by TSA significantly increased the mobility of HP1β in heterochromatic foci (Figure 6f). This supports previous data on the effect of HDAC inhibitors (TSA and sodium butyrate) on heterochromatin; HDAC inhibitors disturb the nuclear arrangement of HP1 subtypes and HP1 colocalisation within centromeric heterochromatin [41, 42]. In contrast, we show here that vorinostat had a more moderate effect on protein mobility than TSA. Thus, the analysis of protein mobility after cytostatic treatments might reveal important side effects of clinically used drugs.
Interestingly, when we increased the transcriptional activity of ribosomal genes along with increased histone acetylation by treating cells with TSA , we observed no effect on RPA194 recovery time, similar to what was observed for HP1β in nucleoli (Figure 10a,b). However, a significantly decreased recovery time was observed for UBF after TSA treatment (Figure 10c). Because RPA194, UBF, and HP1β proteins are all associated with active transcription of ribosomal genes [11, 35], these distinctions may be consequences of protein quantity at particular genomic regions and/or related to protein binding efficiency to other proteins or DNA in this region. Similarly, the local concentrations of proteins can influence the mobilisation and kinetics of particular proteins, as discussed in two papers by Ayoub et al.[43, 44], especially for HP1β recruited to UV-damaged chromatin. Another explanation is that the recovery kinetics after photobleaching may be influenced by cell cycle changes. For RNA polymerase I, the recovery kinetics decrease as transcription increases during the cell cycle . In S-phase, there is a reduced recovery after photobleaching of RPA194 compared with G1 cells . Similarly, the low recovery properties of histone H2B and H3 (, and Figure 2a,b) were significantly reduced when proliferating mESCs were induced to neuronal differentiation, characterised by a cell cycle block in G0 phase . As stated by Melcer and Meshorer , this phenomenon is not without exception. For example, the recovery kinetics of histone H3.3 is similarly hypodynamic in both pluripotent mESCs and after neuronal differentiation. Moreover, no differences in histone H1o kinetics were observed during different cell cycle stages .
Our experiments were directed towards the detailed detection of trajectories of selected protein foci. We also measured a correlation between TRF1 foci trajectory and transcriptional activity, but it was not upheld for HP1β and BMI1 (Figures 3, 4, 5). This is similar to data published for nuclear gene positioning with respect to gene expression. For example Meaburn and Misteli  identified several genes that are spatially repositioned during breast cancer tumorigenesis, but observed gene activity-independent genome repositioning in the early stages of tumour formation. Recently, we showed that the nuclear radial position of the pluripotency gene Oct3/4 did not change when downregulated . Thus, it seems that nuclear radial position is highly gene specific and likely depends on the transcriptional activity of surrounding chromatin and/or, as stated by Küpper et al., radial chromatin positioning is shaped by local gene density, not by gene expression.
Together, the trajectories of HP1β, BMI1, and TRF1 foci were, in many cases, influenced by histone hyperacetylation and by suppression of transcription by actinomycin D (Figures 3, 4, 5). However, the localised dynamics or extended dynamics of other nuclear domains, such as PML bodies, are not changed after inhibition of RNA polymerase II by α-amanitin . Here, we found that a protein trajectory was dependent on foci positioning; it matters whether the foci were located at the nuclear periphery or in the nuclear interior (Figures 3 and 4). Heterogeneity in the mobility of protein foci has also been reported by Guan et al. for α1B-adrenoceptor in living cells. Moreover, increased energy dependent motion was observed by Wang et al. for shorter telomeres when compared with longer and uncapped telomeres. Here, we show that both the trajectory and the average velocity were specific to nuclear positioning for HP1β, BMI1, and TRF1 foci (Table 1). This heterogeneity in foci movement is consistent with the observation of 'calm' and 'jittering' telomeres in live cells . Similarly, Muratani et al. have published velocity differences between individual PML bodies. Trajectories of PML bodies were defined as quasilinear, which is distinct from the trajectories of HP1β, BMI1, and TRF1 foci (Figures 3, 4, 5). These data demonstrate protein-specific kinetics and trajectories that likely reflect several biological events.
We have described protein dynamics that can be influenced by several factors, including binding and release constants, residence times, diffusion coefficients , chromatin condensation , cell differentiation , and cell cycle-dependent transcriptional level . In addition, Odenheimer et al., suggest there are strong effect of chromatin nanostructure on local mobility. In our experimental model, no single factor was directly responsible for slow or fast fluorescence recovery after photobleaching. Thus, we propose that protein dynamics are likely influenced by several factors and cellular processes that combine to determine the kinetic properties of chromatin-related proteins.
Immortalised wild-type (wt) MEFs originated in the Laboratory of Professor Thomas Jenuwein, Max-Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany. HeLa cells expressing histone H2B-GFP were a generous gift from Dr Marion Cremer (Ludwig-Maximilians University, Munich, Germany) and photoconvertible Dendra2 was used to label histone H4 in HepG2 cells (obtained from Professor Ivan Raška, Institute of Cellular Biology and Pathology, First Faculty of Medicine, Charles University, Prague, Czech Republic). U2OS cells expressing BMI1-GFP were a generous gift from Associate Professor Dušan Cmarko (Institute of Cellular Biology and Pathology, First Faculty of Medicine, Charles University, Prague, Czech Republic). Dr Paul Verbruggen (Swammerdam Institute for Life Sciences, University of Amsterdam, Amsterdam, The Netherlands) provided 3T3 cells with stable HP1β expression. The majority of cells were cultivated in high-glucose Dulbecco's modified Eagle medium (DMEM) containing 10% sera, but MEFs were cultivated according to Harničarová Horáková et al. and Dendra-HepG2 were grown in low-glucose DMEM supplemented with 10% sera. When the cultures reached 70% confluence, the cells were treated with a final concentration of 100 nM trichostatin A (TSA; Sigma-Aldrich, St Louis, MO, USA) and 0.5 μg/ml actinomycin D (#A9415, Sigma-Aldrich). The hyperacetylation effects of TSA have been verified elsewhere [41, 52] and actinomycin D is a commonly used reagent to block transcription elongation mediated by RNA polymerases [20, 21, 53].
The mESC line, GOWT1 (a generous gift from Hitoshi Niwa, Laboratory for Pluripotent Stem Cell Studies, RIKEN Center for Developmental Biology, Kobe, Japan), was cultivated in standard mESC medium (Glasgow minimum essential medium (GMEM) + 10% foetal calf serum (FCS)) with leukaemia inhibitory factor (LIF). These cells were maintained in the presence of puromycin (1.5 μg/ml) to select for Oct3/4 positive undifferentiated stem cells. All cell cultures were maintained at 37°C in a humidified atmosphere containing 5% CO2.
The plasmids were transformed into Escherichia coli DH5α for amplification. Plasmid DNA was isolated using the QIAGEN Large-Construct kit (#12462; QIAGEN, Bio-Consult, Prague, Czech Republic). For transfection of MEFs we used the METAFECTENE PRO system (Biontex Laboratories GmbH, Planegg, Germany). The following plasmids were used in this study: JMJD2b-GFP (obtained from Professor Thomas Jenuwein and Dr Nicholas Shukeir, Max-Planck Institute of Immunobiology, Freiburg, Germany); HP1α-GFP and HP1β-GFP (from the laboratory of Dr Tom Misteli, National Institutes of Health, Bethesda, MD, USA); RPA194-GFP (Addgene Inc., Cambridge, MA, USA; #17660); UBF-GFP (Addgene; #17656); Ubiquitin-GFP (Addgene; #11928); p53-GFP (Addgene; #12091); pmCherry-alpha-tubulin-IRES-puro2 (Addgene; #1360); lamin A-GFP (Addgene; #17662); c-myc-GFP (a gift from Dr Hiroyoshi Ariga, Hokkaido University, Graduate School of Pharmaceutical Sciences, Kita-Ku, Sapporo, Japan), β-catenin-enhanced GFP and STAT1-GFP (gifts from Dr Vítězslav Bryja and Dr Jiří Pacherník, Faculty of Sciences, Masaryk University, Brno, Czech Republic). GFP-hTRF1-pS65T-C1 was a gift from Dr C M Counter, Duke University Medical Center, Durham, NC, USA. GFP-PML was obtained from Douglas Durso, Customer Service Representative, Roche NimbleGen, Inc., Madison, WI, USA. pEGFP-C1 cloning plasmid was purchased from Clontech Laboratories, Inc. (Mountain View, CA, USA).
Single particle tracking analysis
Time series of images were acquired on a confocal Leica TSC SP 5× microscope (Leica, Mannheim, Germany). We captured images every 2.6 s for 20 min. The image analysis was performed using Matlab software http://www.mathworks.com/products/matlab/. Each frame was cropped to contain a single cell. Weighted centres of foci were computed for every frame (see below). The centres formed point sets, which were used to compensate for cellular movement, as described in . It was impossible to localise all foci in all frames (mostly due to out-of-focus motion). Therefore, the trajectories were made in two steps: (1) all corresponding points in successive frames were linked together. The correspondence was identified by a point-based alignment algorithm. This step produced many foci subtracks. In this step, we determined the area of the smallest ellipse required to enclose all points in each subtrack. (2) We combined subtracks in one track if their minimal ellipses intersected and followed each other in time. The collection of subtracks defined the trajectory of a focus. We calculated the average velocity of the foci as the length of all trajectory subtracks divided by the sum of the subtracks' time spans. We measured the area of the minimal enclosing ellipse of all trajectory points. We studied the foci trajectories in the interior of nucleus, as well as at the nuclear periphery. The trajectories that had a local radius less of than 60% for the initial frame were taken to be the central trajectory. Local radius has been defined elsewhere .
We calculated the weighted centres of the foci using the following procedure. First, to remove noise, each frame of a sequence was filtered through a median filter with a circular window of 7 pixels. The image was then smoothed by convolution with Gaussian kernel (σ = 3). Next, we computed the morphological h-dome transform with manually selected h value. Only domes containing a pixel intensity exactly equal to h were considered objects of interest. This was achieved using morphological reconstructions of h-dome images from seeds corresponding to pixels of intensity h. Finally, we computed the weighted centres within h-domes of detected objects (foci).
To compute the local radius of points, we had to segment the nucleus, which was performed using the following algorithm. First, the image was smoothed by convolution with a Gaussian kernel σ = 10), and then the threshold was adjusted using the chord method . Morphological closing of a suitable radius was performed to fill the holes.
We calculated MSD and diffusion coefficient as previously described . We used a modified MSD formula for sequences with missing probes as previously described . The MSD curve was calculated for all trajectories up to the time lag of 100 frames (260 s). We present average curves for different cells in Figure 5.
The diffusion coefficient was estimated by linear fitting to the MSD curve. The optimal number of fitting points was determined using the iterative scheme proposed in . We calculated the diffusion coefficient from the average MSD curves.
Live cell imaging, FRAP analysis, and statistics
Cells growing on 50-mm glass-bottomed dishes (MatTek, Corporation, Ashland, MA, USA; #P50G-0-30-F) were placed in a cultivation hood (EMBL, Heidelberg, Germany) with a stable temperature of 37°C and 5% CO2. For live cell studies, we used a confocal Leica TSC SP-5× microscope, equipped with a white light laser (470 to 670 nm); argon laser (488 nm), and UV lasers (355 nm and 405 nm); 64× magnification and numerical aperture NA = 1.4 were used. The movement of individual cells and chromatin-related proteins, either accumulated in distinct foci or dispersed throughout the nucleoplasm, was monitored using Leica LAS AF software (version 2.1.2). For trajectory analysis, the cells were monitored for 20 min; scanning was performed every 2 s. For FRAP, GFP was excited at 488-nm using an argon laser. During FRAP, fluorescent molecules were bleached at defined regions of interest (ROI) (2 μm2). Bleaching was performed with a 488-nm argon laser for 1.2 s. For scanning, we used 10% laser intensity and for FRAP we used 100% laser power. To minimise the laser instability that can be caused by external temperature, the room temperature was stabilised using a highly precise air conditioner (room temperature was 21°C to 22°C). Moreover, a particular set of experiments was analysed by FRAP in a single day to minimise the effects of laser instability. Prebleaching was set to five frames (each frame is 0.66 s), and the rate of fluorescence recovery provided information about the rate at which the fluorescent molecules moved into the bleached region. The cells were monitored over 0.66 s intervals, for 25 s.
To study correlations, we used linear regression analysis as a statistical tool in the SigmaPlot software (version 8.0; SPSS, Chicago, IL, USA). Potential correlation between molecular weight (MW) and fluorescence recovery after photobleaching was evaluated by the Pearson's product-moment correlation coefficient (r). We constructed linear regression curves and a confidence interval around (r) that showed 95% probability of statistically significant results.
This work was supported by the Ministry of Education Youth and Sports of the Czech Republic; the research projects LC535, LC06027, ME919, and by the Academy of Sciences of the Czech Republic; projects AVOZ50040702 and AVOZ50040507. This work was also supported by EU project COST TD09/05 and by Grant Agency of Czech Republic, project P302/10/1022. Many thanks the scientists (mentioned in Methods section) who provided cells and/or plasmids and Dr Jiří Pacherník for GOWT1 mESC cultivation. We are also grateful to BioScience Writers (Houston, TX, USA) for the critical linguistic revision of our manuscript.
- Belmont AS: Visualizing chromosome dynamics with GFP. Trends Cell Biol. 2001, 11: 250-257. 10.1016/S0962-8924(01)02000-1.View ArticlePubMedGoogle Scholar
- Meshorer E: Imaging chromatin in embryonic stem cells. StemBook. 2008, The Stem Cell Research Community, [http://www.ncbi.nlm.nih.gov/books/NBK27072/]Google Scholar
- Misteli T: Beyond the sequence: cellular organization of genome function. Cell. 2007, 128: 787-800. 10.1016/j.cell.2007.01.028.View ArticlePubMedGoogle Scholar
- Lanctot C, Cheutin T, Cremer M, Cavalli G, Cremer T: Dynamic genome architecture in the nuclear space: regulation of gene expression in three dimensions. Nat Rev Genet. 2007, 8: 104-115. 10.1038/nrg2041.View ArticlePubMedGoogle Scholar
- Mattout A, Meshorer E: Chromatin plasticity and genome organization in pluripotent embryonic stem cells. Curr Opin Cell Biol. 2010, 22: 334-341. 10.1016/j.ceb.2010.02.001.View ArticlePubMedGoogle Scholar
- Branco MR, Pombo A: Intermingling of chromosome territories in interphase suggests role in translocations and transcription-dependent associations. PLoS Biol. 2006, 4: e138-10.1371/journal.pbio.0040138.PubMed CentralView ArticlePubMedGoogle Scholar
- Cremer T, Cremer C: Chromosome territories, nuclear architecture and gene regulation in mammalian cells. Nat Rev Genet. 2001, 2: 292-301. 10.1038/35066075.View ArticlePubMedGoogle Scholar
- Probst AV, Almouzni G: Pericentric heterochromatin: dynamic organization during early development in mammals. Differentiation. 2008, 76: 15-23.View ArticlePubMedGoogle Scholar
- Blackburn EH: Telomeres and telomerase: their mechanisms of action and the effects of altering their functions. FEBS Lett. 2005, 579: 859-862. 10.1016/j.febslet.2004.11.036.View ArticlePubMedGoogle Scholar
- Blasco MA: Telomere length, stem cells and aging. Nat Chem Biol. 2007, 3: 640-649. 10.1038/nchembio.2007.38.View ArticlePubMedGoogle Scholar
- Prieto JL, McStay B: Pseudo-NORs: a novel model for studying nucleoli. Biochim Biophys Acta. 2008, 1783: 2116-2123. 10.1016/j.bbamcr.2008.07.004.View ArticlePubMedGoogle Scholar
- Raska I, Shaw PJ, Cmarko D: New insights into nucleolar architecture and activity. Int Rev Cytol. 2006, 255: 177-235. 10.1016/S0074-7696(06)55004-1.View ArticlePubMedGoogle Scholar
- Raska I, Shaw PJ, Cmarko D: Structure and function of the nucleolus in the spotlight. Curr Opin Cell Biol. 2006, 18: 325-334. 10.1016/j.ceb.2006.04.008.View ArticlePubMedGoogle Scholar
- Alcobia I, Quina AS, Neves H, Clode N, Parreira L: The spatial organization of centromeric heterochromatin during normal human lymphopoiesis: evidence for ontogenically determined spatial patterns. Exp Cell Res. 2003, 290: 358-369. 10.1016/S0014-4827(03)00335-5.View ArticlePubMedGoogle Scholar
- Boisvert FM, van Koningsbruggen S, Navascues J, Lamond AI: The multifunctional nucleolus. Nat Rev Mol Cell Biol. 2007, 8: 574-585. 10.1038/nrm2184.View ArticlePubMedGoogle Scholar
- McStay B, Grummt I: The epigenetics of rRNA genes: from molecular to chromosome biology. Annu Rev Cell Dev Biol. 2008, 24: 131-157. 10.1146/annurev.cellbio.24.110707.175259.View ArticlePubMedGoogle Scholar
- Grummt I: Life on a planet of its own: regulation of RNA polymerase I transcription in the nucleolus. Genes Dev. 2003, 17: 1691-1702. 10.1101/gad.1098503R.View ArticlePubMedGoogle Scholar
- Russell J, Zomerdijk JC: RNA-polymerase-I-directed rDNA transcription, life and works. Trends Biochem Sci. 2005, 30: 87-96. 10.1016/j.tibs.2004.12.008.View ArticlePubMedGoogle Scholar
- Gorski SA, Snyder SK, John S, Grummt I, Misteli T: Modulation of RNA polymerase assembly dynamics in transcriptional regulation. Mol Cell. 2008, 30: 486-497. 10.1016/j.molcel.2008.04.021.PubMed CentralView ArticlePubMedGoogle Scholar
- Perry RP, Kelley DE: Inhibition of RNA synthesis by actinomycin D: characteristic dose-response of different RNA species. J Cell Physiol. 1970, 76: 127-139. 10.1002/jcp.1040760202.View ArticlePubMedGoogle Scholar
- Sobell HM: Actinomycin and DNA transcription. Proc Natl Acad Sci USA. 1985, 82: 5328-5331. 10.1073/pnas.82.16.5328.PubMed CentralView ArticlePubMedGoogle Scholar
- Weiland Y, Lemmer P, Cremer C: Combining FISH with localisation microscopy: super-resolution imaging of nuclear genome nanostructures. Chromosome Res. 2011, 19: 5-23. 10.1007/s10577-010-9171-6.View ArticlePubMedGoogle Scholar
- Egner A, Verrier S, Goroshkov A, Soling HD, Hell SW: 4Pi-microscopy of the Golgi apparatus in live mammalian cells. J Struct Biol. 2004, 147: 70-76. 10.1016/j.jsb.2003.10.006.View ArticlePubMedGoogle Scholar
- Schermelleh L, Carlton PM, Haase S, Shao L, Winoto L, Kner P, Burke B, Cardoso MC, Agard DA, Gustafsson MG, Leonhardt H, Sedat JW: Subdiffraction multicolor imaging of the nuclear periphery with 3D structured illumination microscopy. Science. 2008, 320: 1332-1336. 10.1126/science.1156947.PubMed CentralView ArticlePubMedGoogle Scholar
- Hell SW, Wichmann J: Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy. Opt Lett. 1994, 19: 780-782. 10.1364/OL.19.000780.View ArticlePubMedGoogle Scholar
- Misteli T: Protein dynamics: implications for nuclear architecture and gene expression. Science. 2001, 291: 843-847. 10.1126/science.291.5505.843.View ArticlePubMedGoogle Scholar
- Hernandez-Munoz I, Lund AH, van der Stoop P, Boutsma E, Muijrers I, Verhoeven E, Nusinow DA, Panning B, Marahrens Y, van Lohuizen M: Stable X chromosome inactivation involves the PRC1 polycomb complex and requires histone MACROH2A1 and the CULLIN3/SPOP ubiquitin E3 ligase. Proc Natl Acad Sci USA. 2005, 102: 7635-7640. 10.1073/pnas.0408918102.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Kam Z, Carlton PM, Xu L, Sedat JW, Blackburn EH: Rapid telomere motions in live human cells analyzed by highly time-resolved microscopy. Epigenetics Chromatin. 2008, 1: 4-10.1186/1756-8935-1-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Festenstein R, Pagakis SN, Hiragami K, Lyon D, Verreault A, Sekkali B, Kioussis D: Modulation of heterochromatin protein 1 dynamics in primary Mammalian cells. Science. 2003, 299: 719-721. 10.1126/science.1078694.View ArticlePubMedGoogle Scholar
- Cheutin T, McNairn AJ, Jenuwein T, Gilbert DM, Singh PB, Misteli T: Maintenance of stable heterochromatin domains by dynamic HP1 binding. Science. 2003, 299: 721-725. 10.1126/science.1078572.View ArticlePubMedGoogle Scholar
- Meshorer E, Yellajoshula D, George E, Scambler PJ, Brown DT, Misteli T: Hyperdynamic plasticity of chromatin proteins in pluripotent embryonic stem cells. Dev Cell. 2006, 10: 105-116. 10.1016/j.devcel.2005.10.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Minc E, Allory Y, Worman HJ, Courvalin JC, Buendia B: Localization and phosphorylation of HP1 proteins during the cell cycle in mammalian cells. Chromosoma. 1999, 108: 220-234. 10.1007/s004120050372.View ArticlePubMedGoogle Scholar
- Sparmann A, van Lohuizen M: Polycomb silencers control cell fate, development and cancer. Nat Rev Cancer. 2006, 6: 846-856. 10.1038/nrc1991.View ArticlePubMedGoogle Scholar
- Sykorova E, Fajkus J: Structure-function relationships in telomerase genes. Biol Cell. 2009, 101: 375-392. 10.1042/BC20080205.View ArticlePubMedGoogle Scholar
- Horakova AH, Bartova E, Galiova G, Uhlirova R, Matula P, Kozubek S: SUV39h-independent association of HP1 beta with fibrillarin-positive nucleolar regions. Chromosoma. 2010, 119: 227-241. 10.1007/s00412-009-0252-2.View ArticlePubMedGoogle Scholar
- Jacobson K, Wojcieszyn J: The translational mobility of substances within the cytoplasmic matrix. Proc Natl Acad Sci USA. 1984, 81: 6747-6751. 10.1073/pnas.81.21.6747.PubMed CentralView ArticlePubMedGoogle Scholar
- Braga J, Desterro JM, Carmo-Fonseca M: Intracellular macromolecular mobility measured by fluorescence recovery after photobleaching with confocal laser scanning microscopes. Mol Biol Cell. 2004, 15: 4749-4760. 10.1091/mbc.E04-06-0496.PubMed CentralView ArticlePubMedGoogle Scholar
- Houtsmuller AB, Rademakers S, Nigg AL, Hoogstraten D, Hoeijmakers JH, Vermeulen W: Action of DNA repair endonuclease ERCC1/XPF in living cells. Science. 1999, 284: 958-961. 10.1126/science.284.5416.958.View ArticlePubMedGoogle Scholar
- Kruhlak MJ, Lever MA, Fischle W, Verdin E, Bazett-Jones DP, Hendzel MJ: Reduced mobility of the alternate splicing factor (ASF) through the nucleoplasm and steady state speckle compartments. J Cell Biol. 2000, 150: 41-51. 10.1083/jcb.150.1.41.PubMed CentralView ArticlePubMedGoogle Scholar
- Phair RD, Misteli T: High mobility of proteins in the mammalian cell nucleus. Nature. 2000, 404: 604-609. 10.1038/35007077.View ArticlePubMedGoogle Scholar
- Bartova E, Pachernik J, Harnicarova A, Kovarik A, Kovarikova M, Hofmanova J, Skalnikova M, Kozubek M, Kozubek S: Nuclear levels and patterns of histone H3 modification and HP1 proteins after inhibition of histone deacetylases. J Cell Sci. 2005, 118: 5035-5046. 10.1242/jcs.02621.View ArticlePubMedGoogle Scholar
- Taddei A, Maison C, Roche D, Almouzni G: Reversible disruption of pericentric heterochromatin and centromere function by inhibiting deacetylases. Nat Cell Biol. 2001, 3: 114-120. 10.1038/35055010.View ArticlePubMedGoogle Scholar
- Ayoub N, Jeyasekharan AD, Bernal JA, Venkitaraman AR: HP1-beta mobilization promotes chromatin changes that initiate the DNA damage response. Nature. 2008, 453: 682-686. 10.1038/nature06875.View ArticlePubMedGoogle Scholar
- Ayoub N, Jeyasekharan AD, Venkitaraman AR: Mobilization and recruitment of HP1: a bimodal response to DNA breakage. Cell Cycle. 2009, 8: 2945-2950.PubMedGoogle Scholar
- Melcer S, Meshorer E: Chromatin plasticity in pluripotent cells. Essays Biochem. 2010, 48: 245-262. 10.1042/bse0480245.View ArticlePubMedGoogle Scholar
- Meaburn KJ, Misteli T: Locus-specific and activity-independent gene repositioning during early tumorigenesis. J Cell Biol. 2008, 180: 39-50. 10.1083/jcb.200708204.PubMed CentralView ArticlePubMedGoogle Scholar
- Bartova E, Krejci J, Harnicarova A, Kozubek S: Differentiation of human embryonic stem cells induces condensation of chromosome territories and formation of heterochromatin protein 1 foci. Differentiation. 2008, 76: 24-32.View ArticlePubMedGoogle Scholar
- Kupper K, Kolbl A, Biener D, Dittrich S, von Hase J, Thormeyer T, Fiegler H, Carter NP, Speicher MR, Cremer T, Cremer M: Radial chromatin positioning is shaped by local gene density, not by gene expression. Chromosoma. 2007, 116: 285-306. 10.1007/s00412-007-0098-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Muratani M, Gerlich D, Janicki SM, Gebhard M, Eils R, Spector DL: Metabolic-energy-dependent movement of PML bodies within the mammalian cell nucleus. Nat Cell Biol. 2002, 4: 106-110. 10.1038/ncb740.View ArticlePubMedGoogle Scholar
- Guan Y, Xu M, Liang Z, Xu N, Lu Z, Han Q, Zhang Y, Zhao XS: Heterogeneous transportation of alpha1B-adrenoceptor in living cells. Biophys Chem. 2007, 127: 149-154. 10.1016/j.bpc.2007.01.009.View ArticlePubMedGoogle Scholar
- Odenheimer J, Heermann DW, Kreth G: Brownian dynamics simulations reveal regulatory properties of higher-order chromatin structures. Eur Biophys J. 2009, 38: 749-756. 10.1007/s00249-009-0486-1.View ArticlePubMedGoogle Scholar
- Bartova E, Pachernik J, Kozubik A, Kozubek S: Differentiation-specific association of HP1alpha and HP1beta with chromocentres is correlated with clustering of TIF1beta at these sites. Histochem Cell Biol. 2007, 127: 375-388. 10.1007/s00418-006-0259-1.View ArticlePubMedGoogle Scholar
- Nguyen VT, Giannoni F, Dubois MF, Seo SJ, Vigneron M, Kedinger C, Bensaude O: In vivo degradation of RNA polymerase II largest subunit triggered by alpha-amanitin. Nucleic Acids Res. 1996, 24: 2924-2929. 10.1093/nar/24.15.2924.PubMed CentralView ArticlePubMedGoogle Scholar
- Matula P, Kozubek M, Dvorak V: Fast point-based 3-D alignment of live cells. IEEE Trans Image Process. 2006, 15: 2388-2396. 10.1109/TIP.2006.875209.View ArticlePubMedGoogle Scholar
- Uhlírová R, Horáková AH, Galiová G, Legartová S, Matula P, Fojtová M, Varecha M, Amrichová J, Vondrácek J, Kozubek S, Bártová E: SUV39h- and A-type lamin-dependent telomere nuclear rearrangement. J Cell Biochem. 2010, 109: 915-926.PubMedGoogle Scholar
- Zack GW, Rogers WE, Latt SA: Automatic measurement of sister chromatid exchange frequency. J Histochem Cytochem. 1977, 25: 741-753.View ArticlePubMedGoogle Scholar
- Michalet X: Mean square displacement analysis of single-particle trajectories with localization error: Brownian motion in an isotropic medium. Phys Rev E Stat Nonlin Soft Matter Phys. 2011, 82: 041914-10.1103/PhysRevE.82.041914.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.