Analysis Of RNA Sequencing Data Research Articles

Abstract 6839: Comprehensive transcriptomic analysis of acquired resistance after immune-checkpoint inhibitors in non-small cell lung cancer

Abstract Background: Treatment of lung cancer has changed dramatically with the use of immune checkpoint inhibitors (ICI). Nevertheless, most patients underwent drug resistance even though they initially showed a response to ICI, and there is an unmet need to delineate the mechanism. Method: We analyzed RNA sequencing data of paired pre- and post-ICI samples from 31 NSCLC patients. Among them, we identified 17 patients who underwent acquired resistance (AR) who initially showed partial or complete response or stable disease lasting more than six months, and the remaining 14 patients who showed primary resistance (PR) to ICI. The analysis was performed using a tumor microenvironment (TME) classifier, classifying immune-enriched [IE], immune-enriched and fibrotic [IE/F], fibrotic [F] and immune-desert [D] subtype, and CIBERSORTx, the immune cell deconvolutional analysis tool. Results: Most patients with AR were initially classified as IE or IE/F subtypes (N=12, 70.6%), but six out of 14 with PR (42.8%) showed IE or IE/F subtypes. Among 12 patients with initial IE or IE/F subtypes, seven (58.3%) patients turned into the D subtype after AR to ICI. Most patients initially classified as D or F subtypes (N=5, 29.4%) maintained their subtypes (N=4, 80%) after AR.When assessing the pre-ICI samples, memory B cell was enriched, and M2 macrophage and regulatory T cell were depleted in AR patients compared to PR. The memory B cell and resting memory CD4+T cell decreased, and M2 macrophage increased after AR to ICI. These immune cell composition changes were not observed in PR patients. Conclusions: We observed differences in TME phenotype and immune cell composition before ICI between AR & PR patients. The TME change from immune-enriched to desert and cell composition change of increased M2 macrophages and decreased memory B and resting CD4+ T cells after AR to ICI could be significant mechanisms of AR to ICI. Citation Format: Cheol Yong Joe, Yeong Hak Bang, Sehhoon Park, Hyun Ae Jung, Jong-Mu Sun, Yoon-La Choi, Jin Seok Ahn, Myung-Ju Ahn, Se-Hoon Lee. Comprehensive transcriptomic analysis of acquired resistance after immune-checkpoint inhibitors in non-small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6839.

Identification and Validation of T-Cell Exhaustion Signature for Predicting Prognosis and Immune Response in Pancreatic Cancer by Integrated Analysis of Single-Cell and Bulk RNA Sequencing Data.

Pancreatic cancer (PACA) is one of the most fatal malignancies worldwide. Immunotherapy is largely ineffective in patients with PACA. T-cell exhaustion contributes to immunotherapy resistance. We investigated the prognostic potential of T-cell exhaustion-related genes (TEXGs). A single-cell RNA (scRNA) sequencing dataset from Tumor Immune Single-Cell Hub (TISCH) and bulk sequencing datasets from the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) were used to screen differentially expressed TEXGs. Kaplan-Meier survival, LASSO regression, and univariate/multivariate Cox regression analyses were performed to construct a TEXG risk model. This model was used to predict the prognosis, tumor immune microenvironment, and immunotherapy response. The PACA cohorts from the ICGC and GSE71729 datasets were used to validate the risk model. Pan-cancer expression of SPOCK2 was determined using the TISCH database. A six-gene (SPOCK2, MT1X, LIPH, RARRES3, EMP1, and MEG3) risk model was constructed. Patients with low risk had prolonged survival times in both the training (TCGA-PAAD, n = 178) and validation (ICGC-PACA-CA, ICGC-PAAD-US, and GSE71729, n = 412) datasets. Multivariate Cox regression analysis demonstrated that the risk score was an independent prognostic variable for PACA. High-risk patients correlated with their immunosuppressive status. Immunohistochemical staining confirmed the changes in TEXGs in clinical samples. Moreover, pan-cancer scRNA sequencing datasets from TISCH analysis indicated that SPOCK2 may be a novel marker of exhausted CD8+ T-cells. We established and validated a T-cell exhaustion-related prognostic signature for patients with PACA. Moreover, our study suggests that SPOCK2 is a novel marker of exhausted CD8+ T cells.

