scholarly journals Benchmarking of a Bayesian single cell RNAseq differential gene expression test for dose-response study designs

2021 ◽  
Author(s):  
Rance Nault ◽  
Satabdi Saha ◽  
Sudin Bhattacharya ◽  
Jack Dodson ◽  
Samiran Sinha ◽  
...  

AbstractThe application of single-cell RNA sequencing (scRNAseq) for the evaluation of chemicals, drugs, and food contaminants presents the opportunity to consider cellular heterogeneity in pharmacological and toxicological responses. Current differential gene expression analysis (DGEA) methods focus primarily on two group comparisons, not multi-group dose-response study designs used in safety assessments. To benchmark DGEA methods for dose-response scRNAseq experiments, we proposed a multiplicity corrected Bayesian testing approach and compare it against 8 other methods including two frequentist fit-for-purpose tests using simulated and experimental data. Our Bayesian test method outperformed all other tests for a broad range of accuracy metrics including control of false positive error rates. Most notable, the fit-for-purpose and standard multiple group DGEA methods were superior to the two group scRNAseq methods for dose-response study designs. Collectively, our benchmarking of DGEA methods demonstrates the importance in considering study design when determining the most appropriate test methods.

2019 ◽  
Vol 20 (S24) ◽  
Author(s):  
Yu Zhang ◽  
Changlin Wan ◽  
Pengcheng Wang ◽  
Wennan Chang ◽  
Yan Huo ◽  
...  

Abstract Background Various statistical models have been developed to model the single cell RNA-seq expression profiles, capture its multimodality, and conduct differential gene expression test. However, for expression data generated by different experimental design and platforms, there is currently lack of capability to determine the most proper statistical model. Results We developed an R package, namely Multi-Modal Model Selection (M3S), for gene-wise selection of the most proper multi-modality statistical model and downstream analysis, useful in a single-cell or large scale bulk tissue transcriptomic data. M3S is featured with (1) gene-wise selection of the most parsimonious model among 11 most commonly utilized ones, that can best fit the expression distribution of the gene, (2) parameter estimation of a selected model, and (3) differential gene expression test based on the selected model. Conclusion A comprehensive evaluation suggested that M3S can accurately capture the multimodality on simulated and real single cell data. An open source package and is available through GitHub at https://github.com/zy26/M3S.


2019 ◽  
Author(s):  
Kelly M. Bakulski ◽  
John F. Dou ◽  
Robert C. Thompson ◽  
Christopher Lee ◽  
Lauren Y. Middleton ◽  
...  

AbstractBackgroundLead (Pb) exposure is ubiquitous and has permanent developmental effects on childhood intelligence and behavior and adulthood risk of dementia. The hippocampus is a key brain region involved in learning and memory, and its cellular composition is highly heterogeneous. Pb acts on the hippocampus by altering gene expression, but the cell type-specific responses are unknown.ObjectiveExamine the effects of perinatal Pb treatment on adult hippocampus gene expression, at the level of individual cells, in mice.MethodsIn mice perinatally exposed to control water (n=4) or a human physiologically-relevant level (32 ppm in maternal drinking water) of Pb (n=4), two weeks prior to mating through weaning, we tested for gene expression and cellular differences in the hippocampus at 5-months of age. Analysis was performed using single cell RNA-sequencing of 5,258 cells from the hippocampus by 10x Genomics Chromium to 1) test for gene expression differences averaged across all cells by treatment; 2) compare cell cluster composition by treatment; and 3) test for gene expression and pathway differences within cell clusters by treatment.ResultsGene expression patterns revealed 12 cell clusters in the hippocampus, mapping to major expected cell types (e.g. microglia, astrocytes, neurons, oligodendrocytes). Perinatal Pb treatment was associated with 12.4% more oligodendrocytes (P=4.4×10−21) in adult mice. Across all cells, differential gene expression analysis by Pb treatment revealed cluster marker genes. Within cell clusters, differential gene expression with Pb treatment (q<0.05) was observed in endothelial, microglial, pericyte, and astrocyte cells. Pathways up-regulated with Pb treatment were protein folding in microglia (P=3.4×10−9) and stress response in oligodendrocytes (P=3.2×10−5).ConclusionBulk tissue analysis may be confounded by changes in cell type composition and may obscure effects within vulnerable cell types. This study serves as a biological reference for future single cell studies of toxicant or neuronal complications, to ultimately characterize the molecular basis by which Pb influences cognition and behavior.


