Single-cell RNA Research Articles

Loss of Stim2 in zebrafish induces glaucoma-like phenotype.

Calcium is involved in vision processes in the retina and implicated in various pathologies, including glaucoma. Rod cells rely on store-operated calcium entry (SOCE) to safeguard against the prolonged lowering of intracellular calcium ion concentrations. Zebrafish that lacked the endoplasmic reticulum Ca2+ sensor Stim2 (stim2 knockout [KO]) exhibited impaired vision and lower light perception-related gene expression. We sought to understand mechanisms that are responsible for vision impairment in stim2 KO zebrafish. The single-cell RNA (scRNA) sequencing of neuronal cells from brains of 5days postfertilization larvae distinguished 27 cell clusters, 10 of which exhibited distinct gene expression patterns, including amacrine and γ-aminobutyric acid (GABA)ergic retinal interneurons and GABAergic optic tectum cells. Five clusters exhibited significant changes in cell proportions between stim2 KO and controls, including GABAergic diencephalon and optic tectum cells. Transmission electron microscopy of stim2 KO zebrafish revealed decreases in width of the inner plexiform layer, ganglion cells, and their dendrites numbers (a hallmark of glaucoma). GABAergic neuron densities in the inner nuclear layer, including amacrine cells, as well as photoreceptors significantly decreased in stim2 KO zebrafish. Our study suggests a novel role for Stim2 in the regulation of neuronal insulin expression and GABAergic-dependent vision causing glaucoma-like retinal pathology.

RNF144A as a potential risk gene for endometrial carcinoma: Insights from Mendelian randomization, bulk RNA sequencing, single-cell RNA, and experimental analysis.

Endometrial carcinoma (EC) is a prevalent gynecological malignancy that poses a significant threat to women's health worldwide. However, its pathogenesis and underlying mechanisms remains unclear. In this study, expression quantitative trait loci data, Mendelian randomization analysis, and differential expression analysis were performed to identify potential targets. A prognostic risk signature was subsequently constructed for EC patients based on the expression of these genes. Four bioinformatics algorithms, including generalized linear model, extreme gradient boosting, support vector machine, and random forest, were used to identify hub genes in EC. The expression of ring finger protein 144A (RNF144A) was validated using quantitative real-time polymerase chain reaction. Cellular proliferation and migration ability were evaluated using CCK-8 and Transwell assays, respectively. The genes RNF144A, ketohexokinase, and Rab interacting lysosomal protein like 2 were identified as potential targets for EC. Their differential expression was observed in EC patients, and Mendelian randomization analysis revealed a negative correlation between these genes and the development of EC. Mechanistic analyses suggested a strong association between these genes and the tumor immune microenvironment. The constructed risk signature was significantly associated with the prognosis, age, cancer stage, and grade of EC patients. Furthermore, based on interacted model algorithms, RNF144A was identified as a hub gene in EC. It was found to be significantly downregulated in EC samples, and its expression was positively correlated with the stage and grade of EC patients. In vitro experiments showed that overexpression of RNF144A significantly promoted cell growth and migration in EC cells. In conclusion, this study provides insights into the molecular mechanisms underlying EC progression and identifies preliminary candidate biomarkers for the development of EC treatment strategies.

A random survival forest-based pathomics signature classifies immunotherapy prognosis and profiles TIME and genomics in ES-SCLC patients

BackgroundSmall cell lung cancer (SCLC) is a highly aggressive neuroendocrine tumor with high mortality, and only a limited subset of extensive-stage SCLC (ES-SCLC) patients demonstrate prolonged survival under chemoimmunotherapy, which warrants the exploration of reliable biomarkers. Herein, we built a machine learning-based model using pathomics features extracted from hematoxylin and eosin (H&E)-stained images to classify prognosis and explore its potential association with genomics and TIME.MethodsWe retrospectively recruited ES-SCLC patients receiving first-line chemoimmunotherapy at Nanjing Jinling Hospital between April 2020 and August 2023. Digital H&E-stained whole-slide images were acquired, and targeted next-generation sequencing, programmed death ligand-1 staining, and multiplex immunohistochemical staining for immune cells were performed on a subset of patients. A random survival forest (RSF) model encompassing clinical and pathomics features was established to predict overall survival. The function of putative genes was assessed via single-cell RNA sequencing.Results and conclusionDuring the median follow-up period of 12.12 months, 118 ES-SCLC patients receiving first-line immunotherapy were recruited. The RSF model utilizing three pathomics features and liver metastases, bone metastases, smoking status, and lactate dehydrogenase, could predict the survival of first-line chemoimmunotherapy in patients with ES-SCLC with favorable discrimination and calibration. Underlyingly, the higher RSF-Score potentially indicated more infiltration of CD8+ T cells in the stroma as well as a greater probability of MCL-1 amplification and EP300 mutation. At the single-cell level, MCL-1 was associated with TNFA-NFKB signaling and apoptosis-related processes. Hopefully, this noninvasive model could act as a biomarker for immunotherapy, potentially facilitating precision medicine in the management of ES-SCLC.

