Risk factors for portal vein system thrombosis after partial splenic embolisation in cirrhotic patients with hypersplenism

AIM

To determine risk factors for portal venous system thrombosis (PVST) after partial splenic artery embolisation (PSAE) in cirrhotic patients with hypersplenism.

MATERIALS AND METHODS

Between March 2014 and February 2022, 428 cirrhotic patients with hypersplenism underwent partial splenic artery embolisation and from these patients 208 were enrolled and 220 were excluded. Medical records of enrolled patients were collected. Computed tomography (CT) images were reviewed by two blinded, independent radiologists. Statistical analyses were performed by using SPSS.

RESULTS

Progressive PVST was observed in 18.75% (39/208) of cirrhotic patients after PSAE. No significant differences in peripheral blood counts, liver function biomarkers, and renal function were observed between the patients with progressive PVST and the patients without progressive PVST. The imaging data showed significant differences in PVST, the diameters of the portal, splenic, and superior mesenteric veins between the progressive PVST group and non-progressive PVST group. Univariate and multivariate analysis demonstrated portal vein thrombosis, spleen infarction percentage, and the diameter of the splenic vein were independent risk factors for progressive PVST. Seventeen of 173 (9.83%) patients showed new PVST; the growth of PVST was observed in 62.86% (22/35) of the patients with pre-existing PVST. Spleen infarction percentage and the diameter of the splenic vein were independent risk factors for new PVST after PSAE.

CONCLUSION

The present study demonstrated portal vein thrombosis, spleen infarction percentage, and the diameter of the splenic vein were independent risk factors for PVST after PSAE in cirrhotic patients with hypersplenism.

Test of CP Symmetry in Hyperon to Neutron Decays

The quantum entangled J/ψ→Σ^{+}Σ[over ¯]^{-} pairs from (1.0087±0.0044)×10^{10}  J/ψ events taken by the BESIII detector are used to study the nonleptonic two-body weak decays Σ^{+}→nπ^{+} and Σ[over ¯]^{-}→n[over ¯]π^{-}. The CP-odd weak decay parameters of the decays Σ^{+}→nπ^{+} (α_{+}) and Σ[over ¯]^{-}→n[over ¯]π^{-} (α[over ¯]_{-}) are determined to be 0.0481±0.0031_{stat}±0.0019_{syst} and -0.0565±0.0047_{stat}±0.0022_{syst}, respectively. The decay parameter α[over ¯]_{-} is measured for the first time, and the accuracy of α_{+} is improved by a factor of 4 compared to the previous results. The simultaneously determined decay parameters allow the first precision CP symmetry test for any hyperon decay with a neutron in the final state with the measurement of A_{CP}=(α_{+}+α[over ¯]_{-})/(α_{+}-α[over ¯]_{-})=-0.080±0.052_{stat}±0.028_{syst}. Assuming CP conservation, the average decay parameter is determined as ⟨α_{+}⟩=(α_{+}-α[over ¯]_{-})/2=-0.0506±0.0026_{stat}±0.0019_{syst}, while the ratios α_{+}/α_{0} and α[over ¯]_{-}/α[over ¯]_{0} are -0.0490±0.0032_{stat}±0.0021_{syst} and -0.0571±0.0053_{stat}±0.0032_{syst}, where α_{0} and α[over ¯]_{0} are the decay parameters of the decays Σ^{+}→pπ^{0} and Σ[over ¯]^{-}→p[over ¯]π^{0}, respectively.

Measurement of Energy-Dependent Pair-Production Cross Section and Electromagnetic Form Factors of a Charmed Baryon

We study the process e^{+}e^{-}→Λ_{c}^{+}Λ[over ¯]_{c}^{-} at twelve center-of-mass energies from 4.6119 to 4.9509 GeV using data samples collected by the BESIII detector at the BEPCII collider. The Born cross sections and effective form factors (|G_{eff}|) are determined with unprecedented precision after combining the single and double-tag methods based on the decay process Λ_{c}^{+}→pK^{-}π^{+}. Flat cross sections around 4.63 GeV are obtained and no indication of the resonant structure Y(4630), as reported by Belle, is found. In addition, no oscillatory behavior is discerned in the |G_{eff}| energy dependence of Λ_{c}^{+}, in contrast to what is seen for the proton and neutron cases. Analyzing the cross section together with the polar-angle distribution of the Λ_{c}^{+} baryon at each energy point, the moduli of electric and magnetic form factors (|G_{E}| and |G_{M}|) are extracted and separated. For the first time, the energy dependence of the form factor ratio |G_{E}/G_{M}| is observed, which can be well described by an oscillatory function.

Observation of New Baryons in the Ξ_{b}^{-}π^{+}π^{-} and Ξ_{b}^{0}π^{+}π^{-} Systems

