Needle-shaped beam optical coherence tomography angiography visualizing cortical blood perfusion after photothrombotic stroke in vivo

doi:10.21203/rs.3.rs-3087686/v1

Download PDF

Article

Needle-shaped beam optical coherence tomography angiography visualizing cortical blood perfusion after photothrombotic stroke in vivo

https://doi.org/10.21203/rs.3.rs-3087686/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Optical imaging techniques provide low-cost, non-radiative images with high spatiotemporal resolution, making them advantageous for long-term dynamic observation of blood perfusion in stroke research and other brain studies compared to non-optical methods. However, high-resolution imaging in optical microscopy fundamentally requires a tight optical focus, and thus a limited depth of field (DOF). Consequently, large-scale, non-stitched, high-resolution images of curved surfaces, like brains, are difficult to acquire without z-axis scanning. To overcome this limitation, we developed a needle-shaped beam optical coherence tomography angiography (NB-OCTA) system, and for the first time, achieved a volumetric resolution of less than 8 µm in a non-stitched volume space of 6.4 mm × 4 mm × 620 µm in vivo. This system captures the distribution of blood vessels at 3.4-times larger depths than normal OCTA equipped with a Gaussian beam (GB-OCTA). We then employed NB-OCTA to perform long-term observation of cortical blood perfusion after stroke, and quantitatively analyzed the vessel area density (VAD) and the diameters of representative vessels in different regions over 10 days, revealing different spatiotemporal dynamics in the acute, sub-acute and chronic phase of post-ischemic revascularization. Benefiting from our NB-OCTA, we revealed that the recovery process is not only the result of spontaneous reperfusion, but also the formation of new vessels. This study provides visual and mechanistic insights into strokes and helps to deepen our understanding of the spontaneous response of brain after stroke.

Physical sciences/Optics and photonics/Optical techniques/Imaging and sensing

Physical sciences/Optics and photonics/Other photonics/Biophotonics

Stroke is the second highest cause of death globally and a leading cause of disability. It has shown increasing incidence in developing countries, where 70% strokes are of ischemic origin¹. Due to its widespread prevalence and limited therapeutic management options, ongoing stroke research revolves around investigation of several therapeutic interventions in preclinical and clinical settings. Animal models have been extensively used in preclinical research^2,3, and three commonly used methods to model experimental stroke include craniectomy models^4,5, intraluminal suture^6,7, and photothrombotic (PT) stroke model^8,9. Craniectomy models involve invasive clipping, ligation or electrocoagulation techniques that can cause mechanical damage to the cortex, and intraluminal suture models require a high-level of surgical expertise and may not generate reproducible ischemic lesions. Contrastingly, the PT stroke model has gained wide acceptance due to its precise special controllability, which allows targeting of any cortical region of interest in a reproducible and minimally invasive manner¹⁰.

Although the research on the pathological mechanisms of ischemic stroke has made progress^11,12, clinical treatment options and drug development still remain challenging. Quick and effective restoration of cerebral blood flow can arrest pathological events and reduce functional damage¹³. Therefore, thrombolysis is critical in early treatment of ischemic stroke^14,15, with an optimal window of 4.5 hours after ischemia¹⁶, but most patients do not catch this time window. In addition to natural and interventional revascularization, localized spontaneous vascular regeneration, or angiogenesis can occur. Angiogenesis is triggered by hypoxia as oxygen deficiency upregulates proangiogenic factors¹⁷. This results in the subsequent upregulation of hypoxia-inducible factors, followed by the expression of angiopoietin, erythropoietin, nitric oxide synthase, vascular endothelial growth factor (VEGF) and their receptors¹⁸. Angiogenesis is also believed to occur in the adult brain, but this process remains largely unexplored^19,20. Hence, a rapid, high-resolution, in vivo imaging technique is a vital asset in both visually assessing blood perfusion and brain tissue damage, and longitudinal monitoring of the ischemic area to study the mechanisms of post-stroke revascularization.

Magnetic resonance imaging (MRI) and computed tomography (CT) are commonly used brain imaging modalities in clinical practice ^21,22. MRI has a high sensitivity and specificity for detecting pathological changes such as cerebral vessels, lesions, and tumors, but repeated MRI scans can be outrageously expensive and time consuming. CT provides fast and high-quality images, but frequent measurements involve high levels of radiation exposure. While these non-optical imaging methods can achieve a large imaging depth, their resolution is limited to the sub-millimeter level.

Resultingly, optical imaging techniques are becoming a popular choice for brain studies due to their high spatial and temporal resolution^10,23,24. Researchers have demonstrated an increasing impact on brain imaging by using photoacoustic microscopy (PAM)²³, which can measure blood flow and oxygen metabolism by detecting endogenous and exogenous contrasts ²⁵. However, PAM requires water as the coupling medium for ultrasound transmission, making the system huge and complicated. Two-photon (2P) microscopy, allows subcellular image resolution and imaging depths of several hundred micrometers into living brain tissue²⁶, but struggles to assess large-scale blood flow in typical stroke models due to its slow three-dimensional (3D) scan speeds and limited field of view (FOV). Furthermore, 2P microscopy observations of stroke models may permit dye leakage thus impeding microvasculature imaging due to ischemic damage to the cortical tissue²⁷. It would therefore be particularly useful for stroke studies to have intrinsic contrast blood flow signals without introducing exogenous agents into the system. Laser speckle contrast imaging (LSCI) is a label-free technique that utilizes dynamic light scattering to visualize blood flow²⁸, providing measurements of relative flow velocity across a large FOV. However, LSCI suffers from limited spatial resolution, being limited to two-dimensional imaging and lacking depth-resolved information.