BZW2 is a potential regulator of non-small cell lung cancer progression

BackgroundPersonalized targeted therapy has become an important strategy for cancer treatment owing to its remarkable therapeutic efficacy and safety. However, drug resistance remains the primary cause of treatment failure. Basic leucine zipper and W2 domain 2 (BZW2), which is aberrantly expressed in cancer, has been implicated in tumor progression and may serve as a new therapeutic target. Therefore, the role of BZW2 in non-small cell lung cancer (NSCLC) requires further investigation. MethodsThe expression and genetic alterations of BZW2 in pan-cancers were explored using The Cancer Genome Atlas (TCGA) PanCancer databases. The mRNA and protein levels of BZW2 in patients with NSCLC were verified in our cohort. Functional experiments including CCK8, colony formation, and transwell assays were performed to evaluate the impact of BZW2 on the proliferative, migratory, and invasive capacities of SK-MES-1 cells. Gene Set Enrichment Analysis was used to identify underlying biological processes and pathways. Single-cell RNA (scRNA) sequencing data were employed to investigate the tumor microenvironment of NSCLC and the co-expression of BZW2 and stemness-related genes. ResultsDysregulated BZW2 expression was observed in various malignant tumors. BZW2 expression was found to be significantly elevated in NSCLC. BZW2 depletion inhibited the growth, mobility, and invasive abilities of lung squamous cell carcinoma SK-MES-1 cells. BZW2 may be related to signaling pathways such as nucleotide excision repair, ubiquitin-mediated proteolysis, and the P53 signaling pathway. Biological processes, including translational initiation, tRNA processing, and RNA methylation, were observed to be enriched in the high-BZW2 group. Furthermore, there was a positive correlation between BZW2 and the m6A- and m5C-related genes. scRNA analysis revealed a co-expression relationship between BZW2 and stemness-related genes such as CD44, SOX9, and CD133. ConclusionsElevated BZW2 expression is associated with the proliferation, migration, and invasion of NSCLC, and BZW2 may be a potential therapeutic target for NSCLC.

High-dimensional profiling of immune responses to kidney transplant reveals heterogeneous T helper 1 and B cell effectors associated with rejection

Rejection is a primary cause of allograft dysfunction after kidney transplantation. The diversity of immune subpopulations involved in the different endotypes of rejection remains to be delineated at single-cell resolution. In a cohort of 76 kidney transplant recipients, we conducted high-dimensional immune phenotyping of blood CD4 T and B cells, single-cell RNA and T/B cell receptor sequencing, and plasma cytokine profiling. Phenotypic, transcriptional, and clonal states of CD4T and B cells could significantly distinguish stable allograft states from rejection. Patients undergoing T cell-mediated rejection displayed accumulation of clonally expanded cytotoxic T helper (Th)1 cells and Th17-like cells, associated with predominant naive B cell responses. In contrast, antibody-mediated rejection was characterized by clonal expansion of Th1-polarized T follicular helper cells and effector T-bet+ memory B cells, both of which strongly expressed interleukin 12 and tumor necrosis factor-signaling pathways. Plasma cytokine analysis confirmed mixed Th1/Th17 and Th1/T follicular helper cell-driven inflammatory profiles distinguishing T cell-mediated rejection and antibody-mediated rejection, respectively. CD4T and B cell subpopulations and signatures were validated using bulk RNA-seq analysis of matched kidney allografts and using an independent single-cell RNA-seq data set. These data improve mechanistic understanding of the immune pathogenesis of rejection and support the development of more specific immunosuppressive therapies to treat allograft rejection.

