Quantifiable Intravital Light Sheet Microscopy.

Live imaging of zebrafish embryos that maintains normal development can be difficult to achieve due to a combination of sample mounting, immobilization, and phototoxicity issues that, once overcome, often still results in image quality sufficiently poor that computer-aided analysis or even manual analysis is not possible. Here, we describe our mounting strategy for imaging the zebrafish midbrain-hindbrain boundary (MHB) with light sheet fluorescence microscopy (LSFM) and pilot experiments to create a study-specific set of parameters for semiautomatically tracking cellular movements in the embryonic midbrain primordium during zebrafish segmentation.

Navigating the Light-Sheet Image Analysis Software Landscape: Concepts for Driving Cohesion From Data Acquisition to Analysis.

From the combined perspective of biologists, microscope instrumentation developers, imaging core facility scientists, and high performance computing experts, we discuss the challenges faced when selecting imaging and analysis tools in the field of light-sheet microscopy. Our goal is to provide a contextual framework of basic computing concepts that cell and developmental biologists can refer to when mapping the peculiarities of different light-sheet data to specific existing computing environments and image analysis pipelines. We provide our perspective on efficient processes for tool selection and review current hardware and software commonly used in light-sheet image analysis, as well as discuss what ideal tools for the future may look like.

Building a three-dimensional model of early-stage zebrafish embryo brain.

We introduce a computational approach to build three-dimensional (3D) surface mesh models of the early-stage zebrafish brain primordia from time-series microscopy images. The complexity of the early-stage brain primordia and lack of recognizable landmarks pose a distinct challenge for feature segmentation and 3D modeling. Additional difficulty arises because of noise and variations in pixel intensity. We overcome these by using a hierarchical approach in which simple geometric elements, such as "beads" and "bonds," are assigned to represent local features and their connectivity is used to smoothen the surface while retaining high-curvature regions. We apply our method to build models of two zebrafish embryo phenotypes at discrete time points between 19 and 28 h post-fertilization and collect measurements to quantify development. Our approach is fast and applicable to building models of other biological systems, as demonstrated by models from magnetic resonance images of the human fetal brain. The source code, input scripts, sample image files, and generated outputs are publicly available on GitHub.

The Sterol Trafficking Pathway in Arabidopsis thaliana.

In plants, the trafficking mechanisms by which sterols move through the plant and into target cells are unknown. Earlier studies identified endosomes as primary candidates for internalization of sterols in plants, but these results have come into question. Here, we show that in elongating root cells, the internalization of sterol occurs primarily by a non-endocytic mechanism. Added fluorescent sterols [dehydroergosterol (DHE) and BODIPY-cholesterol (BCh)] do not initially label endosomes identified by fluorescent protein markers or by internalized FM4-64. Instead, the nuclear envelope, an organelle not associated with the endocytic pathway but part of the endoplasmic reticulum (ER), becomes labeled. This result is supported by experiments with the inducible overexpression of auxilin-2-like protein (AUX2 line), which blocks most endocytosis upon induction. Internalization and nuclear envelope labeling still occur in induced AUX2 cells. Longer-term incubation labels the oil body, a site involved in sterol storage. Although the first site of localization, the nuclear envelope, is part of the ER, other domains of the ER do not accumulate the label. The trafficking pathway differs from vesicular endocytosis and points toward a different pathway of sterol transport possibly involving other mechanisms, such as ER-plasma membrane contact sites and cytoplasmic transport.

Photothermal modulation of human stem cells using light-responsive 2D nanomaterials.

Two-dimensional (2D) molybdenum disulfide (MoS2) nanomaterials are an emerging class of biomaterials that are photoresponsive at near-infrared wavelengths (NIR). Here, we demonstrate the ability of 2D MoS2 to modulate cellular functions of human stem cells through photothermal mechanisms. The interaction of MoS2 and NIR stimulation of MoS2 with human stem cells is investigated using whole-transcriptome sequencing (RNA-seq). Global gene expression profile of stem cells reveals significant influence of MoS2 and NIR stimulation of MoS2 on integrins, cellular migration, and wound healing. The combination of MoS2 and NIR light may provide new approaches to regulate and direct these cellular functions for the purposes of regenerative medicine as well as cancer therapy.

Wnt signaling regulates neural plate patterning in distinct temporal phases with dynamic transcriptional outputs.

