Imaging embryonic development with ultrashort pulse microscopy

Abstract. We report the application of ultrashort pulse microscopy (UPM) for integrated imaging of embryonic development at the tissue, cell, and molecular length scales. The UPM is a multimodal imaging platform that utilizes the broad-power spectrum and high-peak power of 10-fs pulses to render two-photon excited signals and the short coherence gate of such pulses to render optical coherence signals. We show that ultrashort pulses efficiently excite cellular autofluorescence in developing zebrafish embryos such that tissues are readily visualized and individual cells can be monitored, providing a potential method for label-free cell tracking. We also show the ability of ultrashort pulses, without tuning, to excite a broad spectrum of fluorescent protein variants for tracking genetically labeled cell lineages in live embryos, with no apparent damage to the embryos. Molecular information at the mRNA transcript level can also be obtained from embryos that have been stained to reveal the localization of the expression of a gene using NBT/BCIP, which we show can be detected with three-dimensional resolution using a combination of two-photon and optical coherence signals. From this demonstration, we conclude that UPM is an efficient and a powerful tool for elucidating the dynamic multiparameter and multiscale mechanisms of embryonic development.


Introduction
Embryonic development is the shaping of a single cell into a functioning organism through the interdependent forces of cell proliferation, migration, and genetic specification.Our understanding of embryonic development represents a conceptual synthesis of what has been shown through many different types of direct and indirect experimentations, enabled by the evolving experimental tools available to developmental biologists.Transplantation technology allowed the astounding discovery of the developmental organizer, 1 mutagenesis uncovered many key developmental genes/networks and their roles, 2,3 genomic sequencing revealed unexpected levels of evolutionary conservation in the development of the animal body plan, 4,5 and the discovery of fluorescent proteins 6,7 has enabled the construction of genetic reporters for labeling and tracking specific populations of cells in developing embryos. 85][16] The foundation of this process of discovery is visualization, a function of contrast mechanisms and the ability to meaningfully capture them.Consequently, the progress in understanding embryonic development and downstream clinical application of that basic knowledge is intimately tied to the progress in imaging technology.
To emphasize the demands on an imaging system that would attempt to directly capture them, the events of embryonic development can be briefly summarized as follows.When an egg (a zebrafish egg of approximately 700 μm in diameter 17,18 ) is fertilized, the resulting embryo begins a stage of rapid cell division (15 to 30 min per division 19 ), known as the cleavage stage.Toward the end of this stage, the embryo transits from relying on pre-existing molecular cues deposited in the egg to utilizing its own genomic machinery (thousands of genes 20,21 ) and shortly thereafter, the three germ layers are distinguished during the process of gastrulation, where rapid cell movements (few μm∕ min 22 ) accompanied by spatiotemporal varying patterns of gene expression result in recognizable anterior-posterior and dorsal-ventral embryonic axes.Continued cell proliferation, slowed cell movements (∼0.5 μm∕ min 22 ), and further refinement of gene expression shape the body plan composed of segments of appropriately fated cell lineages (now tens of thousands of cells 22 ).From this point, within each segment cells differentiate and take on specific morphologies and functions, as individual organs mature.Perhaps not surprisingly, despite a detailed framework and list of components, these events as a whole have yet to be fully observed or understood.
Each aspect of embryonic development, then, presents unique challenges.Its dynamism requires imaging to be nondestructive and preferably noninvasive.Its fundamental unit is the cell, but subcellular events, such as cell division and mitotic spindle orientation, are vital to proper embryonic development, as they are collective cell behaviors behind tissue and organ morphogenesis.Adding to the complexity is the daunting number of involved genes that have disparate spatial and temporal expression patterns.All this activity occurs over a relatively large three-dimensional (3-D) space, placing further demands on imaging capabilities.The challenge to modern imaging systems, then, is to rapidly and noninvasively acquire large, high resolution, multichannel volumetric images from a live developing embryo. 23oward this goal, two-photon microscopy (2PM) has become favored over one-photon approaches because of the advantages of near-infrared (NIR) light, including (1) less photodamage to embryos, allowing longer timelapse imaging sessions, (2) less scattering, allowing deeper imaging into embryos, especially more opaque embryos, and (3) the option of simultaneously exciting multiple fluorescent species in the biological window (near IR), without the excitation light (750 to 1000 nm) interfering with imaging signals (400 to 700 nm). 24As 2PM techniques have become more widespread, so has ultrafast laser technology, with most imaging systems typically employing pulses from a Ti:Sapphire oscillator with a temporal duration of 100 to 170 fs. 25,26Despite the advantages gained by adopting 2PM, however, it has remained challenging to acquire high-resolution, 3-D multichannel movies of embryonic development.
Microscopists are well acquainted with the trade-offs encountered when attempting to optimize imaging parameters due to a photon economy, i.e., limited by the organism's tolerance of the probing light radiation.For example, to increase imaging speed to capture more rapid cell movements during embryogenesis, narrowband NIR pulses will generally be tuned to the optimum 2P absorption peak of a single-fluorescent marker and all emitted light will be collected as a single channel. 27This approach provides "in toto" characterization of cell proliferation, death, and migration, but provides little information about the shaping of the embryonic tissue and none about the genetic specification of cells.With a transition from single-point scanning to line scanning, the excitation light source, as in digital scanned light-sheet microscopy (DSLM), can possibly add another channel. 28,29However, imaging two or three genetic markers in four-dimensional (space, time) are the current technological limit, 27,30 while the genetic regulatory networks governing embryonic development contain vastly more components.
Optical coherence tomography/microscopy (OCT/M) is a linear imaging technique that images the morphology of biological samples based on changes in index of refraction that frequently employ NIR laser pulses. 31Because of the high sensitivity of the interferometric detection scheme employed in OCT/M, this approach has been used to image extraordinarily fast movements during embryogenesis such as heart tube pumping and development. 32,33Alternatively, because of the enhanced penetration of both the input and output NIR signals, OCT has also been used for longitudinal imaging of mouse embryos in utero, where the uterus is temporarily dissected from the mother. 34,35Since 2PM and OCT/M can be relatively easily integrated into a multimodal platform, it would be interesting to explore potential applications of OCT/M to additionally render molecularly specific signals in the context of embryonic development.
NBT/BCIP is a robust and strongly scattering precipitate that is commonly used to mark domains of gene expression through a standard immunohistochemical technique called in situ hybridization (ISH). 36Though ISH requires embryos to be fixed, it provides fundamentally unique data from 2PM approaches that track fluorescent protein markers because it provides a more accurate view of gene transcript localization at a specific developmental stage.Such data will be complementary to that obtained by time-lapse 2PM, while also being relatively easily obtained with the same imaging platform.
Here, we present the theoretical benefit of using ultrashort pulses, on the order of 10 fs, for imaging embryonic development and demonstrating the benefit manifest such as unique capabilities for exciting and rendering a wide variety of optical signals from developing embryos simultaneously.In addition, we demonstrate these capabilities in a system, where dispersion compensation is implemented in a relatively simple and user-friendly manner using double-chirped mirrors. 37Our multimodal instrument can be used for 2PM (and other forms of nonlinear optical microscopy such as second-harmonic generation) with standard two-channel detection or with spectral detection, as well as for OCM, and we refer to these techniques collectively as ultrashort pulse microscopy (UPM).Within a single experiment, UPM provides the opportunity to better link the behaviors of genetically labeled cell lineages with the morphogenesis of tissues and organs during development.[40] 2 Theory For imaging embryonic development, ultrashort pulses provide the unique capabilities of simultaneous excitation of a broad range of fluorophores and a short coherence gate for high-resolution optical coherence microscopy of tissues or deposited scattering precipitates.

