Semiparametric Normal Transformation Models for Spatially Correlated Survival Data

There is an emerging interest in modeling spatially correlated survival data in biomedical and epidemiological studies. In this paper, we propose a new class of semiparametric normal transformation models for right censored spatially correlated survival data. This class of models assumes that survival outcomes marginally follow a Cox proportional hazard model with unspecified baseline hazard, and their joint distribution is obtained by transforming survival outcomes to normal random variables, whose joint distribution is assumed to be multivariate normal with a spatial correlation structure. A key feature of the class of semiparametric normal transformation models is that it provides a rich class of spatial survival models where regression coefficients have population average interpretation and the spatial dependence of survival times is conveniently modeled using the transformed variables by flexible normal random fields. We study the relationship of the spatial correlation structure of the transformed normal variables and the dependence measures of the original survival times. Direct nonparametric maximum likelihood estimation in such models is practically prohibited due to the high dimensional intractable integration of the likelihood function and the infinite dimensional nuisance baseline hazard parameter. We hence develop a class of spatial semiparametric estimating equations, which conveniently estimate the population-level regression coefficients and the dependence parameters simultaneously. We study the asymptotic properties of the proposed estimators, and show that they are consistent and asymptotically normal. The proposed method is illustrated with an analysis of data from the East Boston Ashma Study and its performance is evaluated using simulations.

Antitumor activity of an epithelial cell adhesion molecule targeted nanovesicular drug delivery system

Site-specific delivery of anticancer agents to tumors represents a promising therapeutic strategy because it increases efficacy and reduces toxicity to normal tissues compared with untargeted drugs. Sterically stabilized immunoliposomes (SIL), guided by antibodies that specifically bind to well internalizing antigens on the tumor cell surface, are effective nanoscale delivery systems capable of accumulating large quantities of anticancer agents at the tumor site. The epithelial cell adhesion molecule (EpCAM) holds major promise as a target for antibody-based cancer therapy due to its abundant expression in many solid tumors and its limited distribution in normal tissues. We generated EpCAM-directed immunoliposomes by covalently coupling the humanized single-chain Fv antibody fragment 4D5MOCB to the surface of sterically stabilized liposomes loaded with the anticancer agent doxorubicin. In vitro, the doxorubicin-loaded immunoliposomes (SIL-Dox) showed efficient cell binding and internalization and were significantly more cytotoxic against EpCAM-positive tumor cells than nontargeted liposomes (SL-Dox). In athymic mice bearing established human tumor xenografts, pharmacokinetic and biodistribution analysis of SIL-Dox revealed long circulation times in the blood with a half-life of 11 h and effective time-dependent tumor localization, resulting in up to 15% injected dose per gram tissue. These favorable pharmacokinetic properties translated into potent antitumor activity, which resulted in significant growth inhibition (compared with control mice), and was more pronounced than that of doxorubicin alone and nontargeted SL-Dox at low, nontoxic doses. Our data show the promise of EpCAM-directed nanovesicular drug delivery for targeted therapy of solid tumors.

Distribution of amyloid precursor protein and amyloid-beta immunoreactivity in DBA/2J glaucomatous mouse retinas

PURPOSE: Evidence suggests that altered metabolism of amyloid precursor protein (APP) may play a role in the pathophysiology of retinal ganglion cell (RGC) death in the etiology of glaucoma. The authors sought to determine the distribution of APP and amyloid-beta (Abeta) in DBA/2J glaucomatous mouse retinas. METHODS: The retinas of 3- and 15-month-old DBA/2J mice and C57/BL-6 mice (control group) were fixed with 4% paraformaldehyde and processed for immunohistochemistry. Antibodies used included a polyclonal antibody to the C terminus of Abeta 40 and a polyclonal antibody to the APP ectodomain. Immunohistochemically stained tissue was graded using light microscopy. Distribution and semiquantitative expression of APP and Abeta in young and old glaucomatous and normal retinas were determined and compared. RESULTS: Strong APP and Abeta immunoreactivity was found in the RGC layer, optic nerve, and pia/dura of old DBA/2J retinas, with considerably higher intensity found in the old compared with the young DBA/2J mice. In contrast to glaucomatous mice, the control group did not show any notable age-related difference. CONCLUSIONS: Disruption of the homeostatic properties of secreted APP with consecutive Abeta cytotoxicity might be a contributing factor of ganglion cell loss in glaucomatous mouse retinas.

Online distribution management framework by incorporating wind generation

