Physics in medicine and biology最新文献

Objective. The acceleration of magnetic resonance imaging (MRI) acquisition is crucial for both clinical and research applications, particularly in dynamic MRI. Existing compressed sensing (CS) methods, despite being effective for fast imaging, face limitations such as the need for incoherent sampling and residual noise, which restrict their practical use for rapid MRI.Approach. To overcome these challenges, we propose a novel image reconstruction framework that integrates the MRI physical model with a flexible, self-adjusting, decoupling data-driven model. We validated this method through experiments using both simulated andin vivodynamic contrast-enhanced MRI datasets.Main results. The experimental results demonstrate that the proposed framework achieves high spatial and temporal resolution reconstructions. Additionally, when compared to state-of-the-art image reconstruction approaches, our method significantly enhances acceleration capabilities, enabling sparse and rapid imaging with high resolution.Significance. Our proposed framework offers a promising solution for real-time imaging and image-guided radiation therapy applications by providing superior performance in achieving high spatial and temporal resolution reconstructions, thus addressing the limitations of existing CS schemes.

Objective. MRI is the standard imaging modality for high-dose-rate brachytherapy of cervical cancer. Precise contouring of organs at risk (OARs) and high-risk clinical target volume (HR-CTV) from MRI is a crucial step for radiotherapy planning and treatment. However, conventional manual contouring has limitations in terms of accuracy as well as procedural time. To overcome these, we propose a deep learning approach to automatically segment OARs (bladder, rectum, and sigmoid colon) and HR-CTV from female pelvic MRI.Approach. In the proposed pipeline, a coarse multi-organ segmentation model first segments all structures, from which a region of interest is computed for each structure. Then, each organ is segmented using an organ-specific fine segmentation model separately trained for each organ. To account for variable sizes of HR-CTV, a size-adaptive multi-model approach was employed. For coarse and fine segmentations, we designed a dual convolution-transformer UNet (DCT-UNet) which uses dual-path encoder consisting of convolution and transformer blocks. To evaluate our model, OAR segmentations were compared to the clinical contours drawn by the attending radiation oncologist. For HR-CTV, four sets of contours (clinical + three additional sets) were obtained to produce a consensus ground truth as well as for inter/intra-observer variability analysis.Main results. DCT-UNet achieved dice similarity coefficient (mean ± SD) of 0.932 ± 0.032 (bladder), 0.786 ± 0.090 (rectum), 0.663 ± 0.180 (sigmoid colon), and 0.741 ± 0.076 (HR-CTV), outperforming other state-of-the-art models. Notably, the size-adaptive multi-model significantly improved HR-CTV segmentation compared to a single-model. Furthermore, significant inter/intra-observer variability was observed, and our model showed comparable performance to all observers. Computation time for the entire pipeline per subject was 12.59 ± 0.79 s, which is significantly shorter than the typical manual contouring time of >15 min.Significance. These experimental results demonstrate that our model has great utility in cervical cancer brachytherapy by enabling fast and accurate automatic segmentation, and has potential in improving consistency in contouring. DCT-UNet source code is available athttps://github.com/JHU-MICA/DCT-UNet.

Introduction: Ion therapies have an increased relative biological effectiveness (RBE) compared to X-rays, but this remains poorly quantified across different radiation qualities. Mechanistic models that simulate DNA damage and repair after irradiation could be used to help better quantify RBE. However, there is large variation in model design with the simulation detail and number of parameters required to accurately predict key biological endpoints remaining unclear. This work investigated damage models with varying detail to determine how different model features impact the predicted DNA damage. Methods: Damage models of reducing detail were designed in TOPAS-nBio and Medras investigating the inclusion of chemistry, realistic nuclear geometries, single strand break damage, and track structure. The nucleus models were irradiated with 1 Gy of protons across a range of linear energy transfers (LETs). Damage parameters in the models with reduced levels of simulation detail were fit to proton double strand break (DSB) yield predicted by the most detailed model. Irradiation of the optimised models with a range of radiation qualities was then simulated, before undergoing repair in the Medras biological response model. Results: Simplified damage models optimised to proton exposures predicted similar trends in DNA damage across radiation qualities. On average across radiation qualities, the simplified models experienced an 8% variation in double strand break (DSB) yield but a larger 28% variation in chromosome aberrations. Aberration differences became more prominent at higher LETs, with model features having an increasing impact on the distribution and therefore misrepair of DSBs. However, overall trends remained similar with better agreement likely achievable through repair model optimisation. Conclusion: Several model simplifications could be made without compromising key damage yield predictions, although changes in damage complexity and distribution were observed. This suggests simpler, more efficient models may be sufficient for initial radiation damage comparisons, if validated against experimental data. .

