Computerized Determination of Growth Kinetic Curves and

ANALYTICALBIOCHEMISTRY
187,262-267
(1990)
Computerized Determination of Growth Kinetic Curves
and Doubling Times from Cells in Microculture
Herta
Reile,
Institut
fiir
Received
July
Herbert
Pharmazie,
BirnbGck,
Giinther
Sonderforschungsbereich
Bernhardt,
234, Universittit
Sprul3, and Helmut
Regensburg,
D-8400 Regensburg,
Schgnenberger
Federal Republic
of Germany
l&l989
In this paper we describe the microcomputer-aided
determination
of cell proliferation
kinetics
and doubling
times utilizing
a crystal violet assay and a 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazolium
bromide
assay in microtitration
plates. The analysis
of spectrophotometric
data provides
the doubling
times at any
time of incubation.
Plots of doubling
time versus time of
incubation
give reproducible
information
on the exact
duration
of the logarithmic
growth
phase. This method
is applicable
to anchorage-dependent
as well as anchorage-independent
cells when calorimetric
or fluorometric data are accessible.
o isso
Academic
press, h.
A great problem in the testing of new anticancer drugs
by the cell culture technique is the marked genetic instability of tumor cells (1,2) entailing an extreme clonal
variation (3). This becomes evident with growing numbers of passagesfrom a permanent change of parameters
of cell proliferation like the duration of the lag phase and
of the log period, the doubling time, and the saturation
density (plateau level). For example, Reddel et al. (4) observed such a perpetual decrease of doubling time in subcultures of the human breast cancer cell line T-47-D.
Since drugs can show quite different inhibitory effects
dependent of the growth parameters of a given tumor
cell line, the knowledge of the growth parameters in each
test series is of importance for the experimental reproducibility. Although it is generally postulated that the
exposure of the tumor cells to the drug should occur in
the phase of exponential growth, this demand often does
not meet experimental verification.
Conventionally, mean doubling times are determined
graphically from the mid log phase or calculated from
initial and final cell numbers (5) counted with either a
hemocytometer or an electronic particle counter which
requires removing the cells from their substratum. Either method could lead to severe misinterpretation due
262
Thilo
to practical difficulties in distinguishing the lag phase
from the phase of exponential growth.
In this publication we describe a microcomputer technique which allows the registration of growth curves of
cells in monolayer cultures by large scale spectrophotometric measurement after crystal violet staining’ (6) or
MTT incorporation (7-9) using 96-well microtitration
plates. This method enables the overall growth curve especially the lag phase and the log period to be determined in a more exact way. As a consequence it is also
possible to precisely assessthe doubling time (TAU) at
any time of the experiment. In our opinion a plot of TAU
versus time (t) provides the maximum information
available from cell proliferation experiments.
METHODS
Chemicals.
Reagents (A-grade purity) were obtained
from Merck. Crystal violet and MTT were purchased
from Serva; dimethyl sulfoxide (DMSO) spectrophotometric grade was obtained from Aldrich. Millipore-filtered water was used throughout.
Cell
lines
and
culture
conditions.
MDA-MB-231
(ATCC No. HTB 26), a human adenocarcinoma of the
breast, was maintained in McCoy’s 5A medium (Boehringer) containing L-glutamine, NaHCO, (2.3 g/liter),
gentamycin (50 mg/liter), and 10% NCS (GIBCO) in 75
cm2 flasks (Falcon Plastics 3023) as a monolayer.
The cells were serially passagedweekly following trypsinization using trypsin/EDTA
(Boehringer).
P388DI
(ATCC No. CCL 46), a murine lymphoid neoplasm, was grown in a suspension culture in Dulbecco’s
modified Eagle’s medium (Seromed) containing glucose
’ Abbreviations
used: crystal
violet,
N-hexamethylpararosaniline;
MTT,
3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2~-tetrazolium
mide; DMSO,
dimethyl
sulfoxide;
PBS, phosphate-buffered
OD, optical
density;
PC, personal
computer;
NCS, newborn
rum.
brosaline;
calf se-
0003-2697/90 $3.00
Copyright
0 1990 by Academic Press, Inc.
All rights of reproduction
in any form reserved.
COMPUTERIZED
I
DETERMINATION
Cell
Ii&e: MDA-MB-231
cnItum
conditions:
MCCOY’S
+ 10% NCS
humidlfkd
atmosphere.
37O
(96%
air. 5% coa 1
of cells
In
at defined
C
96well
density
Incubation
Fixation
of cells
with
1% glutaraldehyde
in PBS
for 15 min.
after
varying
times
of incubation.
Storage
under
PBS at 4O C
Staining
of 0.02%
of cells
crystal
Extraction
with
with
violet
a solution
(30 min.)
of bound
70% ethanol
GROWTH
KINETIC
263
CURVES
(I) MDA-MB-231
cells were processed as shown in
Flow Chart 1 by modifying the crystal violet staining
procedure described by Gillies et al. (6).
(II) The proliferation
of P388D1 cells was quantitated
according to a microculture
modification
of the tetrazolium assay (7-9, 11): MTT was prepared as a 1 mg/ml
stock solution in medium and filtered through a 0.22-pm
filter to remove undissolved
dye. At the times indicated
in Fig. 1, 50 ~1 of stock MTT solution was added to all
wells, each containing
100 ~1 of cell suspension.
The
plates were gently shaken and incubated at 37°C for another 1.5 h. Supernatant
removal was accomplished
by
carefully inverting
and blotting the tray; no prior centrifugation
was needed, since the macrophage-like
cells
sedimented and adhered to the substratum.
One hundred microliters
of DMSO was added to all wells to dissolve the formazan crystals.
Cytogenetic analysis.
The cells were grown to about
50% confluence on microscopic
slides. The slides were
prepared as described elsewhere
(12). So that spindle
I
1
Plating
nicmplates
OF
MAIN
stain
MENU
J
I
1
Measurement
of optical
(EL 309 Autoreader.
density
BIOTEK)
I
Computer-aided
data
analysis
DATA
FLOW
violet
dye
sensitivity
CHART
1.
Registration
procedure.
Optical
and reliability.
density
of growth
reading
curves
at 578
using
nm
gives
the
T O FILE
crystal
improved
(4.5 g/liter), NaHCO,
(3.7 g/liter), sodium pyruvate (110
mg/liter),
gentamycin
(50 mg/liter),
and 10% horse serum (Boehringer)
in 75cm2 culture flasks. The cells
were passaged every 3-4 days by 1:40 dilution with fresh
medium.
Both cell lines were grown at 37°C in a humidified atmosphere of 95% air and 5% carbon dioxide.
Mycoplasma
contamination
was monitored by routine
assay techniques
(Hoechst
33.258) (lo), and only mycoplasma-free
cultures were used.
Registration
of growth curves.
The assay was carried
out in flat-bottomed
microtitration
trays (Falcon Plastics 3075). Optical density was read using a Biotek EL
309 Autoreader
(Tecnomara)
at 578 nm (crystal violet)
and 540 nm (MTT).
POLYNOMIAL
LEAST SQUARES
FLOW
CHART
FIT
2.
Analysis
of spectrophotometric
data.
264
REILE
ET
AL.
60
1.5
50
0
0
,0
0
0”
10
B
aloo
0
0
0
60
120
180
"
60
TINE [hrs]
FIG.
1. (A) Calorimetric
measurement
of growth
in
were
unknown,
numbered
upon receipt from the ATCC)
MTT was added to all wells, and after another
1.5 h the
cell line MDA-MB-231
was plated in the 38th passage
medium/well
and were allowed
to attach.
Cell densities
microscope
(320X). At the times indicated,
the cells were
end of the experiment
all trays were stained with crystal
at 578 nm. (B) Plots of the corresponding
doubling
times,
= C Uiti;
(i = 0, 1, . . . , n)
Ui
This polynomial
= regression
180
microtitration
plates. (0) P388D1
leukemia
cells from passage 28 (original
passage
seeded at a density of 8 X 10’ cells/ml
(Coulter
Counter
ZsI). At the times indicated,
plates were developed
and read immediately
at 540 nm. (@) The human breast cancer
from origin at a density
of 12 cells per field of vision. The cells were seeded in 100 ~1
were estimated
by counting
50 fields of vision using a Leitz Diavert
phase-contrast
fixed with 1% glutaraldehyde
in PBS for 15 min and stored under PBS at 4°C. At the
violet simultaneously.
OD was measured
after extraction
of the dye with 70% ethanol
calculated
from the data shown in A, versus time of incubation.
formation
could be inhibited, the slide chambers were
inoculated with Colcemid solution (Serva) to a final concentration
of 0.04 pg/ml and incubated for 3 h at 37°C.
The medium was removed by suction and replaced with
0.075 M KCl. After 30 min of incubation at 37°C an equal
volume of cold, freshly made fixative (absolute methanol/glacial acetic acid, 3/l) was added. This hypotonic/
fixative mixture was removed immediately
and replaced
twice with ice-cold, fresh fixative. The slides were removed from the dish and air-dried.
The chromosomes
were stained for 8 min with 10 ml
Giemsa plus 90 ml of 0.025 M KH,PO,,
pH 6.8.
Computer processing.
Readings from the EL 309 autoreader (Biotek) were transferred
directly to an Olivetti
M24 personal computer, using a program that saved the
OD values on a diskette. The program for data processing is graphed in Flow Chart 2.
Doubling time analysis. Curve fitting of experimental
data is accomplished
by a polynomial regression fit applying the least-squares
method. Optical density is described as a function of time and data points are fitted
to an nth order polynomial
(usually n = 5) represented
by
OD = f(t)
120
TIME [hrs]
coefficients.
is defined in the time interval
of points
observed. Once the polynomial
fit has been performed
the doubling time TAU can easily be calculated at any
time of the interval concerned provided that exponential
growth is observed (13),
TAU
= l/cu
a! = l/in(2)
.d ln(OD)/dt
OD = C aiti,
and hence
TAU
= ln(2).(C
Uit”)/(C
i*oit”-‘);
(i = 0, 1, . . . , n).
ABS (TAU) is plotted versus time. From this plot the
time interval where exponential
growth occurs and the
corresponding
doubling time can easily be determined.
The software for evaluating cell kinetic data was developed in FORTRAN
and PASCAL
and implemented
on
an Olivetti
M24 pc.
RESULTS
Figure 1 illustrates
the growth curves and the corresponding doubling times at any time of incubation
for
an anchorage-dependent
(MDA-MB-231)
cell line and
an anchorage-independent
(P388D,)
cell line. The
curves were registered using two different staining tech-
COMPUTERIZED
#
100
DETERMINATION
,
l
80
T
5
60
*
7
P
3
J
x
40
20
0
0
60
120
180
TINE [hrs]
FIG. 2. Reproducibility
of doubling
time determination
by the crystal violet assay. MDA-MB-231
cells were plated at an approximate
density of 15 cells per field of vision (Leitz Diavert,
320X). The calculated doubling
times from four different
experiments
(all processed
separately)
are plotted
versus time of incubation.
niques:
crystal violet for MDA-MB-231
and MTT for
P388D1. Due to their macrophage-like
properties
(tendency to attach) no centrifugation
step was required in
the case of P388DI cells. A general application
of the
MTT
procedure
to anchorage-independent
cells requires a centrifuge equipped with a rotor capable of accepting multiplates.
Under the given experimental
conditions for both cell
lines the logarithmic
phase covers only a small fraction
of the overall growth curve. In a plot, doubling time
(TAU) versus time of incubation (t) (Fig. lB), exponential growth
is characterized
by a parallelism
of the
graphs with the t-axis. Although
in Fig. 1A the time
points are not equidistant
the calculated doubling times
for the exponential
growth phase are not significantly
affected by closely spaced points at either end of the
curve. A close scattering of data at both ends is required
to separate the log phase from the lag and the plateau
phases where cell proliferation
becomes extremely slow
resulting
in infinite
doubling times. Therefore
only
quantification
of the log phase is of biological significance.
MDA-MB-231
grows exponentially
for about two generations with a doubling time of -34 h, whereas the log
phase for P388DI at the given plating density is restricted to one doubling of the population.
By performing
the doubling time calculations
from
four successive crystal violet processing procedures utilizing MDA-MB-231
cells in the same passage seeded at
comparable plating densities, we encountered
an interassay variation of ~13% (Fig. 2).
OF
GROWTH
KINETIC
265
CURVES
The karyotype
information
presented in Fig. 3 gives
an example for the dramatic variability
of a well-defined
cell line after prolonged culture in vitro following cryopreservation.
Genotypic changes do not only reflect phenotypic characteristics
like morphology, but may also alter cell physiology
as well as cell proliferation,
which
directly affects doubling time (14).
Figure 4 illustrates the alterations
in doubling time as
a result of the number of passages. The doubling times of
MDA-MB-231
cells maintained under constant culture
conditions decreased constantly
with prolonged time in
culture. In the crystal violet assay the cells were seeded
at comparable densities. The doubling times vary enormously over a wide range from x10 h for a late (85th) to
=38 h for a relatively early (35th) passage. The variation most probably arises from the selection of faster
growing subclones from an initially heterogeneous population by frequent and incomplete trypsinization
of the
culture.
It is widely accepted that the initial cell number influences the proliferation
kinetics
of eukaryotic
cells.
This marked influence on doubling time is shown in
Fig. 5.
DISCUSSION
While a variety of techniques (6-10, 15) including an
enzyme immunoassay for estimating small cell numbers
(16) have been described recently for the determination
of parameters for cell proliferation, little interest was focused on the construction and analysis of growth curves
in animal cell culture.
/
25
/
I
I
/
80
90
100
/
n
2ot II
50
60
70
110
185
CHROMISOME NUMBER
FIG. 3.
from origin,
in passage
rethawed.
metaphase
Genetic
instability
of MDA-MB-231
cells. (Cl) 27th passage
(H) 35th passage from origin: The cells were cryopreserved
28 for about 6 months
and then passaged
weekly,
when
For each passage the chromsomes
from 50 well-spread
plates were counted.
266
REILE
It has been common practice to assume that cells exhibit ideal exponential
growth regardless of the diversities of experimental
designs without further verification.
In our opinion the description
of the overall growth
characteristics
with an exact determination
of the lag/
log phases, the doubling time (TAU),
and the plateau
level is absolutely
necessary.
In addition to the characterization
of cell lines (e.g., morphology,
karyotype,
receptor pattern, sensitivity
to antitumor
drugs) the description of the exponential growth phase is a prerequisite for any valuable experimental
setup. Especially, the
investigation
of cell lines underlying
extreme genetic
variability,
e.g., transformed
cells and cancer cells, requires a precise description
and verification
of the parameters concerned. Due to time-consuming
and difficult experimental
procedures
this is usually
not
performed.
The major disadvantage
of the methods described so
far (5, 17, 18) is the production
of mean values rather
than exact values. Due to laborious working techniques
only a few data points are generally obtained. This results in poor graphical estimates which make the separation of the lag/log/plateau
phases of growth curves extremely difficult and unsatisfactory.
In this paper we describe a technique which offers a
number
of advantages
compared
with conventional
methods. This inexpensive,
rapid, and technically
easy
microtitration
assay gives access to large sets of data
thus providing
results of high statistical
significance.
The equipment
necessary
(microplate
reader, PC) is
available in almost every analytical laboratory today.
loo-
80 -
ET
AL.
l
I
I
60
120
TIK
180
[hrs]
FIG.
5. Doubling
time as a function
of initial cell number.
MDAMB-231
cells were plated in passage 39 at different
densities:
(0) 10,
(A) 16, and (m) 36 cells per field of vision (Leitz Diavert,
320X).
Because of the mathematical
representation
of the
data an exact determination
of the log period is guaranteed. TAU as a function of time clearly reveals the duration of the lag and plateau phases. This approach overcomes the major drawback
of graphical methods which
only give a poor estimate of cell kinetics. This method
offers universal applicability
to biological investigations
whenever growth kinetics must be monitored. Its use is
by no way restricted
to calorimetric
measurement.
Its
application can also be extended to turbidimetric,
fluorometric,
radiochemical,
etc., monitoring
techniques.
The only prerequisite
is that a large number of data
which can be digitized are available.
ACKNOWLEDGMENTS
This work was supported
by the Deutsche
Forschungsgemeinschaft
and the “Matthias
Lackas-Stiftung
fiir Krebsforschung.”
Thanks
are
also due to the Fonds der Chemischen
Industrie
for financial
support.
The helpful assistance
of F. Birk in solving data transfer
problems
is
gratefully
acknowledged.
60 -
REFERENCES
1. Whang-Peng,
J., Lee, E. C., Kao-Shan,
Lippman,
M. (1983) J. Natl. Cancer Inst.
2. Yunis, J. J. (1983) Science 221,227-236.
0
60
120
TIE
FIG.
180
[hrs]
4.
Doubling
time as a function
of passage using MDA-MB-231
cells in the (0) 35th, (0) 39th, and (0, n ) 85th passage. (m) The cells
were plated at double density
(ca. 30 cells per field of vision,
Leitz
Diavert,
320X).
C.-S., Seibert,
71,687-695.
K., and
3. Seibert,
K., Shafie,
S. M., Triche,
T. J., Whang-Peng,
J. J.,
O’Brien,
S. J., Toney,
J. H., Huff, K. K., and Lippman,
M. E.
(1983) Cancer Res. 43,2223-2239.
4. Reddel, R. R., Alexander,
I. E., Koga, M., Shine, J., and Sutherland, R. L. (1988) Cancer Res. 48,4340-4347.
5. Reddel, R. R., Murphy,
L. C., Hall, R. E., and Sutherland,
R. L.
(1985) Cancer Res. 46,1525-1531.
COMPUTERIZED
6. Gillies,
R. J., Didier,
N., and Denton,
DETERMINATION
M.
(1986)
Anal.
OF
Biochem.
159,109-113.
7. Mosmann,
8. Denizot,
T. (1983)
F., and Lang,
J. Immunol.
R. (1986)
Methods
J. Immunol.
65,55-63.
Methods
89,
271-
277.
9. Alley, M. C., Scudiero,
D. A., Monks,
A., Hursey,
M. L., Czerwinski, M. J., Fine, D. L., Abbott,
B. J., Mayo,
J. G., Shoemaker,
R. H., and Boyd, M. R. (1988) Cancer Res. 48,589-601.
10. Peters, J. H., Baumgarten,
nale Antikorper:
Herstellung
Springer,
Berlin/Heidelberg/New
H., and Schulze, M. (1985) Monoklound Charakterisierung,
pp. 94-101,
York/Tokyo.
11. Carmichael,
J., DeGraff,
W. G., Gazdar,
A. F., Minna,
Mitchell,
J. B. (1987) Cancer Res. 47,936-942.
J. D., and
GROWTH
KINETIC
CURVES
267
12. Rooney,
D. E., and Czepulkowski
B. H. (1986) in Human
Cytogenetics, A Practical
Approach
(Rooney,
D. E., and Czepulkowski
B. H., Eds.), p. 10, IRL Press, Oxford/Washington,
DC.
13. Griffiths,
B. (1986) in Animal
Cell Culture,
Practical
Approach
(Freshney,
R. I., Ed.), p. 38, IRL Press, Oxford/Washington,
DC.
14. Osborne,
C. K., Hobbs, K., and Trent, J. M. (1987) Breast Cancer
Res. Treat. 9, 111-121.
15. Romijn,
J. C., Verkoelen,
C. F., and Schroeder,
F. H. (1988) Prostate 12,99-110.
16. Faraji-Shadan,
F., and Bowman,
P. D. (1989) Anal. Biochem.
177,259-262.
17. Roper, P. R., and Drewinko,
B. (1976) Cancer Res. 36,2182-2188.
18. Leonessa,
F., Coialbu,
T., and Toma, S. (1986) Anticancer
Res. 6,
1291-1296.