|Home | About | Journals | Submit | Contact Us | Français|
We present a 3-D mapping in WI38 human diploid fibroblast cells of chromosome territories (CT) 13,14,15,21, and 22, which contain the nucleolar organizing regions (NOR) and participate in the formation of nucleoli. The nuclear radial positioning of NOR-CT correlated with the size of chromosomes with smaller CT more interior. A high frequency of pairwise associations between NOR-CT ranging from 52% (CT13-21) to 82% (CT15-21) was detected as well as a triplet arrangement of CT15-21-22 (72%). The associations of homologous CT were significantly lower (24–36%). The arrangements of each pairwise CT varied from CT13-14 and CT13-22, which had a majority of cells with single associations, to CT13-15 and CT13-21 where a majority of cells had multiple interactions. In cells with multiple nucleoli, one of the nucleoli (termed “dominant”) always associated with a higher number of CT. Moreover, certain CT pairs more frequently contributed to the same nucleolus than to others. This nonrandom pattern suggests that a large number of the NOR-chromsomes are poised in close proximity during the postmitotic nucleolar recovery and through their NORs may contribute to the formation of the same nucleolus. A global data mining program termed the chromatic median determined the most probable interchromosomal arrangement of the entire NOR-CT population. This interactive network model was significantly above randomized simulation and was composed of 13 connections among the NOR-CT. We conclude that the NOR-CT form a global interactive network in the cell nucleus that may be a fundamental feature for the regulation of nucleolar and other genomic functions.
The human genome is composed of 23 pairs of chromosomes. During interphase, chromosomes are organized in the nucleus as discrete chromosome territories (CT) with well distinguishable boundaries (Bickmore, 2013; Cremer and Cremer, 2010; Cremer et al., 2006; Cremer et al., 2000; Meaburn and Misteli, 2007b; Visser and Aten, 1999). The three dimensional position of chromosomes in the cell nucleus is considered to be important for determining the overall nuclear architecture of eukaryotes, the maintenance of genome stability and the epigenetic regulation of gene expression (Berezney, 2002; Berezney et al., 2005; Bickmore, 2013; Cremer and Cremer, 2010; Kumaran et al., 2008; Malyavantham et al., 2010; Meaburn and Misteli, 2007a; Misteli, 2007; Stein et al., 2003; Stein G, 2008; T, 2005; Zaidi et al., 2007). A generally accepted notion postulates that the three dimensional positions of chromosomes in the cell nucleus are probabilistically non-random. For instance, non-randomness in CT positioning is correlated with gene density, where gene-rich chromosomes are found more interior than gene-poor chromosomes within the nucleus (Boyle et al., 2001; Heride et al., 2010; Kreth et al., 2004). The radial positions of chromosomes also correlate with their sequence length, where longer chromosomes tend to associate with the nuclear periphery, while shorter chromosomes are usually localized more internally relative to the cell nucleus (Bolzer et al., 2005; Marella et al., 2009a; Sun et al., 2000; Zeitz et al., 2009).
In addition to nonrandom radial positioning of CT within the nucleus, specific chromosome neighborhood arrangements have been determined (Bolzer et al., 2005; Nagele et al., 1999) which are different between cell and tissue types (Marella et al., 2009a; Mayer et al., 2005; Neusser et al., 2007; Parada et al., 2004a; Zeitz et al., 2009), and are altered during cell differentiation, development (Kuroda et al., 2004; Marella et al., 2009b), and in cancer cells (Brianna Caddle et al., 2007; Fritz et al., 2014; Parada et al., 2002). Interestingly, increased translocation frequencies have been determined between specific CT which are in closer proximity (Brianna Caddle et al., 2007; Folle, 2008; Parada et al., 2002; Roix et al., 2003; Soutoglou et al., 2007) and those that are actually overlapping or intermingling (Branco and Pombo, 2006). At the same time, the intra-nuclear locations of chromosomes are probabilistic rather than absolute, and, therefore, an understanding of CT distribution patterns is difficult to attain.
In this investigation we determined the interchromosomal and nucleolar associations among the acrocentric CT 13, 14, 15, 21 and 22. During nucleolar biogenesis, these chromosomes congregate in late telophase/early G1 to form large sub-nuclear compartments (Prieto and McStay, 2005; Raska et al., 2006). This clustering of acrocentric chromosomes is thought to occur because they bear the nucleolar organizing regions (NOR), which contain tandem arrays of ribosomal RNA (rRNA) genes. The NOR-chromosomes together harbor approximately 400 ribosomal genes per human diploid genome (Zentner et al., 2011). Synthesis and processing of ribosomal RNA and formation of ribosomes is an exceptionally intense biosynthetic activity which is morphologically expressed in the formation of the nucleolus, the largest structure-function compartment of the cell nucleus (Busch and Smetana, 1970; Koberna et al., 2002; Raska et al., 2006).
Not all the NOR-CT are active at any given time and they can be grouped into transcriptional competent and non-competent NOR (Grob et al., 2014; Kalmarova et al., 2008a; Kalmarova et al., 2007; Smirnov et al., 2006). Moreover, particular NOR-CT have greater numbers of transcriptionally competent NOR and concomitantly a greater number of those CT are adjacent to nucleoli (Kalmarova et al., 2007). Kamarova et. al. (2007, 2008) further found nonrandom patterns of NOR-CT associated with nucleoli in HeLa cells and in daughter cells following mitosis (Kalmarova et al., 2008b; Kalmarova et al., 2007).
Despite this progress, there is still only limited information on how the diverse pairs of NOR-CT cooperate in the formation of nucleoli, as well as the relationships between the distribution of NOR-CT and the number, size and intranuclear positions of the nucleolar domains. Moreover, it is unclear as to whether there are specific positional patterns of arrangement at both the pairwise level of interactions and globally among the entire set of NOR-CT. To a large extent, comprehensive studies on the intranuclear distribution patterns of the NOR-CT and nucleolar associations are scarce because it is technically challenging to perform simultaneous localization of all five pairs of these CT.
In this investigation of the NOR-CT, we took a different approach. The spatial organization and associations of chromosomes 13, 14, 15, 21 and 22 were performed together in the same cell nuclei using repetitive cycles of fluorescence in situ hybridization (re-FISH). We report highly non-random pair-wise association patterns of the NOR-CT and their proximity to nucleoli. Application of a novel computational geometric and data mining approach enables us to demonstrate for the first time a highly preferred, albeit probabilistic, interconnected network of the NOR-CT in the cell nucleus.
Labeling of chromosomes was performed in the WI38 diploid human lung fibroblast cell line (ATTC, Manassas, VA, USA). Cells were grown in Advanced DMEM (Life Technologies, Carlsbad, CA, USA) supplemented with 2.5% fetal bovine serum (FBS) and 1% antibiotic-antimycotic solution (Sigma-Aldrich, St. Louis, MO, USA) in 5% CO2 at 37°C. Prior to chromosome labeling, cells were placed on gridded CeLLocate coverslips (Eppendorf, Hamburg, Germany) and synchronized in G0 by serum starvation for 48 h in medium containing 0.1% FBS. This protocol yielded a nearly 100% synchronization in G0 as judged by the absence of Ki67 and BrdU staining or mitotic cells in the treated culture. The viability of treated cells was monitored by a live-dead cytotoxicity assay based on combined staining with calcein AM and ethidium homodimer (LifeTechnologies, Carlsbad, CA, USA). In this test, the cells with low activity of intracellular esterases and/or compromised integrity of membranes are recognized by a negative calcein AM signal (in the green spectral channel) and a positive ethidium staining (in the red spectral channel). Such cells were negligibly rare in the synchronized cell cultures.
Immunolabeling was performed with mouse anti- digoxigenin, anti-BrdU, anti Ki67 and rabbit anti-nucleolin antibodies (Abcam, Cambridge, England). Chromosome staining was performed by chromosome-specific aquarius paint probes (Cytocell, Windsor, CT), The hybridized segment of rDNA was prepared from a pA plasmid construct (Erickson et al., 1981) which codes for downstream 200 nucleotides of 18 S rRNA, internal transcribed spacer 1, sequence for 5.8 S rRNA, internal transcribed spacer 2 and upstream ~4,500 nucleotides of 28 S rRNA. The probe was labeled with digoxygenin-11-dUTP by nick translation.
An approach for 3-D mapping of CT in the cell nucleus termed “re-FISH” (Walter et al., 2006) was previously adapted by our group (Marella et al., 2009a; Zeitz et al., 2009). This method involves a sequential series of FISH labeling and imaging as follows: (i) two selected pairs of CT are labelled by commercial fluorescence probes, and mapped in a group of cells by high resolution 3D imaging. (ii) The fluorescence probes are striped and two new probes recognizing another two pairs of CT are applied. (iii) Cells studied in previous imaging session are identified and the positions of the new CT are mapped in these cells. Steps (ii) and (iii) can be repeated to map locations of additional CT in the cell nucleus. (iv) Stacks of optical sequences acquired in the same cells at subsequent imaging sessions are aligned in x-, y- and z- planes, CT are segmented and their mutual positions are analyzed. We have validated this technique for confidently mapping up to eight or nine pairs of CT in the cell nucleus (Fritz et al., 2014; Marella et al., 2009a; Zeitz et al., 2009). In order to combine FISH and immunofluorescence detection, we introduced minor alterations to the above approach. The WI38 cells synchronized in G0 were first fixed with 4% paraformaldehyde in PBS for 12 min and permeabilized in 0.5 % Triton in PBS for 5 min. The nucleoli were labeled with anti-nucleolin rabbit antibody (Abcam, Cambridge, UK) followed by Cy-5 anti-rabbit conjugate (Jackson Immuno Research, West Grove, PA, USA). After immunolabeling, the cells were post-fixed overnight with methanol:acetic acid (3:1) at −20°C to crosslink antibodies to the cellular structure. In experiments involving rDNA hybridization, cells were then incubated with 100 μg/ml RNAse A (Roche, Basel, Switzerland) for 2 h at 37°C, to eliminate rRNA and avoid cross hybridization with the rDNA probe. Then cells were incubated in 20% glycerol/PBS (30min), freeze-thawed in liquid nitrogen three times, treated with 0.1N HCl for 5 min and subjected to DNA denaturation (70% deionized formamide/2X SSC, pH 7.0) at 75°C for 3 min. Aquarius paint probes for two chromosome pairs (Cytocell, Windsor, CT) or rDNA probes were then denatured for 8 min at 75°C. Labeling of rDNA was usually performed in the same series of experiments as immunolabeling of nucleolin. Probe solution was applied to coverslips and sealed with rubber cement. Hybridization was carried out at 37°C for 48h. Three consecutive post-hybridization washes each for 30 min at 37°C were: 50% formamide/2X SSC/0.05% Tween-20; 2X SSC/0.05% Tween-20 and 1X SSC. Coverslips were then mounted on slides in Vectashield. Following image collection, chromosome paints were stripped by immersion of the coverslips in 50% formamide/2X SSC for 1 min at 65°C. Another pair of denatured chromosome paint probes were then immediately added to cells and hybridized at 37°C for 48h. This process was repeated one more time to obtain images of all five NOR chromosome pairs. Before mounting of the coverslips on microscopy slides, the genomic DNA was counterstained with 10 μM DAPI (Sigma-Aldrich, St. Louis, MO, USA).
The images were collected on an Olympus BX51 fluorescence microscope equipped with a Sensicam QE (Cooke Corporation, Romulus, MI) digital CCD camera, a motorized z-axis controller (Prior, Rockland, MA) and Slidebook 4.0 software (Intelligent Imaging Innovations, Denver, CO). Stacks of optical sections at 0.5 μm intervals were collected to map nucleoli, CT and rDNA signals. 3-D volume rendering of the collected optical sections for visualization was performed using Slidebook 4.0. The exact position of each cell in coordinates of the CeLLocate grid was documented by 3D phase contrast imaging. These phase-contrast images were used to localize the same cells in subsequent imaging cycles. Finally, z-stacks of the optical sections were aligned in 3D using registration software developed in our laboratory by selecting reference points from corresponding optical sections of phase contrast images as previously described (Zeitz et al., 2009). Following the alignment, optical stacks were merged and subjected to segmentation and quantitative analysis by an in-house developed computer program running on Matlab termed eFISHent (Fritz et al., 2014). This program determines a wide range of spatial distances in 3-D between the segmented objects (Fritz et al., 2014; Zeitz et al., 2009). One feature of this suite of 3-D distance measurements enabled determination of the nearest neighbor 3-D distances (edge to edge distances) among the labeled CT in each image set as well as among the nucleoli and CT. A positive pair-wise association between 2 CT was then scored if the nearest distance was ≤ 8 pixels which corresponds to ≤ 0.56 μm. Similarly, a positive association between a CT and a nucleolus was scored when the nearest edge to edge 3-D distance was ≤ 1.0 μm. For CT-nucleolar interactions, the CT were assigned an association with the nearest nucleolus. In the few rare cases where an NOR-CT was spaced an equal distance between two nucleoli, the neighborhood association was assigned to the nucleolus having the largest surface contact with that particular CT. Thirty nine image sets which contained all 5 pairs of NOR-bearing chromosome territories were analyzed in this study.
Computer simulations were run to determine to what extend the NOR-CT pair-wise interchromosomal associations and associations with the same nucleoli can be accounted for by random interactions of the NOR-CT. A constrained random simulation was performed as previous described for CT that are positioned in close association with the nuclear periphery (Zeitz et al, 2009). In this case, however, the NOR-CT are constrained to be placed at a distance of ≤ 1.0 μm from one of the nucleoli. To best mimic the experimental conditions, each image set was individually simulated. In the first step, the NOR-CT are removed from the each image set while the nuclear boundary and nucleoli remain fixed in position. Then the 5 pairs of NOR-CT are randomly placed back into the nucleus with the constraint that each NOR-CT is placed in close association with one or more of the nucleoli (see Fig 5). The pairwise association patterns and their pairwise associations with individual nucleoli are then calculated from the population of individual simulations. All computer simulations were run using Matlab.
Previously our group developed a computational geometric approach termed the Generalized Median Graph (GMG) to determine the best fit probabilistic model of chromosome territory interactions in G0 WI38 fibroblasts (Mukherjee et al., 2009; Zeitz et al., 2009). The GMG considered all CT associations (i.e., all permutations of the association graphs) and simultaneously identifies the association patterns of all CT. Recently we reported an improved data mining pattern recognition algorithm termed the chromatic median (CM) which uses combinatorial optimization to infer the common chromosome interaction pattern for the cell population (Ding et al., 2013). The CM is capable of determining the corresponding homolog of the same chromosome across all cells based upon the interactions that homolog has with other CT. In the mathematical formulation of the problem, we use a cell representation which is suitable for describing interactions of CT. The interactions within each nucleus are represented as a matrix where a zero indicates no interaction and a one indicates an interaction (Fig 1). The objective is to find the best permutation (relabeling from a to b and vice versa for all chromosome pairs within each cell) which aligns the association matrix of each input cell with that of the common pattern.
This new CM algorithm considers all possible permutations of the interactions and simultaneously optimizes the interactions of all pairs of heterologs. For example, if there is a high frequency of nuclei wherein a homolog of CT13 associates with CT15 and 21 while the other homolog of CT13 associates with CT22, it will classify the first as “CT13a” and the second as “CT13b” across all cells (Fig 1). This process is performed simultaneously for all CT studied to identify similarity across the population. Once this process is complete, we then determine how many input nuclei have an interaction for each of the possible pairwise combinations. After permutation analysis, the CM generates an output matrix in which each excel cell contains the percent of nuclei that have that given interaction (Fig 1b). Using Excel’s conditional formatting, each interaction is filled in with a color ranging from green (high) to red (low). Yellow indicates moderate values. By taking a threshold that is greater than randomizations of the input matrices, we can generate the most probable graph of interactions above random interactions which is the chromatin median graph (Fig 1b). The CM represents the state-of-the-art for determining the most common pattern among a population of cells and yields almost optimal solutions. It has been validated extensively by empirical comparison analysis on both random and real datasets with various data sizes (Ding et al., 2013).
In this communication we characterize the interchromosomal interactions of the NOR-CT and their associations with nucleolar domains in WI38 normal human diploid lung fibroblast cells. To reduce experimental variations, all studies were performed in G0 synchronized cells. Nucleoli were visualized by immunofluorescence labeling of the abundant nucleolar protein nucleolin and the boundaries of the cell nucleus were identified using DAPI staining of genomic DNA. The CT were labeled in a sequential series of re-FISH procedures. In some experiments, in addition to CT, the rDNA genes were labeled, as described in Materials and Methods.
Representative 3-D volume rendering images following this sequence of immunolabeling and re-FISH are shown in Figure 2. Accumulative mapping of fluorescence signals are illustrated where the nucleolus, the five pairs of NOR-CT (# 13, 14, 15, 21 and 22), rDNA and genomic DNA were stained in several steps in the same cell (Fig 2A). Figure 2B displays images of cells with different numbers of nucleoli and all 5 NOR-CT. These multiplex sets of imaging data obtained in 39 cells were then used for computerized analysis of the NOR-CT and nucleoli.
The number of nucleoli in WI38 cells varied from one to six with an average of ~3 nucleoli per cell. As expected of diploid genomes, chromosome painting indicated the presence of only two copies of each labeled NOR-CT in almost all studied cells (> 99%). Consistent with the localization of ribosomal genes within nucleoli (Fig 2A), the majority of the NOR-CT were in close apposition to the nucleoli. A small portion of the NOR-CT, however, were separated from the nucleoli by variable distances that extended up to 4 μm (Fig. 2C and E). Since the rDNA gene clusters on the p-arms of the NOR-CT are completely associated with the nucleolar interior with no detectable extranucleolar signal (Fig 2A, bottom right image), the NOR-CT that are more distal from the nearest nucleoli must also have their rDNA clusters associated with the interior of nucleoli. To allow specificity of chromosome labeling, highly repetitive sequences are removed from the chromosome specific DNA libraries during chromosome paint preparation. Large portions of the NOR-CT p-arms (rDNA clusters and other repetitive sequences) are removed from the paints. These paints therefore predominantly label the q-arms (76–86% of the NOR-CT; see http://genome.ucsc.edu/). It is likely that those NOR-CT more distal to nucleoli are actually in close apposition via connections to the p-arm regions of the CT that are not well visualized with this approach as reported previously (Kalmarova et al., 2007).
For quantitative analysis, the acquired images were segmented based on the signal threshold intensity. The 3D coordinates of the segmented regions and a large number of distance measurements were then determined with our suite of computer programs (see Material and Methods and ((Zeitz et al., 2009). By analyzing the 3D segments, we measured the average volume of the WI38 cell nucleus as ~1300 μm3. Within the nucleus, the average sizes of labeled NOR-CT ranged from ~14 μm3 for the smallest CT 21, to ~25 μm3 for the larger CT 13, 14, and 15 (Fig 3A). Altogether, the NOR-CT occupied about 16% of the overall nuclear volume and there was a strong linear correlation [r2 = 0.94], between mbp DNA lengths and their volumes (Fig 3B).
We documented a direct correlation between the volumes of the nucleoli and the nuclei in WI38 cells. While the total volume of the nucleoli decreased strikingly as the number of nucleoli in each cell increased (~2.4 fold between 1 to 5 nucleoli; Fig 3C), there was a corresponding decrease in total nuclear volume so that the overall nucleolar volume decreased only slightly when expressed as the % of total nuclear volume (8% to 6%, Fig 3C and D). Moreover, the total nucleolar surface area increased only slightly as the number of nucleoli increased (Fig 3E). Thus regardless of the number of nucleoli in different cells, the NOR-CT are associated with approximately the same nucleolar surface area for potential interactions. In addition, the size of the nucleoli did not correlate with the number of associated NOR-CT or with the intensity of the rDNA signal inside the nucleolar domains (Fig 2).
A number of previous studies have shown a linear correlation between chromosome size and their radial position where the smallest chromosomes occupy the nuclear interior and the largest chromosomes are closer to the nuclear periphery (Fritz et al., 2014; Marella et al., 2009a; Sun et al., 2000; Zeitz et al., 2009). It is remarkable to note that although the positions of all NOR-CT are influenced by their associations with nucleoli, a similar linear correlation was preserved for the radial positioning of these CT (Fig 3F). Previous studies of non-NOR-CT demonstrated that gene dense CT are more interior than gene poor CT (Boyle et al., 2001). This was dependent on the shape of nuclei (Cremer et al., 2001). However, we did not detect a correlation between the radial position of CT and gene density for NOR-CT (Fig 3G). For example, despite being of similar mbp length and volume, CT22 is twice the gene density of CT21, yet both share similar radial positioning (Fig 3F–G).
The inherent proximity of NOR-CT to nucleoli provides a convenient experimental model to study chromosomal neighborhood associations using re-FISH in combination with our computerized image analysis and computational geometric approaches. In one approach, the levels of interchromosomal associations as defined by a certain nearest distance threshold are determined for each CT pair. In a second approach the association of the NOR-CT with the nucleoli are measured using the same threshold method.
To determine NOR-CT-nucleolar associations, the nucleoli and NOR-CT were segmented and the nearest edge-to-edge 3-D distances were determined between individual CT and nucleoli. Nearest distances of ≤ 1.0 μm (14 pixels) were counted as positive associations. In cases where a CT was in association with more than one nucleolus, the association was assigned to the nucleolus with the largest contact area. Based on these criteria, we screened the frequency of pair-wise associations of specific NOR-CT with the same nucleolus. A large range of associations was measured for the heterologous CT pair from 51% of cells for the CT13-14 pair to 90% of cells for the CT15-22 pair (Fig 4A). Associations of homologous CT with the same nucleolus also showed a high frequency and wide range from 48% for the CT14 pair to 72 % for the CT21 pair (Fig 4B). In contrast, simulations where the NOR-CT are placed back into the nucleus randomly constrained (see Materials and Methods) to associate with the nucleoli (Fig 5A and B), revealed much lower levels of association and a much narrower range of values (26–38%, Fig 5C). All the individual experimental values were much higher (1.5–3.5 fold) than the values following simulation.
Next, we determined the number of CT associated with each nucleolar domain. Using the same criteria as above, we determined associations between individual NOR-CT and discrete nucleoli. We observed that NOR-CT do not associate with different nucleoli in equal or near-equal amounts as would be implied by a random pattern of associated CT. Instead we report a highly nonrandom phenomenon in the chromosome distribution pattern, where one of the nucleoli in each cell is always association with the greatest number of NOR-CT (Fig 2 and and4C).4C). Nucleoli associated with the greatest number of NOR-CT in each cell were termed “dominant” nucleoli. Those with the lowest and intermediate numbers of CT were termed “minor” and “intermediate” nucleoli, respectively. While the “dominant” nucleolus (indicated by orange arrows in Fig 2) is often the largest nucleolus in the cell (Fig 2D and E), this is not always the case (Figure 2A and F). Figure 2G shows a cell with only one nucleolus and all of the NOR-CT in close proximity.
Depending on the number of nucleoli, a wide range of nucleolar associations were found for individual NOR-CT (Fig 4D–E). In cells containing a single nucleolus, the association frequency ranged from 60% for CT14 to 82% for CT21 (Fig 4D). In cells with two nucleoli, CT13, CT14 and CT22 were predominantly associated with the dominant nucleolus while CT21 was predominantly associated with the intermediate nucleolus. With three or more nucleoli only CT13, CT14 and CT15 were preferential in their association with the dominant nucleolus. Thus while one nucleolus in each cell preferentially contains the greatest number of associated NOR-CT, the distribution of individual NOR-CT among the nucleolar population is more complex.
NOR-CT are defined as “neighbors” and assigned with a positive pair-wise association, when the nearest edge-to-edge 3-D distance between the segmented CT is ≤0.56 μm (8 pixels). This distance is large enough to distinguish boundaries between two homologous CT stained in the same fluorescence channel. A wide range of pairwise associations (Fig 6A,B) were determined for both heterologous (52–82%) and homologous (24–36%) CT pairs. Increasing the threshold distance between 8 to 20 pixels only slightly altered the levels of association (Fig S1). Importantly, the levels of pairwise associations at zero pixel distance (38–72%) were still very high (Fig S1). Moreover, the triplet cluster CT15-21-22 was observed in 72 % of the cells indicating a highly non-random association pattern.
Since the NOR-CT are predominantly confined to the crowded volume of the perinucleolar space, we applied our random constrained simulation program to determine to what extent these high associations values can be explained by a random association of the NOR-CT within this confined space. In this simulation (see Materials and Methods) the 5 pairs of NOR-CT are randomly placed into the nucleus with the constraint that each NOR-CT is positioned ≤ 1.0 μm of one or more of the nucleoli (Fig 5A–B). All of the association frequencies following this simulation were lower than the corresponding experimental results (Fig 5D). 6 of the 10 associations were 1.5–2.8 fold higher than the random simulation values while the remaining 4 associations displayed lower levels of differences (1.2 –1.3 fold).
In most previous studies, distances between centers of gravity have been used to study the organization of CT in the cell nucleus which may not be reflective of the interactions between CT (Fritz et al., 2014). Our analysis demonstrates that NOR-CT which are in close association as measured by the nearest edge-to-edge approach, have widely varying center-to-center distances (Fig 6C). Thus center-to-center distances, while providing important positional information, cannot accurately identify many of the close associations among the CT population in the nucleus.
As illustrated in Figure 6D, a CT pair may produce up to four pairwise interactions with another CT pair. We, therefore, determined the distribution for each NOR-CT pair of single and multiple pairwise interchromosomal associations (Fig 6E). While some CT pairwise associations consisted largely of singular interactions (CT13-14, CT13-22, CT14-15, CT14-21, CT15-21), others had approximately equal amounts of singular and multiple interactions (CT13-15, CT13-21, CT14-22, CT21-22). The great majority of the multiple interactions involved two interactions, while three and four interactions were infrequent and present mainly with the smaller CT21 and CT22 (Fig 6E).
The differences in the frequencies of pairwise associations (Fig 6A and B), the variations in distributions of singular and multiple interactions (Fig 6E), and the association of particular NOR-CT with nucleoli (Fig 4), suggest a complex and nonrandom pattern to the overall organization of NOR-CT in the cell nucleus. To investigate this further we used a novel computational data mining and pattern recognition program termed the chromatic median (CM). The CM is designed to determine the most probable overall pattern of CT associations across the population of nuclei (Ding et al. 2013). A simple illustration of this approach is presented in Figure 1 and is detailed in Materials and Methods. In brief, this program determines the corresponding homologs across all nuclei based upon that homolog’s interactions. It then represents each nucleus as a 10x10 (5 CT, 2 homologs per CT) binary matrix wherein a value of 1 indicates an interaction and a value of 0 indicates the absence of an interaction. Subsequently we determine the percent of nuclei that contain an interaction for each of the 45 positions in the matrices. Within these resulting median matrices, hot and cold-spots were found that ranged from 6 to 47% of input nuclei (Fig 7A) compared to 6–23% following randomizations of the input matrices (Fig 7B). This demonstrates that the process of determining corresponding homologs in a population of nuclei with a high degree of interactions does not artificially create a pattern of those interactions. Chi test analysis determined that median matrices from experimental and randomized input are significantly different from each other (p<0.001).
Using the CM algorithm we were able to determine corresponding homologs based on which other CT each homolog is interacting. Our results demonstrate that homolog a does interact with certain other CT with higher frequency than homolog b. These differences between homologs in experimental nuclei were greater than randomizations. For example, CT15a interacts with CT21a in 47% of nuclei, but CT15b and 21b interact in only 11%. This observation suggests that functional neighborhoods of interaction extend to the homolog level. This could be a result of certain NOR-CT interacting with particular nucleoli.
To determine a preferred probabilistic model of NOR-CT interactions, thresholding was performed on the matrices at 24% association. This enriched for the CT interactions found at the higher levels among the total population. Moreover at 24%, there are no connections in randomizations of input matrices or random simulations. The CM program revealed a global interactive network composed of 13 connections between the 5 pairs of NOR-CT (Fig 7C). Interestingly, CT21 had the largest number of connections (7) despite being the smallest human chromosome. CT14, CT15 and CT22 had 5 connections each, while CT13 had only 3. This lower level was likely due to the absence of any predicted connections at this probabilistic level for the CT13b homolog (Fig 7C). The connections with the highest % associations in the model are indicated as thick connecting lines as CT21-22 and CT15-21 (Fig 7C) Moreover, the two smallest chromosomes in the human genome CT21 and CT22 were the only NOR-CT displaying homologous connections. Sorenson’s analysis (Sorenson, 1948) revealed that each individual nucleus in the population analyzed contains on average 32% of the network connections displayed in this probabilistic model.
The role of the nucleolus in cell function extends beyond its central position as a ribosome factory (Busch and Smetana, 1970; Koberna et al., 2002). A plethora of functionally related properties have been ascribed to this sub-nuclear compartment which have largely baffled our understanding (Boisvert et al., 2007; Martelli et al., 2001; Olson et al., 2000; Pederson, 1998; Raska et al., 2006; Rubbi and Milner, 2003). Proteomics has demonstrated many hundreds of nucleolar proteins the great majority of which are not known to be involved in ribosomal RNA transcription or metabolism (Andersen et al., 2005; Pederson and Tsai, 2009). Genomics has elucidated nucleolar associated sequences (NADs) that are not limited to the NOR-bearing chromosomes but distributed throughout the genome (Nemeth et al., 2010; van Koningsbruggen et al., 2010). It is this latter property that captured our attention and led us to analyze the interchromosomal and nucleolar associations of the NOR-CT.
There is growing emphasis on deciphering the relationship of higher level genomic organization and function to the regulation of gene expression (Berezney, 2002; Berezney et al., 2005; Bickmore, 2013; Cremer and Cremer, 2010; Kumaran et al., 2008; Malyavantham et al., 2010; Meaburn and Misteli, 2007a; Misteli, 2007; Stein et al., 2003; Stein G, 2008; T, 2005; Zaidi et al., 2007). It is now well established that the mitotic chromosomes are arranged in the interphase cell nucleus as discrete 3-D structures termed chromosome territories (CT), (Bickmore, 2013; Cremer and Cremer, 2001; Cremer and Cremer, 2010; Cremer et al., 2006; Meaburn and Misteli, 2007b). Recent findings suggesting an involvement of intra- and inter-chromosomal interactions in gene expression and regulation (Berezney, 2002; Berezney et al., 2005; Bickmore, 2013; Clowney et al., 2012; Kumaran et al., 2008; Lanctot et al., 2007; Malyavantham et al., 2008; Malyavantham et al., 2010; Misteli, 2004; Misteli, 2007; Osborne et al., 2004; Osborne et al., 2007; Parada et al., 2004b; Spilianakis et al., 2005; Stein G, 2003; Stein G, 2008; T, 2005; Zaidi et al., 2007) has brought to the forefront the need for a better understanding of the 3-D organization of CT and their interactions in the cell nucleus. By combining the techniques of re-FISH with 3-D microscopy, computer image analysis and novel computational geometric data mining and pattern recognition algorithms, progress has been made in deciphering the 3-D arrangement of large subsets of human CT (Fritz et al., 2014; Marella et al., 2009a; Zeitz et al., 2009). In this study we use these tools to investigate the interchromosomal arrangements of the five pairs of nucleolar associated NOR-CT (CT 13,14,15, 21 and 22) in WI38 human diploid lung fibroblasts and the relationships of NOR-CT and nucleoli associations.
Analysis of NOR-CT associations with nucleoli in G0 synchronized WI38 fibroblasts demonstrated that one nucleolus was always associated with more NOR-CT than the other nucleoli (Fig 4). Smetena et. al. (1999, 2006) previously observed this in cells from human patients using the silver staining procedure for NORs and termed these the “dominant nucleoli” (Smetana et al., 2006a; Smetana et al., 1999; Smetana et al., 2006b). We further show that while many of the individual NOR-CT are predominantly associated with the dominant nucleolus, others are preferentially associated with nucleoli with intermediate levels of NOR-CT (Fig 4).
A wide range of high level associations (50–90%) were found for heterologous pairs of NOR-CT with the same nucleolus which exceeded those demonstrated by random simulations by 1.5–3.5 fold (Fig 4, ,5).5). There was also an unexpectedly high level of homologous NOR-CT pairs associated with the same nucleolus (48–72%). Random simulations showed a size dependency with the larger NOR-CT (e.g. CT 13–15) at higher association levels than the smaller CT (e.g., CT 21–22, Fig 5C). In contrast, the experimental data had no clear size dependency, e.g., CT 21–22, CT 21–21 and CT 22–22 showing much higher levels of association than CT 14-14 and 13-14 (Fig 4A, B). Together these results indicate the presence of non-random association patterns of NOR-CT with individual nucleoli. Further supporting this conclusion, it was reported that a high percentage of the NOR-CT are associated with nucleoli and that at least a portion of NOR-CT associated with the same nucleolus were preserved in daughter cells (Cvackova et al., 2009; Kalmarova et al., 2008b; Kalmarova et al., 2007).
These findings also have implications for nucleolar biogenesis in the cell cycle. It is well established that the nucleoli are formed in late telophase/ early G1-phase, around the sites of ribosomal genes expression, engulfing or fusing the material of several adjacent NORs protruded from their CT (Prieto and McStay, 2005; Raska et al., 2006). The discovery of high levels of specific NOR-CT pairs associated with the same nucleoli suggests that ordered clusters of certain combinations of NOR-CT may form by the time of exit from mitosis and subsequently contribute their NORs to the formation of the same nucleolus. In some cells this clustering involves all NOR-bearing chromosomes which consequently leads to formation of a single nucleolus. In other cells there is one larger cluster of chromosomes that gives a rise to a “dominant nucleolus” and several smaller sub-clusters or separate chromosomes that produce their own nucleoli (Fig 4). It is important to note that the nucleolar patterns measured in these investigations were performed in cells synchronized in G0. Thus differences in the cell cycle and particularly the observed phenomenon of fusion of smaller nucleoli into single larger ones that occur during the S and G2 phase (Sister Paula and Nardone, 1968; Wachtler et al., 1984; Wachtler and Stahl, 1993), do not contribute to these patterns.
By linking specific NOR-CT with non-NOR CT, the post-mitotic locations of NOR-CT and their associations with different nucleoli could have a major impact on the biogenesis of the overall chromosomal landscape of the cell nucleus. In this regard high levels of interchromosomal associations were detected between NOR-CT 15 and 21 and several non-NOR CT (Fritz et al., 2014). This is consistent with recent genomic findings suggesting that many CT aside from the NOR-CT associate with the nucleolus (Nemeth et al., 2010; van Koningsbruggen et al., 2010) where they presumably contribute to the perinucleolar heterochromatin region (Cvackova et al., 2009).
A high level and wide range (52–82%) were measured for heterologous interchromosomal pairwise associations (Fig 6) while homologous pairwise associations were much lower (24–36%). 6 of the 10 heterologous pairwise combinations were significantly above random simulations when the NOR-CT were rigorously constrained to positions surrounding the nucleoli (Fig 5). There were also major differences in the levels of multiple interactions among the CT pairs. For homologous associations, the smaller CT 21 and 22 pairs had significantly higher levels of association (Fig 4).
The above findings support a high degree of non-randomness to the arrangement of the NOR-CT. We, therefore, applied a novel pattern recognition algorithm termed the chromatic median (Ding et al., 2013) to determine if there was a preferred arrangement of interchromosomal associations for the entire population of NOR-CT. The analysis revealed an interconnected network of associations that encompassed all the individual NOR-CT except one copy of CT13 (Fig 7C). Interestingly, the highest levels of associations within this global interactive network involved the two smallest chromosomes CT 21 and 22. Global interactive networks of CT were previous reported for other CT subsets which showed cell type and tissue specificity (Marella et al., 2009a; Zeitz et al., 2009). Moreover, a recent study revealed major alterations in the global interactive network in malignant breast cancer cells (Fritz et al., 2014).
The analysis of NOR-CT positional associations not only provides valuable information about this subset of CT, but also represents an important approach for investigating the NOR-CT in relation to the rest of the genome. The NOR regions of these CT are anchored inside the nucleoli but the bulk volumes of these CT are located outside the nucleoli and can thus create pair-wise associations with different nucleolar and non-nucleolar CT. In this regard, a study of CT associations among a subset of 9 CT pairs including 2 NOR-CT, showed significant levels of association among the NOR and non-NOR CT which were also present in their global interactive networks (Fritz et al., 2014).
Recent genomic studies demonstrate that many if not all human chromosomes are associated with nucleoli via nucleolar associating domains (NADs) including a host of gene sequences unrelated to ribosomal RNA metabolism and ribosome biogenesis (Nemeth et al., 2010; van Koningsbruggen et al., 2010). This raises the possibility that the nucleolus may be directly involved in the regulation of a variety of genes other than rDNA (Nemeth et al., 2010; van Koningsbruggen et al., 2010). By sequestering inactive genes in the perinucleolar heterochromatin and protrusion of active genes into the extranucleolar regions for expression, the nucleolus may play a fundamental role in the global organization and genomic function of CT (van Koningsbruggen et al., 2010). Interestingly, NADs share partial sequence similarity with lamina-associated domains (LADs) and as a result some of this inactive chromatin may be repositioned to the nuclear periphery (Guelen et al., 2008; van Koningsbruggen et al., 2010).
The nucleoli and nuclear envelope may also be part of an overall architectural system in the cell nucleus for CT positioning and function. Previous studies have demonstrated a possible role of nucleoli and the nuclear periphery in constraining chromatin positions (Chubb et al., 2002). In this regard, the nuclear matrix was originally defined as a structural framework in the cell nucleus composed of residual components of the nucleoli, the peripheral nuclear envelope and a fibrogranular internal matrix (Berezney and Coffey, 1974; Berezney and Coffey, 1977; Berezney et al., 1995). Moreover, appropriate salt extraction of whole cells results in a remarkable preservation of chromosome territory organization but only if nuclear matrix organization is maintained (Ma et al., 1999). Treatments that result in nuclear matrix extraction lead to a corresponding disruption of chromosome territories (Ma et al., 1999).
In summary, our study identifies non-random associations for the subset of NOR-CT and demonstrates their arrangement into a global interactive network. Similar networks in other subsets of CT suggest that the NOR-CT are part of an overall genome wide network extending throughout the cell nucleus. We previously proposed that this interactive CT network may function as a probabilistic chromosome code that mediates genomic function at the global level (Fritz et al., 2014; Marella et al., 2009a; Zeitz et al., 2009). Deciphering the interactions of the NOR-CT in this genome wide network may, therefore, significantly contribute to our understanding of the multi-functional nucleolus.
This research was supported by grants from the National Institutes of Health (GM-072131) to R.B, the National Science Foundation (IIS-0713489 and IIS-1115220) to J.X. and R.B. and the University at Buffalo Foundation (9351115726) to R.B.