Spatial quantitative analysis of fluorescently labeled nuclear structures: problems, methods, pitfalls

The vast majority of microscopic data in biology of the cell nucleus is currently collected using fluorescence microscopy, and most of these data are subsequently subjected to quantitative analysis. The analysis process unites a number of steps, from image acquisition to statistics, and at each of these steps decisions must be made that may crucially affect the conclusions of the whole study. This often presents a really serious problem because the researcher is typically a biologist, while the decisions to be taken require expertise in the fields of physics, computer image analysis, and statistics. The researcher has to choose between multiple options for data collection, numerous programs for preprocessing and processing of images, and a number of statistical approaches. Written for biologists, this article discusses some of the typical problems and errors that should be avoided. The article was prepared by a team uniting expertise in biology, microscopy, image analysis, and statistics. It considers the options a researcher has at the stages of data acquisition (choice of the microscope and acquisition settings), preprocessing (filtering, intensity normalization, deconvolution), image processing (radial distribution, clustering, co-localization, shape and orientation of objects), and statistical analysis.

First steps towards systems radiation biology studies concerned with DNA and chromosome structure within living cells

For the understanding of radiation action on biological systems like cellular macromolecules (e.g., DNA in its higher structures) a synergistic approach of experiments and quantitative modelling of working hypotheses is necessary. Further on, the influence on calculated results of certain assumptions in such working hypotheses must critically be evaluated. In the present work, this issue is highlighted in two aspects for the case of DNA damage in single cells. First, yields of double-strand breaks and frequency distributions of DNA fragment lengths after ion irradiation were calculated using different assumptions on the DNA target model. Compared to a former target model now a moderate effect due to the inclusion of a spherical chromatin domain model has been found. Second, the influence of assumptions on particular geometric chromosome models on calculated chromosome aberration data is illustrated with two target-modelling approaches for this end point.

Nanostructure analysis using spatially modulated illumination microscopy

We describe the usage of the spatially modulated illumination (SMI) microscope to estimate the sizes (and/or positions) of fluorescently labeled cellular nanostructures, including a brief introduction to the instrument and its handling. The principle setup of the SMI microscope will be introduced to explain the measures necessary for a successful nanostructure analysis, before the steps for sample preparation, data acquisition and evaluation are given. The protocol starts with cells already attached to the cover glass. The protocol and duration outlined here are typical for fixed specimens; however, considerably faster data acquisition and in vivo measurements are possible.

Light optical precision measurements of the active and inactive Prader-Willi syndrome imprinted regions in human cell nuclei

Despite the major advancements during the last decade with respect to both knowledge of higher order chromatin organization in the cell nucleus and the elucidation of epigenetic mechanisms of gene control, the true three-dimensional (3D) chromatin structure of endogenous active and inactive gene loci is not known. The present study was initiated as an attempt to close this gap. As a model case, we compared the chromatin architecture between the genetically active and inactive domains of the imprinted Prader-Willi syndrome (PWS) locus in human fibroblast and lymphoblastoid cell nuclei by 3D fluorescence in situ hybridization and quantitative confocal laser scanning microscopy. The volumes and 3D compactions of identified maternal and paternal PWS domains were determined in stacks of light optical serial sections using a novel threshold-independent approach. Our failure to detect volume and compaction differences indicates that possible differences are below the limits of light optical resolution. To overcome this limitation, spectral precision distance microscopy, a method of localization microscopy at the nanometer scale, was used to measure 3D distances between differentially labeled probes located both within the PWS region and in its neighborhood. This approach allows the detection of intranuclear differences between 3D distances down to about 70-90 nm, but again did not reveal clearly detectable differences between active and inactive PWS domains. Despite this failure, a comparison of the experimental 3D distance measurements with computer simulations of chromatin folding strongly supports a non-random higher order chromatin configuration of the PWS locus and argues against 3D configurations based on giant chromatin loops. Our results indicate that the search for differences between endogenous active and inactive PWS domains must be continued at still smaller scales than hitherto possible with conventional light microscopic procedures. The possibilities to achieve this goal are discussed.

