A software solution for recording circadian oscillator features in time-lapse live cell microscopy
© Sage et al; licensee BioMed Central Ltd. 2010
Received: 22 March 2010
Accepted: 6 July 2010
Published: 6 July 2010
Fluorescent and bioluminescent time-lapse microscopy approaches have been successfully used to investigate molecular mechanisms underlying the mammalian circadian oscillator at the single cell level. However, most of the available software and common methods based on intensity-threshold segmentation and frame-to-frame tracking are not applicable in these experiments. This is due to cell movement and dramatic changes in the fluorescent/bioluminescent reporter protein during the circadian cycle, with the lowest expression level very close to the background intensity. At present, the standard approach to analyze data sets obtained from time lapse microscopy is either manual tracking or application of generic image-processing software/dedicated tracking software. To our knowledge, these existing software solutions for manual and automatic tracking have strong limitations in tracking individual cells if their plane shifts.
In an attempt to improve existing methodology of time-lapse tracking of a large number of moving cells, we have developed a semi-automatic software package. It extracts the trajectory of the cells by tracking theirs displacements, makes the delineation of cell nucleus or whole cell, and finally yields measurements of various features, like reporter protein expression level or cell displacement. As an example, we present here single cell circadian pattern and motility analysis of NIH3T3 mouse fibroblasts expressing a fluorescent circadian reporter protein. Using Circadian Gene Express plugin, we performed fast and nonbiased analysis of large fluorescent time lapse microscopy datasets.
Our software solution, Circadian Gene Express (CGE), is easy to use and allows precise and semi-automatic tracking of moving cells over longer period of time. In spite of significant circadian variations in protein expression with extremely low expression levels at the valley phase, CGE allows accurate and efficient recording of large number of cell parameters, including level of reporter protein expression, velocity, direction of movement, and others. CGE proves to be useful for the analysis of widefield fluorescent microscopy datasets, as well as for bioluminescence imaging. Moreover, it might be easily adaptable for confocal image analysis by manually choosing one of the focal planes of each z-stack of the various time points of a time series.
CGE is a Java plugin for ImageJ; it is freely available at: http://bigwww.epfl.ch/sage/soft/circadian/.
Circadian oscillators have been described in virtually all organisms from cyanobacteria to humans. The mammalian circadian timing system has a hierarchical structure in that a master pacemaker residing in the suprachiasmatic nucleus synchronizes slave oscillators existing in most body cells . Moreover, circadian clocks are ticking in mammalian cultured cell lines, like Rat1 or NIH3T3 fibroblasts, and these clocks are self-sustained and cell-autonomous [2, 3]. A negative transcription/translation feedback loop, comprising clock genes repressing their own transcription, was proposed as the universal operational principle for generating circadian rhythm. Posttranslational events, like protein phosphorylation or acetylation, contribute critically to rhythm generation [4, 5].
Recent advances in time-lapse fluorescent imaging have allowed new insights into the mechanisms of circadian rhythms. Luciferase enzymes have been extensively used as reporters for numerous purposes in organisms as diverse as cyanobacteria, plants, fruit flies, and mice . Bioluminescence and fluorescence time lapse microscopy approaches have been successfully used to investigate molecular mechanisms of the mammalian circadian oscillator at a single cell level, the cross talk between individual cell clocks, and the mechanisms of single cell clock synchronization [3, 7]. Transgenic NIH3T3 cell lines stably expressing a short-lived nuclear yellow fluorescent protein (Venus) from circadian regulatory elements of the Rev-erbα locus (Rev-VNP), or luciferase protein driven by circadian Bmal1 promoter (Bmal1-luc), have been established and exploited to unravel different aspects of mammalian circadian clockwork machinery [3, 4, 8].
In spite of remarkable potential of the time lapse microscopy to address various questions of circadian biology, there is a very limited number of data analysis software available. Commercially available software Metamorph (Universal Imaging Corp), Imaris (Bitplane A.G.) and DiaTrack (Semasopht) incorporate modules to track objects and to measure intensity in a region of interest. However, the analysis of the reporter protein level in the described above time lapse microscopy datasets using these software requires a lot of manual interventions. Metamorph interrupts tracking in every valley of the circadian cycle; therefore the user has to manually complete the trace. This is mainly due to the high variation of intensity in the reporter protein level from one frame to another. Approaches based on intensity threshold or on template matching are not able to perform a correct tracking. In addition, a manual analysis is unreasonably time-consuming and subject to errors in observer judgment.
In an attempt to go beyond the tracking capability of conventional software, we tailored our approach towards tracking over longer periods of time. To achieve this, we had to employ advanced image-analysis methods to filter away reliance on a strongly changing fluorescent or bioluminescence reporter signal. We developed this new user-friendly image-analysis software for accurate tracking of individual cells in a living cell population. Tools presented here allow tracking and segmentation of the cells under the conditions of cyclic variations of intensities. The standard approach to track is to decompose the problem into two steps: 1) the segmentation phase which extracts the objects from the background in a frame; 2) the linking phase which tries to find the best match between objects from one frame to the next frame. This is the "frame-to-frame tracking" paradigm, taken by most commercial software and by the majority of the research community . A nearest-neighbor approach fails quickly because the cells have similar appearance in the valley of the circadian cycle. This approach can not be applied for the study of circadian oscillator. The circadian reporter protein expression level oscillates dramatically over 24 hours, to the extent that at the lowest point it comes close to the background level, making it difficult to distinguish the reporter level from the background. To resolve this tracking problem, we propose a solution where there is no explicit detection of the object to track. We formalize the tracking as an optimal process solving the shortest-path problem with the dynamic programming algorithm . Therefore, extracting the cell trajectory consists of finding a path in the spatio-temporal volume 2D + T by optimizing a cost function. Practically, the optimization minimizes error based on the intensity and on the displacement. This procedure which takes into account the past and the future of the particles is efficient for tracking dim particles. The program first tracks the center of gravity of a cell over the whole sequence of images. For the images that have a good enough resolution to precisely delineate the contour of the cells, we have incorporated a segmentation tool to our software based on the active rays' methods. The tool allows contouring of the convex cell shape based on knowledge of the center of the cell. Finally, the program makes measurements of several parameters like size, mean protein expression, displacement, etc. In addition, we have developed a granulometry application, which allows following of protein sub-nuclear localization during the cell cycle. Similarly to endogenous DNA ligase I, RFP-DNA ligase I construct is associated with the site of DNA replication (replication foci) during S phase. Thus RFP-DNA ligase I sub-nuclear distribution changes from punctuated (during S-phase) to dispersed, which makes it possible to easily follow cell cycle progression . Our granulometry application discriminates between condensed and dispersed (homogenous) states of the protein, allowing in the case of RFP-DNA ligase I protein to obtain the accurate time frame for S-phase, and thus to follow cell cycle progression.
The tools are written as Java plugin of the popular ImageJ image-processing software package (ImageJ: National Institutes of Health, Bethesda, MD, USA). Our experimental results show that this semi-automatic analysis method is reliable, reproducible and efficient for individual moving cell tracking, fluorescence or bioluminescent protein expression level quantification and cell trajectory analysis.
1. Tracking of cells with highly variable intensities
The task of automated tracking of moving objects has been studied extensively in the digital image analysis literature. In biology, tracking is of fundamental importance in cell motility studies. Different techniques have been proposed in the past for tracking cells or particles in the context of biological imaging . However, the common tracking paradigm, which consists of pre-processing the image, detecting the objects, and linking them from frame to frame, fails in the application to circadian expressed proteins. Due to extreme variations in reporter protein expression level over 24 hours, there is a high risk of loosing the tracked cell at the lowest phase. At this point the grey-level intensity (corresponding to the level of protein expression) often approaches the background value.
The cells are approximately round.
Each cell leaves a unique trace; in the case of cell division, the requirement is to track a single one of the daughter cells. The circadian rhythm study does not need to recover cell division lineage.
The intensity of a cell can approach that of the background for short periods of time.
Images are down sized so that single cells are reduced to particles of a few pixels. Size reduction provides a significant gain in terms of computational cost but, obviously, one loses in spatial precision. However, at this stage, a rough estimate of cell positions is sufficient.
Choice of start and an end point. DP is performed between user-specified start and end positions P1 (x 1) and Pn (x n). The user can interactively choose the cell to track and provide the necessary positions.
Introduction of a confinement area. It is assumed that the cell cannot exit the area between P1 and PN with a margin.
Maximum displacement constraint. It is additionally assumed that the cell displacement between consequent frames cannot exceed a certain step size Δ.
This cost function (which is to be minimized) includes two standard terms known as external and internal energies in snake terminology. The external energy term depends on the data and favors small variations of intensities, while the internal energy is a spatial regularization term that penalizes large variations of displacements. We additionally included a third term favoring intensities close to the intensity of the starting or end points. The three coefficients λext, λint, and λfix are weighting parameters that also serve as normalization factors. They can be tuned by the user to match the requirements of specific problems. The internal energy term is particularly important in the tracking of particles with low observability; it allows the algorithm to pursue the tracking according to a smooth-trajectory assumption, even when the cells in question have weak protein expressions.
DP is a highly efficient and robust algorithm for extracting spatio-temporal trajectories of particles (Fig. 1C), even in difficult cases where the particles have low intensities. The method is, however, adequate only for tracking individual particles, and does not straightforwardly extend to the joint tracking of several particles. For multiple-trajectory analysis, the user has to choose the cell to track, run the tracking process, and repeat the same operations for other cells.
2. Cell shape delineation
The external term attracts the curve towards the contour via the potential function, while the internal energy term favors smooth curves. The final curve is then represented by a cubic spline that interpolates the nodes Nk. In the implementation, images are preprocessed in order to reduce noise; Gaussian smoothing and a non-linear diffusion filter are applied to each frame in the sequence (Fig. 2B).
3. Cell Feature Measurement
Quantities such as the level of gene expression as well as certain dynamic variables (velocity, direction of movement, and variations in the direction) are among interesting features that can be computed from the information provided by the tracking module, for both bioluminescent and fluorescent images. Other features, such as the size of the cells and mean levels and standard deviations of gene expression, are calculated by also taking into account the shape information obtained using the shape delineation module. Sudden changes in the cell size (from large to small) can be used to detect cell divisions.
For a typical sequence of 200 frames, the software is able to track one cell in less than 200 ms, and to extract the shape in less than 2 s (Apple MacPro DualCore 2.66 MHz). The user has to enter few (around 5) positions by clicking on the cells to track. All the measurements are reported in a spreadsheet format and in graphical form.
Testing: Comparison of cell oscillation pattern analysis by CGE to the existing method
1. Measurement of cell motility, nucleus size and cell division timing
In addition to cell circadian oscillation profile, we assessed other output parameters of CGE which are useful for numerous purposes in cell biology, such as cell nucleus size, cell motility and cell division time. Natural variations in cell size are remarkable between wild type individual NIH3T3 fibroblasts. These variations seem to be even more pronounced in NIH3T3 cells transduced with Large T Antigen of SV40 virus (LTAgSV40; Additional file 2: Movie 2). We used CGE to analyze oscillation pattern in correlation with these cell nucleus size. As depicted in Figure 4D, single cell analysis of NIH3T3 Rev-VNP LTAgSV40 fibroblasts indicated that similarly to parental cell line, these cells exhibit strong circadian oscillation pattern. About 15% of transformed cells had unusually large nucleus size (nucleus square ~1500 pixels2; example cell 1 (1973 pixels2; Additional file 2: Movie 2). Interestingly, increase in the nucleus and cell size correlated in our example experiment with longer oscillation period length: 25.5 h ± 0.8 h in cells with nucleus size smaller then 600 pixels2 (395 ± 90 pixels2, n = 9), in comparison to 27.5 h ± 1.3 h in cells with nucleus larger then 1500 pixels2. Circadian phase might also be delayed in larger cells (see cell 1 profile at the small graph in Figure 4D). Thus LTAgSV40 over expression did not change cell oscillation pattern in most of the cells, however it might affect circadian period length and circadian phase in the cells harboring larger size. Larger datasets have to be analyzed to confirm a possible correlation between cell size and cell circadian pattern.
Cell motility was evaluated using two parameters: the distance and the movement angle change (direction evaluation) between each two subsequent time points. Using a small subset of cells (n = 20), we measured cell circadian phase, cell motility and cell division time (as demonstrated in an example cell in Fig. 4E). In a good agreement with previous publications , cell motility was increased after the cell division. No significant correlation was found between cell circadian phase and cell motility.
Thus, in addition to cell circadian oscillation profile, our software gives accurate evaluation of cell nucleus size, cell motility (distance and angle) and cell division time. It allows semi-automatic tracking and reliable analysis of various parameters during long term experiments implying large data sets of the cells exhibiting highly varying fluorescent protein levels and changing significantly their position during the experiment.
2. Single cell oscillation profile analysis using bioluminescence time lapse microscopy large data sets
Fluorescent time lapse studies expose cells to strong light excitation, and therefore may be prone to toxic effects when using shorter wavelengths over longer periods of time. In our hands treatment of mouse fibroblasts with different drugs in combination with cell excitation every 30 minutes biased the experimental analysis due to cell poor condition, and to significant cell death. Recently developed bioluminescent microscopy, absolutely non-toxic for the cells , might represent an ultimate solution for the toxicity problem in long term experiments. We have successfully applied CGE method in order to analyze large datasets obtained from bioluminescent time lapse microscopy. Using this approach, we have recently unraveled intriguing questions of transcription and temperature compensation of mammalian circadian oscillator . Time lapse microscopy of NIH3T3 Bmal1-luc fibroblasts was performed in this study, and at least 130 cells were tracked with subsequent analysis of their oscillation profiles using CGE. Thus in addition to fluorescent time lapse microscopy analysis, CGE could be used to analyze data sets obtained by bioluminescent time lapse microscopy.
3. Cell cycle analysis
An extra feature we have developed for CGE, granulometry, allows the accurate follow up of protein sub-nuclear localization, specifically the discrimination between condensed and diffused protein expression pattern. In order to evaluate S phase time frame, which corresponds to condensed expression phase of RFP-DNA Ligase I construct , time lapse microscopy of NIH3T3 cells transiently expressing RFP-DNA Ligase I was performed as described in Supplementary Methods. In examples presented at Additional file 3: Movie 3 and Figure 4F experiment, the cell number 1 enters S phase 15 hours after the start of recording, staying in S-phase about 6 hours, and proceeding to G2 phase 21 hours after the beginning recording. Cell 2 stays in S-phase between 16 and 22 hours, cell 3 - between 22 and 28 hours. In the future, we are planning to use this experimental system and CGE granulometry analysis to follow the cell cycle and circadian clock in the same cell.
Discussion and Conclusions
In conclusion, the CGE package allows one to study accurately and efficiently the different features and changes of cells with significantly varying locations and protein expression levels imaged over a period of several days. Our newly developed software proves to be more reliable, reproducible and efficient for individual cell circadian pattern quantification, in comparison to results obtained by different methods (MetaMorph). It is suitable for analysis of oscillation pattern from large data sets obtained by both fluorescence and bioluminescence time lapse microscopy. Experimental results herein and from our previous publications suggest that CGE allows accurate evaluation of fluorescent or bioluminescent reporter protein level and easy correlation between cell oscillation pattern and cell size, cell motility pattern and cell division timing. The granulometry feature provides accurate evaluation of protein sub-nuclear distribution, useful for cell cycle progression analysis, and allowing easy correlation between cell clock and cell cycle by the same experimental approach. More generally, our method offers a wide range of possibilities for cell tracking, cell size and cell motility analyses, for acquiring the protein distribution pattern and evaluation of cell cycle progression, making it useful in various aspects of cell biology research.
Cell lines and constructs
NIH3T3 cell line stably expressing short-lived nuclear fluorescent protein Venus-NLS-Pest1 driven from circadian Reverbα regulatory sequences (NIH3T3 RevVNP) was established as described in Nagoshi et al. . To obtain LargeT Antigen SV40 (LTAg SV40) expressing cells, NIH3T3 Rev-VNP fibroblasts were transduced with LTAgSV40 expressing lentiviral particles as described by P. Salmon and colleagues [16, 17] with subsequent clonal line selection. Parental and transformed NIH3T3 cells were maintained in DMEM supplemented with 10% FCS.
Fluorescence Time lapse Microscopy
For time lapse microscopy experiments, cells were plated in 35-mm glass bottom dishes (WillCo-dish, type 3522, WillCo Wells B.V.) and grown to confluence. After stimulating the cells with 100 nM dexamethasone for 30 minutes, the medium was replaced by 2 ml phenol red-free DMEM supplemented with 10% FCS. The cultures were placed in a 37°C chamber equilibrated with humidified air containing 5% CO2 throughout the microscopy. Time lapse microscopy was performed with a Leica AF6000LX microscope using a 20× objective. For the cell cycle follow up experiments, NIH3T3 cells were plated in 35-mm Falcon plastic dishes, transiently transfected with RFP-DNAI ligase , and microscopy was started 24 hours after transfection using 40× objective. Time lapse images were captured with either a Roper Coolsnap HQ or a Cascade B CCD camera using a YFP filter set (Additional file 1 and 2: Movie 1, 2) or RFP filter set (Additional file 3: Movie 3).
We thank Ueli Schibler for fruitful discussions and funding of this work, Hans Reinke, Gad Asher and David Suter for critical reading of the manuscript, Christoph Bauer and Jerome Bosset for assistance with time lapse microscopy, Nicolas Roggli for the artwork, Richard James and Sara Deakin for the language corrections. This research was supported by the Swiss National Foundation through individual research grants to Ueli Schibler and to CD, Ernst and Lucie Schmidheiny Foundation to CD, and the Louis Jeantet Foundation of Medicine to Ueli Schibler.
- Reppert SM: A colorful model of the circadian clock. Cell 2006, 124: 233–236. 10.1016/j.cell.2006.01.009View ArticlePubMedGoogle Scholar
- Balsalobre A, Damiola F, Schibler U: A serum shock induces circadian gene expression in mammalian tissue culture cells. Cell 1998, 93: 929–937. 10.1016/S0092-8674(00)81199-XView ArticlePubMedGoogle Scholar
- Nagoshi E, Saini C, Bauer C, Laroche T, Naef F, Schibler U: Circadian gene expression in individual fibroblasts cell-autonomous and self-sustained oscillators pass time to daughter cells. Cell 2004, 119: 693–705. 10.1016/j.cell.2004.11.015View ArticlePubMedGoogle Scholar
- Asher G, Gatfield D, Stratmann M, Reinke H, Dibner C, Kreppel F, Mostoslavsky R, Alt FW, Schibler U: SIRT1 regulates circadian clock gene expression through PER2 deacetylation. Cell 2008, 134: 317–328. 10.1016/j.cell.2008.06.050View ArticlePubMedGoogle Scholar
- Reinke H, Saini C, Fleury-Olela F, Dibner C, Benjamin IJ, Schibler U: Differential display of DNA-binding proteins reveals heat-shock factor 1 as a circadian transcription factor. Genes Dev 2008, 22: 331–345. 10.1101/gad.453808PubMed CentralView ArticlePubMedGoogle Scholar
- Welsh DK, Kay SA: Bioluminescence imaging in living organisms. Curr Opin Biotechnol 2005, 16: 73–78. 10.1016/j.copbio.2004.12.006View ArticlePubMedGoogle Scholar
- Liu AC, Welsh DK, Ko CH, Tran HG, Zhang EE, Priest AA, Buhr ED, Singer O, Meeker K, Verma IM, Doyle FJ, Takahashi JS, Kay SA: Intercellular coupling confers robustness against mutations in the SCN circadian clock network. Cell 2007, 129: 605–616. 10.1016/j.cell.2007.02.047PubMed CentralView ArticlePubMedGoogle Scholar
- Dibner C, Sage D, Unser M, Bauer C, d'Eysmond T, Naef F, Schibler U: Circadian gene expression is resilient to large fluctuations in overall transcription rates. EMBO J 2009, 28: 123–134. 10.1038/emboj.2008.262PubMed CentralView ArticlePubMedGoogle Scholar
- Cheezum MK, Walker WF, Guilford WH: Quantitative comparison of algorithms for tracking single fluorescent particles. Biophys J 2001, 81: 2378–2388. 10.1016/S0006-3495(01)75884-5PubMed CentralView ArticlePubMedGoogle Scholar
- Sage D, Neumann FR, Hediger F, Gasser SM, Unser M: Automatic tracking of individual fluorescence particles application to the study of chromosome dynamics. IEEE Trans Image Processing 2005, 14: 1372–1383. 10.1109/TIP.2005.852787View ArticleGoogle Scholar
- Easwaran HP, Leonhardt H, Cardoso MC: Cell cycle markers for live cell analyses. Cell Cycle 2005, 4: 453–455.View ArticlePubMedGoogle Scholar
- Meijering E, Smal I, Danuser G: Tracking in molecular bioimaging. IEEE Signal Processing Magazine 2006, 23: 46–53. 10.1109/MSP.2006.1628877View ArticleGoogle Scholar
- Chen Y, Huang TS, Rui Y: Optimal radial contour tracking by dynamic programming. In IEEE Int Conf Image Processing. Thessaloniki, Greece; 2001:626–629.Google Scholar
- Denzler J, Niemann H: Active rays polar-transformed active contours for real-time contour tracking. Real-Time Imaging 1999, 5: 203–213. 10.1006/rtim.1997.0116View ArticleGoogle Scholar
- Walmod PS, Hartmann-Petersen R, Prag S, Lepekhin EL, Röpke C, Berezin V, Bock E: Cell-cycle-dependent regulation of cell motility and determination of the role of Rac1. Exp Cell Res 2004, 295: 407–420. 10.1016/j.yexcr.2004.01.011View ArticlePubMedGoogle Scholar
- Salmon P, Oberholzer J, Occhiodoro T, Morel P, Lou J, Trono D: Reversible immortalization of human primary cells by lentivector-mediated transfer of specific genes. Mol Ther 2000, 2: 404–414. 10.1006/mthe.2000.0141View ArticlePubMedGoogle Scholar
- Salmon P, Trono D: Production and titration of lentiviral vectors. Curr Protoc Hum Genet 2007.,Chapter 12(Unit 12 10):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.