In distribution system operations, dispatchers at control center closely monitor system operating limits to ensure system reliability and adequacy. This reliability is partly due to the provision of remote controllable tie and sectionalizing switches. While the stochastic nature of wind generation can impact the level of wind energy penetration in the network, an estimate of the output from wind on hourly basis can be extremely useful. Under any operating conditions, the switching actions require human intervention and can be an extremely stressful task. Currently, handling a set of switching combinations with the uncertainty of distributed wind generation as part of the decision variables has been nonexistent. This thesis proposes a three-fold online management framework: (1) prediction of wind speed, (2) estimation of wind generation capacity, and (3) enumeration of feasible switching combinations. The proposed methodology is evaluated on 29-node test system with 8 remote controllable switches and two wind farms of 18MW and 9MW nameplate capacities respectively for generating the sequence of system reconfiguration states during normal and emergency conditions.

Bayesian change-point analysis in linear regression model with scale mixtures of normal distributions

In this thesis, we consider Bayesian inference on the detection of variance change-point models with scale mixtures of normal (for short SMN) distributions. This class of distributions is symmetric and thick-tailed and includes as special cases: Gaussian, Student-t, contaminated normal, and slash distributions. The proposed models provide greater flexibility to analyze a lot of practical data, which often show heavy-tail and may not satisfy the normal assumption. As to the Bayesian analysis, we specify some prior distributions for the unknown parameters in the variance change-point models with the SMN distributions. Due to the complexity of the joint posterior distribution, we propose an efficient Gibbs-type with Metropolis- Hastings sampling algorithm for posterior Bayesian inference. Thereafter, following the idea of [1], we consider the problems of the single and multiple change-point detections. The performance of the proposed procedures is illustrated and analyzed by simulation studies. A real application to the closing price data of U.S. stock market has been analyzed for illustrative purposes.

Subcellular distribution of the insulin-like growth factor (IGF) binding proteins (IGFBPs) 2 and 3 in articular chondrocytes

The insulin-like growth factor (IGF) is a major anabolic regulator in articular cartilage. The IGF-binding proteins (IGFBPs) are increased during osteoarthritis (OA), but the function of the later proteins remains unknown. In general, the IGFBPs are pluripotential effectors capable of IGF regulation and of acting on their own to control key cell functions, including survival and proliferation. The independent functions are often associated with their cell location, and therefore this study explores the distribution of IGFBP-2 and IGFBP-3 in articular chondrocytes. Immunohistochemistry was used to localize IGFBP-2 in normal human articular cartilage. Bovine chondrocytes were used for subcellular fractionation (hypotonic cell lysis) under nonreducing conditions and nuclear purification (centrifugation on sucrose cushions). Cell fraction markers and IGFBPs were assayed in the subcellular fractions by Western immunoblot. The IHC results showed association of IGFBP-2 with chondrocytes, but not with the nuclei. Subcellular fractionation of isolated chondrocytes yielded intact nuclei as assessed at the light microscopic level; the nuclear marker histone H1 was exclusively associated with this fraction. More than 90% of the cytoplasmic marker GAPDH and all the detectable IGFBP-2 were in the cytoplasmic fraction. Immunoreactive IGFBP-3 was found in the cytoplasmic and peri-nuclear/nuclear fractions. Chondrocytes contain intracellular IGFBP-2 and IGFBP-3 but only IGFBP-3 is associated with nuclei. This suggests the hypothesis that the actions of these IGFBPs in articular cartilage extend beyond the classic modulation of IGF receptor action.

Bone mineral density at distal tibia using dual-energy X-ray absorptiometry in normal women and in patients with vertebral osteoporosis or primary hyperparathyroidism

To assess the effect of age and disease on mineral distribution at the distal third of the tibia, bone mineral content (BMC) and bone mineral density (BMD) were measured at lumbar spine (spine), femoral neck (neck), and diaphysis (Dia) and distal epiphysis (Epi) of the tibia in 89 healthy control women of different age groups (20-29, n = 12; 30-39, n = 11; 40-44, n = 12; 45-49, n = 12; 50-54, n = 12; 55-59, n = 10; 60-69, n = 11; 70-79, n = 9), in 25 women with untreated vertebral osteoporosis (VOP), and in 19 women with primary hyperparathyroidism (PHPT) using dual-energy x-ray absorptiometry (DXA; Hologic QDR 1000 and standard spine software). A soft tissue simulator was used to compensate for heterogeneity of soft tissue thickness around the leg. Tibia was scanned over a length of 130 mm from the ankle joint, fibula being excluded from analysis. For BMC and BMD, 10 sections 13 mm each were analyzed separately and then pooled to define the epiphysis (Epi 13-52 mm) and diaphysis area (Dia 91-130 mm). Precision after repositioning was 1.9 and 2.1% for Epi and Dia, respectively. In the control group, at any site there was no significant difference between age groups 20-29 and 30-39, which thus were pooled to define the peak bone mass (PBM).(ABSTRACT TRUNCATED AT 250 WORDS)