UCHL1 Overexpression Is Related to the Aggressive Phenotype of Non-small Cell Lung Cancer.

Ubiquitin C-terminal hydrolase L1 (UCHL1), which encodes thiol protease that hydrolyzes a peptide bond at the C-terminal glycine residue of ubiquitin, regulates cell differentiation, proliferation, transcriptional regulation, and numerous other biological processes and may be involved in lung cancer progression. UCHL1 is mainly expressed in the brain and plays a tumor-promoting role in a few cancer types; however, there are limited reports regarding its role in lung cancer. Single-cell RNA (scRNA) sequencing using 10X chromium v3 was performed on a paired normal-appearing and tumor tissue from surgical specimens of a patient who showed unusually rapid progression. To validate clinical implication of the identified biomarkers, immunohistochemical (IHC) analysis was performed on 48 non-small cell lung cancer (NSCLC) tissue specimens, and the correlation with clinical parameters was evaluated. We identified 500 genes overexpressed in tumor tissue compared to those in normal tissue. Among them, UCHL1, brain expressed X-linked 3 (BEX3), and midkine (MDK), which are associated with tumor growth and progression, exhibited a 1.5-fold increase in expression compared to that in normal tissue. IHC analysis of NSCLC tissues showed that only UCHL1 was specifically overexpressed. Additionally, in 48 NSCLC specimens, UCHL1 was specifically upregulated in the cytoplasm and nuclear membrane of tumor cells. Multivariable logistic analysis identified several factors, including smoking, tumor size, and high-grade dysplasia, to be typically associated with UCHL1 overexpression. Survival analyses using The Cancer Genome Atlas (TCGA) datasets revealed that UCHL1 overexpression is substantially associated with poor survival outcomes. Furthermore, a strong association was observed between UCHL1 expression and the clinicopathological features of patients with NSCLC. UCHL1 overexpression was associated with smoking, tumor size, and high-grade dysplasia, which are typically associated with a poor prognosis and survival outcome. These findings suggest that UCHL1 may serve as an effective biomarker of NSCLC.

Treatment resistance to melanoma therapeutics on a single cell level

Therapy targeting the BRAF-MEK cascade created a treatment revolution for patients with BRAF mutant advanced melanoma. Unfortunately, 80% patients treated will progress by 5 years follow-up. Thus, it is imperative we study mechanisms of melanoma progression and therapeutic resistance. We created a scRNA (single cell RNA) atlas of 128,230 cells from 18 tumors across the treatment spectrum, discovering melanoma cells clustered strongly by transcriptome profiles of patients of origins. Our cell-level investigation revealed gains of 1q and 7q as likely early clonal events in metastatic melanomas. By comparing patient tumors and their derivative cell lines, we observed that PD1 responsive tumor fraction disappears when cells are propagated in vitro. We further established three anti-BRAF-MEK treatment resistant cell lines using three BRAF mutant tumors. ALDOA and PGK1 were found to be highly expressed in treatment resistant cell populations and metformin was effective in targeting the resistant cells. Our study suggests that the investigation of patient tumors and their derivative lines is essential for understanding disease progression, treatment response and resistance.

Tumor-reactive T cell clonotype dynamics underlying clinical response to TIL therapy in melanoma

Adoptive cell therapy (ACT) using invitro expanded tumor-infiltrating lymphocytes (TILs) has inconsistent clinical responses. To better understand determinants of therapeutic success, we tracked TIL clonotypes from baseline tumors to ACT products and post-ACT blood and tumor samples in melanoma patients using single-cell RNA and Tcell receptor (TCR) sequencing. Patients with clinical responses had baseline tumors enriched in tumor-reactive TILs, and these were more effectively mobilized upon invitro expansion, yielding products enriched in tumor-specific CD8+ cells that preferentially infiltrated tumors post-ACT. Conversely, lack of clinical responses was associated with tumors devoid of tumor-reactive resident clonotypes and with cell products mostly composed of blood-borne clonotypes that persisted in blood but not in tumors post-ACT. Upon expansion, tumor-specific TILs lost tumor-associated transcriptional signatures, including exhaustion, and responders exhibited an intermediate exhausted effector state after TIL engraftment in the tumor, suggesting functional reinvigoration. Our findings provide insight into the nature and dynamics of tumor-specific clonotypes associated with clinical response to TIL-ACT, with implications for treatment optimization.

