Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

Use of EORTC Target Definition Guidelines for Dose-Intensified Salvage Radiation Therapy for Recurrent Prostate Cancer: Results of the Quality Assurance Program of the Randomized Trial SAKK 09/10

PURPOSE Different international target volume delineation guidelines exist and different treatment techniques are available for salvage radiation therapy (RT) for recurrent prostate cancer, but less is known regarding their respective applicability in clinical practice. METHODS AND MATERIALS A randomized phase III trial testing 64 Gy vs 70 Gy salvage RT was accompanied by an intense quality assurance program including a site-specific and study-specific questionnaire and a dummy run (DR). Target volume delineation was performed according to the European Organisation for the Research and Treatment of Cancer guidelines, and a DR-based treatment plan was established for 70 Gy. Major and minor protocol deviations were noted, interobserver agreement of delineated target contours was assessed, and dose-volume histogram (DVH) parameters of different treatment techniques were compared. RESULTS Thirty European centers participated, 43% of which were using 3-dimensional conformal RT (3D-CRT), with the remaining centers using intensity modulated RT (IMRT) or volumetric modulated arc technique (VMAT). The first submitted version of the DR contained major deviations in 21 of 30 (70%) centers, mostly caused by inappropriately defined or lack of prostate bed (PB). All but 5 centers completed the DR successfully with their second submitted version. The interobserver agreement of the PB was moderate and was improved by the DR review, as indicated by an increased κ value (0.59 vs 0.55), mean sensitivity (0.64 vs 0.58), volume of total agreement (3.9 vs 3.3 cm(3)), and decrease in the union volume (79.3 vs 84.2 cm(3)). Rectal and bladder wall DVH parameters of IMRT and VMAT vs 3D-CRT plans were not significantly different. CONCLUSIONS The interobserver agreement of PB delineation was moderate but was improved by the DR. Major deviations could be identified for the majority of centers. The DR has improved the acquaintance of the participating centers with the trial protocol.

Allopurinol hypersensitivity is primarily mediated by dose-dependent oxypurinol-specific T cell response

BACKGROUND Allopurinol is a main cause of severe cutaneous adverse reactions (SCAR). How allopurinol induces hypersensitivity remains unknown. Pre-disposing factors are the presence of the HLA-B*58:01 allele, renal failure and possibly the dose taken. OBJECTIVE Using an in vitro model, we sought to decipher the relationship among allopurinol metabolism, HLA-B*58:01 phenotype and drug concentrations in stimulating drug-specific T cells. METHODS Lymphocyte transformation test (LTT) results of patients who had developed allopurinol hypersensitivity were analysed. We generated allopurinol or oxypurinol-specific T cell lines (ALP/OXP-TCLs) from allopurinol naïve HLA-B*58:01(+) and HLA-B*58:01(-) individuals using various drug concentrations. Their reactivity patterns were analysed by flow cytometry and (51) Cr release assay. RESULTS Allopurinol allergic patients are primarily sensitized to oxypurinol in a dose-dependent manner. TCL induction data show that both the presence of HLA-B*58:01 allele and high concentration of drug are important for the generation of drug-specific T cells. The predominance of oxypurinol-specific lymphocyte response in allopurinol allergic patients can be explained by the rapid conversion of allopurinol to oxypurinol in vivo rather than to its intrinsic immunogenicity. OXP-TCLs do not recognize allopurinol and vice versa. Finally, functional avidity of ALP/OXP-TCL is dependent on both the induction dose and HLA-B*58:01 status. CONCLUSIONS AND CLINICAL RELEVANCE This study establishes the important synergistic role of drug concentration and HLA-B*58:01 allele in the allopurinol or oxypurinol-specific T cell responses. Despite the prevailing dogma that Type B adverse drug reactions are dose independent, allopurinol hypersensitivity is primarily driven by oxypurinol-specific T cell response in a dose-dependent manner, particular in the presence of HLA-B*58:01 allele.