The process that partitions the nascent vertebrate central nervous system into forebrain, midbrain, hindbrain, and spinal cord after neural induction is of fundamental interest in developmental biology, and is known to be dependent on Wnt/ß-catenin signaling at multiple steps. Neural induction specifies neural ectoderm with forebrain character that is subsequently posteriorized by graded Wnt signaling: embryological and mutant analyses have shown that progressively higher levels of Wnt signaling induce progressively more posterior fates. However, the mechanistic link between Wnt signaling and the molecular subdivision of the neural ectoderm into distinct domains in the anteroposterior (AP) axis is still not clear. To better understand how Wnt mediates neural AP patterning, we performed a temporal dissection of neural patterning in response to manipulations of Wnt signaling in zebrafish. We show that Wnt-mediated neural patterning in zebrafish can be divided into three phases: (I) a primary AP patterning phase, which occurs during gastrulation, (II) a mes/r1 (mesencephalon-rhombomere 1) specification and refinement phase, which occurs immediately after gastrulation, and (III) a midbrain-hindbrain boundary (MHB) morphogenesis phase, which occurs during segmentation stages. A major outcome of these Wnt signaling phases is the specification of the major compartment divisions of the developing brain: first the MHB, then the diencephalic-mesencephalic boundary (DMB). The specification of these lineage divisions depends upon the dynamic changes of gene transcription in response to Wnt signaling, which we show primarily involves transcriptional repression or indirect activation. We show that otx2b is directly repressed by Wnt signaling during primary AP patterning, but becomes resistant to Wnt-mediated repression during late gastrulation. Also during late gastrulation, Wnt signaling becomes both necessary and sufficient for expression of wnt8b, en2a, and her5 in mes/r1. We suggest that the change in otx2b response to Wnt regulation enables a transition to the mes/r1 phase of Wnt-mediated patterning, as it ensures that Wnts expressed in the midbrain and MHB do not suppress midbrain identity, and consequently reinforce formation of the DMB. These findings integrate important temporal elements into our spatial understanding of Wnt-mediated neural patterning and may serve as an important basis for a better understanding of neural patterning defects that have implications in human health.

Analysis of the wnt1 regulatory chromosomal landscape.

One of the earliest patterning events in the vertebrate neural plate is the specification of mes/r1, the territory comprising the prospective mesencephalon and the first hindbrain rhombomere. Within mes/r1, an interface of gene expression defines the midbrain-hindbrain boundary (MHB), a lineage restriction that separates the mesencephalon and rhombencephalon. wnt1 is critical to mes/r1 development and functions within the MHB as a component of the MHB gene regulatory network (GRN). Despite its importance to these critical and early steps of vertebrate neurogenesis, little is known about the factors responsible for wnt1 transcriptional regulation. In the zebrafish, wnt1 and its neighboring paralog, wnt10b, are expressed in largely overlapping patterns, suggesting co-regulation. To understand wnt1 and wnt10b transcriptional control, we used a comparative genomics approach to identify relevant enhancers. We show that the wnt1-wnt10b locus contains multiple cis-regulatory elements that likely interact to generate the wnt1 and wnt10b expression patterns. Two of 11 conserved enhancers tested show activity restricted to the midbrain and MHB, an activity that is conserved in the distantly related spotted gar orthologous elements. Three non-conserved elements also play a likely role in wnt1 regulation. The identified enhancers display dynamic modes of chromatin accessibility, suggesting controlled deployment during embryogenesis. Our results suggest that the control of wnt1 and wnt10b expression is under complex regulation involving the interaction of multiple enhancers.

Decellularization of Whole Human Heart Inside a Pressurized Pouch in an Inverted Orientation.

The ultimate solution for patients with end-stage heart failure is organ transplant. But donor hearts are limited, immunosuppression is required, and ultimately rejection can occur. Creating a functional, autologous bio-artificial heart could solve these challenges. Biofabrication of a heart comprised of scaffold and cells is one option. A natural scaffold with tissue-specific composition as well as micro- and macro-architecture can be obtained by decellularizing hearts from humans or large animals such as pigs. Decellularization involves washing out cellular debris while preserving 3D extracellular matrix and vasculature and allowing "cellularization" at a later timepoint. Capitalizing on our novel finding that perfusion decellularization of complex organs is possible, we developed a more "physiological" method to decellularize non-transplantable human hearts by placing them inside a pressurized pouch, in an inverted orientation, under controlled pressure. The purpose of using a pressurized pouch is to create pressure gradients across the aortic valve to keep it closed and improve myocardial perfusion. Simultaneous assessment of flow dynamics and cellular debris removal during decellularization allowed us to monitor both fluid inflow and debris outflow, thereby generating a scaffold that can be used either for simple cardiac repair (e.g. as a patch or valve scaffold) or as a whole-organ scaffold.