Revealing the role of SPP1+ macrophages in glioma prognosis and therapeutic targeting by investigating tumor-associated macrophage landscape in grade 2 and 3 gliomas

BackgroundGlioma is a highly heterogeneous brain tumor categorized into World Health Organization (WHO) grades 1–4 based on its malignancy. The suppressive immune microenvironment of glioma contributes significantly to unfavourable patient outcomes. However, the cellular composition and their complex interplays within the glioma environment remain poorly understood, and reliable prognostic markers remain elusive. Therefore, in-depth exploration of the tumor microenvironment (TME) and identification of predictive markers are crucial for improving the clinical management of glioma patients.ResultsOur analysis of single-cell RNA-sequencing data from glioma samples unveiled the immunosuppressive role of tumor-associated macrophages (TAMs), mediated through intricate interactions with tumor cells and lymphocytes. We also discovered the heterogeneity within TAMs, among which a group of suppressive TAMs named TAM-SPP1 demonstrated a significant association with Epidermal Growth Factor Receptor (EGFR) amplification, impaired T cell response and unfavourable patient survival outcomes. Furthermore, by leveraging genomic and transcriptomic data from The Cancer Genome Atlas (TCGA) dataset, two distinct molecular subtypes with a different constitution of TAMs, EGFR status and clinical outcomes were identified. Exploiting the molecular differences between these two subtypes, we developed a four-gene-based prognostic model. This model displayed strong associations with an elevated level of suppressive TAMs and could be used to predict anti-tumor immune response and prognosis in glioma patients.ConclusionOur findings illuminated the molecular and cellular mechanisms that shape the immunosuppressive microenvironment in gliomas, providing novel insights into potential therapeutic targets. Furthermore, the developed prognostic model holds promise for predicting immunotherapy response and assisting in more precise risk stratification for glioma patients.Graphical abstract

Transcriptome sequencing and expression analysis in peanut reveal the potential mechanism response to Ralstonia solanacearum infection

BackgroundBacterial wilt caused by Ralstonia solanacearum severely affects peanut (Arachis hypogaea L.) yields. The breeding of resistant cultivars is an efficient means of controlling plant diseases. Therefore, identification of resistance genes effective against bacterial wilt is a matter of urgency. The lack of a reference genome for a resistant genotype severely hinders the process of identification of resistance genes in peanut. In addition, limited information is available on disease resistance-related pathways in peanut.ResultsFull-length transcriptome data were used to generate wilt-resistant and -susceptible transcript pools. In total, 253,869 transcripts were retained to form a reference transcriptome for RNA-sequencing data analysis. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differentially expressed genes revealed the plant-pathogen interaction pathway to be the main resistance-related pathway for peanut to prevent bacterial invasion and calcium plays an important role in this pathway. Glutathione metabolism was enriched in wilt-susceptible genotypes, which would promote glutathione synthesis in the early stages of pathogen invasion. Based on our previous quantitative trait locus (QTL) mapping results, the genes arahy.V6I7WA and arahy.MXY2PU, which encode nucleotide-binding site-leucine-rich repeat receptor proteins, were indicated to be associated with resistance to bacterial wilt.ConclusionsThis study identified several pathways associated with resistance to bacterial wilt and identified candidate genes for bacterial wilt resistance in a major QTL region. These findings lay a foundation for investigation of the mechanism of resistance to bacterial wilt in peanut.

SOCS3 regulates pathological retinal angiogenesis through modulating SPP1 expression in microglia and macrophages

