Search

Objective To explore the gene expression profile in paraspinal muscle degeneration (PMD) and identify key ferroptosis genes. Methods RNA sequencing was performed on paraspinal muscle tissue of 3 normal and 3 PMD patients respectively to obtain differentially expressed genes. Through protein-protein interaction (PPI) and gene functional enrichment analysis, the intersection of ferroptosis genes was identified to identify key hub genes associated with ferroptosis. The diagnostic value for PMD disease was analyzed by receiver operating characteristic (ROC) curves. Results A total of 292 differentially expressed genes were identified in PMD. Among them, 125 genes were significantly downregulated and 167 genes were significantly upregulated. Bioinformatics analysis revealed that 14 differentially expressed genes were associated with ferroptosis. Among them, ferroptosis genes MUC1, ATF3 and CDKN1A were key hub genes with good specificity and sensitivity for diagnosing PMD. Functional enrichment analysis revealed that they may mediate the occurrence and progression of PMD by regulating cell apoptosis, ferroptosis and skeletal muscle tissue development and differentiation. Conclusion Ferroptosis genes MUC1, ATF3 and CDKN1A can serve as biomarkers for diagnosing PMD, providing theoretical basis for decoding the pathological mechanism of PMD and developing new drugs. [ABSTRACT FROM AUTHOR]

ABSTRACTExtensive research has explored the role of gut microbiota in colorectal cancer (CRC). Nonetheless, metatranscriptomic studies investigating the in situ functional implications of host-microbe interactions in CRC are scarce. Therefore, we characterized the influence of CRC core pathogens and biofilms on the tumor microenvironment (TME) in 40 CRC, paired normal, and healthy tissue biopsies using fluorescence in situ hybridization (FISH) and dual-RNA sequencing. FISH revealed that Fusobacterium spp. was associated with increased bacterial biomass and inflammatory response in CRC samples. Dual-RNA sequencing demonstrated increased expression of pro-inflammatory cytokines, defensins, matrix-metalloproteases, and immunomodulatory factors in CRC samples with high bacterial activity. In addition, bacterial activity correlated with the infiltration of several immune cell subtypes, including M2 macrophages and regulatory T-cells in CRC samples. Specifically, Bacteroides fragilis and Fusobacterium nucleatum correlated with the infiltration of neutrophils and CD4+ T-cells, respectively. The collective bacterial activity/biomass appeared to exert a more significant influence on the TME than core pathogens, underscoring the intricate interplay between gut microbiota and CRC. These results emphasize how biofilms and core pathogens shape the immune phenotype and TME in CRC while highlighting the need to extend the bacterial scope beyond CRC pathogens to advance our understanding and identify treatment targets.

BACKGROUND: Preeclampsia is a serious disease of pregnancy that lacks early diagnosis methods or effective treatment, except delivery. Dysregulated uterine immune cells and spiral arteries are implicated in preeclampsia, but the mechanistic link remains unclear. METHODS: Single-cell RNA sequencing and spatial transcriptomics were used to identify immune cell subsets associated with preeclampsia. Cell-based studies and animal models including conditional knockout mice and a new preeclampsia mouse model induced by recombinant mouse galectin-9 were applied to validate the pathogenic role of a CD11chigh subpopulation of decidual macrophages (dMφ) and to determine its underlying regulatory mechanisms in preeclampsia. A retrospective preeclampsia cohort study was performed to determine the value of circulating galectin-9 in predicting preeclampsia. RESULTS: We discovered a distinct CD11chigh dMφ subset that inhibits spiral artery remodeling in preeclampsia. The proinflammatory CD11chigh dMφ exhibits perivascular enrichment in the decidua from patients with preeclampsia. We also showed that trophoblast-derived galectin-9 activates CD11chigh dMφ by means of CD44 binding to suppress spiral artery remodeling. In 3 independent preeclampsia mouse models, placental and plasma galectin-9 levels were elevated. Galectin-9 administration in mice induces preeclampsia-like phenotypes with increased CD11chigh dMφ and defective spiral arteries, whereas galectin-9 blockade or macrophage-specific CD44 deletion prevents such phenotypes. In pregnant women, increased circulating galectin-9 levels in the first trimester and at 16 to 20 gestational weeks can predict subsequent preeclampsia onset. CONCLUSIONS: These findings highlight a key role of a distinct perivascular inflammatory CD11chigh dMφ subpopulation in the pathogenesis of preeclampsia. CD11chigh dMφ activated by increased galectin-9 from trophoblasts suppresses uterine spiral artery remodeling, contributing to preeclampsia. Increased circulating galectin-9 may be a biomarker for preeclampsia prediction and intervention. [ABSTRACT FROM AUTHOR]

