Search

Background: The use of neoadjuvant chemoradiotherapy (neoCRT) followed by surgery has markedly enhanced the quality of survival in patients suffering from local advanced rectal cancer (LARC). Enhancing this treatment requires a deep understanding of its underlying mechanism. The heterogeneous nature of the tumor microenvironment (TME) significantly impacts therapeutic responses, presenting complex therapeutic challenges., Methods: In this comprehensive study, we explored the intricate cellular and molecular shifts within the TME of LARC after neoCRT administration. Using single-cell transcriptomic analysis, we meticulously examined 32,417 cells sourced from six samples, each representing different tumor regression grades (TRG: 0 versus 2). This detailed analysis enabled us to characterize the various cell subpopulations, encompassing epithelial cells, lymphocytes, myeloid cells, endothelial cells, and fibroblasts. Additionally, we identified their marker genes for deconvolution calculation in the READ cohort of the TCGA project. And we obtain their marker genes for deconvolution calculation in the READ cohort of the TCGA project., Results: Through cluster analysis and pathway comparisons of malignant tumor cells, we discerned that samples with poor tumor regression exhibit enhanced metabolic versatility and adaptability, enabling them to counteract the impacts of both radiotherapy and chemotherapy. Interestingly, within the TRG2 cohort, we observed a predominant immunosuppressive state in the TME, characterized by the activation of CD4 + regulatory T cells, maintained CD8 + T cell functionality, and a heightened M1 to M2 macrophage ratio. Moreover, the differing outcomes of neoCRT were reflected in the varying interaction dynamics between macrophages (M1 and M2) and CD4+/CD8 + T cells. Furthermore, our data reveal that neoCRT intricately modulates fibroblasts and endothelial cells, primarily through the extracellular matrix remodeling pathway, which orchestrates tumor angiogenesis. All changes were validated through immunofluorescence staining on intraoperative samples before and after treatment. To summarize, our investigation presents a comprehensive exploration of the cellular and molecular metamorphoses within the TME post-neoCRT., Conclusions: By unveiling the sophisticated interaction between the multifaceted cells within the TME and their respective reactions to neoCRT, we establish a robust platform for ensuing future investigations. This study paves the way for novel therapeutic strategies that leverage these insights to bolster the efficacy of neoCRT in managing LARC., Competing Interests: Declarations Ethics approval and consent to participate The study was conducted in accordance with the Declaration of Helsinki, as revised in 2013. Conflict of interest The authors declare there are no potential conflicts of interest., (© 2024. The Author(s).)

Gastroesophageal disorders and cancers impose a significant global burden. Particularly, the prevalence of esophageal adenocarcinoma (EAC) has increased dramatically in recent years. Barrett's esophagus, a precursor of EAC, features a unique tissue adaptation at the gastroesophageal squamo-columnar junction (GE-SCJ), where the esophagus meets the stomach. Investigating the evolution of GE-SCJ and understanding dysregulation in its homeostasis are crucial for elucidating cancer pathogenesis. Here, we present the technical quality of the comprehensive single-cell RNA sequencing (scRNA-seq) dataset from mice that captures the transcriptional dynamics during the development of the esophagus, stomach and the GE-SCJ at embryonic, neonatal and adult stages. Through integration with external scRNA-seq datasets and validations using organoid and animal models, we demonstrate the dataset's consistency in identified cell types and transcriptional profiles. This dataset will be a valuable resource for studying developmental patterns and associated signaling networks in the tissue microenvironment. By offering insights into cellular programs during homeostasis, it facilitates the identification of changes leading to conditions like metaplasia and cancer, crucial for developing effective intervention strategies., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Previous studies have demonstrated the clinical relevance of aberrant serum immunoglobulin G (IgG) N-glycomic profiles in liver fibrosis and the pathogenic effects of agalactosyl IgG on activating hepatic stellate cells (HSCs). However, the dynamics of gene expression changes during HSC activation by agalactosyl IgG remain poorly understood. We performed RNA sequencing to analyze the mRNAome of human LX-2 HSCs at multiple time points after treatment with agalactosyl IgG and then compared these results with those obtained after normal IgG and transforming growth factor (TGF)-β1 treatments. Gene expression changes were significantly pronounced on day 5 and subsided by day 11 after HSC activation. A high degree of similarity in gene expression patterns between HSCs treated with agalactosyl IgG and TGF-β1 was observed, of which 1796 and 1785 differentially expressed genes (DEGs) were identified, respectively. Disease ontology analyses revealed that 114 and 105 DEGs in activated HSCs following agalactosyl IgG and TGF-β1 treatments, respectively, were linked to liver cirrhosis, hepatitis, fatty liver disease, hepatitis B, and alcoholic hepatitis, with CCL5 and FAS being the most commonly affected genes. DEGs associated with liver fibrosis or aforementioned liver diseases involved in gene annotation, physiological functions, and signaling pathways regarding secretion of cytokines and chemokines, expression of fibrosis-related growth factors and their receptors, modification of extracellular matrices, and regulation of cell viability in activated HSCs. In conclusion, this study characterized the dynamics of mRNAome and gene networks and identified the liver fibrosis-related DEGs during HSC activation by agalactosyl IgG and TGF-β1., Competing Interests: Declarations Informed consent Not applicable Competing interests The authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Melon thrips, Thrips palmi, represent a significant threat to plants, inducing necrosis and acting as vectors for numerous plant viruses. Entomopathogenic fungi present a promising avenue for the management of melon thrips populations resistant to conventional chemical treatments. In this work, an adult colony of melon thrips was exposed to Beauveria bassiana strain JEF-350, and the ensuing transcriptional response of the infected thrips was scrutinized to elucidate their reactions during fungal pathogenesis. Utilizing Illumina sequencing, RNA samples were extracted from untreated thrips as well as from thrips continuously infected for 2 and 4 days, each with three biological replicates. While no notable alterations in gene expression were observed between the untreated control and thrips infected for 2 days, those infected for 4 days exhibited a plethora of differentially expressed genes. Specifically, in the thrips infected for the extended period, pathways associated with lysosomal function and insect hormone biosynthesis were notably repressed, while others such as serine and glycine metabolism, Toll and Imd, and circadian rhythm pathways displayed heightened activity. Noteworthy downregulation was observed in numerous lysosomal hydrolase genes encoding glycosidases, sulfatases, and lipases, particularly glycosidases. Furthermore, certain genes related to hydrolase precursors within the Golgi apparatus exhibited heightened expression levels but failed to progress toward hydrolase biosynthesis. Upstream regulation of juvenile hormone biosynthesis was augmented, yet downstream genes were significantly downregulated, leading to a disruption in juvenile hormone production. Similarly, while cytochrome P450 genes in the downstream of ecdysone biosynthesis were upregulated, expressions of cholesterol desaturase and cytochrome P450 genes in the upstream were inhibited, consequently dampening ecdysone biosynthesis. The observed differential targeting of organs or pathways by B. bassiana JEF-350, in contrast to conventional chemicals primarily affecting neurotransmission and energy production, suggests its potential efficacy in managing resistant thrips populations. Consequently, integrating JEF-350 into the chemical spray regimen or incorporating it into tank-mix formulations with chemical insecticides emerges as a pragmatic approach within the realm of integrated pest management strategies. KEY POINTS: • Beauveria treatment inhibited lysosomal function and hormone synthesis in thrips. • Thrips serine/glycine metabolism, Toll and Imd, and circadian rhythm pathways were activated. • Upstream functions of thrips hormone biosynthesis increased, while downstream functions were suppressed. • Regarding biosynthesis of metabolites, this fungus targets other pathways with resistance management., Competing Interests: Declarations Conflict of interest The authors declare no competing interests., (© 2024. The Author(s).)

Background: Global warming-induced environmental stresses have diverse effects on gene expression and regulation in the life processes of various aquatic organisms. N6 adenylate methylation (m6A) modifications are known to influence mRNA transcription, localization, translation, stability, splicing, and nuclear export, which are pivotal in mediating stress responses. Apostichopus japonicus is a significant species in aquaculture and a representative of benthic organisms in ecosystems, thus there is a growing need for research on its heat stress mechanism., Results: In this study, m6A-modified whole transcriptome profiles of the respiratory tree tissues of A. japonicus in the control (T18) and high-temperature stress (T32) groups were obtained using MeRIP-seq technology. The results showed that 7,211 common m6A peaks, and 9,459 genes containing common m6A were identified in three replicates T18 and T32 groups. The m6A peaks were found to be highly enriched in the 3' untranslated region, and the common sequence of the m6A peak was also enriched, which was shown as RRACH (R = G or A; H = A, C, or U). A total of 1,200 peaks were identified as significantly differentially enriched in the T32 group compared with the T18 group. Among them, 245 peaks were upregulated and 955 were downregulated, which indicated that high temperature stress significantly altered the methylation pattern of m6A, and there were more demethylation sites in the T32 group. Conjoint analysis of the m6A methylation modification and the transcript expression level (the MeRIP-seq and RNA-seq data) showed co-differentiated 395 genes were identified, which were subsequently divided into four groups with a predominant pattern that more genes with decreased m6A modification and up-regulated expression, including HSP70IV, EIF2AK1, etc. GO enrichment and KEGG analyses of differential m6A peak related genes and co-differentiated genes showed the genes were significantly associated with transcription process and pathways such as protein processing in the endoplasmic reticulum, Wnt signaling pathway, and mTOR signaling pathway, etc. CONCLUSION: The comparisons of m6A modification patterns and conjoint analyses of m6A modification and gene expression profiles suggest that m6A modification was involved in the regulation of heat stress-responsive genes and important functional pathways in A. japonicus in response to high-temperature stress. The study will contribute to elucidate the regulatory mechanism of m6A modification in the response of A. japonicus to environmental stress, as well as the conservation and utilization of sea cucumber resources in the context of environmental changes., Competing Interests: Declarations Ethics approval and consent to participate The animal study was approved by the Institutional Animal Care and Use Committee of the Ludong University (protocol number LDU-IRB20210308NXY). The study was conducted in accordance with the local legislation and institutional requirements. Consent for publication Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Nervous system cancers exhibit diverse transcriptional cell states influenced by normal development, injury response, and growth. However, the understanding of these states' regulation and pharmacological relevance remains limited. Here we present "single-cell regulatory-driven clustering" (scregclust), a method that reconstructs cellular regulatory programs from extensive collections of single-cell RNA sequencing (scRNA-seq) data from both tumors and developing tissues. The algorithm efficiently divides target genes into modules, predicting key transcription factors and kinases with minimal computational time. Applying this method to adult and childhood brain cancers, we identify critical regulators and suggest interventions that could improve temozolomide treatment in glioblastoma. Additionally, our integrative analysis reveals a meta-module regulated by SPI1 and IRF8 linked to an immune-mediated mesenchymal-like state. Finally, scregclust's flexibility is demonstrated across 15 tumor types, uncovering both pan-cancer and specific regulators. The algorithm is provided as an easy-to-use R package that facilitates the exploration of regulatory programs underlying cell plasticity., (© 2024. The Author(s).)