Pathological ocular angiogenesis has long been associated with myeloid cell activation. However, the precise cellular and molecular mechanisms governing the intricate crosstalk between the immune system and vascular changes during ocular neovascularization formation remain elusive. In this study, we demonstrated that the absence of the suppressor of cytokine signaling 3 (SOCS3) in myeloid cells led to a substantial accumulation of microglia and macrophage subsets during the neovascularization process. Our single-cell RNA sequencing data analysis revealed a remarkable increase in the expression of the secreted phosphoprotein 1 (Spp1) gene within these microglia and macrophages, identifying subsets of Spp1-expressing microglia and macrophages during neovascularization formation in angiogenesis mouse models. Notably, the number of Spp1-expressing microglia and macrophages exhibited further elevation during neovascularization in mice lacking myeloid SOCS3. Moreover, our investigation unveiled the Spp1 gene as a direct transcriptional target gene of signal transducer and activator of transcription 3. Importantly, pharmaceutical activation of SOCS3 or blocking of SPP1 resulted in a significant reduction in pathological neovascularization. In conclusion, our study highlights the pivotal role of the SOCS3/STAT3/SPP1 axis in the regulation of pathological retinal angiogenesis.

Sildenafil as a Candidate Drug for Alzheimer's Disease: Real-World Patient Data Observation and Mechanistic Observations from Patient-Induced Pluripotent Stem Cell-Derived Neurons.

Alzheimer's disease (AD) is a chronic neurodegenerative disease needing effective therapeutics urgently. Sildenafil, one of the approved phosphodiesterase-5 inhibitors, has been implicated as having potential effect in AD. To investigate the potential therapeutic benefit of sildenafil on AD. We performed real-world patient data analysis using the MarketScan® Medicare Supplemental and the Clinformatics® databases. We conducted propensity score-stratified analyses after adjusting confounding factors (i.e., sex, age, race, and comorbidities). We used both familial and sporadic AD patient induced pluripotent stem cells (iPSC) derived neurons to evaluate the sildenafil's mechanism-of-action. We showed that sildenafil usage is associated with reduced likelihood of AD across four new drug compactor cohorts, including bumetanide, furosemide, spironolactone, and nifedipine. For instance, sildenafil usage is associated with a 54% reduced incidence of AD in MarketScan® (hazard ratio [HR] = 0.46, 95% CI 0.32- 0.66) and a 30% reduced prevalence of AD in Clinformatics® (HR = 0.70, 95% CI 0.49- 1.00) compared to spironolactone. We found that sildenafil treatment reduced tau hyperphosphorylation (pTau181 and pTau205) in a dose-dependent manner in both familial and sporadic AD patient iPSC-derived neurons. RNA-sequencing data analysis of sildenafil-treated AD patient iPSC-derived neurons reveals that sildenafil specifically target AD related genes and pathobiological pathways, mechanistically supporting the beneficial effect of sildenafil in AD. These real-world patient data validation and mechanistic observations from patient iPSC-derived neurons further suggested that sildenafil is a potential repurposable drug for AD. Yet, randomized clinical trials are warranted to validate the causal treatment effects of sildenafil in AD.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

Key messageBarley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5′ terminal nucleotide of small RNAs, demonstrating functional conservation and divergence.The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5′ adenine residue, while also accepting 5′ guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5′ adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

ACKR3 in olfactory glia cells shapes the immune defense of the olfactory mucosa.

Barrier-forming olfactory glia cells, termed sustentacular cells, play important roles for immune defense of the olfactory mucosa, for example as entry sites for SARS-CoV-2 and subsequent development of inflammation-induced smell loss. Here we demonstrate that sustentacular cells express ACKR3, a chemokine receptor that functions both as a scavenger of the chemokine CXCL12 and as an activator of alternative signaling pathways. Differential gene expression analysis of bulk RNA sequencing data obtained from WT and ACKR3 conditional knockout mice revealed upregulation of genes involved in immune defense. To map the regulated genes to the different cell types of the olfactory mucosa, we employed biocomputational methods utilizing a single-cell reference atlas. Transcriptome analysis, PCR and immunofluorescence identified up-regulation of NF-κB-related genes, known to amplify inflammatory signaling and to facilitate leukocyte transmigration, in the gliogenic lineage. Accordingly, we found a marked increase in leukocyte-expressed genes and confirmed leukocyte infiltration into the olfactory mucosa. In addition, lack of ACKR3 led to enhanced expression and secretion of early mediators of immune defense by Bowman's glands. As a result, the number of apoptotic cells in the epithelium was decreased. In conclusion, our research underlines the importance of sustentacular cells in immune defense of the olfactory mucosa. Moreover, it identifies ACKR3, a druggable G protein-coupled receptor, as a promising target for modulation of inflammation-associated anosmia.