Quantitative 3D strain analysis in analogue experiments: Integration of X-ray computed tomography and digital volume correlation techniques

The combination of scaled analogue experiments, material mechanics, X-ray computed tomography (XRCT) and Digital Volume Correlation techniques (DVC) is a powerful new tool not only to examine the 3 dimensional structure and kinematic evolution of complex deformation structures in scaled analogue experiments, but also to fully quantify their spatial strain distribution and complete strain history. Digital image correlation (DIC) is an important advance in quantitative physical modelling and helps to understand non-linear deformation processes. Optical non-intrusive (DIC) techniques enable the quantification of localised and distributed deformation in analogue experiments based either on images taken through transparent sidewalls (2D DIC) or on surface views (3D DIC). X-ray computed tomography (XRCT) analysis permits the non-destructive visualisation of the internal structure and kinematic evolution of scaled analogue experiments simulating tectonic evolution of complex geological structures. The combination of XRCT sectional image data of analogue experiments with 2D DIC only allows quantification of 2D displacement and strain components in section direction. This completely omits the potential of CT experiments for full 3D strain analysis of complex, non-cylindrical deformation structures. In this study, we apply digital volume correlation (DVC) techniques on XRCT scan data of “solid” analogue experiments to fully quantify the internal displacement and strain in 3 dimensions over time. Our first results indicate that the application of DVC techniques on XRCT volume data can successfully be used to quantify the 3D spatial and temporal strain patterns inside analogue experiments. We demonstrate the potential of combining DVC techniques and XRCT volume imaging for 3D strain analysis of a contractional experiment simulating the development of a non-cylindrical pop-up structure. Furthermore, we discuss various options for optimisation of granular materials, pattern generation, and data acquisition for increased resolution and accuracy of the strain results. Three-dimensional strain analysis of analogue models is of particular interest for geological and seismic interpretations of complex, non-cylindrical geological structures. The volume strain data enable the analysis of the large-scale and small-scale strain history of geological structures.

BAYESIAN PHASE I DOSE FINDING IN CANCER TRIALS

This dissertation explores phase I dose-finding designs in cancer trials from three perspectives: the alternative Bayesian dose-escalation rules, a design based on a time-to-dose-limiting toxicity (DLT) model, and a design based on a discrete-time multi-state (DTMS) model. We list alternative Bayesian dose-escalation rules and perform a simulation study for the intra-rule and inter-rule comparisons based on two statistical models to identify the most appropriate rule under certain scenarios. We provide evidence that all the Bayesian rules outperform the traditional ``3+3'' design in the allocation of patients and selection of the maximum tolerated dose. The design based on a time-to-DLT model uses patients' DLT information over multiple treatment cycles in estimating the probability of DLT at the end of treatment cycle 1. Dose-escalation decisions are made whenever a cycle-1 DLT occurs, or two months after the previous check point. Compared to the design based on a logistic regression model, the new design shows more safety benefits for trials in which more late-onset toxicities are expected. As a trade-off, the new design requires more patients on average. The design based on a discrete-time multi-state (DTMS) model has three important attributes: (1) Toxicities are categorized over a distribution of severity levels, (2) Early toxicity may inform dose escalation, and (3) No suspension is required between accrual cohorts. The proposed model accounts for the difference in the importance of the toxicity severity levels and for transitions between toxicity levels. We compare the operating characteristics of the proposed design with those from a similar design based on a fully-evaluated model that directly models the maximum observed toxicity level within the patients' entire assessment window. We describe settings in which, under comparable power, the proposed design shortens the trial. The proposed design offers more benefit compared to the alternative design as patient accrual becomes slower.

Mechanical constraints imposed by 3D cellular geometry and arrangement modulate growth patterns in the Arabidopsis embryo