Vascular Smooth Muscle Contractile Function Declines With Age in Skeletal Muscle Feed Arteries.

Aging induces a progressive decline in vasoconstrictor responses in central and peripheral arteries. This study investigated the hypothesis that vascular smooth muscle (VSM) contractile function declines with age in soleus muscle feed arteries (SFA). Contractile function of cannulated SFA isolated from young (4 months) and old (24 months) Fischer 344 rats was assessed by measuring constrictor responses of denuded (endothelium removed) SFA to norepinephrine (NE), phenylephrine (PE), and angiotensin II (Ang II). In addition, we investigated the role of RhoA signaling in modulation of VSM contractile function. Structural and functional characteristics of VSM cells were evaluated by fluorescence imaging and atomic force microscopy (AFM). Results indicated that constrictor responses to PE and Ang II were significantly impaired in old SFA, whereas constrictor responses to NE were preserved. In the presence of a Rho-kinase inhibitor (Y27632), constrictor responses to NE, Ang II, and PE were significantly reduced in young and old SFA. In addition, the age-group difference in constrictor responses to Ang II was eliminated. ROCK1 and ROCK2 content was similar in young and old VSM cells, whereas pROCK1 and pROCK2 were significantly elevated in old VSM cells. Aging was associated with a reduction in smooth muscle α-actin stress fibers and recruitment of proteins to cell-matrix adhesions. Old VSM cells presented an increase in integrin adhesion to the matrix and smooth muscle γ-actin fibers that was associated with increased cell stiffness. In conclusion, our results indicate that VSM contractile function declined with age in SFA. The decrement in contractile function was mediated in part by RhoA/ROCK signaling. Upregulation of pROCK in old VSM cells was not able to rescue contractility in old SFA. Collectively, these results indicate that changes at the VSM cell level play a central role in the reduced contractile function of aged SFA.

Nck deficiency is associated with delayed breast carcinoma progression and reduced metastasis.

Although it is known that noncatalytic region of tyrosine kinase (Nck) regulates cell adhesion and migration by bridging tyrosine phosphorylation with cytoskeletal remodeling, the role of Nck in tumorigenesis and metastasis has remained undetermined. Here we report that Nck is required for the growth and vascularization of primary tumors and lung metastases in a breast cancer xenograft model as well as extravasation following injection of carcinoma cells into the tail vein. We provide evidence that Nck directs the polarization of cell-matrix interactions for efficient migration in three-dimensional microenvironments. We show that Nck advances breast carcinoma cell invasion by regulating actin dynamics at invadopodia and enhancing focalized extracellular matrix proteolysis by directing the delivery and accumulation of MMP14 at the cell surface. We find that Nck-dependent cytoskeletal changes are mechanistically linked to enhanced RhoA but restricted spatiotemporal activation of Cdc42. Using a combination of protein silencing and forced expression of wild-type/constitutively active variants, we provide evidence that Nck is an upstream regulator of RhoA-dependent, MMP14-mediated breast carcinoma cell invasion. By identifying Nck as an important driver of breast carcinoma progression and metastasis, these results lay the groundwork for future studies assessing the therapeutic potential of targeting Nck in aggressive cancers.

Midbrain-Hindbrain Boundary Morphogenesis: At the Intersection of Wnt and Fgf Signaling.

A constriction in the neural tube at the junction of the midbrain and hindbrain is a conserved feature of vertebrate embryos. The constriction is a defining feature of the midbrain-hindbrain boundary (MHB), a signaling center that patterns the adjacent midbrain and rostral hindbrain and forms at the junction of two gene expression domains in the early neural plate: an anterior otx2/wnt1 positive domain and a posterior gbx/fgf8 positive domain. otx2 and gbx genes encode mutually repressive transcription factors that create a lineage restriction boundary at their expression interface. Wnt and Fgf genes form a mutually dependent feedback system that maintains their expression domains on the otx2 or gbx side of the boundary, respectively. Constriction morphogenesis occurs after these conserved gene expression domains are established and while their mutual interactions maintain their expression pattern; consequently, mutant studies in zebrafish have led to the suggestion that constriction morphogenesis should be considered a unique phase of MHB development. We analyzed MHB morphogenesis in fgf8 loss of function zebrafish embryos using a reporter driven by the conserved wnt1 enhancer to visualize anterior boundary cells. We found that fgf8 loss of function results in a re-activation of wnt1 reporter expression posterior to the boundary simultaneous with an inactivation of the wnt1 reporter in the anterior boundary cells, and that these events correlate with relaxation of the boundary constriction. In consideration of other results that correlate the boundary constriction with Wnt and Fgf expression, we propose that the maintenance of an active Wnt-Fgf feedback loop is a key factor in driving the morphogenesis of the MHB constriction.

