The increasing use of mouse models for human brain disease studies presents an emerging need for a new functional imaging modality. Using optical excitation and acoustic detection, we developed a functional connectivity photoacoustic tomography system, which allows noninvasive imaging of resting-state functional connectivity in the mouse brain, with a large field of view and a high spatial resolution. Bilateral correlations were observed in eight functional regions, including the olfactory bulb, limbic, parietal, somatosensory, retrosplenial, visual, motor, and temporal regions, as well as in several subregions. The borders and locations of these regions agreed well with the Paxinos mouse brain atlas. By subjecting the mouse to alternating hyperoxic and hypoxic conditions, strong and weak functional connectivities were observed, respectively. In addition to connectivity images, vascular images were simultaneously acquired. These studies show that functional connectivity photoacoustic tomography is a promising, noninvasive technique for functional imaging of the mouse brain.fcPAT | RSFC | mouse brain functional imaging | hyperoxia | hypoxia

We have demonstrated the feasibility of imaging mouse brain metabolism using photoacoustic computed tomography (PACT), a fast, noninvasive and functional imaging modality with optical contrast and acoustic resolution. Brain responses to forepaw stimulations were imaged transdermally and transcranially. 2-NBDG, which diffuses well across the blood-brain-barrier, provided exogenous contrast for photoacoustic imaging of glucose response. Concurrently, hemoglobin provided endogenous contrast for photoacoustic imaging of hemodynamic response. Glucose and hemodynamic responses were quantitatively decoupled by using two-wavelength measurements. We found that glucose uptake and blood perfusion around the somatosensory region of the contralateral hemisphere were both increased by stimulations, indicating elevated neuron activity. While the glucose response area was more homogenous and confined within the somatosensory region, the hemodynamic response area had a clear vascular pattern and spread wider than the somatosensory region. Our results demonstrate that 2-NBDG-enhanced PACT is a promising tool for noninvasive studies of brain metabolism.

Photoacoustic imaging (PAI) is an emerging medical imaging modality capable of providing high spatial resolution of Ultrasound (US) imaging and high contrast of optical imaging. Delay-and-Sum (DAS) is the most common beamforming algorithm in PAI. However, using DAS beamformer leads to low resolution images and considerable contribution of off-axis signals. A new paradigm namely delay-multiply-and-sum (DMAS), which was originally used as a reconstruction algorithm in confocal microwave imaging, was introduced to overcome the challenges in DAS. DMAS was used in PAI systems and it was shown that this algorithm results in resolution improvement and sidelobe degrading. However, DMAS is still sensitive to high levels of noise, and resolution improvement is not satisfying. Here, we propose a novel algorithm based on DAS algebra inside DMAS formula expansion, double stage DMAS (DS-DMAS), which improves the image resolution and levels of sidelobe, and is much less sensitive to high level of noise compared to DMAS. The performance of DS-DMAS algorithm is evaluated numerically and experimentally. The resulted images are evaluated qualitatively and quantitatively using established quality metrics including signal-to-noise ratio (SNR), full-width-half-maximum (FWHM) and contrast ratio (CR). It is shown that DS-DMAS outperforms DAS and DMAS at the expense of higher computational load. DS-DMAS reduces the lateral valley for about 15 dB and improves the SNR and FWHM better than 13% and 30%, respectively. Moreover, the levels of sidelobe are reduced for about 10 dB in comparison with those in DMAS.

Abstract. In Photoacoustic imaging (PA), Delay-and-Sum (DAS) beamformer is a common beamforming algorithm having a simple implementation. However, it results in a poor resolution and high sidelobes. To address these challenges, a new algorithm namely Delay-Multiply-and-Sum (DMAS) was introduced having lower sidelobes compared to DAS. To improve the resolution of DMAS, a novel beamformer is introduced using Minimum Variance (MV) adaptive beamforming combined with DMAS, so-called Minimum Variance-Based DMAS (MVB-DMAS). It is shown that expanding the DMAS equation results in multiple terms representing a DAS algebra. It is proposed to use the MV adaptive beamformer instead of the existing DAS. MVB-DMAS is evaluated numerically and experimentally. In particular, at the depth of 45 mm MVB-DMAS results in about 31 dB, 18 dB and 8 dB sidelobes reduction compared to DAS, MV and DMAS, respectively. The quantitative results of the simulations show that MVB-DMAS leads to improvement in full-width-half-maximum about 96 %, 94 % and 45 % and signal-to-noise ratio about 89 %, 15 % and 35 % compared to DAS, DMAS, MV, respectively. In particular, at the depth of 33 mm of the experimental images, MVB-DMAS results in about 20 dB sidelobes reduction in comparison with other beamformers.