Optical coherence tomography (OCT)²⁹ is an emerging optical imaging modality that enables three-dimensional (3D) volumetric imaging of biological tissue microstructure. Additionally, because OCT is a non-invasive, contrast agent and radiation free technique, it shows high promise across a wide range of medical and biological applications. OCT angiography (OCTA)³⁰ is the application of OCT for measuring and visualizing 3D mapping of blood perfusion, achieved by mathematically analyzing the red blood cells’ motion-induced temporal changes of scattering signals. Compared with other optical techniques, OCTA allows rapid, high-sensitivity, contrast-free imaging with micron-scale resolution, making it a promising choice for large field of view applications like stroke monitoring¹⁰.

However, for high-resolution optical imaging, the pursuit of diffraction-limited resolution is often accompanied by a trade-off in the form of a reduced DOF, resulting in a degradation of lateral resolution beyond the optical focal plane. This constraint impedes rapid acquisition of high-resolution images of samples with irregular surfaces. Imaging such tissues usually involves laborious axial scanning at multiple planes and complicated image processing procedures. Furthermore, for time-sensitive applications like intraoperative histology or cerebral hemodynamics, rapid high-resolution imaging of irregular surfaces along a large DOF is incredibly desirable. Consequently, numerous approaches have been attempted to enable high-resolution, large DOF imaging. Multi-beam structures³¹ can achieve a larger DOF by simultaneously focusing on different depths of the sample, but require complex hardware and software implementations. Additionally, non-diffracting beams such as Bessel^32,33 or Airy³⁴ beams can achieve large DOF and high resolution, but the imaging quality is usually limited by the severe sidelobes and low beam efficiency. In recent years, deep learning models have been used to improve the DOF^35,36, but this requires a large amount of ground truth during training, and obtaining such data may be challenging or infeasible for certain subjects.

To improve 3D resolution in an optical microscope, we developed a needle-shaped beam (NB), demonstrating its superiority in OCT anatomical imaging of both human skin and drosophila larva ³⁷, and furthermore, successfully applied it to PAM systems ²³. In this work, we developed an NB-OCTA system that, for the first time, achieved a volumetric resolution of less than 8 µm in a non-stitched volume space of 6.4 mm × 4 mm × 620 µm. Compared to normal OCTA using Gaussian beam (GB-OCTA), NB-OCTA captures the distribution of blood vessels within a 3.4-times larger depth range. We then employed NB-OCTA to perform long-term in vivo observation of cortical blood perfusion after PT stroke, and quantitatively analyzed the vessel area density (VAD) and the diameters of representative vessels in different regions over 10 days, revealing the dynamic spatial and temporal evolution of post-thrombotic blood flow perfusion. Benefiting from our NB-OCTA, we revealed that the recovery process is not only the result of spontaneous reperfusion, but also the formation of new vessels. This study provides a powerful monitoring tool for therapeutic investigations and drug screening for strokes, and deepens our understanding of the mechanisms of ischemic injury and subsequent revascularization in stroke.

The principle of NB-OCTA

An NB is composed of numerous closely adjacent foci along the axial direction, which are generated by a customized diffractive optical element (DOE) and a corresponding objective³⁷. The beam length can be flexibly adjusted by increasing or decreasing the number of foci. The DOE phase is formulated as

$$\begin{array}{c}{\text{P}}_{\text{DOE}}\left(\text{x}\text{,}\text{ }\text{y}\right)\text{ }\text{=}{\sum }_{\text{m}\text{=1}}^{\text{M}}\left\{\left[\text{-πn}\text{(}{\text{x}}^{\text{2}}\text{+}{\text{y}}^{\text{2}}\text{)(1}\text{ }\text{/}{\text{ }\text{f}}_{\text{m}}\text{-}\text{1}\text{ }\text{/}\text{ }\text{f}\text{ }\text{)}\text{ }\text{/}\text{ }\text{λ}\text{ }\text{-}\text{ }\text{2}\text{π}\text{m}\text{A}\right]{\text{Loc}}_{\text{m}}\text{(}\text{x}\text{,}\text{ }\text{y}\text{)}\right\}\#\left(1\right)\end{array}$$

where (x, y) is the planar coordinate, λ is the wavelength, f is the focal length of the objective, n is the refractive index of the surrounding medium, M is the foci number, m is the focus index, Loc_m (x, y) is a binary matrix whose value is either 0 or 1 to identify whether the pixel is allocated to the focus f_m, 2πmA is the phase regulator for adjusting the beam diameter, and the item $\left[\text{-πn}\text{(}{\text{x}}^{\text{2}}\text{+}{\text{y}}^{\text{2}}\text{)(1}\text{ /}{\text{ f}}_{\text{m}}\text{-}\text{1 }\text{/ f }\text{) }\text{/ λ - }\text{2}\text{π}\text{m}\text{A}\right]$ aims to shift the focus from f to f_m. The DOE pixels are equally and randomly allocated into m subsets (Fig. 1a). The rule for determining the number of foci is to control the average interval between two adjacent foci within one Rayleigh length. f₁ was set to be coincident with the objective lens focus f, the beam length is determined by (f_m– f₁). In this work, we designed two NBs for different scenarios, one being 420 µm long with a diameter between 6.5 ± 0.5 µm and the second being 620 µm long with a diameter between 7.4 ± 0.6 µm. The first NB is composed of 49 foci with a phase modulator A of 0.01, while the second contains 81 foci with A = 0.008.