The nucleus pulposus (NP) plays a vital role in intervertebral disc degeneration (IVDD). Previous studies have revealed cellular heterogeneity in NP tissue during IVDD progression. However, the cellular and molecular alterations of diverse cell clusters during IVDD remain to be fully elucidated. NP tissues were isolated from patients with different grades of IVDD undergoing discectomy, and then subjected to single-cell RNA sequencing (scRNA-seq). Cell subsets were identified based on unbiased clustering of gene expression profiles. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to determine the molecular features of diverse cell clusters. Monocle analysis was used to illustrate the differentiation trajectories of chondrocytes. Additionally, CellPhoneDB analysis revealed potential interactions between chondrocytes and other cells during IVDD. Based on the expression profiles of 47,610 individual cells, eight putative clusters including chondrocytes, endothelial cells, fibroblasts, macrophages, mural cells, osteoclasts, proliferating stromal cells and T cells were identified. The chondrocyte cluster was classified into three subsets, C1-C3, which were associated with stress-resistance, fibrosis and inflammatory responses, respectively. Pseudo-time trajectories suggested that chondrocytes gradually differentiated into fibroblasts during IVDD. Immune cells including cDC2s, macrophages and monocytes were identified. Further analysis showed that chondrocytes might communicate with immune cells via the MIF, TNFSF9, SPP1 and CCL4L2 signaling pathways. In addition, we found that invading endothelial cells might interact with chondrocytes through the COL4A1, CXCL12, VEGFA and SEMA3E signaling pathways. Our results reveal the cellular complexity and phenotypic characteristics of NP tissues at single-cell resolution, which will contribute to the in-depth investigation of preventative and regenerative strategies for IVDD., (© 2024. The Author(s).)

Key Message: The QTLs and candidate genes governing the multifoliolate phenotype were identified by combining linkage mapping with BSR-seq, revealing a possible interplay between genetics and the environment in soybean leaf development. Soybean, as a legume, is typified by trifoliolate leaves. Although multifoliolate leaves (compound leaves with more than three leaflets each) have been reported in soybean, including sporadic appearances in the first compound leaves in a recombinant inbred line (RIL) population from a cross between cultivated soybean C08 and wild soybean W05 from this study, the genetic basis of this phenomenon is still unclear. Here, we integrated quantitative trait locus (QTL) mapping with bulked segregant RNA sequencing (BSR-seq) to identify the genetic loci associated with the multifoliolate phenotype in soybean. Using linkage mapping, ten QTLs related to the multifoliolate trait were identified. Among these, a significant and major QTL, qMF-2-1 on chromosome 2 and consistently detected across biological replicates, explained more than 10% of the phenotypic variation. Together with BSR-seq analyses, which analyzed the RILs with the highest multifoliolate frequencies and those with the lowest frequencies as two distinct bulks, two candidate genes were identified: Glyma.06G204300 encoding the transcription factor TCP5, and Glyma.06G204400 encoding LONGIFOLIA 2 (LNG2). Transcriptome analyses revealed that stress-responsive genes were significantly differentially expressed between high-multifoliolate occurrence lines and low occurrence ones, indicating environmental factors probably influence the appearance of multifoliolate leaves in soybean through stress-responsive genes. Hence, this study offers new insights into the genetic mechanism behind the multifoliolate phenotype in soybean., (© 2024. The Author(s).)

Purpose: Adenoid cystic carcinoma (AdCC) of the head and neck harbors MYB/MYBL1::NFIB fusions in around 60% of cases, with unfavorable long-term survival due to frequent recurrences and metastases, currently lacking effective targeted therapy. The study aims to identify actionable alterations and to elucidate the molecular underpinnings of MYB/MYBL1::NFIB-negative AdCC using a large targeted RNA sequencing panel., Methods and Results: We retrospectively searched our MSK-Solid Fusion clinical sequencing database for head and neck AdCC sequenced between 2016 and 2023. Of a total of 55 cases, 28 showed MYB::NFIB, 7 showed MYBL1::NFIB, and one case each harbored MYB::MPDZ (case 1) and FUS::MYB (case 2). One base of tongue tumor expressed both MYB::NFIB fusion and MET exon 14 skipping transcripts due to concurrent MET splice site mutation, D1010N (case 3). One parotid tumor lacked MYB/MYBL1 rearrangement but instead showed an in-frame SEC16A::NOTCH1 fusion that preserved the secretase cleavage site (case 4). Clinical records on 4 cases with non-canonical sequencing findings were reviewed. Distant metastases were present at the initial diagnosis (case 2) or at recurrence (cases 1, 3, and 4). Disease-related mortality occurred in cases 2 and 4 despite radiotherapy and immunotherapy., Conclusions: The study improved the understanding of AdCC providing the first documentation of tumor clinical behavior associated with MYB::MPDZ and FUS::MYB fusions and reporting potentially actionable SEC16A::NOTCH1 fusion and MET exon 14 skipping mutation. Further research is needed to explore the therapeutic utility of MET inhibition and the efficacy of γ-secretase inhibitors against rare NOTCH1 fusions in AdCC., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Neuroblastoma is the most common extracranial solid tumor in children, and a leading cause of childhood cancer deaths. All neuroblastomas arise from neural crest-derived sympathetic neuronal progenitors, but numerous mutations, the most common of which is MYCN amplification, give rise to these lesions. Epigenetic aberrations also play a role in oncogenesis and tumor progression. To better understand biologic diversity of neuroblastomas, we performed joint single-nucleus ATAC sequencing and single-nucleus RNA sequencing on six neuroblastoma cell lines, three of which are MYCN amplified. After standard filtering for high-quality nuclei, we obtained chromatin accessibility and transcript abundance data from 41,733 neuroblastoma tumor cells. Preliminary analysis reveals significant diversity in chromatin landscape and gene expression across neuroblastoma cell lines. This dataset is a valuable resource for studying the transcriptional and epigenetic mechanisms of this deadly childhood disease., (© 2024. The Author(s).)

Macrophages induce the expression of hundreds of genes in response to immune threats. However, current technology limits our ability to capture single-cell inducible gene expression dynamics. Here, we generated high-resolution time series single-cell RNA sequencing (scRNA-seq) data from mouse macrophages responding to six stimuli, and imputed ensembles of real-time single-cell gene expression trajectories (scGETs). We found that dynamic information contained in scGETs substantially contributes to macrophage stimulus-response specificity (SRS). Dynamic information also identified correlations between immune response genes, indicating biological coordination. Furthermore, we showed that the microenvironmental context of polarizing cytokines profoundly affects scGETs, such that measuring response dynamics offered clearer discrimination of the polarization state of individual macrophage cells than single time-point measurements. Our findings highlight the important contribution of dynamic information contained in the single-cell expression responses of immune genes in characterizing the SRS and functional states of macrophages., Competing Interests: Declaration of interests K.M.S. and A.H. are authors of a patent (publication number 20230408493) that relates to the content of this manuscript., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Autoimmunity underlies many painful disorders, such as enthesopathies, which localize to the enthesis. From infiltration of the synovium and axial skeleton by B cells, to disturbances in the ratio of M1/M2 enthesis macrophages, to CD8 + T cell mediated inflammation, autoimmune dysregulation is becoming increasingly well characterized in enthesopathies. Tissue resident B cells, macrophages, neutrophils, and T cells have also been localized in healthy human entheses. However, the potential developmental origins, presence, and role of immune cells (ICs) in enthesis development is not known. Here, we use single-cell RNA-sequencing analysis to describe IC subtypes present in the enthesis before, during, and after mineralization, and to infer regulatory interactions between ICs and mesenchymal cells (MCs). We report the presence of nine phenotypically distinct IC subtypes, including B cells, macrophages, neutrophils, and T cells. We find that specific IC subtypes may promote MC-proliferation and differentiation, and that MCs may regulate IC phenotype and autoimmunity. Our findings suggest that bidirectional regulatory interactions between ICs and MCs may be important to enthesis mineralization, and suggest that progenitor MCs have a unique ability to limit autoimmunity during development., (© 2024. The Author(s).)