Single-cell integrative analysis reveals consensus cancer cell states and clinical relevance in breast cancer

High heterogeneity and complex interactions of malignant cells in breast cancer has been recognized as a driver of cancer progression and therapeutic failure. However, complete understanding of common cancer cell states and their underlying driver factors remain scarce and challenging. Here, we revealed seven consensus cancer cell states recurring cross patients by integrative analysis of single-cell RNA sequencing data of breast cancer. The distinct biological functions, the subtype-specific distribution, the potential cells of origin and the interrelation of consensus cancer cell states were systematically elucidated and validated in multiple independent datasets. We further uncovered the internal regulons and external cell components in tumor microenvironments, which contribute to the consensus cancer cell states. Using the state-specific signature, we also inferred the abundance of cells with each consensus cancer cell state by deconvolution of large breast cancer RNA-seq cohorts, revealing the association of immune-related state with better survival. Our study provides new insights for the cancer cell state composition and potential therapeutic strategies of breast cancer.

First Detection of Benign Rabbit Caliciviruses in Chile.

Pathogenic lagoviruses (Rabbit hemorrhagic disease virus, RHDV) are widely spread across the world and are used in Australia and New Zealand to control populations of feral European rabbits. The spread of the non-pathogenic lagoviruses, e.g., rabbit calicivirus (RCV), is less well studied as the infection results in no clinical signs. Nonetheless, RCV has important implications for the spread of RHDV and rabbit biocontrol as it can provide varying levels of cross-protection against fatal infection with pathogenic lagoviruses. In Chile, where European rabbits are also an introduced species, myxoma virus was used for localised biocontrol of rabbits in the 1950s. To date, there have been no studies investigating the presence of lagoviruses in the Chilean feral rabbit population. In this study, liver and duodenum rabbit samples from central Chile were tested for the presence of lagoviruses and positive samples were subject to whole RNA sequencing and subsequent data analysis. Phylogenetic analysis revealed a novel RCV variant in duodenal samples that likely originated from European RCVs. Sequencing analysis also detected the presence of a rabbit astrovirus in one of the lagovirus-positive samples.

Potential therapeutic targets for hypotension in duchenne muscular dystrophy

Duchenne Muscular Dystrophy (DMD) is marked by genetic mutations occurring in the DMD gene, which is widely expressed in the cardiovascular system. In addition to developing cardiomyopathy, patients with DMD have been reported to be susceptible to the development of symptomatic hypotension, although the mechanisms are unclear. Analysis of single-cell RNA sequencing data has identified potassium voltage-gated channel subfamily Q member 5 (KCNQ5) and possibly ryanodine receptor 2 (RyR2) as potential candidate hypotension genes whose expression is significantly upregulated in the vascular smooth muscle cells of DMD mutant mice. We hypothesize that heightened KCNQ5 and RyR2 expression contributes to decreased arterial blood pressure in patients with DMD. Exploring pharmacological approaches to inhibit the KCNQ5 and RyR2 channels holds promise in managing the systemic hypotension observed in individuals with DMD. This avenue of investigation presents new prospects for improving clinical outcomes for these patients.

ICELLNET v2: a versatile method for cell-cell communication analysis from human transcriptomic data.

