Search

Significance: Despite the ample progress made toward faster and more accurate Monte Carlo (MC) simulation tools over the past decade, the limited usability and accessibility of these advanced modeling tools remain key barriers to widespread use among the broad user community. Aim: An open-source, high-performance, web-based MC simulator that builds upon modern cloud computing architectures is highly desirable to deliver state-of-the-art MC simulations and hardware acceleration to general users without the need for special hardware installation and optimization. Approach: We have developed a configuration-free, in-browser 3D MC simulation platform—Monte Carlo eXtreme (MCX) Cloud—built upon an array of robust and modern technologies, including a Docker Swarm-based cloud-computing backend and a web-based graphical user interface (GUI) that supports in-browser 3D visualization, asynchronous data communication, and automatic data validation via JavaScript Object Notation (JSON) schemas. Results: The front-end of the MCX Cloud platform offers an intuitive simulation design, fast 3D data rendering, and convenient simulation sharing. The Docker Swarm container orchestration backend is highly scalable and can support high-demand GPU MC simulations using MCX over a dynamically expandable virtual cluster. Conclusion: MCX Cloud makes fast, scalable, and feature-rich MC simulations readily available to all biophotonics researchers without overhead. It is fully open-source and can be freely accessed at http://mcx.space/cloud.

Loop-mediated isothermal amplification or LAMP has been identified to be an efficient technology for point-of-care diagnostics. Paper-based LAMP technique has tremendous potential in replacing the existing tube-based technology as the manufacturing cost of a paper-based device is comparatively lower and easy-to-use. LAMP-based paper diagnostic device for Mycobacterium tuberculosis (MTB) detection is of extreme importance as it will help in early and rapid diagnosis of the affected patients. The fabrication of these devices requires assessment of design parameters on the extent of LAMP amplification reaction. Hence, CFD studies would be extremely beneficial from the design perspective. The current work presents an insight into the CFD simulations for LAMP amplification reaction on a porous paper membrane (nitrocellulose membrane). The convection-diffusion-reaction model is solved on a COMSOL Multiphysics 5.0 platform. Studies on effect of pore size, aspect ratio and initial DNA concentration on the extent of DNA amplification reaction have been carried out. The current paper-based technique is effective in detecting a minimum of 5 copies of DNA contrasting the previous semi-quantitative technique which demonstrated the detection of minimum 98 copies. Overall, the simulation results displayed almost 96% enhancement in the DNA amplification rate on paper membrane. Graphical abstract Graphical abstract for the computational study of DNA amplification reaction via LAMP technique on a porous paper membrane.

Significance Most shape-morphing materials rely on nonrenewable fossil resources or finely extracted biomaterials, which need strict reaction control, elaborate processing equipment, or prefabricated templates to achieve controllable transformation. To circumvent these challenges, we developed an eco-friendly and scalable strategy for programmable shape evolution that integrates easy-to-process pollen biomass with cost-effective digital printing. Using this approach, specific geometrical features and architectures were customized to build complex materials with user-defined, shape-morphing abilities. These fabrication efforts were complemented by computational simulations to build quantitative and mechanistic insights into the biomaterial’s characteristics for creating complex shapes with methods that are suitable for scalable manufacturing., We demonstrate how programmable shape evolution and deformation can be induced in plant-based natural materials through standard digital printing technologies. With nonallergenic pollen paper as the substrate material, we show how specific geometrical features and architectures can be custom designed through digital printing of patterns to modulate hygrophobicity, geometry, and complex shapes. These autonomously hygromorphing configurations can be “frozen” by postprocessing coatings to meet the needs of a wide spectrum of uses and applications. Through computational simulations involving the finite element method and accompanying experiments, we develop quantitative insights and a general framework for creating complex shapes in eco-friendly natural materials with potential sustainable applications for scalable manufacturing.

Significance: Simulations have the potential to be a powerful tool when planning the placement of photoacoustic imaging system components for surgical guidance. While elastic simulations (which include both compressional and shear waves) are expected to more accurately represent the physical transcranial acoustic wave propagation process, these simulations are more time-consuming and memory-intensive than the compressional-wave-only simulations that our group previously used to identify optimal acoustic windows for transcranial photoacoustic imaging. Aim: We present qualitative and quantitative comparisons of compressional and elastic wave simulations to determine which option is more suitable for preoperative surgical planning. Approach: Compressional and elastic photoacoustic k-Wave simulations were performed based on a computed tomography volume of a human cadaver head. Photoacoustic sources were placed in the locations of the internal carotid arteries and likely positions of neurosurgical instrument tips. Transducers received signals from three previously identified optimal acoustic windows (i.e., the ocular, nasal, and temporal regions). Target detectability, image-based target size estimates, and target-to-instrument distances were measured using the generalized contrast-to-noise ratio (gCNR), resolution, and relative source distances, respectively, for each simulation method. Results: The gCNR was equivalent between compressional and elastic simulations. The areas of the −6 dB contours of point spread functions utilized to measure resolution differed by 0.33 to 3.35 mm2. Target-to-instrument distance measurements were within 1.24 mm of the true distances. Conclusions: These results indicate that it is likely sufficient to utilize the less time-consuming, less memory-intensive compressional wave simulations for presurgical planning.

