Development and validation of 'AutoRIF': software for the automated analysis of radiation-induced foci
© McVean et al; licensee BioMed Central Ltd. 2012
Received: 27 September 2011
Accepted: 26 January 2012
Published: 26 January 2012
The quantification of radiation-induced foci (RIF) to investigate the induction and subsequent repair of DNA double strands breaks is now commonplace. Over the last decade systems specific for the automatic quantification of RIF have been developed for this purpose, however to ask more mechanistic questions on the spatio-temporal aspects of RIF, an automated RIF analysis platform that also quantifies RIF size/volume and relative three-dimensional (3D) distribution of RIF within individual nuclei, is required.
A java-based image analysis system has been developed (AutoRIF) that quantifies the number, size/volume and relative nuclear locations of RIF within 3D nuclear volumes. Our approach identifies nuclei using the dynamic Otsu threshold and RIF by enhanced Laplacian filtering and maximum entropy thresholding steps and, has an application 'batch optimisation' process to ensure reproducible quantification of RIF. AutoRIF was validated by comparing output against manual quantification of the same 2D and 3D image stacks with results showing excellent concordance over a whole range of sample time points (and therefore range of total RIF/nucleus) after low-LET radiation exposure.
This high-throughput automated RIF analysis system generates data with greater depth of information and reproducibility than that which can be achieved manually and may contribute toward the standardisation of RIF analysis. In particular, AutoRIF is a powerful tool for studying spatio-temporal relationships of RIF using a range of DNA damage response markers and can be run independently of other software, enabling most personal computers to perform image analysis. Future considerations for AutoRIF will likely include more complex algorithms that enable multiplex analysis for increasing combinations of cellular markers.
The use of markers, such as the phosphorylated variant of histone 2A (γ-H2AX), for the detection of radiation-induced double strand breaks (DSBs) is now standard practice for investigating biologically relevant doses of radiation (for a review of the field see Lobrich et al 2010 ). For instance, the induction and subsequent repair kinetics of DSB has been determined in a range of cell types upon exposure to varying qualities and radiation doses [2–6]. Additionally, the development of antibodies specific for other relevant proteins in the DNA damage response (DDR) pathway have revolutionised our ability to investigate mechanistic aspects of DSB processing in interphase nuclei. For example, 53BP1 has been shown to localise to sites of DSB [7–9] and studies have implicated the localisation and retention of 53BP1 to be γ-H2AX dependant , Mdc1 dependant [11–13] and as a function of chromatin accessibility . Compared to γ-H2AX, 53BP1 has a significantly improved signal to noise ratio and additionally has been applied for live cell visualisation of DSBs through generation of GFP fusion proteins [11, 13, 15, 16]. Radiation-induced foci (RIF) assays are now being applied to assess clinical outcome to radiation  as well as evaluating risk complications associated with radiotherapies and diagnostics [18–20]. RIF are also being exploited as dosimetry biomarkers to identify individuals exposed to unknown levels of radiation  and for the assessment of cellular sensitivity, which could be utilised as to screen patients prior to radiotherapy or diagnostics .
To reliably quantify RIF there must be a clear definition of what actually constitutes a 'focus'. A RIF can be described as a peak of signal intensity, distinct from the background and therefore displaying a strong signal to noise (S:N) ratio. Factors that can reduce this S:N ratio include poor immunofluorescence staining as a consequence of the antibody employed or poor technique, auto-fluorescence within the sample, the number and distribution of RIF (e.g whereby individual RIF overlap) and the optical system used. Counting of RIF has typically been achieved by manual analysis through the optical binoculars of a fluorescence microscope or by taking images and counting the RIF on a digital screen [23–27]. If the latter, then additional factors which may influence the quantification of RIF include (1) characteristics of the objective lens, (2) the camera used to acquire the images, and (3) the way in which the images are subsequently analysed. Thus, the type of image acquisition and the numerical aperture of the objective lens used will directly affect the resolution of RIF, potentially impacting on the number of RIF counted. For instance, if a single image of one focal plane is acquired using a high numerical aperture (N.A. 1.4) lens, then the resulting narrow field of view will mean only a small proportion of the nucleus is actually sampled. This can be resolved by capturing images from multiple focal planes for subsequent analysis as sequential 2D images or as collapsed maximum intensity projections (MIP). The advantage of analysing from an MIP compared to sequential 2D slices are improved S:N, which is particularly useful when manually analysing from digital images, however overlapping RIF from separate focal planes will be visualised as the same single RIF potentially leading to inaccuracies when nuclei contain large numbers of RIF crossing multiple focal planes. With regards to the choice of camera, the ability to resolve RIF will be determined by the resolution of digital images acquired and therefore the sensitivity and range (bit depth) of the camera employed. Accordingly, a minimum of an 8-bit camera is essential to ensure optimal acquisition of all RIF. Regardless of any technological variations for image acquisition however, the principle limiting factor to achieve high-throughput, reproducible and accurate RIF quantification are the difficulties in manually discriminating between background and RIF, particularly when RIF lie on multiple focal planes, and which ultimately lead to varying results between operators and between labs.
A number of automated systems that lend themselves for the quantification of RIF are available commercially. Image analysis packages such as Imaris, ImageJ and CellProfiler, typically rely upon the application of an intensity threshold and minimum size parameters to isolate RIF from the background. Over the last decade systems specific for the automatic quantification of RIF have been developed [27–34]. The majority of these analysis solutions use a single 2D image to count foci, most commonly generated from the MIP of an image stack, but more recently efforts have been made to score multiple focal planes independently. Increasing the axial resolution, i.e. through the nuclear depth, is becoming more important as we ask more mechanistic questions on the spatio-temporal aspects of RIF such as determining the composition of RIF at varying times after irradiation and also, assessing the relevance of RIF in the formation of chromosome exchanges. Thus our aim was to develop a rapid, high-throughput, high fidelity automated RIF analysis platform that required minimum user input for the quantification of the number, size/volume and relative 3D distribution of RIF within individual nuclei.
Material and methods
Cell culture, irradiation and immunofluorescence
Primary human bronchial epithelial (NHBE) cells (Lonza) were cultured in complete medium (Lonza BulletKit CC-3170) which consists of bronchial epithelial basal medium supplemented with bovine pituitary extract (0.2%), insulin (0.1%), hydrocortisone (0.1%), gentamicin sulfate and amphotericin-b (0.1%), retinoic acid (0.1%), transferrin (0.1%), epinephrine (0.1%) and human epithelial growth factor (0.1%). Cells were seeded at a density of 3.5 × 103 cells/cm2 and routinely sub-cultured at~ 80-90% confluence by trypsinising according to the suppliers guide. In brief, the cell sheet was rinsed with pre-warmed HEPES-BSS, incubated with Trypsin/EDTA solution at 37°C until > 90% of cells had rounded up and become detached before the addition of an excess of Trypsin Neutralisation Solution (TNS). For experimentation passage 3-7 NHBE cells were seeded onto sterilised glass microscope slides (Menzel) in quadriPERM® dishes (Sigma) and cultured for between 2 and 3 days until 90% confluent. Cells were then exposed to 60Cobalt γ-rays (~0.33 Gy/min) at 37°C and incubated for varying lengths of time before being washed three times with ice-cold PBS and fixed in 4% paraformaldehyde at room temperature (RT) for 5-10 min. For immunofluorescence, cells were washed three times for 3-5 min each in PBS, permeabilised (5% saponin (w/v) and 5% triton X-100 (v/v) in PBS) for 20 min at RT and incubated in blocking buffer (5% foetal bovine serum in PBS) for 1 hr at RT. Cells were then incubated with primary antibody (mouse monoclonal anti-human 53BP1 (BD Biosciences Clone 19) 1:200 in blocking buffer) for 1 hr at RT, washed three times in PBS for 5 min each with agitation then incubated with secondary antibody (goat anti mouse IgG conjugated with AlexaFluor 555 (Invitrogen)) again for 1 hr at RT. Cells were subsequently washed three times for 5 min each in PBS with agitation, mounted and counterstained with vectashield containing 4',6-diamidino-2-phenylindole (DAPI).
Image acquisition using widefield fluorescence microscopy
Images of fixed NHBE nuclei were acquired on an Axiovert 200 M microscope (Carl Zeiss) equipped with 100× NA 1.3 objective lens and Axiocam HR (14 bit) camera, and controlled with Axiovision software. Image dimensions of 1300 × 1030 at a pixel scale of 0.07 μm/pixel provided a field of view (FOV) of approximately 90 × 70 μm. The exposure time for each fluorescence channel was determined by examining multiple FOV for each slide to identify the highest intensity as to ensure images were not overexposed, whilst balancing the highest signal to noise ratio with the need to preserve the fluorescence signal across the whole slide. This was typically in the range of 100-600 ms for radiation-induced foci (RIF) and 30-80 ms for DAPI. Image stacks of 22-25 slices were acquired at 0.5 μm intervals from the central focal plane with FOV being selected randomly using the DAPI fluorescence channel ensuring every image stack contained at least one nucleus. Acquisition times varied but typically took in the region of 2 min/FOV. For subsequent automated analysis, each image stack was saved within its own folder as tagged image file format (TIFF) series.
Image processing for analysis of radiation-induced foci (RIF)
Manual analysis for the quantification of RIF was conducted on coded slides using Axiovision software (Carl Zeiss) to view the images. To categorise RIF into size integers an acetate sheet with two scaled circles corresponding to 0.5 μm and 1.0 μm was used. RIF smaller than the 0.5 μm circle were categorised as "small", RIF larger than this but smaller than the 1.0 μm circle were categorised as "medium" and any RIF larger than 1.0 μm in diameter were categorised as "large". RIF from multiple focal planes were counted and categorised according to size by scrolling up and down through the image stack, whereby the most in-focus area was used to determine the foci size (diameter). Where RIF had indiscriminate edges or were not sharp enough to discern a boundary and therefore a categorisation of size, they were classified as being "out-of-focus" (OOF).
Development of automated analysis system (AutoRIF)
The filtering stage comprises several steps. First, a modified Hessian operator is applied to remove high-frequency noise (Figure 1D) which gives the effect of smoothing the plots to allow accurate RIF dimensions to be extracted. This step does not affect the mask outline that is made from the DAPI channel but does remove the noise spikes. Next, the Laplacian transform operator converts the smoothed image such that intensity peaks are translated into troughs or concavities (Figure 1E). Although there are fluctuations in curvature of the plot these can easily be distinguished from the RIF. Any noise which escapes the smoothing phase would have been amplified by the application of the Laplacian operator. Therefore, to remove any artificially enhanced noise from the 'Laplacian' image, the Crimmins filter is used. This brightens pixels that are significantly darker than their neighbours, whilst it darkens those that are lighter than their neighbours (Figure 1F). The maximum entropy threshold, derived from the data of the whole image stack, is then used to identify the RIF. The resulting binary image then undergoes the morphological closing operation, using a kernel size of 3x3, which has the effect of 'filling the gaps' in the individual foci (Figure 1G).
This localisation algorithm uses the above to return the number of RIF in a slice, and for each RIF their centre coordinates, area, rotation angle and radii.
Application of automated analysis system
Primary human bronchial epithelial (NHBE) cells were irradiated with sham or 2 Gy 60Co-γ-rays (0.33 Gy/min) and fixed with 4% paraformaldehyde at varying times after exposure. The recruitment of 53BP1 proteins to DNA double strand breaks (DSB) were detected using mouse anti-human 53BP1 (Clone 19 BD Biosystems) antibodies and secondary goat anti-mouse IgG (Invitrogen) antibodies tagged with Alexafluor 555. To test and validate AutoRIF for the analysis of radiation-induced foci (RIF), digital images from a single focal plane (2D) and image stacks from multiple focal planes (3D) were acquired using a widefield microscope (100 nuclei per time point) and used for non-biased analysis using both manual and AutoRIF methods.
The quantification and measurement of RIF is a rapidly developing field. The benefits and limitations of the RIF assay have been reviewed  but there has yet to be a standardisation in the method for RIF analysis. Of the automated systems currently available, the majority utilise filtering processes to enhance contrast between signal (RIF) and noise (background). The 'top hat' filter is commonly used in conjunction with smoothing filters to improve contrast for RIF analysis [34, 37–40] which then allows a threshold to be applied to segment RIF from the background. The application of a fixed threshold alone is possible but there must be consistent contrast between RIF and background for this to be really effective . It is more practical to use a dynamic threshold that changes automatically with each image based on the frequency distribution of pixel intensities (histogram) . Alternatively, the identification of the centres of RIF from peaks of intensity (maxima) used in conjunction with water-shedding to identify the periphery of RIF  can be used. One issue with this is that it can lead to unnecessary segmentation of larger RIF and over-estimates from noisy background which although has been rectified with more complex algorithms , requires significantly more processing power and time.
The approach taken here was to identify nuclei using the dynamic Otsu threshold, which was then used as a mask to exclude non-nuclear signal from the RIF image. The Otsu method was applied to a MIP image produced by a projection of all stack images into a single one. Using such a composite image for nuclei segmentation reduces the chance of false-positives in Otsu methods when some of the images in each stack have only noise and no distinct nuclei. Image contrast of RIF was then enhanced using a Laplacian filter with smoothing processes, which reduces the risk of amplifying noise and generating artefacts. Finally, the maximum entropy threshold, based on the histogram of all the images in a stack (global), is applied to segment RIF from background in what has been described as pseudo 3D processing . The rationale for including all images from the stack for the threshold calculation instead of calculating the threshold for each image separately is based on the following: noise is accumulated into the histogram from all the images and given that it is of the same average value across, becomes more distinct as a class. As well as noise, the signal (RIF) curvature spikes also get accumulated in the global histogram making the separation between noise and signal bigger. This makes the maximum entropy threshold class separation more precise. If the histogram were based on every individual image, in images with no distinct RIF, the maximum entropy method would segment the noise into two classes: strong and weak noise, so the stronger noise would be mistaken for RIF. However when the common threshold value based on the global histogram is applied to the whole stack, the images that do not have distinct RIF would be subject to the same threshold as used in images containing RIF. This prevents the noise from being exaggerated in slices without RIF. There is also the functionality to combine histograms from entire folders of results to further increase the strictness of the threshold, which extends the range of intensities that can be analysed within an experiment. The two filtering methods discussed (top hat and Laplacian) are not mutually exclusive, in fact a combination of Laplacian and top hat filtering has been used to identify similar objects in images from chest X-rays .
At each stage of this development validity checks using binary images and comparisons of overall scores were made, these checks appear to be standard for systems involving the quantification of foci , but detailed validation in relation to size is not . Our results show there to be an excellent correlation between AutoRIF and manual analysis for the average number of RIF per nucleus, trends in induction and decline of RIF over time and, the frequency distributions for nuclei (Figures 7 and 8) for cells exposed to sham and 2 Gy radiation, while work is ongoing to assess these comparisons across a range of radiation doses. There is also good concordance between measurements of RIF size, but comparisons between manual and AutoRIF size proportions should be considered cautiously as the measurements are different. In the manual counting method, each RIF is only considered on a single focal plane and only the diameter is measured, for AutoRIF the volume is measured based on the theoretical 3D periphery of the RIF, (that is an elliptical approximation of the true boundary); the size integers applied to the AutoRIF score assume RIF to have an even and regular outline, which we know to be incorrect. A more accurate method for comparison between the methods would require more complex modelling of foci.
Some differences between manual and AutoRIF quantification were observed (Figure 7, 8, 9, 10) however these were identified to be due to outliers in the data caused by the very close bounding of two nuclei and poor fluorescence signal (data not shown). Thus, where AutoRIF is to be used in conjunction with images that were acquired using an automated image acquisition platform then it would be recommended to screen images e.g. as a thumbnail prior to the analysis processing. In addition to the above consideration, the application of any automated analysis system and subsequent interpretation of the data generated relies on the operators initial input of threshold parameters. This is no different for the application of AutoRIF whereby a robust regime for the experimental batch optimisation of filter settings, as described in Figures 5 and 6, is critical to ensure rigour and reproducibility. This filter optimisation process will be the same irrespective of dose or radiation quality meaning the software is flexible and applicable for a range of exposure types although it is likely that optimal parameters will differ for varying exposure types. For instance, exposure to high-LET radiation will (depending on the time after radiation) result in multiple, closely spaced RIF in those nuclei traversed meaning there is a potential for multiple small RIF to be identified as single larger RIF, if filter settings are not initially validated by comparing with the original raw images for that experiment. This batch optimisation process is by no means a caveat, indeed on-going developments of AutoRIF could incorporate a form of iterative or reinforcement machine learning, whereby the operator is presented with a sample image of what the algorithm believes is correct, the operator can then correct the image and the algorithm automatically adapts.
The main benefit of this analysis tool is the rapid turnaround of data. A single experiment with 10 slides and 100 nuclei per slide takes ~1-2 weeks to categorise all the RIF into size integers if manual analysis methods are used. Although reproducible this method is subject to bias and error. However, the same data can be obtained with more depth of information and consistency within 2 hrs (depending on hardware) with minimal operator input using AutoRIF. Specifically, consistent measurement of RIF size and signal intensity is generated for all RIF detected that is based on the number of pixels and is therefore continuous data rather than pre-defined categories. Further, this system generates x, y and z positional data for each RIF and provides distances between each RIF to all other RIF and also, to the nuclear centre, in each nucleus meaning that the spatial relationship of RIF can be studied at varying times after induction. Thus in addition to rapid quantification of RIF size (volume) and intensity, this automated system is a powerful tool for studying spatio-temporal relationships of RIF using a range of DDR markers.
- Lobrich M, Shibata A, Beucher A, Fisher A, Ensminger M, Goodarzi AA, Barton O, Jeggo PA: gamma H2AX foci analysis for monitoring DNA double-strand break repair Strengths, limitations and optimization. Cell Cycle. 2010, 9 (4): 662-669. 10.4161/cc.9.4.10764.View ArticlePubMedGoogle Scholar
- O'Neill P, Leatherbarrow EL, Harper JV, Cucinotta FA: Induction and quantification of gamma-H2AX foci following low and high LET-irradiation. International Journal of Radiation Biology. 2006, 82 (2): 111-118. 10.1080/09553000600599783.View ArticlePubMedGoogle Scholar
- Svetlova MP, Solovjeva LV, Tomilin NV: Mechanism of elimination of phosphorylated histone H2AX from chromatin after repair of DNA double-strand breaks. Mutation Research-Fundamental and Molecular Mechanisms of Mutagenesis. 2010, 685 (1-2): 54-60. 10.1016/j.mrfmmm.2009.08.001.View ArticlePubMedGoogle Scholar
- Okayasu R, Kato TA, Bedford JS: Comparison of the induction and disappearance of DNA double strand breaks and gamma-H2AX foci after irradiation of chromosomes in G1-phase or in condensed metaphase cells. Mutation Research-Fundamental and Molecular Mechanisms of Mutagenesis. 2008, 639 (1-2): 108-112. 10.1016/j.mrfmmm.2007.11.006.View ArticlePubMedGoogle Scholar
- Karlsson KH, Stenerlow B: Focus formation of DNA repair proteins in normal and repair-deficient cells irradiated with high-LET ions. Radiation Research. 2004, 161 (5): 517-527. 10.1667/RR3171.View ArticlePubMedGoogle Scholar
- Belyaev IY, Markova E, Schultz N: Kinetics and dose-response of residual 53BP1/gamma-H2AX foci: Co-localization, relationship with DSB repair and clonogenic survival. International Journal of Radiation Biology. 2007, 83 (5): 319-329. 10.1080/09553000601170469.View ArticlePubMedGoogle Scholar
- Schultz LB, Chehab NH, Malikzay A, Halazonetis TD: p53 Binding Protein 1 (53BP1) Is an Early Participant in the Cellular Response to DNA Double-Strand Breaks. The Journal of Cell Biology. 2000, 151 (7): 1381-1390. 10.1083/jcb.151.7.1381.PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson L, Henderson C, Adachi Y: Phosphorylation and rapid relocalization of 53BP1 to nuclear foci upon DNA damage. Mol Cell Biol. 2001, 21 (5): 1719-1729. 10.1128/MCB.21.5.1719-1729.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Rappold I, Iwabuchi K, Date T, Chen J: Tumor Suppressor p53 Binding Protein 1 (53BP1) Is Involved in DNA Damage-signaling Pathways. The Journal of Cell Biology. 2001, 153 (3): 613-620. 10.1083/jcb.153.3.613.PubMed CentralView ArticlePubMedGoogle Scholar
- Ward IM, Minn K, Jorda KG, Chen J: Accumulation of checkpoint protein 53BP1 at DNA breaks involves its binding to phosphorylated histone H2AX. The Journal of biological chemistry. 2003, 278 (22): 19579-19582. 10.1074/jbc.C300117200.View ArticlePubMedGoogle Scholar
- Bekker-Jensen S, Lukas C, Melander F, Bartek J, Lukas J: Dynamic assembly and sustained retention of 53BP1 at the sites of DNA damage are controlled by Mdc1/NFBD1. The Journal of cell biology. 2005, 170 (2): 201-211. 10.1083/jcb.200503043.PubMed CentralView ArticlePubMedGoogle Scholar
- Scully R, Xie AY, Hartlerode A, Stucki M, Odate S, Puget N, Kwok A, Nagaraju G, Yan C, Alt FW, Chen J, Jackson SP: Distinct roles of chromatin-associated proteins MDC1 and 53BP1 in mammalian double-strand break repair. Molecular Cell. 2007, 28 (6): 1045-1057. 10.1016/j.molcel.2007.12.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Zgheib O, Pataky K, Brugger J, Halazonetis TD: An oligomerized 53BP1 tudor domain suffices for recognition of DNA double-strand breaks. Molecular and cellular biology. 2009, 29 (4): 1050-1058. 10.1128/MCB.01011-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Huyen Y, Zgheib O, DiTullio RA, Gorgoulis VG, Zacharatos P, Petty TJ, Sheston EA, Mellert HS, Stavridi ES, Halazonetis TD: Methylated lysine 79 of histone H3 targets 53BP1 to DNA double-strand breaks. Nature. 2004, 432: 406-411. 10.1038/nature03114.View ArticlePubMedGoogle Scholar
- Costes SV, Ponomarev A, Chen JL, Nguyen D, Cucinotta FA, Barcellos-Hoff MH: Image-based modeling reveals dynamic redistribution of DNA damage into nuclear sub-domains. PLoS Comput Biol. 2007, 3 (8): e155-10.1371/journal.pcbi.0030155.PubMed CentralView ArticlePubMedGoogle Scholar
- Jakob B, Splinter J, Durante M, Taucher-Scholz G: Live cell microscopy analysis of radiation-induced DNA double-strand break motion. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (9): 3172-3177. 10.1073/pnas.0810987106.PubMed CentralView ArticlePubMedGoogle Scholar
- Banath JP, Klokov D, MacPhail SH, Banuelos CA, Olive PL: Residual gammaH2AX foci as an indication of lethal DNA lesions. BMC cancer. 2010, 10: 4-10.1186/1471-2407-10-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Rothkamm K, Balroop S, Shekhdar J, Fernie P, Goh V: Leukocyte DNA damage after multi-detector row CT: A quantitative biomarker of low-level radiation exposure. Radiology. 2007, 242 (1): 244-251. 10.1148/radiol.2421060171.View ArticlePubMedGoogle Scholar
- Beels L, Bacher K, Smeets P, Verstraete K, Vral A, Thierens H: Dose-length product of scanners correlates with DNA damage in patients undergoing contrast CT. European Journal of Radiology. 2011,Google Scholar
- Lassmann M, Hanscheid H, Gassen D, Biko J, Meineke V, Reiners C, Scherthan H: In vivo formation of gamma-H2AX and 53BP1 DNA repair foci in blood cells after radioiodine therapy of differentiated thyroid cancer. Journal of nuclear medicine: official publication, Society of Nuclear Medicine. 2010, 51 (8): 1318-1325.View ArticleGoogle Scholar
- Roch-Lefevre S, Mandina T, Voisin P, Gaetan G, Mesa JEG, Valente M, Bonnesoeur P, Garcia O, Voisin P, Roy L: Quantification of gamma-H2AX Foci in Human Lymphocytes: A Method for Biological Dosimetry after Ionizing Radiation Exposure. Radiation Research. 2010, 174 (2): 185-194. 10.1667/RR1775.1.View ArticlePubMedGoogle Scholar
- Bourton EC, Plowman PN, Smith D, Arlett CF, Parris CN: Prolonged expression of the γ-H2AX DNA repair biomarker correlates with excess acute and chronic toxicity from radiotherapy treatment. International Journal of Cancer. 2011, 129 (12): 2928-2934. 10.1002/ijc.25953.View ArticleGoogle Scholar
- Paull TT, Rogakou EP, Yamazaki V, Kirchgessner CU, Gellert M, Bonner WM: A critical role for histone H2AX in recruitment of repair factors to nuclear foci after DNA damage. Current Biology. 2000, 10 (15): 886-895. 10.1016/S0960-9822(00)00610-2.View ArticlePubMedGoogle Scholar
- Riballo E, Kuhne M, Rief N, Doherty A, Smith GC, Recio MJ, Reis C, Dahm K, Fricke A, Krempler A, Parker AR, Jackson SP, Gennery A, Jeggo PA, Lobrich M: A pathway of double-strand break rejoining dependent upon ATM, Artemis, and proteins locating to gamma-H2AX foci. Molecular cell. 2004, 16 (5): 715-724. 10.1016/j.molcel.2004.10.029.View ArticlePubMedGoogle Scholar
- Desai N, Davis E, O'Neill P, Durante M, Cucinotta FA, Wu H: Immunofluorescence detection of clustered gamma-H2AX foci induced by HZE-particle radiation. Radiat Res. 2005, 164 (4 Pt 2): 518-522.View ArticlePubMedGoogle Scholar
- Hamada N, Schettino G, Kashino G, Vaid M, Suzuki K, Kodama S, Vojnovic B, Folkard M, Watanabe M, Michael BD, Prise KM: Histone H2AX phosphorylation in normal human cells irradiated with focused ultrasoft × rays: evidence for chromatin movement during repair. Radiat Res. 2006, 166 (1 Pt 1): 31-38.View ArticlePubMedGoogle Scholar
- Su Y, Meador JA, Geard CR, Balajee AS: Analysis of ionizing radiation-induced DNA damage and repair in three-dimensional human skin model system. Experimental dermatology. 2010, 19 (8): e16-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Bocker W, Iliakis G: Computational Methods for Analysis of Foci: Validation for Radiation-Induced g-H2AX Foci in Human Cells. Radiation Research. 2006, 165 (113-124):Google Scholar
- Barber PR, Locke RJ, Pierce GP, Rothkamm K, Vojnovic B: Gamma-H2AX Foci Counting: Image processing and control software for high-content screening. Imaging, Manipulation, and Analysis of Biomolecules, Cells, and Tissues V. 2007, 6441: M4411-M4411.Google Scholar
- Gerashchenko BI, Dynlacht JR: A tool for enhancement and scoring of DNA repair foci. Cytometry. Part A: the journal of the International Society for Analytical Cytology. 2009, 75 (3): 245-252.View ArticleGoogle Scholar
- Costes SV, Ponomarev A, Chen JL, Nguyen D, Cucinotta FA, Barcellos-Hoff MH: Image-based modeling reveals dynamic redistribution of DNA damage into nuclear sub-domains. PLoS computational biology. 2007, 3 (8): e155-10.1371/journal.pcbi.0030155.PubMed CentralView ArticlePubMedGoogle Scholar
- Ishikawa A, Yamauchi M, Suzuki K, Yamashita S: Image-based quantitative determination of DNA damage signal reveals a threshold for G2 checkpoint activation in response to ionizing radiation. Genome integrity. 2010, 1 (1): 10-10.1186/2041-9414-1-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Cai Z, Vallis KA, Reilly RM: Computational analysis of the number, area and density of gamma-H2AX foci in breast cancer cells exposed to (111)In-DTPA-hEGF or gamma-rays using Image-J software. International journal of radiation biology. 2009, 85 (3): 262-271. 10.1080/09553000902748757.View ArticlePubMedGoogle Scholar
- Ivashkevich AN, Martin OA, Smith AJ, Redon CE, Bonner WM, Martin RF, Lobachevsky PN: [gamma]H2AX foci as a measure of DNA damage: A computational approach to automatic analysis. Mutation Research/Fundamental and Molecular Mechanisms of Mutagenesis. 2011, 711 (1-2): 49-60. 10.1016/j.mrfmmm.2010.12.015.View ArticlePubMedGoogle Scholar
- Otsu N: A Threshold Selection Method from Gray-Level Histograms. Systems, Man and Cybernetics, IEEE Transactions on. 1979, 9 (1): 62-66.View ArticleGoogle Scholar
- Löbrich M, Shibata A, Beucher A, Fisher A, Ensminger M, Goodarzi AA, Barton O, Jeggo PA: γH2AX foci analysis for monitoring DNA double-strand break repair: Strengths, limitations and optimization. Cell Cycle. 2010, 9 (4): 662-669. 10.4161/cc.9.4.10764.View ArticlePubMedGoogle Scholar
- Valente M, Voisin P, Laloi P, Roy L, Roch-Lefèvre S: Automated gamma-H2AX focus scoring method for human lymphocytes after ionizing radiation exposure. Radiation Measurements. 2011,Google Scholar
- Qvarnstrom OF, Simonsson M, Johansson KA, Nyman J, Turesson I: DNA double strand break quantification in skin biopsies. Radiother Oncol. 2004, 72 (3): 311-317. 10.1016/j.radonc.2004.07.009.View ArticlePubMedGoogle Scholar
- Jucha A, Wegierek-Ciuk A, Koza Z, Lisowska H, Wojcik A, Wojewodzka M, Lankoff A: FociCounter: A freely available PC programme for quantitative and qualitative analysis of gamma-H2AX foci. Mutation research. 2010, 696 (1): 16-20.View ArticlePubMedGoogle Scholar
- Roch-Lefevre S, Mandina T, Voisin P, Gaetan G, Mesa JE, Valente M, Bonnesoeur P, Garcia O, Roy L: Quantification of gamma-H2AX foci in human lymphocytes: a method for biological dosimetry after ionizing radiation exposure. Radiation research. 2010, 174 (2): 185-194. 10.1667/RR1775.1.View ArticlePubMedGoogle Scholar
- Hou Y, Lavaf A, Huang D, Peters S, Huq R, Friedrich V, Rosenstein BS, Kao J: Development of an Automated g-H2AX Immunocytochemistry Assay. Radiation Research. 2009, 171: 360-367. 10.1667/RR1349.1.View ArticlePubMedGoogle Scholar
- Barber PL, Pierce G, Rothkamm K, Vojnovic B: Gamma-H2AX Foci Counting: Image processing and control software for high-content screening. 2007Google Scholar
- Shi Z, Zhao M, He L, Wang Y: A Combinational Filtering Method for Enhancing Suspicious Structures in Chest X-rays. Journal of Information and Computational Science. 2011, 8 (7): 997-1005.Google 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.