ABSTRACT
Motivation: High-parameter single-cell technologies can reveal novel cell populations of interest, but studying or validating these populations using lower-parameter methods remains challenging. Results: Here, we present GateFinder, an algorithm that enriches high-dimensional cell types with simple, stepwise polygon gates requiring only two markers at a time. A series of case studies of complex cell types illustrates how simplified enrichment strategies can enable more efficient assays, reveal novel biomarkers and clarify underlying biology. Availability and implementation: The GateFinder algorithm is implemented as a free and open-source package for BioConductor: https://nalab.stanford.edu/gatefinder. Supplementary information: Supplementary data are available at Bioinformatics online.
Subject(s)
Algorithms , Biomarkers/analysis , Flow Cytometry , SoftwareABSTRACT
Delayed recovery from surgery causes personal suffering and substantial societal and economic costs. Whether immune mechanisms determine recovery after surgical trauma remains ill-defined. Single-cell mass cytometry was applied to serial whole-blood samples from 32 patients undergoing hip replacement to comprehensively characterize the phenotypic and functional immune response to surgical trauma. The simultaneous analysis of 14,000 phosphorylation events in precisely phenotyped immune cell subsets revealed uniform signaling responses among patients, demarcating a surgical immune signature. When regressed against clinical parameters of surgical recovery, including functional impairment and pain, strong correlations were found with STAT3 (signal transducer and activator of transcription), CREB (adenosine 3',5'-monophosphate response element-binding protein), and NF-κB (nuclear factor κB) signaling responses in subsets of CD14(+) monocytes (R = 0.7 to 0.8, false discovery rate <0.01). These sentinel results demonstrate the capacity of mass cytometry to survey the human immune system in a relevant clinical context. The mechanistically derived immune correlates point to diagnostic signatures, and potential therapeutic targets, that could postoperatively improve patient recovery.
Subject(s)
Arthroplasty, Replacement, Hip/adverse effects , Flow Cytometry , Hip Joint/surgery , Immunophenotyping/methods , Monocytes/immunology , Postoperative Complications/immunology , Signal Transduction/immunology , Aged , Biomarkers/blood , CREB-Binding Protein/blood , Female , Hip Joint/physiopathology , Humans , Lipopolysaccharide Receptors/blood , Male , Middle Aged , Monocytes/metabolism , NF-kappa B/blood , Phenotype , Phosphorylation , Postoperative Complications/blood , Postoperative Complications/physiopathology , Predictive Value of Tests , Protein Interaction Maps , Recovery of Function , STAT3 Transcription Factor/blood , Time Factors , Treatment OutcomeABSTRACT
Elucidation and examination of cellular subpopulations that display condition-specific behavior can play a critical contributory role in understanding disease mechanism, as well as provide a focal point for development of diagnostic criteria linking such a mechanism to clinical prognosis. Despite recent advancements in single-cell measurement technologies, the identification of relevant cell subsets through manual efforts remains standard practice. As new technologies such as mass cytometry increase the parameterization of single-cell measurements, the scalability and subjectivity inherent in manual analyses slows both analysis and progress. We therefore developed Citrus (cluster identification, characterization, and regression), a data-driven approach for the identification of stratifying subpopulations in multidimensional cytometry datasets. The methodology of Citrus is demonstrated through the identification of known and unexpected pathway responses in a dataset of stimulated peripheral blood mononuclear cells measured by mass cytometry. Additionally, the performance of Citrus is compared with that of existing methods through the analysis of several publicly available datasets. As the complexity of flow cytometry datasets continues to increase, methods such as Citrus will be needed to aid investigators in the performance of unbiased--and potentially more thorough--correlation-based mining and inspection of cell subsets nested within high-dimensional datasets.
Subject(s)
Algorithms , Cells/classification , Cells/cytology , Computational Biology/methods , Flow Cytometry/methods , Software , Blood Cells/cytology , Humans , T-Lymphocyte Subsets/cytologyABSTRACT
Mass cytometry facilitates high-dimensional, quantitative analysis of the effects of bioactive molecules on human samples at single-cell resolution, but instruments process only one sample at a time. Here we describe mass-tag cellular barcoding (MCB), which increases mass cytometry throughput by using n metal ion tags to multiplex up to 2n samples. We used seven tags to multiplex an entire 96-well plate, and applied MCB to characterize human peripheral blood mononuclear cell (PBMC) signaling dynamics and cell-to-cell communication, signaling variability between PBMCs from eight human donors, and the effects of 27 inhibitors on this system. For each inhibitor, we measured 14 phosphorylation sites in 14 PBMC types at 96 conditions, resulting in 18,816 quantified phosphorylation levels from each multiplexed sample. This high-dimensional, systems-level inquiry allowed analysis across cell-type and signaling space, reclassified inhibitors and revealed off-target effects. High-content, high-throughput screening with MCB should be useful for drug discovery, preclinical testing and mechanistic investigation of human disease.
Subject(s)
Flow Cytometry/methods , High-Throughput Screening Assays/methods , Leukocytes, Mononuclear/cytology , Systems Biology/methods , Chelating Agents , Heterocyclic Compounds, 1-Ring , Humans , K562 Cells , Leukocytes, Mononuclear/drug effects , Leukocytes, Mononuclear/metabolism , Phosphorylation , Principal Component Analysis , Protein Kinase Inhibitors/pharmacology , Signal TransductionABSTRACT
MOTIVATION: Recent advances in flow cytometry enable simultaneous single-cell measurement of 30+ surface and intracellular proteins. CytoSPADE is a high-performance implementation of an interface for the Spanning-tree Progression Analysis of Density-normalized Events algorithm for tree-based analysis and visualization of this high-dimensional cytometry data. AVAILABILITY: Source code and binaries are freely available at http://cytospade.org and via Bioconductor version 2.10 onwards for Linux, OSX and Windows. CytoSPADE is implemented in R, C++ and Java. CONTACT: michael.linderman@mssm.edu SUPPLEMENTARY INFORMATION: Additional documentation available at http://cytospade.org.
Subject(s)
Algorithms , Flow Cytometry/methods , Software , Computer GraphicsABSTRACT
The ability to analyze multiple single-cell parameters is critical for understanding cellular heterogeneity. Despite recent advances in measurement technology, methods for analyzing high-dimensional single-cell data are often subjective, labor intensive and require prior knowledge of the biological system. To objectively uncover cellular heterogeneity from single-cell measurements, we present a versatile computational approach, spanning-tree progression analysis of density-normalized events (SPADE). We applied SPADE to flow cytometry data of mouse bone marrow and to mass cytometry data of human bone marrow. In both cases, SPADE organized cells in a hierarchy of related phenotypes that partially recapitulated well-described patterns of hematopoiesis. We demonstrate that SPADE is robust to measurement noise and to the choice of cellular markers. SPADE facilitates the analysis of cellular heterogeneity, the identification of cell types and comparison of functional markers in response to perturbations.
Subject(s)
Algorithms , Bone Marrow Cells/cytology , Flow Cytometry/methods , Adult , Animals , Biomarkers/metabolism , Bone Marrow Cells/drug effects , Bone Marrow Cells/metabolism , Computer Simulation , Databases as Topic , Hematopoiesis/drug effects , Humans , Mice , Staining and Labeling , Tumor Necrosis Factor-alpha/pharmacologyABSTRACT
Flow cytometry is an essential tool for dissecting the functional complexity of hematopoiesis. We used single-cell "mass cytometry" to examine healthy human bone marrow, measuring 34 parameters simultaneously in single cells (binding of 31 antibodies, viability, DNA content, and relative cell size). The signaling behavior of cell subsets spanning a defined hematopoietic hierarchy was monitored with 18 simultaneous markers of functional signaling states perturbed by a set of ex vivo stimuli and inhibitors. The data set allowed for an algorithmically driven assembly of related cell types defined by surface antigen expression, providing a superimposable map of cell signaling responses in combination with drug inhibition. Visualized in this manner, the analysis revealed previously unappreciated instances of both precise signaling responses that were bounded within conventionally defined cell subsets and more continuous phosphorylation responses that crossed cell population boundaries in unexpected manners yet tracked closely with cellular phenotype. Collectively, such single-cell analyses provide system-wide views of immune signaling in healthy human hematopoiesis, against which drug action and disease can be compared for mechanistic studies and pharmacologic intervention.
Subject(s)
Bone Marrow Cells/drug effects , Bone Marrow Cells/immunology , Flow Cytometry/methods , Lymphocyte Subsets/drug effects , Lymphocyte Subsets/immunology , Pyrimidines/pharmacology , Signal Transduction , Single-Cell Analysis/methods , Thiazoles/pharmacology , Algorithms , Antibodies , Antigens, Surface/analysis , B-Lymphocytes/drug effects , B-Lymphocytes/immunology , B-Lymphocytes/metabolism , Bone Marrow Cells/cytology , Bone Marrow Cells/metabolism , Cytokines/metabolism , Dasatinib , Hematopoiesis , Humans , Immunophenotyping , Lanthanoid Series Elements , Leukocytes, Mononuclear/drug effects , Leukocytes, Mononuclear/immunology , Leukocytes, Mononuclear/metabolism , Lymphocyte Activation , Lymphocyte Subsets/metabolism , Mass Spectrometry , Phosphorylation , Protein Kinase Inhibitors/pharmacology , Protein-Tyrosine Kinases/antagonists & inhibitors , Signal Transduction/drug effects , T-Lymphocytes/drug effects , T-Lymphocytes/immunology , T-Lymphocytes/metabolism , Transition ElementsABSTRACT
Aberrant intracellular signaling plays an important role in many diseases. The causal structure of signal transduction networks can be modeled as Bayesian Networks (BNs), and computationally learned from experimental data. However, learning the structure of Bayesian Networks (BNs) is an NP-hard problem that, even with fast heuristics, is too time consuming for large, clinically important networks (20-50 nodes). In this paper, we present a novel graphics processing unit (GPU)-accelerated implementation of a Monte Carlo Markov Chain-based algorithm for learning BNs that is up to 7.5-fold faster than current general-purpose processor (GPP)-based implementations. The GPU-based implementation is just one of several implementations within the larger application, each optimized for a different input or machine configuration. We describe the methodology we use to build an extensible application, assembled from these variants, that can target a broad range of heterogeneous systems, e.g., GPUs, multicore GPPs. Specifically we show how we use the Merge programming model to efficiently integrate, test and intelligently select among the different potential implementations.
ABSTRACT
VectorBase (http://www.vectorbase.org) is an NIAID-funded Bioinformatic Resource Center focused on invertebrate vectors of human pathogens. VectorBase annotates and curates vector genomes providing a web accessible integrated resource for the research community. Currently, VectorBase contains genome information for three mosquito species: Aedes aegypti, Anopheles gambiae and Culex quinquefasciatus, a body louse Pediculus humanus and a tick species Ixodes scapularis. Since our last report VectorBase has initiated a community annotation system, a microarray and gene expression repository and controlled vocabularies for anatomy and insecticide resistance. We have continued to develop both the software infrastructure and tools for interrogating the stored data.
Subject(s)
Arthropod Vectors/genetics , Culicidae/genetics , Databases, Genetic , Aedes/genetics , Animals , Anopheles/genetics , Culex/genetics , Culicidae/metabolism , Gene Expression Profiling , Genome, Insect , Genomics , Ixodes/genetics , Pediculus/genetics , Vocabulary, ControlledABSTRACT
We present a draft sequence of the genome of Aedes aegypti, the primary vector for yellow fever and dengue fever, which at approximately 1376 million base pairs is about 5 times the size of the genome of the malaria vector Anopheles gambiae. Nearly 50% of the Ae. aegypti genome consists of transposable elements. These contribute to a factor of approximately 4 to 6 increase in average gene length and in sizes of intergenic regions relative to An. gambiae and Drosophila melanogaster. Nonetheless, chromosomal synteny is generally maintained among all three insects, although conservation of orthologous gene order is higher (by a factor of approximately 2) between the mosquito species than between either of them and the fruit fly. An increase in genes encoding odorant binding, cytochrome P450, and cuticle domains relative to An. gambiae suggests that members of these protein families underpin some of the biological differences between the two mosquito species.
Subject(s)
Aedes/genetics , Genome, Insect , Insect Vectors/genetics , Aedes/metabolism , Animals , Anopheles/genetics , Anopheles/metabolism , Arboviruses , Base Sequence , DNA Transposable Elements , Dengue/prevention & control , Dengue/transmission , Drosophila melanogaster/genetics , Female , Genes, Insect , Humans , Insect Proteins/genetics , Insect Vectors/metabolism , Male , Membrane Transport Proteins/genetics , Molecular Sequence Data , Multigene Family , Protein Structure, Tertiary/genetics , Sequence Analysis, DNA , Sex Characteristics , Sex Determination Processes , Species Specificity , Synteny , Transcription, Genetic , Yellow Fever/prevention & control , Yellow Fever/transmissionABSTRACT
BACKGROUND: The genome of Anopheles gambiae, the major vector of malaria, was sequenced and assembled in 2002. This initial genome assembly and analysis made available to the scientific community was complicated by the presence of assembly issues, such as scaffolds with no chromosomal location, no sequence data for the Y chromosome, haplotype polymorphisms resulting in two different genome assemblies in limited regions and contaminating bacterial DNA. RESULTS: Polytene chromosome in situ hybridization with cDNA clones was used to place 15 unmapped scaffolds (sizes totaling 5.34 Mbp) in the pericentromeric regions of the chromosomes and oriented a further 9 scaffolds. Additional analysis by in situ hybridization of bacterial artificial chromosome (BAC) clones placed 1.32 Mbp (5 scaffolds) in the physical gaps between scaffolds on euchromatic parts of the chromosomes. The Y chromosome sequence information (0.18 Mbp) remains highly incomplete and fragmented among 55 short scaffolds. Analysis of BAC end sequences showed that 22 inter-scaffold gaps were spanned by BAC clones. Unmapped scaffolds were also aligned to the chromosome assemblies in silico, identifying regions totaling 8.18 Mbp (144 scaffolds) that are probably represented in the genome project by two alternative assemblies. An additional 3.53 Mbp of alternative assembly was identified within mapped scaffolds. Scaffolds comprising 1.97 Mbp (679 small scaffolds) were identified as probably derived from contaminating bacterial DNA. In total, about 33% of previously unmapped sequences were placed on the chromosomes. CONCLUSION: This study has used new approaches to improve the physical map and assembly of the A. gambiae genome.
Subject(s)
Anopheles/genetics , Genome, Insect/genetics , Physical Chromosome Mapping , Animals , Bacteria/genetics , Centromere/genetics , Chromosomes/genetics , Euchromatin/genetics , Polymorphism, Genetic , Species SpecificityABSTRACT
VectorBase (http://www.vectorbase.org/) is a web-accessible data repository for information about invertebrate vectors of human pathogens. VectorBase annotates and maintains vector genomes providing an integrated resource for the research community. Currently, VectorBase contains genome information for two organisms: Anopheles gambiae, a vector for the Plasmodium protozoan agent causing malaria, and Aedes aegypti, a vector for the flaviviral agents causing Yellow fever and Dengue fever.