Abstract
In vivo flow cytometry provides a non-invasive way of probing the biology of circulating cells during disease progression and studying cellular response to therapy. However, current methods provide little morphological information which potentially could be new biological marker for early disease diagnosis, and fail to reveal intercellular interactions. Here we report a multi-color, multiphoton in vivo imaging flow cytometry, to image circulating cells within the vasculature of scattering tissues at high spatiotemporal resolution. We apply it in imaging of cellular dynamics in bone marrow through the intact mouse skull, in situ deformability cytometry, distinguishing cellular clusters, and simultaneously monitoring multiple types of trafficking cells based on their morphologies and fluorescence emission colors.
© 2016 Optical Society of America
1. Introduction
Traditional flow cytometry (FC) is used routinely to acquire quantitative information about specific cell populations at high throughput, and has found extensive clinical applications in diagnostic pathology [1, 2]. However, blood extraction in conventional FC is invasive and the extensive processing required by in vitro detection can introduce potential artifacts [3]. Moreover, blood sampling reduces the sensitivity of capturing rare cells (such as circulating tumor cells) and prevents the long-term study of cell metabolism in their native complex biological environment [4]. In fact, detection and quantification of rare circulating cells in vivo are very important for early diagnosis of diseases (such as cancer, stroke, and inflammation), and long-term study is necessary for probing the biology of circulating cells during disease progression, and for studying cellular response to therapy (such as drugs and radiation) [1, 3–9]. To this end, several methods of in vivo FC have been developed to continuously monitor circulating cells in live animals without affecting the physiology of the subject, and found applications in many biomedical studies, such as cancer, immunology, and stem cells [3–5,10–12]. A drawback of in vivo FC is that it generally provides little morphological information (such as the size, shape, morphology, and deformability of circulating cells), which can potentially be new biological marker that is sensitive to early disease development [13, 14].
The recent advent of in vivo imaging FC [8, 15], or flow cytography [14], enables researchers to record morphologies and intercellular interactions of trafficking cells in situ [1, 16]. Wide-field imaging based FC has been demonstrated to distinguish potential cellular clusters [15], and to assess cellular morphology and deformability in the native biological environment [8]. To get depth-resolved images, a three-dimensional (3D) image cytometer based on widefield structured light microscopy and high-speed remote depth scanning was developed, with 3D resolution to improve quantification of cytometric features [13]. However, wide-field imaging based methods are highly susceptible to optical scattering in biological tissue, and are limited to the superficial surface or transparent samples [17]. For deep tissue imaging, multi-photon microscopies, such as two-photon fluorescence microscopy (TPM) [18], remain the methods of choice. In this regard, laser scanning (point scanning) multiphoton imaging FC is ideal [17, 19] due to the intrinsic advantages in the 3D high resolution sectioning, the improved signal from both ballistic and scattered fluorescence emission, and the deep penetration depth of the near infrared laser excitation. Unfortunately, the current imaging speed of TPM is generally too slow to capture in detail the morphology of fast flowing cells [20]. This is primarily a consequence of limited scanning speed, particullary in the axial direction. Generally it is necessary to adopt 1D line scanning, instead of 2D cross-sectional scanning, to ensure adequately fast sampling of vascular dynamics (Fig. 1), which however loses the spatial information in the axial dimension [10]. Meanwhile, methods based on optical coherence tomography [5] and photoacoustics [14] have been developed for in vivo imaging FC.
To improve the axial scanning speed of TPM, several methods have been proposed, such as remote focusing microscopy [21, 22] and random access microscopy [23, 24]. But none of them has shown the capability of high-speed cross-sectional imaging for imaging FC. Recently, we have integrated optical phase-locked ultrasound lens (OPLUL) with a standard TPM for microsecond scale axial scanning [25]. Here we improve the cross-sectional frame rate to 2 kHz, and demonstrate the first multicolor multiphoton in vivo imaging FC, to our best knowledge. We applied it in imaging cellular dynamics in bone marrow through the intact mouse skull, and also showed in situ deformability cytometry [Fig. 1(b)]. With the improved cross-sectional frame rate, we were able to distinguish clusters and resolve morphologies of rapidly flowing cells. We also demonstrated the capability to concurrently monitor multiple types of circulation cells based on fluorescence emission wavelengths.
2. Experimental setup
The multicolor multiphoton in vivo imaging FC system is based on the the intergration of the OPLUL into a standard TPM, as shown in Fig. 2(a). The two-photon excitation femtosecond (fs) beam (80 MHz, 935 nm, 140 fs, Chameleon, Coherent) and the reference continuous-wave (CW) beam at 785 nm, both horizontally polarized, were combined with the dichroic beam splitter (DBS1). After going through the polarization beam splitter (PBS), quarter waveplate (QWP) and ultrasound lenses (UL, TAG lens 2.0, TAG Optics), the beams were imaged onto a mirror (M1) by a pair of relay lens (RL1 and RL1’). The reflected beams from M1 went through the UL and QWP again, and were reflected by the PBS as their polarizations were rotated by 90° after the second pass of the QWP. The second pair of relay lenses (RL2 and RL2’) were used to image the beam onto the x-y galvo scanner. The beams were separated by DBS2, with the reflected CW beam entering the photodiode (PD) and the transmitted fs beam employed for multiphoton excitation. The CW beam was spatially filtered by an iris before entering the PD, whose signal was used as the feedback for the phase-locked loop (PLL, HF2LI-PLL, Zurich Instruments) to maintain stable phase and amplitude of UL resonance. These elements formed the OPLUL, which combined with a conventional TPM enabled volumetric imaging and FC.
We conjugated the oscillating wavefronts [25, 26] generated by the OPLUL to the pupil plane of the objective (Nikon 16 × , 0.8 NA) for high-speed axial scanning [25]. The lateral and axial resolutions were ~0.5-0.6 µm and 2.9-4.6 µm, respectively. We drove the UL at the 450 kHz resonance, and achieved microsecond scale axial scanning. With a galvo scanner for transversal scanning, we obtained the kHz scale cross-sectional frame rate maintaining high spatial resolution. As an example, in Fig. 2(b), we show the raw data acquired by cross-sectional scanning a 1 µm fluorescent bead, in which the reference signal from the PD was used to reconstruct the cross-sectional image [Fig. 2(c)]. We used a high-capacity field-programmable gate array to record, process and display the data stream in real time.
3. Experimental results
We flexibly operated the system in both volumetric imaging (x, y, z, t) and cross-sectional imaging (x, z, t) schemes for FC applications. All procedures involving mice were approved by the Animal Care and Use Committees of Purdue University, HHMI Janelia Research Campus and National Institute of Allergy and Infectious Diseases, National Institutes of Health.
3.1 In vivo volumetric imaging of cellular dynamics in the bone marrow through intact mouse skull
To show the advantage of our system in imaging through scattering tissues, we demonstrated the imaging of cellular dynamics in bone marrow through the intact mouse skull. The bone marrow is an important site, where all blood cells are formed and hematologic malignancies such as leukemia and metastasis of solid tumors emerge [27]. Noninvasive methods to image the cellular dynamics in the bone marrow are particularly important for stem cell studies, where any perturbation would disturb the fundamental biology [27, 28].
We used the CX3CR1gfp/gfp transgenic mice with EGFP-expressing monocytes, dendritic cells, natural killer cells, and brain microglia [29]. After removing the scalp, we directly imaged the dynamics of EGFP-expressing cells in bone marrow through the intact skull. The thickness of the skull above the bone marrow cavity was ~50 µm based on the second harmonic generation signal of the bone matrix [Fig. 3(a), and Visualization 1]. We recorded cellular trafficking and migrations inside the bone marrow [Figs. 3(b) and 3(c), and Visualization 2] at 7 Hz (210 × 105 × 40 µm3, 50-90 µm under surface of the skull, laser power: 70 mW). In Fig. 3(b), we show a cell trafficking through the vascular while simultaneously interacting with the microenvironment (from time point t1 to t3). In Fig. 3(c), we show the traces of a migrating cell. This is the first report on in vivo fast volumetric (4D) imaging of cellular dynamics in the bone marrow without thinning the skull.
3.2 In situ deformability cytometry and imaging of intercellular interactions
We have also applied our system as an in situ deformability cytometer for mechanical phenotyping. Cell deformability indicates changes in the cytoskeleton and nuclear organization, which may provide a label-free biomarker for determining cell states or properties [8, 30, 31]. Current techniques are mostly based on microfluidic devices for in vitro quantification [30, 31], which however may bring artifacts during the invasive sample isolations. In the native biological environment, circulating cells can undergo apparent deformations in geometrically irregular regions of the vasculature, such as bifurcations, curves, and narrowed areas [8]. Taking advantage of our continuous volumetric imaging system, we were able to quantify the deformability of circulating cells in situ. We used the Lyz2gfp/+ B6.Albino mice with GFP-expressing neutrophils [32]. After craniotomy and the installation of optical windows, we stained the brain tissue with sulforhodamine 101 (SR 101) for 2 minutes. We induced a gentle laser ablation of the neocortex to activate neutrophils and to recruit neutrophils from the bone marrow to the damage site through circulation system. We imaged the vasculature around the damage site (151 × 28 × 40 µm3, 10-50 µm under dura, 63 mW), and took snapshots of neutrophils when they passed through a bifurcation of blood vessels at 28 Hz (Fig. 4, and Visualization 3). From the volumetric imaging data, we were able to quantify the flowing speeds, which were 0.79 ± 0.42 (mean ± standard deviation) and 0.38 ± 0.15 mm/s in the main branch and one side branch, respectively. The velocity difference may indicate the physical pressure change that would result in deformation of the neutrophils passing by [Fig. 4(a)]. At the same time, some neutrophils were still able to crawl along the wall of the blood vessels at the bifurcation [Fig. 4(b)]. The capability of quantifying cell deformability and cell-vessel wall interactions may benefit early disease diagnosis and drug screening.
3.3 In vivo imaging FC for distinguishing cell clusters at high spatiotemporal resolution
Using a galvo scanner for transversal scanning, we achieved a cross-sectional (x-z) frame rate up to 2 kHz, which enabled us to distinguish cell clusters at high spatiotemporal resolutions. It was reported that although rare in circulation, clusters of the circulating tumor cells have much higher metastatic potential compared with single ones [33]. Quantifying the clusters of circulating cells and tracking their fates are critical in studying disease progression [15]. We mimicked this process by monitoring activated neutrophils after tissue damage. We used the DsRed+/− Lyz2gfp/+ B6.Albino mice with GFP-expressing neutrophils. After inducing mild tissue damage in the ear pinnae, we imaged the trafficking of GFP-expressing neutrophils in the vasculature of the mouse ear. We took both volumetric imaging and cross-sectional imaging, along and across the blood vessel, respectively. With volumetric imaging (31 Hz, 168 × 28 × 40 µm3, 70-110 µm under surface, 72 mW), we were able to capture the clusters [Fig. 5(a), and Visualization 4] and calculate the flowing speed, while cross-sectional imaging (2 kHz, x × z: 28 × 40 µm2, 70-110 µm under surface, 72 mW) provided cellular images with smaller distortions [inset in Fig. 5(b), and Visualization 5]. In the ~40 minutes of continuous cross-sectional imaging, no tissue damage was observed, and only ~10 out of 764 neutrophils were potential clusters. We were able to statistically analyze the temporal separation between consecutively passing neutrophils [Fig. 5(b)], and capture the morphology details of the clusters flowing at 0.82 ± 0.16 mm/s.
3.4 In vivo imaging FC for distinguishing cells based on morphologies and colors
Simultaneously quantifying multiple populations of circulating cells is important for disease diagnoses [34, 35]. For example, the ratio of neutrophil to lymphocyte is found to be a predictor of stroke [36]. Taking advantages of the high spatial resolution intrinsinc to multiphoton point-scanning and the 2 kHz cross-sectional frame rate of our system, we were able to distinguish different cell types based on cellular morphologies. We collected mito-Dendra2 expressing lymphocytes, whose mitochondria are labeled with green-emitting dendra2, from the spleen of PhAMexcised transgenic mice, and transferred ~100 µL of the diluted cellular solution into the circulation system of Lyz2gfp/+ B6.Albino mice via tail-vein injections. After inducing mild tissue damage in the ear pinnae as described above, we imaged trafficking cells (flowing at 1.35 ± 0.16 mm/s) in a blood vessel of mouse ear with 22.4 × 40 µm2 cross-secional area at 2 kHz (60-100 µm under surface, 72 mW). In Figs. 6(a)-6(c) and Visualization 6,Visualization 7, and Visualization 8, we show the morphologies of the Dendra2-expressing lymphocytes and the GFP-expressing neutrophils. Although they were of the same emission color, we were able to distinguish these two cell types from their morphologies, and quantify their numbers and temporal distributions [Fig. 6(d)].
Moreover, based on the broad multiphoton excitation spectra, we were able to quantify different types of circulating cells with the color-labeling strategy [34]. We injected ~100 µL SR 101 (10 mg/ml concentration in saline, 20 mg/kg) into the Lyz2gfp/+ B6.Albino mice via tail-veins. After ~1 hour, the SR 101 in the plasma was cleared, and certain types of blood cells were stained. We introduced tissue damage in the mouse neocortex, and imaged cell trafficking (at 0.36 ± 0.03 mm/s) in the brain vasculature through cranial windows. We were able to capture the morphologies of the fast flowing SR101 stained cells and the GFP-expressing neutrophils at 2 kHz cross-sectional frame rate (42 × 40 µm2, 10-50 µm under dura, 44 mW), even when they were close in time [Figs. 7(a) and 7(b), and Visualization 9]. We show their temporal distribution in Fig. 7(c).
4. Discussion
Different from the axial scanning methods based on translating objectives in which the vibration of the objective may transmit to the biological specimens, our method achieved microsecond scale axial scanning without introducing such perturbations. With the 2 kHz cross-sectional frame-rate, our system can flexibly perform volumetric imaging and cross-sectional imaging. Volumetric measurement enables deformablity cytometer, imaging intercellular interactions, and qualifying flow speed. Cross-sectional measurement enables imaging cytometry. The intergration with TPM allowed us to image vasculature in deep tissue at high resolution. Moreover, the broad multiphoton excitation spectra makes it possible to monitor multiple circulating cells of different emission colors.
In addition to TPM, our microsecond scale axial scanning method can also be combined with other multiphoton microscopies, such as coherent Raman microscopy [37] and stimulated emission microscopy [38], for lable-free imaging of circulation cells. This would be of great interest for both biological and clinical studies.
Intergrating adaptive optics with our system, we can compensate the focus distortions induced by the inhomegenous refractive index distribution and increase the imaging depth [39–44]. Using a polygon scanner for transversal scanning (x scan), we can potentially improve the cross-sectional frame rate to 30 kHz or even higher value [45].
5. Conclusions
In summary, our multicolor multiphoton in vivo imaging FC system has superior performance in both volumetric imaging for studying cellular deformability/intercellular interactions and quantifying flow speed, and in cross-sectional imaging for capturing precise morphology information of fast circulating cells and distinguishing clusters and cell types. We expect its broad applications in biology and pharmacology.
Acknowledgments
The authors thank Charles P. Lin at Harvard Medical School for the initial idea of blood flow imaging. L.K. and M.C. thank the support from NIH (1U01NS094341-01), Purdue University, and Howard Hughes Medical Institute. This research was also supported in part by the Intramural Research Program of the National Institute of Allergy and Infectious Diseases, National Institutes of Health.
References and links
1. V. V. Tuchin, Advanced Optical Flow Cytometry: Methods and Disease Diagnoses (Wiley-VCH, 2011).
2. S. M. Daly and M. J. Leahy, “‘Go with the flow ’: A review of methods and advancements in blood flow imaging,” J. Biophotonics 6(3), 217–255 (2013). [CrossRef] [PubMed]
3. V. V. Tuchin, A. Tárnok, and V. P. Zharov, “In vivo flow cytometry: A horizon of opportunities,” Cytometry A 79(10), 737–745 (2011). [CrossRef] [PubMed]
4. I. Georgakoudi, N. Solban, J. Novak, W. L. Rice, X. Wei, T. Hasan, and C. P. Lin, “In vivo flow cytometry: a new method for enumerating circulating cancer cells,” Cancer Res. 64(15), 5044–5047 (2004). [CrossRef] [PubMed]
5. S. P. Morgan, “Can new optical techniques for in vivo imaging and flow cytometry of the microcirculation benefit sickle cell disease research?” Cytometry A 79A(10), 766–774 (2011). [CrossRef] [PubMed]
6. Z.-C. Fan, J. Yan, G.-D. Liu, X.-Y. Tan, X.-F. Weng, W.-Z. Wu, J. Zhou, and X.-B. Wei, “Real-time monitoring of rare circulating hepatocellular carcinoma cells in an orthotopic model by in vivo flow cytometry assesses resection on metastasis,” Cancer Res. 72(10), 2683–2691 (2012). [CrossRef] [PubMed]
7. Y. Zeng, B. Yan, Q. Sun, S. K. Teh, W. Zhang, Z. Wen, and J. Y. Qu, “Label-free in vivo imaging of human leukocytes using two-photon excited endogenous fluorescence,” J. Biomed. Opt. 18(4), 040504 (2013). [CrossRef] [PubMed]
8. V. P. Zharov, E. I. Galanzha, Y. Menyaev, and V. V. Tuchin, “In vivo high-speed imaging of individual cells in fast blood flow,” J. Biomed. Opt. 11(5), 054034 (2006). [CrossRef] [PubMed]
9. D. A. Nedosekin, V. V. Verkhusha, A. V. Melerzanov, V. P. Zharov, and E. I. Galanzha, “In vivo photoswitchable flow cytometry for direct tracking of single circulating tumor cells,” Chem. Biol. 21(6), 792–801 (2014). [CrossRef] [PubMed]
10. W. He, H. Wang, L. C. Hartmann, J.-X. Cheng, and P. S. Low, “In vivo quantitation of rare circulating tumor cells by multiphoton intravital flow cytometry,” Proc. Natl. Acad. Sci. U.S.A. 104(28), 11760–11765 (2007). [CrossRef] [PubMed]
11. S. Boutrus, C. Greiner, D. Hwu, M. Chan, C. Kuperwasser, C. P. Lin, and I. Georgakoudi, “Portable two-color in vivo flow cytometer for real-time detection of fluorescently-labeled circulating cells,” J. Biomed. Opt. 12(2), 020507 (2007). [CrossRef] [PubMed]
12. C. Li, R. K. Pastila, C. Pitsillides, J. M. Runnels, M. Puoris’haag, D. Côté, and C. P. Lin, “Imaging leukocyte trafficking in vivo with two-photon-excited endogenous tryptophan fluorescence,” Opt. Express 18(2), 988–999 (2010). [CrossRef] [PubMed]
13. H. Choi, D. N. Wadduwage, T. Y. Tu, P. Matsudaira, and P. T. C. So, “Three-dimensional image cytometer based on widefield structured light microscopy and high-speed remote depth scanning,” Cytometry A 87(1), 49–60 (2015). [CrossRef] [PubMed]
14. G. He, D. Xu, H. Qin, S. Yang, and D. Xing, “In vivo cell characteristic extraction and identification by photoacoustic flow cytography,” Biomed. Opt. Express 6(10), 3748–3756 (2015). [CrossRef] [PubMed]
15. H. Lee, C. Alt, C. M. Pitsillides, M. Puoris’haag, and C. P. Lin, “In vivo imaging flow cytometer,” Opt. Express 14(17), 7789–7800 (2006). [CrossRef] [PubMed]
16. C. M. Pitsillides, J. M. Runnels, J. A. Spencer, L. Zhi, M. X. Wu, and C. P. Lin, “Cell labeling approaches for fluorescence-based in vivo flow cytometry,” Cytometry A 79A(10), 758–765 (2011). [CrossRef] [PubMed]
17. E. R. Tkaczyk and A. H. Tkaczyk, “Multiphoton flow cytometry strategies and applications,” Cytometry A 79(10), 775–788 (2011). [CrossRef] [PubMed]
18. W. Denk, J. H. Strickler, and W. W. Webb, “Two-photon laser scanning fluorescence microscopy,” Science 248(4951), 73–76 (1990). [CrossRef] [PubMed]
19. Y.-C. Chang, J. Y. Ye, T. P. Thomas, Z. Cao, A. Kotlyar, E. R. Tkaczyk, J. R. Baker Jr, and T. B. Norris, “Fiber-optic multiphoton flow cytometry in whole blood and in vivo,” J. Biomed. Opt. 15(4), 047004 (2010). [CrossRef] [PubMed]
20. R. A. Niesner and A. E. Hauser, “Recent advances in dynamic intravital multi-photon microscopy,” Cytometry A 79(10), 789–798 (2011). [CrossRef] [PubMed]
21. E. J. Botcherby, C. W. Smith, M. M. Kohl, D. Débarre, M. J. Booth, R. Juškaitis, O. Paulsen, and T. Wilson, “Aberration-free three-dimensional multiphoton imaging of neuronal activity at khz rates,” Proc. Natl. Acad. Sci. U.S.A. 109(8), 2919–2924 (2012). [CrossRef] [PubMed]
22. F. Anselmi, C. Ventalon, A. Bègue, D. Ogden, and V. Emiliani, “Three-dimensional imaging and photostimulation by remote-focusing and holographic light patterning,” Proc. Natl. Acad. Sci. U.S.A. 108(49), 19504–19509 (2011). [CrossRef] [PubMed]
23. G. Katona, G. Szalay, P. Maák, A. Kaszás, M. Veress, D. Hillier, B. Chiovini, E. S. Vizi, B. Roska, and B. Rózsa, “Fast two-photon in vivo imaging with three-dimensional random-access scanning in large tissue volumes,” Nat. Methods 9(2), 201–208 (2012). [CrossRef] [PubMed]
24. G. Duemani Reddy, K. Kelleher, R. Fink, and P. Saggau, “Three-dimensional random access multiphoton microscopy for functional imaging of neuronal activity,” Nat. Neurosci. 11(6), 713–720 (2008). [CrossRef] [PubMed]
25. L. Kong, J. Tang, J. P. Little, Y. Yu, T. Lämmermann, C. P. Lin, R. N. Germain, and M. Cui, “Continuous volumetric imaging via an optical phase-locked ultrasound lens,” Nat. Methods 12(8), 759–762 (2015). [CrossRef] [PubMed]
26. M. Duocastella, G. Vicidomini, and A. Diaspro, “Simultaneous multiplane confocal microscopy using acoustic tunable lenses,” Opt. Express 22(16), 19293–19301 (2014). [CrossRef] [PubMed]
27. R. Turcotte, C. Alt, L. J. Mortensen, and C. P. Lin, “Characterization of multiphoton microscopy in the bone marrow following intravital laser osteotomy,” Biomed. Opt. Express 5(10), 3578–3588 (2014). [CrossRef] [PubMed]
28. C. Lo Celso, J. W. Wu, and C. P. Lin, “In vivo imaging of hematopoietic stem cells and their microenvironment,” J. Biophotonics 2(11), 619–631 (2009). [CrossRef] [PubMed]
29. S. Jung, J. Aliberti, P. Graemmel, M. J. Sunshine, G. W. Kreutzberg, A. Sher, and D. R. Littman, “Analysis of fractalkine receptor cx(3)crl function by targeted deletion and green fluorescent protein reporter gene insertion,” Mol. Cell. Biol. 20(11), 4106–4114 (2000). [CrossRef] [PubMed]
30. O. Otto, P. Rosendahl, A. Mietke, S. Golfier, C. Herold, D. Klaue, S. Girardo, S. Pagliara, A. Ekpenyong, A. Jacobi, M. Wobus, N. Töpfner, U. F. Keyser, J. Mansfeld, E. Fischer-Friedrich, and J. Guck, “Real-time deformability cytometry: On-the-fly cell mechanical phenotyping,” Nat. Methods 12(3), 199–202 (2015). [CrossRef] [PubMed]
31. D. R. Gossett, H. T. K. Tse, S. A. Lee, Y. Ying, A. G. Lindgren, O. O. Yang, J. Rao, A. T. Clark, and D. Di Carlo, “Hydrodynamic stretching of single cells for large population mechanical phenotyping,” Proc. Natl. Acad. Sci. U.S.A. 109(20), 7630–7635 (2012). [CrossRef] [PubMed]
32. T. Lämmermann, P. V. Afonso, B. R. Angermann, J. M. Wang, W. Kastenmüller, C. A. Parent, and R. N. Germain, “Neutrophil swarms require ltb4 and integrins at sites of cell death in vivo,” Nature 498(7454), 371–375 (2013). [CrossRef] [PubMed]
33. N. Aceto, A. Bardia, D. T. Miyamoto, M. C. Donaldson, B. S. Wittner, J. A. Spencer, M. Yu, A. Pely, A. Engstrom, H. Zhu, B. W. Brannigan, R. Kapur, S. L. Stott, T. Shioda, S. Ramaswamy, D. T. Ting, C. P. Lin, M. Toner, D. A. Haber, and S. Maheswaran, “Circulating tumor cell clusters are oligoclonal precursors of breast cancer metastasis,” Cell 158(5), 1110–1122 (2014). [CrossRef] [PubMed]
34. E. R. Tkaczyk, C. F. Zhong, J. Y. Ye, A. Myc, T. Thomas, Z. Cao, R. Duran-Struuck, K. E. Luker, G. D. Luker, T. B. Norris, and J. R. Baker Jr., “In vivo monitoring of multiple circulating cell populations using two-photon flow cytometry,” Opt. Commun. 281(4), 888–894 (2008). [CrossRef] [PubMed]
35. S. R. Walsh, E. J. Cook, F. Goulder, T. A. Justin, and N. J. Keeling, “Neutrophil-lymphocyte ratio as a prognostic factor in colorectal cancer,” J. Surg. Oncol. 91(3), 181–184 (2005). [CrossRef] [PubMed]
36. S. Tokgoz, M. Kayrak, Z. Akpinar, A. Seyithanoğlu, F. Güney, and B. Yürüten, “Neutrophil lymphocyte ratio as a predictor of stroke,” J. Stroke Cerebrovasc. Dis. 22(7), 1169–1174 (2013). [CrossRef] [PubMed]
37. C. W. Freudiger, W. Min, B. G. Saar, S. Lu, G. R. Holtom, C. He, J. C. Tsai, J. X. Kang, and X. S. Xie, “Label-free biomedical imaging with high sensitivity by stimulated raman scattering microscopy,” Science 322(5909), 1857–1861 (2008). [CrossRef] [PubMed]
38. W. Min, S. Lu, S. Chong, R. Roy, G. R. Holtom, and X. S. Xie, “Imaging chromophores with undetectable fluorescence by stimulated emission microscopy,” Nature 461(7267), 1105–1109 (2009). [CrossRef] [PubMed]
39. J. Tang, R. N. Germain, and M. Cui, “Superpenetration optical microscopy by iterative multiphoton adaptive compensation technique,” Proc. Natl. Acad. Sci. U.S.A. 109(22), 8434–8439 (2012). [CrossRef] [PubMed]
40. L. Kong and M. Cui, “In vivo fluorescence microscopy via iterative multi-photon adaptive compensation technique,” Opt. Express 22(20), 23786–23794 (2014). [CrossRef] [PubMed]
41. L. Kong and M. Cui, “In vivo neuroimaging through the highly scattering tissue via iterative multi-photon adaptive compensation technique,” Opt. Express 23(5), 6145–6150 (2015). [CrossRef] [PubMed]
42. L. Kong and M. Cui, “In vivo deep tissue imaging via iterative multi-photon adaptive compensation technique,” IEEE J. Sel. Top. Quantum Electron. 99, 1 (2015). [CrossRef]
43. L. Kong, J. Tang, and M. Cui, “In vivo volumetric imaging of biological dynamics in deep tissue via wavefront engineering,” Opt. Express 24(2), 1214–1221 (2016). [CrossRef] [PubMed]
44. J. H. Park, W. Sun, and M. Cui, “High-resolution in vivo imaging of mouse brain through the intact skull,” Proc. Natl. Acad. Sci. U.S.A. 112(30), 9236–9241 (2015). [CrossRef] [PubMed]
45. I. Veilleux, J. A. Spencer, D. P. Biss, D. Cote, and C. P. Lin, “In vivo cell tracking with video rate multimodality laser scanning microscopy,” IEEE J. Sel. Top. Quantum Electron. 14(1), 10–18 (2008). [CrossRef]