Significance: Diffuse light is ubiquitous in biomedical optics and imaging. Understanding the process of migration of an initial photon population entering tissue to a completely randomized, diffusely scattered population provides valuable insight to the interpretation and design of optical measurements. Aim: The goal of this perspective is to present a brief, unifying analytical framework to describe how properties of light transition from an initial state to a distributed state as light diffusion occurs. Approach: First, measurement parameters of light are introduced, and Monte Carlo simulations along with a simple analytical expression are used to explore how these individual parameters might exhibit diffusive behavior. Second, techniques to perform optical measurements are considered, highlighting how various measurement parameters can be leveraged to subsample photon populations. Results: Simulation results reinforce the fact that light undergoes a transition from a non-diffuse population to one that is first subdiffuse and then fully diffuse. Myriad experimental methods exist to isolate subpopulations of photons, which can be broadly categorized as source- and/or detector-encoded techniques, as well as methods of tagging the tissue of interest. Conclusions: Characteristic properties of light progressing to diffusion can be described by some form of Gaussian distribution that grows in space, time, angle, wavelength, polarization, and coherence. In some cases, these features can be approximated by simpler exponential behavior. Experimental methods to subsample features of the photon distribution can be achieved or theoretical methods can be used to better interpret the data with this framework.

Over the past decades, the application of single-photon emission computed tomography and positron emission tomography in neuroimaging has markedly increased. In the current study, we used a series of Korean computational head phantoms with detailed cranial structures for 6-, 9-, 12-, 15-y-old children and adult and a Monte Carlo transport code, MCNPX, to calculate age-dependent specific absorbed fraction (SAF) for mono-energetic electrons ranging from 0.01 to 4 MeV and S values for seven radionuclides widely used in nuclear medicine neuroimaging for the combination of ten source and target regions. Compared to the adult phantom, the 6-y phantom showed up to 1.7-fold greater SAF (cerebellum < cerebellum) and up to 1.4-fold greater S values (vitreous body < lens) for 123I. The electron SAF data, combined with our previous photon SAF data, will facilitate absorbed dose calculations for various cranial structures in patients undergoing neuroimaging procedures.

Significance: Coherence, a fundamental property of waves and fields, plays a key role in photoacoustic image reconstruction. Previously, techniques such as short-lag spatial coherence (SLSC) and filtered delay, multiply, and sum (FDMAS) have utilized spatial coherence to improve the reconstructed resolution and contrast with respect to delay-and-sum (DAS). While SLSC uses spatial coherence directly as the imaging contrast, FDMAS employs spatial coherence implicitly. Despite being more robust against noise, both techniques have their own drawbacks: SLSC does not preserve a relative signal magnitude, and FDMAS shows a reduced contrast-to-noise ratio. Aim: To overcome these limitations, our aim is to develop a beamforming algorithm—generalized spatial coherence (GSC)—that unifies SLSC and FDMAS into a single equation and outperforms both beamformers. Approach: We demonstrated the application of GSC in photoacoustic computed tomography (PACT) through simulation and experiments and compared it to previous beamformers: DAS, FDMAS, and SLSC. Results: GSC outperforms the imaging metrics of previous state-of-the-art coherence-based beamformers in both simulation and experiments. Conclusions: GSC is an innovative reconstruction algorithm for PACT, which combines the strengths of FDMAS and SLSC expanding PACT’s applications.

Significance: Single-molecule localization-based super-resolution microscopy has enabled the imaging of microscopic objects beyond the diffraction limit. However, this technique is limited by the requirements of imaging an extremely large number of frames of biological samples to generate a super-resolution image, thus requiring a longer acquisition time. Additionally, the processing of such a large image sequence leads to longer data processing time. Therefore, accelerating image acquisition and processing in single-molecule localization microscopy (SMLM) has been of perennial interest. Aim: To accelerate three-dimensional (3D) SMLM imaging by leveraging a computational approach without compromising the resolution. Approach: We used blind sparse inpainting to reconstruct high-density 3D images from low-density ones. The low-density images are generated using much fewer frames than usually needed, thus requiring a shorter acquisition and processing time. Therefore, our technique will accelerate 3D SMLM without changing the existing standard SMLM hardware system and labeling protocol. Results: The performance of the blind sparse inpainting was evaluated on both simulation and experimental datasets. Superior reconstruction results of 3D SMLM images using up to 10-fold fewer frames in simulation and up to 50-fold fewer frames in experimental data were achieved. Conclusions: We demonstrate the feasibility of fast 3D SMLM imaging leveraging a computational approach to reduce the number of acquired frames. We anticipate our technique will enable future real-time live-cell 3D imaging to investigate complex nanoscopic biological structures and their functions.

