An efficient framework for photon Monte Carlo treatment planning

Currently photon Monte Carlo treatment planning (MCTP) for a patient stored in the patient database of a treatment planning system (TPS) can usually only be performed using a cumbersome multi-step procedure where many user interactions are needed. This means automation is needed for usage in clinical routine. In addition, because of the long computing time in MCTP, optimization of the MC calculations is essential. For these purposes a new graphical user interface (GUI)-based photon MC environment has been developed resulting in a very flexible framework. By this means appropriate MC transport methods are assigned to different geometric regions by still benefiting from the features included in the TPS. In order to provide a flexible MC environment, the MC particle transport has been divided into different parts: the source, beam modifiers and the patient. The source part includes the phase-space source, source models and full MC transport through the treatment head. The beam modifier part consists of one module for each beam modifier. To simulate the radiation transport through each individual beam modifier, one out of three full MC transport codes can be selected independently. Additionally, for each beam modifier a simple or an exact geometry can be chosen. Thereby, different complexity levels of radiation transport are applied during the simulation. For the patient dose calculation, two different MC codes are available. A special plug-in in Eclipse providing all necessary information by means of Dicom streams was used to start the developed MC GUI. The implementation of this framework separates the MC transport from the geometry and the modules pass the particles in memory; hence, no files are used as the interface. The implementation is realized for 6 and 15 MV beams of a Varian Clinac 2300 C/D. Several applications demonstrate the usefulness of the framework. Apart from applications dealing with the beam modifiers, two patient cases are shown. Thereby, comparisons are performed between MC calculated dose distributions and those calculated by a pencil beam or the AAA algorithm. Interfacing this flexible and efficient MC environment with Eclipse allows a widespread use for all kinds of investigations from timing and benchmarking studies to clinical patient studies. Additionally, it is possible to add modules keeping the system highly flexible and efficient.

Testing the GlaaS algorithm for dose measurements on low- and high-energy photon beams using an amorphous silicon portal imager

The GLAaS algorithm for pretreatment intensity modulation radiation therapy absolute dose verification based on the use of amorphous silicon detectors, as described in Nicolini et al. [G. Nicolini, A. Fogliata, E. Vanetti, A. Clivio, and L. Cozzi, Med. Phys. 33, 2839-2851 (2006)], was tested under a variety of experimental conditions to investigate its robustness, the possibility of using it in different clinics and its performance. GLAaS was therefore tested on a low-energy Varian Clinac (6 MV) equipped with an amorphous silicon Portal Vision PV-aS500 with electronic readout IAS2 and on a high-energy Clinac (6 and 15 MV) equipped with a PV-aS1000 and IAS3 electronics. Tests were performed for three calibration conditions: A: adding buildup on the top of the cassette such that SDD-SSD = d(max) and comparing measurements with corresponding doses computed at d(max), B: without adding any buildup on the top of the cassette and considering only the intrinsic water-equivalent thickness of the electronic portal imaging devices device (0.8 cm), and C: without adding any buildup on the top of the cassette but comparing measurements against doses computed at d(max). This procedure is similar to that usually applied when in vivo dosimetry is performed with solid state diodes without sufficient buildup material. Quantitatively, the gamma index (gamma), as described by Low et al. [D. A. Low, W. B. Harms, S. Mutic, and J. A. Purdy, Med. Phys. 25, 656-660 (1998)], was assessed. The gamma index was computed for a distance to agreement (DTA) of 3 mm. The dose difference deltaD was considered as 2%, 3%, and 4%. As a measure of the quality of results, the fraction of field area with gamma larger than 1 (%FA) was scored. Results over a set of 50 test samples (including fields from head and neck, breast, prostate, anal canal, and brain cases) and from the long-term routine usage, demonstrated the robustness and stability of GLAaS. In general, the mean values of %FA remain below 3% for deltaD equal or larger than 3%, while they are slightly larger for deltaD = 2% with %FA in the range from 3% to 8%. Since its introduction in routine practice, 1453 fields have been verified with GLAaS at the authors' institute (6 MV beam). Using a DTA of 3 mm and a deltaD of 4% the authors obtained %FA = 0.9 +/- 1.1 for the entire data set while, stratifying according to the dose calculation algorithm, they observed: %FA = 0.7 +/- 0.9 for fields computed with the analytical anisotropic algorithm and %FA = 2.4 +/- 1.3 for pencil-beam based fields with a statistically significant difference between the two groups. If data are stratified according to field splitting, they observed %FA = 0.8 +/- 1.0 for split fields and 1.0 +/- 1.2 for nonsplit fields without any significant difference.

Analysis and performance of a UPC implementation of a parallel longest common subsequence algorithm