Several methods have been developed in the past years to infer cell-cell communication networks from transcriptomic data based on ligand and receptor expression. Among them, ICELLNET is one of the few approaches to consider the multiple subunits of ligands and receptors complexes to infer and quantify cell communication. In here, we present a major update of ICELLNET. As compared to its original implementation, we (i) drastically expanded the ICELLNET ligand-receptor database from 380 to 1669 biologically curated interactions, (ii) integrated important families of communication molecules involved in immune crosstalk, cell adhesion, and Wnt pathway, (iii) optimized ICELLNET framework for single-cell RNA sequencing data analyses, (iv) provided new visualizations of cell-cell communication results to facilitate prioritization and biological interpretation. This update will broaden the use of ICELLNET by the scientific community in different biological fields. ICELLNET package is implemented in R. Source code, documentation and tutorials are available on GitHub (https://github.com/soumelis-lab/ICELLNET).

Targeted deletion of CD244 on monocytes promotes differentiation into anti-tumorigenic macrophages and potentiates PD-L1 blockade in melanoma

BackgroundIn the myeloid compartment of the tumor microenvironment, CD244 signaling has been implicated in immunosuppressive phenotype of monocytes. However, the precise molecular mechanism and contribution of CD244 to tumor immunity in monocytes/macrophages remains elusive due to the co-existing lymphoid cells expressing CD244.MethodsTo directly assess the role of CD244 in tumor-associated macrophages, monocyte-lineage-specific CD244-deficient mice were generated using cre-lox recombination and challenged with B16F10 melanoma. The phenotype and function of tumor-infiltrating macrophages along with antigen-specific CD8 T cells were analyzed by flow cytometry and single cell RNA sequencing data analysis, and the molecular mechanism underlying anti-tumorigenic macrophage differentiation, antigen presentation, phagocytosis was investigated ex vivo. Finally, the clinical feasibility of CD244-negative monocytes as a therapeutic modality in melanoma was confirmed by adoptive transfer experiments.ResultsCD244fl/flLysMcre mice demonstrated a significant reduction in tumor volume (61% relative to that of the CD244fl/fl control group) 14 days after tumor implantation. Within tumor mass, CD244fl/flLysMcre mice also showed higher percentages of Ly6Clow macrophages, along with elevated gp100+IFN-γ+ CD8 T cells. Flow cytometry and RNA sequencing data demonstrated that ER stress resulted in increased CD244 expression on monocytes. This, in turn, impeded the generation of anti-tumorigenic Ly6Clow macrophages, phagocytosis and MHC-I antigen presentation by suppressing autophagy pathways. Combining anti-PD-L1 antibody with CD244−/− bone marrow-derived macrophages markedly improved tumor rejection compared to the anti-PD-L1 antibody alone or in combination with wild-type macrophages. Consistent with the murine data, transcriptome analysis of human melanoma tissue single-cell RNA-sequencing dataset revealed close association between CD244 and the inhibition of macrophage maturation and function. Furthermore, the presence of CD244-negative monocytes/macrophages significantly increased patient survival in primary and metastatic tumors.ConclusionOur study highlights the novel role of CD244 on monocytes/macrophages in restraining anti-tumorigenic macrophage generation and tumor antigen-specific T cell response in melanoma. Importantly, our findings suggest that CD244-deficient macrophages could potentially be used as a therapeutic agent in combination with immune checkpoint inhibitors. Furthermore, CD244 expression in monocyte-lineage cells serve as a prognostic marker in cancer patients.

Unraveling the connection between Hashimoto's Thyroiditis and non-alcoholic fatty liver disease: exploring the role of CD4+central memory T cells through integrated genetic approaches.

Exploring the connection between Hashimoto's thyroiditis (HT) and non-alcoholic fatty liver disease (NAFLD) through integrated genetic approaches. We utilized integrated genetic approaches, such as single-cell RNA sequencing (scRNA-seq) data analysis, Mendelian Randomization (MR), colocalization analysis, cell communication, and metabolic analyses, to investigate potential correlations between HT and NAFLD. Through the integrated analysis of scRNA-seq data from individuals with HT, NAFLD, and healthy controls, we observed an upregulation in the proportion of CD4+central memory (CD4+CM) T cells among T cells in both diseases. A total of 63 differentially expressed genes (DEGs) were identified in the CD4+CM cells after the differential analysis. By using MR, 8 DEGs (MAGI3, CSGALNACT1, CAMK4, GRIP1, TRAT1, IL7R, ERN1, and MB21D2) were identified to have a causal relationship with HT, and 4 DEGs (MAGI3, RCAN3, DOCK10, and SAMD12) had a causal relationship with NAFLD. MAGI3 was found to be causally linked to both HT and NAFLD. Therefore, MAGI3 was designated as the marker gene. Reverse MR and Steiger filtering showed no evidence of reverse causality. Colocalization analyses further indicated close links between MAGI3 and HT as well as NAFLD. Finally, based on the expression levels of MAGI3, we stratified CD4+CM cells into two subsets: MAGI3+CD4+CM cells and MAGI3-CD4+CM cells. Functional analyses revealed significant differences between the two subsets, potentially related to the progression of the two diseases. This study delves into the potential connections between HT and NAFLD through integrated genetic methods. Our research reveals an elevated proportion of CD4+CM cells within T cells in both HT and NAFLD. Through MR and colocalization analysis, we identify specific genes causally linked to HT and NAFLD, such as MAGI3. Ultimately, based on MAGI3 expression levels, we categorize CD4+CM cells into MAGI3+CD4+CM cells and MAGI3-CD4+CM cells, uncovering significant differences between them through functional analyses.

Comprehensive genomic and clinical analyses identify APOBEC mutational signatures as a brain metastasis risk factor in lung adenocarcinoma patients

BackgroundLung adenocarcinoma is the most common source of brain metastasis (BM), resulting in significant morbidity and mortality. We aimed to identify patients with high BM risk who possibly benefit from brain-penetrant drugs, prophylactic cranial irradiation, or close brain magnetic resonance imaging surveillance. MethodsMetastatic lung adenocarcinoma patients with extracranial tumor samples profiled by a next-generation sequencing panel targeting 425 tumor-related genes were retrospectively enrolled between February 2008 and July 2021. We compared BM and non-BM patients’ genomic and clinical features and studied their associations with BM risk. Two external cohorts were used for result validation and molecular mechanisms investigation, respectively. ResultsWe included 174 eligible patients, including 90 having developed BM by the end of follow-up. Age≤60, EGFR activating mutations, and high-level apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) mutational signatures were associated with elevated BM risk. Similar findings in BM-free survival were obtained by fitting Fine-Gray subdistribution hazard models addressing competing risks. Increased BM risk related to APOBEC mutational signatures was validated in an external cohort (N = 440). RNA sequencing data analyses performed in another external cohort (N = 230) revealed that expressions of metastasis-related pathways such as transforming growth factor (TGF)β and epithelial-mesenchymal transition (EMT) were upregulated in the patients with high-level APOBEC mutational signatures. ConclusionAPOBEC mutational signatures related to upregulated TGFβ and EMT, could serve as an independent risk factor for BM and BM-free survival in metastatic lung adenocarcinoma patients. Further investigations are warranted to tailor personalized treatments to improve the susceptible patient's outcomes.

