Power Analysis in Applied Linear Regression for Cell Type-Specific Differential Expression Detection

The goal of many human disease-oriented studies is to detect molecular mechanisms different between healthy controls and patients. Yet, commonly used gene expression measurements from any tissues suffer from variability of cell composition. This variability hinders the detection of differentially expressed genes and is often ignored. However, this variability may actually be advantageous, as heterogeneous gene expression measurements coupled with cell counts may provide deeper insights into the gene expression differences on the cell type-specific level. Published computational methods use linear regression to estimate cell type-specific differential expression. Yet, they do not consider many artifacts hidden in high-dimensional gene expression data that may negatively affect the performance of linear regression. In this dissertation we specifically address the parameter space involved in the most rigorous use of linear regression to estimate cell type-specific differential expression and report under which conditions significant detection is probable. We define parameters affecting the sensitivity of cell type-specific differential expression estimation as follows: sample size, cell type-specific proportion variability, mean squared error (spread of observations around linear regression line), conditioning of the cell proportions predictor matrix, and the size of actual cell type-specific differential expression. Each parameter, with the exception of cell type-specific differential expression (effect size), affects the variability of cell type-specific differential expression estimates. We have developed a power-analysis approach to cell type by cell type and genomic site by site differential expression detection which relies upon Welch’s two-sample t-test and factors in differences in cell type-specific expression estimate variability and reduces false discovery. To this end we have published an R package, LRCDE, available in GitHub (http://www.github.com/ERGlass/lrcde.dev) which outputs observed statistics of cell type-specific differential expression, including two-sample t- statistic, t-statistic p-value, and power calculated from two-sample t-statistic on a genomic site- by-site basis

B-Cell and Monocyte Contribution to Systemic Lupus Erythematosus Identified by Cell-Type-Specific Differential Expression Analysis in RNA-Seq Data

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by complex interplay among immune cell types. SLE activity is experimentally assessed by several blood tests, including gene expression profiling of heterogeneous populations of cells in peripheral blood. To better understand the contribution of different cell types in SLE pathogenesis, we applied the two methods in cell-type-specific differential expression analysis, csSAM and DSection, to identify cell-type-specific gene expression differences in heterogeneous gene expression measures obtained using RNA-seq technology. We identified B-cell-, monocyte-, and neutrophil-specific gene expression differences. Immunoglobulin-coding gene expression was altered in B-cells, while a ribosomal signature was prominent in monocytes. On the contrary, genes differentially expressed in the heterogeneous mixture of cells did not show any functional enrichment. Our results identify antigen binding and structural constituents of ribosomes as functions altered by B-cell- and monocyte-specific gene expression differences, respectively. Finally, these results position both csSAM and DSection methods as viable techniques for celltype-specific differential expression analysis, which may help uncover pathogenic, cell-type-specific processes in SLE

Detrimental effects of duplicate reads and low complexity regions on RNA- and ChIP-seq data

Background Adapter trimming and removal of duplicate reads are common practices in next-generation sequencing pipelines. Sequencing reads ambiguously mapped to repetitive and low complexity regions can also be problematic for accurate assessment of the biological signal, yet their impact on sequencing data has not received much attention. We investigate how trimming the adapters, removing duplicates, and filtering out reads overlapping low complexity regions influence the significance of biological signal in RNA- and ChIP-seq experiments. Methods We assessed the effect of data processing steps on the alignment statistics and the functional enrichment analysis results of RNA- and ChIP-seq data. We compared differentially processed RNA-seq data with matching microarray data on the same patient samples to determine whether changes in pre-processing improved correlation between the two. We have developed a simple tool to remove low complexity regions, RepeatSoaker, available at https://github.com/mdozmorov/RepeatSoaker, and tested its effect on the alignment statistics and the results of the enrichment analyses. Results Both adapter trimming and duplicate removal moderately improved the strength of biological signals in RNA-seq and ChIP-seq data. Aggressive filtering of reads overlapping with low complexity regions, as defined by RepeatMasker, further improved the strength of biological signals, and the correlation between RNA-seq and microarray gene expression data. Conclusions Adapter trimming and duplicates removal, coupled with filtering out reads overlapping low complexity regions, is shown to increase the quality and reliability of detecting biological signals in RNA-seq and ChIP-seq data

Polymorphic markers related to a single Tcrb-V6 gene segment

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46749/1/251_2004_Article_BF00230507.pd

A hybrid approach to determining cornea mechanical properties in vivo using a combination of nano-indentation and inverse finite element analysis

This article was published in the serial, Journal of the Mechanical Behavior of Biomedical Materials [© Elsevier]. The definitive version is available at: http://dx.doi.org/10.1016/j.jmbbm.2013.05.016An analysis of the material properties of porcine corneas has been performed. A simple stress relaxation test was performed to determine the viscoelastic properties and a rheological model was built based on the Generalized Maxwell (GM) approach. A validation experiment using nano-indentation showed that an isotropic GM model was insufficient for describing the corneal material behaviour when exposed to a complex stress state. A new technique was proposed for determining the properties, using a combination of nano-indentation experiment, an isotropic and orthotropic GM model and inverse finite element method. The good agreement using this method suggests that this is a promising technique for measuring material properties in vivo and further work should focus on the reliability of the approach in practice

Assessing circadian rhythms in propofol PK and PD during prolonged infusion in ICU patients

This study evaluates possible circadian rhythms during prolonged propofol infusion in patients in the intensive care unit. Eleven patients were sedated with a constant propofol infusion. The blood samples for the propofol assay were collected every hour during the second day, the third day, and after the termination of the propofol infusion. Values of electroencephalographic bispectral index (BIS), arterial blood pressure, heart rate, blood oxygen saturation and body temperature were recorded every hour at the blood collection time points. A two-compartment model was used to describe propofol pharmacokinetics. Typical values of the central and peripheral volume of distribution and inter-compartmental clearance were VC = 27.7 l, VT = 801 l, and CLD = 2.73 l/min. The systolic blood pressure (SBP) was found to influence the propofol metabolic clearance according to Cl (l/min) = 2.65·(1 − 0.00714·(SBP − 135)). There was no significant circadian rhythm detected with respect to propofol pharmacokinetics. The BIS score was assessed as a direct effect model with EC50 equal 1.98 mg/l. There was no significant circadian rhythm detected within the BIS scores. We concluded that the light–dark cycle did not influence propofol pharmacokinetics and pharmacodynamics in intensive care units patients. The lack of night–day differences was also noted for systolic blood pressure, diastolic blood pressure and blood oxygenation. Circadian rhythms were detected for heart rate and body temperature, however they were severely disturbed from the pattern of healthy patients

Headache impairs attentional performance: a conceptual replication and extension

Pain is thought to capture our attention. A consequence is that our performance on other tasks may suffer. Research has supported this, showing that pain disrupts our ability to perform various attention tasks. However, the specific nature of the effect of pain on attention is inconsistent, possibly due to different studies investigating different types of pain. Few studies seek to replicate basic findings. Here, we conceptually replicated and extended the headache study by Moore, Keogh & Eccleston (2013), by including two additional attention tasks, a broader sample, and measures of affect and pain cognitions. Participants performed five complex attention tasks and a choice reaction time task with and without a naturally-occurring headache. Headache slowed reaction times to four of the five complex tasks, and this could be attributed to a slower basic processing speed as measured by the choice reaction time task. Our findings differ from those of Moore et al’s headache study, suggesting that the effect of pain on attention is dynamic, even within a given type of pain. While there is growing evidence that pain does disrupt attention, we cannot yet predict the specific nature of disruption in any given case

Adjunctive rifampicin for Staphylococcus aureus bacteraemia (ARREST): a multicentre, randomised, double-blind, placebo-controlled trial.

BACKGROUND: Staphylococcus aureus bacteraemia is a common cause of severe community-acquired and hospital-acquired infection worldwide. We tested the hypothesis that adjunctive rifampicin would reduce bacteriologically confirmed treatment failure or disease recurrence, or death, by enhancing early S aureus killing, sterilising infected foci and blood faster, and reducing risks of dissemination and metastatic infection. METHODS: In this multicentre, randomised, double-blind, placebo-controlled trial, adults (≥18 years) with S aureus bacteraemia who had received ≤96 h of active antibiotic therapy were recruited from 29 UK hospitals. Patients were randomly assigned (1:1) via a computer-generated sequential randomisation list to receive 2 weeks of adjunctive rifampicin (600 mg or 900 mg per day according to weight, oral or intravenous) versus identical placebo, together with standard antibiotic therapy. Randomisation was stratified by centre. Patients, investigators, and those caring for the patients were masked to group allocation. The primary outcome was time to bacteriologically confirmed treatment failure or disease recurrence, or death (all-cause), from randomisation to 12 weeks, adjudicated by an independent review committee masked to the treatment. Analysis was intention to treat. This trial was registered, number ISRCTN37666216, and is closed to new participants. FINDINGS: Between Dec 10, 2012, and Oct 25, 2016, 758 eligible participants were randomly assigned: 370 to rifampicin and 388 to placebo. 485 (64%) participants had community-acquired S aureus infections, and 132 (17%) had nosocomial S aureus infections. 47 (6%) had meticillin-resistant infections. 301 (40%) participants had an initial deep infection focus. Standard antibiotics were given for 29 (IQR 18-45) days; 619 (82%) participants received flucloxacillin. By week 12, 62 (17%) of participants who received rifampicin versus 71 (18%) who received placebo experienced treatment failure or disease recurrence, or died (absolute risk difference -1·4%, 95% CI -7·0 to 4·3; hazard ratio 0·96, 0·68-1·35, p=0·81). From randomisation to 12 weeks, no evidence of differences in serious (p=0·17) or grade 3-4 (p=0·36) adverse events were observed; however, 63 (17%) participants in the rifampicin group versus 39 (10%) in the placebo group had antibiotic or trial drug-modifying adverse events (p=0·004), and 24 (6%) versus six (2%) had drug interactions (p=0·0005). INTERPRETATION: Adjunctive rifampicin provided no overall benefit over standard antibiotic therapy in adults with S aureus bacteraemia. FUNDING: UK National Institute for Health Research Health Technology Assessment