A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

Studying Effects of Primary Care Physicians and Patients on the Trade-Off Between Charges for Primary Care and Specialty Care Using a Hierarchical Multivariate Two-Part Model

Objective. To examine effects of primary care physicians (PCPs) and patients on the association between charges for primary care and specialty care in a point-of-service (POS) health plan. Data Source. Claims from 1996 for 3,308 adult male POS plan members, each of whom was assigned to one of the 50 family practitioner-PCPs with the largest POS plan member-loads. Study Design. A hierarchical multivariate two-part model was fitted using a Gibbs sampler to estimate PCPs' effects on patients' annual charges for two types of services, primary care and specialty care, the associations among PCPs' effects, and within-patient associations between charges for the two services. Adjusted Clinical Groups (ACGs) were used to adjust for case-mix. Principal Findings. PCPs with higher case-mix adjusted rates of specialist use were less likely to see their patients at least once during the year (estimated correlation: –.40; 95% CI: –.71, –.008) and provided fewer services to patients that they saw (estimated correlation: –.53; 95% CI: –.77, –.21). Ten of 11 PCPs whose case-mix adjusted effects on primary care charges were significantly less than or greater than zero (p < .05) had estimated, case-mix adjusted effects on specialty care charges that were of opposite sign (but not significantly different than zero). After adjustment for ACG and PCP effects, the within-patient, estimated odds ratio for any use of primary care given any use of specialty care was .57 (95% CI: .45, .73). Conclusions. PCPs and patients contributed independently to a trade-off between utilization of primary care and specialty care. The trade-off appeared to partially offset significant differences in the amount of care provided by PCPs. These findings were possible because we employed a hierarchical multivariate model rather than separate univariate models.

RANDOM EFFECTS MODELS IN A META-ANALYSIS OF THE ACCURACY OF DIAGNOSTIC TESTS WITHIN A GOLD STANDARD IN THE PRESENCE OF MISSING DATA

In evaluating the accuracy of diagnosis tests, it is common to apply two imperfect tests jointly or sequentially to a study population. In a recent meta-analysis of the accuracy of microsatellite instability testing (MSI) and traditional mutation analysis (MUT) in predicting germline mutations of the mismatch repair (MMR) genes, a Bayesian approach (Chen, Watson, and Parmigiani 2005) was proposed to handle missing data resulting from partial testing and the lack of a gold standard. In this paper, we demonstrate an improved estimation of the sensitivities and specificities of MSI and MUT by using a nonlinear mixed model and a Bayesian hierarchical model, both of which account for the heterogeneity across studies through study-specific random effects. The methods can be used to estimate the accuracy of two imperfect diagnostic tests in other meta-analyses when the prevalence of disease, the sensitivities and/or the specificities of diagnostic tests are heterogeneous among studies. Furthermore, simulation studies have demonstrated the importance of carefully selecting appropriate random effects on the estimation of diagnostic accuracy measurements in this scenario.

Vascular remodeling and antitumoral effects of mTOR inhibition in a rat model of hepatocellular carcinoma

BACKGROUND/AIMS: Hepatocellular carcinoma (HCC) is amenable to only few treatments. Inhibitors of the kinase mTOR are a new class of immunosuppressors already in use after liver transplantation. Their antiproliferative and antiangiogenic properties suggest that these drugs could be considered to treat HCC. We investigated the antitumoral effects of mTOR inhibition in a HCC model. METHODS: Hepatoma cells were implanted into livers of syngeneic rats. Animals were treated with the mTOR inhibitor sirolimus for 4 weeks. Tumor growth was monitored by MR imaging. Antiangiogenic effects were assessed in vivo by microvessel density and corrosion casts and in vitro by cell proliferation, tube formation and aortic ring assays. RESULTS: Treated rats had significantly longer survival and developed smaller tumors, fewer extrahepatic metastases and less ascites than controls. Sirolimus decreased intratumoral microvessel density resulting in extensive necrosis. Endothelial cell proliferation was inhibited at lower drug concentrations than hepatoma cells. Tube formation and vascular sprouting of aortic rings were significantly impaired by mTOR inhibition. Casts revealed that in tumors treated with sirolimus vascular sprouting was absent, whereas intussusception was observed. CONCLUSIONS: mTOR inhibition significantly reduces HCC growth and improves survival primarily via antiangiogenic effects. Inhibitors of mTOR may have a role in HCC treatment.

Effects of an aging vascular model on healthy and diseased hearts

