Single-Cell RNA Sequencing Efficiently Predicts Transcription Factor Targets in Plants
Ontology highlight
ABSTRACT:
SUBMITTER: Xie Y
PROVIDER: S-EPMC7793804 | biostudies-literature | 2020 Jan
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Clonal bacterial populations rely on transcriptional variation across individual cells to produce specialized states that increase fitness. Understanding all cell states requires studying isogenic bacterial populations at the single-cell level. Here we developed probe-based bacterial sequencing (ProBac-seq), a method that uses libraries of DNA probes and an existing commercial microfluidic platform to conduct bacterial single-cell RNA sequencing. We sequenced the transcriptome of thousands of individual bacterial cells per experiment, detecting several hundred transcripts per cell on average. Applied to Bacillus subtilis and Escherichia coli, ProBac-seq correctly identifies known cell states and uncovers previously unreported transcriptional heterogeneity. In the context of bacterial pathogenesis, application of the approach to Clostridium perfringens reveals heterogeneous expression of toxin by a subpopulation that can be controlled by acetate, a short-chain fatty acid highly prevalent in the gut. Overall, ProBac-seq can be used to uncover heterogeneity in isogenic microbial populations and identify perturbations that affect pathogenicity.
| S-EPMC10159851 | biostudies-literature
Project description:Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.
| S-EPMC11222150 | biostudies-literature
Project description:Transcription factors (TFs) exert their regulatory influence through the binding of enhancers, resulting in coordination of gene expression programs. Active enhancers are often characterized by the presence of short, unstable transcripts termed enhancer RNAs (eRNAs). While their function remains unclear, we demonstrate that eRNAs are a powerful readout of TF activity. We infer sites of eRNA origination across hundreds of publicly available nascent transcription data sets and show that eRNAs initiate from sites of TF binding. By quantifying the colocalization of TF binding motif instances and eRNA origins, we derive a simple statistic capable of inferring TF activity. In doing so, we uncover dozens of previously unexplored links between diverse stimuli and the TFs they affect.
| S-EPMC5848612 | biostudies-literature
Project description:Biomedical research often involves conducting experiments on model organisms in the anticipation that the biology learnt will transfer to humans. Previous comparative studies of mouse and human tissues were limited by the use of bulk-cell material. Here we show that transfer learning-the branch of machine learning that concerns passing information from one domain to another-can be used to efficiently map bone marrow biology between species, using data obtained from single-cell RNA sequencing. We first trained a multiclass logistic regression model to recognize different cell types in mouse bone marrow achieving equivalent performance to more complex artificial neural networks. Furthermore, it was able to identify individual human bone marrow cells with 83% overall accuracy. However, some human cell types were not easily identified, indicating important differences in biology. When re-training the mouse classifier using data from human, less than 10 human cells of a given type were needed to accurately learn its representation. In some cases, human cell identities could be inferred directly from the mouse classifier via zero-shot learning. These results show how simple machine learning models can be used to reconstruct complex biology from limited data, with broad implications for biomedical research.
| S-EPMC7718277 | biostudies-literature
Project description:BackgroundMicroglia are innate immune cells that are the only residential macrophages in the central nervous system. They play vital physiological roles in the adult brain and during development. Microglia are particularly in the spotlight because many genetic risk factors recently identified for neurodegenerative diseases are largely expressed in microglia. Rare polymorphisms in these risk alleles lead to abnormal activity of microglia under traumatic or disease conditions.MethodsIn the present study, to investigate the multifaceted functions of human microglia, we established a novel robust protocol to generate microglia from human induced pluripotent stem cells (hiPSCs) using a combination of cytokines and small chemicals essential for microglia ontogeny. Moreover, we highly enhanced the microglial differentiation efficiency by forcing the expression of PU.1, a crucial transcription factor for microglial development, during posterior mesoderm differentiation.ResultsBy our novel method, we demonstrated the generation of a greater number of hiPSC-derived microglia (hiMGLs, approximately 120-folds) than the prior methods (at most 40-folds). Over 90% of the hiMGLs expressed microglia-specific markers, such as CX3CR1 and IBA-1. Whole-transcriptome analysis revealed that these hiMGLs are similar to human primary microglia but differ from monocytes/macrophages. Furthermore, the specific physiological functions of microglia were confirmed through indices of lipopolysaccharide responsiveness, phagocytotic ability, and inflammasome formation. By co-culturing these hiMGLs with mouse primary neurons, we demonstrated that hiMGLs can regulate the activity and maturation of neurons.ConclusionsIn this study, our new simple, rapid, and highly efficient method for generating microglia from hiPSCs will prove useful for future investigations on microglia in both physiological and disease conditions, as well as for drug discovery.
| S-EPMC9248164 | biostudies-literature
Project description:Background: Osteosarcoma (OS) is a common primary tumor with extensive heterogeneity. In this study, we used single-cell RNA sequencing (scRNA-seq) and network pharmacology to analyze effective targets for Osteosarcoma treatment. Methods: The cell heterogeneity of the Osteosarcoma single-cell dataset GSE162454 was analyzed using the Seurat package. The bulk-RNA transcriptome dataset GSE36001 was downloaded and analyzed using the CIBERSORT algorithm. The key targets for OS therapy were determined using Pearson's correlation analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on key targets. The DeepDR algorithm was used to predict potential drugs for Osteosarcoma treatment. Molecular docking analysis was performed to verify the binding abilities of the predicted drugs and key targets. qRT-PCR assay was used to detect the expression of key targets in osteoblasts and OS cells. Results: A total of 21 cell clusters were obtained based on the GSE162454 dataset, which were labeled as eight cell types by marker gene tagging. Four cell types (B cells, cancer-associated fibroblasts (CAFs), endothelial cells, and plasmocytes) were identified in Osteosarcoma and normal tissues, based on differences in cell abundance. In total, 17 key targets were identified by Pearson's correlation analysis. GO and KEGG analysis showed that these 17 genes were associated with immune regulation pathways. Molecular docking analysis showed that RUNX2, OMD, and CD4 all bound well to vincristine, dexamethasone, and vinblastine. The expression of CD4, OMD, and JUN was decreased in Osteosarcoma cells compared with osteoblasts, whereas RUNX2 and COL9A3 expression was increased. Conclusion: We identified five key targets (CD4, RUNX2, OMD, COL9A3, and JUN) that are associated with Osteosarcoma progression. Vincristine, dexamethasone, and vinblastine may form a promising drug-target pair with RUNX2, OMD, and CD4 for Osteosarcoma treatment.
| S-EPMC9853455 | biostudies-literature
Project description:Single-cell RNA sequencing (scRNA-seq) is a widely used method for identifying cell types and trajectories in biologically heterogeneous samples, but it is limited in its detection and quantification of lowly expressed genes. This results in missing important biological signals, such as the expression of key transcription factors (TFs) driving cellular differentiation. We show that targeted sequencing of ∼1000 TFs (scCapture-seq) in iPSC-derived neuronal cultures greatly improves the biological information garnered from scRNA-seq. Increased TF resolution enhanced cell type identification, developmental trajectories, and gene regulatory networks. This allowed us to resolve differences among neuronal populations, which were generated in two different laboratories using the same differentiation protocol. ScCapture-seq improved TF-gene regulatory network inference and thus identified divergent patterns of neurogenesis into either excitatory cortical neurons or inhibitory interneurons. Furthermore, scCapture-seq revealed a role for of retinoic acid signaling in the developmental divergence between these different neuronal populations. Our results show that TF targeting improves the characterization of human cellular models and allows identification of the essential differences between cellular populations, which would otherwise be missed in traditional scRNA-seq. scCapture-seq TF targeting represents a cost-effective enhancement of scRNA-seq, which could be broadly applied to improve scRNA-seq resolution.
| S-EPMC8168586 | biostudies-literature
Project description:High-occupancy target (HOT) regions are compact genome loci occupied by many different transcription factors (TFs). HOT regions were initially defined in invertebrate model organisms, and we here show that they are a ubiquitous feature of the human gene-regulation landscape.We identified HOT regions by a comprehensive analysis of ChIP-seq data from 96 DNA-associated proteins in 5 human cell lines. Most HOT regions co-localize with RNA polymerase II binding sites, but many are not near the promoters of annotated genes. At HOT promoters, TF occupancy is strongly predictive of transcription preinitiation complex recruitment and moderately predictive of initiating Pol II recruitment, but only weakly predictive of elongating Pol II and RNA transcript abundance. TF occupancy varies quantitatively within human HOT regions; we used this variation to discover novel associations between TFs. The sequence motif associated with any given TF's direct DNA binding is somewhat predictive of its empirical occupancy, but a great deal of occupancy occurs at sites without the TF's motif, implying indirect recruitment by another TF whose motif is present.Mammalian HOT regions are regulatory hubs that integrate the signals from diverse regulatory pathways to quantitatively tune the promoter for RNA polymerase II recruitment.
| S-EPMC3826616 | biostudies-literature
Project description:BackgroundIt has been reported that tumor immune microenvironment performs a vital role in tumor progress. However, acting mechanism of immune cell related genes (IRGs) in esophageal squamous cell carcinoma (ESCC) is uncertain.MethodsTCGA-ESCC, GSE23400, GSE26886, GSE75241, and GSE196756 datasets were gained via public databases. First, differentially expressed genes (DEGs) between ESCC and control samples from GSE23400, GSE26886, and GSE75241 were screened out by differential expression analysis, and overlapping DEGs were identified. Single-cell transcriptome data of GSE196756 were applied to explore immune cells that might be involved in regulation of ESCC. Then, weighted gene co-expression network analysis was applied to screen IRGs. Next, differentially expressed IRGs (DE-IRGs) were identified by overlapping IRGs and DEGs, and were incorporated into univariate Cox, least absolute shrinkage and selection operator, and multivariate Cox to acquire prognosis-related genes, and ESCC samples were grouped into high-/low-risk groups on the basis of median risk score. Finally, the role of prognosis model in immunotherapy was analyzed.ResultsTotally 248 DEGs were yielded by overlapping 3,915 DEGs in GSE26886, 459 DEGs in GSE23400, and 1,641 DEGs in GSE75241. Single-cell analysis found that B cells, dendritic cells, monocytes, neutrophils, natural killer cells, and T cells were involved in ESCC development. Besides, MEred, MEblack, MEpink, MEblue and MEbrown modules were considered as key modules because of their highest correlations with immune cell subtypes. A total of 154 DE-IRGs were yielded by taking intersection of DEGs and genes in key modules. Moreover, CTSC, ALOX12, and RMND5B were identified as prognosis-related genes in ESCC. Obviously, Exclusion and TIDE scores were notably lower in high-risk group than in the other one, indicating that high-risk group was more responsive to immunotherapy.ConclusionsThrough bioinformatic analysis, we identified a prognosis model consisting of IRGs (CTSC, ALOX12, and RMND5B) in ESCC, providing new ideas for studies related to treatment and prognosis of ESCC.
| S-EPMC11187079 | biostudies-literature
Project description:The incidence rate and mortality rate of colorectal cancer (CRC) ranks third and second globally. Cancer-associated fibroblasts (CAFs) are the major constituent of the stromal cells in the tumor microenvironment (TME) and are closely associated with patients' prognoses. Our study intended to establish a prognostic model for CRC using hallmark genes of CAFs. The expression values of genes and clinicopathological characteristics of patients were enrolled from the cancer genome atlas database as well as the gene expression omnibus database. The single-cell RNA sequencing data were collected and analyzed in the deeply integrated human single-cell omics database and cancer single-cell expression map databases. The ESTIMATE algorithm was applied to access the infiltration levels of immune and stromal cells. The prognostic genes were selected by the Cox regression analysis and the prognostic signature was constructed by the least absolute shrinkage and selection operator algorithm. gene set enrichment analysis was used to explore the enriched gene sets. In this study, based on bulk RNA sequencing and single-cell RNA sequencing data, and we found that more CAFs were infiltrated in the tumor microenvironment and consisted of 3 subtypes. Then we constructed a prognostic signature for CRC using hallmark genes of CAFs and proved that this signature exhibited high values to predict the overall survival of CRC patients in independent training and validating cohorts. Besides, function enrichment analysis revealed that our prognostic model was significantly associated with immune regulation. Further analysis showed that the infiltrated levels of tumor-suppressing immune cells and the expression of higher immune checkpoint genes in CRC tissues were higher in patients with high-risk scores. Furthermore, immunohistochemistry analysis exhibited that these genes in our prognostic signature were markedly upregulated in CRC tissues. We first constructed a signature based on CAFs hallmark genes to predict the survival of CRC patients and further revealed that the tumor-suppressing microenvironment and dysregulated immune checkpoint genes in CRC tissues were partly responsible for the poor prognosis of patients.
| S-EPMC9997814 | biostudies-literature