Characterizing cell subsets using marker enrichment modeling.
Ontology highlight
ABSTRACT:
SUBMITTER: Diggins KE
PROVIDER: S-EPMC5330853 | biostudies-literature | 2017 Mar
REPOSITORIES: biostudies-literature
Publications
Characterizing cell subsets using marker enrichment modeling.
Diggins Kirsten E KE Greenplate Allison R AR Leelatian Nalin N Wogsland Cara E CE Irish Jonathan M JM
Nature methods 20170130 3
Learning cell identity from high-content single-cell data presently relies on human experts. We present marker enrichment modeling (MEM), an algorithm that objectively describes cells by quantifying contextual feature enrichment and reporting a human- and machine-readable text label. MEM outperforms traditional metrics in describing immune and cancer cell subsets from fluorescence and mass cytometry. MEM provides a quantitative language to communicate characteristics of new and established cytot ...[more]
Similar Datasets
Project description:ATAC-seq has become a leading technology for probing the chromatin landscape of single and aggregated cells. Distilling functional regions from ATAC-seq presents diverse analysis challenges. Methods commonly used to analyze chromatin accessibility datasets are adapted from algorithms designed to process different experimental technologies, disregarding the statistical and biological differences intrinsic to the ATAC-seq technology. Here, we present a Bayesian statistical approach that uses latent space models to better model accessible regions, termed ChromA. ChromA annotates chromatin landscape by integrating information from replicates, producing a consensus de-noised annotation of chromatin accessibility. ChromA can analyze single cell ATAC-seq data, correcting many biases generated by the sparse sampling inherent in single cell technologies. We validate ChromA on multiple technologies and biological systems, including mouse and human immune cells, establishing ChromA as a top performing general platform for mapping the chromatin landscape in different cellular populations from diverse experimental designs.
| S-EPMC7004981 | biostudies-literature
Project description:Our increased understanding of the mechanistic heterogeneity of diseases has pushed the development of targeted therapeutics. We do not expect all patients with a given disease to benefit from a targeted drug; only those in the target population. That is, those with sufficient dysregulation in the biomolecular pathway targeted by treatment. However, due to complexity of the pathway, and/or technical issues with our characterizing assay, it is often hard to characterize the target population until well into large-scale clinical trials. This has stimulated the development of adaptive enrichment trials; clinical trials in which the target population is adaptively learned; and enrollment criteria are adaptively updated to reflect this growing understanding. This paper proposes a framework for group-sequential adaptive enrichment trials. Building on the work of Simon & Simon (2013). Adaptive enrichment designs for clinical trials. Biostatistics 14(4), 613-625), it includes a frequentist hypothesis test at the end of the trial. However, it uses Bayesian methods to optimize the decisions required during the trial (regarding how to restrict enrollment) and Bayesian methods to estimate effect size, and characterize the target population at the end of the trial. This joint frequentist/Bayesian design combines the power of Bayesian methods for decision making with the use of a formal hypothesis test at the end of the trial to preserve the studywise probability of a type I error.
| S-EPMC6075009 | biostudies-literature
Project description:We designed a mathematical model to describe and quantify the mechanisms and dynamics of tumor growth, cell-kill and resistance as they affect durations of benefit after cancer treatment. Our aim was to explore how treatment efficacy may be related to primary tumor characteristics, with the potential to guide future trial design and appropriate selection of therapy. Assuming a log-normal distribution of both resistant disease and tumor doubling times generates disease-free survival (DFS) or invasive DFS curves with specific shapes. Using a multivariate mathematical model, both treatment and tumor characteristics are related to quantified resistant disease and tumor regrowth rates by allowing different mean values for the influence of different treatments or clinical subtypes on these two log-normal distributions. Application of the model to the CALGB 9741 adjuvant breast cancer trial showed that dose-dense therapy was estimated to achieve an extra 3/4 log of cell-kill compared to standard therapy, but only in patients with more rapidly growing ER-negative tumors. Application of the model to the AZURE trial of adjuvant bisphosphonate treatment suggested that the 5-year duration of zoledronic acid was adequate for ER-negative tumors, but may not be so for ER-positive cases, with increased recurrences after ceasing the intervention. Mathematical models can identify different effects of treatment by subgroup and may aid in treatment design, trial analysis, and appropriate selection of therapy. They may provide a more appropriate and insightful tool than the conventional Cox model for the statistical analysis of response durations.
| S-EPMC4729810 | biostudies-other
Project description:We have developed a technique to directly quantify cell-substrate adhesion force using micropipette aspiration. The micropipette is positioned perpendicular to the surface of an adherent cell and a constant-rate aspiration pressure is applied. Since the micropipette diameter and the aspiration pressure are our control parameters, we have direct knowledge of the aspiration force, whereas the cell behavior is monitored either in brightfield or interference reflection microscopy. This setup thus allows us to explore a range of geometric parameters, such as projected cell area, adhesion area, or pipette size, as well as dynamical parameters such as the loading rate. We find that cell detachment is a well-defined event occurring at a critical aspiration pressure, and that the detachment force scales with the cell adhesion area (for a given micropipette diameter and loading rate), which defines a critical stress. Taking into account the cell adhesion area, intrinsic parameters of the adhesion bonds, and the loading rate, a minimal model provides an expression for the critical stress that helps rationalize our experimental results.
| S-EPMC4621874 | biostudies-other
Project description:Decreasing trends in measles mortality have been reported in recent years. However, such estimates of measles mortality have depended heavily on assumed regional measles case fatality risks (CFRs) and made little use of mortality data from low- and middle-income countries in general and India, the country with the highest measles burden globally, in particular.We constructed a dynamic model of measles transmission in India with parameters that were empirically inferred using spectral analysis from a time series of measles mortality extracted from the Million Death Study, an ongoing longitudinal study recording deaths across 2.4 million Indian households and attributing causes of death using verbal autopsy. The model was then used to estimate the measles CFR, the number of measles deaths, and the impact of vaccination in 2000-2015 among under-five children in India and in the states of Bihar and Uttar Pradesh (UP), two states with large populations and the highest numbers of measles deaths in India.We obtained the following estimated CFRs among under-five children for the year 2005: 0.63% (95% confidence interval (CI): 0.40-1.00%) for India as a whole, 0.62% (0.38-1.00%) for Bihar, and 1.19% (0.80-1.75%) for UP. During 2000-2015, we estimated that 607,000 (95% CI: 383,000-958,000) under-five deaths attributed to measles occurred in India as a whole. If no routine vaccination or supplemental immunization activities had occurred from 2000 to 2015, an additional 1.6 (1.0-2.6) million deaths for under-five children would have occurred across India.We developed a data- and model-driven estimation of the historical measles dynamics, CFR, and vaccination impact in India, extracting the periodicity of epidemics using spectral and coherence analysis, which allowed us to infer key parameters driving measles transmission dynamics and mortality.
| S-EPMC5550950 | biostudies-other
Project description:Gait asymmetry is often observed in populations with varying degrees of neuromuscular control. While changes in vertical ground reaction force (vGRF) peak magnitude are associated with altered limb loading that can be observed during asymmetric gait, the challenge is identifying techniques with the sensitivity to detect these altered movement patterns. Autoregressive (AR) modeling has successfully delineated between healthy and pathological gait during running; but has been little explored in walking. Thus, AR modeling was implemented to assess differences in vGRF pattern dynamics during symmetric and asymmetric walking. We hypothesized that the AR model coefficients would better detect differences amongst the symmetric and asymmetric walking conditions than the vGRF peak magnitude mean. Seventeen healthy individuals performed a protocol that involved walking on a split-belt instrumented treadmill at different symmetric (0.75m/s, 1.0 m/s, 1.5 m/s) and asymmetric (Side 1: 0.75m/s-Side 2:1.0 m/s; Side 1:1.0m/s-Side 2:1.5 m/s) gait conditions. Vertical ground reaction force peaks extracted during the weight-acceptance and propulsive phase of each step were used to construct a vGRF peak time series. Then, a second order AR model was fit to the vGRF peak waveform data to determine the AR model coefficients. The resulting AR coefficients were plotted on a stationarity triangle and their distance from the triangle centroid was computed. Significant differences in vGRF patterns were detected amongst the symmetric and asymmetric conditions using the AR modeling coefficients (p = 0.01); however, no differences were found when comparing vGRF peak magnitude means. These findings suggest that AR modeling has the sensitivity to identify differences in gait asymmetry that could aid in monitoring rehabilitation progression.
| S-EPMC7714243 | biostudies-literature
Project description:MotivationGene set enrichment methods are a common tool to improve the interpretability of gene lists as obtained, for example, from differential gene expression analyses. They are based on computing whether dysregulated genes are located in certain biological pathways more often than expected by chance. Gene set enrichment tools rely on pre-existing pathway databases such as KEGG, Reactome, or the Gene Ontology. These databases are increasing in size and in the number of redundancies between pathways, which complicates the statistical enrichment computation.ResultsWe address this problem and develop a novel gene set enrichment method, called pareg, which is based on a regularized generalized linear model and directly incorporates dependencies between gene sets related to certain biological functions, for example, due to shared genes, in the enrichment computation. We show that pareg is more robust to noise than competing methods. Additionally, we demonstrate the ability of our method to recover known pathways as well as to suggest novel treatment targets in an exploratory analysis using breast cancer samples from TCGA.Availability and implementationpareg is freely available as an R package on Bioconductor (https://bioconductor.org/packages/release/bioc/html/pareg.html) as well as on https://github.com/cbg-ethz/pareg. The GitHub repository also contains the Snakemake workflows needed to reproduce all results presented here.
| S-EPMC10471899 | biostudies-literature
Project description:ChipCytometry is a multiplex imaging method that can be used to analyze either cell suspensions or tissue sections. Images are acquired by iterative cycles of immunostaining with fluorescently labeled Abs, followed by photobleaching, which allows the accumulation of multiple markers on a single sample. In this study, we explored the feasibility of using ChipCytometry to identify and phenotype cell subsets, including rare cell types, using a combination of tissue sections and single-cell suspensions. Using ChipCytometry of tissue sections, we successfully demonstrated the architecture of human palatine tonsils, including the B and T cell zones, and characterized subcompartments such as the B cell mantle and germinal center zone, as well as intrafollicular PD1-expressing CD4+ T cells. Additionally, we were able to identify the rare tonsillar T cell subsets, mucosal-associated invariant T (MAIT) and γδ-T cells, within tonsil tissue. Using single-cell suspension ChipCytometry, we further dissected human tonsillar T cell subsets via unsupervised clustering analysis as well as supervised traditional manual gating. We were able to show that PD1+CD4+ T cells are comprised of CXCR5+BCL6high follicular Th cells and CXCR5-BCL6mid pre-follicular Th cells. Both supervised and unsupervised analysis approaches identified MAIT cells in single-cell suspensions, confirming a phenotype similar to that of blood-derived MAIT cells. In this study, we demonstrate that ChipCytometry is a viable method for single-cell suspension cytometry and analysis, with the additional benefit of allowing phenotyping in a spatial context using tissue sections.
| S-EPMC8278278 | biostudies-literature
Project description:Mucosal-associated invariant T (MAIT) cells represent up to 10% of circulating human T cells. They are usually defined using combinations of non-lineage-specific (surrogate) markers such as anti-TRAV1-2, CD161, IL-18Rα and CD26. The development of MR1-Ag tetramers now permits the specific identification of MAIT cells based on T-cell receptor specificity. Here, we compare these approaches for identifying MAIT cells and show that surrogate markers are not always accurate in identifying these cells, particularly the CD4+ fraction. Moreover, while all MAIT cell subsets produced comparable levels of IFNγ, TNF and IL-17A, the CD4+ population produced more IL-2 than the other subsets. In a human ontogeny study, we show that the frequencies of most MR1 tetramer+ MAIT cells, with the exception of CD4+ MAIT cells, increased from birth to about 25 years of age and declined thereafter. We also demonstrate a positive association between the frequency of MAIT cells and other unconventional T cells including Natural Killer T (NKT) cells and Vδ2+ γδ T cells. Accordingly, this study demonstrates that MAIT cells are phenotypically and functionally diverse, that surrogate markers may not reliably identify all of these cells, and that their numbers are regulated in an age-dependent manner and correlate with NKT and Vδ2+ γδ T cells.
| S-EPMC6446826 | biostudies-literature
Project description:Over the past three decades, obesity-related diseases have increased tremendously in China, and are now the leading causes of morbidity and mortality. Patterns of weight change can be used to predict risk of obesity-related diseases, increase understanding of etiology of disease risk, identify groups at particularly high risk, and shape prevention strategies.Latent class trajectory modeling was used to compute weight change trajectories for adults aged 18 to 66 using the China Health and Nutrition Survey (CHNS) data (n = 12,611). Weight change trajectories were computed separately for males and females by age group at baseline due to differential age-related patterns of weight gain in China with urbanization. Generalized linear mixed effects models examined the association between weight change trajectories and baseline characteristics including urbanicity, BMI category, age, and year of study entry.Trajectory classes were identified for each of six age-sex subgroups corresponding to various degrees of weight loss, maintenance and weight gain. Baseline BMI status was a significant predictor of trajectory membership for all age-sex subgroups. Baseline overweight/obesity increased odds of following 'initial loss with maintenance' trajectories. We found no significant association between baseline urbanization and trajectory membership after controlling for other covariates.Trajectory analysis identified patterns of weight change for age by gender groups. Lack of association between baseline urbanization status and trajectory membership suggests that living in a rural environment at baseline was not protective. Analyses identified age-specific nuances in weight change patterns, pointing to the importance of subgroup analyses in future research.
| S-EPMC4336139 | biostudies-literature