Computational prediction of neural progenitor cell fates

Single-shot Fourier ptychographic microscopy with isotropic lateral resolution via polarization-multiplexed LED illumination

Fourier ptychographic microscopy (FPM) has emerged as a new wide-field and high-resolution computational imaging technique in recent years. To ensure data redundancy for a stable convergence solution, conventional FPM requires dozens or hundreds of raw images, increasing the time cost for both data collection and computation. Here, we propose a single-shot Fourier ptychographic microscopy with isotropic lateral resolution via polarization-multiplexed LED illumination, termed SIFPM. Three LED elements covered with 0°/45°/135° polarization films, respectively, are used to provide numerical aperture-matched illumination for the sample simultaneously. Meanwhile, a polarization camera is utilized to record the light field distribution transmitted through the sample. Based on weak object transfer functions, we first obtain the amplitude and phase estimations of the sample by deconvolution, and then we use them as the initial guesses of the FPM algorithm to refine the accuracy of reconstruction. We validate the complex sample imaging performance of the proposed method on quantitative phase target, unstained and stained bio-samples. These results show that SIFPM can realize quantitative imaging for general samples with the resolution of the incoherent diffraction limit, permitting high-speed quantitative characterization for cells and tissues.

The Cluster Structure Function

For each partition of a data set into a given number of parts there is a partition such that every part is as much as possible a good model (an "algorithmic sufficient statistic") for the data in that part. Since this can be done for every number between one and the number of data, the result is a function, the cluster structure function. It maps the number of parts of a partition to values related to the deficiencies of being good models by the parts. Such a function starts with a value at least zero for no partition of the data set and descents to zero for the partition of the data set into singleton parts. The optimal clustering is the one selected by analyzing the cluster structure function. The theory behind the method is expressed in algorithmic information theory (Kolmogorov complexity). In practice the Kolmogorov complexities involved are approximated by a concrete compressor. We give examples using real data sets: the MNIST handwritten digits and the segmentation of real cells as used in stem cell research.

Machine Learning in Tissue Engineering

Machine learning (ML) and artificial intelligence have accelerated scientific discovery, augmented clinical practice, and deepened fundamental understanding of many biological phenomena. ML technologies have now been applied to diverse areas of tissue engineering research, including biomaterial design, scaffold fabrication, and cell/tissue modeling. Emerging ML-empowered strategies include machine-optimized polymer synthesis, predictive modeling of scaffold fabrication processes, complex analyses of structure-function relationships, and deep learning of spatialized cell phenotypes and tissue composition. The emergence of ML in tissue engineering, while relatively recent, has already enabled increasingly complex and multivariate analyses of the relationships between biological, chemical, and physical factors in driving tissue regenerative outcomes. This review highlights the novel methodologies, emerging strategies, and areas of potential growth within this rapidly evolving area of research. Impact statement Machine learning (ML) has accelerated scientific discovery and augmented clinical practice across multiple fields. Now, ML has driven exciting new paradigms in tissue engineering research, including machine-optimized biomaterial design, predictive modeling of scaffold fabrication, and spatiotemporal analysis of cell and tissue systems. The emergence of ML in tissue engineering, while relatively recent, has already enabled increasingly complex analyses of the relationships between biological, chemical, and physical factors in driving tissue regenerative outcomes. This review highlights the novel methodologies, emerging strategies, and areas of potential growth within this rapidly evolving area of research.

Deep Phenotypic Cell Classification using Capsule Neural Network

Recent developments in ultra-high-throughput microscopy have created a new generation of cell classification methodologies focused solely on image-based cell phenotypes. These image-based analyses enable morphological profiling and screening of thousands or even millions of single cells at a fraction of the cost. They have been shown to demonstrate the statistical significance required for understanding the role of cell heterogeneity in diverse biologists. However, these single-cell analysis techniques are slow and require expensive genetic/epigenetic analysis. This treatise proposes an innovative DL system based on the newly created capsule networks (CapsNet) architecture. The proposed deep CapsNet model employs "Capsules" for high-level feature abstraction relevant to the cell category. Experiments demonstrate that our proposed system can accurately classify different types of cells based on phenotypic label-free bright-field images with over 98.06% accuracy and that deep CapsNet models outperform CNN models in the prior art.

A single-cell guide to retinal development: Cell fate decisions of multipotent retinal progenitors in scRNA-seq

Recent advances in high throughput single-cell RNA sequencing (scRNA-seq) technology have enabled the simultaneous transcriptomic profiling of thousands of individual cells in a single experiment. To investigate the intrinsic process of retinal development, researchers have leveraged this technology to quantify gene expression in retinal cells across development, in multiple species, and from numerous important models of human disease. In this review, we summarize recent applications of scRNA-seq and discuss how these datasets have complemented and advanced our understanding of retinal progenitor cell competence, cell fate specification, and differentiation. Finally, we also highlight the outstanding questions in the field that advances in single-cell data generation and analysis will soon be able to answer.

Emerging machine learning approaches to phenotyping cellular motility and morphodynamics

Cells respond heterogeneously to molecular and environmental perturbations. Phenotypic heterogeneity, wherein multiple phenotypes coexist in the same conditions, presents challenges when interpreting the observed heterogeneity. Advances in live cell microscopy allow researchers to acquire an unprecedented amount of live cell image data at high spatiotemporal resolutions. Phenotyping cellular dynamics, however, is a nontrivial task and requires machine learning (ML) approaches to discern phenotypic heterogeneity from live cell images. In recent years, ML has proven instrumental in biomedical research, allowing scientists to implement sophisticated computation in which computers learn and effectively perform specific analyses with minimal human instruction or intervention. In this review, we discuss how ML has been recently employed in the study of cell motility and morphodynamics to identify phenotypes from computer vision analysis. We focus on new approaches to extract and learn meaningful spatiotemporal features from complex live cell images for cellular and subcellular phenotyping.

Deep Convolutional and Recurrent Neural Networks for Cell Motility Discrimination and Prediction

