Correcting the Mean-Variance Dependency for Differential Variability Testing Using Single-Cell RNA Sequencing Data.

Cell-to-cell transcriptional variability in otherwise homogeneous cell populations plays an important role in tissue function and development. Single-cell RNA sequencing can characterize this variability in a transcriptome-wide manner. However, technical variation and the confounding between variability and mean expression estimates hinder meaningful comparison of expression variability between cell populations. To address this problem, we introduce an analysis approach that extends the BASiCS statistical framework to derive a residual measure of variability that is not confounded by mean expression. This includes a robust procedure for quantifying technical noise in experiments where technical spike-in molecules are not available. We illustrate how our method provides biological insight into the dynamics of cell-to-cell expression variability, highlighting a synchronization of biosynthetic machinery components in immune cells upon activation. In contrast to the uniform up-regulation of the biosynthetic machinery, CD4+ T cells show heterogeneous up-regulation of immune-related and lineage-defining genes during activation and differentiation.NE was funded by the European Molecular Biology Laboratory (EMBL) international PhD programme. ACR was funded by the MRC Skills Development Fellowship (MR/P014178/1). SR was funded by MRC grant MC_UP_0801/1. JCM was funded by core support of Cancer Research UK and EMBL. CAV was funded by The Alan Turing Institute, EPSRC grant EP/N510129/1

Detection and removal of barcode swapping in single-cell RNA-seq data.

Barcode swapping results in the mislabelling of sequencing reads between multiplexed samples on patterned flow-cell Illumina sequencing machines. This may compromise the validity of numerous genomic assays; however, the severity and consequences of barcode swapping remain poorly understood. We have used two statistical approaches to robustly quantify the fraction of swapped reads in two plate-based single-cell RNA-sequencing datasets. We found that approximately 2.5% of reads were mislabelled between samples on the HiSeq 4000, which is lower than previous reports. We observed no correlation between the swapped fraction of reads and the concentration of free barcode across plates. Furthermore, we have demonstrated that barcode swapping may generate complex but artefactual cell libraries in droplet-based single-cell RNA-sequencing studies. To eliminate these artefacts, we have developed an algorithm to exclude individual molecules that have swapped between samples in 10x Genomics experiments, allowing the continued use of cutting-edge sequencing machines for these assays

Targeted genomic analysis reveals widespread autoimmune disease association with regulatory variants in the TNF superfamily cytokine signalling network.

BACKGROUND: Tumour necrosis factor (TNF) superfamily cytokines and their receptors regulate diverse immune system functions through a common set of signalling pathways. Genetic variants in and expression of individual TNF superfamily cytokines, receptors and signalling proteins have been associated with autoimmune and inflammatory diseases, but their interconnected biology has been largely unexplored. METHODS: We took a hypothesis-driven approach using available genome-wide datasets to identify genetic variants regulating gene expression in the TNF superfamily cytokine signalling network and the association of these variants with autoimmune and autoinflammatory disease. Using paired gene expression and genetic data, we identified genetic variants associated with gene expression, expression quantitative trait loci (eQTLs), in four peripheral blood cell subsets. We then examined whether eQTLs were dependent on gene expression level or the presence of active enhancer chromatin marks. Using these eQTLs as genetic markers of the TNF superfamily signalling network, we performed targeted gene set association analysis in eight autoimmune and autoinflammatory disease genome-wide association studies. RESULTS: Comparison of TNF superfamily network gene expression and regulatory variants across four leucocyte subsets revealed patterns that differed between cell types. eQTLs for genes in this network were not dependent on absolute gene expression levels and were not enriched for chromatin marks of active enhancers. By examining autoimmune disease risk variants among our eQTLs, we found that risk alleles can be associated with either increased or decreased expression of co-stimulatory TNF superfamily cytokines, receptors or downstream signalling molecules. Gene set disease association analysis revealed that eQTLs for genes in the TNF superfamily pathway were associated with six of the eight autoimmune and autoinflammatory diseases examined, demonstrating associations beyond single genome-wide significant hits. CONCLUSIONS: This systematic analysis of the influence of regulatory genetic variants in the TNF superfamily network reveals widespread and diverse roles for these cytokines in susceptibility to a number of immune-mediated diseases.The Intramural Research Program of the National Institute of Arthritis and Musculoskeletal and Skin Diseases and the National Library of Medicine of the US National Institutes of Health (Intramural Research Program) , Wellcome Trust (080327/Z/06/Z, 087007/Z/08/Z, 094227/Z/10/Z, Clinical PhD Programme, 079895, 076113 and 085475) , Medical Research Council (G0400929) , National Institute for Health Research , National Institutes of Health (Oxford-Cambridge Scholars Program) , Istanbul University Research Fund and UK Behcet’s Syndrome Society.This is the final version of the article. It first appeared from BioMed Central via http://dx.doi.org/10.1186/s13073-016-0329-

Reduced monocyte and macrophage TNFSF15/TL1A expression is associated with susceptibility to inflammatory bowel disease.