When relativistic electrons are used to irradiate tissues, such as during FLASH pre-clinical irradiations, the electron beam energy is one of the critical parameters that determine the dose distribution. Moreover, during such irradiations, linear accelerators (linacs) usually operate with significant beam loading, where a small change in the accelerator output current can lead to beam energy reduction. Optimisation of the tuning of the accelerator’s radio frequency system is often required. We describe here a robust, easy-to-use device for non-interceptive monitoring of potential variations in the electron beam energy during every linac macro-pulse of an irradiation run. Our approach monitors the accelerated electron fringe beam using two unbiased aluminium annular charge collection plates, positioned in the beam path and with apertures (5 cm in diameter) for the central beam. These plates are complemented by two thin annular screening plates to eliminate crosstalk and equalise the capacitances of the charge collection plates. The ratio of the charge picked up on the downstream collection plate to the sum of charges picked up on the both plates is sensitive to the beam energy and to changes in the energy spectrum shape. The energy sensitivity range is optimised to the investigated beam by the choice of thickness of the first plate. We present simulation and measurement data using electrons generated by a nominal 6 MeV energy linac as well as information on the design, the practical implementation and the use of this monitor.

Significance: Line scanning-based temporal focusing multiphoton microscopy (TFMPM) has superior axial excitation confinement (AEC) compared to conventional widefield TFMPM, but the frame rate is limited due to the limitation of the single line-to-line scanning mechanism. The development of the multiline scanning-based TFMPM requires only eight multiline patterns for full-field uniform multiphoton excitation and it still maintains superior AEC. Aim: The optimized parallel multiline scanning TFMPM is developed, and the performance is verified with theoretical simulation. The system provides a sharp AEC equivalent to the line scanning-based TFMPM, but fewer scans are required. Approach: A digital micromirror device is integrated in the TFMPM system and generates the multiline pattern for excitation. Based on the result of single-line pattern with sharp AEC, we can further model the multiline pattern to find the best structure that has the highest duty cycle together with the best AEC performance. Results: The AEC is experimentally improved to 1.7 μm from the 3.5 μm of conventional TFMPM. The adopted multiline pattern is akin to a pulse-width-modulation pattern with a spatial period of four times the diffraction-limited line width. In other words, ideally only four π/2 spatial phase-shift scans are required to form a full two-dimensional image with superior AEC instead of image-size-dependent line-to-line scanning. Conclusions: We have demonstrated the developed parallel multiline scanning-based TFMPM has the multiline pattern for sharp AEC and the least scans required for full-field uniform excitation. In the experimental results, the temporal focusing-based multiphoton images of disordered biotissue of mouse skin with improved axial resolution due to the near-theoretical limit AEC are shown to clearly reduce background scattering.

Significance: Raman spectroscopy (RS) applied to surgical guidance is attracting attention among scientists in biomedical optics. Offering a computational platform for studying depth-resolved RS and probing molecular specificity of different tissue layers is of crucial importance to increase the precision of these techniques and facilitate their clinical adoption. Aim: The aim of this work was to present a rigorous analysis of inelastic scattering depth sampling and elucidate the relationship between sensing depth of the Raman effect and optical properties of the tissue under interrogation. Approach: A new Monte Carlo (MC) package was developed to simulate absorption, fluorescence, elastic, and inelastic scattering of light in tissue. The validity of the MC algorithm was demonstrated by comparison with experimental Raman spectra in phantoms of known optical properties using nylon and polydimethylsiloxane as Raman-active compounds. A series of MC simulations were performed to study the effects of optical properties on Raman sensing depth for an imaging geometry consistent with single-point detection using a handheld fiber optics probe system. Results: The MC code was used to estimate the Raman sensing depth of a handheld fiber optics system. For absorption and reduced scattering coefficients of 0.001 and 1 mm − 1, the sensing depth varied from 105 to 225 μm for a range of Raman probabilities from 10 − 6 to 10 − 3. Further, for a realistic Raman probability of 10 − 6, the sensing depth ranged between 10 and 600 μm for the range of absorption coefficients 0.001 to 1.4 mm − 1 and reduced scattering coefficients of 0.5 to 30 mm − 1. Conclusions: A spectroscopic MC light transport simulation platform was developed and validated against experimental measurements in tissue phantoms and used to predict depth sensing in tissue. It is hoped that the current package and reported results provide the research community with an effective simulating tool to improve the development of clinical applications of RS.

