Improving the accuracy of radiation pneumonitis dose response models

The prognosis for lung cancer patients remains poor. Five year survival rates have been reported to be 15%. Studies have shown that dose escalation to the tumor can lead to better local control and subsequently better overall survival. However, dose to lung tumor is limited by normal tissue toxicity. The most prevalent thoracic toxicity is radiation pneumonitis. In order to determine a safe dose that can be delivered to the healthy lung, researchers have turned to mathematical models predicting the rate of radiation pneumonitis. However, these models rely on simple metrics based on the dose-volume histogram and are not yet accurate enough to be used for dose escalation trials. The purpose of this work was to improve the fit of predictive risk models for radiation pneumonitis and to show the dosimetric benefit of using the models to guide patient treatment planning. The study was divided into 3 specific aims. The first two specifics aims were focused on improving the fit of the predictive model. In Specific Aim 1 we incorporated information about the spatial location of the lung dose distribution into a predictive model. In Specific Aim 2 we incorporated ventilation-based functional information into a predictive pneumonitis model. In the third specific aim a proof of principle virtual simulation was performed where a model-determined limit was used to scale the prescription dose. The data showed that for our patient cohort, the fit of the model to the data was not improved by incorporating spatial information. Although we were not able to achieve a significant improvement in model fit using pre-treatment ventilation, we show some promising results indicating that ventilation imaging can provide useful information about lung function in lung cancer patients. The virtual simulation trial demonstrated that using a personalized lung dose limit derived from a predictive model will result in a different prescription than what was achieved with the clinically used plan; thus demonstrating the utility of a normal tissue toxicity model in personalizing the prescription dose.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Mouse models and quantitative trait loci of susceptibility to bleomycin- and radiation-induced pulmonary fibrosis

Radiotherapy involving the thoracic cavity and chemotherapy with the drug bleomycin are both dose limited by the development of pulmonary fibrosis. From evidence that there is variation in the population in susceptibility to pulmonary fibrosis, and animal data, it was hypothesized that individual variation in susceptibility to bleomycin-induced, or radiation-induced, pulmonary fibrosis is, in part, genetically controlled. In this thesis a three generation mouse genetic model of C57BL/6J (fibrosis prone) and C3Hf/Kam (fibrosis resistant) mouse strains and F1 and F2 (F1 intercross) progeny derived from the parental strains was developed to investigate the genetic basis of susceptibility to fibrosis. In the bleomycin studies the mice received 100 mg/kg (125 for females) of bleomycin, via mini osmotic pump. The animals were sacrificed at eight weeks following treatment or when their breathing rate indicated respiratory distress. In the radiation studies the mice were given a single dose of 14 or 16 Gy (Co$\sp{60})$ to the whole thorax and were sacrificed when moribund. The phenotype was defined as the percent of fibrosis area in the left lung as quantified with image analysis of histological sections. Quantitative trait loci (QTL) mapping was used to identify the chromosomal location of genes which contribute to susceptibility to bleomycin-induced pulmonary fibrosis in C57BL/6J mice compared to C3Hf/Kam mice and to determine if the QTL's which influence susceptibility to bleomycin-induced lung fibrosis in these progenitor strains could be implicated in susceptibility to radiation-induced lung fibrosis. For bleomycin, a genome wide scan revealed QTL's on chromosome 17, at the MHC, (LOD = 11.7 for males and 7.2 for females) accounting for approximately 21% of the phenotypic variance, and on chromosome 11 (LOD = 4.9), in male mice only, adding 8% of phenotypic variance. The bleomycin QTL on chromosome 17 was also implicated for susceptibility to radiation-induced fibrosis (LOD = 5.0) and contributes 7% of the phenotypic variance in the radiation study. In conclusion, susceptibility to both bleomycin-induced and radiation-induced pulmonary fibrosis are heritable traits, and are influenced by a genetic factor which maps to a genomic region containing the MHC. ^

Structural equation modeling of the medical outcome study: Short Form 36 (SF-36)