Cells in culture display diverse motility behaviors that may reflect differences in cell state and function, providing motivation to discriminate between different motility behaviors. Current methods to do so rely upon manual feature engineering. However, the types of features necessary to distinguish between motility behaviors can vary greatly depending on the biological context, and it is not always clear which features may be most predictive in each setting for distinguishing particular cell types or disease states. Convolutional neural networks (CNNs) are machine learning models allowing for relevant features to be learned directly from spatial data. Similarly, recurrent neural networks (RNNs) are a class of models capable of learning long term temporal dependencies. Given that cell motility is inherently spacio-temporal data, we present an approach utilizing both convolutional and long- short-term memory (LSTM) recurrent neural network units to analyze cell motility data. These RNN models provide accurate classification of simulated motility and experimentally measured motility from multiple cell types, comparable to results achieved with hand-engineered features. The variety of cell motility differences we can detect suggests that the algorithm is generally applicable to additional cell types not analyzed here. RNN autoencoders based on the same architecture are capable of learning motility features in an unsupervised manner and capturing variation between myogenic cells in the latent space. Adapting these RNN models to motility prediction, RNNs are capable of predicting muscle stem cell motility from past tracking data with performance superior to standard motion prediction models. This advance in cell motility prediction may be of practical utility in cell tracking applications.

Dynamic Polarization of Rab11a Modulates Crb2a Localization and Impacts Signaling to Regulate Retinal Neurogenesis

Interkinetic nuclear migration (IKNM) is the process in which pseudostratified epithelial nuclei oscillate from the apical to basal surface and in phase with the mitotic cycle. In the zebrafish retina, neuroepithelial retinal progenitor cells (RPCs) increase Notch activity with apical movement of the nuclei, and the depth of nuclear migration correlates with the probability that the next cell division will be neurogenic. This study focuses on the mechanisms underlying the relationships between IKNM, cell signaling, and neurogenesis. In particular, we have explored the role IKNM has on endosome biology within RPCs. Through genetic manipulation and live imaging in zebrafish, we find that early (Rab5-positive) and recycling (Rab11a-positive) endosomes polarize in a dynamic fashion within RPCs and with reference to nuclear position. Functional analyses suggest that dynamic polarization of recycling endosomes and their activity within the neuroepithelia modulates the subcellular localization of Crb2a, consequently affecting multiple signaling pathways that impact neurogenesis including Notch, Hippo, and Wnt activities. As nuclear migration is heterogenous and asynchronous among RPCs, Rab11a-affected signaling within the neuroepithelia is modulated in a differential manner, providing mechanistic insight to the correlation of IKNM and selection of RPCs to undergo neurogenesis.

Nessys: A new set of tools for the automated detection of nuclei within intact tissues and dense 3D cultures

Methods for measuring the properties of individual cells within their native 3D environment will enable a deeper understanding of embryonic development, tissue regeneration, and tumorigenesis. However, current methods for segmenting nuclei in 3D tissues are not designed for situations in which nuclei are densely packed, nonspherical, or heterogeneous in shape, size, or texture, all of which are true of many embryonic and adult tissue types as well as in many cases for cells differentiating in culture. Here, we overcome this bottleneck by devising a novel method based on labelling the nuclear envelope (NE) and automatically distinguishing individual nuclei using a tree-structured ridge-tracing method followed by shape ranking according to a trained classifier. The method is fast and makes it possible to process images that are larger than the computer's memory. We consistently obtain accurate segmentation rates of >90%, even for challenging images such as mid-gestation embryos or 3D cultures. We provide a 3D editor and inspector for the manual curation of the segmentation results as well as a program to assess the accuracy of the segmentation. We have also generated a live reporter of the NE that can be used to track live cells in 3 dimensions over time. We use this to monitor the history of cell interactions and occurrences of neighbour exchange within cultures of pluripotent cells during differentiation. We provide these tools in an open-access user-friendly format.

Separating Touching Cells Using Pixel Replicated Elliptical Shape Models

One of the most important and error-prone tasks in biological image analysis is the segmentation of touching or overlapping cells. Particularly for optical microscopy, including transmitted light and confocal fluorescence microscopy, there is often no consistent discriminative information to separate cells that touch or overlap. It is desired to partition touching foreground pixels into cells using the binary threshold image information only, and optionally incorporating gradient information. The most common approaches for segmenting touching and overlapping cells in these scenarios are based on the watershed transform. We describe a new approach called pixel replication for the task of segmenting elliptical objects that touch or overlap. Pixel replication uses the image Euclidean distance transform in combination with Gaussian mixture models to better exploit practically effective optimization for delineating objects with elliptical decision boundaries. Pixel replication improves significantly on commonly used methods based on watershed transforms, or based on fitting Gaussian mixtures directly to the thresholded image data. Pixel replication works equivalently on both 2-D and 3-D image data, and naturally combines information from multi-channel images. The accuracy of the proposed technique is measured using both the segmentation accuracy on simulated ellipse data and the tracking accuracy on validated stem cell tracking results extracted from hundreds of live-cell microscopy image sequences. Pixel replication is shown to be significantly more accurate compared with other approaches. Variance relationships are derived, allowing a more practically effective Gaussian mixture model to extract cell boundaries for data generated from the threshold image using the uniform elliptical distribution and from the distance transform image using the triangular elliptical distribution.

Induction of rod versus cone photoreceptor-specific progenitors from retinal precursor cells

During development, multipotent progenitors undergo temporally-restricted differentiation into post-mitotic retinal cells; however, the mechanisms of progenitor division that occurs during retinogenesis remain controversial. Using clonal analyses (lineage tracing and single cell cultures), we identify rod versus cone lineage-specific progenitors derived from both adult retinal stem cells and embryonic neural retinal precursors. Taurine and retinoic acid are shown to act in an instructive and lineage-restricted manner early in the progenitor lineage hierarchy to produce rod-restricted progenitors from stem cell progeny. We also identify an instructive, but lineage-independent, mechanism for the specification of cone-restricted progenitors through the suppression of multiple differentiation signaling pathways. These data indicate that exogenous signals play critical roles in directing lineage decisions and resulting in fate-restricted rod or cone photoreceptor progenitors in culture. Additional factors may be involved in governing photoreceptor fates in vivo.

Subwavelength resolution Fourier ptychography with hemispherical digital condensers

Fourier ptychography (FP) is a promising computational imaging technique that overcomes the physical space-bandwidth product (SBP) limit of a conventional microscope by applying angular-varied illuminations. However, to date, the effective imaging numerical aperture (NA) achievable with a commercial LED board is still limited to the range of 0.3-0.7 with a 4 × /0.1NA objective due to the geometric constraint with the declined illumination intensities and attenuated signal-to-noise ratio (SNR). Thus the highest achievable half-pitch resolution is usually constrained between 500-1000 nm, which cannot meet the requirements of high-resolution biomedical imaging applications. Although it is possible to improve the resolution by using a high-NA objective lens, the FP approach is less appealing as the decrease of field-of-view (FOV) will far exceed the improvement of spatial resolution in this case. In this paper, we initially present a subwavelength resolution Fourier ptychography (SRFP) platform with a hemispherical digital condenser to provide high-angle programmable plane-wave illuminations of 0.95NA, attaining a 4 × /0.1NA objective with the final effective imaging performance of 1.05NA at a half-pitch resolution of 244 nm with the incident wavelength of 465 nm across a wide FOV of 14.60 mm², corresponding to a SBP of 245 megapixels. Our work provides an essential step of FP towards high-throughput imaging applications.