Chromosome neighborhood composition determines translocation outcomes after exposure to high-dose radiation in primary cells

Radiation exposure is an occupational hazard for military personnel, some health care professionals, airport security screeners, and medical patients, with some individuals at risk for acute, high-dose exposures. Therefore, the biological effects of radiation, especially the potential for chromosome damage, are major occupational and health concerns. However, the biophysical mechanisms of chromosome instability subsequent to radiation-induced DNA damage are poorly understood. It is clear that interphase chromosomes occupy discrete structural and functional subnuclear domains, termed chromosome territories (CT), which may be organized into 'neighborhoods' comprising groups of specific CTs. We directly evaluated the relationship between chromosome positioning, neighborhood composition, and translocation partner choice in primary lymphocytes, using a cell-based system in which we could induce multiple, concentrated DNA breaks via high-dose irradiation. We critically evaluated mis-rejoining profiles and tested whether breaks occurring nearby were more likely to fuse than breaks occurring at a distance. We show that CT neighborhoods comprise heterologous chromosomes, within which inter-CT distances directly relate to translocation partner choice. These findings demonstrate that interphase chromosome arrangement is a principal factor in genomic instability outcomes in primary lymphocytes, providing a structural context for understanding the biological effects of radiation exposure, and the molecular etiology of tumor-specific translocation patterns.

Cell type-specific quantitative predictions of radiation-induced chromosome aberrations: a computer model approach

A quantitative computer model was applied to simulate the three-dimensional (3D) spatial organization of chromatin in human cell nuclei under defined conditions of virtual irradiation to explore the implications of spatial organization on chromosome aberrations. To calibrate the virtual irradiation algorithm, a dose-dependent spectrum of radiation-induced chromosome aberrations such as dicentrics, translocations and centric rings was calculated for low-LET radiation doses ranging from 0.5 to 5 Gy. This was compared with the results from experimental studies. While the dose-response curves calculated from model simulations agree well with experimental dose-response curves for dicentrics and translocations, centric rings are significantly more frequent in the model simulation than in experiments despite taking into account exclusive arm territories in the applied Spherical 1 Mbp Chromatin Domain (SCD) computer model explicitly. Taking into account the non-random positioning of chromosome territories observed in lymphocyte cell nuclei (a so-called gene density-correlated arrangement of chromosome territories), aberration frequencies were calculated with the calibrated irradiation algorithm to investigate the impact of chromosome territory neighborhood effects (proximity effects). The absolute frequencies of pairwise exchanges agree well with those found in an experimental study. In conclusion, the results obtained using the computer model approach presented here based on only a few adjustable parameters correlated well with those of experimental studies of chromosome aberration frequencies. Thus the model may be a useful tool in radiation-induced cancer risk estimates in combination with epidemiological studies.

Comparison of triple helical COMBO-FISH and standard FISH by means of quantitative microscopic image analysis of abl/bcr positions in cell nuclei

In this study, a novel DNA fluorescence labelling technique, called triple helical COMBO-FISH (Combinatorial Oligo Fluorescence In Situ Hybridisation), was compared to the standard FISH (Fluorescence In Situ Hybridisation by means of commercially available probe kits) by quantitative evaluation of the nuclear position of the hybridisation signals of the Abelson murine leukaemia (abl) region and the breakpoint cluster region (bcr) in 3D-conserved cell nuclei of lymphocytes and CML blood cells. Two sets of 31 homopyrimidine oligonucleotides each, corresponding to co-localising sequences in the abl region of chromosome 9 and in the bcr region of chromosome 22 were synthesised. Probe types and sizes (in bases) as well as the binding mechanisms of both FISH techniques were completely different. In accordance to established findings that cell type specific radial positioning of chromosomes and sub-chromosomal elements is evolutionarily conserved, no significant difference was found between the two FISH techniques for the radial localisation of the barycentre of the analysed genomic loci. Thermal denaturation and hypotonic treatment of cell nuclei subjected to standard FISH, however, led to different absolute radii and volumes of the cell nuclei, in comparison to the quantities determined for the triple helical COMBO-FISH technique; the chromatin appears to shrink in laterally enlarged, flat nuclei. Consequently, the absolute distances of the homologous labelled sites shifted to greater values. For precise quantitative microscopic analysis of genomic loci, fluorescence labelling procedures are recommended that well maintain the native chromatin topology. Triple helical COMBO-FISH may offer such an approach.