Curcumin suppresses metastasis of triple-negative breast cancer cells by modulating EMT signaling pathways: An integrated study of bioinformatics analysis.

This study aimed to use bioinformatics approaches for predicting the anticancer mechanisms of curcumin on triple-negative breast cancer (TNBC) and to verify these predictions through in vitro experiments. Initially, the Cell Counting Kit-8 (CCK8) assay was employed to rigorously investigate the influence of curcumin on the proliferative capacity of TNBC cells. Subsequently, flow cytometry was employed to meticulously assess the impact of curcumin on cellular apoptosis and the cell cycle regulation. Transwell assays were employed to meticulously evaluate the effect of curcumin on the motility of TNBC cells. RNA sequencing was conducted, followed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of differentially expressed genes, aiming to elucidate the potential anticancer mechanisms underlying curcumin's effects. To thoroughly elucidate the interactions among multiple proteins, we constructed a protein-protein interaction (PPI) network. Finally, the expression levels of several key proteins, including fibronectin, mTOR, β-Catenin, p-Akt, Akt, N-Cadherin, p-S6, and S6, were assessed using the western blot. The CCK8 assay results showed that curcumin significantly inhibited the proliferation of Hs578T and MDA-MB-231 cells. Flow cytometry results showed that curcumin induced apoptosis in these cells and arrested the cell cycle at the G2/M phase. Additionally, Transwell assay results showed that curcumin effectively reduced the motility of Hs578T and MDA-MB-231 cells. Enrichment analysis of RNA sequencing data showed that the mechanism of action of curcumin was significantly associated with signaling pathways such as pathways in cancer, focal adhesion, and PI3K-Akt signaling pathways. Subsequently, we constructed a protein-protein interaction network to elucidate the interactions among multiple proteins. Finally, Western blotting analysis showed that curcumin significantly decreased the expression levels of key proteins including Fibronectin, mTOR, β-Catenin, p-Akt, Akt, N-Cadherin, p-S6, and S6. Curcumin exhibits its therapeutic potential in TNBC by modulating multiple signaling pathways. It may inhibit the epithelial-mesenchymal transition process by downregulating the expression of proteins involved in the mTOR and PI3K-Akt signaling pathways, thereby suppressing the motility of TNBC cells. These findings provide experimental evidence for considering curcumin as a potential therapeutic strategy in the treatment of TNBC.