Background: Metabolic distress is often associated with heart failure with preserved ejection fraction (HFpEF) and represents a therapeutic challenge. Metabolism-induced systemic inflammation links comorbidities with HFpEF. How metabolic changes affect myocardial inflammation in the context of HFpEF is not known., Methods: We found that ApoE knockout mice fed a Western diet recapitulate many features of HFpEF. Single-cell RNA sequencing was used for expression analysis of CD45 + cardiac cells to evaluate the involvement of inflammation in diastolic dysfunction. We focused bioinformatics analysis on macrophages, obtaining high-resolution identification of subsets of these cells in the heart, enabling us to study the outcomes of metabolic distress on the cardiac macrophage infiltrate and to identify a macrophage-to-cardiomyocyte regulatory axis. To test whether a clinically relevant sodium glucose cotransporter-2 inhibitor could ameliorate the cardiac immune infiltrate profile in our model, mice were randomized to receive the sodium glucose cotransporter-2 inhibitor dapagliflozin or vehicle for 8 weeks., Results: ApoE knockout mice fed a Western diet presented with reduced diastolic function, reduced exercise tolerance, and increased pulmonary congestion associated with cardiac lipid overload and reduced polyunsaturated fatty acids. The main immune cell types infiltrating the heart included 4 subpopulations of resident and monocyte-derived macrophages, determining a proinflammatory profile exclusively in ApoE knockout-Western diet mice. Lipid overload had a direct effect on inflammatory gene activation in macrophages, mediated through endoplasmic reticulum stress pathways. Investigation of the macrophage-to-cardiomyocyte regulatory axis revealed the potential effects on cardiomyocytes of multiple inflammatory cytokines secreted by macrophages, affecting pathways such as hypertrophy, fibrosis, and autophagy. Finally, we describe an anti-inflammatory effect of sodium glucose cotransporter-2 inhibition in this model., Conclusions: Using single-cell RNA sequencing in a model of diastolic dysfunction driven by hyperlipidemia, we have determined the effects of metabolic distress on cardiac inflammatory cells, in particular on macrophages, and suggest sodium glucose cotransporter-2 inhibitors as potential therapeutic agents for the targeting of a specific phenotype of HFpEF., Competing Interests: None.

Background: Recent studies have suggested that cell death may be involved in bone loss or the resolution of inflammation in periodontitis. Immunogenic cell death (ICD), a recently identified cell death pathway, may be involved in the development of this disease., Methods: By analyzing single-cell RNA sequencing (scRNA-seq) for periodontitis and scoring gene set activity, we identified cell populations associated with ICD, which were further verified by qPCR, enzyme linked immunosorbent assay (ELISA) and immunofluorescence (IF) staining. By combining the bulk transcriptome and applying machine learning methods, we identified several potential ICD-related hub genes, which were then used to build diagnostic models. Subsequently, consensus clustering analysis was performed to identify ICD-associated subtypes, and multiple bioinformatics algorithms were used to investigate differences in immune cells and pathways between subtypes. Finally, qPCR and immunohistochemical staining were performed to validate the accuracy of the models., Results: Single-cell gene set activity analysis found that in non-immune cells, fibroblasts had a higher ICD activity score, and KEGG results showed that fibroblasts were enriched in a variety of ICD-related pathways. qPCR, Elisa and IF further verified the accuracy of the results. From the bulk transcriptome, we identified 11 differentially expressed genes (DEGs) associated with ICD, and machine learning methods further identified 5 hub genes associated with ICD. Consensus cluster analysis based on these 5 genes showed that there were differences in immune cells and immune functions among subtypes associated with ICD. Finally, qPCR and immunohistochemistry confirmed the ability of these five genes as biomarkers for the diagnosis of periodontitis., Conclusion: Fibroblasts may be the main cell source of ICD in periodontitis. Adaptive immune responses driven by ICD may be one of the pathogenesis of periodontitis. Five key genes associated with ICD (ENTPD1, TLR4, LY96, PRF1 and P2RX7) may be diagnostic biomarkers of periodontitis and future therapeutic targets., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Wu, Yin, Liang, Zhou, Hu, Yuan, Gu, Zhao, Gao, Cheng, Yang, Zhang, Wang, Sun and Shao.)

Introduction: Kinesin family member 20A (KIF20A) is essential for cell proliferation and is implicated in promoting tumor progression, but its role in hepatocellular carcinoma (HCC) remains poorly studied., Methods: Through the analysis of bulk RNA-sequencing (bulk RNA-seq) and single-cell RNA sequencing (scRNA-seq) data, the expression of KIF20A and its relationship with diagnosis, prognosis, and the immune microenvironment were examined. The association between KIF20A and the malignant progression and metastasis of HCC was confirmed through in vitro and in vivo experiments. Furthermore, patient re-staging was performed using Recursive Partitioning Analysis (RPA) to enhance clinical benefit., Results: In this study, we firstly found KIF20A was overexprerssed in HCC both by bulk RNA-seq and scRNA-seq, and then the overexpression of KIF20A significantly promoted the proliferation, invasion, and metastasis in vitro. In vivo, the overexpression of KIF20A promoted the growth and lung metastasis of HCC. Furthermore, gene set variation analysis of bulk RNA-seq and scRNA-seq revealed that KIF20A might be associated with cell cycle related signaling pathways of E2F and G2M, and overexpression of KIF20A inhibited the activity of p21 and bax, as well as shortened G2 phase. Importantly, we found that KIF20A could induce T cell exhaustion via the SPP1-CD44 axe using scRNA-seq. Additionally, KIF20A was also correlated with the expression of immune checkpoint inhibitors (ICIs), and KIF20Ahigh subgroup might be benefited from the ICIs therapy., Conclusion: KIF20A emerges as a pivotal driver of HCC progression, intricately regulating cell cycle pathways and modulating immune responses, which position KIF20A as a promising target for HCC management., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Su, Zhong, He, You, Xin, Wang and Liu.)

Aim: This study was to explore role of immune landscape and the immune cells in acute lymphoblastic leukemia (ALL) progression., Background: The most prevalent genetic alteration in childhood ALL is the ETV6-RUNX1 fusion. The increased proliferation of B progenitor cells could expedite the disease's progression due to irregularities in the cell cycle. Nevertheless, the mechanisms by which particular cell clusters influence the cell cycle and promote the advancement of ALL are still not well understood., Objective: This study was to explore role of immune landscape and the immune cells in ALL progression., Methods: Single-cell RNA sequencing (scRNA-seq) data of ETV6-RUNX1 and healthy pediatric samples obtained from GSE132509 were clustered and annotated using the Seurat package, and differentially highly expressed genes identified in each cluster were analyzed using DAVID for pathway annotation. Chromosome amplification and deletion were analyzed using the inferCNV package. SCENIC evaluated the regulation of transcription factors and target gene formation in cells. cellphoneDB and CellChat were served to infer ligand-receptor pairs that mediate interactions between subpopulations. The role of the target gene in regulating ALL progression was assessed using RT-qPCR, Transwell and scratch healing assays., Results: The bone marrow mononuclear cells (BMMCs) from ETV6-RUNX1 and healthy pediatric samples in GSE132509 were divided into 11 clusters, and B cell cluster 1 was identified as B progenitor cell, which was amplified on chromosome 6p. B progenitor cells were divided into seven clusters. Expression levels of amplified genes in chromosome 6p of B progenitor cell cluster 5 were the highest, and its specific highly expressed genes were annotated to pathways promoting cell cycle progression. Regulons formed in B progenitor cell cluster 5 were all involved in promoting cell cycle progression, so it was regarded as the B progenitor cell cluster that drives cell cycle progression. The key regulator of the B progenitor cell is E2F1, which promotes the migration and invasion ability of the cell line HAP1. The major ligand-receptor pairs that mediate the communication of B progenitor cell cluster 5 with cytotoxic NK/T cells or naive T cells included FAM3C-CLEC2D, CD47-SIRPG, HLAE-KLRC2, and CD47-KLRC2. HLAE-KLRC1 and TGFB1-(TGFBR1+TGFBR2)., Conclusion: This study outlined the immune cell landscape of ETV6-RUNX1 ALL and identified chromosome 6p amplification in B progenitor cells, described the major B progenitor cell cluster driving cell cycle progression and its potential regulatory mechanisms on NK cells and T cells, providing cellular and molecular insights into ETV6-RUNX1 ALL., Competing Interests: The authors declare that they have no competing interests., (© 2024 Qu et al.)

Although metastatic disease is the leading cause of cancer-related deaths, its tumor microenvironment remains poorly characterized due to technical and biospecimen limitations. In this study, we assembled a multi-modal spatial and cellular map of 67 tumor biopsies from 60 patients with metastatic breast cancer across diverse clinicopathological features and nine anatomic sites with detailed clinical annotations. We combined single-cell or single-nucleus RNA sequencing for all biopsies with a panel of four spatial expression assays (Slide-seq, MERFISH, ExSeq and CODEX) and H&E staining of consecutive serial sections from up to 15 of these biopsies. We leveraged the coupled measurements to provide reference points for the utility and integration of different experimental techniques and used them to assess variability in cell type composition and expression as well as emerging spatial expression characteristics across clinicopathological and methodological diversity. Finally, we assessed spatial expression and co-localization features of macrophage populations, characterized three distinct spatial phenotypes of epithelial-to-mesenchymal transition and identified expression programs associated with local T cell infiltration versus exclusion, showcasing the potential of clinically relevant discovery in such maps., Competing Interests: Competing interests A. Regev, J.K. and D.L.A. are co-inventors on patent application number 17/156,392, filed by the Broad Institute, for inventions relating to work in this manuscript. A. Regev is a co-founder and equity holder of Celsius Therapeutics, an equity holder in Immunitas and, until 31 July 2020, a scientific advisory board (SAB) member of Thermo Fisher Scientific, Syros Pharmaceuticals, Neogene Therapeutics and Asimov. From 1 August 2020, A. Regev is an employee of Genentech and has equity in Roche. A. Rotem is an employee of AstraZeneca since September 2020 and has equity in AstraZeneca. E.S.B. co-founded a company that is exploring commercial applications of expansion microscopy. G.N. holds equity in and consults for Akoya Biosciences. J.G.T.Z. owns stocks in the biotechnology exchange-traded funds CNCR, IDNA, IBB and XBI, owns stock in Novo Nordisk and owns stock in Adaptive Biotechnologies, 2seventy bio and bluebird bio. J.J-V. is a contractor at Genentech since 2021. M.S. is a contractor at Genentech since November 2020. O.R.-R. is a co-inventor on patent applications filed by the Broad Institute for inventions related to single-cell genomics. She has given numerous lectures on the subject of single-cell genomics to a wide variety of audiences and, in some cases, has received remuneration to cover time and costs. O.R.-R. is an employee of Genentech since 19 October 2020 and has equity in Roche. S.J.R. is a member of the SAB of Immunitas Therapeutics and receives research funding from Bristol Myers Squibb and Kite/Gilead. X.Z. is a co-founder of and consultant for Vizgen. The other authors declare no competing interests., (© 2024. The Author(s).)