The factorial validity of the SF-36 was evaluated using confirmatory factor analysis (CFA) methods, structural equation modeling (SEM), and multigroup structural equation modeling (MSEM). First, the measurement and structural model of the hypothesized SF-36 was explicated. Second, the model was tested for the validity of a second-order factorial structure, upon evidence of model misfit, determined the best-fitting model, and tested the validity of the best-fitting model on a second random sample from the same population. Third, the best-fitting model was tested for invariance of the factorial structure across race, age, and educational subgroups using MSEM.^ The findings support the second-order factorial structure of the SF-36 as proposed by Ware and Sherbourne (1992). However, the results suggest that: (a) Mental Health and Physical Health covary; (b) general mental health cross-loads onto Physical Health; (c) general health perception loads onto Mental Health instead of Physical Health; (d) many of the error terms are correlated; and (e) the physical function scale is not reliable across these two samples. This hierarchical factor pattern was replicated across both samples of health care workers, suggesting that the post hoc model fitting was not data specific. Subgroup analysis suggests that the physical function scale is not reliable across the "age" or "education" subgroups and that the general mental health scale path from Mental Health is not reliable across the "white/nonwhite" or "education" subgroups.^ The importance of this study is in the use of SEM and MSEM in evaluating sample data from the use of the SF-36. These methods are uniquely suited to the analysis of latent variable structures and are widely used in other fields. The use of latent variable models for self reported outcome measures has become widespread, and should now be applied to medical outcomes research. Invariance testing is superior to mean scores or summary scores when evaluating differences between groups. From a practical, as well as, psychometric perspective, it seems imperative that construct validity research related to the SF-36 establish whether this same hierarchical structure and invariance holds for other populations.^ This project is presented as three articles to be submitted for publication. ^

Bending Shear: The Rate-Controlling Mechanism for Calving Ice Walls

Bending shear was observed to produce nearly vertical shear bands in a calving ice wall standing on dry land on Deception Island (Iat. 63.0 oS., long. 60.6 W.), and slabs calved straight downward when shear rupture occurred along these shear bands (Hughes, 1989). A formula for the calving rate was developed from the Deception Island data, and we have attempted to justify generalizing this formula to include ice walls standing along beaches or in water. These are environments in which a wave-washed groove develops along the base of the ice wall or along a water line above the base. The rate of wave erosion provides an alternative mechanism for controlling the calving rate in these environments. We have determined that the rate at which bending creep produces nearly vertical shear bands, along which shear r upture occurs, controls the calving rate in all environments. Shear rupture occurs at a calving shear stress of about I bar. Our results justify using the calving formula to compute the calving rate of ice walls in computer models of ice-sheet dynamics. This is especially important in simulating retreat of Northern Hemisphere ice sheets during the last deglaciation, when marine and lacustrine environments were common along retreating ice margins. These margins would have been ice walls standing along beaches or in water, because floating ice shelves are not expected in the ablation zone of retreating ice sheets.

Heart rate is a prognostic risk factor for myocardial infarction: a post hoc analysis in the PERFORM (Prevention of cerebrovascular and cardiovascular Events of ischemic origin with teRutroban in patients with a history oF ischemic strOke or tRansient ischeMic attack) study population

BACKGROUND Elevated resting heart rate is known to be detrimental to morbidity and mortality in cardiovascular disease, though its effect in patients with ischemic stroke is unclear. We analyzed the effect of baseline resting heart rate on myocardial infarction (MI) in patients with a recent noncardioembolic cerebral ischemic event participating in PERFORM. METHODS We compared fatal or nonfatal MI using adjusted Cox proportional hazards models for PERFORM patients with baseline heart rate <70 bpm (n=8178) or ≥70 bpm (n=10,802). In addition, heart rate was analyzed as a continuous variable. Other cerebrovascular and cardiovascular outcomes were also explored. RESULTS Heart rate ≥70 bpm was associated with increased relative risk for fatal or nonfatal MI (HR 1.32, 95% CI 1.03-1.69, P=0.029). For every 5-bpm increase in heart rate, there was an increase in relative risk for fatal and nonfatal MI (11.3%, P=0.0002). Heart rate ≥70 bpm was also associated with increased relative risk for a composite of fatal or nonfatal ischemic stroke, fatal or nonfatal MI, or other vascular death (excluding hemorrhagic death) (P<0001); vascular death (P<0001); all-cause mortality (P<0001); and fatal or nonfatal stroke (P=0.04). For every 5-bpm increase in heart rate, there were increases in relative risk for fatal or nonfatal ischemic stroke, fatal or nonfatal MI, or other vascular death (4.7%, P<0.0001), vascular death (11.0%, P<0.0001), all-cause mortality (8.0%, P<0.0001), and fatal and nonfatal stroke (2.4%, P=0.057). CONCLUSION Elevated heart rate ≥70 bpm places patients with a noncardioembolic cerebral ischemic event at increased risk for MI.