This article discusses the limitations of using rodent models in heart transplantation research and highlights the shift towards using human specimens for discovery research. The authors focus on a study that uses single-nucleus RNA sequencing to investigate the cellular landscape of transplanted pediatric hearts. The study identifies important candidate mediators of heart transplant failure, such as graft-infiltrating macrophages and natural killer cells. The authors also explore the dynamics of immune cells within the failing transplanted heart and highlight the potential therapeutic targets that could be investigated further. However, the study has limitations due to the small sample size and the difficulty in obtaining healthy heart transplant tissue. Overall, the article emphasizes the potential of single-cell sequencing and other molecular pathology techniques in advancing our understanding of heart transplantation. [Extracted from the article]

Objective To summarize the clinical features and treatment of genetic epilepsy with febrile seizures plus (GEFS+) caused by SCN1A gene mutations, to improve the understanding of SCN1A gene mutations, and to provide a basis for the clinical diagnosis and treatment of this disease. Methods Clinical data were collected from a child with GEFS+ and her family members who attended Department of Neurology in our hospital, and the third-generation gene sequencing technique was used to perform genetic analysis and validation of peripheral blood DNA of the child and her parents and aunt. Results A girl, aged 7 years, was admitted to the hospital due to "paroxysmal seizures lasting for more than 4 years", with loss of consciousness during seizures and more seizures during pyrexia. The seizures increased after the use of oxcarbazepine and were controlled after the admi-nistration of sodium valproate. Her father had a history of seizures, while her aunt had no clinical symptoms. Third-generation gene sequencing showed the same heterozygous mutation, c.4787G>A(p.R1596H), in the SCN1A gene in the child and her father and aunt. Conclusion The onset of GEFS+ is associated with SCN1A gene mutations, but there is heterogeneity in clinical manifestations among patients in one family, and sodium channel blockers as anti-epileptic drugs will aggravate epileptic seizures in GEFS+.

Background: Irritable bowel syndrome (IBS) is a functional bowel disorder (FBD). Objectives: To assess the therapeutic effects of paeoniflorin (PF) on IBS in rats. Method: Sixty male Sprague–Dawley rats were randomly divided into normal, model, positive drug, low-dose PF, medium-dose PF and high-dose PF groups (n = 10). After gavage for 2 consecutive weeks, the effect of PF on abdominal pain symptoms was assessed based on the abdominal withdrawal reflex (AWR) score, fecal water content and pathological changes in colon tissues. D-lactate, interleukin-1β (IL-1β), transforming growth factor-β (TGF-β) and tumor necrosis factor-α (TNF-α) were detected by enzyme-linked immunosorbent assay, and phosphorylated nuclear factor kappa B (p-NF-κB) p65 was detected by Western blotting. The abundance and diversity changes of intestinal flora were explored using 16S ribosomal RNA sequencing. Result: In PF groups, the mucosal morphology of colon tissues was intact, and the glands were arranged neatly and structured clearly, without obvious inflammatory cell infiltration. Compared with the model group, PF groups had significantly elevated pain threshold, and mRNA and protein levels of zonula occludens-1 (ZO-1) and occludin, decreased AWR score at 20 mmHg pressure, fecal water content, mRNA levels of IL-1β, TGF-β, and TNF-α, protein level of p-NF-κB p65 and level of serum D-lactate, and reduced levels of serum IL-1β, TGF-β, and TNF-α (p < 0.05, p < 0.01). PF groups had higher abundance of Lactobacillus, Akkermansia, Alistipes, and Bacteroides, but lower abundance of Desulfovibrio, Parasutterella, and Enterococcus than those of the model group. Conclusions: PF exerts therapeutic effects on IBS in rats probably by regulating the intestinal flora, and then up-regulating the expressions of ZO-1 and occludin in colon tissue while down-regulating the levels of IL-1β, TGF-β, TNF-α, D-lactate and p-NF-κB p65. [ABSTRACT FROM AUTHOR]