Spiral wave classification using normalized compression distance: Towards atrial tissue spatiotemporal electrophysiological behavior characterization

Analysis of electrical activation patterns such as re-entries during atrial fibrillation (Afib) is crucial in understanding arrhythmic mechanisms and assessment of diagnostic measures. Spiral waves are a phenomena that provide intuitive basis for re-entries occurring in cardiac tissue. Distinct spiral wave behaviors such as stable spiral waves, meandering spiral waves, and spiral wave break-up may have distinct electrogram manifestations on a mapping catheter. Hence, it is desirable to have an automated classification of spiral wave behavior based on catheter recordings for a qualitative characterization of spatiotemporal electrophysiological activity on atrial tissue. In this study, we propose a method for classification of spatiotemporal characteristics of simulated atrial activation patterns in terms of distinct spiral wave behaviors during Afib using two different techniques: normalized compressed distance (NCD) and normalized FFT (NFFTD). We use a phenomenological model for cardiac electrical propagation to produce various simulated spiral wave behaviors on a 2D grid and labeled them as stable, meandering, or breakup. By mimicking commonly used catheter types, a star shaped and a circular shaped both of which do the local readings from atrial wall, monopolar and bipolar intracardiac electrograms are simulated. Virtual catheters are positioned at different locations on the grid. The classification performance for different catheter locations, types and for monopolar or bipolar readings were also compared. We observed that the performance for each case differed slightly. However, we found that NCD performance is superior to NFFTD. Through the simulation study, we showed the theoretical validation of the proposed method. Our findings suggest that a qualitative wavefront activation pattern can be assessed during Afib without the need for highly invasive mapping techniques such as multisite simultaneous electrogram recordings.

ATLANTIS - Attractor Landscape Analysis Toolbox for Cell Fate Discovery and Reprogramming

Boolean modelling of biological networks is a well-established technique for abstracting dynamical biomolecular regulation in cells. Specifically, decoding linkages between salient regulatory network states and corresponding cell fate outcomes can help uncover pathological foundations of diseases such as cancer. Attractor landscape analysis is one such methodology which converts complex network behavior into a landscape of network states wherein each state is represented by propensity of its occurrence. Towards undertaking attractor landscape analysis of Boolean networks, we propose an Attractor Landscape Analysis Toolbox (ATLANTIS) for cell fate discovery, from biomolecular networks, and reprogramming upon network perturbation. ATLANTIS can be employed to perform both deterministic and probabilistic analyses. It has been validated by successfully reconstructing attractor landscapes from several published case studies followed by reprogramming of cell fates upon therapeutic treatment of network. Additionally, the biomolecular network of HCT-116 colorectal cancer cell line has been screened for therapeutic evaluation of drug-targets. Our results show agreement between therapeutic efficacies reported by ATLANTIS and the published literature. These case studies sufficiently highlight the in silico cell fate prediction and therapeutic screening potential of the toolbox. Lastly, ATLANTIS can also help guide single or combinatorial therapy responses towards reprogramming biomolecular networks to recover cell fates.

Inferring cell state by quantitative motility analysis reveals a dynamic state system and broken detailed balance

Cell populations display heterogeneous and dynamic phenotypic states at multiple scales. Similar to molecular features commonly used to explore cell heterogeneity, cell behavior is a rich phenotypic space that may allow for identification of relevant cell states. Inference of cell state from cell behavior across a time course may enable the investigation of dynamics of transitions between heterogeneous cell states, a task difficult to perform with destructive molecular observations. Cell motility is one such easily observed cell behavior with known biomedical relevance. To investigate heterogenous cell states and their dynamics through the lens of cell behavior, we developed Heteromotility, a software tool to extract quantitative motility features from timelapse cell images. In mouse embryonic fibroblasts (MEFs), myoblasts, and muscle stem cells (MuSCs), Heteromotility analysis identifies multiple motility phenotypes within the population. In all three systems, the motility state identity of individual cells is dynamic. Quantification of state transitions reveals that MuSCs undergoing activation transition through progressive motility states toward the myoblast phenotype. Transition rates during MuSC activation suggest non-linear kinetics. By probability flux analysis, we find that this MuSC motility state system breaks detailed balance, while the MEF and myoblast systems do not. Balanced behavior state transitions can be captured by equilibrium formalisms, while unbalanced switching between states violates equilibrium conditions and would require an external driving force. Our data indicate that the system regulating cell behavior can be decomposed into a set of attractor states which depend on the identity of the cell, together with a set of transitions between states. These results support a conceptual view of cell populations as dynamical systems, responding to inputs from signaling pathways and generating outputs in the form of state transitions and observable motile behaviors.

Identifying suitable reference genes for developing and injured mouse CNS tissues

Accurate quantification of gene expression is fundamental for understanding the molecular, genetic and functional bases of tissue development and diseases. Quantitative real-time PCR (qPCR) is now the most widely used method of quantifying gene expression due to its simplicity, specificity, sensitivity, and wide quantification range. The use of appropriate reference genes to ensure accurate normalization is crucial for the correct quantification of gene expression from the early development, maturation, aging to injury processes in the central nervous system (CNS). In this study, we have determined the expression profiles of 12 candidate housekeeping genes (ACTB, CYC1, HMBS, GAPDH, HPRT1, RPL13A, YWHAZ, PPIA, RPLP0, TFRC, GUS, and 18S rRNA) in developing mouse brain and spinal cord. Throughout development, there was a significant degree of fluctuations in their expression levels, indicating the importance and complexity of finding appropriate reference genes. Three software including BestKeeper, geNorm and NormFinder were used to evaluate the stability of potential reference genes. GUS was the most stable gene and GUS/YWHAZ were the most stable reference gene pair across different developmental stages in different CNS regions, whereas HPRT1 and GAPDH were the most variable genes and thus inappropriate to use as reference genes. Therefore, our results identified GUS and YWHAZ as the best combination of two reference genes for expression data normalization in CNS developmental studies. © 2017 Wiley Periodicals, Inc. Develop Neurobiol 78: 39-50, 2018.