The first observation and study of two new baryonic structures in the final state Ξ_{b}^{0}π^{+}π^{-} and the confirmation of the Ξ_{b}(6100)^{-} state in the Ξ_{b}^{-}π^{+}π^{-} decay mode are reported using proton-proton collision data collected by the LHCb experiment, corresponding to an integrated luminosity of 9  fb^{-1}. In addition, the properties of the known Ξ_{b}^{*0}, Ξ_{b}^{'-} and Ξ_{b}^{*-} resonances are measured with improved precision. The new decay mode of the Ξ_{b}^{0} baryon to the Ξ_{c}^{+} π^{-} π^{+} π^{-} final state is observed and exploited for the first time in these measurements.

Determination of Spin-Parity Quantum Numbers for the Narrow Structure near the pΛ[over ¯] Threshold in e^{+}e^{-}→pK^{-}Λ[over ¯]+c.c

A narrow structure in the pΛ[over ¯] system near the mass threshold, named as X(2085), is observed in the process e^{+}e^{-}→pK^{-}Λ[over ¯] with a statistical significance greater than 20σ. Its spin and parity are determined for the first time to be J^{P}=1^{+} in an amplitude analysis, with a statistical significance greater than 5σ over other quantum numbers (0^{-},1^{-} and 2^{+}). The pole positions of X(2085) are measured to be M_{pole}=(2084_{-2}^{+4}±9)  MeV and Γ_{pole}=(58_{-3}^{+4}±25)  MeV, where the first uncertainties are statistical and the second ones are systematic. The analysis is based on the study of the process e^{+}e^{-}→pK^{-}Λ[over ¯] and uses the data samples collected with the BESIII detector at the center-of-mass energies sqrt[s]=4.008, 4.178, 4.226, 4.258, 4.416, and 4.682 GeV with a total integrated luminosity of 8.35  fb^{-1}.

Measurement of the Λ_{b}^{0}→Λ(1520)μ^{+}μ^{-} Differential Branching Fraction

The branching fraction of the rare decay Λ_{b}^{0}→Λ(1520)μ^{+}μ^{-} is measured for the first time, in the squared dimuon mass intervals q^{2}, excluding the J/ψ and ψ(2S) regions. The data sample analyzed was collected by the LHCb experiment at center-of-mass energies of 7, 8, and 13 TeV, corresponding to a total integrated luminosity of 9  fb^{-1}. The result in the highest q^{2} interval, q^{2}>15.0  GeV^{2}/c^{4}, where theoretical predictions have the smallest model dependence, agrees with the predictions.

Precise Measurement of the e^{+}e^{-}→D_{s}^{*+}D_{s}^{*-} Cross Sections at Center-of-Mass Energies from Threshold to 4.95 GeV

The process e^{+}e^{-}→D_{s}^{*+}D_{s}^{*-} is studied with a semi-inclusive method using data samples at center-of-mass energies from threshold to 4.95 GeV collected with the BESIII detector operating at the Beijing Electron Positron Collider. The Born cross sections of the process are measured for the first time with high precision in this energy region. Two resonance structures are observed in the energy-dependent cross sections around 4.2 and 4.4 GeV. By fitting the cross sections with a coherent sum of three Breit-Wigner amplitudes and one phase-space amplitude, the two significant structures are assigned masses of (4186.8±8.7±30) and (4414.6±3.4±6.1)  MeV/c^{2}, widths of (55±15±53) and (122.5±7.5±8.1)  MeV, where the first errors are statistical and the second ones are systematic. The inclusion of a third Breit-Wigner amplitude is necessary to describe a structure around 4.79 GeV.

[Analysis of factors affecting the quality of donor corneal endothelial cells]

Objective: To explore the factors influencing the quality of donor corneal endothelium. Methods: A retrospective case series study was conducted. Data from 568 donor corneas obtained from the Shandong Eye Bank between July 1, 2020, and June 30, 2021, were collected for analysis. The corneal endothelium of the donor corneas was observed using corneal endothelial microscopy to assess corneal endothelial cell density (ECD), coefficient of variation, and hexagonal cell ratio (HEX). Relevant factors of corneal donors were collected, including gender, age, cause of death, season of death, time from death to corneal retrieval, and methods of corpse preservation, to investigate their impact on the quality of donor corneal endothelium. The age factor was divided into five age groups: 0-20 years, 21-40 years, 41-60 years, 61-80 years, and >80 years. The time of corneal retrieval was divided into three periods based on the time elapsed since the donor's death: <6 hours, 6-12 hours, and >12 hours. The relationship between these factors and corneal endothelial conditions was analyzed. Results: The 568 donor corneas were obtained from 288 donors, including 225 males (78.13%) and 63 females (21.87%). The mean age was 51.77±18.48 years. The causes of death among donors were as follows: cardiovascular diseases 54.58% (275 individuals), cancer 17.96% (74 individuals), organ failure 14.26% (49 individuals), and accidents 13.20% (64 individuals). The mean time of corneal retrieval after donor death was 140 (76, 400) minutes (ranging from 30 minutes to 45 hours). Among the 145 corneas (25.53%) that had their initial corneal endothelial microscopy examination, the images were not clear, and after thorough rewarming, 106 corneas (18.7%) still had unclear images and could not be analyzed. Among the 462 corneas (81.3%) with clear images, the ECD was (2 602.23±318.40) cells/mm², the coefficient of variation was 36.61%±4.81%, and the HEX was 52.73%±7.15%. The ECD of corneas from older donors was lower compared to younger donors, and the differences between age groups were statistically significant (P<0.001). Corneas from donors who died due to accidents had a higher ECD [(2 829.88±313.90) cells/mm²] compared to those who died from cancer, cardiovascular diseases, and organ failure, and the differences were statistically significant (P<0.001). The ECD was highest when corneas were retrieved within 6 hours after death, and the difference was statistically significant (P<0.001). Older donors had higher coefficients of variation but lower HEX values (both P<0.05). Corneas retrieved after a longer time from death had higher coefficients of variation, and the difference was statistically significant (P<0.05), but there was no statistically significant difference in HEX (P>0.05). Organ failure, cryopreservation, and corneal retrieval time >12 hours were risk factors for unclear corneal endothelial imaging (all P<0.001). Among the 136 corneal endothelial images (23.94%), circular, oval, or band-shaped dark areas were observed, and corneas with dark areas had lower ECD (P<0.05). The longer the time elapsed from death to corneal retrieval, the more dark areas were observed (P<0.001). The presence of dark areas did not affect the coefficient of variation and HEX (P>0.05). Conclusion: Advanced donor age, death due to chronic diseases, longer time elapsed from death to corneal retrieval, and cryopreservation of the body lead to a decrease in the quality of donor corneal endothelium.

First Experimental Study of the Purely Leptonic Decay D_{s}^{*+}→e^{+}ν_{e}

Using 7.33  fb^{-1} of e^{+}e^{-} collision data taken with the BESIII detector at the BEPCII collider, we report the first experimental study of the purely leptonic decay D_{s}^{*+}→e^{+}ν_{e}. Our data contain a signal of this decay with a statistical significance of 2.9σ. The branching fraction of D_{s}^{*+}→e^{+}ν_{e} is measured to be (2.1_{-0.9_{stat}}^{+1.2}±0.2_{syst})×10^{-5}, corresponding to an upper limit of 4.0×10^{-5} at the 90% confidence level. Taking the total width of the D_{s}^{*+} [(0.070±0.028)  keV] predicted with the radiative D_{s}^{*+} decay from the lattice QCD calculation as input, the decay constant of the D_{s}^{*+} is determined to be f_{D_{s}^{*+}}=(214_{-46_{stat}}^{+61}±44_{syst})  MeV, corresponding to an upper limit of 354 MeV at the 90% confidence level.

Towards Clinical Proton Minibeam Radiation Therapy (pMBRT): Development of Clinical pMBRT System Prototype and pMBRT-Specific Treatment Planning Method

PURPOSE/OBJECTIVE(S)

Proton minibeam radiation therapy (pMBRT) is an emerging spatially fractionated RT (SFRT) modality that can provide very high therapeutic index compared to conventional radiotherapy methods and clinically-available SFRT methods (GRID and LATTICE). The biological data collected thus far encourage the preparation of clinical trials in pMBRT. This work is to facilitate the clinical translation of pMBRT by developing (1) the first clinical pMBRT system prototype worldwide, readily available for both small-animal biology studies and large-animal pMBRT trials; (2) pMBRT-specific treatment planning method with peak-valley dose ratio (PVDR) optimization capability for large animal and patient pMBRT trials, which is currently unavailable.

MATERIALS/METHODS

The pMBRT system is based on clinically-used pencil-beam-scanning proton machine equipped with large-field clinical-size pMBRT collimator, pMBRT-dedicated treatment planning system, and KV/CBCT imaging guidance. The multi-slit brass collimator has 10 × 10 cm field size, 0.4mm width per slit and 4 mm center-to-center distance. The divergence of slits is tailed to the divergence of the proton beam. A unique universal collimator design is implemented, so that we can keep the outer fitting to the snout and conveniently inter-change collimators as needed. The pMBRT-specific treatment planning method jointly optimizes PVDR and dose objectives, to meet a minimal PVDR threshold and maximize PVDR, to avoid the situations where meeting dose objectives can compromise PVDR when PVDR were not optimized. In addition, the survival fraction for organs at risk is also optimized. The dose calculation engine is based on the Monte Carlo method using TOPAS. The optimization algorithm utilizes total variation and L1 sparsity regularization to maximize PVDR and iterative convex relaxation method to solve the optimization problem.

RESULTS

Monte Carlo simulations via TOPAS were performed to design this large-field multi-slit collimator, using the beam structure and the beam data specific to our proton system, with mean dose rate of 8 Gy/min under clinical condition with the collimator in place. The feasibility of using this pMBRT system for small-animal studies has been demonstrated, with customized 3D printed holder for immobilizing small animals, on-board KV imaging system for accurate small-animal positioning, and the GAFchromic film for verifying radiation dose and PVDR. On the other hand, the efficacy of pMBRT-specific treatment planning method (with PVDR optimization capability) to improve PVDR has been demonstrated using retrospective patient planning studies in comparison with standard proton treatment planning method (without PVDR optimization capability).

CONCLUSION

The initial development of a clinical pMBRT system prototype and pMBRT-specific treatment planning method of PVDR optimization capability has been completed with ongoing efforts to make this system ready for large-animal pMBRT studies.

SEA-Net: Structure-Enhanced Attention Network for Limited-Angle CBCT Reconstruction of Clinical Projection Data

PURPOSE/OBJECTIVE(S)

Limited-angle CBCT (LA-CBCT) is of great clinical interest, because the scanning time and the patient dose are proportional to the scanning range of gantry rotation angles of CBCT. However, the image reconstruction for LA-CBCT remains technically challenging, which suffers from severe wedge artifacts and image distortions. This work aims to improve LA-CBCT by developing deep learning (DL) methods for real clinical CBCT projection data, which is the first feasibility study of clinical-projection-data-based LA-CBCT, to the best of our knowledge.

MATERIALS/METHODS

Targeting at real clinical projection data, we have explored various DL methods such as image/data/hybrid-domain methods and finally developed a so-called Structure-Enhanced Attention Network (SEA-Net) method that has the best image quality from clinical projection data among the DL methods we have implemented. Specifically, the proposed SEA-Net employs a specialized structure enhancement sub-network to promote texture preservation. Based on the observation that the distribution of wedge artifacts in reconstruction images is non-uniform, the spatial attention module is utilized to emphasize the relevant regions while ignores the irrelevant ones, which leads to more accurate texture restoration.

RESULTS

SEA-Net was validated in comparison with analytic (FDK), iterative (TV), image-domain DL (DDNet and FED-INet, data-domain DL (DCAR), dual-domain DL (Sam'Net), and various unrolling DL (hdNet, CTNet, FSR-Net, CasRedSCAN) methods. Among all methods, the SEA-Net had the best image reconstruction quality as quantified by root-mean-square error (RMSE), peak signal-to-noise ratio (PSNR), and structural similarity index (SSIM), for various LA-CBCT problems of 90°-180° projection data. In addition, LA-CBCT via SEA-Net provided comparable accuracy for both patient setup (quantified by image registration accuracy from planning CT (pCT) to CBCT) and dose calculation (see the table), with full-view CBCT.

CONCLUSION

We explored various DL methods and developed an image-domain-based method termed SEA-Net that provided the best image quality for clinical projection data. To the best of our knowledge, this is the first feasibility study of the real clinical-projection-data-based LA-CBCT. Moreover, LA-CBCT via SEA-Net can potentially provide comparable accuracy for patient setup and dose calculation, with full-view CBCT.

Investigating the Role of MK2 in Head and Neck Squamous Cell Carcinoma Growth, Metastasis and STING Pathway Activation

PURPOSE/OBJECTIVE(S)

Our prior work demonstrated that inhibition of MAPKAPK2 (MK2) can enhance radiation (RT)-mediated in vivo head and neck squamous cell carcinoma (HNSCC) tumor control and survival in preclinical immune incompetent models. The cytosolic DNA sensor cyclic GMP-AMP synthetase (cGAS) and its downstream adaptor protein, stimulator of interferon genes (STING), are conserved proteins within the innate immune signaling pathways and are important for mediating host defense against microbial infection and can play a role in anti-cancer immunity. We hypothesized that loss of MK2 enhances radiation-induced cGAS-STING pathway activation leading to improved tumor control and survival.

MATERIALS/METHODS

MK2 shRNA knockdown human (Tu167, CAL27) and MK2 Cas9/CRISPR knockout (KO) syngeneic murine (Ly2, MLM3) HNSCC cell lines were treated with 10 Gy irradiation. Micronuclei were quantitated by DAPI-immunofluorescence (IF). Protein changes in cGAS-STING were evaluated by immunoblot. Inflammatory cytokine production including Type I IFNβ1 were evaluated by RT-qPCR. Ly2 and MLM3 cells were orthotopically or flank engrafted into immune competent mice (Balb/c, C57Bl/6, respectively) for animal tumor control-survival studies. Tumor immune cell infiltrate was examined using FACS and immunohistochemistry. Selected drug studies using the MK2 inhibitor, ATI-450, were performed with RT.

RESULTS

Loss of MK2 in HNSCC (Tu167, CAL27, Ly2, MLM3) treated with RT led to a significant increase in micronuclei formation compared to control cells. MK2-enhanced micronuclei generation following RT could be inhibited with the actin filament polymerization inhibitor, cytochalasin B. RT treatment of MK2 shRNA cells led to increased cGAS and phospho-STING levels compared to either treatment alone. IFNβ1 levels were significantly higher in HNSCC cell lines treated with RT and with MK2 inhibited by an MK2 inhibitor (ATI-450) or genetic reduction compared to either treatment alone. In-vivo implantation of MLM3 cells into C57Bl/6 comparing control vs MK2 KO tumors treated with ±8 Gy RT demonstrated improved mouse survival favoring RT+MK2 KO over RT, MK2 KO or parental (63, 58, 58.5, 35 days, respectively). FACS analysis of MLM3 WT v KO tumors 3 days post RT showed an overall increase in the number of CD3/CD8 T-cells infiltrating into the tumor in all groups except for parental tumors. Further analysis demonstrated that loss of MK2 reversed CD8 T-cell exhaustion and when combined with radiotherapy led to increased CD8 T-cell activation. Furthermore, activated CD4 and CD8 T cells were reduced in WT+RT cells compared to WT tumors whereas no reduction was seen in the KO or KO+RT.

CONCLUSION

HNSCC tumor MK2 inhibition enhances RT-mediated micronuclei formation and subsequent cGAS-STING-IFNβ1 levels. Loss of HNSCC MK2 leads to increased CD4-CD8 T-cell infiltration into the tumor and this effect is enhanced following RT. Targeting tumor MK2 may facilitate improved tumor control.

Proton LET Optimization Via Iterative Convex Relaxation Method

PURPOSE/OBJECTIVE(S)

If not optimized, the LET distribution can greatly impact normal tissue toxicity near tumor targets, because the LET often peaks at the distal edge of Bragg peak. LET optimization can account for biological effectiveness of protons during treatment planning, for minimizing biological proton dose and hot spots to normal tissues. However, the LET optimization is nonlinear and nonconvex to solve, which poses a great challenge in optimization. This work will develop an effective LET optimization method via iterative convex relaxation (ICR).

MATERIALS/METHODS

In contrast to the generic nonlinear optimization method, such as Quasi-Newton (QN) method, that does not account for specific characteristics of LET optimization, ICR is tailored to LET modeling and optimization in order to effectively and efficiently solve the LET problem. Specifically, nonlinear dose-averaged LET term is iteratively linearized and becomes convex during ICR, while nonconvex dose-volume constraint and minimum-monitor-unit constraint are also handled by ICR, so that the solution for LET optimization is obtained by solving a sequence of convex and linearized convex subproblems. Since the high LET mostly occurs near the target, a 1cm normal-tissue expansion of clinical target volume (CTV) (excluding CTV), i.e., CTV1cm, is defined to as an auxiliary structure during treatment planning, where LET is minimized.

RESULTS

ICR was validated in comparison with QN for abdomen, lung, and head-and-neck (HN) cases. ICR was effective for LET optimization, as ICR substantially reduced the LET and biological dose in CTV1cm the ring, with preserved dose conformality to CTV. Compared to QN, ICR had smaller LET, physical and biological dose in CTV1cm, and higher conformity index values; ICR was also computationally more efficient, which was about 3 times faster than QN. A lung case is presented in the table, where the quantities from top to bottom are computational time T (unit: second); total objective F, dose objective Fd and LET objective FL (unit: 10-3); conformity index for physical dose CId and biological dose CIb; mean LET L (unit: keV/μm), mean physical dose d and mean biological dose b (in ratio to prescription dose) for CTV1cm; mean LET L (unit: keV/μm), mean physical dose d and mean biological dose b (in ratio to prescription dose; unit: 10-1) for the heart.

CONCLUSION

A LET-specific optimization method based on ICR has been developed for solving proton LET optimization, which has been shown to be more computationally efficient than generic nonlinear optimizer via QN, with better plan quality in terms of LET, biological and physical dose conformality.

Virtual-Collimator Based Spatial Dose Modulation for Proton GRID Therapy

PURPOSE/OBJECTIVE(S)

Compared to conventional proton therapy, the proton GRID therapy can substantially improve normal tissue protection (with the delivery of spatially-modulated peak-valley dose pattern to normal tissues) while maintaining the tumor control efficacy (with the delivery of uniform dose pattern to tumor targets). The realization of proton GRID often relies on the use of physical collimators to shape the spatial dose distribution. However, the physical collimator may increase neutron dose, decrease delivery efficiency, and limit the freedom for patient positioning. Here we propose a virtual-collimator (VC) method for proton GRID. This new approach can generate peak-to-valley pattern with high peak-to-valley dose ratio (PVDR), without using a physical collimator.

MATERIALS/METHODS

The principle behind the VC method to modulate the spatial dose distribution consists of two major steps: (1) the primary beam is essentially halved, i.e., the beamlets are interleaved, so that the organ-at-risk (OAR) plane has the peak-valley dose pattern, while the target plane also has the valley dose; (2) the complementary beam is added with half complementary beamlets to fill in the previously valley-dose positions at the target plane, so that the target dose is uniform, while on the other hand, the complementary beam is angled slightly from the primary beam, so that the OAR still has the peak-valley dose pattern. Moreover, on top of VC, we also utilize sparsity regularization method using total variation and L1 sparsity (TVL1) to further jointly optimize PVDR and dose objectives, namely VC-TVL1.

RESULTS

VC and VC-TVL1 were validated in comparison with conventional proton GRID treatment planning method via IMPT ("CONV") and TVL1-based proton GRID treatment planning method without VC ("TVL1"), for a prostate case with single-beam (270° only) or two-beam (90° and 270°) scenarios. As shown in the table, the results show that VC can indeed modulate spatial dose with higher PVDR than CONV or even TVL1. VC had higher spatial modulation frequency with smaller peak-to-peak distance than TVL1. Moreover, VC+TVL1, as the synergy of VC and TVL1, further improved PVDR from VC or TVL1 alone.

CONCLUSION

A new way to deliver proton GRID therapy without a physical collimator is developed using the VC method. The VC method can be synergized with TVL1 optimization algorithm to further jointly optimize PVDR and dose objectives.

2V-CBCT: Two-Orthogonal-Projection Based CBCT Reconstruction and Dose Calculation from Real CBCT Projection Data

PURPOSE/OBJECTIVE(S)

Not all radiation therapy (RT) treatments/fractions have CBCT acquired, but two orthogonal projections (i.e., KV radiography) are always available. This work demonstrates the feasibility of two-orthogonal-projection-based CBCT (2V-CBCT) reconstruction and dose calculation for RT from real CBCT projection data, which is the first 2V-CBCT feasibility study using real projection data, to the best of our knowledge.

MATERIALS/METHODS

2V-CBCT is a severely ill-posed inverse problem for which we propose a coarse-to-fine learning strategy. First, a 3D deep neural network that can extract and exploit the inter-slice and intra-slice information is adopted to predict the initial 3D volumes. Then, a 2D deep neural network is utilized to fine-tune the initial 3D volumes slice-by-slice. During the fine-tuning stage, a perceptual loss based on multi-frequency features is employed to enhance the image reconstruction. Dose calculation results from both photon and proton RT demonstrate that 2V-CBCT provides comparable accuracy with full-view CBCT based on real projection data.

RESULTS

The proposed method was evaluated on real HN data acquired from on-board CBCT scanners rather than the low-resolution simulated data or down-sampled data. Both visual assessment and quantitative analysis demonstrate that the proposed coarse-to-fine learning strategy has the potential to produce satisfactory volumetric images from two orthogonal projections. Furthermore, we assessed the utility of 2V-CBCT in RT. The results show that the dose distribution maps, dose-volume histograms, and dose parameters calculated using 2V-CBCT have comparable accuracy with the counterparts calculated using the corresponding full-view CBCT for both photon and proton RT. In the table, the methods under comparison are pCT (planning CT), FV-CBCT (CBCT reconstructed with full-view projection data), and 2V-CBCT (CBCT reconstructed with two orthogonal projections).

CONCLUSION

A new effective 2V-CBCT reconstruction method is proposed and validated using real CBCT projection data, which can potentially provide comparable dose calculation accuracy for both photon and proton RT.

A New Treatment Planning Method for Efficient Proton ARC Therapy with Direct Minimization of Number of Energy Jumps

PURPOSE/OBJECTIVE(S)

The optimization of energy layer distributions is crucial for efficient proton ARC therapy: on one hand, a sufficient number of energy layers is needed to ensure the plan quality; on the other hand, an excess number of energy jumps can substantially slow down the treatment delivery. This work will develop a new treatment plan optimization method with direct minimization of number of energy jumps (NEJ), which will be shown to outperform state-of-the-art methods in both plan quality and delivery efficiency.

MATERIALS/METHODS

The proposed method jointly optimizes the plan quality and minimizes the NEJ. To minimize NEJ, (1) the proton spots x is summed per energy layer to form the energy vector y; (2) y is binarized via sigmoid transform into y1; (3) y1 is multiplied with a predefined energy order vector via dot product into y2; (4) y2 is filtered through the finite-differencing kernel into y3 in order to identify NEJ; (5) only the NEJ of y3 is penalized, while x is optimized for plan quality. The solution algorithm to this new method is based on iterative convex relaxation.

RESULTS

The new method is validated in comparison with state-of-the-art methods called energy sequencing (ES) method and energy matrix (EM) method. In terms of delivery efficiency, the new method had fewer NEJ, less energy switching time, and generally less total delivery time. In terms of plan quality, the new method had smaller optimization objective values, lower normal tissue dose, and generally better target coverage. A head-and-neck case is provided in the table with the following dosimetric parameters: planning objective value F; conformity index CI; homogeneity index HI; mean dose of larynx DOAR; mean body dose Dbody; the unit of dose is Gy.

CONCLUSION

We have developed a new treatment plan optimization method with direct minimization of NEJ, and demonstrated that this new method outperformed state-of-the-art methods (ES and EM) in both plan quality and delivery efficiency.

Vertex Position Optimization for LATTICE Therapy

PURPOSE/OBJECTIVE(S)

LATTICE radiation therapy (RT) aims to deliver 3D heterogenous dose of high peak-to-valley dose ratio (PVDR) to the tumor target, with peak dose at lattice vertices inside the target and valley dose for the rest of the target. In current clinical practice the lattice vertex positions are constant during treatment planning. This work proposes a new LATTICE plan optimization method that can optimize lattice vertex positions as plan variables, which is the first lattice vertex position optimization study to the best of our knowledge.

MATERIALS/METHODS

The new LATTICE treatment planning method optimizes lattice vertex positions as well as other plan variables (e.g., photon fluences or proton spot weights), with optimization objectives for target PVDR and organs-at-risk (OAR) sparing. To satisfy mathematical differentiability, the lattice vertices are approximated in sigmoid functions. For geometric feasibility, proper geometry constraints are enforced onto the lattice vertex positions. The lattice vertex position optimization problem is solved by iterative convex relaxation method, where lattice vertex positions and photon/proton plan variables are jointly updated via the Quasi-Newton method.

RESULTS

Both photon and proton LATTICE RT were considered, and the optimal lattice vertex positions in terms of plan objectives were found by solving all possible combinations on given discrete positions via heuristic searching based on standard IMRT/IMPT, which served as the ground truth for validating the new LATTICE method ("NEW"). That is, the plan with the smallest optimization objective ("BEST"), the plan with the median optimization objective ("MID"), and the plan with the largest optimization objective ("WORST") were selected as the reference plans to be compared with NEW. The table was for an abdomen case with the large bowel as the OAR, where the parameters are total optimization objective f, the mean valley dose of target Dvalley, the mean peak dose of target Dpeak, PVDR = Dpeak/Dvalley, and the mean dose of large bowel Dbowel. The unit of doses is Gy. The results in the table show that the new method indeed provided the optimal lattice vertex positions with the smallest optimization objective, the largest target PVDR, and the best OAR sparing.

CONCLUSION

A new LATTICE treatment planning method is proposed and validated that can optimize lattice vertex positions as well as other photon or proton plan variables for improving target PVDR and OAR sparing.

An Effective Dose Rate Optimization Algorithm for Efficient Conventional-Dose-Rate Proton Therapy and Ultra-High-Dose-Rate FLASH Proton Therapy

PURPOSE/OBJECTIVE(S)

The treatment delivery with high dose rate is desirable for efficient conventional-dose-rate proton therapy (IMPT and proton ARC) and essential for achieving biological sparing effect to normal tissues for ultra-high-dose-rate FLASH proton therapy. The minimum-monitor-unit (MMU) problem is mathematically fundamental to these high-dose-rate proton therapy problems, since the dose rate is proportionally associated with the MMU threshold. However, the MMU problem is nonconvex and larger MMU threshold that leads to higher dose rate makes the nonconvex problem more difficult to solve, for which the conventional optimization methods (e.g., alternating direction method of multipliers (ADMM) or proximal gradient descent (PGD) usually cannot satisfactorily solve. This work will develop an effective dose rate optimization algorithm for the MMU problem.

MATERIALS/METHODS

The new method is based on orthogonal matching pursuit (OMP), and consists of two essential components. First, the iterative convex relaxation method is used to determine the active sets for dose-volume planning constraints and decouple the MMU constraint from the rest. Second, a modified OMP optimization algorithm is used to handle the MMU constraint: the non-zero spots are greedily selected via OMP to form the solution set to be optimized, and then a convex constrained subproblem is formed and can be conveniently solved to optimize the spot weights restricted to this solution set via OMP. During this iterative process, the new non-zero spots localized via OMP will be adaptively added to or removed from the optimization objective.

RESULTS

The new method via OMP is validated in comparison with conventional methods (ADMM and PGD) and a recently proposed new method (called stochastic coordinate descent (SCD) method) for high-dose-rate IMPT, ARC, and FLASH problems of large MMU thresholds, and the results suggest that OMP substantially improved the plan quality from PGD, ADMM and SCD in terms of both target dose conformality (e.g., quantified by max target dose and conformity index) and normal tissue sparing (e.g., mean and max dose). A brain case is provided in the table: compared to PGD/ADMM/SCD, OMP improved the conformity index from 0.42/0.52/0.33 to 0.65 for IMPT and 0.46/0.60/0.61 to 0.83 for ARC.

CONCLUSION

A new dose rate optimization algorithm via OMP is developed to solve the MMU problem with large MMU thresholds, and validated using examples of IMPT, ARC, and FLASH with substantially improved plan quality from conventional methods (ADMM and PGD) and a recently proposed new method SCD.

[Application of intraoperative optical coherence tomography in deep lamellar keratoplasty]

Objective: To evaluate the clinical application value of intraoperative optical coherence tomography (iOCT) in deep anterior lamellar keratoplasty (DALK) using the big-bubble technique to bare Descemet's membrane. Methods: Retrospective case series. Clinical data of 92 patients (92 eyes) with monocular stromal corneal diseases who underwent big-bubble DALK in the Eye Hospital of Shandong First Medical University from January 2020 to August 2021 were collected. There were 53 males and 39 females. The average age was (53.2±16.0) years old. All patients underwent iOCT scanning to determine the location and depth of the injection needle after initial removal of the corneal lesion, to observe the integrity of the recipient bed, Descemet's membrane, after complete lesion removal, and to observe the adhesion between the corneal graft and the recipient bed and check folds on the recipient bed after suturing of the corneal graft. The intraoperative perforation of Descemet's membrane, postoperative thickness of the cornea and the recipient bed, visual acuity, and corneal astigmatism were recorded. Results: By iOCT, the thickness of the recipient bed was found to be about 1/2 of the corneal thickness and relatively uniform in all directions in 62 eyes (67.4%), so the sterile air was injected from the center of the recipient bed to separate it from the stromal layer. In 30 eyes (32.6%) with an uneven thickness of the recipient bed, the sterile air was injected from the paracentral area of the recipient bed. Under the guidance of iOCT scanning, 89 eyes (96.7%) did not experience any perforation of Descemet's membrane during surgery. The Descemet's membrane folds in the central 5-mm area of the recipient bed was observed and flattened in 20 eyes with the assistance of iOCT scanning. The postoperative corneal thickness was (578.95±108.26) μm, and the recipient bed thickness was (36.06±23.11) μm. The best corrected visual acuity of all patients at 6 months after surgery was 0.57±0.25 logMAR, which was significantly better than that before surgery (1.61±1.27 logMAR; P<0.001). The average corneal astigmatism at 6 months after surgery was (2.72±2.44) diopters. Conclusions: The application of iOCT scanning in DALK surgery assisted by the big-bubble method can provide safe guidance for surgeons to adopt correct surgical procedures, decrease the risk of Descemet's membrane perforation, reduce the recipient bed folds, and facilitate corneal interlayer adhesion, thereby improving the visual prognosis.

[Predictive value of low-density lipoprotein cholesterol/remnant cholesterol coordination on cardiovascular and cerebrovascular events in elderly patients with hypertension]

Objective: To explore the predictive value of serum low-density lipoprotein cholesterol/residual cholesterol (LDL-C/RC) coordination on major adverse cardiovascular events (MACE) in elderly hypertensive patients. Methods: This is a prospective cohort study. Elderly hypertensive patients hospitalized in Beijing Anzhen Hospital from June 2018 to June 2020 were prospectively enrolled. According to the coordination of baseline LDL-C/RC, patients were divided into four groups: low LDL-C and low RC (LDL-C<2.6 mmol/L, and RC<0.62 mmol/L), low LDL-C and high RC (LDL-C<2.6 mmol/L, and RC≥0.62 mmol/L), high LDL-C and low RC (LDL-C≥2.6 mmol/L, and RC<0.62 mmol/L), and high LDL-C and high RC (LDL-C≥2.6 mmol/L, and RC≥0.62 mmol/L). Patients were followed up until June 2022. The primary outcome was MACE, including non-fatal acute coronary syndrome, non-fatal acute heart failure, non-fatal stroke and death. Kaplan-Meier survival analysis was used to evaluate MACE among the four groups, and Cox regression analysis was performed to evaluate the related factors of MACE. Results: A total of 847 hypertensive patients were enrolled. There were 453 males (53.5%), mean age was (72.4±8.8) years. There were 207, 162, 263, 215 patients in the low LDL-C and low RC, low LDL-C and high RC, high LDL-C and low RC and high LDL-C and high RC group, respectively. During a median follow-up of 37 months, 196 patients (23.1%) had MACE, including 77 cases (9.1%) of ACS, 56 cases (6.6%) of acute heart failure, 34 cases (4.0%)of non-fatal stroke, 29 cases(3.4%) of all-cause mortality, and 16 cases (1.9%) of cardiovascular death. Risk of MACE of the four groups was significantly different (log-rank P<0.001). Compared with low LDL-C and low RC groups, high LDL-C and high RC groups had the highest incidence of MACE (HR=2.237, 95%CI 1.328-3.783, P=0.004), followed by low LDL-C and high RC groups (HR=1.745, 95%CI 1.220-2.527, P=0.003) and high LDL-C and low RC groups (HR=1.393, 95%CI1.048-1.774, P=0.022). In addition, the risk of nonfatal ACS among the four groups was also statistically significant (P=0.037), while the risk of nonfatal acute heart failure, nonfatal stroke, all-cause mortality and cardiovascular death were similar (all P>0.05). Multivariate Cox regression analysis showed that age (HR=1.271), duration of hypertension (HR=1.339), diabetes (HR=1.415), hyperlipidemia (HR=1.348), serum creatinine (HR=1.263), N-terminal pro-B-type natriuretic peptide (HR=1.316), LDL-C (HR=1.205), RC (HR=1.302), low LDL-C and high RC (HR=1.745), high LDL-C and low RC (HR=1.393), high LDL-C and high RC (HR=2.237) were independently associated with the occurrence of MACE. Conclusion: The coordination of LDL-C/RC affects the risk of MACE in elderly hypertensive patients, and the risk of MACE is the highest in patients with high LDL-C and high RC.

First Study of Reaction Ξ^{0}n→Ξ^{-}p Using Ξ^{0}-Nucleus Scattering at an Electron-Positron Collider

Using (1.0087±0.0044)×10^{10}  J/ψ events collected with the BESIII detector at the BEPCII storage ring, the process Ξ^{0}n→Ξ^{-}p is studied, where the Ξ^{0} baryon is produced in the process J/ψ→Ξ^{0}Ξ[over ¯]^{0} and the neutron is a component of the ^{9}Be, ^{12}C, and ^{197}Au nuclei in the beam pipe. A clear signal is observed with a statistical significance of 7.1σ. The cross section of the reaction Ξ^{0}+^{9}Be→Ξ^{-}+p+^{8}Be is determined to be σ(Ξ^{0}+^{9}Be→Ξ^{-}+p+^{8}Be)=(22.1±5.3_{stat}±4.5_{sys})  mb at the Ξ^{0} momentum of 0.818  GeV/c, where the first uncertainty is statistical and the second is systematic. No significant H-dibaryon signal is observed in the Ξ^{-}p final state. This is the first study of hyperon-nucleon interactions in electron-positron collisions and opens up a new direction for such research.

A steady-state kinetic interface-sensitive tracer (KIS-SST) method to measure capillary associated interfacial area in a simultaneous co-flow condition

We propose a novel method to measure the specific capillary-associated interfacial area (a_wn) between non-wetting and wetting fluids by applying kinetic interface-sensitive (KIS) tracers under steady-state two-phase co-flow conditions. Seven column experiments were conducted with a column filled with glass beads (diameter d₅₀ = 170 μm), serving as the solid grain network of a porous granular material. The experiments were performed for two flow scenarios, i.e., five experiments were performed for drainage conditions (increasing non-wetting saturation) and two experiments for imbibition conditions (increasing wetting saturation). To obtain different saturation levels in the column and, consequently, different capillarity-induced interfacial areas between the fluids, the experiments were performed for different fractional flow ratios (i.e., the ratio between injection rate of the wetting phase and total injection rate). The concentrations of the KIS tracer reaction by-product were recorded at each saturation level and the corresponding interfacial area was calculated. As a result of the fractional flow condition a wide range of wetting phase saturation is created (0.3 < S_w < 0.8). The measured a_wn increases with decreasing wetting phase saturation for the range 0.55 < S_w < 0.8, and then is followed by a drop in wetting phase saturation of 0.3 < S_w < 0.55. A good fit for our calculated a_wn is obtained (RMSE <0.16) using a polynomial model. Additionally, the results of the proposed method are compared to published experimental data and the method's main advantages and limitations are discussed.

Measurements of Normalized Differential Cross Sections of Inclusive π^{0} and K_{S}^{0} Production in e^{+}e^{-} Annihilation at Energies from 2.2324 to 3.6710 GeV

Based on electron positron collision data collected with the BESIII detector operating at the BEPCII storage rings, the differential cross sections of inclusive π^{0} and K_{S}^{0} production as a function of hadron momentum, normalized by the total cross section of the e^{+}e^{-}→hadrons process, are measured at six center-of-mass energies from 2.2324 to 3.6710 GeV. Our results, which cover a relative hadron energy range from 0.1 to 0.9, significantly deviate from several theoretical calculations based on existing fragmentation functions.

Precision Measurement of the Decay Σ^{+}→pγ in the Process J/ψ→Σ^{+}Σ[over ¯]^{-}

Using (10 087±44)×10^{6} J/ψ events collected with the BESIII detector, the radiative hyperon decay Σ^{+}→pγ is studied at an electron-positron collider experiment for the first time. The absolute branching fraction is measured to be (0.996±0.021_{stat}±0.018_{syst})×10^{-3}, which is lower than its world average value by 4.2 standard deviations. Its decay asymmetry parameter is determined to be -0.652±0.056_{stat}±0.020_{syst}. The branching fraction and decay asymmetry parameter are the most precise to date, and the accuracies are improved by 78% and 34%, respectively.