Single cell-transcriptomic analysis informs the lncRNA landscape in metastatic castration resistant prostate cancer

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal form of prostate cancer. Although long-noncoding RNAs (lncRNAs) have been implicated in mCRPC, past studies have relied on bulk sequencing methods with low depth and lack of single-cell resolution. Hence, we performed a lncRNA-focused analysis of single-cell RNA-sequencing data (n = 14) from mCRPC biopsies followed by integration with bulk multi-omic datasets. This yielded 389 cell-enriched lncRNAs in prostate cancer cells and the tumor microenvironment (TME). These lncRNAs demonstrated enrichment with regulatory elements and exhibited alterations during prostate cancer progression. Prostate-lncRNAs were correlated with AR mutational status and response to treatment with enzalutamide, while TME-lncRNAs were associated with RB1 deletions and poor prognosis. Finally, lncRNAs identified between prostate adenocarcinomas and neuroendocrine tumors exhibited distinct expression and methylation profiles. Our findings demonstrate the ability of single-cell analysis to refine our understanding of lncRNAs in mCRPC and serve as a resource for future mechanistic studies.

A multi-faceted approach to unravel coding and non-coding gene fusions and target chimeric proteins in ataxia

Ataxia represents a heterogeneous group of neurodegenerative disorders characterized by a loss of balance and coordination, often resulting from mutations in genes vital for cerebellar function and maintenance. Recent advances in genomics have identified gene fusion events as critical contributors to various cancers and neurodegenerative diseases. However, their role in ataxia pathogenesis remains largely unexplored. Our study Hdelved into this possibility by analyzing RNA sequencing data from 1443 diverse samples, including cell and mouse models, patient samples, and healthy controls. We identified 7067 novel gene fusions, potentially pivotal in disease onset. These fusions, notably in-frame, could produce chimeric proteins, disrupt gene regulation, or introduce new functions. We observed conservation of specific amino acids at fusion breakpoints and identified potential aggregate formations in fusion proteins, known to contribute to ataxia. Through AI-based protein structure prediction, we identified topological changes in three high-confidence fusion proteins—TEN1-ACOX1, PEX14-NMNAT1, and ITPR1-GRID2—which could potentially alter their functions. Subsequent virtual drug screening identified several molecules and peptides with high-affinity binding to fusion sites. Molecular dynamics simulations confirmed the stability of these protein-ligand complexes at fusion breakpoints. Additionally, we explored the role of non-coding RNA fusions as miRNA sponges. One such fusion, RP11-547P4-FLJ33910, showed strong interaction with hsa-miR-504-5p, potentially acting as its sponge. This interaction correlated with the upregulation of hsa-miR-504-5p target genes, some previously linked to ataxia. In conclusion, our study unveils new aspects of gene fusions in ataxia, suggesting their significant role in pathogenesis and opening avenues for targeted therapeutic interventions. Communicated by Ramaswamy H. Sarma

