Patient-specific Motion Research Articles

Soft robotic artificial left ventricle simulator capable of reproducing myocardial biomechanics.

The heart's intricate myocardial architecture has been called the Gordian knot of anatomy, an impossible tangle of intricate muscle fibers. This complexity dictates equally complex cardiac motions that are difficult to mimic in physical systems. If these motions could be generated by a robotic system, then cardiac device testing, cardiovascular disease studies, and surgical procedure training could reduce their reliance on animal models, saving time, costs, and lives. This work introduces a bioinspired soft robotic left ventricle simulator capable of reproducing the minutiae of cardiac motion while providing physiological pressures. This device uses thin-filament artificial muscles to mimic the multilayered myocardial architecture. To demonstrate the device's ability to follow the cardiac motions observed in the literature, we used canine myocardial strain data as input signals that were subsequently applied to each artificial myocardial layer. The device's ability to reproduce physiological volume and pressure under healthy and heart failure conditions, as well as effective simulation of a cardiac support device, were experimentally demonstrated in a left-sided mock circulation loop. This work also has the potential to deliver faithful simulated cardiac motion for preclinical device and surgical procedure testing, with the potential to simulate patient-specific myocardial architecture and motion.

Stability analysis of patient-specific 4DCT- and 4DCBCT-based correspondence models.