The profiles of Gaussian beam and 620 µm NB are compared in Fig. 1b. We define the beam length as the range in which the spot size is smaller than 8 µm. The Gaussian beam has a 4 µm (full-width-at-half-maximum, FWHM) diameter size at the focal plane, expanding to 8 µm in the axial range of 180 µm. The NB achieves a minimal spot size of 6.8 µm, remaining smaller than 8 µm across a depth range of 620 µm. The NB's consistent spot size is demonstrated across three axial positions (z = -300 µm, z = 0 µm, z = 300 µm). Note the NB is tapered, with narrower edges and a thicker center. Contrastingly, the Gaussian beam has the smallest spot size centered at the focal plane. With a DOE, we achieved a maximum efficiency of 30%, which means that the energy enclosed in the main lobe is 30% of the GB focal spot energy. Compared to Bessel beam whose sidelobe-to-main-lobe ratio can be up to 20%, our NB can be optimized to be less than 3%, eliminating the need for complicated image processing in NB-OCTA. The DOE was easily installed in the sample arm before the galvo mirror and objective lens (Fig. 1c), without requiring modifications to the original system (details shown in Fig. S1).

To acquire OCTA signals, intensity information from consecutive B-scans at the same spatial position is compared. Static tissue produces consistent intensity, while the movement of red blood cells through vessels produces intensity variations. By subtracting the intensity values of the consecutive B-scans, the OCTA algorithm selectively enhances the changes caused by blood flow while suppressing the static signals.

NB-OCTA system performance

The focused Gaussian beam and two NBs were tested with 5 µm polystyrene (PS) beads (XFJ100, XFnano) embedded in ultrasound gel (Ultrasonic Coupling, Cofoe). The phantom with beads concentration around 2×10⁵ /mm³ was degassed with a centrifuge (10min at 15 000 rpm, H1850, Cence). The B-scan image acquired by the GB only showed recognizable beads only near the focal plane, while areas far from the focus displayed significant blurring, Fig. 2a. For 420µm NB, particles can be distinguished across a wider range, while the 620 µm NB could distinguish particles wider still. Figure. 2b shows x-y planes at five depths. At z = 0, all three beams can distinguish beads clearly. For 620 µm NB, a 21% sidelobe ratio appears, but in the dense bio-samples, the sidelobe effect can be limited, as demonstrated in the following experiments. At z= ± 210 µm, Gaussian beam beads profiles deform significantly, while 420 µm and 620 µm NBs maintain clarity. At z= ± 310 µm, only the 620 µm NB provides high-quality images. Complete data of large scanning area and at more given depth can be found in Supplementary information Fig. S3-S4. The contrast-to-noise ratios (CNRs) along the depth and the signal-to-noise ratios (SNRs) along the depth are shown in Fig. 2c and 2d. Here, CNR= (average signal intensity- average background intensity)/standard deviation of background intensity, SNR = average signal intensity/ standard deviation of background intensity. We demonstrate the contrast and sensitivity of NB stay relatively stable within the designed beam lengths, and furthermore, NBs outperform the focused Gaussian beam in the depths far from the focal plane. The beam diameter and relative intensity shown in Fig. 2e and 2f were measured by a microscope system together with a beam profiler (BP209IR1, Thorlabs) with 10 µm as the axial scanning step size (experimental setup shown in Fig. S2).

OCTA has been widely used for retinal and cerebral vascular imaging. However, the large curvature of the mouse brain cortex makes it challenging for conventional GBs to cover the entire region due to its limited DOF. Consequently, generating a large-field brain vascular image^10,38 has previously required stitching multiple small FOVs, which fundamentally prohibits simultaneous imaging of the entire cortical vasculature, and introduces significant complexity and uncertainty in FOV merging during image processing. Using NB-OCTA, we collected in vivo, non-stitched, large field of view (6.4 mm × 4 mm) images of the mouse brain (Fig. 2g and 2h). To better visualize the brain’s curvature, we encoded the depth with colors. Conventional GB-OCTA could image the small vessels within its focal range, but missed many vessels in the surrounding lower periphery. In contrast, NB-OCTA clearly resolved more vessels around the margins and could image the entire 800 µm depth range. The brain was scanned five times using 58.5 s, with a fast-axis scanning step of 8 µm, and a slow-axis scanning step of 5 µm.

Quantitative comparison between GB-OCTA and NB-OCTA

Before monitoring cerebral stroke, we imaged the selected area using both GB-OCTA and NB-OCTA. This region, situated on the left side of the mouse brain's sagittal suture, exhibited a depth difference of 450 µm between the left and right sides due to the brain's natural curvature. This height difference exceeded the DOF of most high-resolution OCTA systems. Traditional solution is to rotate the sample or the probe to align the Gaussian beam perpendicular to one subregion of sample surface³⁸, capture image of that region, and then rotate to the next. This process is repeated multiple times and the images are stitched together to minimize defocusing. However, this method is time-consuming and requires software and hardware coordination, thus being inconvenient during rapid imaging. In contrast, our proposed NB has a focal depth larger than the height difference of the sample surface, enabling direct high-resolution imaging of the selected region without the need for rotation. Figure 3a and Fig. 3b show OCTA images obtained with GB and NB, respectively. Comparing these images, the GB-OCTA image exhibits blurriness on the left and right sides due to defocusing, while the NB-OCTA image maintains high resolution. We zoomed in on three selected regions at depths of z=-200 µm, 0 µm, and 200 µm, Fig. 3c, 3d, and 3e. Profiles of selected single vessels (indicated by blue arrows) were drawn in Fig. 3f, 3g, and 3h, respectively. Outside the DOF of GB (yellow and green boxes), vessels in the GB-OCTA image appear widened, whereas the NB-OCTA image remains unaffected. In the central region (red boxes), both GB-OCTA and NB-OCTA were able to provide a high-resolution image, Fig. 3d, and the selected blood vessel diameter was measured to be the same in both modalities, Fig. 3g. However, as light penetrates deeper into the tissue, the Gaussian beam apparently defocused, resulting in weaker light intensity and larger beam size at the current scanning point, which increased noise. Consequently, some micro blood vessels can only be visualized in NB-OCTA, red arrows in Fig. 3d.