Factors other than the glomerular filtration rate that determine the serum beta-2-microglobulin level.

BACKGROUND β2-microglobulin has been increasingly investigated as a diagnostic marker of kidney function and a prognostic marker of adverse outcomes. To date, non-renal determinants of β2-microglobulin levels have not been well described. Non-renal determinants are important for the interpretation and appraisal of the diagnostic and prognostic value of any endogenous kidney function marker. METHODS This cross-sectional analysis was performed within the framework of the www.seniorlabor.ch study, which includes subjectively healthy individuals aged ≥ 60 years. Factors known or suspected to have a non-renal association with kidney function markers were investigated for a non-renal association with serum β2-microglobulin. As a marker of kidney function, the Berlin Initiative Study equation 2 for the estimation of the estimated glomerular filtration rate (eGFR(BIS2)) in the elderly was employed. RESULTS A total of 1302 participants (714 females and 588 males) were enrolled in the study. The use of a multivariate regression model adjusting for age, gender and kidney function (eGFR(BIS2)) revealed age, male gender, and C-reactive protein level to be positively associated with β2-microglobulin levels. In addition, there was an inverse non-renal relationship between systolic blood pressure, total cholesterol and current smoking status. No association with markers of diabetes mellitus, body stature, nutritional risk, thyroid function or calcium and phosphate levels was observed. CONCLUSIONS Serum β2-microglobulin levels in elderly subjects are related to several non-renal factors. These non-renal factors are not congruent to those known from other markers (i.e. cystatin C and creatinine) and remind of classical cardiovascular risk factors.

Organ donation in Switzerland - an analysis of factors associated with consent rate

BACKGROUND AND AIM Switzerland has a low post mortem organ donation rate. Here we examine variables that are associated with the consent of the deceased's next of kin (NOK) for organ donation, which is a prerequisite for donation in Switzerland. METHODS AND ANALYSIS During one year, we registered information from NOK of all deceased patients in Swiss intensive care units, who were approached for consent to organ donation. We collected data on patient demographics, characteristics of NOK, factors related to the request process and to the clinical setting. We analyzed the association of collected predictors with consent rate using univariable logistic regression models; predictors with p-values <0.2 were selected for a multivariable logistic regression. RESULTS Of 266 NOK approached for consent, consent was given in 137 (51.5%) cases. In multivariable analysis, we found associations of consent rates with Swiss nationality (OR 3.09, 95% CI: 1.46-6.54) and German language area (OR 0.31, 95% CI: 0.14-0.73). Consent rates tended to be higher if a parent was present during the request (OR 1.76, 95% CI: 0.93-3.33) and if the request was done before brain death was formally declared (OR 1.87, 95% CI: 0.90-3.87). CONCLUSION Establishing an atmosphere of trust between the medical staff putting forward a request and the NOK, allowing sufficient time for the NOK to consider donation, and respecting personal values and cultural differences, could be of importance for increasing donation rates. Additional measures are needed to address the pronounced differences in consent rates between language regions.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

Incidence Rate of Kaposi Sarcoma in HIV-Infected Patients on Antiretroviral Therapy in Southern Africa: A Prospective Multicohort Study.

