Браузер не даёт сохранять данные (например, приватный режим): закладки пропадут после закрытия вкладки.
Поиск статей PubMed
Ключевые слова, автор, журнал, годы — результаты из PubMed, аннотации с переводом на русский, закладки в браузере.
Найдено: 4 928
Страница 235 из 493-
scTWAS Atlas: an integrative knowledgebase of single-cell transcriptome-wide association studies.
Single-cell transcriptome-wide association studies (scTWAS) is a new method for conducting TWAS analysis at the cellular level to identify gene-trait associations with higher precision. This approach helps overcome the challenge of interpreting cell-type heterogeneity in traditional TWAS results. As the field of scTWAS rapidly advances, there is a growing need for additional database platforms to integrate this wealth of data and knowledge effectively. To address this gap, we present scTWAS Atlas (https://ngdc.cncb.ac.cn/sctwas/), a comprehensive database of scTWAS information integrating literature curation and data analysis. The current version of scTWAS Atlas amasses 2,765,211 associations encompassing 34 traits, 30 cell types, 9 cell conditions and 16,470 genes. The database features visualization tools, including an interactive knowledge graph that integrates single-cell expression quantitative trait loci (sc-eQTL) and scTWAS associations to build a multi-omics level regulatory network at the cellular level. Additionally, scTWAS Atlas facilitates cross-cell-type analysis, highlighting cell-type-specific and shared TWAS genes. The database is designed with user-friendly interfaces and allows for easy browsing, searching, and downloading of relevant information. Overall, scTWAS Atlas is instrumental in exploring the genetic regulatory mechanisms at the cellular level and shedding light on the role of various cell types in biological processes, offering novel insights for human health research.
-
Gene-regulatory programs that specify age-related differences during thymocyte development.
T cell development is fundamental to immune system establishment, but how this development changes with age remains poorly understood. Here, we construct a transcriptional and chromatin accessibility atlas of T cell developmental programs in neonatal and adult mice, revealing the ontogeny of divergent gene-regulatory programs and their link to age-related differences. Specifically, we identify a gene module that diverges with age from the earliest stages of genesis and includes programs that govern the effector response and cell cycle. Moreover, we reveal that neonates possess more accessible chromatin during early thymocyte development, likely establishing poised gene expression programs that manifest later in thymocyte development. Finally, we leverage this atlas, employing a CRISPR-based perturbation approach coupled with single-cell RNA sequencing readout, to uncover a conserved transcriptional regulator, Zbtb20, that contributes to age-dependent differences in T cell development. In summary, our study defines gene-regulatory programs that regulate age-specific differences in T cell development.
-
Inferring gene regulatory networks from single-cell multiome data using atlas-scale external data.
Existing methods for gene regulatory network (GRN) inference rely on gene expression data alone or on lower resolution bulk data. Despite the recent integration of chromatin accessibility and RNA sequencing data, learning complex mechanisms from limited independent data points still presents a daunting challenge. Here we present LINGER (Lifelong neural network for gene regulation), a machine-learning method to infer GRNs from single-cell paired gene expression and chromatin accessibility data. LINGER incorporates atlas-scale external bulk data across diverse cellular contexts and prior knowledge of transcription factor motifs as a manifold regularization. LINGER achieves a fourfold to sevenfold relative increase in accuracy over existing methods and reveals a complex regulatory landscape of genome-wide association studies, enabling enhanced interpretation of disease-associated variants and genes. Following the GRN inference from reference single-cell multiome data, LINGER enables the estimation of transcription factor activity solely from bulk or single-cell gene expression data, leveraging the abundance of available gene expression data to identify driver regulators from case-control studies.
-
SCMO: a deep learning model integrating the single-cell resolution TME ecosystem and multi-omics for survival prediction in CRC patients.
BACKGROUND: Colorectal cancer (CRC) remains a leading cause of global cancer mortality, highlighting the need for precise survival prediction to guide clinical decisions. Although tissue-level multi-omics is widely utilized for survival prediction, its limited resolution cannot capture tumor heterogeneity. Single-cell RNA sequencing (scRNA-seq) enables dissection of the tumor microenvironment (TME) at cellular resolution, supporting personalized prognostic assessment. METHODS: We collected 213 CRC scRNA-seq samples and established a CRC-specific TME atlas comprising 339,060 cells. Using this atlas as a reference, we deconvolved bulk RNA-seq data from TCGA-CRC cohort with the EcoTyper algorithm to reconstruct TME features. Clinical, genomic, and transcriptomic data were obtained from the Xena platform; microbial data were sourced from the BIC database. We integrated TME and multi-omics features through a self-normalizing neural network to construct a deep learning model (single-cell resolution TME ecosystem with multi-omics data [SCMO]) for survival prediction. To enhance interpretability, we utilized the Integrated Gradients algorithm and spatial transcriptomic data to analyze multi-omics and TME features. We performed anticancer drug screening with tumor necrosis factor receptor-associated protein 1 (TRAP1), a critical feature according to the Integrated Gradients algorithm, as a potential target. RESULTS: We identified 13 survival-related TME features from the CRC-specific atlas: 12 cell states and one multi-cellular ecosystem. SCMO, which combined TME and multi-omics features, improved survival prediction and outperformed existing methods, achieving a concordance index of 0.762. The SCMO demonstrated robust performance for long-term predictions, achieving areas under the curve (AUCs) of 0.752, 0.772, and 0.869 for 1-, 3-, and 5-year predictions in the training set, with corresponding test set AUCs of 0.639, 0.756, and 0.772. TME features from the SCMO model revealed that ecosystem density increased with CRC malignancy. Multi-omics features included TRAP1 as a potential drug target. Drug screening identified saikosaponin A as a novel TRAP1 inhibitor, and its anticancer activity was validated in vitro. We developed SCMO-Lite, a simplified model incorporating 12 high-attribution-weight multi-omics features, which demonstrated robust risk stratification. CONCLUSIONS: SCMO combines analytical precision with biological interpretability, offering novel insights for oncology survival prediction.
-
Single-Cell Analysis of Regional Differences in Adult V-SVZ Neural Stem Cell Lineages.
The ventricular-subventricular zone (V-SVZ) harbors adult neural stem cells. V-SVZ neural stem cells exhibit features of astrocytes, have a regional identity, and depending on their location in the lateral or septal wall of the lateral ventricle, generate different types of neuronal and glial progeny. We performed large-scale single-cell RNA sequencing to provide a molecular atlas of cells from the lateral and septal adult V-SVZ of male and female mice. This revealed regional and sex differences among adult V-SVZ cells. We uncovered lineage potency bias at the single-cell level among lateral and septal wall astrocytes toward neurogenesis and oligodendrogenesis, respectively. Finally, we identified transcription factor co-expression modules marking key temporal steps in neurogenic and oligodendrocyte lineage progression. Our data suggest functionally important spatial diversity in neurogenesis and oligodendrogenesis in the adult brain and reveal molecular correlates of adult NSC dormancy and lineage specialization.
-
A Vascular-Extracellular Matrix Molecular Program Identifies High-Risk Diffuse Glioma Across Independent Multi-Omics.
BACKGROUND: Gliomas are characterized by a high degree of molecular heterogeneity, which impairs the reproducibility of predictive biomarkers derived from bulk-based molecular profiling due to immune/stromal contamination of tumors and the high prevalence of the IDH mutation signature. METHODS: In this study, we used MOFA+ to derive intrinsic molecular signatures from transcriptional, methylation, and genomic profiles of a cohort of 667 diffuse gliomas in the Cancer Genome Atlas database. Thereafter, factor scores were derived for two separate Chinese Glioma Genome Atlas batches (Batch 1, n = 325; Batch 2, n = 693) without any retraining on the model. The prognostic independence of identified molecular signatures was assessed using multivariable Cox regression adjusted for IDH mutation status and tumor purity; purity-residualized survival analyses; IDH-stratified Cox regression in each cohort; validation by concordance index against established molecular signatures; and survival extreme profiling. To characterize the biological significance of factor signatures, we projected gene set signatures corresponding to each factor signature onto a single-cell RNA-seq dataset of GBM (GSE131928). RESULTS: MOFA+ identified 12 latent factors, of which a vascular-extracellular matrix (ECM) remodeling axis (Factor 1) explained the highest multi-omics variance (24.9%) and was the strongest independent prognostic factor. In multivariable Cox regression adjusting for IDH status and tumor purity, Factor 1 remained independently prognostic (HR = 1.67, 95% CI 1.27-2.20, p = 0.0002); in a fully-adjusted model additionally including age, WHO grade, MGMT methylation, and 1p/19q codeletion (plus radiotherapy and chemotherapy status in the CGGA cohorts), Factor 1 remained prognostic in both CGGA cohorts (CGGA1: HR = 1.50, p = 3.8 × 10-5; CGGA2: HR = 1.18, p = 0.003) but lost significance in TCGA (HR = 1.04, p = 0.83), consistent with the cohort-dependent magnitude reported in the IDH-stratified and meta-regression analyses below. Purity-residualized survival analysis showed negligible attenuation of the Factor 1 signal (raw HR = 3.57 vs. residualized HR = 3.72; concordance 96.5%). Within IDH-wildtype gliomas, Factor 1 was significant in both external validation cohorts (CGGA1: HR = 1.64, FDR = 4.6 × 10-6; CGGA2: HR = 1.20, FDR = 0.02), though the TCGA IDH-wildtype subgroup showed a trend that did not survive FDR correction (FDR = 0.060). All validation was performed without model retraining. Within IDH-mutant gliomas, Factor 1 was strongly prognostic in both CGGA cohorts but was not significant in TCGA (HR = 1.17, FDR = 0.33). These findings should therefore be interpreted as consistent in directionality across cohorts but not uniformly replicated at the FDR-adjusted significance threshold in the TCGA discovery dataset. Concordance index benchmarking on a matched subset (n = 503) showed Factor 1 achieved discrimination comparable to the Mesenchymal signature (C = 0.797 vs. 0.801; ΔC = -0.004) while outperforming four other established classifiers. Factor 1 consistently separated patients with extreme survival phenotypes (OS < 6 vs. >15 months) across all three cohorts (all log-rank p < 0.001). Projection onto a single-cell GBM atlas (GSE131928), supported by inferCNV-based malignant-cell classification, localized the Vascular-ECM program to malignant cells and the Immune-ECM axis to myeloid compartments. CONCLUSIONS: The Vascular-ECM axis is a consistent, prognostic program robust to purity adjustment for diffuse gliomas that remains relevant across IDH-defined subgroups in three independent datasets comprising 1685 patients. The Vascular-ECM axis is a reproducible, purity-robust prognostic program in diffuse glioma, with directionally consistent adverse effects across TCGA, CGGA Batch 1, and CGGA Batch 2 (pooled n = 1685). Given the strong co-loading of endothelial, ECM, and myeloid genes observed in the single-cell projection, Factor 1 is best interpreted as a vascular/ECM-associated tumor-microenvironment ecosystem program rather than a malignant-cell-autonomous signature. Its FDR-adjusted significance within IDH-stratified subgroups is cohort-dependent and robust in both CGGA cohorts but attenuated in the TCGA IDH-wildtype (FDR = 0.060) and TCGA IDH-mutant (FDR = 0.33) strata. The pooled signal should therefore be interpreted as evidence of a generalizable biological program rather than a uniformly replicated subgroup-specific biomarker. It is possible to calculate factor scores based on RNA sequencing alone using fixed loadings (Z = XWᵀ), which may have implications for future translational applications. All findings are correlative; a causal role for the Vascular-ECM program in glioma progression, invasion, or therapy resistance remains to be established through functional perturbation experiments.
-
Processing and cryopreservation of human ureter tissues for single-cell and spatial transcriptomics assays.
Characterizing the cellular heterogeneity of human ureter tissues using single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics provides a detailed atlas of cell types, signaling networks, and potential cell-cell cross talk underlying developmental and regenerative pathways. We describe an optimized protocol for generating, cryopreserving, and thawing single-cell suspensions from ureter tissues isolated post-cystectomy for scRNA-seq. In addition, we describe an optimized protocol for cryopreserving human ureter tissues for 10x Genomics Visium spatial gene expression platform. For complete details on the use and execution of this protocol, please refer to Fink et al. (2022).1.
-
Critical Differential Expression Assessment for Individual Bulk RNA-Seq Projects.
Finding the right balance of quality and quantity can be important, and it is essential that project quality does not drop below the level where important main conclusions are missed or misstated. We use knock-out and over-expression studies as a simplification to test recovery of a known causal gene in RNA-Seq cell line experiments. When single-end RNA-Seq reads are aligned with STAR and quantified with htseq-count, we found potential value in testing the use of the Generalized Linear Model (GLM) implementation of edgeR with robust dispersion estimation more frequently for either single-variate or multi-variate 2-group comparisons (with the possibility of defining criteria less stringent than |fold-change| > 1.5 and FDR < 0.05). When considering a limited number of patient sample comparisons with larger sample size, there might be some decreased variability between methods (except for DESeq1). However, at the same time, the ranking of the gene identified using immunohistochemistry (for ER/PR/HER2 in breast cancer samples from The Cancer Genome Atlas) showed as possible shift in performance compared to the cell line comparisons, potentially highlighting utility for standard statistical tests and/or limma-based analysis with larger sample sizes. If this continues to be true in additional studies and comparisons, then that could be consistent with the possibility that it may be important to allocate time for potential methods troubleshooting for genomics projects. Analysis of public data presented in this study does not consider all experimental designs, and presentation of downstream analysis is limited. So, any estimate from this simplification would be an underestimation of the true need for some methods testing for every project. Additionally, this set of independent cell line experiments has a limitation in being able to determine the frequency of missing a highly important gene if the problem is rare (such as 10% or lower). For example, if there was an assumption that only one method can be tested for "initial" analysis, then it is not completely clear to the extent that using edgeR-robust might perform better than DESeq2 in the cell line experiments. Importantly, we do not wish to cause undue concern, and we believe that it should often be possible to define a gene expression differential expression workflow that is suitable for some purposes for many samples. Nevertheless, at the same time, we provide a variety of measures that we believe emphasize the need to critically assess every individual project and maximize confidence in published results.
-
A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis.
Complex biological processes such as embryogenesis require precise coordination of cell differentiation programs across both space and time. Using protein-fusion fluorescent reporters and four-dimensional live imaging, we present a protein expression atlas of transcription factors (TFs) mapped onto developmental cell lineages during Caenorhabditis elegans embryogenesis, at single-cell resolution. This atlas reveals a spatiotemporal combinatorial code of TF expression, and a cascade of lineage-specific, tissue-specific and time-specific TFs that specify developmental states. The atlas uncovers regulators of embryogenesis, including an unexpected role of a skin specifier in neurogenesis and the critical function of an uncharacterized TF in convergent muscle differentiation. At the systems level, the atlas provides an opportunity to model cell state-fate relationships, revealing a lineage-dependent state diversity within functionally related cells and a winding trajectory of developmental state progression. Collectively, this single-cell protein atlas represents a valuable resource for elucidating metazoan embryogenesis at the molecular and systems levels.
-
C5aR1 blockade reshapes immunosuppressive tumor microenvironment and synergizes with immune checkpoint blockade therapy in high-grade serous ovarian cancer.
High-grade serous ovarian cancer (HGSC), with a modest response to immune checkpoint blockade (ICB) targeting PD-1/PD-L1 monotherapy, is densely infiltrated by M2-polarized tumor-associated macrophages (TAMs) and regulatory T (Treg) cells. The complement C5a/C5aR1 axis contributes to the programming of the immunosuppressive phenotype of TAMs in solid tumors and represents a promising immunomodulatory target for treating HGSCs. Here, we aimed to identify the relevance of C5aR1 in prognosis, immune microenvironment, and immunotherapy response in HGSCs. The expression and relationship of C5aR1 with tumor-infiltrating immune cells were assessed by immunohistochemistry and flow cytometry in the training cohort (n = 120) and fresh HGSC tissues (n = 36). Transcriptomic analyses of the xenografts delineated the mechanisms driving the immunomodulatory activity of PMX53, an orally bioavailable C5aR1 inhibitor. Therapeutic relevance was confirmed in ex vivo tumor cultures and The Cancer Genome Atlas (TCGA) datasets. C5aR1 expression independently predicted dismal prognosis and was linked to the immunoevasive subtype of HGSC, characterized by increased infiltration of pro-tumor cells (Treg cells, M2-polarized macrophages, and neutrophils) and impaired CD8+T functions. PMX53 antagonized subcutaneous tumor growth, modulated immunosuppressive mechanisms and synergized with aPD-1 in several tumor types. Single-cell RNA-seq analysis revealed predominant C5aR1 expression in TAMs, with an immunosuppressive-related expression signature in C5aR1+TAMs. Furthermore, the combination of C5aR1 and PD-L1 was associated with specific molecular characteristics and matched clinical response annotations. Therefore, the abundance of C5aR1 could predict an inferior prognosis in HGSCs, and incorporating PD-L1 may serve as a novel predictive biomarker to guide therapeutic options.