In addition, we also calculated CNR, SNR, and vessel connectivity to quantitatively evaluate the image quality of Fig. 3a and Fig. 3b. In Fig. 3a, the CNR, SNR, and vessel connectivity were 16.8 dB, 17.7 dB, and 54.8%, respectively, while in Fig. 3b, they were 20.6 dB, 22.1 dB, and 73.7%, respectively. The connectivity is defined as the ratio of the biggest number of connected pixels and the total number of pixels that refers to vessels. Higher CNR helps to accurately distinguish blood vessels from surrounding tissues, while higher SNR improves image clarity and visibility of details, aiding in the accurate assessment of vessel connectivity.

To further demonstrate that the GB is not suitable for this scenario, we present experimental results of focusing at multiple z-positions, see Fig. S5-S6 in Supplementary. The results show that no matter where the focus is, the DOF of the GB cannot cover the height difference of the sample surface, so there would always be positions out of focus. Based on the experimental results above, NB-OCTA was employed to monitor the whole process of the recovery from PT stroke.

Long-term observation of blood perfusion after stroke with NB-OCTA

NB-OCTA can overcome the natural curvature of the mouse brain and achieve fast high-resolution imaging, making it an ideal tool for monitoring blood thrombosis. Figure 4 illustrates the entire process in mouse cerebral vasculature from acute ischemia to chronic recovery following PT stroke induction. The first two images illustrate a baseline prior to PT induction with GB and NB. After the PT stroke was induced, a focal ischemic area was formed. On the first day after PT induction (30 minutes post-induction), partial loss of blood flow signals was observed. In the following two days, the ischemic area expanded significantly, with the vascular network in the core region disappearing remarkably. Around the fourth day, the vascular network began to recover and the ischemic region began to shrink. Some blood vessels that had disappeared due to blockage reappeared due to spontaneous recanalization. In addition, the diameters of large blood vessels in the core area (indicated by red arrows) expanded significantly, possibly to compensate for other blocked vessels. On the fifth and sixth days, the area of the core ischemic region continued to shrink, and more blocked small blood vessels spontaneously reperfused. Starting from the seventh day the ischemic area tended to stabilize, and the expanded blood vessels due to compensation effects also began to regress. However, compared with the baseline image before the thrombosis, some deep capillaries did not reperfuse.

In addition, we quantitatively analyzed the entire process as shown in Fig. 5. The core ischemic area, denoted by the red circle, encapsulated the vast majority of the ischemic tissue, while the area outside the red circle was defined as the peripheral area. We conducted a long-term analysis of VAD in two regions, and, as shown in Fig. 5b, the VAD in the core area dropped sharply at first, then began to recover on the fourth day and stabilized after the seventh day. The VAD of the peripheral area was relatively stable with small fluctuations, which was consistent with the intuitive results of Fig. 4. We also selected three representative blood vessels of three groups to analyze the diameter changes during the whole process, shown in Fig. 5b. Group A is composed of vessels from the core area that did not disappear or were not blocked. The diameter of these vessels changed over the course of the experiment. Vessel A2 was partially blocked and then recovered. Vessel A1 and A3 got widened, which indicates they did more transportation work to compensate for the blocked vessels. Following the reperfusion of these blocked vessels, A1 and A3 returned back to their normal state. Group B consists of vessels that were blocked and disappeared from the core, but then reperfused spontaneously. Group C is comprised of the vessels from the periphery. These vessels were relatively stable, and do not have obvious diameter changes, indicating that the peripheral vessels were not significantly affected. This further supports the advantages of the PT model in terms of its precision and controllability.

Additionally, by comparing the spatial location of the reperfused flow with the original vessels (Fig. 5c and 5d), we believe that the recovery process is not only the result of spontaneous reperfusion (indicated by the red lines), but also the formation of new vessels, as indicated by the blue oval area in Fig. 5d. While currently the difference is manually identified, more compelling evidence about angiogenesis can be obtained through molecular experiments such as the Western blotting for angiogenic markers such as VEGF and its fms-like tyrosine (FLT) kinase receptor³⁹.

In this work we have designed and demonstrated a novel NB-OCTA system featuring a substantially improved field of view and DOF compared to conventional GB-OCTA systems. This system enables, for the first time, a volumetric resolution of 8 µm across a non-stitched volume space of 6.4 mm × 4 mm × 620 µm. The NB's elongated DOF enables rapid high-resolution imaging across the brain's curved surface, something previously impossible without time-consuming Z-stacking and laborious image compiling. Note that the current improvements of DOF are not the theoretical limit, which can be further adjusted according to the laser energy and required efficiency of the system. In addition, we demonstrated 10 days observation of blood perfusion after stroke with NB-OCTA, revealing distinct spatiotemporal dynamics in the acute, sub-acute and chronic phase of post-ischemic revascularization. Our experiment displays the initial three days following PT occlusion can be considered as the acute phase of vessel thrombosis, during which the VAD decreased gradually, resulting in peak ischemia area on day 3. Subsequently, the blood flow began to recover from the peripheral region towards the ischemic core, and from the seventh day onwards, the ischemic area tended to stabilize. Furthermore, the NB-OCTA has demonstrated the recovery process is not only the product of spontaneous reperfusion, but also the formation of new vessels.