BACKGROUND The risk of Kaposi sarcoma (KS) among HIV-infected persons on antiretroviral therapy (ART) is not well defined in resource-limited settings. We studied KS incidence rates and associated risk factors in children and adults on ART in Southern Africa. METHODS We included patient data of 6 ART programs in Botswana, South Africa, Zambia, and Zimbabwe. We estimated KS incidence rates in patients on ART measuring time from 30 days after ART initiation to KS diagnosis, last follow-up visit, or death. We assessed risk factors (age, sex, calendar year, WHO stage, tuberculosis, and CD4 counts) using Cox models. FINDINGS We analyzed data from 173,245 patients (61% female, 8% children aged <16 years) who started ART between 2004 and 2010. Five hundred and sixty-four incident cases were diagnosed during 343,927 person-years (pys). The overall KS incidence rate was 164/100,000 pys [95% confidence interval (CI): 151 to 178]. The incidence rate was highest 30-90 days after ART initiation (413/100,000 pys; 95% CI: 342 to 497) and declined thereafter [86/100,000 pys (95% CI: 71 to 105), >2 years after ART initiation]. Male sex [adjusted hazard ratio (HR): 1.34; 95% CI: 1.12 to 1.61], low current CD4 counts (≥500 versus <50 cells/μL, adjusted HR: 0.36; 95% CI: 0.23 to 0.55), and age (5-9 years versus 30-39 years, adjusted HR: 0.20; 95% CI: 0.05 to 0.79) were relevant risk factors for developing KS. INTERPRETATION Despite ART, KS risk in HIV-infected persons in Southern Africa remains high. Early HIV testing and maintaining high CD4 counts is needed to further reduce KS-related morbidity and mortality.

All-loop anomalous dimensions in integrable λ-deformed σ-models

We calculate the all-loop anomalous dimensions of current operators in λ-deformed σ-models. For the isotropic integrable deformation and for a semi-simple group G we compute the anomalous dimensions using two different methods. In the first we use the all-loop effective action and in the second we employ perturbation theory along with the Callan–Symanzik equation and in conjunction with a duality-type symmetry shared by these models. Furthermore, using CFT techniques we compute the all-loop anomalous dimension of bilinear currents for the isotropic deformation case and a general G . Finally we work out the anomalous dimension matrix for the cases of anisotropic SU(2) and the two couplings, corresponding to the symmetric coset G/H and a subgroup H, splitting of a group G.

Plasma environment of a weak comet – Predictions for Comet 67P/Churyumov–Gerasimenko from multifluid-MHD and Hybrid models

The interaction of a comet with the solar wind undergoes various stages as the comet’s activity varies along its orbit. For a comet like 67P/Churyumov–Gerasimenko, the target comet of ESA’s Rosetta mission, the various features include the formation of a Mach cone, the bow shock, and close to perihelion even a diamagnetic cavity. There are different approaches to simulate this complex interplay between the solar wind and the comet’s extended neutral gas coma which include magnetohydrodynamics (MHD) and hybrid-type models. The first treats the plasma as fluids (one fluid in basic single fluid MHD) and the latter treats the ions as individual particles under the influence of the local electric and magnetic fields. The electrons are treated as a charge-neutralizing fluid in both cases. Given the different approaches both models yield different results, in particular for a low production rate comet. In this paper we will show that these differences can be reduced when using a multifluid instead of a single-fluid MHD model and increase the resolution of the Hybrid model. We will show that some major features obtained with a hybrid type approach like the gyration of the cometary heavy ions and the formation of the Mach cone can be partially reproduced with the multifluid-type model.

Comparison of 3D kinetic and hydrodynamic models to ROSINA-COPS measurements of the neutral coma of 67P/Churyumov-Gerasimenko

67P/Churyumov-Gerasimenko (67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P’s neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 10²⁶ molecules s⁻¹.

Exclusion of leptophilic dark matter models using XENON100 electronic recoil data

We have searched for periodic variations of the electronic recoil event rate in the (2-6) keV energy range recorded between February 2011 and March 2012 with the XENON100 detector, adding up to 224.6 live days in total. Following a detailed study to establish the stability of the detector and its background contributions during this run, we performed an un-binned profile likelihood analysis to identify any periodicity up to 500 days. We find a global significance of less than 1 sigma for all periods suggesting no statistically significant modulation in the data. While the local significance for an annual modulation is 2.8 sigma, the analysis of a multiple-scatter control sample and the phase of the modulation disfavor a dark matter interpretation. The DAMA/LIBRA annual modulation interpreted as a dark matter signature with axial-vector coupling of WIMPs to electrons is excluded at 4.8 sigma.