Browsing by Author "Castillo, Edward"
Now showing 1 - 10 of 10
Results Per Page
Sort Options
Item Association of anticoagulation dose and survival in hospitalized COVID‐19 patients: A retrospective propensity score‐weighted analysis(Wiley, 2020) Ionescu, Filip; Jaiyesimi, Ishmael; Petrescu, Ioana; Lawler, Patrick R.; Castillo, Edward; Munoz‐Maldonado, Yolanda; Imam, Zaid; Narasimhan, Mangala; Abbas, Amr E.; Konde, Anish; Nair, Girish B.Background: Hypercoagulability may contribute to COVID‐19 pathogenicity. The role of anticoagulation (AC) at therapeutic (tAC) or prophylactic doses (pAC) is unclear. Objectives: We evaluated the impact on survival of different AC doses in COVID‐19 patients. Methods: Retrospective, multi‐center cohort study of consecutive COVID‐19 patients hospitalized between March 13 and May 5, 2020. Results: A total of 3480 patients were included (mean age, 64.5 years [17.0]; 51.5% female; 52.1% black and 40.6% white). 18.5% (n = 642) required intensive care unit (ICU) stay. 60.9% received pAC (n = 2121), 28.7% received ≥3 days of tAC (n = 998), and 10.4% (n = 361) received no AC. Propensity score (PS) weighted Kaplan‐Meier plot demonstrated different 25‐day survival probability in the tAC and pAC groups (57.5% vs 50.7%). In a PS–weighted multivariate proportional hazards model, AC was associated with reduced risk of death at prophylactic (hazard ratio [HR] 0.35 [95% confidence interval {CI} 0.22‐0.54]) and therapeutic doses (HR 0.14 [95% CI 0.05‐0.23]) compared to no AC. Major bleeding occurred more frequently in tAC patients (81 [8.1%]) compared to no AC (20 [5.5%]) or pAC (46 [2.2%]) subjects. Conclusions: Higher doses of AC were associated with lower mortality in hospitalized COVID‐19 patients. Prospective evaluation of efficacy and risk of AC in COVID‐19 is warranted.Item Extubation Failure in Critically Ill COVID-19 Patients: Risk Factors and Impact on In-Hospital Mortality(Sage, 2021) Ionescu, Filip; Zimmer, Markie S.; Petrescu, Ioana; Castillo, Edward; Bozyk, Paul; Abbas, Amr; Abplanalp, Lauren; Dogra, Sanjay; Nair, Girish B.Purpose:We sought to identify clinical factors that predict extubation failure (reintubation) and its prognostic implications in critically ill COVID-19 patients.Materials and Methods:Retrospective, multi-center cohort study of hospitalized COVID-19 patients. Multivariate competing risk models were employed to explore the rate of reintubation and its determining factors.Results:Two hundred eighty-one extubated patients were included (mean age, 61.0 years [±13.9]; 54.8% male). Reintubation occurred in 93 (33.1%). In multivariate analysis accounting for death, reintubation risk increased with age (hazard ratio [HR] 1.04 per 1-year increase, 95% confidence interval [CI] 1.02 -1.06), vasopressors (HR 1.84, 95% CI 1.04-3.60), renal replacement (HR 2.01, 95% CI 1.22-3.29), maximum PEEP (HR 1.07 per 1-unit increase, 95% CI 1.02 -1.12), paralytics (HR 1.48, 95% CI 1.08-2.25) and requiring more than nasal cannula immediately post-extubation (HR 2.19, 95% CI 1.37-3.50). Reintubation was associated with higher mortality (36.6% vs 2.1%; P < 0.0001) and risk of inpatient death after adjusting for multiple factors (HR 23.2, 95% CI 6.45-83.33). Prone ventilation, corticosteroids, anticoagulation, remdesivir and tocilizumab did not impact the risk of reintubation or death.Conclusions:Up to 1 in 3 critically ill COVID-19 patients required reintubation. Older age, paralytics, high PEEP, need for greater respiratory support following extubation and non-pulmonary organ failure predicted reintubation. Extubation failure strongly predicted adverse outcomes.Item Functional avoidance-based intensity modulated proton therapy with 4DCT derived ventilation imaging for lung cancer(Wiley, 2021) Dougherty, Jingjing M.; Castillo, Edward; Castillo, Richard; Faught, Austin M.; Pepin, Mark; Park, Sean S.; Beltran, Chris J.; Guerrero, Thomas; Grills, Inga; Vinogradskiy, YevgeniyThe primary objective is to evaluate the potential dosimetric gains of performing functional avoidance-based proton treatment planning using 4DCT derived ventilation imaging. 4DCT data of 31 patients from a prospective functional avoidance clinical trial were evaluated with intensity modulated proton therapy (IMPT) plans and compared with clinical volumetric modulated arc therapy (VMAT) plans. Dosimetric parameters were compared between standard and functional plans with IMPT and VMAT with one-way analysis of variance and post hoc paired student t-test. Normal Tissue Complication Probability (NTCP) models were employed to estimate the risk of two toxicity endpoints for healthy lung tissues. Dose degradation due to proton motion interplay effect was evaluated. Functional IMPT plans led to significant dose reduction to functional lung structures when compared with functional VMAT without significant dose increase to Organ at Risk (OAR) structures. When interplay effect is considered, no significant dose degradation was observed for the OARs or the clinical target volume (CTV) volumes for functional IMPT. Using fV20 as the dose metric and Grade 2+ pneumonitis as toxicity endpoint, there is a mean 5.7% reduction in Grade 2+ RP with the functional IMPT and as high as 26% in reduction for individual patient when compared to the standard IMPT planning. Functional IMPT was able to spare healthy lung tissue to avoid excess dose to normal structures while maintaining satisfying target coverage. NTCP calculation also shows that the risk of pulmonary complications can be further reduced with functional based IMPT.Item Medical image registration: A review of existing methods and preliminary numerical results(2005) Castillo, Edward; Zhang, YinRegistration of medical images has become an important area of research. In particular, registration of computed tomography (CT) lung images is of great interest to radiation oncologists planning radiation treatment for patients with lung cancer. A review of existing image registration methods, as well as preliminary numerical results, indicate that methods based on a constant pixel intensity assumption, such as traditional optical flow methods, cannot be expected to produce accurate registration of lung CT images. Nonlinear methods allowing variations in pixel intensities, though more costly than linear methods, promise to be more accurate for this application.Item Optical Flow Methods for the Registration of Compressible Flow Images and Images Containing Large Voxel Displacements or Artifacts(2008-10) Castillo, EdwardThree optical flow image registration (IR) methods referred to as Combined Compressible Local Global (CCLG) optical flow, Large Displacements Optical Flow (LDOF), and Large Displacement Compressible Optical Flow (LDCOF) are introduced. The three novel methods are designed to account for difficulties raised by 4D throacic Computed Tomography (CT) image registration problem, which currently cannot be effectively addressed by existing methods. The 4D CT image registration problem is more challenging than typical IR problems for three key reasons. First, voxel intensities for CT images are proportional to the density of the material imaged. Given that the density of lung tissue changes with respiration, the constant voxel intensity assumption employed by most IR methods is invalid for thoracic CT images. Second, due to the image acquisition procedure, 4D CT image sets are known to suffer from image noise, blurring, and artifacts. Finally, the large size of the image sets requires a computationally efficient and parallelizable algorithm. The CCLG method models compressible image flow with the mass conservation equation coupled with a local-global strategy that alleviates the effects of image noise, and incorporates local image information into the voxel motion model. After a finite element discretization, the resulting large scale linear system is solved using a parallelizable, multi-grid preconditioned conjugate gradient algorithm. The LDOF and LDCOF methods are designed for image sets containing large voxel displacements or erroneous image artifacts. Both methods incorporate unknown image information into the IR problem formulation, which results in a nonlinear least squares problem for both the pixel displacement components and the unknown image values. An alternating linear least squares algorithm is introduced for solving the LDOF and LDCOF nonlinear least squares problems efficiently. After Chapter 1 introduces the basics of IR, the main body of the thesis is divided into two parts. Part 1 is a review of existing IR methodologies. Part 2 derives the three aforementioned new approaches, and presents the results of numerical experiments testing each of the three methods. The computational experiments are carried out on both synthetic and genuine image data. Finally, the thesis concludes in Chapter 8 with a discussion on areas of future research.Item Optical flow methods for the registration of compressible flow images and images containing large voxel displacements or artifacts(2007) Castillo, Edward; Zhang, YinThree optical flow image registration (IR) methods referred to as Combined Compressible Local Global (CCLG) optical flow, Large Displacements Optical Flow (LDOF), and Large Displacement Compressible Optical Flow (LDCOF) are introduced. The three novel methods are designed to account for difficulties raised by 4D throacic Computed Tomography (CT) image registration problems, which currently cannot be effectively addressed by existing methods. The 4D CT image registration problem is more challenging than typical IR problems for three key reasons. First, voxel intensities for CT images are proportional to the density of the material imaged. Given that the density of lung tissue changes with respiration, the constant voxel intensity assumption employed by most IR methods is invalid for thoracic CT images. Second, due to the image acquisition procedure, 4D CT image sets are known to suffer from image noise, blurring, and artifacts. Finally, the large size of the image sets requires a computationally efficient and parallelizable algorithm. The CCLG method models compressible image flow with the mass conservation equation coupled with a local-global strategy that alleviates the effects of image noise, and incorporates local image information into the voxel motion model. After a finite element discretization, the resulting large scale linear system is solved using a parallelizable, multi-grid preconditioned conjugate gradient algorithm. The LDOF and LDCOF methods are designed for image sets containing large voxel displacements or erroneous image artifacts. Both methods incorporate unknown image information into the IR problem formulation, which results in a nonlinear least squares problem for both the pixel displacement components and the unknown image values. An alternating linear least squares algorithm is introduced for solving the LDOF and LDCOF nonlinear least squares problems efficiently. After Chapter 1 introduces the basics of IR, the main body of the thesis is divided into two parts. Part 1 is a review of existing IR methodologies. Part 2 derives the three aforementioned new approaches and presents testing results for the three methods, respectively. The computational experiments are carried out on both synthetic and genuine image data. Finally, the thesis concludes in Chapter 8 with a discussion of possible areas of future research.Item Quantifying pulmonary perfusion from noncontrast computed tomography(Wiley, 2021) Castillo, Edward; Nair, Girish; Turner‐Lawrence, Danielle; Myziuk, Nicholas; Emerson, Scott; Al‐Katib, Sayf; Westergaard, Sarah; Castillo, Richard; Vinogradskiy, Yevgeniy; Quinn, Thomas; Guerrero, Thomas; Stevens, CraigPurpose: Computed tomography (CT)-derived ventilation methods compute respiratory induced volume changes as a surrogate for pulmonary ventilation. Currently, there are no known methods to derive perfusion information from noncontrast CT. We introduce a novel CT-Perfusion (CT-P) method for computing the magnitude mass changes apparent on dynamic noncontrast CT as a surrogate for pulmonary perfusion. Methods: CT-Perfusion is based on a mass conservation model which describes the unknown mass change as a linear combination of spatially corresponding inhale and exhale HU estimated voxel densities. CT-P requires a deformable image registration (DIR) between the inhale/exhale lung CT pair, a preprocessing lung volume segmentation, and an estimate for the Jacobian of the DIR transformation. Given this information, the CT-P image, which provides the magnitude mass change for each voxel within the lung volume, is formulated as the solution to a constrained linear least squares problem defined by a series of subregional mean magnitude mass change measurements. Similar to previous robust CT-ventilation methods, the amount of uncertainty in a subregional sample mean measurement is related to measurement resolution and can be characterized with respect to a tolerance parameter. Spatial Spearman correlation between single photon emission CT perfusion (SPECT-P) and the proposed CT-P method was assessed in two patient cohorts via a parameter sweep of . The first cohort was comprised of 15 patients diagnosed with pulmonary embolism (PE) who had SPECT-P and 4DCT imaging acquired within 24 h of PE diagnosis. The second cohort was comprised of 15 nonsmall cell lung cancer patients who had SPECT-P and 4DCT images acquired prior to radiotherapy. For each test case, CT-P images were computed for 30 different uncertainty parameter values, uniformly sampled from the range [0.01, 0.125], and the Spearman correlation between the SPECT-P and the resulting CT-P images were computed. Results: The median correlations between CT-P and SPECT-P taken over all 30 test cases ranged between 0.49 and 0.57 across the parameter sweep. For the optimal tolerance τ = 0.0385, the CT-P and SPECT-P correlations across all 30 test cases ranged between 0.02 and 0.82. A one-sample sign test was applied separately to the PE and lung cancer cohorts. A low Spearmen correlation of 15% was set as the null median value and two-sided alternative was tested. The PE patients showed a median correlation of 0.57 (IQR = 0.305). One-sample sign test was statistically significant with 96.5 % confidence interval: 0.20–0.63, P < 0.00001. Lung cancer patients had a median correlation of 0.57(IQR = 0.230). Again, a one-sample sign test for median was statistically significant with 96.5 percent confidence interval: 0.45–0.71, P < 0.00001. Conclusion: CT-Perfusion is the first mechanistic model designed to quantify magnitude blood mass changes on noncontrast dynamic CT as a surrogate for pulmonary perfusion. While the reported correlations with SPECT-P are promising, further investigation is required to determine the optimal CT acquisition protocol and numerical method implementation for CT-P imaging.Item Radial MILO: A 4D Image Registration Algorithm Based on Filtering Block Match Data via l1-minimization(2015-04-21) Vargas, Arturo; Zhang, Yin; Castillo, Edward; Tapia, Richard; Warburton, TimMinimal l1 Perturbation to Block Match Data (MILO) is a spatially accurate image registration algorithm developed for thoracic CT inhale/exhale images. The MILO algorithm consists of three components: (1) creating an initial estimate for voxel displacement via a Mutual Minimizing Block Matching Algorithm (MMBM), (2) a filtering step based on l1 minimization and a uniform B-spline parameterization, and (3) recovering a full displacement field based on the filtered estimates. This thesis presents a variation of MILO for 4DCT images. In practice, the use of uniform B-splines has led to rank deficient linear systems due to the spline's inability to conform to non-structured MMBM estimates. In order to adaptively conform to the data an octree is paired with radial functions. The l1 minimization problem had previously been addressed by employing QR factorization, which required substantial storage. As an alternative a block coordinate descent algorithm is employed, relieving the need for QR factorization. Furthermore, by modeling voxel trajectories as quadratic functions in time, the proposed method is able to register multiple images.Item Robust CT ventilation from the integral formulation of the Jacobian(Wiley, 2019) Castillo, Edward; Castillo, Richard; Vinogradskiy, Yevgeniy; Dougherty, Michele; Solis, David; Myziuk, Nicholas; Thompson, Andrew; Guerra, Rudy; Nair, Girish; Guerrero, ThomasComputed tomography (CT) derived ventilation algorithms estimate the apparent voxel volume changes within an inhale/exhale CT image pair. Transformation-based methods compute these estimates solely from the spatial transformation acquired by applying a deformable image registration (DIR) algorithm to the image pair. However, approaches based on finite difference approximations of the transformation's Jacobian have been shown to be numerically unstable. As a result, transformation-based CT ventilation is poorly reproducible with respect to both DIR algorithm and CT acquisition method. PURPOSE: We introduce a novel Integrated Jacobian Formulation (IJF) method for estimating voxel volume changes under a DIR-recovered spatial transformation. The method is based on computing volume estimates of DIR mapped subregions using the hit-or-miss sampling algorithm for integral approximation. The novel approach allows for regional volume change estimates that (a) respect the resolution of the digital grid and (b) are based on approximations with quantitatively characterized and controllable levels of uncertainty. As such, the IJF method is designed to be robust to variations in DIR solutions and thus overall more reproducible. METHODS: Numerically, Jacobian estimates are recovered by solving a simple constrained linear least squares problem that guarantees the recovered global volume change is equal to the global volume change obtained from the inhale and exhale lung segmentation masks. Reproducibility of the IJF method with respect to DIR solution was assessed using the expert-determined landmark point pairs and inhale/exhale phases from 10 four-dimensional computed tomographies (4DCTs) available on www.dir-lab.com. Reproducibility with respect to CT acquisition was assessed on the 4DCT and 4D cone beam CT (4DCBCT) images acquired for five lung cancer patients prior to radiotherapy. RESULTS: The ten Dir-Lab 4DCT cases were registered twice with the same DIR algorithm, but with different smoothing parameter. Finite difference Jacobian (FDJ) and IFJ images were computed for both solutions. The average spatial errors (300 landmarks per case) for the two DIR solution methods were 0.98 (1.10) and 1.02 (1.11). The average Pearson correlation between the FDJ images computed from the two DIR solutions was 0.83 (0.03), while for the IJF images it was 1.00 (0.00). For intermodality assessment, the IJF and FDJ images were computed from the 4DCT and 4DCBCT of five patients. The average Pearson correlation of the spatially aligned FDJ images was 0.27 (0.11), while it was 0.77 (0.13) for the IFJ method. CONCLUSION: The mathematical theory underpinning the IJF method allows for the generation of ventilation images that are (a) computed with respect to DIR spatial accuracy on the digital voxel grid and (b) based on DIR-measured subregional volume change estimates acquired with quantifiable and controllable levels of uncertainty. Analyses of the experiments are consistent with the mathematical theory and indicate that IJF ventilation imaging has a higher reproducibility with respect to both DIR algorithm and CT acquisition method, in comparison to the standard finite difference approach.Item Severity of radiation pneumonitis, from clinical, dosimetric and biological features: a pilot study(Springer Nature, 2020) Aso, Samantha; Navarro‑Martin, Arturo; Castillo, Richard; Padrones, Susana; Castillo, Edward; Montes, Ana; Martínez, José Ignacio; Cubero, Noelia; López, Rosa; Rodríguez, Laura; Palmero, Ramon; Manresa, Federico; Guerrero, Thomas; Molina, MaríaBackground and objective: Radiation pneumonitis (RP) could be a lethal complication of lung cancer treatment. No reliable predictors of RP severity have been recognized. This prospective pilot study was performed to identify early predictors of high grade lung toxicity and to evaluate clinical, biological or dosimetric features associated with different grades of toxicity. Method: Sixteen patients with non-small cell lung cancer with indication of concurrent chemoradiotherapy using 60 Gy/2 Gy/fraction starting at cycle one of platinum based chemotherapy were included. Bronchoalveolar lavage (BAL), pulmonary function testing (PFT), and 18F-2-fluoro-2-deoxy-D-glucose positron-emission tomography was performed before radiotherapy (RT), after three weeks of treatment, and two months post-RT. For analysis, patients were grouped by grade (low [G1-G2] vs. high [G3-G5]). The two groups were compared to identify predictors of RP. Protein expression BAL and lung tissue metabolism was evaluated in two patients (RP-G1 vs. RP-G3). Categorical variables such as comorbidities, stages and locations were summarized as percentages. Radiation doses, pulmonary function values and time to RP were summarized by medians with ranges or as means with standard deviation. Longitudinal analysis PFT was performed by a T-test. Results: All 16 patients developed RP, as follows: G1 (5 pts; 31.3%); G2 (5 pts; 31.3%); G3 (5 pts; 31.3%); and G5 (1 pts; 6.1%). Patients with high grade RP presented significant decrease (p = 0.02) in diffusing lung capacity for carbon monoxide (DLCO) after three weeks of RT. No correlation between dosimetric values and RP grades was observed. BAL analysis of the selected patients showed that CXCL-1, CD154, IL-1ra, IL-23, MIF, PAI-1 and IFN-γ were overexpressed in the lungs of the RP-G3 patient, even before treatment. The pre-RT SUVmax value in the RP-G3 patient was non-significantly higher than in the patient with RP-G1. Conclusions: RT induces some degree of RP. Our data suggest that decrease in DLCO% is the most sensitive parameter for the early detection of RP. Moreover, we detect biological differences between the two grades of pneumonitis, highlighting the potential value of some cytokines as a prognostic marker for developing high grade lung toxicity. Further multicenter studies with larger sample size are essential to validate these findings.