Inhibitory neuron links the causal relationship from air pollution to psychiatric disorders: a large multi-omics analysis

Psychiatric disorders are severe health challenges that exert a heavy public burden. Air pollution has been widely reported as related to psychiatric disorder risk, but their casual association and pathological mechanism remained unclear. Herein, we systematically investigated the large genome-wide association studies (6 cohorts with 1,357,645 samples), single-cell RNA (26 samples with 157,488 cells), and bulk-RNAseq (1595 samples) datasets to reveal the genetic causality and biological link between four air pollutants and nine psychiatric disorders. As a result, we identified ten positive genetic correlations between air pollution and psychiatric disorders. Besides, PM2.5 and NO2 presented significant causal effects on schizophrenia risk which was robust with adjustment of potential confounders. Besides, transcriptome-wide association studies identified the shared genes between PM2.5/NO2 and schizophrenia. We then discovered a schizophrenia-derived inhibitory neuron subtype with highly expressed shared genes and abnormal synaptic and metabolic pathways by scRNA analyses and confirmed their abnormal level and correlations with the shared genes in schizophrenia patients in a large RNA-seq cohort. Comprehensively, we discovered robust genetic causality between PM2.5, NO2, and schizophrenia and identified an abnormal inhibitory neuron subtype that links schizophrenia pathology and PM2.5/NO2 exposure. These discoveries highlight the schizophrenia risk under air pollutants exposure and provide novel mechanical insights into schizophrenia pathology, contributing to pollutant-related schizophrenia risk control and therapeutic strategies development.Graphical

ThymoSpheres culture: A model to study human polyclonal unconventional T cells.

In vitro cultures remain crucial for studying the fundamental mechanisms of human T-cell development. Here, we introduce a novel in vitro cultivation system based on ThymoSpheres (TS): dense spheroids consisting of DLL4-expressing stromal cells and human hematopoietic precursor cells, in the absence of thymic epithelial cells. These spheroids are subsequently cultured at the air-liquid interphase. TS generate large numbers of mature T cells, are easy to manipulate, scalable, and can be repeatably sampled to monitor T-cell differentiation. The mature T cells generated from primary human hematopoietic precursor cells were extensively characterized using single-cell RNA and combined T-cell receptor (TCR) sequencing. These predominantly CD8α T cells exhibit transcriptional and TCR CDR3 characteristics similar to the recently described human polyclonal αβ unconventional T cell (UTC) lineage. This includes the expression of hallmark genes associated with agonist selection, such as IKZF2 (Helios), and the expression of various natural killer receptors. The TCR repertoire of these UTCs is polyclonal and enriched for CDR3-associated autoreactive features and early rearrangements of the TCR-α chain. In conclusion, TS cultures offer an intriguing platform to study the development of this human polyclonal UTC lineage and its inducing selection mechanisms.

Acyl-coenzyme a binding protein (ACBP) - a risk factor for cancer diagnosis and an inhibitor of immunosurveillance

BackgroundThe plasma concentrations of acyl coenzyme A binding protein (ACBP, also known as diazepam-binding inhibitor, DBI, or ‘endozepine’) increase with age and obesity, two parameters that are also amongst the most important risk factors for cancer.MethodsWe measured ACBP/DBI in the plasma from cancer-free individuals, high-risk patients like the carriers of TP53 or BRCA1/2 mutations, and non-syndromic healthy subjects who later developed cancer. In mice, the neutralization of ACBP/DBI was used in models of non-small cell lung cancer (NSCLC) and breast cancer development and as a combination treatment with chemoimmunotherapy (chemotherapy + PD-1 blockade) in the context of NSCLC and sarcomas. The anticancer T cell response upon ACBP/DBI neutralization was characterized by flow cytometry and single-cell RNA sequencing.ResultsCirculating levels of ACBP/DBI were higher in patients with genetic cancer predisposition (BRCA1/2 or TP53 germline mutations) than in matched controls. In non-syndromic cases, high ACBP/DBI levels were predictive of future cancer development, and especially elevated in patients who later developed lung cancer. In preclinical models, ACBP/DBI neutralization slowed down breast cancer and NSCLC development and enhanced the efficacy of chemoimmunotherapy in NSCLC and sarcoma models. When combined with chemoimmunotherapy, the neutralizing monoclonal antibody against ACBP/DBI reduced the frequency of regulatory T cells in the tumor bed, modulated the immune checkpoint profile, and increased activation markers.ConclusionThese findings suggest that ACBP/DBI acts as an endogenous immune suppressor. We conclude that elevation of ACBP/DBI constitutes a risk factor for the development of cancer and that ACBP/DBI is an actionable target for improving cancer immunosurveillance.