Non-monotonic Changes in Progenitor Cell Behavior and Gene Expression during Aging of the Adult V-SVZ Neural Stem Cell Niche

Neural stem cell activity in the ventricular-subventricular zone (V-SVZ) decreases with aging, thought to occur by a unidirectional decline. However, by analyzing the V-SVZ transcriptome of male mice at 2, 6, 18, and 22 months, we found that most of the genes that change significantly over time show a reversal of trend, with a maximum or minimum expression at 18 months. In vivo, MASH1⁺ progenitor cells decreased in number and proliferation between 2 and 18 months but increased between 18 and 22 months. Time-lapse lineage analysis of 944 V-SVZ cells showed that age-related declines in neurogenesis were recapitulated in vitro in clones. However, activated type B/type C cell clones divide slower at 2 to 18 months, then unexpectedly faster at 22 months, with impaired transition to type A neuroblasts. Our findings indicate that aging of the V-SVZ involves significant non-monotonic changes that are programmed within progenitor cells and are observable independent of the aging niche.

•

RNA sequencing analysis of the adult V-SVZ NSC niche at 2, 6, 18, and 22 months

•

During aging, most V-SVZ niche genes show max/min expression at 18 months

•

In vivo MASH1⁺ cells cycle slowest at 18 months but at 22 months return to 2-month rate

•

Time-lapse analyses of isolated SVZ cells show that age-associated changes are programmed

ECG data compression using a neural network model based on multi-objective optimization

Electrocardiogram (ECG) data analysis is of great significance to the diagnosis of cardiovascular disease. ECG compression should be processed in real time, and the data should be based on lossless compression and have high predictability. In terms of the real time aspect, short-time Fourier transformation is applied to the processing of signal wave for reducing computational time. For the lossless compression requirement, wavelet-transformation that is a coding algorithm can be used to avoid loss of data. In practice, compression is required to avoid storing redundant recording data that are not useful in the diagnosis platform. The obtained data can be preprocessed to remove noise by using wavelet transform, and then a multi-objective optimize neural network model is used to extract feature information. Compared with the existing traditional methods such as direct data processing method and transform method, our proposed compression model has self-learning ability to achieve high data compression ratio at 1:19 without losing important ECG information and compromising quality. Upon testing, we demonstrated that the proposed ECG data compression method based on multi-objective optimization neural network is effective and efficient in clinical practice.

Molecular Dynamics Validation of Crizotinib Resistance to ALK Mutations (L1196M and G1269A) and Identification of Specific Inhibitors

Anaplastic lymphoma kinase (ALK) positive non-small cell lung cancer (NSCLC) patients are mostly treated with ALK tyrosine kinase inhibitors (TKIs). Crizotinib is the first generation ALK inhibitor practiced as a primary chemo to combat cancer cells followed by second generation inhibitor ceritinib which are effective against crizotinib resistant ALK mutations. However, patients treated with these drugs invariably relapsed because of the development of new drug resistance mutations. In this study we explored the crizotinib resistance in the presence of ALK mutations L1196M and G1269A through molecular dynamics simulation studies. Further mutation specific inhibitors CID 71748211 and CID 71728095 were identified to potentially inhibit ALK with mutations L1196M and G1269A, respectively. This computational investigation in-sighted the molecular factors involved in crizotinib resistance which enhanced in the identification of new ALK drugs that brings individualized medicine to treat ALK positive NSCLC patients with specific mutations. J. Cell. Biochem. 118: 3462-3471, 2017. © 2017 Wiley Periodicals, Inc.

Resolution-enhanced Fourier ptychographic microscopy based on high-numerical-aperture illuminations

High-resolution and wide field-of-view (FOV) microscopic imaging plays a central role in diverse applications such as high-throughput screening and digital pathology. However, conventional microscopes face inherent trade-offs between the spatial resolution and FOV, which are fundamental limited by the space-bandwidth product (SBP) of the optical system. The resolution-FOV tradeoff can be effectively decoupled in Fourier ptychography microscopy (FPM), however, to date, the effective imaging NA achievable with a typical FPM system is still limited to the range of 0.4–0.7. Herein, we report, for the first time, a high-NA illumination based resolution-enhanced FPM (REFPM) platform, in which a LED-array-based digital oil-immersion condenser is used to create high-angle programmable plane-wave illuminations, endowing a 10×, 0.4 NA objective lens with final effective imaging performance of 1.6 NA. With REFPM, we present the highest-resolution results with a unprecedented half-pitch resolution of 154 nm at a wavelength of 435 nm across a wide FOV of 2.34 mm², corresponding to an SBP of 98.5 megapixels (~50 times higher than that of the conventional incoherent microscope with the same resolution). Our work provides an important step of FPM towards high-resolution large-NA imaging applications, generating comparable resolution performance but significantly broadening the FOV of conventional oil-immersion microscopes.

Machine learning applications in cell image analysis

Machine learning (ML) refers to a set of automatic pattern recognition methods that have been successfully applied across various problem domains, including biomedical image analysis. This review focuses on ML applications for image analysis in light microscopy experiments with typical tasks of segmenting and tracking individual cells, and modelling of reconstructed lineage trees. After describing a typical image analysis pipeline and highlighting challenges of automatic analysis (for example, variability in cell morphology, tracking in presence of clutters) this review gives a brief historical outlook of ML, followed by basic concepts and definitions required for understanding examples. This article then presents several example applications at various image processing stages, including the use of supervised learning methods for improving cell segmentation, and the application of active learning for tracking. The review concludes with remarks on parameter setting and future directions.

Prospective identification of hematopoietic lineage choice by deep learning

Differentiation alters molecular properties of stem and progenitor cells, leading to changing shape and movement characteristics. We present a deep neural network that prospectively predicts lineage choice in differentiating primary hematopoietic progenitors, using image patches from brightfield microscopy and cellular movement. Surprisingly, lineage choice can be detected up to three generations before conventional molecular markers are observable. Our approach allows identifying cells with differentially expressed lineage-specifying genes without molecular labeling.

Lineage mapper: A versatile cell and particle tracker