Women typically live longer than men, and constitute the majority of centenarians. We applied RNA-sequencing (RNA-seq) of blood-derived lymphoblastoid cell lines (LCLs) from women aged 60-80 years and centenarians (100-105 years), validated the RNA-seq findings by real-time PCR, and additionally measured the differentially expressed genes in LCLs from young women aged 20-35 years. Top RNA-seq genes with differential expression between the age groups included three selenoproteins (GPX1, SELENOW, SELENOH) and three heat shock proteins (HSPA6, HSPA1A, HSPA1B), with the highest expression in LCLs from young women, indicating that young women are better protected from oxidative stress. The expression of two additional genes, BID encoding BH3-interacting domain death agonist and CD99 encoding CD99 antigen, showed unique age dependence, with similar expression levels in young and centenarian women while exhibiting higher and lower expression levels, respectively, in LCLs from women aged 60-80 years compared with the two other age groups. This age-related differential expression of BID and CD99 suggests elevated inflammation susceptibility in middle-aged women compared with either young or centenarian women. Our findings, once validated with human peripheral blood mononuclear cells and further cell types, may lead to novel healthy aging diagnostics and therapeutics., (© 2024 The Author(s). Drug Development Research published by Wiley Periodicals LLC.)

The growth of the ovarian antral follicle is a complex process that is difficult to study, especially in human and nonhuman primates. Understanding the antral stage of development is key to new approaches to regulating reproduction. This study analyzed cohorts of three sizes of developing antral follicles obtained from adult rhesus macaque females using RNA sequencing of oocytes and cumulus and granulosa cells. The overall objective of this study was to identify key developmental changes in gene expression in oocytes, granulosa, and cumulus cells, as nonhuman primate antral stage follicles transition through progressively larger sizes in the absence of exogenous hormonal stimulation. Only a relatively small number of genes displayed altered mRNA expression levels in any of the three cell types during this period. Most of the identified differentially expressed genes (DEGs) decreased in the granulosa cells or increased in the cumulus cells. Although the number of DEGs observed was small, these DEGs indicate predicted effects on distinct upstream regulators in the cumulus and granulosa cells. This study is particularly important because it shows for the first time the gene expression changes during antral follicle growth in a medically relevant model. NEW & NOTEWORTHY Changes in gene expression in oocytes, granulosa, and cumulus cells were determined in nonhuman primate antral stage ovarian follicles transitioning through progressively larger sizes without exogenous hormonal stimulation. Only a small number of genes displayed altered mRNA expression levels in any of the three cell types. Most of the differentially expressed genes (DEGs) decreased in granulosa cells or increased in cumulus cells. These results identified upstream regulators of antral follicle development.

This study aimed to evaluate the efficacy and therapeutic mechanism of parthenolide (PTL) in breast cancer (BC) through a comprehensive strategy integrating network pharmacology, single-cell RNA sequencing (scRNA-seq) and metabolomics. In network pharmacology, 70 therapeutic targets were identified, of which 16 core targets were filtered out through seven classical algorithms of Cytohubba plugin. Additionally, the hub module of PPI network was extracted using MCODE plugin. Molecular docking and molecular dynamics simulation showed a potent binding affinity between PTL and JNK, subsequently validated by MST and SPR assays. Further, Mendelian randomization analysis indicated that JNK was causally associated with BC. GO and KEGG enrichment analyses revealed that PTL counteracted BC via promoting ROS generation, inducing apoptosis and suppressing proliferation, which potentially involved the coordinated regulation of MAPK and FoxO1 pathways. Moreover, ssGSEA and scRNA-seq analysis suggested that PTL may act on T cell immune microenvironment of BC. Subsequently, these bioinformatics-based predictions were experimentally validated using in-vitro and in-vivo models. Finally, metabolome profiling unveiled that PTL remodeled the glycine, serine and threonine metabolism as well as biosynthesis of unsaturated fatty acids, and thereby contributed to BC inhibition. From molecular, immune and metabolic perspectives, this study not only provided a unique insight into the mechanistic details of PTL against BC, but also proposed a novel promising therapeutic strategy for BC., Competing Interests: Declaration of competing interest The authors have no competing interest., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Techniques for triggering neural differentiation of embryonic and induced pluripotent stem cells into neural stem cells and neurons have been established. However, neural induction of mesenchymal stem cells, including dental pulp stem cells (DPSCs), has been assessed primarily based on neural-related gene regulation, and detailed studies into the characteristics and differentiation status of cells are lacking. Therefore, this study was aimed at evaluating the cellular components and differentiation pathways of neural lineage cells obtained via neural induction of human DPSCs. Human DPSCs were induced to neural cells in monolayer culture and examined for gene expression and mechanisms underlying differentiation using microarray-based ingenuity pathway analysis. In addition, the neural lineage cells were subjected to single-cell RNA sequencing (scRNA-seq) to classify cell populations based on gene expression profiles and to elucidate their differentiation pathways. Ingenuity pathway analysis revealed that genes exhibiting marked overexpression, post-neuronal induction, such as FABP7 and ZIC1, were associated with neurogenesis. Furthermore, in canonical pathway analysis, axon guidance signals demonstrated maximum activation. The scRNA-seq and cell type annotations revealed the presence of neural progenitor cells, astrocytes, neurons, and a small number of non-neural lineage cells. Moreover, trajectory and pseudotime analyses demonstrated that the neural progenitor cells initially engendered neurons, which subsequently differentiated into astrocytes. This result indicates that the aforementioned neural induction strategy generated neural stem/progenitor cells from DPSCs, which might differentiate and proliferate to constitute neural lineage cells. Therefore, neural induction of DPSCs may present an alternative approach to pluripotent stem cell-based therapeutic interventions for nervous system disorders., (© 2024. The Author(s) under exclusive licence to Japan Human Cell Society.)

Although next-generation sequencing has enabled diagnoses for many patients with Mendelian disorders, the majority remain undiagnosed. Here, we present a sibling pair who were clinically diagnosed with Escobar syndrome, however targeted gene testing was negative. Exome sequencing (ES), and later genome sequencing (GS), revealed compound heterozygous TTN variants in both siblings, a maternally inherited frameshift variant [(NM_133378.4):c.36812del; p.(Asp12271Valfs*10)], and a paternally inherited missense variant [(NM_133378.4):c.12322G > A; p.(Asp4108Asn)]. This result was considered nondiagnostic due to poor clinical fit and limited pathogenicity evidence for the missense variant of uncertain significance (VUS). Following initial nondiagnostic RNA sequencing (RNAseq) on muscle and further pursuit of other variants detected on the ES/GS, a reanalysis of noncanonical splice sites in the muscle transcriptome identified an out-of-frame exon retraction in TTN, near the known VUS. Interim literature included reports of patients with similar TTN variants who had phenotypic concordance with the siblings, and a diagnosis of a congenital titinopathy was given 4 years after the TTN variants had been initially reported. This report highlights the value of reanalysis of RNAseq with a different approach, expands the phenotypic spectrum of congenital titinopathy and also illustrates how a perceived phenotypic mismatch, and failure to consider known variants, can result in a prolongation of the diagnostic journey., (© 2024 Wiley Periodicals LLC.)