Two-Photon Excitation by Ultrashort Pulses
Temporally ultrashort pulses, such as in the 10-fs regime, are advantageous compared with conventional 100 to 170-fs pulses used in 2PM applications, because (1) for a given pulse energy, the peak power of an ultrashort pulse will be higher because the energy is compressed in time, and (2) an ultrashort pulse necessarily contains a broader range of wavelengths that can together create two-photon transitions. 41,42Using a semi-classical model of a two-photon transition, where the electric field is treated classically and the molecule is treated quantum mechanically as a twolevel system, this benefit can be quantified in terms of relative two-photon transition probabilities.
To compare the theoretical transition probabilities of our sub-10-fs pulses with conventional pulses, we first measured the spectral intensity, jEðωÞj 2 , of sub-10-fs pulses output from our Ti:Sapphire oscillator.To model the electric field, EðωÞ, of alternative transform-limited pulses (all frequencies arrive at the focus in phase) of varying temporal durations, τ p , we used the Gaussian relationship shown in Eq. (1).
Two-photon transitions can occur through degenerate and nondegenerate mixings of frequency components, so if the pulse is transform limited, then all possible combinations of the frequencies within the pulse can be given by the two-photon excitation power spectrum, 43 TðωÞ, which is described by Eq. (2).
The two-photon excitation spectrum of our sub-10-fs pulse is compared with that of a 100-fs pulse in Fig. 1(a).In the context of imaging embryonic development, molecules of interest generally have a frequency-dependent response, sometimes called a two-photon action cross-section. 44To meaningfully compare the theoretical performance of these pulses, this response must be known.The yield or transition probability, 45,46 Γ, is then the overlap integral of TðωÞ and the frequency response of the molecule, γðωÞ, as shown in Eq. (3).
From Eq. ( 3), it can be seen that the probability of a twophoton transition depends on the magnitude of the area under the two-photon power spectrum curve and the degree to which this overlaps with the absorption lineshape of a given fluorophore.We used γðωÞ measured from the commonly used enhanced green fluorescent protein (eGFP) 47 to calculate transition probabilities for our sub-10-fs pulses centered at 800 nm and for 100-fs pulses centered at 940 nm, the maximum of the eGFP frequency response, illustrated in Fig. 1(b).We calculated a fivefold enhancement of yield in the case of un-tuned, sub-10-fs pulses compared with optimally tuned 100-fs pulses given the same pulse energies.In terms of imaging embryonic development, this enhancement could be used to image more rapidly or, alternatively, to decrease pulse energy to reduce stress on the developing embryo.
For imaging multiple fluorescent protein variants simultaneously, the advantage of using sub-10-fs pulses increases.A single narrowband source will necessarily be tuned to a compromised wavelength that interacts with the frequency responses of the different molecules, decreasing the yield for each, while the sub-10-fs pulse interacts substantially with each simultaneously. 48,49Thus, with ultrashort pulses, it may be more feasible to image multiple genetically marked lineages during embryonic development.