Single-cell DNA methylation analysis tool Amethyst reveals distinct noncanonical methylation patterns in human glial cells.

Single-cell sequencing technologies have revolutionized biomedical research by enabling deconvolution of cell type-specific properties in highly heterogeneous tissue. While robust tools have been developed to handle bioinformatic challenges posed by single-cell RNA and ATAC data, options for emergent modalities such as methylation are much more limited, impeding the utility of results. Here we present Amethyst, a comprehensive R package for atlas-scale single-cell methylation sequencing data analysis. Amethyst begins with base-level methylation calls and expedites batch integration, doublet detection, dimensionality reduction, clustering, cell type annotation, differentially methylated region calling, and interpretation of results, facilitating rapid data interaction in a local environment. We introduce the workflow using published single-cell methylation human peripheral blood mononuclear cell (PBMC) and human cortex data. We further leverage Amethyst on an atlas-scale brain dataset to describe a noncanonical methylation pattern in human astrocytes and oligodendrocytes, challenging the notion that this form of methylation is principally relevant to neurons in the brain. Tools such as Amethyst will increase accessibility to single-cell methylation data analysis, catalyzing research progress across diverse contexts.

Single-cell and spatial sequencing identifies senescent and germinal tumor cells in adamantinomatous craniopharyngiomas

Adamantinomatous craniopharyngioma (ACP) is a clinically aggressive tumor without effective treatment method. Previous studies proposed a paracrine tumorigenesis model, in which oncogenic β-catenin induces senescence in pituitary stem cells and the senescent cells lead the formation of paracrine tumors through secretion of pro-tumorigenic factors. However, there lacks characterization on senescent cells in ACPs. Here, we profiled 12 ACPs with single-cell RNA and TCR-sequencing to elucidate the cellular atlas in ACPs and 3 of them were also subject to spatial sequencing to localize different subpopulations of the tumor cells. In total, we obtained the transcriptome profiles of 70,682 cells. Tumor cells, which were unambiguously identified through the cellular mutation status of the driver CTNNB1 mutations, were clustered into 6 subsets. The whorl-like cluster (WC) cells show distinct molecular features from the other tumor cells and the palisading epithelium (PE) cells consists of a proliferating subset. Other than typical PE and WC, we identified two novel subpopulations of the tumor cells. In one subpopulation, the cells express a high level of cytokines, e.g., FDCSP and S100A8/A9, and are enriched with the senescence-associated secretory phenotype (SASP) factors. Hematoxylin and eosin staining reveals that these SASP cells lack an ordered structures and their nuclei are elongated. In the other subpopulation, the cell sizes are small and they are tightly packed together with an unusual high density expressing a high level of mitochondrial genes (median 10.9%). These cells are the origin of the tumor developmental trajectories revealed by RNA velocity and pseudo-time analysis. Single-cell RNA and TCR analysis reveals that some ACPs are infiltrated with clonally expanded cytotoxic T cells. We propose a hypothesis that WC and PE are formed via different negative regulation mechanisms of the overactivated WNT/β-catenin signaling which provides a new understanding on the tumorigenesis of ACPs. The study lays a foundation for future studies on targeting senescent cells in ACPs with senolytic compounds or other therapeutic agents.

Generation of dual-attribute iTNK cells from hPSCs for cancer immunotherapy