Background: Transplant arteriosclerosis is a major complication in long-term survivors of heart transplantation. Increased lymph flow from donor heart to host lymph nodes has been reported to play a role in transplant arteriosclerosis, but how lymphangiogenesis affects this process is unknown. Methods: Vascular allografts were transplanted among various combinations of mice, including wild-type, Lyve1 -CreERT2;R26-tdTomato, CAG -Cre-tdTomato, severe combined immune deficiency, Ccr2 KO, Foxn1 KO, and lghm / lghd KO mice. Whole-mount staining and 3-dimensional reconstruction identified lymphatic vessels within the grafted arteries. Lineage tracing strategies delineated the cellular origin of lymphatic endothelial cells. Adeno-associated viral vectors and a selective inhibitor were used to regulate lymphangiogenesis. Results: Lymphangiogenesis within allograft vessels began at the anastomotic sites and extended from preexisting lymphatic vessels in the host. Tertiary lymphatic organs were identified in transplanted arteries at the anastomotic site and lymphatic vessels expressing CCL21 (chemokine [C-C motif] ligand 21) were associated with these immune structures. Fibroblasts in the vascular allografts released VEGF-C (vascular endothelial growth factor C), which stimulated lymphangiogenesis into the grafts. Inhibition of VEGF-C signaling inhibited lymphangiogenesis, neointima formation, and adventitial fibrosis of vascular allografts. These studies identified VEGF-C released from fibroblasts as a signal stimulating lymphangiogenesis extending from the host into the vascular allografts. Conclusions: Formation of lymphatic vessels plays a key role in the immune response to vascular transplantation. The inhibition of lymphangiogenesis may be a novel approach to prevent transplant arteriosclerosis. [ABSTRACT FROM AUTHOR]

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&#95;133378.4):c.36812del; p.(Asp12271Valfs*10)], and a paternally inherited missense variant [(NM&#95;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.)

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 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).)

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).)

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.

Background: Myelodysplastic syndromes (MDS) are heterogeneous and clonal hematological disorders. The role and mechanism of necroptosis in MDS remain poorly understood., Methods: mRNA expression profiles and single-cell RNA-sequencing (scRNA-seq) data were sourced from the GEO database. ScRNA-seq data were processed using the "Seurat" package. After cell annotation, necroptosis-related scores (NRscores) for each cell were calculated using the "UCell" package. Differentially expressed genes (DEGs) and their associated biological functions in NRscore-related cell populations were identified. Additionally, DEGs and necroptosis-related genes (DE-NRGs) between MDS patients and healthy controls were identified. Consensus clustering was employed to classify MDS patients into distinct subclusters based on DE-NRGs. The biological functions and immune characteristics of these classifications were analyzed. Prognostic gene signatures were determined using LASSO and SVM-RFE analyses, and a nomogram was constructed based on the prognostic gene signature., Results: A total of 12 cell types were identified in MDS and healthy controls. NRscore was found to be elevated in monocytes and common lymphoid precursors (CLPs). Enrichment analysis revealed that monocytes and CLPs with high NRscore were associated with mitochondria-related and immune-related pathways. Eleven DEGs in monocytes and CLPs between MDS patients and healthy controls were identified. Additionally, 13 DE-NRGs were identified from 951 DEGs between MDS and healthy controls. MDS patients were classified into two distinct subclusters based on these 13 DE-NRGs, revealing several immune-related processes and signaling pathways. Differences in immune subpopulations between the two subclusters were observed. A necroptosis-related diagnostic gene signature (IRF9, PLA2G4A, MLKL, BAX, JAK2, and STAT3) was identified as predictive of MDS prevalence., Conclusion: Necroptosis plays a role in MDS progression by inducing inflammation. A novel necroptotic gene signature has been developed to distinguish and diagnose MDS at early stages of the disease., (© 2024. The Author(s).)