Purpose: Breast cancer is a common malignancy in women, and its metastasis is a leading cause of cancer-related deaths. Single-cell RNA sequencing (scRNA-seq) can distinguish the molecular characteristics of metastasis and identify predictor genes for patient prognosis. This article explores gene expression in primary breast cancer tumor tissue against metastatic cells in the lymph node and liver using scRNA-seq., Methods: Breast cancer scRNA-seq data from the Gene Expression Omnibus were used. The data were processed using R and the Seurat package. The cells were clustered and identified using Metascape. InferCNV is used to analyze the variation in copy number. Differential expression analysis was conducted for the cancer cells using Seurat and was enriched using Metascape., Results: We identified 18 distinct cell clusters, 6 of which were epithelial. CNV analysis identified significant alterations with duplication of chromosomes 1, 8, and 19. Differential gene analysis resulted in 17 upregulated and 171 downregulated genes for the primary tumor in the primary tumor vs. liver metastasis comparison and 43 upregulated and 4 downregulated genes in the primary tumor in the primary tumor vs lymph node metastasis comparison. Several enriched terms include Ribosome biogenesis, NTP synthesis, Epithelial dedifferentiation, Autophagy, and genes associated with epithelial-to-mesenchymal transitions., Conclusion: No single gene or pathway can clearly explain the mechanisms behind tumor metastasis. Several mechanisms contribute to lymph node and liver metastasis, such as the loss of differentiation, epithelial-to-mesenchymal transition, and autophagy. These findings necessitate further study of metastatic tissue for effective drug development., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Objective: Synovial fluid (SF)-derived T cells are frequently studied as a proxy for investigating the synovial tissue (ST) T cell infiltrate in inflammatory arthritis. However, because ST is the primary site of inflammatory activity, there is debate as to whether SF provides a true reflection of the ST T cell population., Methods: In this study, we used single-cell RNA sequencing paired with single-cell T cell receptor (TCR) sequencing to directly compare memory T cells from paired samples of SF and ST from six patients with inflammatory arthritis to investigate their similarity in terms of TCR repertoire and T cell subset composition., Results: The TCR repertoires of SF and ST T cells were strikingly similar, particularly for CD8+ T cells. A median of 49% of the total CD8+ TCR repertoire in SF was shared with ST, compared with 20% shared with blood. Similarly, 47% of the ST CD8+ TCR repertoire was shared with SF compared to 25% with blood. Furthermore, once the effect of collagenase digestion on gene expression by ST T cells had been accounted for, the frequencies of specific CD8+ and CD4+ T cell subsets were, in general, similar in SF and ST and were distinct from blood., Conclusion: Our results suggest that T cells migrate and equilibrate between the SF and ST and maintain similar phenotypes in both sites. We conclude that SF is an appropriate proxy for investigating the T cell infiltrate in inflamed synovium, particularly in terms of investigating the TCR repertoire., (© 2024 The Author(s). Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

Idiopathic carpal tunnel syndrome is the most common entrapment neuropathy in hand surgery, and it is characterized by Noninflammatory fibrosis of subsynovial connective tissues. The prevalence and incidence differ between male and female individuals, and the mechanism underlying this difference remains largely unclear. In the present study, we collected subsynovial connective tissues from six male and six female patients diagnosed with idiopathic carpal tunnel syndrome during surgery. We performed a comprehensive gene expression analysis using RNA sequencing to compare the gene expression profiles between male and female patients with idiopathic carpal tunnel syndrome. We identified 26 genes with significantly different expressions between male and female patients, in which POSTN, COL1A1, and COL3A1, which are involved in extracellular matrix organization, and IGF1, an important fibrotic factor, were significantly upregulated in male patients. Immunohistochemistry confirmed the expression of proteins encoded by these genes in tissues, and male patients tended to show increased POSTN expression. Our results indicate that fibrosis of subsynovial connective tissues is induced by different mechanisms in male and female patients, and genes involved in extracellular matrix organization, especially POSTN, might be important factors in male patients. This study provides insight into the pathogenesis of idiopathic carpal syndrome and might contribute to the development of new treatment strategies., (© 2024 Orthopaedic Research Society.)

Factors that determine nonresponse to immune checkpoint inhibitor (ICI) remain unclear. The protumor activities of cancer-associated fibroblasts (CAFs) suggest that they are potential therapeutic targets for cancer treatment. There is, however, a lack of CAF-related signature in predicting response to immunotherapy in gastric cancer (GC). Single-cell RNA sequencing (scRNA-seq) and RNA sequencing (RNA-seq) data of GC immunotherapy were downloaded from the Gene Expression Omnibus database. Bulk RNA-seq data were obtained from The Cancer Genome Atlas. The R package 'Seurat' was used for scRNA-seq data processing. Cellular infiltration, receptor-ligand interactions, and evolutionary trajectory analysis were further explored. Differentially expressed genes affecting overall survival were obtained using the limma package. Weighted Gene Correlation Network Analysis was used to identify key modules of immunotherapy nonresponder. Prognostic model was constructed by univariate Cox and least absolute contraction and selection operator analysis using the intersection of activated fibroblast genes (AFGs) with key module genes. The differences in clinicopathological features, immune microenvironment, immunotherapy prediction, and sensitivity to small molecule agents between the high- and low-risk groups were further investigated. Based on scRNA-seq, we finally identified 20 AFGs associations with the prognosis of GC patients. AFGs' high expression levels were correlated with both poor prognosis and tumor progression. Three genes ( FRZB , SPARC , and FKBP10 ) were identified as immunotherapy nonresponse-related fibroblast genes and used to construct the prognostic signature. This signature is an independent significant risk factor affecting the clinical outcomes of GC patients. Remarkably, there were more CD4 memory T cells, resting mast cells, and M2 macrophages infiltrating in the high-risk group, which was characterized by higher tumor immune exclusion. Moreover, patients with higher risk scores were more prone to not respond to immunotherapy but were more sensitive to various small molecule agents, such as memantine. In conclusion, this study constructed a fibroblast-associated ICI nonresponse gene signature, which could predict the response to immunotherapy. This study potentially revealed a novel way to overcome immune resistance in GC., (Copyright © 2024 Wolters Kluwer Health, Inc. All rights reserved.)

The intestinal helminth Ascaris lumbricoides infects over 800 million people. Infections are often chronic and immunity is not sterilizing due to host-immune modulation, therefore reinfection is common after antihelmintic treatment. We have previously demonstrated a role for Ascaris spp. extracellular vesicles (EVs) in host immune modulation but whether EVs are recognized by the adaptive immune system and are present systemically in the host remains unknown. Therefore, we employed a well-established trickle infection model in pigs to mimic natural Ascaris infection in humans. EVs were isolated from adult Ascaris suum followed by immunoblotting of EV and EV-depleted secretory fractions using plasma from infected and uninfected pigs. Next, EVs were isolated from pig plasma at day 56 post first infection and subjected to deep small RNAseq analysis. RNAs were aligned to A. suum and Sus scrofa miRNA complements to detect A. suum EVs and elucidate the host EV micro RNA (miRNA) response to infection, respectively. Infection generates robust antibody responses against A. suum EVs that is distinct from EV-depleted fractions. However, A. suum miRNAs were not detectable in EVs from the peripheral blood. Notably, host plasma-derived EV miRNA profiles showed significant changes between infected and uninfected pigs, indicating that Ascaris infection drives systemic changes in host EV composition., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Forensic genetic analyses aim to retrieve as much information as possible from biological trace material recovered from crime scenes. While standard short tandem repeat (STR) profiling is essential to individualize biological traces, its significance is diminished in crime scenarios where the presence of a suspect's DNA is acknowledged by all parties. In such cases, forensic (m)RNA analysis can provide crucial contextualizing information on the source level about a trace's composition, i.e., body fluids/tissues, and has therefore emerged as a powerful tool for modern forensic investigations. However, the question which of several suspects contributed a specific component (body fluid) to a mixed trace cannot be answered by RNA analysis using conventional methods. This individualizing information is stored within the sequence of the mRNA transcripts. Massively parallel sequencing (MPS) represents a promising alternative, offering not only higher multiplex capacity, but also the typing of individual coding region SNPs (cSNPs) to enable the assignment of contributors to mixture components, thereby reducing the risk of association fallacies. Herein, we describe the development of an extensive mRNA/cSNP panel for targeted sequencing on the IonTorrent S5 platform. Our panel comprises 30 markers for the detection of six body fluids/tissues (blood, saliva, semen, skin, vaginal and menstrual secretion), along with 70 linkage-controlled cSNPs for contributor assignment. It exhibited high reliable detection sensitivity with RNA inputs down to 0.75 ng and a conservatively calculated probability of identity of 0.03 - 6 % for individual body fluid-specific cSNP profiles. Limitations and areas for future work include RNA-related allele imbalances, inclusion of markers to correctly identify rectal mucosa and the optimization of specific markers. In summary, our new panel is intended to be a major step forward to interpret biological evidence at sub-source and source level based on cSNP attribution of a body fluid component to a suspect and victim, respectively., Competing Interests: Declaration of Competing Interest All authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

The occurrence of poisoning incidents caused by cyanobacterial blooms has aroused wide public concern. Microcystin-leucine arginine (MC-LR) is a well-established toxin produced by cyanobacterial blooms, which is widely distributed in eutrophic waters. MC-LR is not only hazardous to the water environment but also exerts multiple toxic effects including liver toxicity in both humans and animals. However, the underlying mechanisms of MC-LR-induced liver toxicity are unclear. Herein, we used advanced single-cell RNA sequencing technology to characterize MC-LR-induced liver injury in mice. We established the first single-cell atlas of mouse livers in response to MC-LR. Our results showed that the differentially expressed genes and pathways in diverse cell types of liver tissues of mice treated with MC-LR are highly heterogeneous. Deep analysis showed that MC-LR induced an increase in a subpopulation of hepatocytes that highly express Gstm3, which potentially contributed to hepatocyte apoptosis in response to MC-LR. Moreover, MC-LR increased the proportion and multiple subtypes of Kupffer cells with M1 phenotypes and highly expressed proinflammatory genes. Furthermore, the MC-LR increased several subtypes of CD8 + T cells with highly expressed multiple cytokines and chemokines. Overall, apart from directly inducing hepatocytes apoptosis, MC-LR activated proinflammatory Kupffer cell and CD8 + T cells, and their interaction may constitute a hostile microenvironment that contributes to liver injury. Our findings not only present novel insight into underlying molecular mechanisms but also provide a valuable resource and foundation for additional discovery of MC-LR-induced liver toxicity., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article., (Copyright © 2023. Published by Elsevier B.V.)

The pathogenesis of sepsis, with a high mortality rate and often poor prognosis, has not been fully elucidated. Therefore, an in-depth study on the pathogenesis of sepsis at the molecular level is essential to identify key sepsis-related genes. The aim of this study was to explore the key genes and potential molecular mechanisms of sepsis using a bioinformatics approach. In addition, key genes with miRNA network correlation analysis and immune infiltration correlation analysis were investigated. The scRNA dataset (GSE167363) and RNA-seq dataset (GSE65682, GSE134347) from GEO database were used for screening out differentially expressed genes using single-cell sequencing and transcriptome sequencing. The analysis of immune infiltration was evaluated by the CIBERSORT method. Key genes and possible mechanisms were identified by WGCNA analysis, GSVA analysis, GSEA enrichment analysis and regulatory network analysis, and miRNA networks associated with key genes were constructed. Nine key genes associated with the development of sepsis, namely IL7R, CD3D, IL32, GPR183, HLA-DPB1, CD81, PEBP1, NCL, and ETS1 were screened, and the specific signaling mechanisms associated with the key genes causing sepsis were predicted. Immune profiling showed immune heterogeneity between control and sepsis samples. A regulatory network of 82 miRNAs, 266 pairs of mRNA-miRNA relationship pairs was also constructed. These nine key genes have the potential to become biomarkers for the diagnosis of sepsis and provide new targets and research directions for the treatment of sepsis.

Tri (2-Ethylhexyl) phosphate (TEHP), widely used as a fire retardant and plasticizer, has been commonly found in the environment. Its potential health-related risks, especially reproductive toxicity, have aroused concern. However, the potential cellular mechanisms remain unexplored. In this study, we aimed to investigate the molecular mechanisms underlying TEHP-caused cell damage in Sertoli cells, which play a crucial role in supporting spermatogenesis. Our findings indicate that TEHP induces apoptosis in 15P-1 mouse Sertoli cells. Subsequently, we conducted RNA sequencing analyses, which suggested that ER stress, autophagy, and MAPK-related pathways may participate in TEHP-induced cytotoxicity. Furthermore, we demonstrated that TEHP triggers ER stress, activates p38 MAPK, and inhibits autophagy flux. Then, we showed that the inhibition of ER stress or p38 MAPK activation attenuates TEHP-induced apoptosis, while the inhibition of autophagy flux is responsible for TEHP-induced apoptosis. These results collectively reveal that TEHP induces ER stress, activates p38, and inhibits autophagy flux, ultimately leading to apoptosis in Sertoli cells. These shed light on the molecular mechanisms underlying TEHP-associated testicular toxicity., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Single-cell RNA sequencing (scRNA-seq) has facilitated a deeper comprehension of the molecular mechanisms behind eye diseases and has prompted the selection of precise therapeutic targets by examining the cellular and molecular intricacies at the single-cell level. This review delineates the pivotal role of scRNA-seq in elucidating the functions of innate immune cells within the context of ocular pathologies. Recent advancements in scRNA-seq have revealed that innate immune cells, both from the periphery and resident in the retina, are actively engaged in various stages of multiple eye diseases. Notably, resident microglia and infiltrating neutrophils exhibit swift responses during the initial phase of injury, while peripheral monocyte-derived macrophages exhibit transcriptomic profiles akin to those of activated microglia, suggesting their potential for long-term residence within the retina. The scRNA-seq analyses have underscored the cellular heterogeneity and gene expression alterations within innate immune cells, which, while sharing commonalities, exhibit disease-specific variations. These insights have not only broadened our understanding of the cellular and molecular mechanisms in eye diseases but also paved the way for the identification of candidate targets for targeted therapeutic interventions. The application of scRNA-seq technology has heralded a new era in the study of ocular pathologies, enabling a more detailed appreciation of the roles that innate immune cells play across a spectrum of eye diseases., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Lu, Mao and Yuan.)