Dual-attribute immune cells possess advantageous features of cytotoxic Tcells and natural killer (NK) cells and hold promise for advancing immunotherapy. Dual-attribute cell types such as invariant natural killer Tcells, induced T-to-NK cells, and cytokine-induced killer cells have demonstrated efficacy and safety in preclinical and clinical studies. However, their limited availability hinders their widespread application. Human pluripotent stem cells (hPSCs) offer an ideal source. Here, we generate dual-attribute induced T-NK (iTNK) cells from hPSCs, expressing markers of both cytotoxic T and NK cells. Single-cell RNA and Tcell receptor (TCR) sequencing analyses reveal that iTNK cells expressed signature genes associated with both NK and Tcells and displayed a diverse TCR repertoire. iTNK cells release cytotoxic mediators, exert cytotoxicity against diverse tumor cell lines, and inhibit tumor growth invivo. By harnessing adaptive and innate immune responses, hPSC-derived iTNK cells offer promising strategies for cancer immunotherapy.

Dynamic regulation of tissue fluidity controls skin repair during wound healing

During wound healing, different pools of stem cells (SCs) contribute to skin repair. However, how SCs become activated and drive the tissue remodeling essential for skin repair is still poorly understood. Here, by developing a mouse model allowing lineage tracing and basal cell lineage ablation, we monitor SC fate and tissue dynamics during regeneration using confocal and intravital imaging. Analysis of basal cell rearrangements shows dynamic transitions from a solid-like homeostatic state to a fluid-like state allowing tissue remodeling during repair, as predicted by a minimal mathematical modeling of the spatiotemporal dynamics and fate behavior of basal cells. The basal cell layer progressively returns to a solid-like state with re-epithelialization. Bulk, single-cell RNA, and epigenetic profiling of SCs, together with functional experiments, uncover a common regenerative state regulated by the EGFR/AP1 axis activated during tissue fluidization that is essential for skin SC activation and tissue repair.

Dissecting the Distinct Tumor Microenvironments of HRD and HRP Ovarian Cancer: Implications for Targeted Therapies to Overcome PARPi Resistance in HRD Tumors and Refractoriness in HRP Tumors.

High-grade serous tubo-ovarian cancer (HGSTOC) is an aggressive gynecological malignancy including homologous recombination deficient (HRD) and homologous recombination proficient (HRP) groups. Despite the therapeutic potential of poly (ADP-ribose) polymerase inhibitors (PARPis) and anti-PDCD1 antibodies, acquired resistance in HRD and suboptimal response in HRP patients necessitate more precise treatment. Herein, single-cell RNA and single-cell T-cell receptor sequencing on 5 HRD and 3 HRP tumors are performed to decipher the heterogeneous tumor immune microenvironment (TIME), along with multiplex immunohistochemistry staining and animal experiments for validation. HRD tumors are enriched with immunogenic epithelial cells, FGFR1+PDGFRβ+ myCAFs, M1 macrophages, tumor reactive CD8+/CD4+ Tregs, whereas HRP tumors are enriched with HDAC1-expressing epithelial cells, indolent CAFs, M2 macrophages, and bystander CD4+/CD8+ T cells. Significantly, customized therapies are proposed. For HRD patients, targeting FGFR1+PDGFRβ+ myCAFs via tyrosine kinase inhibitors, targeting Tregs via anti-CCR8 antibodies/TNFRSF4 stimulation, and targeting CXCL13+ exhausted T cells by blocking PDCD1/CTLA-4/LAG-3/TIGIT are proposed. For HRP patients, targeting indolent CAFs, targeting M2 macrophages via CSF-1/CSF-1R inhibitors, targeting bystander T cells via tumor vaccines, and targeting epithelial cells via HDAC inhibitors. The study provides comprehensive insights into HRD and HRP TIME and tailored therapeutic approaches, addressing the challenges of PARPi-resistant HRD and refractory HRP tumors.

Unveiling Cellular Identities and Gene Expression Pathways in Overlapping Myocarditis and Myositis.

Immune checkpoint therapies can drive antitumor responses and benefit patients but can also induce life-threatening immune-related adverse events such as myocarditis and myositis. These immune-related adverse events are rare but carry substantial morbidity and mortality. In this issue, Siddiqui and colleagues use single-cell RNA and T-cell receptor sequencing to identify novel cellular subsets and propose various mechanisms that could contribute to the pathogenesis of immune checkpoint inhibitor-associated myocarditis and myositis. These new insights should help move the field toward the development of improved treatment and prevention options, ultimately improving patient outcomes. See related article by Siddiqui et al., p. 964 (1).

Comparison of uridine and N1-methylpseudouridine mRNA platforms in development of an Andes virus vaccine