The vitamin D(3) and nicotine (VDN) model is a model of isolated systolic hypertension (ISH) due to arterial calcification raising arterial stiffness and vascular impedance similar to an aged and stiffened arterial tree. We therefore analyzed the impact of this aging model on normal and diseased hearts with myocardial infarction (MI). Wistar rats were treated with VDN (n = 9), subjected to MI by coronary ligation (n = 10), or subjected to a combination of both MI and VDN treatment (VDN/MI, n = 14). A sham-treated group served as control (Ctrl, n = 10). Transthoracic echocardiography was performed every 2 wk, whereas invasive indexes were obtained at week 8 before death. Calcium, collagen, and protein contents were measured in the heart and the aorta. Systolic blood pressure, pulse pressure, thoracic aortic calcium, and end-systolic elastance as an index of myocardial contractility were highest in the aging model group compared with MI and Ctrl groups (P(VDN) < 0.05, 2-way ANOVA). Left ventricular wall stress and brain natriuretic peptide (P(VDNxMI) = not significant) were highest, while ejection fraction, stroke volume, and cardiac output were lowest in the combined group versus all other groups (P(VDNxMI) < 0.05). The combination of ISH due to this aging model and MI demonstrates significant alterations in cardiac function. This model mimics several clinical phenomena of cardiovascular aging and may thus serve to further study novel therapies.

Drug-induced respiratory depression: an integrated model of drug effects on the hypercapnic and hypoxic drive

Drug-induced respiratory depression is a common side effect of the agents used in anesthesia practice to provide analgesia and sedation. Depression of the ventilatory drive in the spontaneously breathing patient can lead to severe cardiorespiratory events and it is considered a primary cause of morbidity. Reliable predictions of respiratory inhibition in the clinical setting would therefore provide a valuable means to improve the safety of drug delivery. Although multiple studies investigated the regulation of breathing in man both in the presence and absence of ventilatory depressant drugs, a unified description of respiratory pharmacodynamics is not available. This study proposes a mathematical model of human metabolism and cardiorespiratory regulation integrating several isolated physiological and pharmacological aspects of acute drug-induced ventilatory depression into a single theoretical framework. The description of respiratory regulation has a parsimonious yet comprehensive structure with substantial predictive capability. Simulations relative to the synergistic interaction of the hypercarbic and hypoxic respiratory drive and the global effect of drugs on the control of breathing are in good agreement with published experimental data. Besides providing clinically relevant predictions of respiratory depression, the model can also serve as a test bed to investigate issues of drug tolerability and dose finding/control under non-steady-state conditions.

Exogenous surfactant application in a rat lung ischemia reperfusion injury model: effects on edema formation and alveolar type II cells

BACKGROUND: Prophylactic exogenous surfactant therapy is a promising way to attenuate the ischemia and reperfusion (I/R) injury associated with lung transplantation and thereby to decrease the clinical occurrence of acute lung injury and acute respiratory distress syndrome. However, there is little information on the mode by which exogenous surfactant attenuates I/R injury of the lung. We hypothesized that exogenous surfactant may act by limiting pulmonary edema formation and by enhancing alveolar type II cell and lamellar body preservation. Therefore, we investigated the effect of exogenous surfactant therapy on the formation of pulmonary edema in different lung compartments and on the ultrastructure of the surfactant producing alveolar epithelial type II cells. METHODS: Rats were randomly assigned to a control, Celsior (CE) or Celsior + surfactant (CE+S) group (n = 5 each). In both Celsior groups, the lungs were flush-perfused with Celsior and subsequently exposed to 4 h of extracorporeal ischemia at 4 degrees C and 50 min of reperfusion at 37 degrees C. The CE+S group received an intratracheal bolus of a modified natural bovine surfactant at a dosage of 50 mg/kg body weight before flush perfusion. After reperfusion (Celsior groups) or immediately after sacrifice (Control), the lungs were fixed by vascular perfusion and processed for light and electron microscopy. Stereology was used to quantify edematous changes as well as alterations of the alveolar epithelial type II cells. RESULTS: Surfactant treatment decreased the intraalveolar edema formation (mean (coefficient of variation): CE: 160 mm3 (0.61) vs. CE+S: 4 mm3 (0.75); p < 0.05) and the development of atelectases (CE: 342 mm3 (0.90) vs. CE+S: 0 mm3; p < 0.05) but led to a higher degree of peribronchovascular edema (CE: 89 mm3 (0.39) vs. CE+S: 268 mm3 (0.43); p < 0.05). Alveolar type II cells were similarly swollen in CE (423 microm3(0.10)) and CE+S (481 microm3(0.10)) compared with controls (323 microm3(0.07); p < 0.05 vs. CE and CE+S). The number of lamellar bodies was increased and the mean lamellar body volume was decreased in both CE groups compared with the control group (p < 0.05). CONCLUSION: Intratracheal surfactant application before I/R significantly reduces the intraalveolar edema formation and development of atelectases but leads to an increased development of peribronchovascular edema. Morphological changes of alveolar type II cells due to I/R are not affected by surfactant treatment. The beneficial effects of exogenous surfactant therapy are related to the intraalveolar activity of the exogenous surfactant.