Objective. Optical fibre dosimeters (OFDs) offer great promise for real-timein vivodose measurement in radiation-based treatment modalities such as radiotherapy and brachytherapy. This is attributed to their many useful qualities such as high spatial resolution and sensitivity. However, there are several requirements that an optical fibre dosimeter must meet to be acceptable for dose measurement in a specified treatment modality. In this work, the dosimetric performance of a novel optical fibre dosimeter for use in external beam radiotherapy is presented.Approach. The dosimeter was characterised for photon beam energies between 6-15 MV using a Varian TrueBeam Linac at dose rates between 100-2400 MU/min and assessed based on its repeatability, dose dependence, dose rate dependence, energy dependence and dose-per-pulse dependence.Main Results. The results demonstrated excellent short-term repeatability of 0.3%, good linearity in response (R2>0.9997), and minor dose rate dependence between 0.53%-2.49% for all beam qualities investigated. As the scintillator of the OFD is non-water equivalent, Monte-Carlo-TOPAS simulations were used to calculate the absorbed dose energy dependence. A dose-per-pulse dependence was also investigated and compared with dosimetry measurements made with an ionisation chamber and simulated from the treatment planning system. An over-response of 20%was found at the lowest investigated dose-per-pulse, and an under-response of 34%was found at the highest investigated dose-per-pulse. This is believed to be due to an intrinsic energy dependence making this type of OFD unsuitable for external beam radiotherapy dosimetry.Significance. The OFD evaluated in this work was primarily designed for high-dose-rate brachytherapy whereas this study includes the first measurements made in external beam radiotherapy and highlights the challenges of transferability of the dosimeter to a different radiation source.

Cancer has a high incidence and lethality rate, which is a significant threat to human health. With the development of high-throughput technologies, different types of cancer genomics data have been accumulated, including genomics, epigenomics, transcriptomics, proteomics, and metabolomics. A comprehensive analysis of various omics data is needed to understand the underlying mechanisms of tumor development. However, integrating such a massive amount of data is one of the main challenges today. Artificial intelligence (AI) techniques such as machine learning are now becoming practical tools for analyzing and understanding multi-omics data on diseases. Enabling great optimization of existing research paradigms for cancer screening, diagnosis, and treatment. In addition, intelligent healthcare has received widespread attention with the development of healthcare informatization. As an essential part of innovative healthcare, practical, intelligent prognosis analysis and personalized treatment for cancer patients are also necessary. This paper introduces the advanced multi-omics data analysis technology in recent years, presents the cases and advantages of the combination of both omics data and AI applied to cancer diseases, and finally briefly describes the challenges faced by multi-omics analysis and AI at the current stage, aiming to provide new perspectives for oncology research and the possibility of personalized cancer treatment.

Objective. In magneto-acousto-electrical tomography (MAET), linearly frequency-modulated (LFM) signal stimulation uses much lower peak voltage than the spike pulse stimulation, lengthening the operation life of the transducer. However, due to the uneven frequency responses of the transducer, the low-noise amplifier (LNA), and the bandpass filter (BPF), MAET using LFM signal stimulation suffers from longitudinal resolution loss. In this paper, frequency response compensated linearly frequency-modulated (FRC-LFM) signal stimulation is investigated to resolve the problem.Approach. The physical model of measurement of the frequency responses of the transducer and the cascading module of the detection electrodes, the LNA, and the BPF is constructed. The frequency responses are approximated by fitting a curve to the measurement data. The frequency response compensation function is set to the reciprocal of the product of the frequency responses. The digital FRC-LFM signal is generated in MATLAB and converted to analog signal through an arbitrary waveform generator. Two groups of MAET experiments are designed to confirm the performance of the FRC-LFM signal stimulation. Pure agar phantom with rectangular through-holes and agar phantom with pork tissue inclusion serve as the samples.Main results. The pulse-compressed magneto-acousto-electrical signal obtained using FRC-LFM stimulation has narrower main-lobe than that obtained using LFM excitation, although the signal to noise pulse interference ratio of the former is little lower than that of the latter, which is due to the limitation of the power amplifier. FRC-LFM also proves to be an effective method to utilize the frequency outside the working band of the transducer in MAET.Significance. The method in this study compensates for the longitudinal resolution loss due to the uneven frequency responses. Combining with high-capability power amplifier and high-performance LNA, the MAET using FRC-LFM signal stimulation can potentially achieve high longitudinal resolution and high sensitivity, advancing MAET toward the clinical application.