An important problem in computational biology is finding the longest common subsequence (LCS) of two nucleotide sequences. This paper examines the correctness and performance of a recently proposed parallel LCS algorithm that uses successor tables and pruning rules to construct a list of sets from which an LCS can be easily reconstructed. Counterexamples are given for two pruning rules that were given with the original algorithm. Because of these errors, performance measurements originally reported cannot be validated. The work presented here shows that speedup can be reliably achieved by an implementation in Unified Parallel C that runs on an Infiniband cluster. This performance is partly facilitated by exploiting the software cache of the MuPC runtime system. In addition, this implementation achieved speedup without bulk memory copy operations and the associated programming complexity of message passing.

Experimental model to evaluate in vivo and in vitro cartilage MR imaging by means of histological analyses

OBJECTIVES: Implementation of an experimental model to compare cartilage MR imaging by means of histological analyses. MATERIAL AND METHODS: MRI was obtained from 4 patients expecting total knee replacement at 1.5 and/or 3T prior surgery. The timeframe between pre-op MRI and knee replacement was within two days. Resected cartilage-bone samples were tagged with Ethi((R))-pins to reproduce the histological cutting course. Pre-operative scanning at 1.5T included following parameters for fast low angle shot (FLASH: TR/TE/FA=33ms/6ms/30 degrees , BW=110kHz, 120mmx120mm FOV, 256x256 matrix, 0.65mm slice-thickness) and double echo steady state (DESS: TR/TE/FA=23.7ms/6.9ms/40 degrees , BW=130kHz, 120x120mm FOV, 256x256 matrix, 0.65mm slice-thickness). At 3T, scan parameters were: FLASH (TR/TE/FA=12.2ms/5.1ms/10 degrees , BW=130kHz, 170x170mm FOV, 320x320, 0.5mm slice-thickness) and DESS (TR/TE/FA=15.6ms/4.5ms/25 degrees , BW=200kHz, 135mmx150mm FOV, 288x320matrix, 0.5mm slice-thickness). Imaging of the specimens was done the same day at 1.5T. MRI (Noyes) and histological (Mankin) score scales were correlated using the paired t-test. Sensitivity and specificity for the detection of different grades of cartilage degeneration were assessed. Inter-reader and intra-reader reliability was determined using Kappa analysis. RESULTS: Low correlation (sensitivity, specificity) was found for both sequences in normal to mild Mankin grades. Only moderate to severe changes were diagnosed with higher significance and specificity. The use of higher field-strengths was advantageous for both protocols with sensitivity values ranging from 13.6% to 93.3% (FLASH) and 20.5% to 96.2% (DESS). Kappa values ranged from 0.488 to 0.944. CONCLUSIONS: Correlating MR images with continuous histological slices was feasible by using three-dimensional imaging, multi-planar-reformat and marker pins. The capability of diagnosing early cartilage changes with high accuracy could not be proven for both FLASH and DESS.

Control design and genetic algorithm optimization for electrostatic MEMS

This dissertation discusses structural-electrostatic modeling techniques, genetic algorithm based optimization and control design for electrostatic micro devices. First, an alternative modeling technique, the interpolated force model, for electrostatic micro devices is discussed. The method provides improved computational efficiency relative to a benchmark model, as well as improved accuracy for irregular electrode configurations relative to a common approximate model, the parallel plate approximation model. For the configuration most similar to two parallel plates, expected to be the best case scenario for the approximate model, both the parallel plate approximation model and the interpolated force model maintained less than 2.2% error in static deflection compared to the benchmark model. For the configuration expected to be the worst case scenario for the parallel plate approximation model, the interpolated force model maintained less than 2.9% error in static deflection while the parallel plate approximation model is incapable of handling the configuration. Second, genetic algorithm based optimization is shown to improve the design of an electrostatic micro sensor. The design space is enlarged from published design spaces to include the configuration of both sensing and actuation electrodes, material distribution, actuation voltage and other geometric dimensions. For a small population, the design was improved by approximately a factor of 6 over 15 generations to a fitness value of 3.2 fF. For a larger population seeded with the best configurations of the previous optimization, the design was improved by another 7% in 5 generations to a fitness value of 3.0 fF. Third, a learning control algorithm is presented that reduces the closing time of a radiofrequency microelectromechanical systems switch by minimizing bounce while maintaining robustness to fabrication variability. Electrostatic actuation of the plate causes pull-in with high impact velocities, which are difficult to control due to parameter variations from part to part. A single degree-of-freedom model was utilized to design a learning control algorithm that shapes the actuation voltage based on the open/closed state of the switch. Experiments on 3 test switches show that after 5-10 iterations, the learning algorithm lands the switch with an impact velocity not exceeding 0.2 m/s, eliminating bounce.

