<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Adaptive Mitosis Detection in Large in vitro Stem Cell Populations using Timelapse Microscopy</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Tim Becker</string-name>
          <email>tim.becker@emb.fraunhofer.de</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Daniel H. Rapoport</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Amir Madany Mamlouk</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Fraunhofer Research Institution for Marine Biotechnology</institution>
          ,
          <addr-line>Lu ̈beck</addr-line>
          ,
          <country country="DE">Germany</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Graduate School for Computing in Medicine and Life Science, University of Lu ̈beck</institution>
          ,
          <country country="DE">Germany</country>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>Institute for Neuroand Bioinformatics, University of Lu ̈beck</institution>
          ,
          <country country="DE">Germany</country>
        </aff>
      </contrib-group>
      <fpage>49</fpage>
      <lpage>53</lpage>
      <abstract>
        <p>Reliable analysis of adult stem cell populations in in vitro experiments still poses a problem on the way to fully understand the regulating mechanism of these cultures. However, it is essential in the use of cultivated endogenous cells in stem cell therapies. One crucial feature during automated analysis is clearly the robust detection of mitotic events. In this work, we use the fully labeled stem cell benchmark data set CeTReS I in order to evaluate different approaches of mitosis detection: a purely time line based approach; a feature-based motility detector; and a detector based on the cell morphology changes, for which we also propose an adaptive version. We demonstrate that the approach based on morphological change outperforms the static detectors. However, the set of optimal features is changing over time, and thus it is not surprising that a feature set adapted to the systems confluency shows the best performance.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        In order to fully understand the behavior and regulating mechanisms of in vitro