Blood ◽  
2018 ◽  
Vol 132 (Supplement 1) ◽  
pp. 5087-5087
Author(s):  
Francis Guitart ◽  
Moonjung Jung ◽  
Stefan Cordes ◽  
Shiqin Yu ◽  
Jizhong Zou ◽  
...  

Abstract GATA2 deficiency is a rare, inherited or sporadic genetic disorder characterized by variable onset of a pleomorphic constellation of immune, hematologic and lymphatic abnormalities linked to heterozygous mutations in the Gata2 gene. Patients develop monocyte, B cell, NK cell and dendritic cell deficiencies resulting in vulnerabilities to unusual infections. Patients with GATA2 deficiency also frequently progress to bone marrow failure, myelodysplastic syndrome and/or acute myelogenous leukemia. GATA proteins are transcription factors with central roles in early embryonic development and lineage specification. GATA2 is a master regulator of hematopoiesis, implicated in the initial generation and maintenance of hematopoietic stem cells (HSC). Murine models recapitulate the human phenotype incompletely: GATA2 heterozygous knockout mice do not manifest loss of monocyte, B cells or NK cells; however, serial repopulation assays show decreased engraftment potential. Direct studies of primary HSC from patients with GATA2 deficiency are challenging due to the generally hypocellular marrow. We hypothesized that human pluripotent stem cells, particularly patient-specific iPSC, could be used to study potential developmental defects in GATA2 deficiency, overcoming a lack of primary HSC. In order to gain insights into the impact of human GATA2 deficiency on hematopoietic differentiation, we compared the single cell transcriptomes of HSPC differentiated from (i) iPSCs from a patient with GATA2 deficiency due to a mutation p.R337X (c.1009C>T) (ii) isogenic iPSCs created via homology-directed repair of Gata2 p.R337X, (iii) iPSCs from a healthy control and (iv) isogenic Gata2 heterozygous mutant iPSCs with monoallelic frameshift mutations in the second zinc finger domain. Mesodermal and hematopoietic differentiation was performed under feeder-free, defined media conditions. At day 0, iPSCs were plated in mesodermal induction media containing VEGF, SCF, Activin A and Y27632 in STEMdiff APEL media. Mesodermal induction was continued until day 4, when embryoid bodies were cultured in hematopoietic specification media with SCF, FLT3L, IL3, IL6, G-CSF and BMP4 until day 16, when CD34+CD45+ iPSC-derived hematopoietic stem and progenitor cells (iHSC) were enumerated and sorted by fluorescence-activated cell sorting. Single cell RNA-seq was performed using the 10XGenomics Chromium platform and primary analysis via CellRanger. Scater was used to filter outlier cells. Seurat as used to compute multiple manifold alignment and differential gene expression. Cell classification, pseudotemporal ordering and branch point analysis were performed with monocle. URD was used to calculate confirmatory diffusion maps and pseudotemporal ordering. We analyzed 7,855 iHSPC (2952 from GATA2-deficient patient, 241 isogenic iHSPCs after repair of Gata2 mutation, 2,605 from a healthy volunteer and 2,057 from isogenic heterozygous Gata2 knockout iHSPCs) after filtering of outliers. We computed multiple manifold alignment to mitigate batch effects. Differential gene expression across Gata2 mutation status found that 42 out of 102 (42%) target genes of GATA2 (c.f. TRANSFAC database of curated transcription factor targets) were differentially expressed with adjusted p-values less than 0.05. Semi-supervised classification of cell-types and pseudotemporal ordering via monocle revealed two branch points, consistent with developmental branchings at the level of CLP and CMP multipotent progenitors. The numbers of cells along each branch was found to be statistically different (χ2=30.07, p-value = 3e-7) with the biggest differences noted in the lymphoid branch (state 4). Differential gene expression in this branch revealed a differential up-regulation of Notch1, CD69 and FKBPs and differential down-regulation of CD14. In conclusion, iPSC/iHSPC differentiation models combined with single cell transcriptome analysis may be a valuable tool to identify pathways responsible for impaired hematopoietic/lymphatic development in GATA2 deficiency. Figure. Figure. Disclosures Dunbar: National Institute of Health: Research Funding. Winkler:National Institute of Health: Research Funding.