Morphogenesis occurs in 3D space over time and is guided by coordinated gene expression programs. Here we use postembryonic development in Arabidopsis plants to investigate the genetic control of growth. We demonstrate that gene expression driving the production of the growth-stimulating hormone gibberellic acid and downstream growth factors is first induced within the radicle tip of the embryo. The center of cell expansion is, however, spatially displaced from the center of gene expression. Because the rapidly growing cells have very different geometry from that of those at the tip, we hypothesized that mechanical factors may contribute to this growth displacement. To this end we developed 3D finite-element method models of growing custom-designed digital embryos at cellular resolution. We used this framework to conceptualize how cell size, shape, and topology influence tissue growth and to explore the interplay of geometrical and genetic inputs into growth distribution. Our simulations showed that mechanical constraints are sufficient to explain the disconnect between the experimentally observed spatiotemporal patterns of gene expression and early postembryonic growth. The center of cell expansion is the position where genetic and mechanical facilitators of growth converge. We have thus uncovered a mechanism whereby 3D cellular geometry helps direct where genetically specified growth takes place.

Relationship between survival and increased radiation dose to subventricular zone in glioblastoma is controversial.

To test the hypothesis on prolonged survival in glioblastoma cases with increased subventricular zone (SVZ) radiation dose. Sixty glioblastoma cases were previously treated with adjuvant radiotherapy and Temozolamide. Ipsilateral, contralateral and bilateral SVZs were contoured and their doses were retrospectively evaluated. Median follow-up, progression free survival (PFS) and overall survival (OS) were 24.5, 8.5 and 19.3 months respectively. Log-rank tests showed a statistically significant correlation between contralateral SVZ (cSVZ) dose > 59.2 Gy (75th percentile) and poor median PFS (10.37 [95% CI 8.37-13.53] vs 7.1 [95% CI 3.5-8.97] months, p = 0.009). cSVZ dose > 59.2 Gy was associated with poor OS in the subgroup with subtotal resection/biopsy (HR: 4.83 [95% CI 1.71-13.97], p = 0.004). High ipsilateral SVZ dose of > 62.25 Gy (75th percentile) was associated with poor PFS in both subgroups of high performance status (HR: 2.58 [95% CI 1.03-6.05], p = 0.044) and SVZ without tumoral contact (HR: 10.57 [95% CI 2.04-49], p = 0.008). The effect of high cSVZ dose on PFS lost its statistical significance in multivariate Cox regression analysis. We report contradictory results compared to previous publications. Changing the clinical practice based on retrospective studies which even do not indicate consistent results among each other will be dangerous. We need carefully designed prospective randomized studies to evaluate any impact of radiation to SVZ in glioblastoma.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Adaptive Beaconless Opportunistic Routing for Multimedia Distribution

User experience on watching live videos must be satisfactory even under the inuence of different network conditions and topology changes, such as happening in Flying Ad-Hoc Networks (FANETs). Routing services for video dissemination over FANETs must be able to adapt routing decisions at runtime to meet Quality of Experience (QoE) requirements. In this paper, we introduce an adaptive beaconless opportunistic routing protocol for video dissemination over FANETs with QoE support, by taking into account multiple types of context information, such as link quality, residual energy, buffer state, as well as geographic information and node mobility in a 3D space. The proposed protocol takes into account Bayesian networks to define weight vectors and Analytic Hierarchy Process (AHP) to adjust the degree of importance for the context information based on instantaneous values. It also includes a position prediction to monitor the distance between two nodes in order to detect possible route failure.

Dose escalated intensity modulated radiotherapy in the treatment of cervical cancer.