Coherence Gating with Ultrashort Pulses
The benefits of utilizing the short coherence length of ultrashort pulses to improve axial resolution of OCT/M are well known. 50Using Fourier-domain detection, the frequency components of the interferogram produced by the sample and reference arms encode the spatial profile of the biological sample along the optical axis (A-line), within the confocal parameter of the imaging objective. 51Rapid 3-D imaging can be achieved by scanning this A-line across the sample in two dimensions.Thus, for OCT imaging, a low-numerical aperture (NA) objective resulting in a large depth of field is desirable, as long as the diffraction-limited lateral resolution remains comparable with the axial resolution, Δz, 51,52 which depends on the source bandwidth at full-widthhalf-maximum (FWHM), as shown in Eq. ( 4).
Our sub-10-fs pulses have a FWHM of Δλ ¼ 133 nm centered at λ 0 ¼ 800 nm corresponding to a theoretical axial resolution of 2.12 μm, whereas a more conventional pulse with a bandwidth of 10 nm at this central wavelength will Fig. 1 Advantage of using ultrashort pulses in two-photon microscopy (2PM).(a) Two-photon power spectra, T ðωÞ, for 10 fs (black line) versus 100 fs (blue line) pulses with the same pulse energy centered at 800 nm.The area under the curve represents the ability to create two-photon transitions in materials with frequency-independent response, γðωÞ.(b) T ðωÞ for 10 fs (black line) pulses at 800 nm and 100 fs pulses (blue line) tuned to 940 nm, the maxima of γðωÞ for green fluorescent protein (GFP) (dotted green line).The ability to create two-photon transitions for GFP, where the response is now frequency dependent, is the overlap integral of T ðωÞ and γðωÞ, which is five times greater in the case of untuned, transform-limited pulses at 10-fs pulses at 800 nm versus optimally tuned 100-fs pulses.
have an axial resolution of 28 μm.This resolution difference is significant when imaging individual cells and small-scale tissue structures in developing embryos.In the context of combining OCM with 2PM, the resolution offered by the use of ultrashort pulses also better matches the two-photon interaction volume resulting from a tightly focusing high NA objective, so simultaneously acquired 2PM and OCM images are co-registered. 53Methods

Instrumentation
Our multimodal UPM system, summarized in Fig. 2, combines 2PM and Fourier-domain OCT/M in a single platform with the option of either two-channel (with OCM) or spectral detection (stand alone) capabilities for 2PM.Detailed more in these reports, 37,53,54 we describe the instrument here briefly.The sub-10-fs pulses from a Kerr-lens mode-locked Ti:Sapphire oscillator (Femtolasers, Vienna, Austria) at a 800-nm central wavelength with 133-nm bandwidth (FWHM) are precompensated using double-chirped mirrors (GSM 270, Femtolasers) and coupled to the multimodal imaging system.For combined 2PM-OCM, a 5% beamsplitter splits the beam into signal and reference arms in a Michelson interferometer configuration.The signal arm is coupled using galvanometer-driven mirrors (Cambridge Technology, Cambridge, Massachusetts) into an upright microscope (Axioskop2, MAT, Carl Zeiss, Thornwood, New York), expanded, and directed by a 635-nm short-pass dichroic mirror (Chroma Technologies, Bellows Falls, New York) to the back aperture of the water-immersion imaging objective (20×, 1.0 NA or 40×, 0.8 NA, Carl Zeiss), which focuses the pulses to a diffraction-limited spot that is scanned in the x-y plane across the sample.The backscattered near-infrared light is collected by the imaging objective and reflected back to the beamsplitter, where it recombines with the reference arm and is coupled into a single-mode fiber and sent to a home-built spectrometer consisting of a collimating lens, reflective grating (1200 lines∕mm), focusing lens, and CCD array (Basler, Exton, Pennsylvania).Dispersion is matched in the reference arm using a combination of a prism pair and BK7 glass.For OCM, the interferogram obtained from each pixel is integrated into a single intensity value to build up a 256 × 256 two-dimensional (2-D) image rendered en face.Simultaneously, 2PM signals are collected by the imaging objective and directed onto the two PMT detectors (Hammamatsu, Bridgewater, NJ) using appropriate dichroic mirrors and bandpass filters (Chroma).Photon counts are saved as intensity values.
For spectral 2PM, the 2PM signals are collected by the imaging objective, coupled into a multimode fiber, and sent to a home-built spectrometer consisting of a collimating lens, diffraction grating (700 lines∕mm), focusing lens, and a 16-channel multianode PMT array (Hammamatsu).The 3-D images are built up by translating the sample in the z-direction along the optical axis.Instrumentation control and data acquisition are performed using custom LabVIEW software (National Instruments, Austin, Texas).