Inverted orientation improves decellularization of whole porcine hearts.

In structurally heterogeneous organs, such as heart, it is challenging to retain extracellular matrix integrity in the thinnest regions (eg, valves) during perfusion decellularization and completely remove cellular debris from thicker areas. The high inflow rates necessary to maintain physiologic pressure can distend or damage thin tissues, but lower pressures prolong the process and increase the likelihood of contamination. We examined two novel retrograde decellularization methods for porcine hearts: inverting the heart or venting the apex to decrease inflow rate. We measured flow dynamics through the aorta (Ao) and pulmonary artery (PA) at different Ao pressures and assessed the heart's appearance, turbidity of the outflow solutions, and coronary perfusion efficiency. We used rectangle image fitting of decellularized heart images to obtain a heart shape index. Using nonlinear optical microscopy, we determined the microstructure of collagen and elastin fibers of the aortic valve cusps. DNA, glycosaminoglycan, and residual detergent levels were compared. The inverted method was superior to the vented method, as shown by a higher coronary perfusion efficiency, more cell debris outflow, higher collagen and elastin content inside the aortic valve, lower DNA content, and better retention of the heart shape after decellularization. To our knowledge, this is the first study to use flow dynamics in a whole heart throughout the decellularization procedure to provide real-time information about the success of the process and the integrity of the vulnerable regions of the matrix. Heart orientation was important in optimizing decellularization efficiency and maintaining extracellular matrix integrity. STATEMENT OF SIGNIFICANCE: The use of decellularized tissue as a suitable scaffold for engineered tissue has emerged over the past decade as one of the most promising biofabrication platforms. The decellularization process removes all native cells, leaving the natural biopolymers, extracellular matrix materials and native architecture intact. This manuscript describes heart orientation as important in optimizing decellularization efficiency and maintaining extracellular matrix integrity. To our knowledge, this is the first study to assess flow dynamics in a whole heart throughout the decellularization procedure. Our findings compared to currently published methods demonstrate that continuous complex real-time measurements and analyses are required to produce an optimal scaffold for cardiac regeneration.

Combined lineage mapping and gene expression profiling of embryonic brain patterning using ultrashort pulse microscopy and image registration.

During embryogenesis, presumptive brain compartments are patterned by dynamic networks of gene expression. The spatiotemporal dynamics of these networks, however, have not been characterized with sufficient resolution for us to understand the regulatory logic resulting in morphogenetic cellular behaviors that give the brain its shape. We have developed a new, integrated approach using ultrashort pulse microscopy [a high-resolution, two-photon fluorescence (2PF)-optical coherence microscopy (OCM) platform using 10-fs pulses] and image registration to study brain patterning and morphogenesis in zebrafish embryos. As a demonstration, we used time-lapse 2PF to capture midbrain-hindbrain boundary morphogenesis and a wnt1 lineage map from embryos during brain segmentation. We then performed in situ hybridization to deposit NBT/BCIP, where wnt1 remained actively expressed, and reimaged the embryos with combined 2PF-OCM. When we merged these datasets using morphological landmark registration, we found that the mechanism of boundary formation differs along the dorsoventral axis. Dorsally, boundary sharpening is dominated by changes in gene expression, while ventrally, sharpening may be accomplished by lineage sorting. We conclude that the integrated visualization of lineage reporter and gene expression domains simultaneously with brain morphology will be useful for understanding how changes in gene expression give rise to proper brain compartmentalization and structure.

Sequential multimodal microscopic imaging and biaxial mechanical testing of living multicomponent tissue constructs.