PURPOSE Standard dose of external beam radiotherapy seems to be insufficient for satisfactory control of loco-regionally advanced cervical cancer. Aim of our study is to evaluate the outcome as well as early and chronic toxicities in patients with loco-regionally advanced cervical cancer, treated with dose escalated intensity modulated radiotherapy (IMRT) combined with cisplatin chemotherapy. MATERIAL AND METHODS Thirty-nine patients with cervical carcinoma FIGO stage IB2 - IVA were treated with curative intent between 2006 and 2010. The dose of 50.4 Gy was prescribed to the elective pelvic nodal volume. Primary tumors < 4 cm in diameter (n = 6; 15.4 %) received an external beam radiotherapy (EBRT) boost of 5.4 Gy, primary tumors > 4 cm in diameter (n = 33; 84.6 %) received an EBRT boost of 9 Gy. Patients with positive lymph nodes detected with (18)FDG-PET/CT (n = 22; 56.4 %) received a boost to a total dose of 59.4 - 64.8 Gy. The para-aortic region was included in the radiation volume in 8 (20.5 %) patients and in 5 (12.8 %) patients the para-aortic macroscopic lymph nodes received an EBRT boost. IMRT was followed with a 3D planned high dose rate intrauterine brachytherapy given to 36 (92.3 %) patients with a total dose ranging between 15-18 Gy in three fractions (single fraction: 4-6.5 Gy). Patients without contraindications (n = 31/79.5 %) received concomitantly a cisplatin-based chemotherapy (40 mg/kg) weekly. Toxicities were graded according to the common terminology criteria for adverse events (CTCAE v 4.0). RESULTS Mean overall survival for the entire cohort was 61.1 months (±3.5 months). Mean disease free survival was 47.2 months (±4.9 months) and loco-regional disease free survival was 55.2 months (±4.4 months). 65 % of patients developed radiotherapy associated acute toxicities grade 1, ca. 30 % developed toxicities grade 2 and just two (5.2 %) patients developed grade 3 toxicities, one acute diarrhea and one acute cystitis. 16 % of patients had chronic toxicities grade 1, 9 % grade 2 and one patient (2.6 %) toxicities grade 3 in the form of vaginal dryness. CONCLUSION Dose escalated IMRT appears to have a satisfactory outcome with regards to mean overall survival, disease free and loco-regional disease free survival, whereas the treatment-related toxicities remain reasonably low.

The assumptions and effects of conservative procedures in low -dose risk assessment for setting safety standards

Conservative procedures in low-dose risk assessment are used to set safety standards for known or suspected carcinogens. However, the assumptions upon which the methods are based and the effects of these methods are not well understood.^ To minimize the number of false-negatives and to reduce the cost of bioassays, animals are given very high doses of potential carcinogens. Results must then be extrapolated to much smaller doses to set safety standards for risks such as one per million. There are a number of competing methods that add a conservative safety factor into these calculations.^ A method of quantifying the conservatism of these methods was described and tested on eight procedures used in setting low-dose safety standards. The results using these procedures were compared by computer simulation and by the use of data from a large scale animal study.^ The method consisted of determining a "true safe dose" (tsd) according to an assumed underlying model. If one assumed that Y = the probability of cancer = P(d), a known mathematical function of the dose, then by setting Y to some predetermined acceptable risk, one can solve for d, the model's "true safe dose".^ Simulations were generated, assuming a binomial distribution, for an artificial bioassay. The eight procedures were then used to determine a "virtual safe dose" (vsd) that estimates the tsd, assuming a risk of one per million. A ratio R = ((tsd-vsd)/vsd) was calculated for each "experiment" (simulation). The mean R of 500 simulations and the probability R $<$ 0 was used to measure the over and under conservatism of each procedure.^ The eight procedures included Weil's method, Hoel's method, the Mantel-Byran method, the improved Mantel-Byran, Gross's method, fitting a one-hit model, Crump's procedure, and applying Rai and Van Ryzin's method to a Weibull model.^ None of the procedures performed uniformly well for all types of dose-response curves. When the data were linear, the one-hit model, Hoel's method, or the Gross-Mantel method worked reasonably well. However, when the data were non-linear, these same methods were overly conservative. Crump's procedure and the Weibull model performed better in these situations. ^

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.