Chromatin domains and the interchromatin compartment form structurally defined and functionally interacting nuclear networks

In spite of strong evidence that the nucleus is a highly organized organelle, a consensus on basic principles of the global nuclear architecture has not so far been achieved. The chromosome territory-interchromatin compartment (CT-IC) model postulates an IC which expands between chromatin domains both in the interior and the periphery of CT. Other models, however, dispute the existence of the IC and claim that numerous chromatin loops expand between and within CTs. The present study was undertaken to resolve these conflicting views. (1) We demonstrate that most chromatin exists in the form of higher-order chromatin domains with a compaction level at least 10 times above the level of extended 30 nm chromatin fibers. A similar compaction level was obtained in a detailed analysis of a particularly gene-dense chromosome region on HSA 11, which often expanded from its CT as a finger-like chromatin protrusion. (2) We further applied an approach which allows the experimental manipulation of both chromatin condensation and the width of IC channels in a fully reversible manner. These experiments, together with electron microscopic observations, demonstrate the existence of the IC as a dynamic, structurally distinct nuclear compartment, which is functionally linked with the chromatin compartment.

Rise, fall and resurrection of chromosome territories: a historical perspective. Part II. Fall and resurrection of chromosome territories during the 1950s to 1980s. Part III. Chromosome territories and the functional nuclear architecture: experiments and models from the 1990s to the present

Part II of this historical review on the progress of nuclear architecture studies points out why the original hypothesis of chromosome territories from Carl Rabl and Theodor Boveri (described in part I) was abandoned during the 1950s and finally proven by compelling evidence forwarded by laser-uv-microbeam studies and in situ hybridization experiments. Part II also includes a section on the development of advanced light microscopic techniques breaking the classical Abbe limit written for readers with little knowledge about the present state of the theory of light microscopic resolution. These developments have made it possible to perform 3D distance measurements between genes or other specifically stained, nuclear structures with high precision at the nanometer scale. Moreover, it has become possible to record full images from fluorescent structures and perform quantitative measurements of their shapes and volumes at a level of resolution that until recently could only be achieved by electron microscopy. In part III we review the development of experiments and models of nuclear architecture since the 1990s. Emphasis is laid on the still strongly conflicting views about the basic principles of higher order chromatin organization. A concluding section explains what needs to be done to resolve these conflicts and to come closer to the final goal of all studies of the nuclear architecture, namely to understand the implications of nuclear architecture for nuclear functions.

4Pi microscopy deconvolution with a variable point-spread function

To remove the axial sidelobes from 4Pi images, deconvolution forms an integral part of 4Pi microscopy. As a result of its high axial resolution, the 4Pi point spread function (PSF) is particularly susceptible to imperfect optical conditions within the sample. This is typically observed as a shift in the position of the maxima under the PSF envelope. A significantly varying phase shift renders deconvolution procedures based on a spatially invariant PSF essentially useless. We present a technique for computing the forward transformation in the case of a varying phase at a computational expense of the same order of magnitude as that of the shift invariant case, a method for the estimation of PSF phase from an acquired image, and a deconvolution procedure built on these techniques.

Folding and organization of a contiguous chromosome region according to the gene distribution pattern in primary genomic sequence