NB-OCTA holds a fundamental advantage over traditional GB-OCTA by enabling rapid imaging without scanning and stitching of subregions, which was previously impossible for samples with uneven surfaces. Real-time imaging is beneficial for operation navigation and monitoring during the surgical process⁴⁰. With the assistance of GPU⁴¹, NB-OCTA can provide high-resolution real-time cross-section images and rapid 3D volumetric data during surgery, assisting doctors in accurately locating lesion areas and organ structures, such as determining the range of tumor resection, which helps to maximize the preservation of normal tissue and enables doctors to make quick and precise decisions during surgery.

NB-OCTA shows promise across multiple applications, potentially enabling novel functions and providing improved imaging performance. For example, by combining NB-OCTA with Doppler technology⁴², it is possible to measure blood flow velocity in deep vessels. Furthermore, contrast agents⁴³ can alleviate the energy attenuation caused by the efficiency and elongation of the beam in NB-OCTA. NB-OCTA can also be used to guide and observe micro and nano robots⁴⁴, making it potentially suitable for targeted drug delivery and minimally invasive surgery.

In addition, technological improvements of NB-OCTA can extend its applications, expanding its relevancy across medical and biological research. NB-OCTA can benefit from adaptive optics technology to correct wavefront distortions caused by tissue heterogeneity, thereby obtaining higher-contrast and higher-resolution images⁴⁵. Combining NB-OCTA with other imaging modalities can provide complementary information⁴⁶, allowing analysis beyond the sum of the individual parts. For example, combining NB-OCT(A) with photoacoustic imaging can provide functional information about tissue oxygenation and blood flow as well as high axial resolution cross-section images⁴⁷, making diagnosis and monitoring of disease progression more accurate. Furthermore, artificial intelligence (AI) algorithms can both improve image quality and make analysis and diagnosis⁴⁸, reducing the workload of clinical doctors and improving the accuracy and reliability of diagnosis.

Our study provides a new approach for developing optical microscopy systems with consistently high spatial resolution and contributes to a better understanding of the mechanisms of ischemic injury and revascularization in stroke. Furthermore, we believe that NB-OCTA has great potential for a broader range of applications in the future by combining with other technologies.

DOE fabrication

The phase mask was created on a DOE which was fabricated on a 500-µm-thick fused silica wafer via four rounds of lithography. The relationship between the phase modulation P (x, y) and the height H (x, y) of the DOE can be expressed as

$$\begin{array}{c}\text{P}\text{ }\text{(}\text{x}\text{,}\text{ }\text{y}\text{)}\text{=}\frac{\text{2}\text{π}\text{(}{\text{n}}_{\text{DOE}}\text{-}\text{1)}\text{∙H}\text{ }\text{(}\text{x}\text{,}\text{ }\text{y}\text{)}}{\text{λ}} \# \left(2\right)\end{array}$$

where n_DOE is the refractive index of the material and λ is the wavelength. For fused silica, n_DOE =1.453 at 850nm, and the incremental thickness is 117nm (π / 8 in phase) for the 16 height levels. The first round of etching is to finish 8×117nm, the second for 4×117nm, the third for 2×117nm and the fourth for 117nm, combing to generate 16 heights. Each DOE has 1024×1024 pixels, with a 10-µm pixel resolution. The fabrication process was completed at the Stanford Nanofabrication Facility, and the detailed procedure can be found in our previous work^23,37. It takes approximately 10 hours to manufacture one batch of DOEs, considering all the time costs. The fabrication process achieves a height accuracy better than 95%, and the alignment error for the four rounds of photolithography is around 1 µm.

NB-OCTA system

A homebuilt SD-OCTA, as shown in Fig. 1c was used for imaging with a light source of broadband super luminescent diode (SLD-371HP3, Superlum, 850 ± 25nm, 15mW), giving an axial resolution of 6.4 µm. The light from the SLD passes through and a 2×2 coupler and is divided into the reference arm and sample arm. We placed the DOE in the sample arm before the objective lens (LSM02BB, Thorlabs, focal size = 5.8 µm), converting the conventional Gaussian beam-OCTA system to an NB-OCTA system. A dispersion compensation unit (LSM02DC, Thorlabs) was used in the sample arm to minimize the dispersion mismatch between the sample arm and the reference arm. The reflected reference light and sample light are combined through the coupler, and a home-built spectrometer with 0.05nm resolution, 2048 pixel, and 70kHz maximum A-scan rate is used to record the interference fringes.

Animal Preparation

All procedures involving mice were approved by the Institutional Animal Care and Use Committee of Tsinghua University. Wild-type mice (C57BL/6, 9 weeks to 3 months old; weight, 20-30g) were used in all the in vivo experiments.