The Cardica C-Port System: clinical and angiographic evaluation of a new device for automated, compliant distal anastomoses in coronary artery bypass grafting surgery--a multicenter prospective clinical trial

OBJECTIVES: The C-Port System (Cardica, Inc, Redwood City, Calif) integrates in one tool all functions necessary to enable rapid automated distal coronary anastomoses. The goal of this prospective, nonrandomized, and multicenter study is to determine the safety and efficacy of this novel anastomotic system. METHODS: Five centers enrolled 133 patients awaiting elective coronary artery bypass grafting surgery. Outcome variables were intraoperative device performance, incidence of device-related adverse events, predischarge and 6-month angiographic graft patency, and 12-month clinical outcome. Independent core laboratories performed qualitative and quantitative angiographic and computed tomographic assessments. RESULTS: The C-Port was used to perform a vein-to-coronary anastomosis in 130 patients. Intraoperative conversion to a hand-sewn anastomosis was necessary in 11 patients because of inadequate target site preparation, inappropriate target vessel selection, or both. Inadequate blood flow related to poor runoff required conversion in 3 additional patients. Three patients died before discharge of causes unrelated to the device. At discharge, 113 patients had a C-Port implant in place, and 104 C-Port anastomoses were studied by means of angiography, resulting in 100 FitzGibbon A, 3 FitzGibbon B, and 1 FitzGibbon 0 classifications. At 6 months, one additional patient died of a device-unrelated cause, and 98 patients were evaluated by means of angiography (n = 89). Overall patency (FitzGibbon A) was 92.1%. Three C-Port anastomoses were rated FitzGibbon B, and 4 were rated FitzGibbon 0. At 12 months, 107 (98.2%) of 109 alive patients were followed up, without any reports of device-related major adverse cardiac events. CONCLUSIONS: The C-Port System allows for a rapid, reliable, and compliant distal anastomosis and yields favorable 6-month angiographic and 12-month clinical results when compared with published studies.

Radiation metabolomics. 3. Biomarker discovery in the urine of gamma-irradiated rats using a simplified metabolomics protocol of gas chromatography-mass spectrometry combined with random forests machine learning algorithm

Abstract Radiation metabolomics employing mass spectral technologies represents a plausible means of high-throughput minimally invasive radiation biodosimetry. A simplified metabolomics protocol is described that employs ubiquitous gas chromatography-mass spectrometry and open source software including random forests machine learning algorithm to uncover latent biomarkers of 3 Gy gamma radiation in rats. Urine was collected from six male Wistar rats and six sham-irradiated controls for 7 days, 4 prior to irradiation and 3 after irradiation. Water and food consumption, urine volume, body weight, and sodium, potassium, calcium, chloride, phosphate and urea excretion showed major effects from exposure to gamma radiation. The metabolomics protocol uncovered several urinary metabolites that were significantly up-regulated (glyoxylate, threonate, thymine, uracil, p-cresol) and down-regulated (citrate, 2-oxoglutarate, adipate, pimelate, suberate, azelaate) as a result of radiation exposure. Thymine and uracil were shown to derive largely from thymidine and 2'-deoxyuridine, which are known radiation biomarkers in the mouse. The radiation metabolomic phenotype in rats appeared to derive from oxidative stress and effects on kidney function. Gas chromatography-mass spectrometry is a promising platform on which to develop the field of radiation metabolomics further and to assist in the design of instrumentation for use in detecting biological consequences of environmental radiation release.

Severe autoimmune cytopenias in treatment-naive hepatitis C virus infection: clinical description of 16 cases

OBJECTIVE: To describe the prevalence, main characteristics, and treatment of severe autoimmune cytopenias [autoimmune hemolytic anemia (AIHA), autoimmune thrombocytopenic purpura (AITP)] in patients with chronic hepatitis C virus (HCV) infection. METHODS: Retrospective chart review of patients with chronic HCV infection seen at our institution. Two additional departments contributed eight more patients to assess therapy of HCV-related autoimmune cytopenias. RESULTS: Eight patients (seven AITP, one AIHA) fulfilled the inclusion criteria in our population of 4345 HCV-infected patients. The number of patients with AITP was much greater than would be expected by chance (P<0.0001). Patients with HCV-related AITP were older and demonstrated more immunological markers than a group of 40 controls. Eight additional patients (six AITP, two Evans syndrome) were included. We only assessed the response for AITP patients because of the single case of AIHA. Patients with AITP had a poor response to initial corticosteroids [one complete response (CR), three partial response (PR), and four failures]. Intravenous immunoglobulins led to transient efficacy in three of four patients. In second-line therapy, five of seven patients responded to splenectomy. Rituximab proved effective in increasing platelets in two patients. Of eight patients treated with antiviral therapy (IFN-alpha+/-ribavirin), five responded (three CR, two PR). CONCLUSION: AITP occurs more commonly in patients with chronic HCV infection than would be expected by chance. HCV-positive AITP requires a treatment strategy different from that used in HCV-negative AITP. On the basis of the results from our study and a literature analysis, we propose an algorithm for treatment of severe HCV-related autoimmune cytopenias.