The ability to accurately track cells and particles from images is critical to many biomedical problems. To address this, we developed Lineage Mapper, an open-source tracker for time-lapse images of biological cells, colonies, and particles. Lineage Mapper tracks objects independently of the segmentation method, detects mitosis in confluence, separates cell clumps mistakenly segmented as a single cell, provides accuracy and scalability even on terabyte-sized datasets, and creates division and/or fusion lineages. Lineage Mapper has been tested and validated on multiple biological and simulated problems. The software is available in ImageJ and Matlab at isg.nist.gov.

Measuring Process Dynamics and Nuclear Migration for Clones of Neural Progenitor Cells

Neural stem and progenitor cells (NPCs) generate processes that extend from the cell body in a dynamic manner. The NPC nucleus migrates along these processes with patterns believed to be tightly coupled to mechanisms of cell cycle regulation and cell fate determination. Here, we describe a new segmentation and tracking approach that allows NPC processes and nuclei to be reliably tracked across multiple rounds of cell division in phase-contrast microscopy images. Results are presented for mouse adult and embryonic NPCs from hundreds of clones, or lineage trees, containing tens of thousands of cells and millions of segmentations. New visualization approaches allow the NPC nuclear and process features to be effectively visualized for an entire clone. Significant differences in process and nuclear dynamics were found among type A and type C adult NPCs, and also between embryonic NPCs cultured from the anterior and posterior cerebral cortex.

Automated Measurement of Cobblestone Morphology for Characterizing Stem Cell Derived Retinal Pigment Epithelial Cell Cultures

PURPOSE

Assessing the morphologic properties of cells in microscopy images is an important task to evaluate cell health, identity, and purity. Typically, subjective visual assessments are accomplished by an experienced researcher. This subjective human step makes transfer of the evaluation process from the laboratory to the cell manufacturing facility difficult and time consuming.

METHODS

Automated image analysis can provide rapid, objective measurements of cultured cells, greatly aiding manufacturing, regulatory, and research goals. Automated algorithms for classifying images based on appearance characteristics typically either extract features from the image and use those features for classification or use the images directly as input to the classification algorithm. In this study we have developed both feature and nonfeature extraction methods for automatically measuring "cobblestone" structure in human retinal pigment epithelial (RPE) cell cultures.

RESULTS

A new approach using image compression combined with a Kolmogorov complexity-based distance metric enables robust classification of microscopy images of RPE cell cultures. The automated measurements corroborate determinations made by experienced cell biologists. We have also developed an approach for using steerable wavelet filters for extracting features to characterize the individual cellular junctions.

CONCLUSIONS

Two image analysis techniques enable robust and accurate characterization of the cobblestone morphology that is indicative of viable RPE cultures for therapeutic applications.

Live Imaging of Adult Neural Stem Cells in Rodents

The generation of cells of the neural lineage within the brain is not restricted to early development. New neurons, oligodendrocytes, and astrocytes are produced in the adult brain throughout the entire murine life. However, despite the extensive research performed in the field of adult neurogenesis during the past years, fundamental questions regarding the cell biology of adult neural stem cells (aNSCs) remain to be uncovered. For instance, it is crucial to elucidate whether a single aNSC is capable of differentiating into all three different macroglial cell types in vivo or these distinct progenies constitute entirely separate lineages. Similarly, the cell cycle length, the time and mode of division (symmetric vs. asymmetric) that these cells undergo within their lineage progression are interesting questions under current investigation. In this sense, live imaging constitutes a valuable ally in the search of reliable answers to the previous questions. In spite of the current limitations of technology new approaches are being developed and outstanding amount of knowledge is being piled up providing interesting insights in the behavior of aNSCs. Here, we will review the state of the art of live imaging as well as the alternative models that currently offer new answers to critical questions.

Single step synthesized sulfur and nitrogen doped carbon nanodots from whey protein: nanoprobes for longterm cell tracking crossing the barrier of photo-toxicity

Long-term cell tracking via whey protein derived carbon nanodots.

Adult stem cell lineage tracing and deep tissue imaging

Lineage tracing is a widely used method for understanding cellular dynamics in multicellular organisms during processes such as development, adult tissue maintenance, injury repair and tumorigenesis. Advances in tracing or tracking methods, from light microscopy-based live cell tracking to fluorescent label-tracing with two-photon microscopy, together with emerging tissue clearing strategies and intravital imaging approaches have enabled scientists to decipher adult stem and progenitor cell properties in various tissues and in a wide variety of biological processes. Although technical advances have enabled time-controlled genetic labeling and simultaneous live imaging, a number of obstacles still need to be overcome. In this review, we aim to provide an in-depth description of the traditional use of lineage tracing as well as current strategies and upcoming new methods of labeling and imaging.

Survey statistics of automated segmentations applied to optical imaging of mammalian cells

BACKGROUND

The goal of this survey paper is to overview cellular measurements using optical microscopy imaging followed by automated image segmentation. The cellular measurements of primary interest are taken from mammalian cells and their components. They are denoted as two- or three-dimensional (2D or 3D) image objects of biological interest. In our applications, such cellular measurements are important for understanding cell phenomena, such as cell counts, cell-scaffold interactions, cell colony growth rates, or cell pluripotency stability, as well as for establishing quality metrics for stem cell therapies. In this context, this survey paper is focused on automated segmentation as a software-based measurement leading to quantitative cellular measurements.

METHODS

We define the scope of this survey and a classification schema first. Next, all found and manually filteredpublications are classified according to the main categories: (1) objects of interests (or objects to be segmented), (2) imaging modalities, (3) digital data axes, (4) segmentation algorithms, (5) segmentation evaluations, (6) computational hardware platforms used for segmentation acceleration, and (7) object (cellular) measurements. Finally, all classified papers are converted programmatically into a set of hyperlinked web pages with occurrence and co-occurrence statistics of assigned categories.

RESULTS

The survey paper presents to a reader: (a) the state-of-the-art overview of published papers about automated segmentation applied to optical microscopy imaging of mammalian cells, (b) a classification of segmentation aspects in the context of cell optical imaging, (c) histogram and co-occurrence summary statistics about cellular measurements, segmentations, segmented objects, segmentation evaluations, and the use of computational platforms for accelerating segmentation execution, and (d) open research problems to pursue.

CONCLUSIONS

The novel contributions of this survey paper are: (1) a new type of classification of cellular measurements and automated segmentation, (2) statistics about the published literature, and (3) a web hyperlinked interface to classification statistics of the surveyed papers at https://isg.nist.gov/deepzoomweb/resources/survey/index.html.

Inference of spatiotemporal effects on cellular state transitions from time-lapse microscopy

BACKGROUND