Significance: Photoacoustic imaging (PAI) has been greatly developed in a broad range of diagnostic applications. The efficiency of light to sound conversion in PAI is limited by the ubiquitous noise arising from the tissue background, leading to a low signal-to-noise ratio (SNR), and thus a poor quality of images. Frame averaging has been widely used to reduce the noise; however, it compromises the temporal resolution of PAI. Aim: We propose an approach for photoacoustic (PA) signal denoising based on a combination of low-pass filtering and sparse coding (LPFSC). Approach: LPFSC method is based on the fact that PA signal can be modeled as the sum of low frequency and sparse components, which allows for the reduction of noise levels using a hybrid alternating direction method of multipliers in an optimization process. Results: LPFSC method was evaluated using in-silico and experimental phantoms. The results show a 26% improvement in the peak SNR of PA signal compared to the averaging method for in-silico data. On average, LPFSC method offers a 63% improvement in the image contrast-to-noise ratio and a 33% improvement in the structural similarity index compared to the averaging method for objects located at three different depths, ranging from 10 to 20 mm, in a porcine tissue phantom. Conclusions: The proposed method is an effective tool for PA signal denoising, whereas it ultimately improves the quality of reconstructed images, especially at higher depths, without limiting the image acquisition speed.

The high-Tc superconducting (HTS) dynamo is a promising device that can inject large DC supercurrents into a closed superconducting circuit. This is particularly attractive to energise HTS coils in NMR/MRI magnets and superconducting rotating machines without the need for connection to a power supply via current leads. It is only very recently that quantitatively accurate, predictive models have been developed which are capable of analysing HTS dynamos and explain their underlying physical mechanism. In this work, we propose to use the HTS dynamo as a new benchmark problem for the HTS modelling community. The benchmark geometry consists of a permanent magnet rotating past a stationary HTS coated-conductor wire in the open-circuit configuration, assuming for simplicity the 2D (infinitely long) case. Despite this geometric simplicity the solution is complex, comprising time-varying spatially-inhomogeneous currents and fields throughout the superconducting volume. In this work, this benchmark problem has been implemented using several different methods, including H-formulation-based methods, coupled H-A and T-A formulations, the Minimum Electromagnetic Entropy Production method, and integral equation and volume integral equation-based equivalent circuit methods. Each of these approaches show excellent qualitative and quantitative agreement for the open-circuit equivalent instantaneous voltage and the cumulative time-averaged equivalent voltage, as well as the current density and electric field distributions within the HTS wire at key positions during the magnet transit. A critical analysis and comparison of each of the modelling frameworks is presented, based on the following key metrics: number of mesh elements in the HTS wire, total number of mesh elements in the model, number of degrees of freedom (DOFs), tolerance settings and the approximate time taken per cycle for each model., 22 pages, 7 figures, 2 tables

Significance: A scene-based adaptive-optics (AO) system is developed and a method for investigating its imaging performance is proposed. The system enables derivation of Strehl ratios from observed images via collaboration with computer simulations. The resultant Strehl ratios are comparable with those of other current AO systems. Aim: For versatile and noninvasive AO microscopy, a scene-based wavefront-sensing technique working on a Shack–Hartmann wavefront sensor is developed in a modal control system. The purpose of the research is to clarify the imaging performance of the AO system via the derivation of Strehl ratios from observed images toward applications in microscopy of living cells and tissues. Approach: Two imaging metrics that can be directly measured from observed images (i.e., an energy concentration ratio and unbiased maximum ratio) are defined and related to the Strehl ratio via computer simulations. Experiments are conducted using artificial targets to measure the imaging metrics, which are then converted to Strehl ratios. Results: The resultant Strehl ratios are >0.7 and 0.5 in the cases of defocus and higher aberrations, respectively. The half-widths at half-maximum of the AO-corrected bead images are favorably comparable to those of on-focus images under simple defocus aberration, and the AO system works both under bright-field illumination and on fluorescent bead images. Conclusions: The proposed scene-based AO system is expected to work with a Strehl ratio of more than 0.5 when applied to high-resolution live imaging of cells and tissues under bright-field and fluorescence microscopies.

The hybrid trapped field magnet lens (HTFML) is a promising device that is able to concentrate a magnetic field higher than an applied background field continuously, even after removing a background field, which was conceptually proposed by the authors in 2018. We have numerically investigated the HTFML performance, consisting of a REBaCuO cylindrical magnetic lens and REBaCuO trapped field magnet (TFM) cylinder, magnetized by pulsed fields. Single magnetic pulses were applied ranging from B app = 1.5 T to 5.0 T at the operating temperature of T s = 30, 40 and 50 K, and the performance was compared with that of the single REBaCuO TFM cylinder. The HTFML effect was clearly confirmed for the lower B app values. However, for the higher B app values, the trapped field in the magnetic lens bore was nearly equal to or slightly lower than that for the single TFM cylinder because of a weakened lens effect due to magnetic flux penetration into the lens. A temperature rise in the REBaCuO magnetic lens and TFM cylinder was also observed. These results strongly suggest that lowering the temperature of the REBaCuO magnetic lens could enhance the HTFML effect even for higher B app., Adaptable and Seamless Technology Transfer Program through Target-driven R&D (A-STEP), Japan Science and Technology Agency (JST), Grant No. VP30218088419 JSPS KAKENHI Grant No.19K05240