We performed the craniotomy surgery as follows: First, mice were anesthetized using isoflurane (3% for induction, 1.5% during surgery), then a toe pinch procedure was conducted to verify the status of anesthesia. After confirming the animal's anesthesia status, the mouse was transferred onto a stereotactic frame with its body temperature maintained by using a heating pad. Before the surgery, ophthalmic ointment was topically applied to the corneal surfaces to prevent corneal drying and damage from prolonged exposure to the illuminating light source. An incision was made in the scalp and craniotomy was performed with a drill⁴⁹. After the skull was removed, the exposed brain was protected during the procedure from drying and infection. For acute experiments (Fig. 2g and 2h), the brain was covered with 2% low gelling temperature agarose (A4718, Sigma-Aldrich), and Phosphate-buffered saline (PBS) was used instead of water to dissolve the agarose to achieve the same osmotic pressure. A custom aluminum head post was then fixed to the skull using dental cement. For chronic imaging experiments (Fig. 3-Fig. 5), the brain was covered with 2% agarose and then a cover glass (0.17mm thickness, 6mm diameter) was glued to the skull to reduce motion of the exposed brain, and also a custom aluminum head post was fixed to the skull using dental cement. Imaging experiments were started ~ 1–2 weeks after chronic window implantation.

Photothrombotic vascular occlusion

After the mice recovered from the craniotomy surgery, PT occlusion was induced as follows: Rose Bengal (Sigma-Aldrich) solution (2.5mg/100g, 7.5mg/mL) was intraperitoneally injected into the mouse, and then the mouse was allowed to rest for 30 minutes to allow the solution to be absorbed into the systemic circulation. Then the occlusion was induced by focal illumination (1mm diameter focal spot, 30mW/mm²) with 532nm laser (MSL-III-532, CNI) in the somatosensory cortex region of the left cerebral hemisphere of the mouse brain for 30 minutes¹⁰. After PT occlusion was induced, we acquired OCTA data and then returned the mouse to its cage, where it had free access to food and water. After 24 hours, the mouse was anesthetized again for OCTA data acquisition. The experiment lasted for ten days, during which data collection was performed at a consistent time every day.

Data acquisition and processing

In this work, a stepwise raster scanning protocol was used to acquire volumetric dataset, the slow scanner (y) was driven by a step waveform with a total of 800 steps, at each step, five repeated B-scan were successively acquired for analyzing dynamic flow signals, with each B-scan composed of 1024 A-lines in the fast-scanning direction. Adjacent B-scans were paired for blood flow extraction at each scanning step y. The raw spectral interference fringe signal S (k, x, t) was Fourier transformed along the wavenumber direction to produce the backscattered profile C (z, x, t) in the depth space, which consists of amplitude and phase information:

$$\begin{array}{c}{\text{C}}_{\text{OCT}}\left(\text{x}\text{,}\text{ }\text{z}\text{,}\text{ }\text{t}\right)\text{ }\text{=}\text{ }\text{I}\left(\text{x}\text{,}\text{ }\text{z}\text{,}\text{ }\text{t}\right){\text{e}}^{\text{-}\text{iΦ}\left(\text{x}\text{,}\text{ }\text{z}\text{,}\text{ }\text{t}\right)}\#\#\left(3\right)\end{array}$$

The cross-sectional angiogram was then created by subtracting the amplitude component I (z, x, t) between paired B-scans:

$$\begin{array}{c}\text{O}\text{C}\text{T}\text{A}\text{ }\text{=}\text{ }\left|\text{I}\left(\text{z}\text{, }\text{x}\text{, }\text{i}\right)\text{ - }\text{I}\left(\text{z}\text{, }\text{x}\text{,}\text{ i+1}\right)\right| \# \left(4\right)\end{array}$$

Then the OCTA signals at each step y were averaged to improve flow contrast. Performing this process for all y values, OCTA volumetric data containing microvasculature flow information can be obtained. The en face angiograms are obtained by preforming maximum intensity projection (MIP) along the axial direction.

The depth-resolved image was color-encoded by assigning depth information to the hue channel, the 2D en face image to the value channel, and a constant value to the saturation channel. Subsequently, the color-encoded image was converted from HSV space to RGB space.

In order to assess the quality of the en face OCTA images acquired by GB and NB, several parameters were evaluated, including contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), and vessel connectivity. Prior to analysis, the en face OCTA images were normalized and converted to binary format by applying a threshold to differentiate blood vessels from static tissues, a skeleton map (M) was then generated. The CNR and SNR were calculated according to the following definitions:

$$\begin{array}{c}{\text{CNR}}_{\text{en face OCTA}}\text{ }\text{=}\text{ }\frac{\stackrel{-}{{\text{I}\left(\text{x}\text{,}\text{ }\text{y}\right)\text{|}}_{\text{M}\left(\text{x}\text{,}\text{ }\text{y}\right)\text{=1}}} \text{-}\text{ }\stackrel{-}{{\text{I}\left(\text{x}\text{,}\text{ }\text{y}\right)\text{|}}_{\text{Bck}}}}{{\text{σ}}_{\text{I}\left(\text{x,}\text{ }\text{y}\right)\text{|Bck}}}\#\left(5\right)\end{array}$$

$$\begin{array}{c}{\text{SNR}}_{\text{en face OCTA}}\text{=}\frac{\stackrel{-}{{\text{I}\left(\text{x,}\text{ }\text{y}\right)\text{|}}_{\text{M}\left(\text{x,}\text{ }\text{y}\right)\text{=}\text{1}}}}{{\text{σ}}_{\text{I}\left(\text{x,}\text{ }\text{y}\right)\text{|Bck}}}\#\left(6\right)\end{array}$$

where $\stackrel{-}{{\text{I}\left(\text{x,}\text{ }\text{y}\right)\text{|}}_{\text{M}\left(\text{x,}\text{ }\text{y}\right)\text{=}\text{1}}}$is the mean value of the intensities of the pixels that corresponded to the skeleton map in the en face OCTA image, $\stackrel{-}{{\text{I}\left(\text{x}\text{,}\text{ }\text{y}\right)\text{|}}_{\text{Bck}}}$represents the average intensity of pixels in the background (tissue) area of the en face image, and ${\text{σ}}_{\text{I}\left(\text{x,}\text{ }\text{y}\right)\text{|Bck}}$is the standard deviation of the background signals.