Label-Free Imaging of Fixed and Live Wild-Type Embryos
Zebrafish were maintained and bred according to the standard protocols. 55For imaging, wild-type embryos were treated with 0.003% (w/v) phenylthioluria, to suppress the formation of pigment cells (except for label-free imaging of melanocytes), and mounted in 1.2% low-melt agarose hydrated with Instant Ocean (United Pet Group, Blacksburg, VA).Live wild-type embryos were imaged at 25°C with time-lapse imaging performed at 27°C for 15to 30-min intervals.In some cases, embryos were fixed in 4% (w/v) paraformaldehyde overnight at 4°C and stored in phosphate-buffered saline (PBS) before imaging.Average power of the excitation beam was attenuated to 25 to 30 mW on the sample (not accounting for scattering due to the agarose embedding) and pixel dwell time ranged from 120 to 480 μs.In most cases, lateral resolution was limited by digital resolution but kept to the order of 1 to 2 μm, while axial step size was generally 3 μm, as recommended by Keller et al., 22 to capture all cells in the imaging volume.Signals were collected into a single channel using a BG39 filter (Schott, Elmsford, NY).

Imaging Transgenic Embryos
Stable transgenic lines expressing eGFP under a wnt1 regulatory element that recapitulates wnt1 expression in the midbrain and at the midbrain-hindbrain boundary (Lilie et al., unpublished) were used to demonstrate time-lapse imaging of genetically labeled cell lineages in live zebrafish embryos.Transgenic fish were intercrossed, and embryos were imaged during the segmentation period of development as described above, with stacks acquired at 15-min intervals.Autofluorescence and eGFP signals were separated with a 490nm long-pass dichroic mirror (Chroma) and further discriminated with 450∕60 and 525∕50 nm band-pass filters (Chroma), respectively.For demonstrating excitation and spectral detection of a variety of fluorescent protein variants within zebrafish embryos using ultrashort pulses, wild-type embryos were injected at the one-cell stage with plasmid DNA containing a fluorescent protein gene (eBFP2, CFP, mOrange, and mRFP) or a fluorescent protein fusion (citrine-H2B) under a cmv promoter, resulting in mosaic expression throughout cells of the embryo at 24 hours post-fertilization (hpf).For acquiring spectral images of eGFP, the aforementioned stable transgenic lines were used.Embryos were imaged at 24 hpf, and fluorescent signals were collected with the spectral detection scheme described in Sec.3.1.
Fig. 2 Multiodal ultrashort pulse microscopy (UPM).A single oscillator outputting 10-fs pulses centered at 800 nm at 75 MHz is coupled to a microscopy system with several detection modes.Pulses are precompensated for dispersion introduced by the microscopy system with double-chirped mirrors.

In Situ Hybridization and Imaging NBT/BCIP with 3-D Resolution
The ISH was performed as previously described 56 for either krox20 or egfp mRNAs.Embryos stained with NBT/BCIP and hydrated in PBS were mounted as described above and imaged in the combined 2PM-OCM modality with an average power of 35 mW on the sample (not accounting for scattering due to the agarose embedding), a pixel dwell time of 480 μs, and a step size of 3 μm.

Image Processing
Images were processed using custom MATLAB routines and openly available image-processing software ImageJ and FIJI. 573-D multichannel reconstructions were rendered in V3D. 58Nuclear segmentation from autofluorescence signal was performed manually.Spectral unmixing was performed with custom MATLAB routines that utilized spectral standards of the components to be unmixed (obtained online from Invitrogen and Molecular Probes or measured experimentally).The algorithm creates a set of fits, FðλÞ, that are the linear-weighted sum of n scalar coefficients, c i , and 16-element reference spectra, R i ðλÞ, of n presumed constituents, as shown in Eq. ( 5) The solution is the set of positive coefficients that sum to 1 and produces a fit that minimizes the mean squared error between the fit and measured spectra.

Label-Free Imaging of Embryonic Tissues
To determine whether endogenous fluorescence could be used to image vertebrate development, we imaged wildtype zebrafish embryos free of exogenous labels with 10fs pulses centered at 800 nm.The zebrafish, also known as Danio rerio, is a well characterized and established model organism, especially useful for time-lapse imaging analysis of developmental mechanisms. 17e first examined embryonic development at the tissue level.At 24 hpf, the major organ systems in zebrafish embryos have taken on a recognizable morphology according to the basic vertebrate body plan, so we imaged a fixed embryo in its entirety using overlapping fields of view to observe those features visible from endogenous fluorescence.The combined image rendered using FIJI volume viewer is shown in Fig. 3(a).
We observed a strong autofluorescence signal from most of the tissues in the embryo, and developing structures were easily visible in optical sections, including structures of the eye, ear, developing brain and ventricular system, notochord, blood islands, and muscle segments [Figs.3(b)-3(f)].These data illustrate that ultrashort pulses efficiently excite autofluorescence from the constituent cells of the embryo, such that the morphogenesis of a wide variety of organ systems can be characterized quantitatively without the need for exogenous labels.Such data will be useful for the re-emerging field of physical biology that seeks to model the mechanical forces that give rise to embryonic shape 59,60 and potentially connect these forces to specific cell behaviors and patterns of gene expression.