Over the past five decades, DNA restriction enzymes have revolutionized biotechnology. While these enzymes are widely used in DNA research and DNA engineering, the emerging field of RNA and mRNA therapeutics requires sequence-specific RNA endoribonucleases. Here, we describe EcoToxN1, a member of the type III toxin-antitoxin family of sequence-specific RNA endoribonucleases, and its use in RNA and mRNA analysis. This enzyme recognizes a specific pentamer in a single-stranded RNA and cleaves the RNA within this sequence. The enzyme is neither dependent on annealing of guide RNA or DNA oligos to the template nor does it require magnesium. Furthermore, it performs over a wide range of temperatures. With its unique functions and characteristics, EcoToxN1 can be classified as an RNA restriction enzyme. EcoToxN1 enables new workflows in RNA analysis and biomanufacturing, meeting the demand for faster, cheaper, and more robust analysis methods., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The pathogenesis of Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) remains unclear, though increasing evidence suggests inflammatory processes play key roles. In this study, single-cell RNA sequencing (scRNA-seq) of peripheral blood mononuclear cells (PBMCs) was used to decipher the immunometabolic profile in 4 ME/CFS patients and 4 heathy controls. We analyzed changes in the composition of major PBMC subpopulations and observed an increased frequency of total T cells and a significant reduction in NKs, monocytes, cDCs and pDCs. Further investigation revealed even more complex changes in the proportions of cell subpopulations within each subpopulation. Gene expression patterns revealed upregulated transcription factors related to immune regulation, as well as genes associated with viral infections and neurodegenerative diseases.CD4 + and CD8 + T cells in ME/CFS patients show different differentiation states and altered trajectories, indicating a possible suppression of differentiation. Memory B cells in ME/CFS patients are found early in the pseudotime, indicating a unique subtype specific to ME/CFS, with increased differentiation to plasma cells suggesting B cell overactivity. NK cells in ME/CFS patients exhibit reduced cytotoxicity and impaired responses, with reduced expression of perforin and CD107a upon stimulation. Pseudotime analysis showed abnormal development of adaptive immune cells and an enhanced cell-cell communication network converging on monocytes in particular. Our analysis also identified the estrogen-related receptor alpha (ESRRA)-APP-CD74 signaling pathway as a potential biomarker for ME/CFS in peripheral blood. In addition, data from the GSE214284 database confirmed higher ESRRA expression in the monocyte cell types of male ME/CFS patients. These results suggest a link between immune and neurological symptoms. The results support a disease model of immune dysfunction ranging from autoimmunity to immunodeficiency and point to amyloidotic neurodegenerative signaling pathways in the pathogenesis of ME/CFS. While the study provides important insights, limitations include the modest sample size and the evaluation of peripheral blood only. These findings highlight potential targets for diagnostic biomarkers and therapeutic interventions. Further research is needed to validate these biomarkers and explore their clinical applications in managing ME/CFS., (© 2024. The Author(s).)

Rosacea is a chronic inflammatory skin disorder, whose underlying cellular and molecular mechanisms remain obscure. Here, we generate a single-cell atlas of facial skin from female rosacea patients and healthy individuals. Among keratinocytes, a subpopulation characterized by IFNγ-mediated barrier function damage is found to be unique to rosacea lesions. Blocking IFNγ signaling alleviates rosacea-like phenotypes and skin barrier damage in mice. The papulopustular rosacea is featured by expansion of pro-inflammatory fibroblasts, Schwann, endothelial and macrophage/dendritic cells. The frequencies of type 1/17 and skin-resident memory T cells are increased, and vascular mural cells are characterized by activation of inflammatory pathways and impaired muscle contraction function in rosacea. Most importantly, fibroblasts are identified as the leading cell type producing pro-inflammatory and vasodilative signals in rosacea. Depletion of fibroblasts or knockdown of PTGDS, a gene specifically upregulated in fibroblasts, blocks rosacea development in mice. Our study provides a comprehensive understanding of the aberrant alterations of skin-resident cell populations and identifies fibroblasts as a key determinant in rosacea development., (© 2024. The Author(s).)

Angiogenesis, metastasis, and resistance to therapy are all facilitated by cancer-associated fibroblasts (CAFs). A CAF-based risk signature can be used to predict patients' prognoses for Lung adenocarcinoma (LUAD) based on CAF characteristics. The Gene Expression Omnibus (GEO) database was used to gather signal-cell RNA sequencing (scRNA-seq) data for this investigation. The GEO and TCGA databases were used to gather bulk RNA-seq and microarray data for LUAD. The scRNA-seq data were analyzed using the Seurat R program based on the CAF markers. Our goal was to use differential expression analysis to discover differentially expressed genes (DEGs) across normal and tumor samples in the TCGA dataset. Pearson correlation analysis was utilized to discover prognostic genes related with CAF, followed by univariate Cox regression analysis. Using Lasso regression, a risk signature based on CAF-related prognostic genes was created. A nomogram model was created based on the clinical and pathological aspects. 5 CAF clusters were identified in LUAD, 4 of which were associated with prognosis. From 2811 DEGs, 1002 genes were found to be significantly correlated with CAF clusters, which led to the creation of a risk signature with 8 genes. These 8 genes were primarily connected with 41 pathways, such as antigen paocessing and presentation, apoptosis, and cell cycle. Meanwhile, the risk signature was significantly associated with stromal and immune scores, as well as some immune cells. Multivariate analysis revealed that risk signature was an independent prognostic factor for LUAD, and its value in predicting immunotherapeutic outcomes was confirmed. A novel nomogram integrating the stage and CAF-based risk signature was constructed, which exhibited favorable predictability and reliability in the prognosis prediction of LUAD. CAF-based risk signatures can be effective in predicting the prognosis of LUAD, and they may provide new strategies for cancer treatments by interpreting the response of LUAD to immunotherapy., (© 2024. The Author(s).)

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