The rodent-borne Andes virus (ANDV) causes a severe disease in humans. We developed an ANDV mRNA vaccine based on the M segment of the viral genome, either with regular uridine (U-mRNA) or N1-methylpseudouridine (m1Ψ-mRNA). Female mice immunized by m1Ψ-mRNA developed slightly greater germinal center (GC) responses than U-mRNA-immunized mice. Single cell RNA and BCR sequencing of the GC B cells revealed similar levels of activation, except an additional cluster of cells exhibiting interferon response in animals vaccinated with U-mRNA but not m1Ψ-mRNA. Similar immunoglobulin class-switching and somatic hypermutations were observed in response to the vaccines. Female Syrian hamsters were immunized via a prime–boost regimen with two doses of each vaccine. The titers of glycoprotein-binding antibodies were greater for U-mRNA construct than for m1Ψ-mRNA construct; however, the titers of ANDV-neutralizing antibodies were similar. Vaccinated animals were challenged with a lethal dose of ANDV, along with a naïve control group. All control animals and two animals vaccinated with a lower dose of m1Ψ-mRNA succumbed to infection whereas other vaccinated animals survived without evidence of virus replication. The data demonstrate the development of a protective vaccine against ANDV and the lack of a substantial effect of m1Ψ modification on immunogenicity and protection in rodents.

Analysis of IDH and EGFR as biomarkers in glioblastoma multiforme: A case-control study

BackgroundGlioblastoma multiforme (GBM) is a very aggressive primary central nervous system (CNS) tumor with limited therapeutic options and poor prognosis. This study aimed to analyze the association between single nucleotide polymorphisms (SNPs), including IDH1 rs121913500C > T, IDH2 rs11540478G > A, and EGFR rs1468727C > T, and their association on the risk and overall survival of GBM patients in Jordan. MethodsUsing a case-control study design involving 63 GBM patients and 226 healthy controls was conducted at King Abdullah University Hospital in Jordan. DNA extraction was performed using formalin-fixed and paraffin-embedded tissue for GBM samples and blood samples for controls. SNPs analysis was performed using the Sequenom iPLEX assay sequencing technique. Survival outcomes were assessed using Cox models and hazard ratios (HR), and single-cell RNA (scRNA) analysis was performed from GSE70630. ResultsThe study showed a significant association between genotype frequency in GBM cases and controls for specific SNPs, including IDH1 rs121913500C > T, and EGFR rs1468727C > T. The G/G genotype of rs11540478 (IDH2) was associated with better prognostic outcomes in GBM patients. The scRNA analysis demonstrated the differential expression of IDH1, IDH2, and EGFR in GBM, with enrichment in central carbon metabolism in cancer. ConclusionOur findings suggest that SNPs, particularly in IDH1 and IDH2 genes and EGFR, may serve as diagnostic and prognostic biomarkers for GBM. While the study underscores the clinical relevance of these genetic variants, further investigations with larger and more diverse populations are essential to validate and extend these associations.

Single-cell total-RNA profiling unveils regulatory hubs of transcription factors

Recent development of RNA velocity uses master equations to establish the kinetics of the life cycle of RNAs from unspliced RNA to spliced RNA (i.e., mature RNA) to degradation. To feed this kinetic analysis, simultaneous measurement of unspliced RNA and spliced RNA in single cells is greatly desired. However, the majority of single-cell RNA-seq chemistry primarily captures mature RNA species to measure gene expressions. Here, we develop a one-step total-RNA chemistry-based single-cell RNA-seq method: snapTotal-seq. We benchmark this method with multiple single-cell RNA-seq assays in their performance in kinetic analysis of cell cycle by RNA velocity. Next, with LASSO regression between transcription factors, we identify the critical regulatory hubs mediating the cell cycle dynamics. We also apply snapTotal-seq to profile the oncogene-induced senescence and identify the key regulatory hubs governing the entry of senescence. Furthermore, from the comparative analysis of unspliced RNA and spliced RNA, we identify a significant portion of genes whose expression changes occur in spliced RNA but not to the same degree in unspliced RNA, indicating these gene expression changes are mainly controlled by post-transcriptional regulation. Overall, we demonstrate that snapTotal-seq can provide enriched information about gene regulation, especially during the transition between cell states.