The connectivity is defined as the ratio of the biggest number of connected pixels and the total number of pixels that refers to vessels.

Acknowledgements

We thank Yang Lin and Xiaochen Li for helping with bio-sample preparation and discussing for data processing. This work was supported by the National Key R&D Program of China (No.2022YFB4702902), National Natural Science Foundation of China (Nos. 61831014, 62275023, and 32021002), Beijing Municipal Natural Science Foundation (No. 4232077), Overseas Expertise Introduction Project for Discipline Innovation (No. B18005), STI2030-Major Projects (No. 2022ZD0212000).

Author contributions

X.G. and J.Z. designed the experiments. X.G. and L.S. performed system integration and calibration. J.Z. designed and fabricated the DOEs. L.D. contributed to the mask preparation and wafer dividing. X.G. and J.Z conducted performance optimization. L.K. prepared the animals and supervised the biological experiments. X.G performed the imaging experiments. V.G. and K.S. contributed to stroke model processing and data analysis. X.G., J.Z., L.S. and Y.H. contributed to image processing. A.V.V. and K.L. contributed to editing of the manuscript. L.C. and Y.Y. contributed to the optimization of system design. G.J., A.d.l.Z. and Y.H. supervised the project. All authors provided critical revision feedback to the manuscript.

Data availability

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from corresponding authors.

Conflict of interest

The authors declare no competing interests.

Supplementary Information