Specific mammalian genes functionally and dynamically associate together within the nucleus. Yet, how an array of many genes along the chromosome sequence can be spatially organized and folded together is unknown. We investigated the 3D structure of a well-annotated, highly conserved 4.3-Mb region on mouse chromosome 14 that contains four clusters of genes separated by gene “deserts.” In nuclei, this region forms multiple, nonrandom “higher order” structures. These structures are based on the gene distribution pattern in primary sequence and are marked by preferential associations among multiple gene clusters. Associating gene clusters represent expressed chromatin, but their aggregation is not simply dependent on ongoing transcription. In chromosomes with aggregated gene clusters, gene deserts preferentially align with the nuclear periphery, providing evidence for chromosomal region architecture by specific associations with functional nuclear domains. Together, these data suggest dynamic, probabilistic 3D folding states for a contiguous megabase-scale chromosomal region, supporting the diverse activities of multiple genes and their conserved primary sequence organization.

Rise, fall and resurrection of chromosome territories: a historical perspective. Part I. The rise of chromosome territories

It is now generally accepted that chromosomes in the cell nucleus are organized in distinct domains, first called chromosome territories in 1909 by the great cytologist Theodor Boveri. Yet, even today chromosomes have remained enigmatic individuals, whose structures, arrangements and functions in cycling and post-mitotic cells still need to be explored in full detail. Whereas numerous recent reviews describe present evidence for a dynamic architecture of chromosome territories and discuss the potential significance within the functional compartmentalization of the nucleus, a comprehensive historical account of this important concept of nuclear organization was lacking so far. Here, we describe the early rise of chromosome territories within the context of the discovery of chromosomes and their fundamental role in heredity, covering a period from the 1870th to the early 20th century (part I, this volume). In part II (next volume) we review the abandonment of the chromosome territory concept during the 1950th to 1980th and the compelling evidence, which led to its resurrection during the 1970th to 1980th.

COMBO-FISH for focussed fluorescence labelling of gene domains: 3D-analysis of the genome architecture of abl and bcr in human blood cells

Structural analysis and nanosizing of gene domains requires not only high-resolution microscopy but also improved techniques of fluorescence labelling strongly focussed on the gene domains. To investigate the architecture of abl and bcr in blood cell nuclei forming the Philadelphia chromosome in CML, we applied COMBO-FISH using specifically colocalising combinations of triple strand forming oligonucleotide probes for abl on chromosome 9 and bcr on chromosome 22. Each probe set consisting of 31 homopyrimidine oligonucleotides was computer selected from the human genome database. Measurements by 3D microscopy were compared to results obtained after standard FISH using commercially available abl/bcr BAC probes. The relative radial fluorescence distributions in lymphocyte cell nuclei of healthy donors in comparison to cell nuclei of blood cells of CML patients showed a strong correlation in the location of abl and bcr for both labelling techniques. The absolute distances of the homologous bcr domains and the abl domain-nuclear center-abl domain angles in cell nuclei of CML donors differed significantly from those of healthy donors only when COMBO-FISH was applied. These results indicate that COMBO-FISH may be more sensitive than standard FISH in case of slight modifications in the genome architecture.

Nanostructure of specific chromatin regions and nuclear complexes

Spatially modulated illumination (SMI) microscopy is a method of widefield fluorescence microscopy featuring interferometric illumination, which delivers structural information about nanoscale features in fluorescently labeled cells. Using this approach, structural changes in the context of gene activation and chromatin remodeling may be revealed. In this paper we present the application of SMI microscopy to size measurements of the 7q22 gene region, giving us a size estimate of 105+/-16 nm which corresponds to an average compaction ratio of 1:324. The results for the 7q22 domain are compared with the previously measured sizes of other fluorescently labeled gene regions, and to those obtained for transcription factories. The absence of a correlation between the measured and genomic sizes of the various gene regions indicate that a high variability in chromatin folding is present, with factors other than the sequence length contributing to the chromatin compaction. Measurements of the 7q22 region in different preparations and at different excitation wavelengths show a good agreement, thus demonstrating that the technique is robust when applied to biological samples.