Despite the scientific relevance of circular RNAs (circRNAs), the study of these RNAs in non-model organisms, especially in sheep, is still in its infancy. On the other hand, while some studies have focused on sheep circRNA identification in a limited number of tissues, there is a lack of comprehensive analysis that profile circRNA expression patterns across the tissues not yet investigated. In this study, 61 public RNA sequencing datasets from 12 different tissues were uniformly analyzed to identify circRNAs, profile their expression and investigate their various characteristics. We reported for the first time a circRNA expression landscape with functional annotation in sheep tissues not yet investigated including hippocampus, BonMarrowMacrophage, left-ventricle, thymus, ileum, reticulum and 23-day-embryo. A stringent computational pipeline was employed and 8919 exon-derived circRNAs with high confidence were identified, including 88 novel circRNAs. Tissue-specificity analysis revealed that 3059 circRNAs were tissue-specific, which were also more specific to the tissues than linear RNAs. The highest number of tissue-specific circRNAs was found in kidney, hippocampus and thymus, respectively. Co-expression analysis revealed that expression of circRNAs may not be affected by their host genes. While most of the host genes produced more than one isoform, only one isoform had dominant expression across the tissues. The host genes of the tissue-specific circRNAs were significantly enriched in biological/pathways terms linked to the important functions of their corresponding tissues, suggesting potential roles of circRNAs in modulating physiological activity of those tissues. Interestingly, functional terms related to the regulation and various signaling pathways were significantly enriched in all tissues, suggesting some common regulatory mechanisms of circRNAs to modulate the physiological functions of tissues. Finding of the present study provide a valuable resource for depicting the complexity of circRNAs expression across tissues of sheep, which can be useful for the field of sheep genomic and veterinary research., (© 2024. The Author(s).)

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive lung disease with limited effective therapies. Interstitial macrophages (IMs), especially those derived from monocytes, play an unknown role in IPF pathogenesis. By using single-cell RNA sequencing (scRNA-seq), bleomycin (BLM)-induced pulmonary fibrosis mouse lungs were analyzed to characterize the cellular landscape and heterogeneity of macrophages in this model. scRNA-seq was used to identify distinct interstitial macrophage subpopulations in fibrotic lungs, with monocyte-derived macrophages exhibiting a pro-fibrotic gene expression profile enriched in wound healing, extracellular matrix (ECM) remodeling, and pro-fibrotic cytokine production functions. A pseudotime analysis revealed that IMs originated from monocytes and differentiated along a specific trajectory. A cell-cell communication analysis demonstrated strong interactions between monocyte-derived interstitial macrophages (Mo-IMs) and fibroblasts through the transforming growth factor beta (TGFβ), secreted phosphoprotein 1 (SPP1), and platelet-derived growth factor (PDGF) signaling pathways. Flow cytometry validated the presence and expansion of Mo-IMs subpopulations in BLM-treated mice. This study reveals the cellular heterogeneity and developmental trajectory of lung macrophages in early BLM-induced pulmonary fibrosis, highlighting the crucial role of Mo-IMs with a pro-fibrotic phenotype in IPF pathogenesis via interactions with fibroblasts. Targeting these specific macrophage subpopulations and associated signaling pathways may provide novel therapeutic strategies for IPF.

Single-cell RNA sequencing has enabled the study of aging at a molecular scale. While substantial progress has been made in measuring age-related gene expression, the underlying patterns and mechanisms of aging transcriptomes remain poorly understood. To address this gap, we propose a physics-inspired, data-analysis approach to extract additional insights from single-cell RNA sequencing data. By considering the genome as a many-body interacting system, we leverage central idea of the Renormalization Group to construct an approach to hierarchically describe aging across a spectrum of scales for the gene expresion. This framework provides a quantitative language to study the multiscale patterns of aging transcriptomes. Overall, our study demonstrates the value of leveraging theoretical physics concepts like the Renormalization Group to gain new biological insights from complex high-dimensional single-cell data., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Li, Mani. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The quality of Recombination signal sequences (RSSs), location, and genetics of mammalian V, D, and J genes synergistically affect the recombination frequency of genes; however, the specific regulatory mechanism and efficiency have not been elucidated. By taking advantage of single-cell RNA-sequencing (scRNA-seq) and high-throughput sequencing (HTS) to investigate V(D)J rearrangement characteristics in the CDR3 repertoire, we found that the distal and proximal V genes (or J genes) "to D" gene were involved in rearrangement significantly more frequently than the middle V genes (or J genes) in the TRB locus among various species, including Primates (human and rhesus monkey), Rodentia (BALB/c, C57BL/6, and Kunming mice), Artiodactyla (buffalo), and Chiroptera (Rhinolophus affinis). The RSS quality of the V and J genes affected their frequency in rearrangement to varying degrees, especially when the V-RSSs with recombination signal information content (RIC) score < -45 significantly reduced the recombination frequency of the V gene. The V and J genes that were "away from D" had the dual advantages of recombinant structural accessibility and relatively high-quality RSSs, which promoted their preferential utilization in rearrangement. The quality of J-RSSs formed during mammalian evolution was apparently greater than that of V-RSSs, and the D-J distance was obviously shorter than that of V-D, which may be one of the reasons for guaranteeing that the "D-to-J preceding V-to-DJ rule" occurred when rearranged. This study provides a novel perspective on the mechanism and efficiency of V-D-J rearrangement in the mammalian TRB locus, as well as the biased utilization characteristics and application of V and J genes in the initial CDR3 repertoire., (© 2024. The Author(s).)

Objective: This study aimed to explore the heterogeneity of tumor endothelial cells (TECs) in hepatocellular carcinoma (HCC) and their role in tumor progression, with the goal of identifying new therapeutic targets and strategies to improve patient prognosis., Methods: Single-cell RNA sequencing data from nine primary liver cancer samples were analyzed, obtained from the Gene Expression Omnibus (GEO) database. Data preprocessing, normalization, dimensionality reduction, and batch effect correction were performed based on the Seurat package. HCC cell types were identified using uniform manifold approximation and projection (UMAP) and cluster analysis, and the different cell types were annotated using the CellMarker database. Pseudotime trajectory analysis was conducted with Monocle to explore the differentiation trajectory of TECs. MAPK signaling pathway activity and copy number variations (CNV) in TECs were analyzed in conjunction with data from The Cancer Genome Atlas (TCGA), the trans-well and wound healing assay was used for cell invasion and migration activity assessment., Results: Two subgroups of TECs (TECs 1 and TECs 2) were identified, exhibiting distinct functional activities and signaling pathways. Specifically, TECs 1 may be involved in tumor cell proliferation and inflammatory responses, whereas TECs 2 is not only involved in cell proliferation pathways, but also enriched in pathways such as metabolic synthesis. Pseudotime analysis revealed dynamic changes in TECs subgroups during HCC progression, correlating specific gene expressions (such as PDGFRB, PGF, JUN, and NR4A1). Subsequently, the JUN gene was predicted by performing binding sites and was shown to act as a transcription factor that may regulate the expression of the PGF gene. CNV analysis highlighted key genes and pathways in TECs that might influence HCC progression, and the PGF as key regulatory factor mediated cell proliferation and migration., Conclusion: The study revealed the heterogeneity of TECs in HCC and their potential roles in tumor progression, offering new perspectives and potential therapeutic targets for HCC molecular mechanisms. The findings emphasize the importance of further exploring TECs heterogeneity for understanding HCC pathogenesis and developing personalized treatment strategies., Competing Interests: The authors declare that they have no competing interests., (© 2024 Sun et al.)