Label-Free Imaging of Individual Cells and Their Movements
We were also interested to determine whether cellular autofluorescence signal was sufficient for tracking individual cells; so, we performed time-lapse imaging of the cranial region of zebrafish embryos during the segmentation period when neural crest cells are derived from the roof plate of the neural tube and migrate to give rise to different structures and cell types in the embryo including craniofacial cartilage and bone, melanocytes, and peripheral sensory neurons. 61As seen in Figs.4(a) and 4(b), these neural crest cells generally originate from the dorsal side of the densely packed neuroepithelium and migrate laterally and ventrally into a less-dense mesenchyme outside the neural tube.In this mesenchyme, individual neural crest cells were easily visible from autofluorescence signal and demonstrate that UPM could potentially be utilized to track the dynamic migration of different populations of neural crest, i.e., produce a highresolution neural crest fate map independent of exogenous cell labels.We also decreased the field-of-view to match our digital resolution to the diffraction-limited performance of the system and to assess the possibility of imaging individual cells within the tightly packed neuroepithelium at 24 hpf.At this time in development, the neuroepithelium is a contorted tube, highly stratified in the apical side in contact with the brain ventricular system and the basal side facing outward.Figure 4(c) shows a cartoon of this arrangement from a dorsal view, where the anterior portion of the head is directed to the left.The neural tube has been segmented into morphologically distinct regions including the forebrain (fb), midbrain (mb), and hindbrain (hb).The boundary between the midbrain and hindbrain is referred to as the midbrain-hindbrain boundary (MHB) organizer [the term organizer denotes the importance of this region for correct development of structures in the midbrain and hindbrain such as the tectum (tc) and cerebellum (cb), respectively].Figure 4(c)' shows autofluorescence from an optical section taken from this boundary region on the right side of the embryo.The large dark portion of the image in the bottom left corresponds to the fluid-filled brain ventricle, whereas the other dark spots correspond to individual cell nuclei.Figure 4(c)" is the same image inverted and marked to show the arrangement of individual cell nuclei in yellow and the outline of the neuroepithelium in blue.We were able to image a unique region of the epithelium at the boundary, where a group of cells seems to have recaptured the character of a simple onelayer epithelium with the nuclei aligned rather than staggered.This particular arrangement may provide some structural support to the MHB constriction, as the tectum and cerebellum rapidly develop on either side.We have also observed that this group of cells can be identified by their strong expression of wnt1 (Figs. 6 and 8).While there is a significant literature concerning the genetic basis underlying the structure and specification of this region, 62,63 there is relatively little known concerning details about how these genes direct individual or collective cell behaviors that achieve this complex tissue architecture.
The autofluorescence signal we detect arises predominantly from pyridine nucleotides and flavins found in the mitochondria and cytoplasm, whose 2P excitation spectra have been measured. 64While narrowband pulses can simultaneously excite key autofluorescence molecules NADH, NAD(P)H, FAD, and lipamide dehydrogenase, we have not seen reports of these signals used to measure tissue morphology and cellular movements over time in embryonic development.It may be that for a given pulse energy, ultrashort pulses provide a significant enhancement over conventional pulses such that the autofluorescence is transformed from unwanted background noise to a label-free method to integrate the cellular and tissue scale dynamics of development.

Label-Free Imaging of Pigmented Cells
Melanocytes, derived from the multipotent neural crest cell lineage, have been a prominent model for the study of cell differentiation and establishment of adult stem cell populations during development. 65We have observed strong (up to 10-fold stronger than cellular autofluorescence) 2P excited fluorescence from melanocytes (commonly called melanophores in fish), as shown in a volumetric maximum intensity projection in Fig. 5(a).Live embryos that were not treated to block melanin production were imaged at 24 hpf.We observed brightly fluorescing cells distributed throughout the cranial region underneath the periderm displaying expected differentiated cranial melanocyte morphology and positioning. 61,65Neural crest delamination and migration occur with varying dynamics along the anterior-posterior axis and, at a specific spatial location, neural crest may delaminate in temporally sequential waves.For example, in the zebrafish trunk, at least two different temporal waves give rise to neural crest migration along a medial and, later, a more lateral pathway. 66In the cranial region, these dynamics are not well characterized and could be significantly more complex.Optical sections in Figs.5(b)-5(d) reveal a striking population in the optic tectum localized at the dorsolateral hinge points of the neural tube, which appear pigmented while still integrated in the neural epithelium.This finding is unusual, because melanocytes are generally thought to differentiate from melanoblast precursors that have already delaminated from the neuroepithelium. 67Some neurons, such as dopaminergic neurons, are known to be pigmented, 68 but no pigmented neurons have been reported in this region of the optic tectum in zebrafish.One interpretation of this data is that ultrashort pulses provide significantly more sensitive detection of melanin than can be achieved with standard light microscopy and that, at least for a subpopulation of melanocytes, differentiation occurs prior to or during the epithelialto-mesenchyme transition.Efficient 2P excitation of melanin using conventional pulses and resulting in bright fluorescent signals has been reported in cell culture 69 and in a concentration-dependent manner in vivo. 70It would be interesting to use a combination of autofluorescence and melanin fluorescence to track the origin and eventual fate of this unique population, although the overwhelming melanin signal may impede the addition of any fluorescent protein gene reporters simultaneously.