Influence of left ventricular relaxation on the pressure half time of aortic regurgitation

BACKGROUND The severity of aortic regurgitation can be estimated using pressure half time (PHT) of the aortic regurgitation flow velocity, but the correlation between regurgitant fraction and PHT is weak. AIM To test the hypothesis that the association between PHT and regurgitant fraction is substantially influenced by left ventricular relaxation. METHODS In 63 patients with aortic regurgitation, subdivided into a group without (n = 22) and a group with (n = 41) left ventricular hypertrophy, regurgitant fraction was calculated using the difference between right and left ventricular cardiac outputs. Left ventricular relaxation was assessed using the early to late diastolic Doppler tissue velocity ratio of the mitral annulus (E/ADTI), the E/A ratio of mitral inflow (E/AM), and the E deceleration time (E-DT). Left ventricular hypertrophy was assessed using the M mode derived left ventricular mass index. RESULTS The overall correlation between regurgitant fraction and PHT was weak (r = 0.36, p < 0.005). In patients without left ventricular hypertrophy, there was a significant correlation between regurgitant fraction and PHT (r = 0.62, p < 0.005), but not in patients with left ventricular hypertrophy. In patients with a left ventricular relaxation abnormality (defined as E/ADTI< 1, E/AM< age corrected lower limit, E-DT >/= 220 ms), no associations between regurgitant fraction and PHT were found, whereas in patients without left ventricular relaxation abnormalities, the regurgitant fraction to PHT relations were significant (normal E/AM: r = 0.57, p = 0.02; E-DT< 220 ms: r = 0.50, p < 0.001; E/ADTI < 1: r = 0.57, p = 0.02). CONCLUSIONS Only normal left ventricular relaxation allows a significant decay of PHT with increasing aortic regurgitation severity. In abnormal relaxation, which is usually present in left ventricular hypertrophy, wide variation in prolonged backward left ventricular filling may cause dissociation between the regurgitant fraction and PHT. Thus the PHT method should only be used in the absence of left ventricular relaxation abnormalities.

Synchrotron microbeam radiation therapy induces hypoxia in intracerebral gliosarcoma but not in the normal brain

PURPOSE Synchrotron microbeam radiation therapy (MRT) is an innovative irradiation modality based on spatial fractionation of a high-dose X-ray beam into lattices of microbeams. The increase in lifespan of brain tumor-bearing rats is associated with vascular damage but the physiological consequences of MRT on blood vessels have not been described. In this manuscript, we evaluate the oxygenation changes induced by MRT in an intracerebral 9L gliosarcoma model. METHODS Tissue responses to MRT (two orthogonal arrays (2 × 400Gy)) were studied using magnetic resonance-based measurements of local blood oxygen saturation (MR_SO2) and quantitative immunohistology of RECA-1, Type-IV collagen and GLUT-1, marker of hypoxia. RESULTS In tumors, MR_SO2 decreased by a factor of 2 in tumor between day 8 and day 45 after MRT. This correlated with tumor vascular remodeling, i.e. decrease in vessel density, increases in half-vessel distances (×5) and GLUT-1 immunoreactivity. Conversely, MRT did not change normal brain MR_SO2, although vessel inter-distances increased slightly. CONCLUSION We provide new evidence for the differential effect of MRT on tumor vasculature, an effect that leads to tumor hypoxia. As hypothesized formerly, the vasculature of the normal brain exposed to MRT remains sufficiently perfused to prevent any hypoxia.

Osteopontin (OPN) and neoplastic disease: Circulating levels and tissue distribution