Significance: It is commonly assumed that using the objective lens to create a tightly focused light spot for illumination provides a twofold resolution improvement over the Rayleigh resolution limit and that resolution improvement is independent of object properties. Nevertheless, such an assumption has not been carefully examined. We examine this assumption by analyzing the performance of two super-resolution methods, known as image scanning microscopy (ISM) and illumination-enhanced sparsity (IES). Aim: We aim to identify the fundamental differences between the two methods, and to provide examples that help researchers determine which method to utilize for different imaging conditions. Approach: We input the same image datasets into the two methods and analyze their restorations. In numerical simulations, we design objects of distinct brightness and sparsity levels for imaging. We use biological imaging experiments to verify the simulation results. Results: The resolution of IES often exceeds twice the Rayleigh resolution limit when imaging sparse objects. A decrease in object sparsity negatively affects the resolution improvement in both methods. Conclusions: The IES method is superior for imaging sparse objects with its main features being bright and small against a dark, large background. For objects that are largely bright with small dark features, the ISM method is favorable.

To detect small-scale changes in tissue with optical techniques, small sampling volumes are required. Single fiber reflectance (SFR) spectroscopy has a sampling depth of a few hundred micrometers. SFR spectroscopy uses a single fiber to emit and collect light. The only available model to determine optical properties with SFR spectroscopy was derived for tissues with modified Henyey-Greenstein phase functions. Previously, we demonstrated that this model is inadequate for other tissue phase functions. We develop a model to relate SFR measurements to scattering properties for a range of phase functions, in the absence of absorption. Since the source and detector overlap, the reflectance cannot be accurately described by diffusion theory alone: SFR measurements are subdiffuse. Therefore, we describe the reflectance as a combination of a diffuse and a semiballistic component. We use the model of Farrell et al. for the diffuse component, solved for an overlapping source and detector fiber. For the semiballistic component, we derive a new parameter, psb, which incorporates the integrals of the phase function over 1 deg in the backward direction and 23 deg in the forward direction. Our model predicts the reflectance with a median error of 2.1%, compared to 9.0% for the currently available model..

Since our previous publication of organ dose for the pediatric CT cohort in the UK, there have been questions about the magnitude of uncertainty in our dose estimates. We therefore quantified shared and unshared uncertainties in empirical CT parameters extracted from 1073 CT films (1978–2008) from 36 hospitals in the study and propagated these uncertainties into organ doses using Monte Carlo random sampling and NCICT organ dose calculator. The average of 500 median brain and marrow doses for the full cohort was 35 (95% confidence interval: 30–40) mGy and 6 (5–7) mGy, respectively. We estimated that shared uncertainty contributed ~99% of coefficient of variation of median brain doses in brain scans compared to unshared uncertainty (1% contribution). We found that the previous brain doses were slightly underestimated for 1990 compared to the results in the current study due to the revised CTDI models based on CT films.

Significance: Ultrasound-assisted optical imaging techniques, such as ultrasound-modulated optical tomography, allow for imaging deep inside scattering media. In these modalities, a fraction of the photons passing through the ultrasound beam is modulated. The efficiency by which the photons are converted is typically referred to as the ultrasound modulation’s “tagging efficiency.” Interestingly, this efficiency has been defined in varied and discrepant fashion throughout the scientific literature. Aim: The aim of this study is the ultrasound tagging efficiency in a manner consistent with its definition and experimentally verify the contributive (or noncontributive) relationship between the mechanisms involved in the ultrasound optical modulation process. Approach: We adopt a general description of the tagging efficiency as the fraction of photons traversing an ultrasound beam that is frequency shifted (inclusion of all frequency-shifted components). We then systematically studied the impact of ultrasound pressure and frequency on the tagging efficiency through a balanced detection measurement system that measured the power of each order of the ultrasound tagged light, as well as the power of the unmodulated light component. Results: Through our experiments, we showed that the tagging efficiency can reach 70% in a scattering phantom with a scattering anisotropy of 0.9 and a scattering coefficient of 4 mm⁻¹ for a 1-MHz ultrasound with a relatively low (and biomedically acceptable) peak pressure of 0.47 MPa. Furthermore, we experimentally confirmed that the two ultrasound-induced light modulation mechanisms, particle displacement and refractive index change, act in opposition to each other. Conclusion: Tagging efficiency was quantified via simulation and experiments. These findings reveal avenues of investigation that may help improve ultrasound-assisted optical imaging techniques.