Background: Esophageal cancer (EC) is a major health issue, ranking seventh in incidence and sixth in mortality worldwide. Despite advancements in multidisciplinary treatment approaches, the 5-year survival rate for EC remains low at 21%. Challenges in EC treatment arise from late-stage diagnosis, high malignancy, and poor prognosis. Understanding the tumor microenvironment is critical, as it includes various cellular and extracellular components that influence tumor behavior and treatment response. Mast cells (MCs), as tissue-resident immune cells, play dual roles in tumor dynamics. High-throughput single-cell RNA sequencing offers a powerful tool for analyzing tumor heterogeneity and immune interactions, although its application in EC is limited., Methods: In this study, we investigated the immune microenvironment of EC using single-cell RNA sequencing and established a comprehensive immune profile. We also performed analysis of upstream transcription factors and downstream pathway enrichment to further comprehensively decipher MCs in EC. Besides, we performed knockdown experiments to explore the role of epidermal growth factor receptor ( EGFR ) signaling pathway in MCs-tumor cell interactions, highlighting its potential as a prognostic marker. Finally, we constructed a prognostic model for EC, which provided valuable suggestions for the diagnosis and prognosis of EC., Results: Our analysis identified 11 major cell types, of which MCs were particularly present in pericarcinoma tissues. Further grouping of the 5,001 MCs identified 8 distinct subtypes, including SRSF7 -highly expressed MCs, which showed strong tumor preference and potential tumor-promoting properties. Moreover, we identified the key signaling receptor EGFR and validated it by in vitro knockdown experiments, demonstrating its cancer-promoting effects. In addition, we established an independent prognostic indicator, SRSF7 + MCs risk score (SMRS), which showed a correlation between high SMRS group and poor prognosis., Conclusion: These findings illuminate the complex interactions within the tumor microenvironment of EC and suggest that targeting specific MCs subtypes, particularly via the EGFR signaling pathway, may present novel therapeutic strategies. This study establishes a comprehensive immune map of EC, offering insights for improved treatment approaches., 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 Zhang, Zhang, Xiahou, Zuo, Xue and Zhang.)

Background: Huntington's disease (HD) is a hereditary neurological disorder caused by mutations in HTT, leading to neuronal degeneration. Traditionally, HD is associated with the misfolding and aggregation of mutant huntingtin due to an extended polyglutamine domain encoded by an expanded CAG tract. However, recent research has also highlighted the role of global transcriptional dysregulation in HD pathology. However, understanding the intricate relationship between mRNA expression and HD at the cellular level remains challenging. Our study aimed to elucidate the underlying mechanisms of HD pathology using single-cell sequencing data., Results: We used single-cell RNA sequencing analysis to determine differential gene expression patterns between healthy and HD cells. HD cells were effectively modeled using a residual neural network (ResNet), which outperformed traditional and convolutional neural networks. Despite the efficacy of our approach, the F1 score for the test set was 96.53%. Using the SHapley Additive exPlanations (SHAP) algorithm, we identified genes influencing HD prediction and revealed their roles in HD pathobiology, such as in the regulation of cellular iron metabolism and mitochondrial function. SHAP analysis also revealed low-abundance genes that were overlooked by traditional differential expression analysis, emphasizing its effectiveness in identifying biologically relevant genes for distinguishing between healthy and HD cells. Overall, the integration of single-cell RNA sequencing data and deep learning models provides valuable insights into HD pathology., Conclusion: We developed the model capable of analyzing HD at single-cell transcriptomic level., (© 2024. The Author(s).)