Understanding relationships between mechanical stimuli and cellular responses require measurements of evolving tissue structure and mechanical properties. We developed a 3D tissue bioreactor that couples to both the stage of a custom multimodal microscopy system and a biaxial mechanical testing platform. Time dependent changes in microstructure and mechanical properties of fibroblast seeded cruciform fibrin gels were investigated while cultured under either anchored (1.0:1.0 stretch ratio) or strip biaxial (1.0:1.1) conditions. A multimodal nonlinear optical microscopy-optical coherence microscopy (NLOM-OCM) system was used to delineate noninvasively the relative spatial distributions of original fibrin, deposited collagen, and fibroblasts during month long culture. Serial in-culture mechanical testing was also performed to track the evolution of bulk mechanical properties under sterile conditions. Over the month long time course, seeded cells and deposited collagen were randomly distributed in equibiaxially anchored constructs, but exhibited preferential alignment parallel to the direction of the 10% stretch in constructs cultured under strip biaxial stretch. Surprisingly, both anchored and strip biaxial stretched constructs exhibited isotropic mechanical properties (including progressively increasing stiffness) despite developing a very different collagen microstructural organization. In summary, our biaxial bioreactor system integrating both NLOM-OCM and mechanical testing provided complementary information on microstructural organization and mechanical properties and, thus, may enable greater fundamental understanding of relationships between engineered soft tissue mechanics and mechanobiology.

Dynamic multicomponent engineered tissue reorganization and matrix deposition measured with an integrated nonlinear optical microscopy-optical coherence microscopy system.

Multicomponent tissue models are viable tools to better understand cell responses in complex environments, but present challenges when investigated with live cell microscopy noninvasively. In this study, integrated nonlinear optical microscopy-optical coherence microscopy (NLOM-OCM) was used to characterize cell interactions within three-dimensional (3-D), multicomponent extracellular matrices. In fibrin-collagen mixtures, 3T3 fibroblasts were observed to recruit both fibrin and collagen fibers while remodeling matrices. Also, NLOM-OCM was used to observe collagen deposition by neonatal human dermal fibroblasts within originally fibrin matrices over an extended time. It was observed that preferentially aligned collagen deposition could be achieved with aligned fibroblasts but that cell alignment could be achieved without aligning the extant extracellular matrix. In summary, this multimodel imaging system has potential for both real-time and longitudinal imaging of living 3-D cultures, which is particularly important for evaluating cell microenvironments in composite scaffolds or serial characterization of engineered tissue constructs during culture.

Time-lapse ultrashort pulse microscopy of infection in three-dimensional versus two-dimensional culture environments reveals enhanced extra-chromosomal virus replication compartment formation.

The mechanisms that enable viruses to harness cellular machinery for their own survival are primarily studied in cell lines cultured in two-dimensional (2-D) environments. However, there are increasing reports of biological differences between cells cultured in 2-D versus three-dimensional (3-D) environments. Here we report differences in host-virus interactions based on differences in culture environment. Using ultrashort pulse microscopy (UPM), a form of two-photon microscopy that utilizes sub-10-fs pulses to efficiently excite fluorophores, we have shown that de novo development of extra-chromosomal virus replication compartments (VRCs) upon murine cytomegalovirus (mCMV) infection is markedly enhanced when host cells are cultured in 3-D collagen gels versus 2-D monolayers. In addition, time-lapse imaging revealed that mCMV-induced VRCs have the capacity to grow by coalescence. This work supports the future potential of 3-D culture as a useful bridge between traditional monolayer cultures and animal models to study host-virus interactions in a more physiologically relevant environment for the development of effective anti-viral therapeutics. These advances will require broader adoption of modalities, such as UPM, to image deep within scattering tissues.

Passive pressure-diameter relationship and structural composition of rat mesenteric lymphangions.