Tissue properties of the human vertebral body sub-structures evaluated by means of microindentation

Purpose The better understanding of vertebral mechanical properties can help to improve the diagnosis of vertebral fractures. As the bone mechanical competence depends not only from bone mineral density (BMD) but also from bone quality, the goal of the present study was to investigate the anisotropic indentation moduli of the different sub-structures of the healthy human vertebral body and spondylophytes by means of microindentation. Methods Six human vertebral bodies and five osteophytes (spondylophytes) were collected and prepared for microindentation test. In particular, indentations were performed on bone structural units of the cortical shell (along axial, circumferential and radial directions), of the endplates (along the anterio-posterior and lateral directions), of the trabecular bone (along the axial and transverse directions) and of the spondylophytes (along the axial direction). A total of 3164 indentations down to a maximum depth of 2.5 µm were performed and the indentation modulus was computed for each measurement. Results The cortical shell showed an orthotropic behavior (indentation modulus, Ei, higher if measured along the axial direction, 14.6±2.8 GPa, compared to the circumferential one, 12.3±3.5 GPa, and radial one, 8.3±3.1 GPa). Moreover, the cortical endplates (similar Ei along the antero-posterior, 13.0±2.9 GPa, and along the lateral, 12.0±3.0 GPa, directions) and the trabecular bone (Ei= 13.7±3.4 GPa along the axial direction versus Ei=10.9±3.7 GPa along the transverse one) showed transversal isotropy behavior. Furthermore, the spondylophytes showed the lower mechanical properties measured along the axial direction (Ei=10.5±3.3 GPa). Conclusions The original results presented in this study improve our understanding of vertebral biomechanics and can be helpful to define the material properties of the vertebral substructures in computational models such as FE analysis.

Thermal acclimation to 4 or 10 degrees C imparts minimal benefit on swimming performance in Atlantic cod (Gadus morhua L.).

Thermal acclimation is frequently cited as a means by which ectothermic animals improve their Darwinian fitness, i.e. the beneficial acclimation hypothesis. As the critical swimming speed (U (crit)) test is often used as a proxy measure of fitness, we acclimated Atlantic cod (Gadus morhua) to 4 and 10 degrees C and then assessed their U (crit) swimming performance at their respective acclimation temperatures and during acute temperature reversal. Because phenotypic differences exist between different populations of cod, we undertook these experiments in two different populations, North Sea cod and North East Arctic cod. Acclimation to 4 or 10 degrees C had a minimal effect on swimming performance or U (crit), however test temperature did, with all groups having a 10-17% higher U (crit) at 10 degrees C. The swimming efficiency was significantly lower in all groups at 4 degrees C arguably due to the compression of the muscle fibre recruitment order. This also led to a reduction in the duration of "kick and glide" swimming at 4 degrees C. No significant differences were seen between the two populations in any of the measured parameters, due possibly to the extended acclimation period. Our data indicate that acclimation imparts little benefit on U (crit) swimming test in Atlantic cod. Further efforts need to identify the functional consequences of the long-term thermal acclimation process.

Are Greenhouse Gas Signals of Northern Hemisphere winter extra-tropical cyclone activity dependent on the identification and tracking algorithm?

For Northern Hemisphere extra-tropical cyclone activity, the dependency of a potential anthropogenic climate change signal on the identification method applied is analysed. This study investigates the impact of the used algorithm on the changing signal, not the robustness of the climate change signal itself. Using one single transient AOGCM simulation as standard input for eleven state-of-the-art identification methods, the patterns of model simulated present day climatologies are found to be close to those computed from re-analysis, independent of the method applied. Although differences in the total number of cyclones identified exist, the climate change signals (IPCC SRES A1B) in the model run considered are largely similar between methods for all cyclones. Taking into account all tracks, decreasing numbers are found in the Mediterranean, the Arctic in the Barents and Greenland Seas, the mid-latitude Pacific and North America. Changing patterns are even more similar, if only the most severe systems are considered: the methods reveal a coherent statistically significant increase in frequency over the eastern North Atlantic and North Pacific. We found that the differences between the methods considered are largely due to the different role of weaker systems in the specific methods.