Background: Batten disease is a group of rare inherited neurodegenerative diseases. Juvenile CLN3 disease is the most prevalent type, and the most common pathogenic variant shared by most patients is the "1-kb" deletion which removes two internal coding exons (7 and 8) in CLN3. Previously, we identified two transcripts in patient fibroblasts homozygous for the 1-kb deletion: the 'major' and 'minor' transcripts. To understand the full variety of disease transcripts and their role in disease pathogenesis, it is necessary to first investigate CLN3 transcription in "healthy" samples without juvenile CLN3 disease., Methods: We leveraged PacBio long-read RNA sequencing datasets from ENCODE to investigate the full range of CLN3 transcripts across various tissues and cell types in human control samples. Then we sought to validate their existence using data from different sources., Results: We found that a readthrough gene affects the quantification and annotation of CLN3. After taking this into account, we detected over 100 novel CLN3 transcripts, with no dominantly expressed CLN3 transcript. The most abundant transcript has median usage of 42.9%. Surprisingly, the known disease-associated 'major' transcripts are detected. Together, they have median usage of 1.5% across 22 samples. Furthermore, we identified 48 CLN3 ORFs, of which 26 are novel. The predominant ORF that encodes the canonical CLN3 protein isoform has median usage of 66.7%, meaning around one-third of CLN3 transcripts encode protein isoforms with different stretches of amino acids. The same ORFs could be found with alternative UTRs. Moreover, we were able to validate the translational potential of certain transcripts using public mass spectrometry data., Conclusion: Overall, these findings provide valuable insights into the complexity of CLN3 transcription, highlighting the importance of studying both canonical and non-canonical CLN3 protein isoforms as well as the regulatory role of UTRs to fully comprehend the regulation and function(s) of CLN3. This knowledge is essential for investigating the impact of the 1-kb deletion and rare pathogenic variants on CLN3 transcription and disease pathogenesis., (© 2024. The Author(s).)

Objectives: Monocytes are plastic cells that assume different polarization states that can either promote inflammation or tissue repair and inflammation resolution. Polarized monocytes are partially defined by their transcriptional profiles that are influenced by environmental stimuli. The airway monocyte response in pediatric acute respiratory distress syndrome (PARDS) is undefined. To identify differentially expressed genes and networks using a novel transcriptomic reporter assay with donor monocytes exposed to the airway fluid of intubated children with and at-risk for PARDS. To determine differences in gene expression at two time points using the donor monocyte assay exposed to airway fluid from intubated children with PARDS obtained 48-96 hours following initial tracheal aspirate sampling., Design: In vitro pilot study carried out using airway fluid supernatant., Setting: Academic 40-bed PICU., Participants: Fifty-seven children: 44 children with PARDS and 13 children at-risk for PARDS., Interventions: None., Measurements and Main Results: We performed bulk RNA sequencing using a transcriptomic reporter assay of monocytes exposed to airway fluid from intubated children to discover gene networks differentiating PARDS from at-risk for PARDS and those differentiating mild/moderate from severe PARDS. We also report differences in gene expression in children with PARDS 48-96 hours following initial tracheal aspirate sampling. We found that interleukin (IL)-10, IL-4, and IL-13, cytokine/chemokine signaling, and the senescence-associated secretory phenotype are upregulated in monocytes exposed to airway fluid from intubated children with PARDS compared with those at-risk for PARDS. Signaling by NOTCH, histone deacetylation/acetylation, DNA methylation, chromatin modifications (B-WICH complex), and RNA polymerase I transcription and its associated regulatory apparatus were upregulated in children with PARDS 48-96 hours following initial tracheal aspirate sampling., Conclusions: We identified gene networks important to the PARDS airway immune response using bulk RNA sequencing from a monocyte reporter assay that exposed monocytes to airway fluid from intubated children with and at-risk for PARDS. Mechanistic investigations are needed to validate our findings., Competing Interests: The authors have disclosed that they do not have any potential conflicts of interest., (Copyright © 2024 The Authors. Published by Wolters Kluwer Health, Inc. on behalf of the Society of Critical Care Medicine.)