Objective. This work aims to investigate the iso-effectiveness of conventional and FLASH radiotherapy on tumors through in-silico mathematical models. We focused on the role of radiolytic oxygen depletion (ROD), which has been argued as a possible factor to explain the FLASH effect.Approach. We used a spatiotemporal reaction-diffusion model, including ROD, to simulate tumor oxygenation and response. From those oxygen distributions we obtained surviving fractions (SFs) using the linear-quadratic (LQ) model with the oxygen enhancement ratios (OERs). We then employed the calculated SFs to describe the evolution of preclinical tumor volumes through a mathematical model of tumor response, and we also extrapolated those results to calculate tumor control probabilities (TCPs) using the Poisson-LQ approach.Main results. Our study suggests that the ROD effect may cause differences in SF between FLASH and conventional radiotherapy, especially in lowα/βandpoorly oxygenatedcells. However, a statistical analysis showed that these changes in SF generally do not result in significant differences in the evolution of preclinical tumor growth curves when the sample size is small, because such differences in SF may not be noticeable in the heterogeneity of the population of animals. Nonetheless, when extrapolating this effect to TCP curves, we observed important differences between both techniques (TCP is lower in FLASH radiotherapy). When analyzing the response of tumors with heterogeneous oxygenations, differences in TCP are more important forwell oxygenatedtumors. This apparent contradiction with the results obtained for homogeneously oxygenated cells is explained by the complex interplay between the heterogeneity of tumor oxygenation, the OER effect, and the ROD effect.Significance. This study supports the experimentally observed iso-effectiveness of FLASH and conventional radiotherapy when analyzing the volume evolution of preclinical tumors (that are far from control). However, this study also hints that tumor growth curves may be less sensitive to small variations in SF than tumor control probability: ROD may lead to increased SF in FLASH radiotherapy, which while not large enough to cause significant differences in tumor growth curves, could lead to important differences in clinical TCPs. Nonetheless, it cannot be discarded that other effects not modeled in this work, like radiation-induced immune effects, can contribute to tumor control and maintain the iso-effectiveness of FLASH radiotherapy. The study of tumor growth curves may not be the ideal experiment to test the iso-effectiveness of FLASH, and experiments reporting TCP orD₅₀may be preferred.

Objective.Fluorescence molecular tomography (FMT) holds promise for early tumor detection by mapping fluorescent agents in three dimensions non-invasively with low cost. However, since ill-posedness and ill-condition due to strong scattering effects in biotissues and limited measurable data, current FMT reconstruction is still up against unsatisfactory accuracy, including location prediction and morphological preservation.Approach.To strike the above challenges, we propose a novel Sparse-Laplace hybrid graph manifold (SLHGM) model. This model integrates a hybrid Laplace norm-based graph manifold learning term, facilitating a trade-off between sparsity and preservation of morphological features. To address the non-convexity of the hybrid objective function, a fixed-point equation is designed, which employs two successive resolvent operators and a forward operator to find a converged solution.Main results.Through numerical simulations andin vivoexperiments, we demonstrate that the SLHGM model achieves an improved performance in providing accurate spatial localization while preserving morphological details.Significance.Our findings suggest that the SLHGM model has the potential to advance the application of FMT in biological research, not only in simulation but also inin vivostudies.

Objective Achieving ultra-precise temporal resolution in ionizing radiation detection is essential, particularly in positron emission tomography, where precise timing enhances signal-to-noise ratios and may enable reconstruction-less imaging. A promising approach involves utilizing ultrafast modulation of the complex refractive index, where sending probe pulses to the detection crystals will result in changes in picoseconds (ps), and thus a sub - 10 ps coincidence time resolution can be realized. Towards this goal, here, we aim to first measure the ps changes in probe pulses using an ionizing radiation source with high time resolution. Approach We used relativistic, ultrafast electrons to induce complex refractive index and use probe pulses in the near-infrared (800 nm) and terahertz (THz, 300 µm) regimes to test the hypothesized wavelength-squared increase in absorption coefficient in the Drude free-carrier absorption model. We measured BGO, ZnSe, BaF2, ZnS, PBG, and PWO with 1 mm thickness to control the deposited energy of the 3 MeV electrons, simulating ionization energy of the 511 keV photons. Main results Both with the 800 nm and THz probe pulses, transmission decreased across most samples, indicating the free carrier absorption, with an induced signal change of 11% in BaF2, but without the predicted Drude modulation increase. To understand this discrepancy, we simulated ionization tracks and examined the geometry of the free carrier distribution, attributing the mismatch in THz modulations to the sub-wavelength diameter of trajectories, despite the lengths reaching 500 µm to 1 mm. Additionally, thin samples truncated the final segments of the ionization tracks, and the measured initial segments have larger inter-inelastic collision distances due to lower stopping power (dE/dx) for high-energy electrons, exacerbating diffraction-limited resolution. Significance Our work offers insights into ultrafast radiation detection using complex refractive index modulation and highlights critical considerations in sample preparation, probe wavelength, and probe-charge carrier coupling scenarios.

Diffusion magnetic resonance imaging (dMRI) currently stands as the foremost noninvasive method for quantifying brain tissue microstructure and reconstructing white matter fiber pathways. However, the inherent free diffusion motion of water molecules in dMRI results in signal decay, diminishing the signal-to-noise ratio (SNR) and adversely affecting the accuracy and precision of microstructural data. In response to this challenge, we propose a novel method known as the Multiscale Fast Attention-Multibranch Irregular Convolutional Neural Network for dMRI image denoising. In this work, we introduce Multiscale Fast Channel Attention, a novel approach for efficient multiscale feature extraction with attention weight computation across feature channels. This enhances the model's capability to capture complex features and improves overall performance. Furthermore, we propose a multi-branch irregular convolutional architecture that effectively disrupts spatial noise correlation and captures noise features, thereby further enhancing the denoising performance of the model. Lastly, we design a novel loss function, which ensures excellent performance in both edge and flat regions. Experimental results demonstrate that the proposed method outperforms other state-of-the-art deep learning denoising methods in both quantitative and qualitative aspects for dMRI image denoising with fewer parameters and faster operational speed.