Campbell, B. C. et al. Ischaemic stroke. Nature Reviews Disease Primers 5, 70 (2019).
Fluri, F., Schuhmann, M. K. & Kleinschnitz, C. Animal models of ischemic stroke and their application in clinical research. Drug design, development and therapy 9, 3445 (2015).
Sommer, C. J. Ischemic stroke: experimental models and reality. Acta neuropathologica 133, 245–261 (2017).
Robinson, R. G., Shoemaker, W. J., Schlumpf, M., Valk, T. & Bloom, F. E. Effect of experimental cerebral infarction in rat brain on catecholamines and behaviour. Nature 255, 332–334 (1975).
Llovera, G., Roth, S., Plesnila, N., Veltkamp, R. & Liesz, A. Modeling stroke in mice: permanent coagulation of the distal middle cerebral artery. Journal of visualized experiments: JoVE (2014).
Belayev, L., Alonso, O. F., Busto, R., Zhao, W. & Ginsberg, M. D. Middle cerebral artery occlusion in the rat by intraluminal suture: neurological and pathological evaluation of an improved model. Stroke 27, 1616–1623 (1996).
UlucÌʹ, K. Focal cerebral ischemia model by endovascular suture occlusion of the middle cerebral artery in the rat. (MyJoVE Corporation, 2016).
Kuroiwa, T. et al. Development of a rat model of photothrombotic ischemia and infarction within the caudoputamen. Stroke 40, 248–253 (2009).
Lu, H. et al. Induction and imaging of photothrombotic stroke in conscious and freely moving rats. Journal of Biomedical Optics 19, 096013–096013 (2014).
Yang, S. et al. Longitudinal in vivo intrinsic optical imaging of cortical blood perfusion and tissue damage in focal photothrombosis stroke model. Journal of Cerebral Blood Flow & Metabolism 39, 1381–1393 (2019).
Zhang, H. et al. Molecularly self-fueled nano-penetrator for nonpharmaceutical treatment of thrombosis and ischemic stroke. Nature Communications 14, 255 (2023).
He, L. et al. Highly bioactive zeolitic imidazolate framework-8–capped nanotherapeutics for efficient reversal of reperfusion-induced injury in ischemic stroke. Science Advances 6, eaay9751 (2020).
Murphy, T. H. & Corbett, D. Plasticity during stroke recovery: from synapse to behaviour. Nature reviews neuroscience 10, 861–872 (2009).
Molina, C. A. Reperfusion therapies for acute ischemic stroke: current pharmacological and mechanical approaches. Stroke 42, S16-S19 (2011).
Bhaskar, S., Stanwell, P., Cordato, D., Attia, J. & Levi, C. Reperfusion therapy in acute ischemic stroke: dawn of a new era? BMC neurology 18, 1–26 (2018).
Thomalla, G. et al. MRI-guided thrombolysis for stroke with unknown time of onset. New England Journal of Medicine 379, 611–622 (2018).
Yanev, P. & Dijkhuizen, R. M. In vivo imaging of neurovascular remodeling after stroke. Stroke 43, 3436–3441 (2012).
Carmeliet, P. Mechanisms of angiogenesis and arteriogenesis. Nature medicine 6, 389–395 (2000).
Risau, W. Mechanisms of angiogenesis. Nature 386, 671–674 (1997).
Ruan, L., Wang, B., ZhuGe, Q. & Jin, K. Coupling of neurogenesis and angiogenesis after ischemic stroke. Brain research 1623, 166–173 (2015).
Wang, X. et al. Functional MRI reveals brain-wide actions of thalamically-initiated oscillatory activities on associative memory consolidation. Nature Communications 14, 2195 (2023).
Tayyebi, S. et al. Diagnostic value of non-contrast brain computed tomography in the evaluation of acute cerebral venous thrombosis. Scientific reports 10, 883 (2020).
Cao, R. et al. Optical-resolution photoacoustic microscopy with a needle-shaped beam. Nature Photonics 17, 89–95 (2023).
Wang, T. & Xu, C. Three-photon neuronal imaging in deep mouse brain. Optica 7, 947–960 (2020).
Zhu, X. et al. Real-time whole-brain imaging of hemodynamics and oxygenation at micro-vessel resolution with ultrafast wide-field photoacoustic microscopy. Light: Science & Applications 11, 138 (2022).
Zong, W. et al. Fast high-resolution miniature two-photon microscopy for brain imaging in freely behaving mice. Nature methods 14, 713–719 (2017).
Burgess, A., Nhan, T., Moffatt, C., Klibanov, A. & Hynynen, K. Analysis of focused ultrasound-induced blood–brain barrier permeability in a mouse model of Alzheimer's disease using two-photon microscopy. Journal of controlled release 192, 243–248 (2014).
Li, D.-Y., Xia, Q., Yu, T.-T., Zhu, J.-T. & Zhu, D. Transmissive-detected laser speckle contrast imaging for blood flow monitoring in thick tissue: from Monte Carlo simulation to experimental demonstration. Light: Science & Applications 10, 241 (2021).
Drexler, W. & Fujimoto, J. G. Optical coherence tomography: technology and applications. Vol. 2 (Springer, 2015).
Spaide, R. F., Fujimoto, J. G., Waheed, N. K., Sadda, S. R. & Staurenghi, G. Optical coherence tomography angiography. Progress in retinal and eye research 64, 1–55 (2018).
Standish, B. A. et al. In vivo endoscopic multi-beam optical coherence tomography. Physics in Medicine & Biology 55, 615 (2010).
Lorenser, D., Singe, C. C., Curatolo, A. & Sampson, D. D. Energy-efficient low-Fresnel-number Bessel beams and their application in optical coherence tomography. Optics letters 39, 548–551 (2014).
Vairagi, K. et al. Common-path optical coherence tomography using the Bessel beam from negative axicon optical fiber tip. IEEE Journal of Selected Topics in Quantum Electronics 25, 1–6 (2018).
Zhang, M., Ren, Z. & Yu, P. Improve depth of field of optical coherence tomography using finite energy Airy beam. Optics letters 44, 3158–3161 (2019).
Jin, L. et al. Deep learning extended depth-of-field microscope for fast and slide-free histology. Proceedings of the National Academy of Sciences 117, 33051–33060 (2020).
Wu, Y. et al. Three-dimensional virtual refocusing of fluorescence microscopy images using deep learning. Nature methods 16, 1323–1331 (2019).
Zhao, J. et al. Flexible method for generating needle-shaped beams and its application in optical coherence tomography. Optica 9, 859–867 (2022).
He, B. et al. Whole Brain Micro-Vascular Imaging Using Robot Assisted Optical Coherence Tomography Angiography. IEEE Journal of Selected Topics in Quantum Electronics 29, 1–9 (2022).
Kovács, Z., Ikezaki, K., Samoto, K., Inamura, T. & Fukui, M. VEGF and flt: expression time kinetics in rat brain infarct. Stroke 27, 1865–1873 (1996).
Kut, C. et al. Detection of human brain cancer infiltration ex vivo and in vivo using quantitative optical coherence tomography. Science translational medicine 7, 292ra100-292ra100 (2015).
Wei, X. et al. Real-time cross-sectional and en face OCT angiography guiding high-quality scan acquisition. Optics letters 44, 1431–1434 (2019).
Li, Y., Chen, J. & Chen, Z. Advances in Doppler optical coherence tomography and angiography. Translational biophotonics 1, e201900005 (2019).
Nguyen, V. P. et al. Chain-like gold nanoparticle clusters for multimodal photoacoustic microscopy and optical coherence tomography enhanced molecular imaging. Nature Communications 12, 34 (2021).
Wu, Z. et al. A microrobotic system guided by photoacoustic computed tomography for targeted navigation in intestines in vivo. Science robotics 4, eaax0613 (2019).
Camino, A. et al. Depth-resolved optimization of a real-time sensorless adaptive optics optical coherence tomography. Optics letters 45, 2612–2615 (2020).
Nguyen, V. P. et al. Long-term, noninvasive in vivo tracking of progenitor cells using multimodality photoacoustic, optical coherence tomography, and fluorescence imaging. ACS nano 15, 13289–13306 (2021).
Hosseinaee, Z. et al. Functional and structural ophthalmic imaging using noncontact multimodal photoacoustic remote sensing microscopy and optical coherence tomography. Scientific Reports 11, 1–11 (2021).
Varadarajan, A. V. et al. Predicting optical coherence tomography-derived diabetic macular edema grades from fundus photographs using deep learning. Nature communications 11, 130 (2020).
Chen, T.-W. et al. Ultrasensitive fluorescent proteins for imaging neuronal activity. Nature 499, 295–300 (2013).

(Not answered)

Supplementaryinformation.docx

Download PDF

Version 1

posted

You are reading this latest preprint version

Needle-shaped beam optical coherence tomography angiography visualizing cortical blood perfusion after photothrombotic stroke in vivo

Status:

Version 1

Abstract

Figures

Introduction

Results

The principle of NB-OCTA

NB-OCTA system performance

Quantitative comparison between GB-OCTA and NB-OCTA

Long-term observation of blood perfusion after stroke with NB-OCTA

Discussion

Materials and methods

DOE fabrication

NB-OCTA system

Animal Preparation

Photothrombotic vascular occlusion

Data acquisition and processing

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1