Introduction: Rheumatoid arthritis (RA) is a chronic autoimmune disorder characterized by persistent synovitis, systemic inflammation, and autoantibody production. This study aims to explore the role of lactylation in plasma cells and its impact on RA pathogenesis., Methods: We utilized single-cell RNA sequencing (scRNA-seq) data and applied bioinformatics and machine learning techniques. A total of 10,163 cells were retained for analysis after quality control. Clustering analysis identified 13 cell clusters, with plasma cells displaying the highest lactylation scores. We performed pathway enrichment analysis to examine metabolic activity, such as oxidative phosphorylation and glycolysis, in highly lactylated plasma cells. Additionally, we employed 134 machine learning algorithms to identify seven core lactylation-promoting genes and constructed a diagnostic model with an average AUC of 0.918., Results: The RA lactylation score (RAlac&#95;score) was significantly elevated in RA patients and positively correlated with immune cell infiltration and immune checkpoint molecule expression. Differential expression analysis between two plasma cell clusters revealed distinct metabolic and immunological profiles, with cluster 2 demonstrating increased immune activity and extracellular matrix interactions. qRT-PCR validation confirmed that NDUFB3, NGLY1, and SLC25A4 are highly expressed in RA., Conclusion: This study highlights the critical role of lactylation in plasma cells for RA pathogenesis and identifies potential biomarkers and therapeutic targets, which may offer insights for future therapeutic strategies., 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 Fu, Wang, Lu, Shi and Yang.)

Background: Myocardial infarction (MI) induces complex transcriptional changes across diverse cardiac cell types. Single-cell RNA sequencing (scRNA-seq) provides an unparalleled ability to discern cellular diversity during infarction, yet the veracity of these discoveries necessitates confirmation. This investigation sought to elucidate MI mechanisms by integrating scRNA-seq and bulk RNA-seq data., Methods: Publicly available scRNA-seq (GSE136088) and bulk RNA-seq (GSE153485) data from mice MI models were analyzed. Cell types were annotated, and differential expression analysis conducted. Bulk RNA-seq underwent quality control, principal component analysis, and differential expression analysis., Results: In scRNA-seq data, the comparison between MI and sham groups unveiled a reduction in endothelial cell populations, but macrophages and monocytes increased. Within fibroblast subgroups, three distinct categories were discerned, with two exhibiting upregulation in MI. Notably, endothelial cells exhibited an elevated expression of genes associated with apoptosis and ferroptosis. In bulk RNA-seq analysis, distinct patterns emerged when comparing MI and sham groups. Specifically, six genes linked to endothelial ferroptosis exhibited heightened expression in MI group, thereby corroborating the scRNA-seq findings. Moreover, the examination of isolated cardiac macrophages from mice MI model revealed increased expression of Spp1, Col1a2, Col3a1, Ctsd, and Lgals3 compared to sham group, thus substantiating the dysregulation of macrophage apoptosis-related proteins following MI., Conclusion: MI altered the transcriptomic landscapes of cardiac cells with increased expression of apoptotic genes. Moreover, the upregulation of macrophage apoptosis marker was confirmed within MI models. The presence of endothelial cell depletion and ferroptosis in MI has been demonstrated., (© 2024. The Author(s).)

Purpose: This study aimed to define the heterogeneity, spatial localization, and functional roles of immune cells in the mouse cornea using single-cell RNA sequencing (scRNA-seq) and immunofluorescent staining., Methods: Enriched mouse corneal immune cells (C57BL/6 strain, age 16-20 weeks) underwent single-cell RNA sequencing library preparation, sequencing, and analysis with Seurat, Monocle 3, and CellChat packages in R. Pathway analysis used Qiagen Ingenuity Pathway Analysis software. Immunostaining confirmed cell distribution., Results: We identified 14 distinct immune cell clusters (56% myeloid and 44% lymphoid). Myeloid populations included resident macrophages, conventional dendritic cells (cDC2s), Langerhans cells, neutrophils, monocytes, and mast cells. Additionally, lymphocyte subsets (B, CD8, CD4, γδT, natural killer, natural killer T, and group 2 innate lymphoid cells) were found. We also found three new subtypes of resident macrophages in the cornea. Trajectory analysis suggested a differentiation pathway from monocytes to conventional dendritic cells, resident macrophages, and LCs. Intercellular communication network analysis using cord diagram identified amyloid precursor protein, chemokine (C-C motif) ligands (2, 3, 4, 6, 7, 9, and 12), Cxcl2, Mif, Tnf, Tgfb1, Igf1, and Il10 as prominent ligands involved in these interactions. Sexually dimorphic gene expression patterns were observed, with male myeloid cells expressing genes linked to immune regulation (Egr1, Foxp1, Mrc1, and Il1rn) and females showing higher expression of antigen presentation genes (Clic1, Psmb8, and Psmb9). Finally, immunostaining confirmed the spatial distribution of these cell populations within the cornea., Conclusions: This study unveils a diverse immune landscape in the mouse cornea, with evidence for cell differentiation and sex-based differences. Immunostaining validates the spatial distribution of these populations, furthering our knowledge of corneal immune function.

Purpose: Microglia-like cells derived from stem cells (iMG) provide a plentiful cell source for studying the functions of microglia in both normal and pathological conditions. Our goal is to establish a simplified and effective method for generating iMG in a precisely defined system. Additionally, we aim to achieve functional maturation of iMG through coculture with retinal organoids., Methods: In this study, iMG were produced under precisely defined conditions. They were subjected to LPS and poly IC stimulation. Additionally, we examined distinct phenotypic and functional variances between iMG and HMC3, a commonly used human microglia cell line. To investigate how the retinal cell interaction enhances microglial properties, iMG were cocultured with retinal organoids, producing CC-iMG. We performed RNA sequencing, electrophysiological analysis, and transmission electron microscope (TEM) to examine the maturation of CC-iMG compared to iMG., Results: Our results demonstrated that iMG performed immune-responsive profiles closely resembling those of primary human microglia. Compared to HMC3, iMG expressed a higher level of typical microglial markers and exhibited enhanced phagocytic activity. The transcriptomic analysis uncovered notable alterations in the ion channel profile of CC-iMG compared to iMG. Electrophysiological examination demonstrated a heightened intensity of inward- and outward-rectifying K+ currents in CC-iMG. Furthermore, CC-iMG displayed elevated numbers of lysosomes and mitochondria, coupled with increased phagocytic activity., Conclusions: These findings contribute to advancing our understanding of human microglial biology, specifically in characterizing and elucidating the functions of CC-iMG, thereby offering an in vitro microglial model for future scientific research and potential clinical applications in cell therapy.

Understanding of molecular mechanisms contributing to the pathophysiology of endometriosis, and upstream drivers of lesion formation, remains limited. Using a C57Bl/6 mouse model in which decidualized endometrial tissue is injected subcutaneously in the abdomen of recipient mice, we generated a comprehensive profile of gene expression in decidualized endometrial tissue (n=4), and in endometriosis-like lesions at Day 7 (n=4) and Day 14 (n=4) of formation. High-throughput mRNA sequencing allowed identification of genes and pathways involved in the initiation and progression of endometriosis-like lesions. We observed distinct patterns of gene expression with substantial differences between the lesions and the decidualized endometrium that remained stable across the two lesion timepoints, and showed similarity to transcriptional changes implicated in human endometriosis lesion formation. Pathway enrichment analysis revealed several immune and inflammatory response-associated canonical pathways, multiple potential upstream regulators, and involvement of genes not previously implicated in endometriosis pathogenesis, including IRF2BP2 and ZBTB10, suggesting novel roles in disease progression. Collectively, the provided data will be a useful resource to inform research on the molecular mechanisms contributing to endometriosis-like lesion development in this mouse model., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2024. Published by The Company of Biologists Ltd.)

