Mathematical modelling of the automated FADU assay for the quantification of DNA strand breaks and their repair in human peripheral mononuclear blood cells
© Junk et al.; licensee BioMed Central Ltd. 2014
Received: 8 May 2014
Accepted: 19 August 2014
Published: 9 September 2014
Cells continuously undergo DNA damage from exogenous agents like irradiation or genotoxic chemicals or from endogenous radicals produced by normal cellular metabolic activities. DNA strand breaks are one of the most common genotoxic lesions and they can also arise as intermediates of DNA repair activity. Unrepaired DNA damage can lead to genomic instability, which can massively compromise the health status of organisms. Therefore it is important to measure and quantify DNA damage and its repair.
We have previously published an automated method for measuring DNA strand breaks based on fluorimetric detection of alkaline DNA unwinding , and here we present a mathematical model of the FADU assay, which enables to an analytic expression for the relation between measured fluorescence and the number of strand breaks.
Assessment of the formation and also the repair of DNA strand breaks is a crucial functional parameter to investigate genotoxicity in living cells. A reliable and convenient method to quantify DNA strand breakage is therefore of significant importance for a wide variety of scientific fields, e.g. toxicology, pharmacology, epidemiology and medical sciences.
KeywordsFADU DNA repair DNA strand breaks Mathematical model
The fact that the DNA denatures (unwinds) under alkaline conditions has been used to establish FADU assay for measuring DNA strand breaks . This method is based on the limited denaturation of the DNA under precisely controlled conditions of pH, temperature and time. Under these conditions only chromosome ends plus open sites in the DNA serve as starting points for the unwinding, which proceeds bidirectionally from internal damage sites. After stopping the unwinding process at a specified time, a dye is added that interacts with the remaining double stranded DNA and during contact with DNA emits a fluorescence signal whose intensity is inversely related to the number of DNA strand breaks originally present: the less fluorescence signal, the more DNA strand breaks.
When X-radiation is applied to living cells, the number of strand breaks is known to increase linearly with the dose applied ,. However, the fluorescence signal intensity depends on the dose in a nonlinear fashion and displays saturation. Previously, the percentage of double stranded DNA, remaining after the alkaline unwinding, relative to undenatured control or more precisely the decadic logarithm of the intensity ratio with and without unwinding has been used for that purpose, but this does not fully result in a linear function with radiation dose. In this article, we present a mathematical model that describes the effect of DNA unwinding on the resulting strength of the fluorescence signal and thereby captures the non-linear relationship more precisely.
The main assumptions of the model are large cell numbers, homogeneously distributed damage sites, and small variation of cell parameters like initial damage, susceptibility to applied damage, and potential unwinding distance from the site of damage. The resulting formula relates the relative fluorescence intensity to the applied X-ray dose and fits well with experiments. To obtain this result, it suffices to use elementary arguments from probability theory. In particular, a detailed biochemical model of the unwinding process as in  is not required.
In other words, the relative fluorescence intensity I depends on the average relative length of double stranded DNA and the relative background fluorescence B under total unwinding. Later, the fact that N is large allows us to apply the mathematical law of large numbers to safely replace the average by an expected value.
It remains to model the relative length L i of double stranded DNA after the unwinding process and its relation to the damage D i . Since we assume that the process happens identically and independently in each cell, we consider a generic situation and drop the cell index i in the following discussion.
Our basic assumption is that the break points x k are uniformly distributed along the DNA strands and that the DNA unwinds a distance δ to the left and to the right of x k unless there is interference with the unwinding process going on at a neighboring break.
where Δ k = z k + 1 -z k is the distance between consecutive break points. For m we use the splitting m + 1 = D 0 + D x into the number of breaks D x = βd due to X-ray damage and the number D 0 accounting for breaks at zero damage, resulting from normal cell metabolism. Note that in our model, D 0 is at least the number of chromosomes because we have added their terminal points as artificial breaks. Due to the stochastic nature of the damaging process, the cell metabolism and the unwinding, the number of strand breaks per dose β, the zero dose value D 0 and the unwinding distance δ should be considered variables whose value may vary randomly from cell to cell.
To estimate the resulting distribution of the relative fluorescence signal ƒ with some manageable expression it seems worthwhile to simplify the model further, even if it requires additional assumptions. For example, the situation simplifies a lot if we assume that all z 1, , z m are independently and uniformly distributed in the interval [0,1], which means that the particular role of the original terminal chromosome points is dropped. Since in a normal diploid human cell, the number of chromosomes is only 46, this may be acceptable as long as m is comparably large.
In order to compute the unconditional expectation E(L), reasonable assumptions on the probability distributions of D 0, β and δ are required. Finally, E(L) depends in a complicated way on the dose d, the mean values , the corresponding variances and maybe other parameters of the distributions.
which is a reasonable assumption when D 0, β and δ have negligible variance among all cells in the assay. Then, the effective values and δ eff are close to the mean values .
In order to test our model, we use experimental results from the experiment specified in Material and Methods. Here, three replicates of the relative intensities at seven dose values are available for eleven donors. We denote them P jkl where the index j = 0,..,6 labels the X-ray dose, k = 1,..., 11 represents the donor and l = 1,2,3 indicates the independent repetitions of the experiment.
where is the residual in case of the optimal parameters. The optimization is carried out with the Levenberg-Marquardt algorithm of the Matlab software package with a resulting value of R 2 = 0.99.
Confidence intervals for the fit parameters are computed with the bootstrap method  exploiting the information contained in the threefold repetition of each experiment. For each index combination (j, k, l), N = 10000 bootstrap samples are generated by drawing intensities randomly from the values P jk1, P jk2, P jk3 with replacement. Then the computation of the optimal parameters is repeated for each n = 1,..., N and the resulting distribution is used to construct 95%-confidence intervals.
Fit parameters for the relative fluorescence intensity
0.79 ± 0.01
0.18 ± 0.02
0.74 ± 0.06
0.24 ± 0.04
0.64 ± 0.04
0.21 ± 0.03
0.65 ± 0.02
0.21 ± 0.02
0.51 ± 0.04
0.20 ± 0.05
0.65 ± 0.05
0.20 ± 0.04
0.60 ± 0.03
0.30 ± 0.04
0.63 ± 0.02
0.25 ± 0.02
0.66 ± 0.04
0.31 ± 0.05
0.74 ± 0.03
0.28 ± 0.02
0.66 ± 0.08
0.18 ± 0.06
DNA strand breaks are one of the most common genotoxic lesions. One of the methods used for measuring DNA strand breaks and their repair is the FADU assay. Under the selected alkaline conditions a large number of base pairs will be unwound to the right and left starting from each DNA "open side". Due to this amplification it is possible to detect very low number of DNA strand breaks. The main advantage of the FADU assay is that it can be performed in a fully automated fashion. However there are few methods that can be performed in a semi-automated version such as the COMET and γH2AX assays. Even though these techniques provide useful tools for measuring DNA strand breaks in a high-throughput fashion, they have some disadvantages. A comprehensive comparison between these assays has been published before . Briefly, γH2AX assay is considered a very sensitive method that specifically detects double DNA strand breaks however in some cases the presence of γH2AX foci in the absence of DNA damage has been demonstrated. The comet assay protocol includes many different steps making difficult the development of a fully automated version. Furthermore these also affect both intra-assay variability and inter-assay reproducibility.
We have previously published a modified and automated version of the FADU assay ,,. Compared with the original FADU assay , the automated FADU assay shows comparable sensitivity, yet increased robustness and throughput, and decreased operator time. However a suitable mathematical model, which allows translating the fluorescence intensity in DNA damage increases accuracy. In order to test whether the model is capable of describing the experimentally observed dependence between applied damage and resulting fluorescence intensity, we have measured the DNA strand breaks in peripheral mononuclear blood cells of 11 individuals. It turns out that the model is able to fit the dose-effect relation very accurately.
As shown in Table1, the resulting model parameters λ and P 0, which describe the cell's susceptibility to DNA damage and the level of endogenous DNA strand breaks respectively, show a certain variation among the individuals which is stronger for the P 0 value. An individual endogenous level of DNA strand breaks is expected; in general the importance of endogenous DNA damage has been reported before . For example, it is known that DNA double-strand breaks (DSBs) accumulate in senescing human cell cultures and in ageing mice . This is in accordance with our recent data showing that senescing human T lymphocytes have acquired more DNA strand breaks in vivo and an impaired DNA repair capacity in vitro. Additionally, endogenous DNA strand breaks accumulation in human cells has been associated with nutrition , cancer , diabetes ,, Down's syndrome , rheumatoid arthritis  and psychological stress disorders , among others. In other words, differences in the level of endogenous DNA strand breaks might reflect aging processes or health status of individuals.
The susceptibility to DNA damage λ also shows variation between individuals. This finding is not surprising and has been described before. Odarigi and colleagues reported an inter-individual variance in chromosomal damage after x-irradiation . Since radiation induces free oxygen radical formation, an individual damage response to x-ray might result from the different capacity of cells to scavenge these free radicals.
The investigation of DNA strand breaks and their repair is of great interest to toxicology, pharmacology, epidemiology and medical sciences at large. In general impaired DNA repair can lead to embryonic lethality, impaired growth, accelerated aging, shortened life span, and increased incidence of a variety of diseases, including a pronounced manifestation of cancer -. Therefore robust methods for quantification of DNA damage and its repair are necessary.
Finally, in our model, induced DNA strand breaks are assumed to be uniformly distributed. However chromatin state could have an impact on the distribution of DNA strand breaks. Chromatin state differs depending on cell type and cell status (cell cycle phase and/or transcription) therefore it is important to take these aspects into consideration. A mathematical model of DNA unwinding considering the chromatin state would be an interesting task for future research, especially when comparing different cell types or cell status. A more refined model could be based on a probability distribution, which describes the susceptibility to breaking under X-ray exposure along the DNA strand.
We have applied the above mathematical model for the quantification of DNA strand breaks in irradiated human PBMCs. This model can be also applied to other cell types and other DNA damaging agents as long as the damaging agent leads to uniformly distributed DNA strand breaks. Furthermore measuring the residual DNA strand breaks in damaged cells after several periods of time allow us to calculate the rate of repaired DNA strand breaks. Last but not least the amount of endogenous DNA strand breaks provides important information about the individual "cell-damage-status" and fluorescence signal after induced damage serves as biomarker of cell sensitivity to DNA strand breaks inducers, therefore an accurate model for DNA breakage quantification could be used for diagnostic and therapy purposes towards a personalized medicine, e.g. in cancer therapy.
Selection of volunteers was carried out in accordance with the Declaration of Helsinki and ethical approval was obtained from the Ethics Committee of University of Constance. A signed Informed Consent was obtained from each subject. Blood from healthy volunteers aged 23-55 was collected in sodium citrat tubes (Sarstedt, Nümbrecht, Germany). The automated FADU assay has previously been described in detail . Briefly, peripheral blood mononuclear cells (PBMC) were isolated by Biocoll (Biochrom, Berlin, Germany) density gradient centrifugation. Cells were counted using a cell counting device (Casy® counter), pelleted (10 min, 200 g), and resuspended in isotonic buffer (0.25 M meso-inositol; 10 mM sodium phosphate, pH 7.4; 1 mM magnesium chloride) or RPMI-1640 medium (Invitrogen) containing 100 U/ml penicillin (Invitrogen), 100 Mg/ml streptomycin (Invitrogen) and 10% fetal calf serum (FCS) at 4 × 106 cells per ml. In order to induce DNA strand breaks several aliquots of 100 l cell suspension were irradiated on ice at a dose rate of 1.9 Gy/min using an X-ray generator (CHF Müller, Hamburg, Germany, 70 keV, 1 mM Al-filter). Qualitative behaviour has been modelled, but it still depends on few quantitative parameters, which where determined based on the results obtained from unwinding experiments after various (0.3, 1, 1.9, 3.8, 7.5 and 14.9 Gy) irradiation doses (Figure 1).
The detection of DNA strand breaks is based on progressive DNA unwinding (denaturation) under highly controlled conditions of alkaline pH, time and temperature. The starting points for the unwinding process are DNA "open sites" like replication forks or chromosome ends, but also DNA strand breaks. For monitoring DNA damage a commercially available fluorescent dye (SYBR Green®) is used as marker for double stranded DNA, and a decrease in the fluorescence intensity of SYBR Green® indicates an increase of DNA unwinding and consequently a greater number of DNA strand breaks.
T and P0 are controls to be run in parallel with the experimentally treated cells. In P0 samples alkaline unwinding is allowed and represents the DNA strand breaks under physiological conditions (i.e. without exogenous DNA damage). In T samples the neutralisation follows the lysis and therefore unwinding cannot take place; T-samples provide a measure of total DNA content and yield a fluorescence signal defined as 100%. Px samples (P1,P2,P3,P4....Px) are the different extent of damage to be measured.
The work presented here was carried out in collaboration between all authors. MMV defined the research theme and designed methods and experiments. JS carried out the laboratory experiments and analyzed the data. TS co-worked on data interpretation and discussed analyses. MJ developed the mathematical model and carried out the parameter fitting. AB co-designed experiments, discussed analyses, interpretation, and presentation. MMV and MJ wrote the paper. All authors have contributed to, seen and approved the manuscript.
This work was supported by the European Union FP7 Health research, grant number HEALTH-F4-2008-200880 [MARK-AGE]. The authors would like to acknowledge Ms. Gudrun von Scheven for organizing the blood collection.
- Moreno-Villanueva M, Pfeiffer R, Sindlinger T, Leake A, Müller M, Kirkwood T, Bürkle A: A modified and automated version of the `Fluorimetric Detection of Alkaline DNA Unwinding' method to quantify formation and repair of DNA strand breaks. BMC Biotechnol. 2009, 9: 39-10.1186/1472-6750-9-39.View ArticleGoogle Scholar
- Birnboim HC, Jevcak JJ: Fluorometric method for rapid detection of DNA strand breaks in human white blood cells produced by low doses of radiation. Cancer Res. 1981, 41: 1889-1892.Google Scholar
- Wood DK, Weingeist DM, Bhatia SN, Engelward BP: Single cell trapping and DNA damage analysis using microwell arrays. PNAS. 2010, 1;107 (22): 10008-10013. 10.1073/pnas.1004056107.View ArticleGoogle Scholar
- Olive PL, Banáth JP: The comet assay: a method to measure DNA damage in individual cells. Nat Protoc. 2006, 1 (1): 23-29. 10.1038/nprot.2006.5.View ArticleGoogle Scholar
- Gaver DP, Jacobs PA, Carpenter RL, Burkhart JG: A mathematical model for intracellular effects of toxins on DNA adduction and repair. Bull Math Biol. 1997, 59 (1): 89-106. 10.1007/BF02459472.MATHView ArticleGoogle Scholar
- Sokhansanj BA, Rodrigue GR, Fitch JP, Wilson DM: A quantitative model of human DNA base excision repair. I. Mechanistic insights. Nucleic Acids Res. 2002, 30 (8): 1817-1825. 10.1093/nar/30.8.1817.View ArticleGoogle Scholar
- Kesseler KJ, Kaufmann WK, Reardon JT, Elston TC, Sancar A: A mathematical model for human nucleotide excision repair: damage recognition by random order assembly and kinetic proofreading. J Theor Biol. 2007, 249 (2): 361-375. 10.1016/j.jtbi.2007.07.025.MathSciNetView ArticleGoogle Scholar
- Crooke PS, Parl FF: A mathematical model for DNA damage and repair. J Nucleic Acids. 2010, 2010: 352603-10.4061/2010/352603.View ArticleGoogle Scholar
- Devroye L: Non-Uniform Random Variate Generation. 1987, Springer Verlag, New YorkGoogle Scholar
- Efron B, Tibshirani RJ: An Introduction to the Bootstrap. Monographs on Statistics and Applied Probability 57. 1994, Chapman&Hall/CRC,Google Scholar
- Moreno-Villanueva M, Bürkle A: High-throughput Assays to Quantify the Formation of DNA Strand Breaks. High Throughput Screening Methods in Toxicity Testing. Edited by: Steinberg P. 2013, John Wiley & Sons, Inc, Hoboken, NJ, USA, 285-294. 10.1002/9781118538203.ch14.View ArticleGoogle Scholar
- Moreno-Villanueva M, Eltze T, Dressler D, Bernhardt J, Hirsch C, Wick P, von Scheven G, Lex K, Bürkle A: The automated FADU-assay, a potential high-throughput in vitro method for early screening of DNA breakage. ALTEX. 2011, 28 (4): 295-303. 10.14573/altex.2011.4.295.View ArticleGoogle Scholar
- De Bont R, van Larebeke N: Endogenous DNA damage in humans: a review of quantitative data. Mutagenesis. 2004, 19 (3): 169-185. 10.1093/mutage/geh025.View ArticleGoogle Scholar
- Sedelnikova OA, Horikawa I, Zimonjic DB, Popescu NC, Bonner WM, Barrett JC: Senescing human cells and ageing mice accumulate DNA lesions with unrepairable double-strand breaks. Nat Cell Biol. 2004, 6 (2): 168-170. 10.1038/ncb1095.View ArticleGoogle Scholar
- Brunner S, Herndler-Brandstetter D, Arnold CR, Wiegers GJ, Villunger A, Hackl M, Grillari J, Moreno-Villanueva M, Bürkle A, Grubeck-Loebenstein B: Upregulation of miR-24 is associated with a decreased DNA damage response upon etoposide treatment in highly differentiated CD8(+) T cells sensitizing them to apoptotic cell death. Aging Cell. 2012, 11 (4): 579-587. 10.1111/j.1474-9726.2012.00819.x.View ArticleGoogle Scholar
- Krajcovicová-Kudlácková M, Dusinská M: Oxidative DNA damage in relation to nutrition. Neoplasma. 2004, 51 (1): 30-33.Google Scholar
- Vilenchik MM, Knudson AG: Endogenous DNA double-strand breaks: production, fidelity of repair, and induction of cancer. Proc Natl Acad Sci. 2003, 100 (22): 12871-12876. 10.1073/pnas.2135498100.View ArticleGoogle Scholar
- Lorenzi M, Montisano DF, Toledo S, Wong HC: Increased single strand breaks in DNA of lymphocytes from diabetic subjects. J Clin Invest. 1987, 79 (2): 653-656. 10.1172/JCI112863.View ArticleGoogle Scholar
- Blasiak J, Arabski M, Krupa R, Wozniak K, Zadrozny M, Kasznicki J, Zurawska M, Drzewoski J: DNA damage and repair in type 2 diabetes mellitus. Mutat Res. 2004, 554 (1-2): 297-304. 10.1016/j.mrfmmm.2004.05.011.View ArticleGoogle Scholar
- Zana M, Szécsény A, Czibula A, Bjelik A, Juhász A, Rimanóczy A, Szabó K, Vetró A, Szucs P, Várkonyi A, Pákáski M, Boda K, Raskó I, Janka Z, Kálmán J: Age-dependent oxidative stress-induced DNA damage in down's lymphocytes. Biochem Biophys Res Commun. 2006, 345 (2): 726-733. 10.1016/j.bbrc.2006.04.167. Epub 2006 May 5View ArticleGoogle Scholar
- Bhusate LL, Herbert KE, Scott DL, Perrett D: Increased DNA strand breaks in mononuclear cells from patients with rheumatoid arthritis. Ann Rheum Dis. 1992, 51 (1): 8-12. 10.1136/ard.51.1.8.View ArticleGoogle Scholar
- Morath J, Moreno-Villanueva M, Gilava H, Kolassa S, Ruf-Leuschner M, Schauer M, Elbert T, Bürkle A, Kolassa I-T: Effects of psychotherapy on DNA strand break accumulation originating from traumatic stress. Psychother Psychosom. 2014, 83 (5): 289-297. 10.1159/000362739.View ArticleGoogle Scholar
- Odagiri Y, Uchida H, Shibazaki S: Interindividual variation in cytogenetic response to X-ray and colchicine measured with the cytokinesis-block micronucleus assay. Mutat Res. 1997, 381 (1): 1-13. 10.1016/S0027-5107(97)00138-3.View ArticleGoogle Scholar
- Hakem R: DNA-damage repair; the good, the bad, and the ugly. EMBO J. 2008, 27: 589-605. 10.1038/emboj.2008.15.View ArticleGoogle Scholar
- Diderich K, Alanazi M, Hoeijmakers JH: Premature aging and cancer in nucleotide excision repair-disorders. DNA Repair. 2011, 10 (7): 772-780. 10.1016/j.dnarep.2011.04.025.View ArticleGoogle Scholar
- Keimling M, Deniz M, Varga D, Stahl A, Schrezenmeier H, Kreienberg R, Hoffmann I, König J, WiesMüller L: The power of DNA double-strand break (DSB) repair testing to predict breast cancer susceptibility. FASEB J. 2012, 26 (5): 2094-2104. 10.1096/fj.11-200790.View ArticleGoogle Scholar
- Varga D, Deniz M, Schwentner L, WiesMüller L: Ovarian cancer: in search of better marker systems based on DNA repair defects. Int J Mol Sci. 2013, 14 (1): 640-673. 10.3390/ijms14010640.View ArticleGoogle Scholar
- Toulany M, Rodemann HP: Membrane receptor signaling and control of DNA repair after exposure to ionizing radiation. Nuklearmedizin. 2010, 2010 (49 Suppl 1): S26-S30.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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.