stem cell populations, time-lapse phase contrast microscopy plays a major role.
This technique enables the observation and analysis of a cell population in a
completely label-free and noninvasive manner and allows a simultaneous
processing of thousands of cells in vitro. An automated cell detection and tracking
system facilitates the analysis of these data and represents a big improvement
to the actual gold standard method: the tedious work of manual cell tracking.
The key event in analyzing this data is cell division or mitosis. The mitosis is
the most difficult part to detect as it implies a strong non-linearity or, in other
words, new cells and tracking candidates are added to the scenario. However, an
accurate knowledge of these events enables the reconstruction of cell lineages [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]
as well as the complete characterization of cell populations including life time
distributions, change of cell cycle length and (a)symmetry of cell divisions. Up
to now, only a few approaches explicitly address the detection of mitotic events
in unlabeled cell populations. Those that do typically use a spatio-temporal
pattern that arises during mitosis. As an example, the radius of a cell can be used
to determine such an event [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. Other approaches use a characteristic change of
different cell parameters, such as size, brightness and circularity [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ].
2
2.1
      </p>
    </sec>
    <sec id="sec-2">
      <title>Material and Methods</title>
      <sec id="sec-2-1">
        <title>Data</title>
        <p>
          All detection and tracking tasks were performed using the freely available data
set CeTReS I [
          <xref ref-type="bibr" rid="ref4">4</xref>
          ]. This reference set contains two image series of adult stem
cells from rat pancreas [
          <xref ref-type="bibr" rid="ref5">5</xref>
          ]. The first series consists of 209 images acquired
every 15 minutes and shows the complete growth curve in detail: the initial lag,
exponential and stationary phase (Fig. 1). This image series shows the biggest
variety of cell characteristics and was therefore chosen for testing and evaluation.
2.2
        </p>
      </sec>
      <sec id="sec-2-2">
        <title>Cell Detection and Tracking</title>
        <p>
          Within this study, we use a straightforward algorithm to detect the cells
automatically on the given microscope images. To get a consistent input, all images
were preprocessed in two steps: first, an image normalization was used to expand
the dynamic range of the data; second, the contrast was increased by performing
a histogram equalization. To perform cell detection, each image is segmented
using the Otsu intensity threshold. The result of this operation is a binary
image that is further refined by morphological filters. This algorithm achieves a
false acceptance rate (FAR) and false rejection rate (FRR) of 5% for cell
detection and thus give a detection quality comparable with the state-of-the-art
approaches [
          <xref ref-type="bibr" rid="ref1">1</xref>
          ]. The next task consists of following the cells over time, i.e. to
track the cells. This task is performed using the overlap of the detected cell
shapes in consecutive frames. As long as a unique cell successor can be
determined, the cell shapes are used to construct cell paths p, i.e. a cell path consists
of a series of a uniquely merged time series of cells. To correctly connect these
paths and to reconstruct complete genealogical trees, it is essential that a path
starts and ends with a mitosis.
2.3
        </p>
      </sec>
      <sec id="sec-2-3">
        <title>Mitosis Detection</title>
        <p>
          In general, there are two main approaches to detect mitoses: first, through
temporal context; second, through a spatio-temporal pattern. Using the
temporal context, a mitosis can be found as one cell (mother) divides into two cells
(daughters). In a binary tree, this gives the typical Y-shape in a family tree.
Furthermore, mitotic cells feature a characteristic spatio-temporal pattern
(illustrated in Fig. 2). In particular, they contract, become smaller, rounder, appear
brighter and decrease their velocity. Based on these criteria, four detectors were
implemented for comparison:
{ Y detector: The naive approach classifies each path p as mitotic if it can be
linked to two successive cell paths and both have p as their solely predecessor.
{ Motility detector The second detector is motivated by Debeir et al. [
          <xref ref-type="bibr" rid="ref2">2</xref>
          ]. It
simply classifies a cell path p as being mitotic if both, size S and circularity
C lie below a predefined threshold .
{ Morphology detector: The following set of cell features is computed for all
detected single cells: perimeter P , cell area A, circularity C, average intensity
(brightness) B and the length L. The changes of these parameters, denoted
as d, are calculated for all time points of a cell path p. A path p is classified
as being mitotic if the multidimensional normal distribution Prob(d) exceeds
a threshold ; Prob(d) is given by
        </p>
        <p>Prob(d) =</p>
        <p>1
√(2 )4det( )
exp
(
(d
1 (d</p>
        <p>) )
)T
2
(1)
where denotes the covariance matrix and the mean feature difference
vector. Values for and are trained using the reference data. These
features combine to a set of potential classifiers c(!; d), with ! 2 f0; 1gN
indicating which of the N features are turned on or off.
{ Adaptive morphology detector Finally, we propose a detector that substitutes
! in c(!; d) by a time-dependent feature parameter !(t). This appears to be
a reasonable extension as single cell morphology is strongly changing over
time due to population effects, e.g., confluency.
3</p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>Results</title>
      <p>The reference data CeTReS I was used to evaluate the proposed detectors. This
set contains a total of 4689 cell paths and 2019 labeled mitotic events, of which
1903 (95%) are successfully found by our proposed detection and tracking
algorithm (2.2). To compare the different detection results, the false acceptance rate
(FAR=specificity) and false rejection rate (FRR=sensitivity) for each detector
was calculated with varying thresholds. In Figure 3(a), the ROC (Receiver
Operator Characteristic) curves are shown for the different classifiers. The motility
detection clearly performs worst. The motility detector is outperformed by the
morphology detector using only size as a feature. Note that the motility variant
use cell size as a static measure while the morphology detector is working on size
differences. As shown, the combination of brightness and size parameter resulted
in the best performance for the morphology detector. Figure 3(b) illustrates the
behavior of the different detection strategies over time. The motility detection
shows a strong decrease in performance leading to the worst result. Some of the
morphology features show a contrary behavior: size has a decreasing performance
while compactness has an increasing performace over time. This illustrates why
the adaptive detection strategy performs best. From the biological point of view,
(a) ROC for different detections.
(b) EER over time
detection algorithm
complete paths
reliability
the most important measure is the rate of correct detected complete paths, i.e.
paths that describe complete cell cycles from the “birth” mitosis to the “end”
mitosis. The results presented in table 1 show the number of detected paths and
their reliability. To achieve the most reliable results, the Y detection was used
to preselect possible mitosis candidates.
4</p>
    </sec>
    <sec id="sec-4">
      <title>Discussion</title>
      <p>
        In this paper, we present a mitosis detection based on the morphological change
of the cells during division. The novelty of this algorithm lies in the fact that it
adapts to the changes of cell appearance during culturing. Therefore, it
outperforms other, non adaptive approaches like those of [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ] or [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. One has to bear
in mind, that the different mitosis detection approaches have been tried on
different time-lapse data. Therefore, a real quantitative comparison is not possible
yet. We suggest to establish a common benchmark data set, like the CeTReS [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]
timelapse series. Summarizing, the presented data illustrate the varying and
inhomogeneous character of adult stem cell populations. This character has to
be taken into account, when a fully automated cell tracking shall be used for
future automated cell culture systems.
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Li</surname>
            <given-names>K</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chen</surname>
            <given-names>M</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kanade</surname>
            <given-names>T</given-names>
          </string-name>
          , et al.
          <article-title>Cell population tracking and lineage construction with spatiotemporal context</article-title>
          .
          <source>Med Image Anal</source>
          .
          <year>2008</year>
          ;
          <volume>12</volume>
          (
          <issue>1</issue>
          ):
          <fpage>546</fpage>
          -
          <lpage>66</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Debeir</surname>
            <given-names>O</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Milojevic</surname>
            <given-names>D</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Leloup</surname>
            <given-names>T</given-names>
          </string-name>
          , et al.
          <article-title>Mitotic tree construction by computer in vitro cell tracking: a tool for proliferation and motility features extraction</article-title>
          .
          <source>In: Proc EUROCON; 2005</source>
          . p.
          <fpage>951</fpage>
          -
          <lpage>4</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>yang</surname>
            <given-names>F</given-names>
          </string-name>
          ,
          <article-title>michael a mackey, fiorenza ianzini</article-title>
          , et al.
          <article-title>Cell segmentation, tracking, and mitosis detection using temporal context</article-title>
          .
          <source>In: Proc MICCAI; 2005</source>
          . p.
          <fpage>302</fpage>
          -
          <lpage>9</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Madany-Mamlouk</surname>
            <given-names>A</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Becker</surname>
            <given-names>T</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Rapoport</surname>
            <given-names>D.</given-names>
          </string-name>
          <article-title>A novel benchmark data set for adult stem cell tracking in time-lapse microscopy</article-title>
          . In: Front Neur Conference Abstract: Neuroinformatics;
          <year>2009</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Kruse</surname>
            <given-names>C</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kajahn</surname>
            <given-names>J</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Petschnik</surname>
            <given-names>AE</given-names>
          </string-name>
          , et al.
          <article-title>Adult pancreatic stem/progenitor cells spontaneously differentiate in vitro into multiple cell lineages and form teratomalike structures</article-title>
          .
          <source>Ann Anat</source>
          .
          <year>2006</year>
          ;
          <volume>188</volume>
          :
          <fpage>503</fpage>
          -
          <lpage>17</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>