OPN is a secreted phosphate containing protein which is expressed by osteoblasts and a variety of other cells in vivo. Data from in vitro studies has accumulated which relates OPN to cellular transformation. We hypothesize that OPN expression is associated with neoplastic disease in humans as suggested by cell culture models. The overall objective of the current study was to determine the tissue distribution of OPN in human malignancy and to determine whether or not a correlation exists between OPN serum levels and malignancy. At the inception of this project, no study had been made demonstrating the relevance of OPN expression with naturally occurring neoplastic disease in humans. To date, few studies have reported OPN distribution in human neoplasia and are limited by either the number of specimens analyzed or the technique used in analysis. In this dissertation study, OPN was purified from human milk and $\alpha$-OPN antiserum developed and characterized. Following antibody development, the distribution and prevalence of OPN in human oral squamous cell carcinoma and human prostate carcinoma was evaluated using immunohistochemical localization. OPN immunolocalization was found in a high percentage of oral epithelial dysplasia and oral squamous cell carcinoma in humans. One oral squamous cell carcinoma cells line, UMSCC-1, was found to express OPN mRNA using Northern blotting. OPN localized to a high percentage of primary prostate adenocarcinomas. OPN localized to 52% of androgen dependent cases and 100% of androgen independent cases. Androgen dependent cell lines such as LNCap and NbE showed minimal OPN mRNA expression while the androgen independent lines C4-2 and PC3 produced ample OPN mRNA. An OPN sandwich assay was developed and used to determine the serum level of OPN in normal males, patients with BPH (benign prostate hypertrophy), and patients with prostate carcinoma. No statistically significant difference was found in OPN serum levels among the three groups. However, a trend of increasing OPN in the serum was noted in patients with BPH and prostate cancer. ^

Statistical characterization and development of summary measures of non-normal distributions of nuclear morphometry data from prostate cancer cells for use as prognostic indicators

Nuclear morphometry (NM) uses image analysis to measure features of the cell nucleus which are classified as: bulk properties, shape or form, and DNA distribution. Studies have used these measurements as diagnostic and prognostic indicators of disease with inconclusive results. The distributional properties of these variables have not been systematically investigated although much of the medical data exhibit nonnormal distributions. Measurements are done on several hundred cells per patient so summary measurements reflecting the underlying distribution are needed.^ Distributional characteristics of 34 NM variables from prostate cancer cells were investigated using graphical and analytical techniques. Cells per sample ranged from 52 to 458. A small sample of patients with benign prostatic hyperplasia (BPH), representing non-cancer cells, was used for general comparison with the cancer cells.^ Data transformations such as log, square root and 1/x did not yield normality as measured by the Shapiro-Wilks test for normality. A modulus transformation, used for distributions having abnormal kurtosis values, also did not produce normality.^ Kernel density histograms of the 34 variables exhibited non-normality and 18 variables also exhibited bimodality. A bimodality coefficient was calculated and 3 variables: DNA concentration, shape and elongation, showed the strongest evidence of bimodality and were studied further.^ Two analytical approaches were used to obtain a summary measure for each variable for each patient: cluster analysis to determine significant clusters and a mixture model analysis using a two component model having a Gaussian distribution with equal variances. The mixture component parameters were used to bootstrap the log likelihood ratio to determine the significant number of components, 1 or 2. These summary measures were used as predictors of disease severity in several proportional odds logistic regression models. The disease severity scale had 5 levels and was constructed of 3 components: extracapsulary penetration (ECP), lymph node involvement (LN+) and seminal vesicle involvement (SV+) which represent surrogate measures of prognosis. The summary measures were not strong predictors of disease severity. There was some indication from the mixture model results that there were changes in mean levels and proportions of the components in the lower severity levels. ^

Impact of synchronous metastasis distribution on cancer specific survival in renal cell carcinoma after radical nephrectomy with tumor thrombectomy

PURPOSE Metastatic renal cell carcinoma can be clinically diverse in terms of the pattern of metastatic disease and response to treatment. We studied the impact of metastasis and location on cancer specific survival. MATERIALS AND METHODS The records of 2,017 patients with renal cell cancer and tumor thrombus who underwent radical nephrectomy and tumor thrombectomy from 1971 to 2012 at 22 centers in the United States and Europe were analyzed. Number and location of synchronous metastases were compared with respect to patient cancer specific survival. Multivariable Cox regression models were used to quantify the impact of covariates. RESULTS Lymph node metastasis (155) or distant metastasis (725) was present in 880 (44%) patients. Of the patients with distant disease 385 (53%) had an isolated metastasis. The 5-year cancer specific survival was 51.3% (95% CI 48.6-53.9) for the entire group. On univariable analysis patients with isolated lymph node metastasis had a significantly worse cancer specific survival than those with a solitary distant metastasis. The location of distant metastasis did not have any significant effect on cancer specific survival. On multivariable analysis the presence of lymph node metastasis, isolated distant metastasis and multiple distant metastases were independently associated with cancer specific survival. Moreover higher tumor thrombus level, papillary histology and the use of postoperative systemic therapy were independently associated with worse cancer specific survival. CONCLUSIONS In our multi-institutional series of patients with renal cell cancer who underwent radical nephrectomy and tumor thrombectomy, almost half of the patients had synchronous lymph node or distant organ metastasis. Survival was superior in patients with solitary distant metastasis compared to isolated lymph node disease.