Ethnopharmacological Relevance: Nonalcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver disease worldwide. Nonalcoholic steatohepatitis (NASH) is a crucial component of this disease spectrum. The Yanxiao Di'naer formula (YXDNE) is an Uyghur medical extract that has been used in folk medicine to treat hepatitis for a long time. However, the role and mechanism of action of YXDNE in NASH treatment remains unclear., Objective: The objective of this study was to assess the effectiveness of YXDNE in treating NASH induced by injections of carbon tetrachloride combined with a high-fat high-cholesterol diet (HFHCD), and to clarify the underlying mechanisms., Methods: The compounds in the YXDNE extract were analysed for classification and proportions using ultra-performance liquid chromatography-mass spectrometry. The efficacy of YXDNE in treating abnormal lipid metabolism was evaluated in L02 cells in vitro. In addition, a C57BL/6 mouse model of NASH was established to evaluate the therapeutic efficacy of YXDNE in vivo. Metabolomics and RNA sequencing were used to analyse the therapeutic effects of YXDNE on the liver. The corresponding signalling pathways were found to target AMPKα1, PPARα, and NF-κB. The efficacy of YXDNE was validated using inhibitors or silencing RNA (siRNA) against AMPKα1 and PPARα., Results: This study confirmed that YXDNE treatment ameliorated NASH in a murine model of this disease. Metabolomics analysis suggested that YXDNE efficacy was associated with fatty acid catabolism and AMPK signalling pathways. RNA sequencing results showed that YXDNE efficacy in treating NASH was highly correlated with the AMPK, PPARα and NF-κB pathways. Both in vitro and in vivo experimental data demonstrated that YXDNE affected the expression of p-AMPKα1, PPARα, p-NF-κB, IκB, and p-IκB. The efficacy of YXDNE in treating NASH in vitro was cancelled when AMPK was inhibited with Compound C or PPARα was modulated via siRNA., Conclusions: YXDNE may have a therapeutic effect on abnormal lipid metabolism in L02 cells and in a murine model of NASH by affecting the AMPKα1/PPARα/NF-κB signalling pathway. Therefore, YXDNE has the potential for clinical application in the prevention and treatment of NASH., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Background: Chronic obstructive pulmonary disease (COPD) is a complex respiratory condition influenced by environmental and genetic factors. Using next-generation sequencing, we aimed to identify dysregulated genes and potential therapeutic targets for COPD., Methods: Peripheral blood leukocyte RNA profiles from COPD patients and healthy controls were analyzed using next-generation sequencing. Key genes involved in COPD pathogenesis were identified through protein-protein interaction network analysis. In vitro, bronchial epithelial cells treated with cigarette smoke extract (CSE) were used to study the effects on gene expression, cell viability, apoptosis, and ferroptosis. Additionally, Lipocalin 2 (LCN2) inhibition experiments were conducted to elucidate its role in COPD-related cellular processes., Results: Analysis of RNA profiles revealed consistent downregulation of 17 genes and upregulation of 21 genes across all COPD groups. Among these, Cathelicidin Antimicrobial Peptide(CAMP), Defensin Alpha 4(DEFA4), Neutrophil Elastase(ELANE), LCN2 and Lactotransferrin(LTF) were identified as potentially important players in COPD pathogenesis. Particularly, LCN2 exhibited a close association with COPD and was found to be involved in cellular processes. In vitro experiments demonstrated that CSE treatment significantly increased LCN2 expression in bronchial epithelial cells in a concentration-dependent manner. Moreover, CSE-induced apoptosis and ferroptosis were observed, along with alterations in cell viability, Glutathione content, Fe2 + accumulation, ROS: Reactive Oxygen Species and Malondialdehyde levels, Lactate Dehydrogenase(LDH) release and Glutathione Peroxidase 4(GPX4) expression. Inhibition of LCN2 expression partially reversed these effects, indicating the pivotal role of LCN2 in COPD-related cellular processes., Conclusion: Our study identified six candidate genes: CAMP, DEFA4, ELANE, LCN2, and LTF were upregulated, HSPA1B was downregulated. Notably, LCN2 emerges as a significant biomarker in COPD pathogenesis, exerting its effects by promoting apoptosis and ferroptosis in bronchial epithelial cells., (© 2024. The Author(s).)

Purpose: MYB has been shown to play a central role in oncogenesis in a majority of adenoid cystic carcinomas (ACC). Testing for MYB expression via immunohistochemistry (IHC) or testing for the MYB gene fusion by next-generation sequencing (NGS) have become useful tools for the diagnosis of ACC. In addition, detection of MYB expression may have implications for patient management., Methods: A cohort of 35 ACC cases was identified from the archival pathology files of the Massachusetts General Hospital. Cases were tested for MYB expression using a panel of 4 different commercially available MYB antibodies and scored using a modified Allred system. RNA-based NGS for MYB gene fusion detection was also performed., Results: Among 4 different MYB antibodies, the sensitivity for MYB detection ranged from 26 to 97%. When a 30% threshold for determination of MYB immunohistochemical positivity was used, the AB&#95;10900735 IHC clone showed the maximum sensitivity (97%). RNA sequencing revealed 19 (54%) cases positive for MYB fusions, and expression analysis derived from the sequencing data confirmed a significant association between MYB expression and fusion status (p = 0.036). Although less sensitive, the AB&#95;778878 MYB clone showed a significant positive association between IHC staining and MYB RNA expression (R 2  = 0.15, p = 0.023)., Conclusion: The detection of MYB expression using immunohistochemistry varies significantly depending on the antibody used. Comparison with MYB fusion and transcription levels, as determined by NGS, reveals that MYB has a complex relationship between genetic alterations, transcript levels, and protein abundance., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

The mammalian cerebral cortex is composed of a rich diversity of cell types. Sensory cortical cells are organized into networks that rely on their functional diversity to ultimately carry out a variety of sophisticated cognitive functions for perception, learning, and memory. The auditory cortex (AC) has been most extensively studied for its experience-dependent effects, including for perceptual learning and associative memory. Here, we used single-nucleus RNA sequencing (snRNA-seq) in the AC of the adult rat to investigate the breadth of transcriptionally diverse cell types that likely support the role of AC in experience-dependent functions. A variety of unique excitatory and inhibitory neuron subtypes were identified that harbor unique transcriptional profiles of genes with putative relevance for the adaptive neuroplasticity of cortical microcircuits. In addition, we report for the first time a diversity of astrocytes in AC that may represent functionally unique subtypes, including those that could integrate experience-dependent adult neuroplasticity at cortical synapses. Together, these results pave the way for building models of how cortical neurons work in concert with astrocytes to fulfill dynamic and experience-dependent cognitive functions., (© 2024. The Author(s).)

Background: Ovarian cancer (OC) is one of the cancers with a high incidence at present, which poses a severe threat to women's health. This study focused on identifying the heterogeneity among malignant epithelial cell OC and constructing an effective prognostic signature to predict prognosis and immunotherapy according to a multidisciplinary study., Methods: The InterCNV algorithm was used to identify the heterogeneity of OC based on the scRNA-seq and bulk RNA-seq data. Six algorithms selected EMTscore. An effective prognostic signature was conducted using the COX and Least Absolute Shrinkage and Selection Operator (LASSO) regression algorithms. The texting datasets were used to assess the accuracy of the prognostic signature. We evaluated different immune characteristics and immunotherapy response differences among other risk groups., Results: A prognostic signature including 14 genes was established. The patients in the high-risk group have poor survival outcomes. We also found that the patients in the low-risk group have higher immune cell infiltration, enrichment of immune checkpoints, and immunotherapy response, suggesting that the patients in the low-risk group may be more sensitive to immunotherapy. Finally, the laboratory test results showed that KREMEN2 was identified as a novel biomarker and therapeutic target for OC patients., Conclusions: Our study established a GRG signature consisting of 16 genes based on the scRNA-seq and bulk RNA-seq data, which provides a new perspective on the prediction of prognosis and treatment strategy for OC., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Background: Lung cancer is a highly aggressive and prevalent disease worldwide. By the time it is first diagnosed, distant metastases have usually already occurred. Among them, the prognosis of patients with brain metastasis from lung cancer is very poor. Therefore, it is particularly important to identify the evolutionary status of tumor cells during lung cancer brain metastases and discover the underlying mechanisms of lung cancer brain metastases., Methods: In this study, we analysed three types of data: single-cell RNA sequencing, bulk RNA sequencing, and spatial transcriptome. Firstly, we identified early metastatic epithelial cell clusters (EMEC) using CNV and trajectory analysis in scRNA-seq data. Secondly, we integrated scRNA-seq and spatial transcriptome data with the help of MIA (Multimodal intersection analysis) to explore the biological characteristics of EMEC. Finally, we used bulk RNA-seq data to validate the molecular characteristics of EMEC., Result: A total of 55,763 single cells were obtained and divided into 9 cell types. In brain metastasis, we found a significantly higher proportion of epithelial cells. In addition, we identified a specific subpopulation of epithelial cells, which was named as "early metastatic epithelial cell clusters (EMEC)". It is enriched in oxidative phosphorylation, coagulation, complement. Moreover, we also found that EMEC underwent cellular communication with other immune cells through ligand-receptor pairs such as MIF-(CD74 + CXCR4) and MIF-(CD74 + CD44). Next, we validated that EMEC were associated with poor clinical prognosis using three independent external datasets. Finally, spatial transcriptome analysis revealed specificity in the spatial distribution of EMEC, which shifted from the peripheral regions to the central regions of the tumour as the depth of tumor invasion progressed., Conclusion: This study reveals the potential molecular mechanisms of lung cancer brain metastasis from both single-cell and spatial transcriptomic perspectives, providing biological insights and clinical reference value for detecting patients suffering from lung cancer brain metastasis., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Pathogenic variants in the MED13L gene are associated with the autosomal dominant MED13L syndrome, which is characterised by global developmental delay and cardiac malformations. We investigated two heterozygous MED13L variants located at the canonical donor splice site motif of exon 7: c.1009+1G>C and c.1009+5G>C. We report that in silico predictions suggested two possible outcomes: exon 7 skipping, resulting in loss of the phosphodegron motif essential for MED13L regulation, or activation of a cryptic donor site in intron 7, leading to intron retention. RNA analysis confirmed that both variants affected the exon 7 splice donor site, resulting in the retention of 73 bp of intron 7. This retention caused a frameshift and premature translation termination, consistent with haploinsufficiency. Our results highlight the importance of combining predictive and experimental approaches to understand the functional impact of splice site variants. These insights into the molecular consequences of MED13L variants provide a deeper understanding of the genetic basis of MED13L syndrome., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2024. No commercial re-use. See rights and permissions. Published by BMJ.)