Effects of EDP-420 on penicillin-resistant and quinolone- and penicillin-resistant pneumococci in the rabbit meningitis model

OBJECTIVES: To test the efficacy of EDP-420, a new ketolide, in experimental pneumococcal meningitis and to determine its penetration into the CSF. METHODS: The experimental rabbit model was used in this study and EDP-420 was tested against a penicillin-resistant and a penicillin- and quinolone-resistant mutant. EDP-420 was also tested against both strains in time-killing assays over 8 h in vitro. RESULTS: In experimental meningitis, EDP-420 produced a bactericidal activity comparable to the standard regimen based on a combination of vancomycin with ceftriaxone against a penicillin-resistant Streptococcus pneumoniae and a penicillin- and quinolone-resistant S. pneumoniae isolate. The penetration of EDP-420 into inflamed meninges was 38% after an i.v. injection of 10 mg/kg. The bactericidal activity of EDP-420 was also confirmed in in vitro time-killing assays. CONCLUSIONS: EDP-420 is an efficacious alternative treatment in pneumococcal meningitis, especially when resistant strains are suspected.

Cytotoxic effects of camptothecin and cisplatin combined with tumor necrosis factor-related apoptosis-inducing ligand (Apo2L/TRAIL) in a model of primary culture of non-small cell lung cancer

The cytokine tumor-necrosis factor-related apoptosis-inducing ligand (Apo2L/TRAIL) has been shown to preferentially induce apoptosis in cancer cells. A previous study of our group demonstrated that non-small cell lung cancer cell lines can be sensitized to Apo2L/TRAIL-induced apoptosis by chemotherapeutic agents. The aim of the present study was the evaluation of these results in a model of primary culture of non-small cell lung cancer.

Effects of GDNF pretreatment on function and survival of transplanted fetal ventral mesencephalic cells in the 6-OHDA rat model of Parkinson's disease

Transplantation of fetal dopaminergic (DA) neurons offers an experimental therapy for Parkinson's disease (PD). The low availability and the poor survival and integration of transplanted cells in the host brain are major obstacles in this approach. Glial cell line-derived neurotrophic factor (GDNF) is a potent neurotrophic factor with growth- and survival-promoting capabilities for developing DA neurons. In the present study, we examined whether pretreatment of ventral mesencephalic (VM) free-floating roller tube (FFRT) cultures with GDNF would improve graft survival and function. For that purpose organotypic cultures of E14 rat VM were grown for 2, 4 or 8 days in the absence (control) or presence of GDNF [10 ng/ml] and transplanted into the striatum of 6-hydroxydopamine-lesioned rats. While all groups of rats showed a significant reduction in d-amphetamine-induced rotations at 6 weeks posttransplantation a significantly improved graft function was observed only in the days in vitro (DIV) 4 GDNF pretreated group compared to the control group. In addition, no statistical significant differences between groups were found in the number of surviving tyrosine hydroxylase-immunoreactive (TH-ir) neurons assessed at 9 weeks posttransplantation. However, a tendency for higher TH-ir fiber outgrowth from the transplants in the GDNF pretreated groups as compared to corresponding controls was observed. Furthermore, GDNF pretreatment showed a tendency for a higher number of GIRK2 positive neurons in the grafts. In sum, our findings demonstrate that GDNF pretreatment was not disadvantageous for transplants of embryonic rat VM with the FFRT culture technique but only marginally improved graft survival and function.

Assessment of the Matrix Degenerative Effects of MMP-3, ADAMTS-4 and HTRA1 injected into a bovine Intervertebral Disc Organ Culture Model