The mesh-based Monte Carlo (MMC) algorithm is increasingly used as the gold-standard for developing new biophotonics modeling techniques in 3-D complex tissues, including both diffusion-based and various Monte Carlo (MC)-based methods. Compared to multilayered and voxel-based MCs, MMC can utilize tetrahedral meshes to gain improved anatomical accuracy but also results in higher computational and memory demands. Previous attempts of accelerating MMC using graphics processing units (GPUs) have yielded limited performance improvement and are not publicly available. We report a highly efficient MMC—MMCL—using the OpenCL heterogeneous computing framework and demonstrate a speedup ratio up to 420× compared to state-of-the-art single-threaded CPU simulations. The MMCL simulator supports almost all advanced features found in our widely disseminated MMC software, such as support for a dozen of complex source forms, wide-field detectors, boundary reflection, photon replay, and storing a rich set of detected photon information. Furthermore, this tool supports a wide range of GPUs/CPUs across vendors and is freely available with full source codes and benchmark suites at http://mcx.space/#mmc.

As photoacoustic imaging (PAI) technology matures, computational modeling will increasingly represent a critical tool for facilitating clinical translation through predictive simulation of real-world performance under a wide range of device and biological conditions. While modeling currently offers a rapid, inexpensive tool for device development and prediction of fundamental image quality metrics (e.g., spatial resolution and contrast ratio), rigorous verification and validation will be required of models used to provide regulatory-grade data that effectively complements and/or replaces in vivo testing. To address methods for establishing model credibility, we developed an integrated computational model of PAI by coupling a previously developed three-dimensional Monte Carlo model of tissue light transport with a two-dimensional (2D) acoustic wave propagation model implemented in the well-known k-Wave toolbox. We then evaluated ability of the model to predict basic image quality metrics by applying standardized verification and validation principles for computational models. The model was verified against published simulation data and validated against phantom experiments using a custom PAI system. Furthermore, we used the model to conduct a parametric study of optical and acoustic design parameters. Results suggest that computationally economical 2D acoustic models can adequately predict spatial resolution, but metrics such as signal-to-noise ratio and penetration depth were difficult to replicate due to challenges in modeling strong clutter observed in experimental images. Parametric studies provided quantitative insight into complex relationships between transducer characteristics and image quality as well as optimal selection of optical beam geometry to ensure adequate image uniformity. Multidomain PAI simulation tools provide high-quality tools to aid device development and prediction of real-world performance, but further work is needed to improve model fidelity, especially in reproducing image noise and clutter.

In this work, we characterize the impact of large-volume processing upon the analytical sensitivity of flow-through paper-based immunoassays. Larger sample volumes feature greater molar quantities of available analyte, but the assay design principles which would enable the rapid collection of this dilute target are ill-defined. We developed a finite-element model to explore the operating conditions under which processing large sample volumes via pressure-driven convective flow would yield an improved binding signal. Our simulation results underscore the importance of establishing a high local concentration of the analyte-binding species within the porous substrate. This elevated abundance serves to enhance the binding kinetics, matching the time scale of target capture to the period during which the sample is in contact with the test zone (i.e., the effective residence time). These findings were experimentally validated using the rcSso7d-cellulose-binding domain (CBD) fusion construct, a bifunctional binding protein which adsorbs to cellulose in high abundance. As predicted by our modeling efforts, the local concentration achieved using the rcSso7d-CBD species is uniquely enabling for sensitivity enhancement through large-volume processing. The rapid analyte depletion which occurs at this high surface density also permits the processing of large sample volumes within practical time scales and flow regimes. Using these findings, we present guidance for the optimal means of processing large sample volumes for enhanced assay sensitivity.

Machine learning methods and in particular random forests are promising approaches for prediction based on high dimensional omics data sets. They provide variable importance measures to rank predictors according to their predictive power. If building a prediction model is the main goal of a study, often a minimal set of variables with good prediction performance is selected. However, if the objective is the identification of involved variables to find active networks and pathways, approaches that aim to select all relevant variables should be preferred. We evaluated several variable selection procedures based on simulated data as well as publicly available experimental methylation and gene expression data. Our comparison included the Boruta algorithm, the Vita method, recurrent relative variable importance, a permutation approach and its parametric variant (Altmann) as well as recursive feature elimination (RFE). In our simulation studies, Boruta was the most powerful approach, followed closely by the Vita method. Both approaches demonstrated similar stability in variable selection, while Vita was the most robust approach under a pure null model without any predictor variables related to the outcome. In the analysis of the different experimental data sets, Vita demonstrated slightly better stability in variable selection and was less computationally intensive than Boruta. In conclusion, we recommend the Boruta and Vita approaches for the analysis of high-dimensional data sets. Vita is considerably faster than Boruta and thus more suitable for large data sets, but only Boruta can also be applied in low-dimensional settings.