The large group of negative-strand RNA viruses (NSVs) comprises many important pathogens. To identify conserved patterns in host responses, we systematically compared changes in the cellular RNA levels after infection of human hepatoma cells with nine different NSVs of different virulence degrees. RNA sequencing experiments indicated that the amount of viral RNA in host cells correlates with the number of differentially expressed host cell transcripts. Time-resolved differential gene expression analysis revealed a common set of 178 RNAs that are regulated by all NSVs analyzed. A newly developed open access web application allows downloads and visualizations of all gene expression comparisons for individual viruses over time or between several viruses. Most of the genes included in the core set of commonly differentially expressed genes (DEGs) encode proteins that serve as membrane receptors, signaling proteins and regulators of transcription. They mainly function in signal transduction and control immunity, metabolism, and cell survival. One hundred sixty-five of the DEGs encode host proteins from which 47 have already been linked to the regulation of viral infections in previous studies and 89 proteins form a complex interaction network that may function as a core hub to control NSV infections.IMPORTANCEThe infection of cells with negative-strand RNA viruses leads to the differential expression of many host cell RNAs. The differential spectrum of virus-regulated RNAs reflects a large variety of events including anti-viral responses, cell remodeling, and cell damage. Here, these virus-specific differences and similarities in the regulated RNAs were measured in a highly standardized model. A newly developed app allows interested scientists a wide range of comparisons and visualizations., Competing Interests: The authors declare no conflict of interest.

Background: Pancreatic ductal adenocarcinoma (PDAC) is one of the main types of malignant tumor of the digestive system, and patient prognosis is affected by difficulties in early diagnosis, poor treatment response, and a high postoperative recurrence rate. Carbohydrate antigen 19-9 (CA19-9) has been widely used as a biomarker for the diagnosis and postoperative follow-up of PDAC patients. Nevertheless, the production mechanism and potential role of CA19-9 in PDAC progression have not yet been elucidated., Methods: We performed single-cell RNA sequencing on six samples pathologically diagnosed as PDAC (three CA19-9-positive and three CA19-9-negative PDAC samples) and two paracarcinoma samples. We also downloaded and integrated PDAC samples (each from three CA19-9-positive and CA19-9-negative patients) from an online database. The dynamics of the proportion and potential function of each cell type were verified through immunofluorescence. Moreover, we built an in vitro coculture cellular model to confirm the potential function of CA19-9., Results: Three subtypes of cancer cells with a high ability to produce CA19-9 were identified by the markers TOP2A , AQP5 , and MUC5AC . CA19-9 production bypass was discovered on antigen-presenting cancer-associated fibroblasts (apCAFs). Importantly, the proportion of immature ficolin-1 positive (FCN1+) macrophages was high in the CA19-9-negative group, and the proportion of mature M2-like macrophages was high in the CA19-9-positive group. High proportions of these two macrophage subtypes were associated with an unfavourable clinical prognosis. Further experiments indicated that CA19-9 could facilitate the transformation of M0 macrophages into M2 macrophages in the tumor microenvironment., Conclusions: Our study described CA19-9 production at single-cell resolution and the dynamics of the immune atlas in CA19-9-positive and CA19-9-negative PDAC. CA19-9 could promote M2 polarization of macrophage in the pancreatic tumor microenvironment., (Copyright © 2024 The Chinese Medical Association, produced by Wolters Kluwer, Inc. under the CC-BY-NC-ND license.)

Background: Chronic Kidney Disease (CKD) impacts over 10% of the global population, and recent advancements in high-throughput analytical technologies are uncovering the complex physiology underlying this condition. By integrating Genome-Wide Association Studies (GWAS), RNA sequencing (RNA-seq/RNA array), and single-cell RNA sequencing (scRNA-seq) data, our study aimed to explore the genes and cell types relevant to CKD traits., Methods: GWAS summary data for end-stage renal failure (ESRD) and decreased eGFR (CKD) with or without diabetes and (micro)proteinuria were obtained from the GWAS Catalog and the UK Biobank (UKB) database. Two gene Expression Omnibus (GEO) transcriptome datasets were used to establish glomerular and tubular gene expression differences between CKD patients and healthy individuals. Two scRNA-seq datasets were utilized to obtain the expression of key genes at the single-cell level. The expression profile, differentially expressed genes (DEGs), gene-gene interaction, and pathway enrichment were analysed for these CKD risk genes., Results: A total of 779 distinct SNPs were identified from GWAS across different CKD traits, involving 681 genes. While many of these risk genes are specific to the CKD traits of renal failure, decreased eGFR, and (micro)proteinuria, they share common pathways, including extracellular matrix (ECM). ECM modeling was enriched in upregulated glomerular and tubular DEGs from CKD kidneys compared to healthy controls, with the expression of relevant collagen genes, such as COL1A2, prevalent in fibroblasts/myofibroblasts. Additionally, immune responses, including T cell differentiation, were dysregulated in CKD kidneys. The late podocyte signature gene THSD7A was enriched in podocytes but downregulated in CKD. We also highlighted that the regulated risk genes of CKD are mainly expressed in tubular cells and immune cells in the kidney., Conclusions: Our integrated analysis highlight the genes, pathways, and relevant cell types associational with the pathogenesis of kidney traits, as a basis for further mechanistic studies to understand the pathogenesis of CKD., (© 2024. The Author(s).)

It has been demonstrated that CXCL12 inhibits hair growth via CXCR4, and its neutralizing antibody (Ab) increases hair growth in alopecia areata (AA). However, the molecular mechanisms have not been fully elucidated. In the present study, we further prepared humanized CXCL12 Ab for AA treatment and investigated underlying molecular mechanisms using single-cell RNA sequencing. Subcutaneous injection of humanized CXCL12 Ab significantly delayed AA onset in mice, and dorsal skin was analyzed. T cells and dendritic cells/macrophages were increased in the AA model, but decreased after CXCL12 Ab treatment. Pseudobulk RNA sequencing identified 153 differentially expressed genes that were upregulated in AA model and downregulated after Ab treatment. Gene ontology analysis revealed that immune cell chemotaxis and cellular response to type II interferon were upregulated in AA model but downregulated after Ab treatment. We further identified key immune cell-related genes such as Ifng , Cd8a , Ccr5 , Ccl4 , Ccl5 , and Il21r , which were colocalized with Cxcr4 in T cells and regulated by CXCL12 Ab treatment. Notably, CD8+ T cells were significantly increased and activated via Jak/Stat pathway in the AA model but inactivated after CXCL12 Ab treatment. Collectively, these results indicate that humanized CXCL12 Ab is promising for AA treatment via immune modulatory effects., Competing Interests: Authors MZ and J-HS were employed by Epi Biotech Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 An, Zheng, Park, Park, Noh and Sung.)

Exchangeable aluminum (Al) ions released from acidic soils with pH < 5.5 inhibit root elongation of crops, ultimately leading to yield reduced. It is necessary to identify the quantitative trait locus (QTLs) and candidate genes that confer toxicity resistance to understand the mechanism and improve tolerance of rapeseed. In this study, an F 2 segregating population was derived from a cross between Al-tolerance inbred line FDH188 (R178) and -sensitive inbred line FDH152 (S169), and the F 2:3 were used as materials to map QTLs associated with the relative elongation of taproot (RET) under Al toxicity stress. Based on bulked segregant analysis sequencing (BSA-seq), three QTLs ( qAT-A07-1 , qAT-A07-2 , and qAT-A09-1 ) were detected as significantly associated with RET, and 656 candidate genes were screened. By combined BSA and RNA-seq analysis, 55 candidate genes showed differentially expressed, including genes encoding ABC transporter G (ABCG), zinc finger protein, NAC, ethylene-responsive transcription factor (ERF), etc. These genes were probably positive factors in coping with Al toxicity stress in rapeseed. This study provides new insight into exploring the QTLs and candidate genes' response to Al toxicity stress by combined BSA-seq and RNA-seq and is helpful to further research on the mechanism of Al resistance in rapeseed.