Surrogate-based motion compensation in stereotactic body radiation therapy (SBRT) strongly relies on a constant relationship between an external breathing signal and the internal tumor motion over the course of treatment, that is, a stable patient-specific correspondencemodel. This study aims to develop methods for analyzing the stability of correspondence models by integrating planning 4DCT and pretreatment 4D cone-beam computed tomography (4DCBCT) data and assessing the relation to patient-specific clinicalparameters. For correspondence modeling, a regression-based approach is applied, correlating patient-specific internal motion (vector fields computed by deformable image registration) and external breathing signals (recorded by Varian's RPM and RGSC system). To analyze correspondence model stability, two complementary methods are proposed. (1) Target volume-based analysis: 4DCBCT-based correspondence models predict clinical target volumes (GTV and internal target volume [ITV]) within the planning 4DCT, which are evaluated by overlap and distance measures (Dice similarity coefficient [DSC]/average symmetric surface distance [ASSD]). (2) System matrix-based analysis: 4DCBCT-based regression models are compared to 4DCT-based models using mean squared difference (MSD) and principal component analysis of the system matrices. Stability analysis results are correlated with clinical parameters. Both methods are applied to a dataset of 214 pretreatment 4DCBCT scans (Varian TrueBeam) from a cohort of 46 lung tumor patients treated with ITV-based SBRT (planning 4DCTs acquired with Siemens AS Open and SOMATOM go.OPEN Pro CT scanners). Consistent results across the two complementary analysis approaches (Spearman correlation coefficient of between system matrix-based MSD and GTV-based DSC/ASSD) were observed. Analysis showed that stability was not predominant, with 114/214 fraction-wise models not surpassing a threshold of for the GTV, and only 14/46 patients demonstrating a in all fractions. Model stability did not degrade over the course of treatment. The mean GTV-based DSC is (mean ASSD of ) and the respective ITV-based DSC is (mean ASSD of ). The clinical parameters showed a strong correlation between smaller tumor motion ranges and increasedstability. The proposed methods identify patients with unstable correspondence models prior to each treatment fraction, serving as direct indicators for the necessity of replanning and adaptive treatment approaches to account for internal-external motion variations throughout the course oftreatment.

An interdisciplinary treatment to reshape upper anterior displaced teeth by using dynamic occlusal recording with an intraoral scanner.

This article describes an interdisciplinary treatment that helped a patient with displaced upper anterior teeth. A gingivectomy, root canal therapies, digital smile design, digital wax-up, and guided tooth preparations were applied. A patient with pathologically migrated teeth asked for treatment without orthodontic involvement due to a primary failed orthodontic treatment history. A smile photo was taken and superimposed with the dentition in a CAD software to accomplish a digital smile design. The jaw movements were recorded with two different methods, a mechanical articulator and an intraoral scanner with Patient-Specific-Motion function. The occlusal contacts during protrusive and lateral movements were compared and the digital wax-up was designed according to the later occlusal data. An aesthetic crown lengthening and pre-op root canal treatment were carried out in advance accordingly. After guided tooth preparation with a silicone index, the final fixed restorations were manufactured and cemented. A 2-year follow-up showed that our prosthesis functions well. This clinical report revealed that an intraoral scanner with Patient-Specific-Motion function can effectively record individual dynamic occlusal patterns and these data can be integrated into the CAD/CAM process to enhance the fulfillment of clinical requirements. This clinical procedure with a 2-year follow-up demonstrated that a prosthodontic-based interdisciplinary treatment of pathologically migrated teeth using dynamic occlusal recording with an intraoral scanner could achieve satisfactory esthetics in a relatively short treatment period. The Patient Specific Motion module may be used to record a personalized functional movement and the data can be integrated into the design process of the final restorations.

Poster 208: Predicting Patient-Specific Hip Motion During Functional Activities After Periacetabular Osteotomy

Objectives: Developmental dysplasia of the hip (DDH) is a condition in which the hypoplastic acetabulum provides inadequate coverage of the femoral head. This abnormal coverage, while pathological, may be advantageous in some sporting and recreational activities needing significant hip range of motion (ROM), including dance, ballet, and performing arts. Although DDH allows patients to perform activities requiring supraphysiologic hip motion, subluxation of the hip joint generates excessive articular stresses, leading ultimately to pain, osteoarthritis (OA), and total joint arthroplasty. In some DDH patients, this degenerative pathway can be interrupted by performing a periacetabular osteotomy (PAO) to reorient the acetabular fragment, which increases femoral head coverage and reduces intra-articular stresses. This improved coverage is believed to come at the sacrifice of range of motion, particularly noticed while performing high-demand activities requiring supraphysiologic end-range motion. This study addresses the questions: 1. Does a PAO reduce the range of motion of the hip enough to reduce a patient’s ability to perform daily living, recreational and sporting activities? 2. If so, what activities are most impacted by these changes? Methods: We created an inventory of daily living, recreational and sporting activities commonly performed by patients with hip dysplasia. A systematic review of the literature provided the maximum hip range of motion required (required ROM) to perform eight activities of daily living (ADL) and thirteen sporting movements. The concurrent peak value for each anatomic component of hip motion, (i.e flexion/extension; adduction/abduction; internal/external rotation) was classified for every activity in terms of physiologic demand relative to the average ROM of the adult hip. Based on these values, we created a 7-level ordinal scale according to the peak kinematic demand required to perform each activity (Table 1), ranging from low to extreme difficulty. To examine the impact of a PAO on the limits of hip motion, we selected, with IRB approval, a cohort of 24 dysplastic hips (12 left and 12 right). Patient-specific 3D models of each pelvis and hip joint were created from preoperative CT scans, and each activity in our hip inventory was simulated using a custom Matlab hip simulation routine. For each activity, the ROM of the hip was recorded as the point at which either bony impingement was detected, or the specified movement was successfully completed. The percentage of the required motion successfully completed (percent function) was then calculated (100*measured ROM/required ROM) for each movement. This procedure was repeated on the same hips after completion of a computer-simulated PAO according to the technique described by Ganz and colleagues in which the LCEA was increased to 35° and reoriented anteriorly to appropriately restore head coverage. Wilcoxon Signed Rank tests were performed to compare preoperative and postoperative percent function for each movement. Results: When all activities were combined, there was no significant difference (p > 0.05) between preoperative function (99.2 ± 0.6%) and postoperative function (98.0 ± 1.1%). The flexion portions of the right and left grand écart latéral, commonly used ballet movements, were the only activities where percent function was significantly decreased after PAO. Before surgery, the left hips performed 96.8 ± 1.8% of the motion required for the left grand écart latéral whereas after surgery hip function decreased to 82.5 ± 7.6% (p < 0.05). Preoperatively, the right hips performed the right grand écart latéral at 96.7 ± 3.3% compared to 77.4 ± 8.3% postoperatively (p < 0.05). There were no significant differences between preoperative and postoperative function for any other activities within the hip inventory. Conclusions: After testing a broad range of movements found in both daily life and sporting and recreational activities, just two of the twenty-one activities were significantly impaired after a computer-simulated PAO. A lack of significant change in percent function for most movements between preoperative and postoperative hips indicates that even high performing patients who undergo a PAO can expect little to no change in their ability to perform many of the activities they enjoyed before surgery. A significant change, however, was seen in the ballet positions right and left grand écart latéral and patients requiring performance of these activities after PAO should be counseled preoperatively regarding the demands of their sporting and recreational activities and their postoperative expectations. [Figure: see text]

Dosimetric impact of respiratory motion on proximal bronchial tree during lung cancer stereotactic body radiation therapy: A patient-specific phantom study

IntroductionThe respiratory motion of the proximal bronchial tree (PBT) adjacent to the tumor is not taken into consideration in the treatment plan, and there is a risk that the PBT may be irradiated with a higher dose. Therefore, this study aimed to measure the PBT dose directly using patient-specific phantoms of patients treated with SBRT for a centrally located lung tumor and to determine the dose variation with the respiratory motion of the PBT. Materials and methodsThis study included 10 patients who received SBRT for central lung tumors. We utilized patient-specific tumor phantoms fabricated with a three-dimensional printer using the patient's computed tomography (CT) data to analyze the dose variation in PBT owing to respiratory motion. Two SBRT plans were generated: a full-phase plan and a 30−70 gating plan. Dose measurements were performed using glass dosimeters. We measured the dose three times for each plan, with the phantom either in motion or fixed at the 50% phase. ResultsWe analyzed the target volumes and measured doses for each SBRT plan. Compared with the gross tumor volume on the 50% phase CT image, the internal target volume was larger by approximately 13.1% and 45.9% during the 30−70% phase and full phases CT images, respectively. Similarly, the planning target volumes increased by 9.3% and 34.4%, respectively. The average variations in the measured dose of the glass dosimeter owing to patient-specific motion were 3.85% for the full-phase plans and 1.92% for the 30−70 gating plans. ConclusionRespiratory motion can affect the absorbed dose of the PBT, and this effect is particularly pronounced when the PBT is near the tumor. The gating method reduces the degree of dose variation. This study highlights the importance of understanding the respiratory motion of the PBT in SBRT for centrally located lung tumors.

Dynamic CBCT imaging using prior model-free spatiotemporal implicit neural representation (PMF-STINR)

Objective. Dynamic cone-beam computed tomography (CBCT) can capture high-spatial-resolution, time-varying images for motion monitoring, patient setup, and adaptive planning of radiotherapy. However, dynamic CBCT reconstruction is an extremely ill-posed spatiotemporal inverse problem, as each CBCT volume in the dynamic sequence is only captured by one or a few x-ray projections, due to the slow gantry rotation speed and the fast anatomical motion (e.g. breathing). Approach. We developed a machine learning-based technique, prior-model-free spatiotemporal implicit neural representation (PMF-STINR), to reconstruct dynamic CBCTs from sequentially acquired x-ray projections. PMF-STINR employs a joint image reconstruction and registration approach to address the under-sampling challenge, enabling dynamic CBCT reconstruction from singular x-ray projections. Specifically, PMF-STINR uses spatial implicit neural representations to reconstruct a reference CBCT volume, and it applies temporal INR to represent the intra-scan dynamic motion of the reference CBCT to yield dynamic CBCTs. PMF-STINR couples the temporal INR with a learning-based B-spline motion model to capture time-varying deformable motion during the reconstruction. Compared with the previous methods, the spatial INR, the temporal INR, and the B-spline model of PMF-STINR are all learned on the fly during reconstruction in a one-shot fashion, without using any patient-specific prior knowledge or motion sorting/binning. Main results. PMF-STINR was evaluated via digital phantom simulations, physical phantom measurements, and a multi-institutional patient dataset featuring various imaging protocols (half-fan/full-fan, full sampling/sparse sampling, different energy and mAs settings, etc). The results showed that the one-shot learning-based PMF-STINR can accurately and robustly reconstruct dynamic CBCTs and capture highly irregular motion with high temporal (∼ 0.1 s) resolution and sub-millimeter accuracy. Significance. PMF-STINR can reconstruct dynamic CBCTs and solve the intra-scan motion from conventional 3D CBCT scans without using any prior anatomical/motion model or motion sorting/binning. It can be a promising tool for motion management by offering richer motion information than traditional 4D-CBCTs.

Computational Lower Limb Simulator Boundary Conditions to Reproduce Measured TKA Loading in a Cohort of Telemetric Implant Patients.

Recent advancements in computational modeling offer opportunities to refine total knee arthroplasty (TKA) design and treatment strategies. This study developed patient-specific simulator external boundary conditions (EBCs) using a PID-controlled lower limb finite element (FE) model. Calibration of the external actuation required to achieve measured patient-specific joint loading and motion was completed for nine patients with telemetric implants during gait, stair descent, and deep knee bend. The study also compared two EBC scenarios: activity-specific hip AP motion and pelvic rotation (that was averaged across all patients for an activity) and patient-specific hip AP motion and pelvic rotation. Including patient-specific data significantly improved reproduction of joint-level loading, reducing root mean squared error between the target and achieved loading by 28.7% and highlighting the importance of detailed patient data in replicating joint kinematics and kinetics. The principal component analysis (PCA) of the EBCs for the patient dataset showed that one component represented 77.8% of the overall variation, while the first three components represented 97.8%. Given the significant loading variability within the patient cohort, this group of patient-specific models can be run individually to provide insight into expected TKA mechanics variability, and the PCA can be utilized to further create reasonable EBCs that expand the variability evaluated.

A fully digital workflow to design anterior guidance for an implant-supported single crown using a modified patient-specific motion technique

Current methods for designing anterior guidance of anterior fixed prostheses are either complicated or lack accuracy. The article describes a fully digital workflow to design individualized anterior guidance of an implant-supported single crown by using a modified patient-specific motion technique. The technique aims to optimize the digital occlusal design workflow, thereby improving the occlusal fit and long-term stability of anterior fixed prostheses.

Digital Protocol to Record Occlusal Analysis in Prosthodontics: A Pilot Study.

Background: Digital technologies enable the accurate replication of occlusion, which is pivotal for stability in maximum intercuspation and dynamic occlusion. CAD softwares generates standardized occlusal morphologies requiring significant adjustments. The consideration of individual mandibular movements during restoration leads to better functional integration. This pilot study evaluates the efficacy of a novel, fully digital protocol for occlusal analysis recording in prosthodontics. Methods: Patients needing single or multiple metal-free restorations were included. Teeth underwent horizontal finish line preparation, while restorations on implants were either directly screwed or used multi-unit abutments. A digital impression (Trios 3 Intraoral Scanner) captured the mouth's elements. Dynamic occlusion was recorded via Patient Specific Motion (PSM). After the placement and functionalization of temporary restorations, subsequent scans included various elements, and CAD software (Dental system) was used for the restoration design. Restorations were milled in monolithic zirconia, pressed from CAD/CAM-milled wax, and sintered. Results: An evaluation of 52 restorations in 37 patients indicated high accuracy in restorations manufactured via the fully digital workflow. Monolithic zirconia was predominantly used. Subtractive (17.3%) and additive (7.7%) occlusal adjustments were mainly chairside. Conclusion: This study underscores the efficacy of meticulous verification measures and a centric contact system in reducing the need for clinical occlusal refinements in prosthetic restorations.

Is dynamic occlusal design necessary for anterior guidance recovery in the computer-aided design process? An in vitro study

ObjectivesThis study aimed to assess the effectiveness of patient-specific motion in restoring anterior guidance and to investigate the influence of occlusal plane position within a virtual articulator on the design of the anterior guide slope for incisors. MethodsTwenty participants' intraoral scan, occlusal plane position, and jaw motion data were recorded. The maxillary anterior teeth were virtually prepared, and the crowns were designed based on average virtual articulator (AVR), personalized virtual articulator (ART), and patient-specific motion (PSM). The anterior guide slope of maxillary central incisors (S1, S2, Sc, Sp) and the mesio-distal angle (MDA) of the canine of prostheses were compared to that of natural teeth (NAT). One-Way ANOVA was utilized to evaluate the effectiveness of the three methods in restoring the anterior guidance of maxillary anterior teeth. ResultsThe comparison of Sp and MDA showed no significant difference between the PSM and NAT groups (p > 0.05). However, Sp of the ART group was significantly smaller, while MDA was higher than that of the NAT group (p < 0.05). Sp did not differ significantly (p > 0.05) when the angle of the occlusal plane (AOP) was small. As the AOP increased, Sp of the ART and AVR groups were significantly smaller than that of the NAT group (p < 0.05). With a large AOP, Sp of the ART group was notably smaller than that of the NAT group (p < 0.05), while there was no significant difference between the AVR and NAT groups (p > 0.05). ConclusionsOcclusal design based on patient-specific motion proved more effective in restoring natural anterior guidance. The anterior guidance of prostheses designed using a virtual articulator was influenced by occlusal plane position. Clinical SignificanceThe utilization of a jaw motion tracer facilitated the transfer of personalized occlusal plane positions and recorded jaw motion, which can be integrated into the digital prosthetic workflow for virtual occlusion adjustment. Occlusal design based on patient-specific motion more effectively restored lingual guidance of maxillary anterior crowns.

Dynamic Chest Radiograph Simulation Technique with Deep Convolutional Neural Networks: A Proof-of-Concept Study.

In this study, we present an innovative approach that harnesses deep neural networks to simulate respiratory lung motion and extract local functional information from single-phase chest X-rays, thus providing valuable auxiliary data for early diagnosis of lung cancer. A novel radiograph motion simulation (RMS) network was developed by combining a U-Net and a long short-term memory (LSTM) network for image generation and sequential prediction. By utilizing a spatial transformer network to deform input images, our proposed network ensures accurate image generation. We conducted both qualitative and quantitative assessments to evaluate the effectiveness and accuracy of our proposed network. The simulated respiratory motion closely aligns with pulmonary biomechanics and reveals enhanced details of pulmonary diseases. The proposed network demonstrates precise prediction of respiratory motion in the test cases, achieving remarkable average Dice scores exceeding 0.96 across all phases. The maximum variation in lung length prediction was observed during the end-exhale phase, with average deviation of 4.76 mm (±6.64) for the left lung and 4.77 mm (±7.00) for the right lung. This research validates the feasibility of generating patient-specific respiratory motion profiles from single-phase chest radiographs.

Interplay Effect of Splenic Motion for Total Lymphoid Irradiation in Pediatric Proton Therapy.

(1) Background: The most significant cause of an unacceptable deviation from the planned dose during respiratory motion is the interplay effect. We examined the correlation between the magnitude of splenic motion and its impact on plan quality for total lymphoid irradiation (TLI); (2) Methods: Static and 4D CT images from ten patients were used for interplay effect simulations. Patients' original plans were optimized based on the average CT extracted from the 4D CT and planned with two posterior beams using scenario-based optimization (±3 mm of setup and ±3% of range uncertainty) and gradient matching at the level of mid-spleen. Dynamically accumulated 4D doses (interplay effect dose) were calculated based on the time-dependent delivery sequence of radiation fluence across all phases of the 4D CT. Dose volume parameters for each simulated treatment delivery were evaluated for plan quality; (3) Results: Peak-to-peak splenic motion (≤12 mm) was measured from the 4D CT of ten patients. Interplay effect simulations revealed that the ITV coverage of the spleen remained within the protocol tolerance for splenic motion, ≤8 mm. The D100% coverage for ITV spleen decreased from 95.0% (nominal plan) to 89.3% with 10 mm and 87.2% with 12 mm of splenic motion; (4) Conclusions: 4D plan evaluation and robust optimization may overcome problems associated with respiratory motion in proton TLI treatments. Patient-specific respiratory motion evaluations are essential to confirming adequate dosimetric coverage when proton therapy is utilized.

Evaluation of Patient-Specific Quality Control (QC) for Markerless Dynamic Tumor Tracking (MDTT) Deliveries

To report on the patient-specific QC for MDTT for liver and lung tumors treated using the dome-of-liver/diaphragm as a surrogate for the target position and compare to the QC of conventional fiducial-based tracking. Vero4DRT is a linear accelerator that has a gimbal-mounted waveguide and collimation system allowing the radiation beam to perform dynamic tumor tracking based on a patient's respiratory motion. A correlation model is built between an external IR marker placed on the chest and an internal structure, which may be implanted fiducial surrogates or in the case of MDTT, an anatomic landmark. During treatment an orthogonal kV image pair is taken every second. The auto detected internal structure position versus the predicted position is recorded in treatment log files. Five IMRT tracking plans were delivered to two commercial motion phantoms with in-house additions for patient-specific tracking QC. For point dose measurements, a 0.6cc farmer chamber was placed inside the dynamic tracking phantom which contains fiducials for fiducial-based tracking delivery and an anatomic landmark for MDTT. For 2D dose distribution measurements, a radiographic film insert for the Quasar Respiratory Motion Phantom was constructed. The phantom contains both fiducials and has a liver-dome shape at one end to facilitate MDTT. Both motion platforms were interfaced with software that enabled patient-specific respiratory motion traces acquired during a 4DCT scan. Each plan was delivered three times on each phantom: 1) in static mode, 2) during fiducial-based tracking, and 3) during MDTT. Chamber measurements taken during static and tracked delivery were compared to the average chamber dose from the treatment planning software (TPS). Film distributions measured during tracked deliveries were compared to the static film distribution using gamma analysis with FILMQA PRO software. Finally, dynamic tracking treatment statistics, extracted from log files, were compared between fiducial and markerless tracking deliveries. All chamber measurements resulted in dose differences <2.5% compared with the TPS. Dose differences between fiducial-based and MDTT were 0.26% on average (range 0.03-0.62%). For film dose map analysis, gamma pass rates were >95% for all plans for all tracking methods. The average gamma pass rate for 3%/3mm was 99.4% (fiducial-based) vs 98.4% (markerless). For 2%/2mm, the average gamma pass rate was 96.4% (fiducial-based tracking) vs 95.6% (MDTT). Dynamic treatment statistics from logs files reported an average 3D absolute deviation from predicted position of 0.52mm (±0.24) for fiducial tracking and 1.19mm (±0.50) for MDTT. Both fiducial-based and MDTT plans meet passing criteria for patient-specific QC. 3D absolute deviation of detected versus predicted position are larger for MDTT compared to fiducial tracking, however no significant impact on dose delivery was measured.

Extracting lung contour deformation features with deep learning for internal target motion tracking: a preliminary study

Objective. To propose lung contour deformation features (LCDFs) as a surrogate to estimate the thoracic internal target motion, and to report their performance by correlating with the changing body using a cascade ensemble model (CEM). LCDFs, correlated to the respiration driver, are employed without patient-specific motion data sampling and additional training before treatment. Approach. LCDFs are extracted by matching lung contours via an encoder–decoder deep learning model. CEM estimates LCDFs from the currently captured body, and then uses the estimated LCDFs to track internal target motion. The accuracy of the proposed LCDFs and CEM were evaluated using 48 targets’ motion data, and compared with other published methods. Main results. LCDFs estimated the internal targets with a localization error of 2.6 ± 1.0 mm (average ± standard deviation). CEM reached a localization error of 4.7 ± 0.9 mm and a real-time performance of 256.9 ± 6.0 ms. With no internal anatomy knowledge, they achieved a small accuracy difference (of 0.34∼1.10 mm for LCDFs and of 0.43∼1.75 mm for CEM at the 95% confidence level) with a patient-specific lung biomechanical model and the deformable image registration models. Significance. The results demonstrated the effectiveness of LCDFs and CEM on tracking target motion. LCDFs and CEM are non-invasive, and require no patient-specific training before treatment. They show potential for broad applications.

Elevated atrial blood stasis in paroxysmal atrial fibrillation during sinus rhythm: a patient-specific computational fluid dynamics study.

Atrial fibrillation (AF) is associated with an increased risk of stroke, often caused by thrombi that form in the left atrium (LA), and especially in the left atrial appendage (LAA). The underlying mechanism is not fully understood but is thought to be related to stagnant blood flow, which might be present despite sinus rhythm. However, measuring blood flow and stasis in the LAA is challenging due to its small size and low velocities. We aimed to compare the blood flow and stasis in the left atrium of paroxysmal AF patients with controls using computational fluid dynamics (CFD) simulations. The CFD simulations were based on time-resolved computed tomography including the patient-specific cardiac motion. The pipeline allowed for analysis of 21 patients with paroxysmal AF and 8 controls. Stasis was estimated by computing the blood residence time. Residence time was elevated in the AF group (p < 0.001). Linear regression analysis revealed that stasis was strongest associated with LA ejection ratio (p < 0.001, R2 = 0.68) and the ratio of LA volume and left ventricular stroke volume (p < 0.001, R2 = 0.81). Stroke risk due to LA thrombi could already be elevated in AF patients during sinus rhythm. In the future, patient specific CFD simulations may add to the assessment of this risk and support diagnosis and treatment.

Are Abnormal Muscle Biomechanics and Patient-reported Outcomes Associated in Patients With Hip Dysplasia?

Developmental dysplasia of the hip (DDH) is a major risk factor for the early development of hip osteoarthritis. Recent studies have demonstrated how DDH alters hip muscle moment arms and elevates muscle-induced biomechanical variables such as joint reaction forces and acetabular edge loads. Understanding the link between abnormal biomechanics and patient-reported outcome measures (PROMs) is important for evidence-based clinical interventions that improve patient symptoms and functional outcomes. To our knowledge, there are no reports of the relationships between muscle-induced biomechanics and PROMs. (1) Are there associations between PROMs and muscle-induced hip biomechanics during gait for patients with DDH and controls? (2) Are there associations among PROMs and separately among biomechanical variables? Participants in this prospective cross-sectional comparative study included 20 female patients with DDH who had no prior surgery or osteoarthritis and 15 female individuals with no evidence of hip pathology (controls) (age: median 23 years [range 16 to 39 years]; BMI: median 22 kg/m 2 [range 17 to 27 kg/m 2 ]). Muscle-induced biomechanical variables for this cohort were reported and had been calculated from patient-specific musculoskeletal models, motion data, and MRI. Biomechanical variables included joint reaction forces, acetabular edge loads, hip center lateralization, and gluteus medius muscle moment arm lengths. PROMs included the Hip Disability and Osteoarthritis Outcome Score (HOOS), the WOMAC, International Hip Outcome Tool-12, National Institutes of Health Patient-Reported Outcome Measure Information System (PROMIS) Pain Interference and Physical Function subscales, and University of California Los Angeles activity scale. Associations between PROMs and biomechanical variables were tested using Spearman rank-order correlations and corrected for multiple comparisons using the Benjamini-Yekutieli method. For this study, associations between variables were considered to exist when correlations were statistically significant (p < 0.05) and were either strong (ρ ≥ 0.60) or moderate (ρ = 0.40 to 0.59). Acetabular edge load impulses (the cumulative acetabular edge load across the gait cycle), medially directed joint reaction forces, and hip center lateralization most commonly demonstrated moderate or strong associations with PROMs. The strongest associations were a negative correlation between acetabular edge load impulse on the superior acetabulum and the HOOS function in daily living subscale (ρ = -0.63; p = 0.001), followed by a negative correlation between hip center lateralization and the HOOS pain subscale (ρ = -0.6; p = 0.003), and a positive correlation between hip center lateralization and the PROMIS pain subscale (ρ = 0.62; p = 0.002). The University of California Los Angeles activity scale was the only PROM that did not demonstrate associations with any biomechanical variable. All PROMs, aside from the University of California Los Angeles activity scale, were associated with one another. Although most of the biomechanical variables were associated with one another, these relationships were not as consistent as those among PROMs. The associations with PROMs detected in the current study suggest that muscle-induced biomechanics may have wide-reaching effects not only on loads within the hip, but also on patients' perceptions of their health and function. As the treatment of DDH evolves, patient-specific joint preservation strategies may benefit from targeting the underlying causes of biomechanical outcomes associated with PROMs. Level III, prognostic study.

A motion model-guided 4D dose reconstruction for pencil beam scanned proton therapy

Objective. 4D dose reconstruction in proton therapy with pencil beam scanning (PBS) typically relies on a single pre-treatment 4DCT (p4DCT). However, breathing motion during the fractionated treatment can vary considerably in both amplitude and frequency. We present a novel 4D dose reconstruction method combining delivery log files with patient-specific motion models, to account for the dosimetric effect of intra- and inter-fractional breathing variability. Approach. Correlation between an external breathing surrogate and anatomical deformations of the p4DCT is established using principal component analysis. Using motion trajectories of a surface marker acquired during the dose delivery by an optical tracking system, deformable motion fields are retrospectively reconstructed and used to generate time-resolved synthetic 4DCTs (‘5DCTs’) by warping a reference CT. For three abdominal/thoracic patients, treated with respiratory gating and rescanning, example fraction doses were reconstructed using the resulting 5DCTs and delivery log files. The motion model was validated beforehand using leave-one-out cross-validation (LOOCV) with subsequent 4D dose evaluations. Moreover, besides fractional motion, fractional anatomical changes were incorporated as proof of concept. Main results. For motion model validation, the comparison of 4D dose distributions for the original 4DCT and predicted LOOCV resulted in 3%/3 mm gamma pass rates above 96.2%. Prospective gating simulations on the p4DCT can overestimate the target dose coverage V95% by up to 2.1% compared to 4D dose reconstruction based on observed surrogate trajectories. Nevertheless, for the studied clinical cases treated with respiratory-gating and rescanning, an acceptable target coverage was maintained with V95% remaining above 98.8% for all studied fractions. For these gated treatments, larger dosimetric differences occurred due to CT changes than due to breathing variations. Significance. To gain a better estimate of the delivered dose, a retrospective 4D dose reconstruction workflow based on motion data acquired during PBS proton treatments was implemented and validated, thus considering both intra- and inter-fractional motion and anatomy changes.

A CARDIAC MOTION MODEL TO EVALUATE INTRA-FRACTION DOSIMETRIC VARIATIONS IN RADIOTHERAPY TREATMENTS

Purpose: A cardiac cycle model was implemented to simulate cardiac motion during radiotherapy to evaluate the intra-fraction dosimetric impact on cardiac sub-structures comparing different planning techniques. Methods: Cardiac sub-structures were automatically contoured in 10 CTs acquired in deep inspiration breath hold (DIBH) by using a recently developed hierarchical-clustering atlas-based algorithm. A deformable image registration algorithm was used to simulate the cardiac motion cycle based on volume variations available in the literature. Two synthetic CTs were created and contoured simulating contraction and expansion during the cardiac cycle. Ninety radiotherapy plans were calculated using three radiotherapy paradigms: tangential fields planned as Linac-3D-CRT with a steep linear dose gradient toward the heart-modulated therapy with an intermediately steep concave gradient of intermediate-to-high doses toward the heart, planned as Linac-VMAT; modulated therapy with a steep concave gradient of intermediate-to-high doses toward the heart, planned as helical tomotherapy. Python scripts were developed for autocontouring, automatic creation of synthetic CTs and data extraction. Results: Comparison between paradigms shows that different constraints (maximal gradient toward heart/lung versus maximal sparing of contralateral breast/axilla) do not necessarily result in preferred or reduced heart sparing, but this depends more on individual anatomy. A planning paradigm with an intermediate-steepness dose gradient showed the best robustness against intra-fraction organ motion. Conclusions: Patient-specific organ motion models may reduce differences between planned and delivered RT and may thus help to refine dose–volume–toxicity models for cardiac sub-structures and, as a consequence, clinical constraints. Automatized plan recalculation on synthetic image sets might be used for robustness optimization and evaluation.

Evaluating motion of pancreatic tumors and anatomical surrogates using cine MRI in 0.35T MRgRT under free breathing conditions.

Treatment tolerability is a significant limitation to pancreatic cancer treatment with radiotherapy due to proximity to highly radiosensitive organs and respiratory motion necessitating expanded target margins. Further, pancreatic tumors are difficult to visualize on conventional radiotherapy systems. Surrogates are often used to locate the tumor but are often inconsistent and do not provide strong positional relations throughout the respiratory cycle. This work utilizes a retrospective dataset of 45 pancreatic cancer patients treated on an MR-Linac system with cine MRI acquired for real-time target tracking. We investigated intra-fraction motion of tumors and two abdominal surrogates, leading to prediction models between the tumor and surrogate. Patient specific motion evaluation and prediction models were generated from 225 cine MRI series acquired during treatment. Tumor contours were used to evaluate the pancreatic tumor motion. Linear regression and principal component analysis (PCA) based models were used to predict tumor position from the anterior-posterior (AP) motion of the abdominal surface, the superior-inferior (SI) motion of the diaphragm, or a combination. Models were evaluated using mean squared error (MSE) and mean absolute error (MAE). Contour analysis showed the average pancreatic tumor motion range was 7.4±2.7mm and 14.9±5.8mm in the AP and SI directions, respectively. The PCA model had MSE of 1.4mm2 and 0.6mm2 , for the SI and AP directions, respectively, with both surrogates as inputs for the models. When only the abdomen surrogate was used, MSE was 1.3mm2 and 0.4mm2 in the SI and AP directions, while it was 0.4mm2 and 1.3mm2 when only the diaphragm surrogate was used. We evaluated intra-fraction pancreatic tumor motion and demonstrated prediction models between the tumor and surrogate. The models calculated the pancreatic tumor position from diaphragm, abdominal, or both contours within standard pancreatic cancer target margin, and the process could be applied to other disease sites in the abdominothoracic cavity.

4DCT is long overdue for improvement.

4DCT is long overdue for improvement.