Three-dimensional maps of all chromosomes in human male fibroblast nuclei and prometaphase rosettes

Studies of higher-order chromatin arrangements are an essential part of ongoing attempts to explore changes in epigenome structure and their functional implications during development and cell differentiation. However, the extent and cell-type-specificity of three-dimensional (3D) chromosome arrangements has remained controversial. In order to overcome technical limitations of previous studies, we have developed tools that allow the quantitative 3D positional mapping of all chromosomes simultaneously. We present unequivocal evidence for a probabilistic 3D order of prometaphase chromosomes, as well as of chromosome territories (CTs) in nuclei of quiescent (G0) and cycling (early S-phase) human diploid fibroblasts (46, XY). Radial distance measurements showed a probabilistic, highly nonrandom correlation with chromosome size: small chromosomes-independently of their gene density-were distributed significantly closer to the center of the nucleus or prometaphase rosette, while large chromosomes were located closer to the nuclear or rosette rim. This arrangement was independently confirmed in both human fibroblast and amniotic fluid cell nuclei. Notably, these cell types exhibit flat-ellipsoidal cell nuclei, in contrast to the spherical nuclei of lymphocytes and several other human cell types, for which we and others previously demonstrated gene-density-correlated radial 3D CT arrangements. Modeling of 3D CT arrangements suggests that cell-type-specific differences in radial CT arrangements are not solely due to geometrical constraints that result from nuclear shape differences. We also found gene-density-correlated arrangements of higher-order chromatin shared by all human cell types studied so far. Chromatin domains, which are gene-poor, form a layer beneath the nuclear envelope, while gene-dense chromatin is enriched in the nuclear interior. We discuss the possible functional implications of this finding.

Human archival tissues provide a valuable source for the analysis of spatial genome organization

Sections from archival formalin-fixed, paraffin wax-embedded human tissues are a valuable source for the study of the nuclear architecture of specific tissue types in terms of the three-dimensional spatial positioning and architecture of chromosome territories and sub-chromosomal domains. Chromosome painting, centromeric, and locus-specific probes were hybridized to tissue microarrays prepared from formalin-fixed paraffin wax-embedded samples of pancreas and breast. The cell nuclei were analyzed using quantitative three-dimensional image microscopy. The results obtained from non-neoplastic pancreatic cells of randomly selected individuals indicated that the radial arrangement of the chromosome 8 territories as well as their shape (roundness) did not significantly differ between the individuals and were in accordance with assumptions of a probabilistic model for computer simulations. There were considerable differences between pancreatic tumor and non-neoplastic cells. In non-neoplastic ductal epithelium of the breast there was a larger, but insignificant, variability in the three-dimensional positioning of the centromere 17 and HER2 domains between individuals. In neoplastic epithelial breast cells, however, the distances between centromere and gene domains were, on average, smaller than in non-neoplastic cells. In conclusion, our results demonstrate the feasibility of studying the genome architecture in archival, formalin-fixed, paraffin wax-embedded human tissues, opening new directions in tumor research and cell classification.

Nano-sizing of specific gene domains in intact human cell nuclei by spatially modulated illumination light microscopy

Although light microscopy and three-dimensional image analysis have made considerable progress during the last decade, it is still challenging to analyze the genome nano-architecture of specific gene domains in three-dimensional cell nuclei by fluorescence microscopy. Here, we present for the first time chromatin compaction measurements in human lymphocyte cell nuclei for three different, specific gene domains using a novel light microscopic approach called Spatially Modulated Illumination microscopy. Gene domains for p53, p58, and c-myc were labeled by fluorescence in situ hybridization and the sizes of the fluorescence in situ hybridization "spots" were measured. The mean diameters of the gene domains were determined to 103 nm (c-myc), 119 nm (p53), and 123 nm (p58) and did not correlate to the genomic, labeled sequence length. Assuming a spherical domain shape, these values would correspond to volumes of 5.7 x 10(-4) microm(3) (c-myc), 8.9 x 10(-4) microm(3) (p53), and 9.7 x 10(-4) microm(3) (p58). These volumes are approximately 2 orders of magnitude smaller than the diffraction limited illumination or observation volume, respectively, in a confocal laser scanning microscope using a high numerical aperture objective lens. By comparison of the labeled sequence length to the domain size, compaction ratios were estimated to 1:129 (p53), 1:235 (p58), and 1:396 (c-myc). The measurements demonstrate the advantage of the SMI technique for the analysis of gene domain nano-architecture in cell nuclei. The data indicate that chromatin compaction is subjected to a large variability which may be due to different states of genetic activity or reflect the cell cycle state.