Time-Lapse Imaging of Genetically Labeled Cell Lineages
We next demonstrated the imaging of the dynamics of a genetically labeled cell lineage simultaneously with morphogenetic events in the embryonic brain.Stable transgenic embryos expressing eGFP under the control of regulatory elements from the wnt1 gene were imaged during segmentation over a period of 12 to 13 h.During segmentation, wnt1 is expressed in the midbrain and at the MHB and is important for proper development of the region.In mice lacking Wnt1, the entire midbrain and hindbrain regions (known as mes/r1) are deleted. 71In zebrafish, other Wnts have redundant roles to wnt1; however, loss of these Wnts in combination also results in the deletion of midbrain and hindbrain structures. 56,72Thus, wnt1 has an important role in the specification of these embryonic brain tissues.Morphogenesis of the MHB region imaged with UPM progressed in step with previous reports. 73,74We show in Fig. 6 the dynamics of the wnt1-expressing lineage in this complex structural environment of the mes/r1 region.When the MHB is initially formed, genetic determinants from either side of the boundary initially overlap, and subsequently their expression sharpens. 63This sharpening coincides with the development of a morphologically visible constriction at the boundary by the 10-somite stage (t ¼ 0∶00).We observe eGFP positive cells in the hindbrain at this stage whose fluorescence gradually diminishes, suggesting that the boundary sharpens at least in part when cells alter their pattern of gene expression in response to a still undetermined mechanism of spatial identity acquisition.At t ¼ 0∶00, we also observe the wnt1 lineage to have a longer extent along the anterior-posterior axis in the dorsal neural tube than in ventral neural tube and find that this arrangement persists throughout the periods of segmentation and ventricle formation, suggesting that this initially labeled lineage is proliferative but relatively stationary along the anterior-posterior axis within the neuroepithelium.Intensity of the eGFP signal increases over time in cells at the boundary itself and anteriorly along the dorsal midline up to the epiphysis, indicating continuous production of reporter transcripts in these regions, which maintain wnt1 expression.At t ¼ 3∶53, we also captured the epithelialto-mesenchyme transition of a portion of the wnt1 lineage comprising neural crest cells, marked by arrowheads.After the time-lapse imaging concluded, we observed disfiguration of the tail and slight edema as expected when immobilizing zebrafish embryos with a combination of tricaine and agarose, 75 but otherwise embryos continued to develop normally through larval stages.
Photobleaching by ultrashort pulses appears to be minimal as subpopulations of the wnt1 lineage that have stopped transcribing the wnt1 reporter (as indicated by ISH at previous stages) can be tracked for several hours, despite the dilution of available eGFP molecules by cell division.As mentioned previously, ultrashort pulses do not appear to cause undue stress to the developing embryo as evidenced by normal morphological progression of the irradiated region of the embryo throughout the period of imaging and up to 7 days later.However, in the future, it may be prudent to assess the damage threshold for these organisms using ultrashort pulses similar to what has been reported for Drosophila. 76Using a combination of markers, including lethality, necrosis, and TUNEL staining, Saytashev et al. showed that for a given signal level, shorter pulses (37 fs) resulted in less thermal damage and fewer DNA strand breaks than longer pulses (100 fs).

Imaging a Variety of Labels with a Single Pulsed Source
To demonstrate the potential for capturing the dynamics of multiple fluorescently labeled genetic lineages simultaneously during embryonic development, we obtained several fluorescent protein (FP) variants and expressed them mosaically in independent zebrafish embryos through transient transgenesis.These embryos were imaged at 24 hpf using our home-built 16-channel spectral detector, spanning 350 to 630 nm for capturing 2P-excited fluorescence signals.
For each pixel in a 2-D image or voxel in a 3-D image, a 16-channel lambda stack is acquired from live embryos.The average of the spectral profile for different regions in the embryos, as shown in Fig. 7, was compared with delineated autofluorescence contributions (acquired from areas of the embryo not expressing the injected fluorescent protein variant) and spectral signatures from eBFP2, CFP, eGFP, citrine, mOrange, and mRFP.Except for eGFP, the amount of plasmid injected was the same for each FP variant, and thus, the relative intensity of fluorescent protein to autofluorescence can be interpreted using the 2P overlap integral among the variants with our two-photon excitation power spectrum, as in Eq. ( 3).As expected, for pulses centered at 800 nm, the overlap appears to be more significant for blue variants than for red, although the other red variants such as TagRFP1 77 have been reported to have a significant 2P action cross-section near 800 nm and, thus, may be better suited for this application.These spectral signatures show that we can indeed efficiently excite a wide spectrum of fluorescent protein variants simultaneously, as all spectra were obtained with the same untuned ultrashort pulses.Once a suitable transgenic line is created, it may be then possible to image the dynamics of up to five different genetically labeled lineages.
To demonstrate the potential to accurately delineate these simultaneously excited overlapping spectral emissions, we performed spectral time-lapse imaging with simultaneous excitation of closely overlapping eGFP and citrine (a yellow fluorescent protein) emissions and separated them with linear unmixing.The unmixing algorithm assumed three constituents: (1) autofluorescence, (2) eGFP, and (3) citrine.Autofluorescence reference spectra were measured directly from wild-type embryos, whereas eGFP and citrine reference spectra are available online (Molecular Probes and Invitrogen).The stable transgenic line marking the wnt1 lineage with eGFP was injected with plasmid DNA encoding citrine-H2B at the one-cell stage, resulting in differential distribution of these two fluorescent markers across the MHB at 24 hpf.Summed images are shown in the top row of Fig. 8. Individual cells on either side of the boundary were tracked (green and yellow arrowheads) and their spectral profiles unmixed revealing the expected higher levels of eGFP expression in the cell tracked in the midbrain (bottom row) compared with a cell tracked in the hindbrain (middle row) reflected by the respective areas underneath the unmixed eGFP reference profiles.Within the neuroepithelium, cells undergo interkinetic migration that is characterized by cellular and nuclear movements between the basal and apical surfaces, with cell division occurring at the apical side.It would be interesting to determine whether the fluorescence contribution from the citrine-H2B fusion protein, which labels nucleosomes, could be quantified as a means to track cell cycle dynamics.Intensity from the marker would presumably fluctuate as nuclear chromatin replicates and subsequently condenses prior to chromosome segregation.