RNA sequencing (RNA-seq) has become a standard procedure to investigate transcriptional changes between conditions and is routinely used in research and clinics. While standard differential expression (DE) analysis between two conditions has been extensively studied, and improved over the past decades, RNA-seq time course (TC) DE analysis algorithms are still in their early stages. In this study, we compare, for the first time, existing TC RNA-seq tools on an extensive simulation data set and validated the best performing tools on published data. Surprisingly, TC tools were outperformed by the classical pairwise comparison approach on short time series (, Briefings in Bioinformatics, 20 (1), ISSN:1467-5463, ISSN:1477-4054

Rank aggregation (RA), the process of combining multiple ranked lists into a single ranking, has played an important role in integrating information from individual genomic studies that address the same biological question. In previous research, attention has been focused on aggregating full lists. However, partial and/or top ranked lists are prevalent because of the great heterogeneity of genomic studies and limited resources for follow-up investigation. To be able to handle such lists, some ad hoc adjustments have been suggested in the past, but how RA methods perform on them (after the adjustments) has never been fully evaluated. In this article, a systematic framework is proposed to define different situations that may occur based on the nature of individually ranked lists. A comprehensive simulation study is conducted to examine the performance characteristics of a collection of existing RA methods that are suitable for genomic applications under various settings simulated to mimic practical situations. A non-small cell lung cancer data example is provided for further comparison. Based on our numerical results, general guidelines about which methods perform the best/worst, and under what conditions, are provided. Also, we discuss key factors that substantially affect the performance of the different methods.

Despite decades of epidemiological research, it remains uncertain whether ionizing radiation can cause lymphomas. Most epidemiological studies of lymphoma risk following non-uniform exposure used dose to red bone marrow (RBM), constituting a small fraction of the lymphocytes, as a surrogate of dose to the lymphocytes. We developed a method to estimate dose to the lymphocytes using the reference distribution of lymphocytes throughout the body and Monte Carlo simulations of computational human phantoms. We applied our method to estimating lymphocyte doses for a pediatric CT patient cohort in the United Kingdom. Estimated dose to the RBM was greater than lymphocyte dose for most scan types (up to 2.6-fold higher, a 5-year-old brain scan) except abdomen scan (RBM dose was about half the lymphocyte dose, a 5-year-old abdomen scan). The lymphocyte dose in the UK cohort showed that T-spine and whole body scans delivered the highest lymphocyte doses (up to 22.4 mGy).

DNA methylation is one of the most extensively studied epigenetic modifications of genomic DNA. In recent years, sequencing of bisulfite-converted DNA, particularly via next-generation sequencing technologies, has become a widely popular method to study DNA methylation. This method can be readily applied to a variety of species, dramatically expanding the scope of DNA methylation studies beyond the traditionally studied human and mouse systems. In parallel to the increasing wealth of genomic methylation profiles, many statistical tools have been developed to detect differentially methylated loci (DMLs) or differentially methylated regions (DMRs) between biological conditions. We discuss and summarize several key properties of currently available tools to detect DMLs and DMRs from sequencing of bisulfite-converted DNA. However, the majority of the statistical tools developed for DML/DMR analyses have been validated using only mammalian data sets, and less priority has been placed on the analyses of invertebrate or plant DNA methylation data. We demonstrate that genomic methylation profiles of non-mammalian species are often highly distinct from those of mammalian species using examples of honey bees and humans. We then discuss how such differences in data properties may affect statistical analyses. Based on these differences, we provide three specific recommendations to improve the power and accuracy of DML and DMR analyses of invertebrate data when using currently available statistical tools. These considerations should facilitate systematic and robust analyses of DNA methylation from diverse species, thus advancing our understanding of DNA methylation.

In recent decades, applications of single photon emission computed tomography and positron emission tomography in clinical neuroimaging have markedly increased. In this study, we developed a series of Korean computational head phantoms with detailed cranial substructures for 6-, 9-, 12- and 15-year-old children and adult by non-uniformly adjusting a template head phantom to match the Korean standard head dimensions. The Korean head phantoms were coupled with a Monte Carlo transport code to calculate age-dependent specific absorbed fraction (SAF) for the combination of 10 source and target regions and mono-energetic photons ranging from 0.01 to 4 MeV. Compared to the adult phantom, the 6-y phantom showed up to 1.4-fold greater self-absorption SAF (cerebellum) and up to 1.8-fold greater cross-irradiation SAF (cerebellum < eye balls). With addition of electron SAFs in the future, our photon SAF data will facilitate dose calculations for various cranial substructures in patients undergoing cranial neuroimaging procedures.

Since it was first demonstrated more than a decade ago, the single-pixel camera concept has been used in numerous applications in which it is necessary or advantageous to reduce the channel count, cost, or data volume. Here, three-dimensional (3-D), compressed-sensing photoacoustic tomography (PAT) is demonstrated experimentally using a single-pixel camera. A large area collimated laser beam is reflected from a planar Fabry–Pérot ultrasound sensor onto a digital micromirror device, which patterns the light using a scrambled Hadamard basis before it is collected into a single photodetector. In this way, inner products of the Hadamard patterns and the distribution of thickness changes of the FP sensor—induced by the photoacoustic waves—are recorded. The initial distribution of acoustic pressure giving rise to those photoacoustic waves is recovered directly from the measured signals using an accelerated proximal gradient-type algorithm to solve a model-based minimization with total variation regularization. Using this approach, it is shown that 3-D PAT of imaging phantoms can be obtained with compression rates as low as 10%. Compressed sensing approaches to photoacoustic imaging, such as this, have the potential to reduce the data acquisition time as well as the volume of data it is necessary to acquire, both of which are becoming increasingly important in the drive for faster imaging systems giving higher resolution images with larger fields of view.

Interstitial photodynamic therapy (iPDT) with 5-aminolevulinic acid (ALA) is a possible alternative treatment for malignant brain tumors. Further evaluation is, however, required before it can be clinically applied. Computational simulation of the photophysical process in ALA-iPDT can offer a quantitative tool for understanding treatment outcomes, which depend on various variables related to clinical treatment conditions. We propose a clinical simulation method of ALA-iPDT for malignant brain tumors using a singlet oxygen (O12) model and O12 threshold to induce cell death. In this method, the amount of O12 generated is calculated using a photosensitizer photobleaching coefficient and O12 quantum yield, which have been measured in several previous studies. Results of the simulation using clinical magnetic resonance imaging data show the need to specify the insertion positions of cylindrical light diffusers and the level of light fluence. Detailed analysis with a numerical brain tumor model demonstrates that ALA-iPDT treatment outcomes depend on combinations of photobleaching and threshold values. These results indicate that individual medical procedures, including pretreatment planning and treatment monitoring, will greatly benefit from simulation of ALA-iPDT outcomes.

Multiphoton microscopy provides a suitable technique for imaging biological tissues with submicrometer resolution. Usually a Gaussian beam (GB) is used for illumination, leading to a reduced power efficiency in the multiphoton response and vignetting for a square-shaped imaging area. A flat-top beam (FTB) provides a uniform spatial intensity distribution that equalizes the probability of a multiphoton effect across the imaging area. We employ a customized widefield multiphoton microscope to compare the performance of a square-shaped FTB illumination with that based on using a GB, for both two-photon fluorescence (TPF) and second-harmonic generation (SHG) imaging. The variation in signal-to-noise ratio across TPF images of fluorescent dyes spans ∼5.6 dB for the GB and ∼1.2 dB for the FTB illumination, respectively. For the GB modality, TPF images of mouse colon and Convallaria root, and SHG images of chicken tendon and human breast biopsy tissue showcase ∼20% area that are not imaged due to either insufficient or lack of illumination. For quantitative analysis that depends on the illuminated area, this effect can potentially lead to inaccuracies. This work emphasizes the applicability of FTB illumination to multiphoton applications.

Near-infrared spectroscopy (NIRS) is widely used in biomedical optics with applications ranging from basic science, such as in functional neuroimaging, to clinical, as in pulse oximetry. Despite the relatively low absorption of tissue in the near-infrared, there is still a significant amount of optical attenuation produced by the highly scattering nature of tissue. Because of this, designers of NIRS systems have to balance source optical power and source–detector separation to maximize the signal-to-noise ratio (SNR). However, theoretical estimations of SNR neglect the effects of speckle. Speckle manifests as fluctuations of the optical power received at the detector. These fluctuations are caused by interference of the multiple random paths taken by photons in tissue. We present a model for the NIRS SNR that includes the effects of speckle. We performed experimental validations with a NIRS system to show that it agrees with our model. Additionally, we performed computer simulations based on the model to estimate the contribution of speckle noise for different collection areas and source–detector separations. We show that at short source–detector separation, speckle contributes most of the noise when using long coherence length sources. Considering this additional noise is especially important for hybrid applications that use NIRS and speckle contrast simultaneously, such as in diffuse correlation spectroscopy.

X-ray luminescence computed tomography (XLCT) is an emerging hybrid imaging modality which has the potential for achieving both high sensitivity and spatial resolution simultaneously. For the narrow x-ray beam-based XLCT imaging, based on previous work, a spatial resolution of about double the x-ray beam size can be achieved using a translate/rotate scanning scheme, taking step sizes equal to the x-ray beam width. To break the current spatial resolution limit, we propose a scanning strategy achieved by reducing the scanning step size to be smaller than the x-ray beam size. We performed four sets of numerical simulations and a phantom experiment using cylindrical phantoms and have demonstrated that our proposed scanning method can greatly improve the XLCT-reconstructed image quality compared with the traditional scanning approach. In our simulations, by using a fixed x-ray beam size of 0.8 mm, we were able to successfully reconstruct six embedded targets as small as 0.5 mm in diameter and with the same edge-to-edge distances by using a scanning step as small as 0.2 mm which is a 1.6 times improvement in the spatial resolution compared with the traditional approach. Lastly, the phantom experiment further demonstrated the efficacy of our proposed method in improving the XLCT image quality, with all image quality metrics improving as the step size decreased.