2020 ◽  
Vol 17 (6) ◽  
pp. 583-585
Author(s):  
Etienne Becht ◽  
Edward Zhao ◽  
Robert Amezquita ◽  
Raphael Gottardo

mBio ◽  
2020 ◽  
Vol 11 (4) ◽  
Author(s):  
Rosiris León-Rivera ◽  
Brenda Morsey ◽  
Meng Niu ◽  
Howard S. Fox ◽  
Joan W. Berman

ABSTRACT HIV reservoirs persist despite successful antiretroviral therapy (ART) and are a major obstacle to the eradication and cure of HIV. The mature monocyte subset, CD14+CD16+, contributes to viral reservoirs and HIV-associated comorbidities. Only a subset of monocytes harbors HIV (HIV+), while the rest remain uninfected, exposed cells (HIVexp). We developed an innovative single cell RNA sequencing (scRNAseq) pipeline that detects HIV and host transcripts simultaneously, enabling us to examine differences between HIV+ and HIVexp mature monocytes. Using this, we characterized uninfected, HIV+, and HIVexp primary human mature monocytes with and without ART. We showed that HIV+ mature monocytes do not form their own cluster separately from HIVexp but can be distinguished by significant differential gene expression. We found that ART decreased levels of unspliced HIV transcripts potentially by modulating host transcriptional regulators shown to decrease viral infection and replication. We also identified and characterized mature monocyte subpopulations differentially impacted by HIV and ART. We identified genes dysregulated by ART in HIVexp monocytes compared to their uninfected counterpart and, of interest, the junctional protein ALCAM, suggesting that ART impacts monocyte functions. Our data provide a novel method for simultaneous detection of HIV and host transcripts. We identify potential targets, such as those genes whose expression is increased in HIV+ mature monocytes compared to HIVexp, to block their entry into tissues, preventing establishment/replenishment of HIV reservoirs even with ART, thereby reducing and/or eliminating viral burden and HIV-associated comorbidities. Our data also highlight the heterogeneity of mature monocyte subsets and their potential contributions to HIV pathogenesis in the ART era. IMPORTANCE HIV enters tissues early after infection, leading to establishment and persistence of HIV reservoirs despite antiretroviral therapy (ART). Viral reservoirs are a major obstacle to the eradication and cure of HIV. CD14+CD16+ (mature) monocytes may contribute to establishment and reseeding of reservoirs. A subset of monocytes, consisting mainly of CD14+CD16+ cells, harbors HIV (HIV+), while the rest remain uninfected, exposed cells (HIVexp). It is important to identify cells harboring virus to eliminate reservoirs. Using an innovative single-cell RNA sequencing (scRNAseq) pipeline to detect HIV and host transcripts simultaneously, we characterized HIV+ and HIVexp primary human mature monocytes with and without ART. HIV+ mature monocytes are not a unique subpopulation but rather can be distinguished from HIVexp by differential gene expression. We characterized mature monocyte subpopulations differently impacted by HIV and ART, highlighting their potential contributions to HIV-associated comorbidities. Our data propose therapeutic targets to block HIV+ monocyte entry into tissues, preventing establishment and replenishment of reservoirs even with ART.


Author(s):  
Peng He ◽  
Brian A. Williams ◽  
Diane Trout ◽  
Georgi K. Marinov ◽  
Henry Amrhein ◽  
...  

AbstractIn mammalian embryogenesis differential gene expression gradually builds the identity and complexity of each tissue and organ system. We systematically quantified mouse polyA-RNA from embryo day E10.5 to birth, sampling 17 whole tissues, enhanced with single-cell measurements for the developing limb. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets, characterized by their promoters’ transcription factor (TF) motif codes. We decomposed the tissue-level transcriptome using scRNA-seq and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for 1/3 of differential gene expression and over 40% of identified cell types. Integrating promoter sequence motifs with companion ENCODE epigenomic profiles identified a promoter de-repression mechanism unique to neuronal expression clusters and attributable to known and novel repressors. Focusing on the developing limb, scRNA-seq identified 25 known and candidate novel cell types, including progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell type TF networks and complementary sets of candidate enhancer elements by de-convolving whole-tissue IDEAS epigenome chromatin state models. These ENCODE reference data, computed network components and IDEAS chromatin segmentations, are companion resources to the matching epigenomic developmental matrix, available for researchers to further mine and integrate.


Sign in / Sign up

Export Citation Format

Share Document