Two-dimensional and three-dimensional oxygen mapping by 3He-MRI validation in a lung phantom

The aim of this study was to validate oxygen-sensitive 3He-MRI in noninvasive determination of the regional, two- and three-dimensional distribution of oxygen partial pressure. In a gas-filled elastic silicon ventilation bag used as a lung phantom, oxygen sensitive two- and three-dimensional 3He-MRI measurements were performed at different oxygen concentrations which had been equilibrated in a range of normal and pathologic values. The oxygen partial pressure distribution was determined from 3He-MRI using newly developed software allowing for mapping of oxygen partial pressure. The reference bulk oxygen partial pressure inside the phantom was measured by conventional respiratory gas analysis. In two-dimensional measurements, image-based and gas-analysis results correlated with r=0.98; in three-dimensional measurements the between-methods correlation coefficient was r=0.89. The signal-to-noise ratio of three-dimensional measurements was about half of that of two-dimensional measurements and became critical (below 3) in some data sets. Oxygen-sensitive 3He-MRI allows for noninvasive determination of the two- and three-dimensional distribution of oxygen partial pressure in gas-filled airspaces.

Spatial Distribution of Stem Cell-Like Keratinocytes in Dissected Compound Hair Follicles of the Dog

Hair cycle disturbances are common in dogs and comparable to some alopecic disorders in humans. A normal hair cycle is maintained by follicular stem cells which are predominately found in an area known as the bulge. Due to similar morphological characteristics of the bulge area in humans and dogs, the shared particularity of compound hair follicles as well as similarities in follicular biomarker expression, the dog is a promising model to study human hair cycle and stem cell disorders. To gain insight into the spatial distribution of follicular keratinocytes with stem cell potential in canine compound follicles, we microdissected hair follicles in anagen and telogen from skin samples of freshly euthanized dogs. The keratinocytes isolated from different locations were investigated for their colony forming efficiency, growth and differentiation potential as well as clonal growth. Our results indicate that i) compound and single hair follicles exhibit a comparable spatial distribution pattern with respect to cells with high growth potential and stem cell-like characteristics, ii) the lower isthmus (comprising the bulge) harbors most cells with high growth potential in both, the anagen and the telogen hair cycle stage, iii) unlike in other species, colonies with highest growth potential are rather small with an irregular perimeter and iv) the keratinocytes derived from the bulbar region exhibit characteristics of actively dividing transit amplifying cells. Our results now provide the basis to conduct comparative studies of normal dogs and those with hair cycle disorders with the possibility to extend relevant findings to human patients.

Body fat distribution and serum lipids in children and adolescents

Previous studies of normal children have linked body fat but not body fat distribution (BFD), to higher blood pressures, lipids, and insulin resistance (Berenson et al., 1988) BFD is a well-established risk factor for cardiovascular disease in adults (Björntorp, 1988). This study investigates the relation of BFD and serum lipids at baseline in children from Project HeartBeat!, a study of the growth and development of cardiovascular risk factors in 678 children in three cohorts measured initially at ages 8, 11, and 14 years. Initially, two of four indices of BFD were significantly related to the lipids: ratio of upper to lower body skinfolds (ln US:LS) and conicity (C Index). A factor analysis reduced the information in the serum lipids to two vectors: (1) total cholesterol + LDL-cholesterol and (2) HDL-cholesterol − triglycerides, which together accounted for 85% of the lipid variation. Using each serum lipid and vector as separate dependent variables, linear and quadratic regression models were constructed to examine the predictive ability of the two BFD variables, controlling for total body fat, gender, ethnicity (Black, non-Black) and maturation. Linear models provided an acceptable fit. Percent body fat (%BF) was a significant predictor in each and every lipid model, independent of age, maturation, or ethnicity (p ≤ 0.05). No BFD variable entered the equation for total or LDL-cholesterol, although there was a significant maturity by BFD interaction for LDL (ln US:LS was a significant predictor in more mature individuals). Both %BF and BFD (by way of Conicity) were significant predictors of HDL-cholesterol and triglycerides (p ≤ 0.01). All models were statistically significant at a high level (p ≤ 0.01), but adjusted R 2's for all models were low (0.05–0.15). Body fat distribution is a significant predictor of lipids in normal children, but secondarily to %BF, and for LDL-cholesterol in particular, the relation is dependent on maturity status. ^