Optical coherence tomography (OCT) has the potential for skin tissue characterization due to its high axial and transverse resolution and its acceptable depth penetration. In practice, OCT cannot reach the theoretical resolutions due to imperfections of some of the components used. One way to improve the quality of the images is to estimate the point spread function (PSF) of the OCT system and deconvolve it from the output images. In this paper, we investigate the use of solid phantoms to estimate the PSF of the imaging system. We then utilize iterative Lucy-Richardson deconvolution algorithm to improve the quality of the images. The performance of the proposed algorithm is demonstrated on OCT images acquired from a variety of samples, such as epoxy-resin phantoms, fingertip skin and basaloid larynx and eyelid tissues.

With the growing application of photoacoustic imaging (PAI) in medical fields, there is a need to make them more compact, portable, and affordable. Therefore, we designed very low-cost PAI systems by replacing the expensive and sophisticated laser with a very low-energy laser diode. We implemented photoacoustic (PA) microscopy, both reflection and transmission modes, as well as PA computed tomography systems. The images obtained from tissue-mimicking phantoms and biological samples determine the feasibility of using a very low-energy laser diode in these configurations.

The current gold standard for clinical diagnosis of melanoma is excisional biopsy and histopathologic analysis. Approximately 15-30 benign lesions are biopsied to diagnose each melanoma. In addition, biopsies are invasive and result in pain, anxiety, scarring, and disfigurement of patients, which can add additional burden to the health care system. Among several imaging techniques developed to enhance melanoma diagnosis, optical coherence tomography (OCT), with its highresolution and intermediate penetration depth, can potentially provide required diagnostic information noninvasively. Here, we present an image analysis algorithm, "optical properties extraction (OPE)," which improves the specificity and sensitivity of OCT by identifying unique optical radiomic signatures pertinent to melanoma detection. We evaluated the performance of the algorithm using several tissue-mimicking phantoms and then tested the OPE algorithm on 69 human subjects. Our data show that benign nevi and melanoma can be differentiated with 97% sensitivity and 98% specificity. These findings suggest that the adoption of OPE algorithm in the clinic can lead to improvements in melanoma diagnosis and patient experience. Significance: This study describes a noninvasive, safe, simple-to-implement, and accurate method for the detection and differentiation of malignant melanoma versus benign nevi.

Spontaneous fluctuations of resting state functional MRI (rsfMRI) have been widely used to understand the macro-connectome of the human brain. However, these fluctuations are not synchronized among subjects, which leads to limitations and makes utilization of first-level model-based methods challenging. Considering this limitation of rsfMRI data in the time domain, we propose to transfer the spatiotemporal information of the rsfMRI data to another domain, the connectivity domain, in which each value represents the same effect across subjects. Using a set of seed networks and a connectivity index to calculate the functional connectivity for each seed network, we transform data into the connectivity domain by generating connectivity weights for each subject. Comparison of the two domains using a data-driven method suggests several advantages in analyzing data using data-driven methods in the connectivity domain over the time domain. We also demonstrate the feasibility of applying model-based methods in the connectivity domain, which offers a new pathway for the use of first-level model-based methods on rsfMRI data. The connectivity domain, furthermore, demonstrates a unique opportunity to perform first-level feature-based data-driven and model-based analyses. The connectivity domain can be constructed from any technique that identifies sets of features that are similar across subjects and can greatly help researchers in the study of macro-connectome brain function by enabling us to perform a wide range of model-based and data-driven approaches on rsfMRI data, decreasing susceptibility of analysis techniques to parameters that are not related to brain connectivity information, and evaluating both static and dynamic functional connectivity of the brain from a new perspective.