Study Design. In vitro study to develop an intervertebral disc degeneration (IDD) organ culture model, using coccygeal bovine intervertebral discs (IVDs) and injection of proteolytic enzymes MMP-3, ADAMTS-4 and HTRA1.Objective. This study aimed to develop an in-vitro model of enzyme-mediated IDD to mimic the clinical outcome in humans for investigation of therapeutic treatment options.Summary of Background Data. Bovine IVDs are comparable to human IVDs in terms of cell composition and biomechanical behavior. Researchers injected papain and trypsin into them to create an IDD model with a degenerated nucleus pulposus (NP) area. They achieved macroscopic cavities as well as a loss of glycosaminoglycans (GAGs). However, none of these enzymes are clinically relevant.Methods. Bovine IVDs were harvested maintaining the endplates. Active forms of MMP-3, ADAMTS-4 and HTRA1 were injected at a dose of 10μg/ml each. Phosphate buffered saline (PBS) was injected as a control. Discs were cultured for 8 days and loaded diurnally (day 1 to day 4 with 0.4 MPa for 16 h) and left under free swelling condition from day 4 to day 8 to avoid expected artifacts due to dehydration of the NP. Outcome parameters included disc height, metabolic cell activity, DNA content, glycosaminoglycan (GAG) content, total collagen content, relative gene expression and histological investigation.Results. The mean metabolic cell activity was significantly lower in the NP area of discs injected with ADAMTS-4 compared to the day 0 control discs. Disc height was decreased following injection with HTRA1, and was significantly correlated with changes in GAG/DNA of the NP tissue. Total collagen content tended to be lower in groups injected with ADAMTS4 and MMP-3.Conclusion. MMP-3, ADAMTS-4 and HTRA1 neither provoked visible matrix degradation nor major shifts in gene expression. However, cell activity was significantly reduced and HTRA1 induced loss of disc height which positively correlated with changes in GAG/DNA content. The use of higher doses of these enzymes or a combination thereof may therefore be necessary to induce disc degeneration

Imaging of pulmonary embolism and t-PA therapy effects using MDCT and liposomal iohexol blood pool agent – preliminary results in a rabbit model

Hypothesis and Objectives PEGylated liposomal blood pool contrast agents maintain contrast enhancement over several hours. This study aimed to evaluate (long-term) imaging of pulmonary arteries, comparing conventional iodinated contrast with a liposomal blood pool contrast agent. Secondly, visualization of the (real-time) therapeutic effects of tissue-Plasminogen Activator (t-PA) on pulmonary embolism (PE) was attempted. Materials and Methods Six rabbits (approximate 4 kg weight) had autologous blood clots injected through the superior vena cava. Imaging was performed using conventional contrast (iohexol, 350 mg I/ml, GE HealthCare, Princeton, NJ) at a dose of 1400 mgI per animal and after wash-out, animals were imaged using an iodinated liposomal blood pool agent (88 mg I/mL, dose 900 mgI/animal). Subsequently, five animals were injected with 2mg t-PA and imaging continued for up to 4 ½ hours. Results Both contrast agents identified PE in the pulmonary trunk and main pulmonary arteries in all rabbits. Liposomal blood pool agent yielded uniform enhancement, which remained relatively constant throughout the experiments. Conventional agents exhibited non uniform opacification and rapid clearance post injection. Three out of six rabbits had mistimed bolus injections, requiring repeat injections. Following t-PA, Pulmonary embolus volume (central to segmental) decreased in four of five treated rabbits (range 10–57%, mean 42%). One animal showed no response to t-PA. Conclusions Liposomal blood pool agents effectively identified acute PE without need for re-injection. PE resolution following t-PA was quantifiable over several hours. Blood pool agents offer the potential for repeated imaging procedures without need for repeated (nephrotoxic) contrast injections

Imaging of pulmonary embolism and t-PA therapy effects using MDCT and liposomal iohexol blood pool agent: preliminary results in a rabbit model.

RATIONALE AND OBJECTIVES: Polyethylene glycol-coated liposomal blood pool contrast agents maintain contrast enhancement over several hours. This study aimed to evaluate (long-term) imaging of pulmonary arteries, comparing conventional iodinated contrast with a liposomal blood pool contrast agent. Also, visualization of the (real-time) therapeutic effects of tissue plasminogen activator (t-PA) on pulmonary embolism (PE) was attempted. MATERIALS AND METHODS: Six rabbits (weight approximately 4 kg) had autologous blood clots injected through the superior vena cava. Imaging was performed using conventional contrast (iohexol, 350 mg I/ml; GE HealthCare, Princeton, NJ) at a dose of 1400 mg I per animal, and after wash-out, animals were imaged using an iodinated liposomal blood pool agent (88 mg I/mL, dose 900 mg I/animal). Subsequently, five animals were injected with 2 mg of t-PA and imaging continued for up to 4(1/2) hours. RESULTS: Both contrast agents identified PE in the pulmonary trunk and main pulmonary arteries in all rabbits. Liposomal blood pool agent yielded uniform enhancement, which remained relatively constant throughout the experiments. Conventional agents exhibited nonuniform opacification and rapid clearance postinjection. Three of six rabbits had mistimed bolus injections, requiring repeat injections. Following t-PA, pulmonary embolus volume (central to segmental) decreased in four of five treated rabbits (range 10-57%, mean 42%). One animal showed no response to t-PA. CONCLUSIONS: Liposomal blood pool agents effectively identified acute PE without need for reinjection. PE resolution following t-PA was quantifiable over several hours. Blood pool agents offer the potential for repeated imaging procedures without need for repeated (nephrotoxic) contrast injections.