Imaging In Situ Hybridization with Three-Dimensional Resolution
Having established that the broad 2P power spectrum of ultrashort pulses is useful for integrating molecular, cellular, and morphogenetic information in live organisms over time,  To demonstrate how this OCM technique can be used to complement 2PM data, we performed ISH to detect egfp transcripts in the stable transgenic line, where eGFP is driven by wnt1 enhancer elements, and imaged those embryos with combined 2PM-OCM at 24 hpf [Figs.9(e)-9(h)].Comparing the distribution of egfp transcripts versus eGFP protein revealed that a substantial portion of the wnt1 lineage stopped expressing the wnt1 reporter, as suggested by tracking eGFP protein signal over time.Although this type of analysis can only be performed for a single time point in a given embryo, it will be useful for characterization of the genetic networks governing cell fate decisions during embryogenesis, as they involve not only the turning on of certain genes, but also the silencing of others.

Conclusion
We have provided both theoretical arguments and data in favor of the use of ultrashort pulses for imaging the dynamic multiparameter mechanisms of embryonic development that begin at the molecular scale and culminate in the shaping of an entire organism.In particular, we have shown (1) the ability to image tissue morphogenesis and individual cells and  their movements within developing embryos from intrinsic autofluorescence, i.e., label free, (2) the ability to simultaneously track tissue morphogenesis and genetically labeled cell lineages, (3) the ability to excite a wide range of fluorescent protein markers simultaneously for tracking multiple genetically labeled lineages in their morphological context, and (4) the ability to image mRNA transcript localization marked by NBT/BCIP staining with 3-D resolution.These capabilities, summarized in Fig. 10, are available with a simple implementation for dispersion control, double-chirped mirrors which are relatively easy to align and produce less unwanted third-order dispersion (TOD) than prism pairs. 80n reality, it is highly unlikely we have achieved the theoretical enhancement predicted for transform-limited pulses.It will be interesting in the future to implement more sophisticated pulse-control methods for managing TOD such as MIIPs pulse shaping, 81 MEMs, 82 or a grating/prism compressor 83 that could be used to tune second-and third-order dispersions simultaneously.Such improvements may be used to more precisely account for varying imaging conditions and biological samples, further increase imaging speed, or if needed, decrease irradiation to the embryo while maintaining the same signal level.

Fig. 3
Fig. 3 Label-free imaging of embryo and organ morphology with UPM.(a) Lateral view of a volume rendering of endogenous signals excited by ultrashort pulses in a zebrafish embryo at 24 hours post fertilization (hpf), also referred to as the prim-5 stage.Image is the product of fusing nine different volumes.(b-f) Optical sections of the same embryo detailing various forming organs including (b) the eye, composed of the lens (l) and retina (r), (c) the developing otic cup (oc), (d) different brain segments, including the forebrain (fb), midbrain (mb), and hindbrain (hb), and ventricular spaces (v), (e) trunk region including notochord (n) and blood islands (bi), and (f) more laterally in the trunk region, myotome (m) or muscle segments.Scale bars are 100 μm.