Chronic inflammation in inflammatory bowel disease (IBD) results from a breakdown of intestinal immune homeostasis and compromise of the intestinal barrier. Genome-wide association studies have identified over 200 genetic loci associated with risk for IBD, but the functional mechanisms of most of these genetic variants remain unknown. Polymorphisms at the TNFSF15 locus, which encodes the TNF superfamily cytokine commonly known as TL1A, are associated with susceptibility to IBD in multiple ethnic groups. In a wide variety of murine models of inflammation including models of IBD, TNFSF15 promotes immunopathology by signaling through its receptor DR3. Such evidence has led to the hypothesis that expression of this lymphocyte costimulatory cytokine increases risk for IBD. In contrast, here we show that the IBD-risk haplotype at TNFSF15 is associated with decreased expression of the gene by peripheral blood monocytes in both healthy volunteers and IBD patients. This association persists under various stimulation conditions at both the RNA and protein levels and is maintained after macrophage differentiation. Utilizing a "recall-by-genotype" bioresource for allele-specific expression measurements in a functional fine-mapping assay, we localize the polymorphism controlling TNFSF15 expression to the regulatory region upstream of the gene. Through a T cell costimulation assay, we demonstrate that genetically regulated TNFSF15 has functional relevance. These findings indicate that genetically enhanced expression of TNFSF15 in specific cell types may confer protection against the development of IBD

Stimulation strength controls the rate of initiation but not the molecular organisation of TCR-induced signalling

Funder: European Molecular Biology Organization; FundRef: http://dx.doi.org/10.13039/100004410Millions of naïve T cells with different TCRs may interact with a peptide-MHC ligand, but very few will activate. Remarkably, this fine control is orchestrated using a limited set of intracellular machinery. It remains unclear whether changes in stimulation strength alter the programme of signalling events leading to T cell activation. Using mass cytometry to simultaneously measure multiple signalling pathways during activation of murine CD8+ T cells, we found a programme of distal signalling events that is shared, regardless of the strength of TCR stimulation. Moreover, the relationship between transcription of early response genes Nr4a1 and Irf8 and activation of the ribosomal protein S6 is also conserved across stimuli. Instead, we found that stimulation strength dictates the rate with which cells initiate signalling through this network. These data suggest that TCR-induced signalling results in a coordinated activation program, modulated in rate but not organization by stimulation strength

Prohibitin 1 Modulates Mitochondrial Stress-Related Autophagy in Human Colonic Epithelial Cells

Autophagy is an adaptive response to extracellular and intracellular stress by which cytoplasmic components and organelles, including damaged mitochondria, are degraded to promote cell survival and restore cell homeostasis. Certain genes involved in autophagy confer susceptibility to Crohn's disease. Reactive oxygen species and pro-inflammatory cytokines such as tumor necrosis factor α (TNFα), both of which are increased during active inflammatory bowel disease, promote cellular injury and autophagy via mitochondrial damage. Prohibitin (PHB), which plays a role in maintaining normal mitochondrial respiratory function, is decreased during active inflammatory bowel disease. Restoration of colonic epithelial PHB expression protects mice from experimental colitis and combats oxidative stress. In this study, we investigated the potential role of PHB in modulating mitochondrial stress-related autophagy in intestinal epithelial cells.We measured autophagy activation in response to knockdown of PHB expression by RNA interference in Caco2-BBE and HCT116 WT and p53 null cells. The effect of exogenous PHB expression on TNFα- and IFNγ-induced autophagy was assessed. Autophagy was inhibited using Bafilomycin A(1) or siATG16L1 during PHB knockdown and the affect on intracellular oxidative stress, mitochondrial membrane potential, and cell viability were determined. The requirement of intracellular ROS in siPHB-induced autophagy was assessed using the ROS scavenger N-acetyl-L-cysteine.TNFα and IFNγ-induced autophagy inversely correlated with PHB protein expression. Exogenous PHB expression reduced basal autophagy and TNFα-induced autophagy. Gene silencing of PHB in epithelial cells induces mitochondrial autophagy via increased intracellular ROS. Inhibition of autophagy during PHB knockdown exacerbates mitochondrial depolarization and reduces cell viability.Decreased PHB levels coupled with dysfunctional autophagy renders intestinal epithelial cells susceptible to mitochondrial damage and cytotoxicity. Repletion of PHB may represent a therapeutic approach to combat oxidant and cytokine-induced mitochondrial damage in diseases such as inflammatory bowel disease

Comparison of gene expression microarray data with count-based RNA measurements informs microarray interpretation.