Time-lapse microscopy allows to monitor cell state transitions in a spatiotemporal context. Combined with single cell tracking and appropriate cell state markers, transition events can be observed within the genealogical relationship of a proliferating population. However, to infer the correlations between the spatiotemporal context and cell state transitions, statistical analysis with an appropriately large number of samples is required.

RESULTS

Here, we present a method to infer spatiotemporal features predictive of the state transition events observed in time-lapse microscopy data. We first formulate a generative model, simulate different scenarios, such as time-dependent or local cell density-dependent transitions, and illustrate how to estimate univariate transition rates. Second, we formulate the problem in a machine-learning language using regularized linear models. This allows for a multivariate analysis and to disentangle indirect dependencies via feature selection. We find that our method can accurately recover the relevant features and reconstruct the underlying interaction kernels if a critical number of samples is available. Finally, we explicitly use the tree structure of the data to validate if the estimated model is sufficient to explain correlated transition events of sister cells.

CONCLUSIONS

Using synthetic cellular genealogies, we prove that our method is able to correctly identify features predictive of state transitions and we moreover validate the chosen model. Our approach allows to estimate the number of cellular genealogies required for the proposed spatiotemporal statistical analysis, and we thus provide an important tool for the experimental design of challenging single cell time-lapse microscopy assays.

Computational Image Analysis Reveals Intrinsic Multigenerational Differences between Anterior and Posterior Cerebral Cortex Neural Progenitor Cells

Time-lapse microscopy can capture patterns of development through multiple divisions for an entire clone of proliferating cells. Images are taken every few minutes over many days, generating data too vast to process completely by hand. Computational analysis of this data can benefit from occasional human guidance. Here we combine improved automated algorithms with minimized human validation to produce fully corrected segmentation, tracking, and lineaging results with dramatic reduction in effort. A web-based viewer provides access to data and results. The improved approach allows efficient analysis of large numbers of clones. Using this method, we studied populations of progenitor cells derived from the anterior and posterior embryonic mouse cerebral cortex, each growing in a standardized culture environment. Progenitors from the anterior cortex were smaller, less motile, and produced smaller clones compared to those from the posterior cortex, demonstrating cell-intrinsic differences that may contribute to the areal organization of the cerebral cortex.

•

Open-source automated software designed to track stem/progenitor clones in time-lapse movies

•

Software tools for easy data validation and visualization greatly improve efficiency

•

Lineage tree reconstruction from hundreds of embryonic mouse forebrain clones

•

Intrinsic differences in progenitor behavior from anterior/posterior cerebral cortex

Normalized Compression Distance of Multisets with Applications

Pairwise normalized compression distance (NCD) is a parameter-free, feature-free, alignment-free, similarity metric based on compression. We propose an NCD of multisets that is also metric. Previously, attempts to obtain such an NCD failed. For classification purposes it is superior to the pairwise NCD in accuracy and implementation complexity. We cover the entire trajectory from theoretical underpinning to feasible practice. It is applied to biological (stem cell, organelle transport) and OCR classification questions that were earlier treated with the pairwise NCD. With the new method we achieved significantly better results. The theoretic foundation is Kolmogorov complexity.

Extracting meaning from biological imaging data

Biological imaging continues to improve, capturing continually longer-term, richer, and more complex data, penetrating deeper into live tissue. How do we gain insight into the dynamic processes of disease and development from terabytes of multidimensional image data? Here I describe a collaborative approach to extracting meaning from biological imaging data. The collaboration consists of teams of biologists and engineers working together. Custom computational tools are built to best exploit application-specific knowledge in order to visualize and analyze large and complex data sets. The image data are summarized, extracting and modeling the features that capture the objects and relationships in the data. The summarization is validated, the results visualized, and errors corrected as needed. Finally, the customized analysis and visualization tools together with the image data and the summarization results are shared. This Perspective provides a brief guide to the mathematical ideas that rigorously quantify the notion of extracting meaning from biological image, and to the practical approaches that have been used to apply these ideas to a wide range of applications in cell and tissue optical imaging.

Global linking of cell tracks using the Viterbi algorithm

Automated tracking of living cells in microscopy image sequences is an important and challenging problem. With this application in mind, we propose a global track linking algorithm, which links cell outlines generated by a segmentation algorithm into tracks. The algorithm adds tracks to the image sequence one at a time, in a way which uses information from the complete image sequence in every linking decision. This is achieved by finding the tracks which give the largest possible increases to a probabilistically motivated scoring function, using the Viterbi algorithm. We also present a novel way to alter previously created tracks when new tracks are created, thus mitigating the effects of error propagation. The algorithm can handle mitosis, apoptosis, and migration in and out of the imaged area, and can also deal with false positives, missed detections, and clusters of jointly segmented cells. The algorithm performance is demonstrated on two challenging datasets acquired using bright-field microscopy, but in principle, the algorithm can be used with any cell type and any imaging technique, presuming there is a suitable segmentation algorithm.

Quantitative single-cell approaches to stem cell research

Understanding the molecular control of cell fates is central to stem cell research. Such insight requires quantification of molecular and cellular behavior at the single-cell level. Recent advances now permit high-throughput molecular readouts from single cells as well as continuous, noninvasive observation of cell behavior over time. Here, we review current state-of-the-art approaches used to query stem cell fate at the single-cell level, including advances in lineage tracing, time-lapse imaging, and molecular profiling. We also offer our perspective on the advantages and drawbacks of available approaches, key technical limitations, considerations for data interpretation, and future innovation.

Visualization and correction of automated segmentation, tracking and lineaging from 5-D stem cell image sequences

BACKGROUND

Neural stem cells are motile and proliferative cells that undergo mitosis, dividing to produce daughter cells and ultimately generating differentiated neurons and glia. Understanding the mechanisms controlling neural stem cell proliferation and differentiation will play a key role in the emerging fields of regenerative medicine and cancer therapeutics. Stem cell studies in vitro from 2-D image data are well established. Visualizing and analyzing large three dimensional images of intact tissue is a challenging task. It becomes more difficult as the dimensionality of the image data increases to include time and additional fluorescence channels. There is a pressing need for 5-D image analysis and visualization tools to study cellular dynamics in the intact niche and to quantify the role that environmental factors play in determining cell fate.

RESULTS