Integrated analysis of single-cell and bulk RNA sequencing data reveals the association between hypoxic tumor cells and exhausted T cells in predicting immune therapy response

Continuous stimulation of tumor neoantigens and various cytokines in the tumor microenvironment leads to T cell dysfunction, but the specific mechanisms by which these key factors are distributed among different cell subpopulations and how they affect patient outcomes and treatment response are incompletely characterized.By integrating single-cell and bulk sequencing data of non-small cell lung cancer patients, we constructed a clinical outcome-associated T cell exhaustion signature. We discovered a significant association between the T cell exhaustion state and tumor cell hypoxia. Hypoxic malignant cells were significantly correlated with the proportion of exhausted T cells, and they co-occurred in patients at advanced stage. By analyzing the ligand-receptor interactions between these two cell states, we observed that T cells were recruited towards tumor cells through production of chemokines such as CXCL16-CXCR6 axis and CCL3/CCL4/CCL5-CCR5 axis. Based on 15 immune checkpoint blockade (ICB)-treatment cohorts, we constructed an interaction signature that can be used to predict the response to immune checkpoint blockade therapy. Among genes composed of the signature, CXCR6 alone has similarly high prediction efficacy (Area Under Curve (AUC) = 1, 0.89 and 0.73 for GSE126044, GSE135222 and GSE93157, respectively) with the signature and thus could serve as a potential biomarker for predicting immunotherapy response.Together, we have discovered and validated a significant association between exhausted T cells and hypoxic malignant cells, elucidating key interaction factors that significantly associated with response to immunotherapy.

High CD133 expression in proximal tubular cells in diabetic kidney disease: good or bad?

BackgroundProximal tubular cells (PTCs) play a critical role in the progression of diabetic kidney disease (DKD). As one of important progenitor markers, CD133 was reported to indicate the regeneration of dedifferentiated PTCs in acute kidney disease. However, its role in chronic DKD is unclear. Therefore, we aimed to investigate the expression patterns and elucidate its functional significance of CD133 in DKD.MethodsData mining was employed to illustrate the expression and molecular function of CD133 in PTCs in human DKD. Subsequently, rat models representing various stages of DKD progression were established. The expression of CD133 was confirmed in DKD rats, as well as in human PTCs (HK-2 cells) and rat PTCs (NRK-52E cells) exposed to high glucose. The immunofluorescence and flow cytometry techniques were utilized to determine the expression patterns of CD133, utilizing proliferative and injury indicators. After overexpression or knockdown of CD133 in HK-2 cells, the cell proliferation and apoptosis were detected by EdU assay, real-time cell analysis and flow analysis. Additionally, the evaluation of epithelial, progenitor cell, and apoptotic indices was performed through western blot and quantitative RT-PCR analyses.ResultsThe expression of CD133 was notably elevated in both human and rat PTCs in DKD, and this expression increased as DKD progressed. CD133 was found to be co-expressed with CD24, KIM-1, SOX9, and PCNA, suggesting that CD133+ cells were damaged and associated with proliferation. In terms of functionality, the knockdown of CD133 resulted in a significant reduction in proliferation and an increase in apoptosis in HK-2 cells compared to the high glucose stimulus group. Conversely, the overexpression of CD133 significantly mitigated high glucose-induced cell apoptosis, but had no impact on cellular proliferation. Furthermore, the Nephroseq database provided additional evidence to support the correlation between CD133 expression and the progression of DKD. Analysis of single-cell RNA-sequencing data revealed that CD133+ PTCs potentially play a role in the advancement of DKD through multiple mechanisms, including heat damage, cell microtubule stabilization, cell growth inhibition and tumor necrosis factor-mediated signaling pathway.ConclusionOur study demonstrates that the upregulation of CD133 is linked to cellular proliferation and protects PTC from apoptosis in DKD and high glucose induced PTC injury. We propose that heightened CD133 expression may facilitate cellular self-protective responses during the initial stages of high glucose exposure. However, its sustained increase is associated with the pathological progression of DKD. In conclusion, CD133 exhibits dual roles in the advancement of DKD, necessitating further investigation.