Generic features of tertiary chromatin structure as detected in natural chromosomes

Knowledge of tertiary chromatin structure in mammalian interphase chromosomes is largely derived from artificial tandem arrays. In these model systems, light microscope images reveal fibers or beaded fibers after high-density targeting of transactivators to insertional domains spanning several megabases. These images of fibers have lent support to chromonema fiber models of tertiary structure. To assess the relevance of these studies to natural mammalian chromatin, we identified two different approximately 400-kb regions on human chromosomes 6 and 22 and then examined light microscope images of interphase tertiary chromatin structure when the regions were transcriptionally active and inactive. When transcriptionally active, these natural chromosomal regions elongated, yielding images characterized by a series of adjacent puncta or "beads", referred to hereafter as beaded images. These elongated structures required transcription for their maintenance. Thus, despite marked differences in the density and the mode of transactivation, the natural and artificial systems showed similarities, suggesting that beaded images are generic features of transcriptionally active tertiary chromatin. We show here, however, that these images do not necessarily favor chromonema fiber models but can also be explained by a radial-loop model or even a simple nucleosome affinity, random-chain model. Thus, light microscope images of tertiary structure cannot distinguish among competing models, although they do impose key constraints: chromatin must be clustered to yield beaded images and then packaged within each cluster to enable decondensation into adjacent clusters.

Radial arrangement of chromosome territories in human cell nuclei: a computer model approach based on gene density indicates a probabilistic global positioning code

Numerous investigations in the last years focused on chromosome arrangements in interphase nuclei. Recent experiments concerning the radial positioning of chromosomes in the nuclear volume of human and primate lymphocyte cells suggest a relationship between the gene density of a chromosome territory (CT) and its distance to the nuclear center. To relate chromosome positioning and gene density in a quantitative way, computer simulations of whole human cell nuclear genomes of normal karyotype were performed on the basis of the spherical 1 Mbp chromatin domain model and the latest data about sequence length and gene density of chromosomes. Three different basic assumptions about the initial distribution of chromosomes were used: a statistical, a deterministic, and a probabilistic initial distribution. After a simulated decondensation in early G1, a comparison of the radial distributions of simulated and experimentally obtained data for CTs Nos. 12, 18, 19, and 20 was made. It was shown that the experimentally observed distributions can be fitted better assuming an initial probabilistic distribution. This supports the concept of a probabilistic global gene positioning code depending on CT sequence length and gene density.

Virtual radiation biophysics: implications of nuclear structure

The non-random positioning of chromosome territories (CTs) in lymphocyte cell nuclei has raised the question whether systematic chromosome-chromosome associations exist which have significant influence on interchange rates. In such a case the spatial proximity of certain CTs or even of clusters of CTs is expected to increase the respective exchange yields significantly, in comparison to a random association of CTs. In the present study we applied computer simulated arrangements of CTs to calculate interchange frequencies between all heterologous CT pairs, assuming a uniform action of the molecular repair machinery. For the positioning of CTs in the virtual nuclear volume we assumed a) a statistical, and b) a gene density-correlated arrangement. The gene density-correlated arrangement regards the more experimentally observed interior localization of gene-rich and the more peripheral positioning of gene-poor CTs. Regarding one-chromosome yields, remarkable differences for single CTs were observed taking into account the gene density-correlated distribution of CTs.