We present an application that integrates visualization and quantitative analysis of 5-D (x,y,z,t,channel) and large montage confocal fluorescence microscopy images. The image sequences show stem cells together with blood vessels, enabling quantification of the dynamic behaviors of stem cells in relation to their vascular niche, with applications in developmental and cancer biology. Our application automatically segments, tracks, and lineages the image sequence data and then allows the user to view and edit the results of automated algorithms in a stereoscopic 3-D window while simultaneously viewing the stem cell lineage tree in a 2-D window. Using the GPU to store and render the image sequence data enables a hybrid computational approach. An inference-based approach utilizing user-provided edits to automatically correct related mistakes executes interactively on the system CPU while the GPU handles 3-D visualization tasks.

CONCLUSIONS

By exploiting commodity computer gaming hardware, we have developed an application that can be run in the laboratory to facilitate rapid iteration through biological experiments. We combine unsupervised image analysis algorithms with an interactive visualization of the results. Our validation interface allows for each data set to be corrected to 100% accuracy, ensuring that downstream data analysis is accurate and verifiable. Our tool is the first to combine all of these aspects, leveraging the synergies obtained by utilizing validation information from stereo visualization to improve the low level image processing tasks.

Induction of the ganglion cell differentiation program in human retinal progenitors before cell cycle exit

BACKGROUND

Despite the disease relevance, understanding of human retinal development lags behind that of other species. We compared the kinetics of gene silencing or induction during ganglion cell development in human and murine retina.

RESULTS

Induction of POU4F2 (BRN3B) marks ganglion cell commitment, and we detected this factor in S-phase progenitors that had already silenced Cyclin D1 and VSX2 (CHX10). This feature was conserved in human and mouse retina, and the fraction of Pou4f2+ murine progenitors labeled with a 30 min pulse of BrdU matched the fraction of ganglion cells predicted to be born in a half-hour period. Additional analysis of 18 markers revealed many with conserved kinetics, such as the POU4F2 pattern above, as well as the surprising maintenance of "cell cycle" proteins KI67, PCNA, and MCM6 well after terminal mitosis. However, four proteins (TUBB3, MTAP1B, UCHL1, and RBFOX3) showed considerably delayed induction in human relative to mouse retina, and two proteins (ISL1, CALB2) showed opposite kinetics, appearing on either side of terminal mitosis depending on the species.

CONCLUSION

With some notable exceptions, human and murine ganglion cell differentiation show similar kinetics, and the data add weight to prior studies supporting the existence of biased ganglion cell progenitors.

Labour-efficient in vitro lymphocyte population tracking and fate prediction using automation and manual review

Interest in cell heterogeneity and differentiation has recently led to increased use of time-lapse microscopy. Previous studies have shown that cell fate may be determined well in advance of the event. We used a mixture of automation and manual review of time-lapse live cell imaging to track the positions, contours, divisions, deaths and lineage of 44 B-lymphocyte founders and their 631 progeny in vitro over a period of 108 hours. Using this data to train a Support Vector Machine classifier, we were retrospectively able to predict the fates of individual lymphocytes with more than 90% accuracy, using only time-lapse imaging captured prior to mitosis or death of 90% of all cells. The motivation for this paper is to explore the impact of labour-efficient assistive software tools that allow larger and more ambitious live-cell time-lapse microscopy studies. After training on this data, we show that machine learning methods can be used for realtime prediction of individual cell fates. These techniques could lead to realtime cell culture segregation for purposes such as phenotype screening. We were able to produce a large volume of data with less effort than previously reported, due to the image processing, computer vision, tracking and human-computer interaction tools used. We describe the workflow of the software-assisted experiments and the graphical interfaces that were needed. To validate our results we used our methods to reproduce a variety of published data about lymphocyte populations and behaviour. We also make all our data publicly available, including a large quantity of lymphocyte spatio-temporal dynamics and related lineage information.

Segmentation of occluded hematopoietic stem cells from tracking

Image sequences of live proliferating cells often contain visual ambiguities that are difficult even for human domain experts to resolve. Here we present a new approach to analyzing image sequences that capture the development of clones of hematopoietic stem cells (HSCs) from live cell time lapse microscopy. The HSCs cannot survive long term imaging unless they are cultured together with a secondary cell type, OP9 stromal cells. The HSCs frequently disappear under the OP9 cell layer, making segmentation difficult or impossible from a single image frame, even for a human domain expert. We have developed a new approach to the segmentation of HSCs that captures these occluded cells. Starting with an a priori segmentation that uses a Monte Carlo technique to estimate the number of cells in a clump of touching cells, we proceed to track and lineage the image data. Following user validation of the lineage information, an a posteriori resegmentation step utilizing tracking results delineates the HSCs occluded by the OP9 layer. Resegmentation has been applied to 3031 occluded segmentations from 77 tracks, correctly recovering over 84% of the occluded segmentations.

Making of a retinal cell: insights into retinal cell-fate determination

Understanding the process by which an uncommitted dividing cell produces particular specialized cells within a tissue remains a fundamental question in developmental biology. Many tissues are well suited for cell-fate studies, but perhaps none more so than the developing retina. Traditionally, experiments using the retina have been designed to elucidate the influence that individual environmental signals or transcription factors can have on cell-fate decisions. Despite a substantial amount of information gained through these studies, there is still much that we do not yet understand about how cell fate is controlled on a systems level. In addition, new factors such as noncoding RNAs and regulators of chromatin have been shown to play roles in cell-fate determination and with the advent of "omics" technology more factors will most likely be identified. In this chapter we summarize both the traditional view of retinal cell-fate determination and introduce some new ideas that are providing a challenge to the older way of thinking about the acquisition of cell fates.

An analysis-synthesis approach for neurosphere modelisation under phase-contrast microscopy

The study of stem cells is one of the most important biomedical research. Understanding their development could allow multiple applications in regenerative medicine. For this purpose, automated solutions for the observation of stem cell development process are needed. This study introduces an on-line analysis method for the modelling of neurosphere evolution during the early time of their development under phase contrast microscopy. From the corresponding phase contrast time-lapse sequences, we extract information from the neurosphere using a combination of phase contrast physics deconvolution and curve detection for locate the cells inside the neurosphere. Then, based on prior biological knowledge, we generate possible and optimal 3-dimensional configuration using 2D to 3D registration methods and evolutionary optimisation algorithm.

Personalized prediction of EGFR mutation-induced drug resistance in lung cancer