Fig. 4
Fig. 4 Label-free imaging of individual cells and their movements with UPM.(a) For reference, a cartoon of a dorsal view of the embryonic brain including the forebrain (fb), eye (e), neural crest (nc), and midbrain (mb).Left-right symmetry can be assumed.The rest of (a) are the optical sections ðx; yÞ depicting coordinated migration of neural crest cells in a low-density mesenchyme tracked during the segmentation period of development (t ¼ 0∶00 is the 10 somite stage).Images were rendered en face in this orientation.(b) Cartoon depicting a transverse view of the embryonic neural tube (nt).The rest of (b) are the optical sections ðy; zÞ of neural crest cells migrating lateral to the neural tube.(c) Cartoon depicting morphology of the midbrain-hindbrain boundary (MHB) at 24 hpf.(c') Optical section of the area outlined in gray.At this stage, the neuroepithelium is stratified and densely packed, but autofluorescence excited by ultrashort pulses can still be used to identify individual cells.Nuclei, lacking endogenous fluorophores, remain dark.(c") Different cellular organizations in different brain compartments, such as the tectum (tc) and cerebellum (cb), including a special nonstratified epithelium at the MHB.Scale bars are 100 μm.

Fig. 5
Fig. 5 Label-free imaging of melanocytes and other pigment cells.(a) Dorsal view of a three-dimensional (3-D) maximum intensity projection of the embryonic zebrafish brain at 24 hpf.At this stage, cells derived from the neural crest have differentiated into melanocytes in a characteristic pattern in the cranial region.A strong signal (10-fold brighter than tissue autofluorescence) excited by ultrashort pulses matches this characteristic pattern.(b-d) Optical sections (b-lateral, c-dorsal, and d-transverse) from volume in (a) also reveal a population of pigmented cells in the neural epithelium that is localized at dorsolateral hinge points in the tectal primordium.Cells appear to be undergoing an epithelial-to-mesenchyme transition.Scale bars are 100 μm.

Fig. 7
Fig.7Excitation and spectral detection of fluorescent protein variants transiently expressed in zebrafish embryos.First panel shows a representative image from an embryo mosaically expressing a fluorescent protein variant (citrine-H2B, in this case, segmented and superimposed in yellow).Dotted line represents the embryonic midline, and boxes show areas from which pixel-by-pixel 16-channel spectra have been summed.Subsequent graphs show the spectral data from embryos injected with different variants, spanning a range from blue to red fluorescent proteins.One spectral profile obtained from a region of the embryo with little or no exogenous label (black line) is compared with an area where the injected fluorescent protein is present (appropriately colored line).The autofluorescence (AF) contributes a broad and somewhat indistinct spectral profile, while peaks from individual fluorescent proteins are readily identifiable at their expected locations.Scale bar is 50 μm.

Fig. 6
Fig. 6 Imaging genetically labeled cell lineages in their morphological context.First panel shows a cartoon of the region of the brain depicted in subsequent 3-D maximum intensity projections from a time-lapse experiment on a stable transgenic line expressing eGFP protein under the control of a wnt1 regulatory element.Time-lapse follows the midbrain (mb) and hindbrain (hb) regions during the formation of the midbrain-hindbrain boundary (MHB) constriction and brain ventricle formation.t ¼ 0∶00 corresponds to the 10 to 12 somite stages, when the neural tube is still closed.The wnt1 population (green) is tracked with 3-D resolution alongside the dramatic reorganization of neural morphology visible by autofluorescence (white).Arrows point to the migration of neural crest cells that are derived from a wnt1 lineage.Scale bar is 100 μm.

Fig. 8
Fig.8Spectral unmixing of overlapping fluorescence spectra from live cells at the midbrain-hindbrain boundary.Stable transgenic embryos expressing eGFP in the wnt1 lineage were injected at the one-cell stage with plasmid DNA containing a citrine-H2B fusion protein to mosaically mark inidividual nuclei.Top row shows summed spectral image, and arrowheads point to cells tracked in the midbrain (green arrowhead) and hindbrain (yellow arrohead).The spectral profiles of these cells were unmixed by varying weighting parameters for a linear sum of reference spectra (solid curves) to minimize the mean squared error between the proposed fit (dotted line) and measured data (solid line).Scale bar is 50 μm.

Figure 9 (
a) shows NBT/ BCIP staining marking the location of krox20, which is expressed in rhombomeres 3 and 5 of the developing hindbrain, in an easily distinguisable pattern.Figure9(b)shows the quenching that occurs in that region.Figure9(c)shows the OCM signal from both the tissue and NBT/BCIP stain.The overlay, shown in Fig.9(d), highlights regions stained by NBT/BCIP (low autofluorescence and high OCM signal) in red, while tissue architecture (marked by comparable OCM and autofluorescence signal) is shown in green to yellow.

Fig. 9
Fig. 9 Imaging gene expression domains with multimodal UPM.(a) Light microscopy image of an in situ hybridization (ISH), where the expression domain of the krox20 gene is marked by NBT/BCIP staining (arrowheads).Combined 2PM-OCM from the embryo shows that NBT/BCIP quenches tissue autofluorescence (b), while it enhances scattering detected by OCM (c).When these signals are combined (d), the NBT/BCIP positive domains are marked by red, while the rest of the tissue is visible as either green (low OCM signal in tissue relative to autofluorescence) or yellow (comparable OCM signal in tissue relative to autofluorescence).(e-g) Optical sections of combined 2PM-OCM using ultrashort pulses from the midbrain-hindbrain boundary, showing where egfp transcripts are present in NBT/BCIP-stained embryos with 3-D resolution.Sections correspond to placement of cross-hairs.(e) En-face rendered image dependent on lateral resolution.Transcripts are localized to the MHB [(f) transverse and (g) lateral].Optical sections depend on axial resolution.Asterisks point to characteristic lessening of staining intensity in the basal portion of the neural tube at the MHB.Arrow points to NBT/BCIP staining in the most ventral portion of the neural tube, at a depth of about 180 μm.(h) The 3-D maximum intensity projection of the same dataset as (e-g).Scale bars are 50 μm.

Fig. 10
Fig. 10 Integrated imaging of multiscale embryonic development using multimodal UPM.Left diagram depicts the components that drive embryogenesis in a dynamic, responsive manner.Right diagram highlights the components that can be visualized simultaneously with multimodal UPM.