BACKGROUND: Lymph flow depends on both the rate of lymph production by tissues and the extent of passive and active pumping. Here we aim to characterize the passive mechanical properties of a lymphangion in both mid-lymphangion and valve segments to assess regional differences along a lymphangion, as well as evaluating its structural composition. METHODS AND RESULTS: Mesenteric lymphatic vessels were isolated and cannulated in a microchamber for pressure-diameter (P-D) testing. Vessels were inflated from 0 to 20 cmH(2)O at a rate of 4 cmH(2)O/min, and vessel diameter was continuously tracked, using an inverted microscope, video camera, and custom LabVIEW program, at both mid-lymphangion and valve segments. Isolated lymphatic vessels were also pressure-fixed at 2 and 7 cmH(2)O and imaged using a nonlinear optical microscope (NLOM) to obtain collagen and elastin structural information. We observed a highly nonlinear P-D response at low pressures (3-5 cmH(2)O), which was modeled using a three-parameter constitutive equation. No significant difference in the passive P-D response was observed between mid-lymphangion and valve regions. NLOM imaging revealed an inner elastin layer and outer collagen layer at all locations. Lymphatic valve leaflets were predominantly elastin with thick axially oriented collagen bands at the insertion points. CONCLUSIONS: We observed a highly nonlinear P-D response at low pressures (3-5 cmH(2)O) and developed the first constitutive equation to describe the passive P-D response for a lymphangion. The passive P-D response did not vary among regions, in agreement with the composition of elastin and collagen in the lymphatic wall.

Two-photon excited fluorescence enhancement with broadband versus tunable femtosecond laser pulse excitation.

The inverse relationship between two-photon excited fluorescence (TPEF) and laser pulse duration suggests that two-photon microscopy (TPM) performance may be improved by decreasing pulse duration. However, for ultrashort pulses of sub-10 femtosecond (fs) in duration, its spectrum contains the effective gain bandwidth of Ti:Sapphire and its central wavelength is no longer tunable. An experimental study was performed to explore this apparent tradeoff between untuned sub-10 fs transform-limited pulse (TLP) and tunable 140 fs pulse for TPEF. Enhancement factors of 1.6, 6.7, and 5.2 are measured for Indo-1, FITC, and TRITC excited by sub-10 fs TLP compared with 140 fs pulse tuned to the two-photon excitation (TPE) maxima at 730 nm, 800 nm, and 840 nm, respectively. Both degenerate (v(1) = v(2)) and nondegenerate (v(1) ≠ v(2)) mixing of sub-10 fs TLP spectral components result in its broad second-harmonic (SH) power spectrum and high spectral density, which can effectively compensate for the lack of central wavelength tuning and lead to large overlap with dye TPE spectra for TPEF enhancements. These pulse properties were also exploited for demonstrating its potential applications in multicolor imaging with TPM.

Mechanical assessment of elastin integrity in fibrillin-1-deficient carotid arteries: implications for Marfan syndrome.

AIMS: Elastin is the primary component of elastic fibres in arteries, which contribute significantly to the structural integrity of the wall. Fibrillin-1 is a microfibrillar glycoprotein that appears to stabilize elastic fibres mechanically and thereby to delay a fatigue-induced loss of function due to long-term repetitive loading. Whereas prior studies have addressed some aspects of ageing-related changes in the overall mechanical properties of arteries in mouse models of Marfan syndrome, we sought to assess for the first time the load-carrying capability of the elastic fibres early in maturity, prior to the development of ageing-related effects, dilatation, or dissection. METHODS AND RESULTS: We used elastase to degrade elastin in common carotid arteries excised, at 7-9 weeks of age, from a mouse model (mgR/mgR) of Marfan syndrome that expresses fibrillin-1 at 15-25% of normal levels. In vitro biaxial mechanical tests performed before and after exposure to elastase suggested that the elastic fibres exhibited a nearly normal load-bearing capability. Observations from nonlinear optical microscopy suggested further that competent elastic fibres not only contribute to load-bearing, they also increase the undulation of collagen fibres, which endows the normal arterial wall with a more compliant response to pressurization. CONCLUSION: These findings support the hypothesis that it is an accelerated fatigue-induced damage to or protease-related degradation of initially competent elastic fibres that render arteries in Marfan syndrome increasingly susceptible to dilatation, dissection, and rupture.

Cornea microstructure and mechanical responses measured with nonlinear optical and optical coherence microscopy using sub-10-fs pulses.

A combined nonlinear optical microscopy (NLOM) and optical coherence microscopy (OCM) imaging system has been assembled in order to simultaneously capture co-registered volumetric images of corneal morphology and biochemistry. Tracking of cell nuclei visible in the OCM volume enabled the calculation of strain depth profile in response to changes in intraocular pressure for rabbit cornea stroma. Results revealed nonlinear responses with a depth dependent strain distribution, exhibiting smaller strains in the anterior and larger strains in the posterior stroma. Cross-sectional images of collagen lamellae, visible in NLOM, showed inhomogeneous collagen structure along the anterior-posterior direction that correlated well with the noted heterogeneous corneal mechanical responses.