EGFR mutation-induced drug resistance has significantly impaired the potency of small molecule tyrosine kinase inhibitors in lung cancer treatment. Computational approaches can provide powerful and efficient techniques in the investigation of drug resistance. In our work, the EGFR mutation feature is characterized by the energy components of binding free energy (concerning the mutant-inhibitor complex), and we combine it with specific personal features for 168 clinical subjects to construct a personalized drug resistance prediction model. The 3D structure of an EGFR mutant is computationally predicted from its protein sequence, after which the dynamics of the bound mutant-inhibitor complex is simulated via AMBER and the binding free energy of the complex is calculated based on the dynamics. The utilization of extreme learning machines and leave-one-out cross-validation promises a successful identification of resistant subjects with high accuracy. Overall, our study demonstrates advantages in the development of personalized medicine/therapy design and innovative drug discovery.

An automatic method for robust and fast cell detection in bright field images from high-throughput microscopy

Background

In recent years, high-throughput microscopy has emerged as a powerful tool to analyze cellular dynamics in an unprecedentedly high resolved manner. The amount of data that is generated, for example in long-term time-lapse microscopy experiments, requires automated methods for processing and analysis. Available software frameworks are well suited for high-throughput processing of fluorescence images, but they often do not perform well on bright field image data that varies considerably between laboratories, setups, and even single experiments.

Results

In this contribution, we present a fully automated image processing pipeline that is able to robustly segment and analyze cells with ellipsoid morphology from bright field microscopy in a high-throughput, yet time efficient manner. The pipeline comprises two steps: (i) Image acquisition is adjusted to obtain optimal bright field image quality for automatic processing. (ii) A concatenation of fast performing image processing algorithms robustly identifies single cells in each image. We applied the method to a time-lapse movie consisting of ∼315,000 images of differentiating hematopoietic stem cells over 6 days. We evaluated the accuracy of our method by comparing the number of identified cells with manual counts. Our method is able to segment images with varying cell density and different cell types without parameter adjustment and clearly outperforms a standard approach. By computing population doubling times, we were able to identify three growth phases in the stem cell population throughout the whole movie, and validated our result with cell cycle times from single cell tracking.

Conclusions

Our method allows fully automated processing and analysis of high-throughput bright field microscopy data. The robustness of cell detection and fast computation time will support the analysis of high-content screening experiments, on-line analysis of time-lapse experiments as well as development of methods to automatically track single-cell genealogies.

Synergistic combination of clinical and imaging features predicts abnormal imaging patterns of pulmonary infections

We designed and tested a novel hybrid statistical model that accepts radiologic image features and clinical variables, and integrates this information in order to automatically predict abnormalities in chest computed-tomography (CT) scans and identify potentially important infectious disease biomarkers. In 200 patients, 160 with various pulmonary infections and 40 healthy controls, we extracted 34 clinical variables from laboratory tests and 25 textural features from CT images. From the CT scans, pleural effusion (PE), linear opacity (or thickening) (LT), tree-in-bud (TIB), pulmonary nodules, ground glass opacity (GGO), and consolidation abnormality patterns were analyzed and predicted through clinical, textural (imaging), or combined attributes. The presence and severity of each abnormality pattern was validated by visual analysis of the CT scans. The proposed biomarker identification system included two important steps: (i) a coarse identification of an abnormal imaging pattern by adaptively selected features (AmRMR), and (ii) a fine selection of the most important features from the previous step, and assigning them as biomarkers, depending on the prediction accuracy. Selected biomarkers were used to classify normal and abnormal patterns by using a boosted decision tree (BDT) classifier. For all abnormal imaging patterns, an average prediction accuracy of 76.15% was obtained. Experimental results demonstrated that our proposed biomarker identification approach is promising and may advance the data processing in clinical pulmonary infection research and diagnostic techniques.

Live imaging and analysis of postnatal mouse retinal development

Background

The explanted, developing rodent retina provides an efficient and accessible preparation for use in gene transfer and pharmacological experimentation. Many of the features of normal development are retained in the explanted retina, including retinal progenitor cell proliferation, heterochronic cell production, interkinetic nuclear migration, and connectivity. To date, live imaging in the developing retina has been reported in non-mammalian and mammalian whole-mount samples. An integrated approach to rodent retinal culture/transfection, live imaging, cell tracking, and analysis in structurally intact explants greatly improves our ability to assess the kinetics of cell production.

Results

In this report, we describe the assembly and maintenance of an in vitro, CO₂-independent, live mouse retinal preparation that is accessible by both upright and inverted, 2-photon or confocal microscopes. The optics of this preparation permit high-quality and multi-channel imaging of retinal cells expressing fluorescent reporters for up to 48h. Tracking of interkinetic nuclear migration within individual cells, and changes in retinal progenitor cell morphology are described. Follow-up, hierarchical cluster screening revealed that several different dependent variable measures can be used to identify and group movement kinetics in experimental and control samples.

Conclusions

Collectively, these methods provide a robust approach to assay multiple features of rodent retinal development using live imaging.

Deciphering the single-cell omic: innovative application for translational medicine

Traditional technologies to investigate system biology are limited by the detection of parameters resulting from the averages of large populations of cells, missing cells produced in small numbers, and attempting to uniform the heterogeneity. The advent of proteomics and genomics at a single-cell level has set the basis for an outstanding improvement in analytical technology and data acquisition. It has been well demonstrated that cellular heterogeneity is closely related to numerous stochastic transcriptional events leading to variations in patterns of expression among single genetically identical cells. The new-generation technology of single-cell analysis is able to better characterize a cell's population, identifying and differentiating outlier cells, in order to provide both a single-cell experiment and a corresponding bulk measurement, through the identification, quantification and characterization of all system biology aspects (genomics, transcriptomics, proteomics, metabolomics, degradomics and fluxomics). The movement of omics into single-cell analysis represents a significant and outstanding shift.

Neural progenitor cells labeling with microbubble contrast agent for ultrasound imaging in vivo

Tracking neuroprogenitor cells (NPCs) that are used to target tumors, infarction or inflammation, is paramount for cell-based therapy. We employed ultrasound imaging that can detect a single microbubble because it can distinguish its unique signal from those of surrounding tissues. NPCs efficiently internalized positively charged microbubbles allowing a clinical ultrasound system to detect a single cell at 7 MHz. When injected intravenously, labeled NPCs traversed the lungs to be imaged in the left ventricle and the liver where they accumulated. Internalized microbubbles were not only less sensitive to destruction by ultrasound, but remained visible in vivo for days as compared to minutes when given free. The extended longevity provides ample time to allow cells to reach their intended target. We were also able to transfect NPCs in vitro when microbubbles were preloaded with GFP plasmid only when cells were insonated. Transfection efficiency and cell viability were both greater than 90%.