Transmembrane protein 176B (TMEM176B), localized mainly on the endosomal membrane, has been reported as an immune regulatory factor in malignant diseases. However, the biological function of this molecule remains undetermined during respiratory viral infections. To investigate the functions and prognostic value of this gene, six gene sets were selected from the Gene Expression Omnibus database for research. First, the function of TMEM176B and its co-expressed genes were evaluated at different levels (cell, peripheral blood, lung tissue). Afterwards, a machine learning algorithm was utilized to analyze the relationship between TMEM176B and its interacting genes with prognosis. After importance evaluation and variable screening, a prognostic model was established. Finally, the reliability of the model was further verified through external data sets. In vitro experiments were conducted to validate the function of TMEM176B. TMEM176B and its co-expressed genes are involved in multiple processes such as inflammasome activation, myeloid immune cell development, and immune cell infiltration. Machine learning further screened 27 interacting gene modules including TMEM176B as prognostic models for severe respiratory viral infections, with the area under the ROC curve (AUCs) of 0.986 and 0.905 in derivation and external validation sets, respectively. We further confirmed that viral load as well as NLRP3 activation and cell death were significantly enhanced in TMEM176B -/- THP-1-differentiated macrophages via in vitro experiments. Our study revealed that TMEM176B is involved in a wide range of biological functions in respiratory viral infections and has potential prognostic value, which is expected to bring new insights into the clinical management of severe respiratory viral infection hosts., (© 2024 Wiley Periodicals LLC.)