BACKGROUND: Although numerous investigations have compared gene expression microarray platforms, preprocessing methods and batch correction algorithms using constructed spike-in or dilution datasets, there remains a paucity of studies examining the properties of microarray data using diverse biological samples. Most microarray experiments seek to identify subtle differences between samples with variable background noise, a scenario poorly represented by constructed datasets. Thus, microarray users lack important information regarding the complexities introduced in real-world experimental settings. The recent development of a multiplexed, digital technology for nucleic acid measurement enables counting of individual RNA molecules without amplification and, for the first time, permits such a study. RESULTS: Using a set of human leukocyte subset RNA samples, we compared previously acquired microarray expression values with RNA molecule counts determined by the nCounter Analysis System (NanoString Technologies) in selected genes. We found that gene measurements across samples correlated well between the two platforms, particularly for high-variance genes, while genes deemed unexpressed by the nCounter generally had both low expression and low variance on the microarray. Confirming previous findings from spike-in and dilution datasets, this "gold-standard" comparison demonstrated signal compression that varied dramatically by expression level and, to a lesser extent, by dataset. Most importantly, examination of three different cell types revealed that noise levels differed across tissues. CONCLUSIONS: Microarray measurements generally correlate with relative RNA molecule counts within optimal ranges but suffer from expression-dependent accuracy bias and precision that varies across datasets. We urge microarray users to consider expression-level effects in signal interpretation and to evaluate noise properties in each dataset independently

Homer2 within the nucleus accumbens core bidirectionally regulates alcohol intake by both P and Wistar rats

In murine models of alcoholism, the glutamate receptor scaffolding protein Homer2 bidirectionally regulates alcohol intake. Although chronic alcohol drinking increases Homer2 expression within the core subregion of the nucleus accumbens (NAc) of alcohol-preferring P rats, the relevance of this neuroadaptation for alcohol intake has yet to be determined in rats. Thus, the present study employed an adeno-associated viral vector (AAV) strategy to over-express and knock down the major rodent isoform Homer2b within the NAc of both P and outbred Wistar rats to examine for changes in alcohol preference and intake (0-30% v/v) under continuous-access procedures. The generalization of AAV effects to non-drug, palatable, sweet solutions was also determined in tests of sucrose (0-5% w/v) and saccharin (0-0.125% w/v) intake/preference. No net-flux in vivo microdialysis was conducted for glutamate in the NAc to relate Homer2-dependent changes in alcohol intake to extracellular levels of glutamate. Line differences were noted for sweet solution preference and intake, but these variables were not affected by intra-NAc AAV infusion in either line. In contrast, Homer2b over-expression elevated, while Homer2b knock-down reduced, alcohol intake in both lines, and this effect was greatest at the highest concentration. Strikingly, in P rats there was a direct association between changes in Homer2b expression and NAc extracellular glutamate levels, but this effect was not seen in Wistar rats. These data indicate that NAc Homer2b expression actively regulates alcohol consumption by rats, paralleling this previous observation in mice. Overall, these findings underscore the importance of mesocorticolimbic glutamate activity in alcohol abuse/dependence and suggest that Homer2b and/or its constituents may serve as molecular targets for the treatment of these disorders

Insight into Genotype-Phenotype Associations through eQTL Mapping in Multiple Cell Types in Health and Immune-Mediated Disease

Genome-wide association studies (GWAS) have transformed our understanding of the genetics of complex traits such as autoimmune diseases, but how risk variants contribute to pathogenesis remains largely unknown. Identifying genetic variants that affect gene expression (expression quantitative trait loci, or eQTLs) is crucial to addressing this. eQTLs vary between tissues and following in vitro cellular activation, but have not been examined in the context of human inflammatory diseases. We performed eQTL mapping in five primary immune cell types from patients with active inflammatory bowel disease (n = 91), anti-neutrophil cytoplasmic antibody-associated vasculitis (n = 46) and healthy controls (n = 43), revealing eQTLs present only in the context of active inflammatory disease. Moreover, we show that following treatment a proportion of these eQTLs disappear. Through joint analysis of expression data from multiple cell types, we reveal that previous estimates of eQTL immune cell-type specificity are likely to have been exaggerated. Finally, by analysing gene expression data from multiple cell types, we find eQTLs not previously identified by database mining at 34 inflammatory bowel disease-associated loci. In summary, this parallel eQTL analysis in multiple leucocyte subsets from patients with active disease provides new insights into the genetic basis of immune-mediated diseases.This research was funded by a Wellcome Trust Clinical PhD Programme Fellowship (JEP), the NIH-Oxford-Cambridge Scholars Program (ACR), Wellcome Trust Grant 083650/Z/07/Z and MRC Grant MR/L19027/1 (KGCS), and the National Institute for Health Research Cambridge Biomedical Research Centre. KGCS is a National Institute for Health Research Senior Investigator. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

Testing for differential abundance in mass cytometry data.

When comparing biological conditions using mass cytometry data, a key challenge is to identify cellular populations that change in abundance. Here, we present a computational strategy for detecting 'differentially abundant' populations by assigning cells to hyperspheres, testing for significant differences between conditions and controlling the spatial false discovery rate. Our method (http://bioconductor.org/packages/cydar) outperforms other approaches in simulations and finds novel patterns of differential abundance in real data.This work was supported by Cancer Research UK (core funding to J.C.M., award no. A17197), the University of Cambridge and Hutchison Whampoa Limited. J.C.M. was also supported by core funding from EMBL