Search

Purpose: The purpose of this study was to identify conjunctival transcriptome differences in patients with Acanthamoeba keratitis compared with keratitis with no known associated pathogen., Methods: The host conjunctival transcriptome of 9 patients with Acanthamoeba keratitis (AK) is compared with the host conjunctival transcriptome of 13 patients with pathogen-free keratitis. Culture and/or confocal confirmed Acanthamoeba in 8 of 9 participants with AK who underwent metagenomic RNA sequencing as the likely pathogen. Cultures were negative in all 13 cases where metagenomic RNA sequencing did not identify a pathogen., Results: Transcriptome analysis identified 36 genes differently expressed between patients with AK and patients with presumed sterile, or pathogen-free, keratitis. Gene enrichment analysis revealed that some of these genes participate in several biologic pathways important for cellular signaling, ion transport and homeostasis, glucose transport, and mitochondrial metabolism. Notable relatively differentially expressed genes with potential relevance to Acanthamoeba infection included CPS1 , SLC35B4 , STEAP2 , ATP2B2 , NMNAT3 , and AKAP12 ., Conclusions: This research suggests that the local transcriptome in Acanthamoeba keratitis may be sufficiently robust to be detected in the conjunctiva and that corneas infected with Acanthamoeba may be distinguished from the inflamed cornea where no pathogen was identified. Given the low sensitivity for corneal cultures, identification of differentially expressed genes may serve as a suggestive transcriptional signature allowing for a complementary diagnostic technique to identify this blinding parasite. Knowledge of differentially expressed genes may also direct investigation of disease pathophysiology and suggest novel pathways for therapeutic targets., Competing Interests: The authors have no funding or conflicts of interest to disclose., (Copyright © 2024 Wolters Kluwer Health, Inc. All rights reserved.)

Background: Eosinophils are elusive cells involved in allergic inflammation. Single-cell RNA-sequencing (scRNA-seq) is an emerging approach to deeply characterize cellular properties, heterogeneity, and functionality., Objectives: We sought to comprehensively characterize the transcriptome and biological functions of human eosinophils at a site of severe allergic inflammation in the esophagus (ie, eosinophilic esophagitis [EoE])., Methods: We employed a gravity-based scRNA-seq methodology to sequence blood eosinophils from patients with EoE and control individuals compared to a reanalyzed public scRNA-seq dataset of human esophageal eosinophils of EoE patients. We used flow cytometry, immunostaining, and a stimulation assay to verify mRNA findings., Results: In total, scRNA-seq was obtained from 586 eosinophils (188 from blood [n = 6 individuals] and 398 from esophagus [n = 6 individuals]). The esophageal eosinophils were composed of a population of activated eosinophils (enriched in 659 genes compared with peripheral blood-associated eosinophils) and a small population of eosinophils resembling peripheral blood eosinophils (enriched in 62 genes compared with esophageal eosinophils). Esophageal eosinophils expressed genes involved in sensing and responding to diverse stimuli, most notably IFN-γ, IL-10, histamine and leukotrienes, and succinate. Esophageal eosinophils were most distinguished from other esophageal populations by gene expression of the receptors CCR3, HRH4, SUCNR1, and VSTM1; transcription factors CEBPE, OLIG1, and OLIG2; protease PRSS33; and the hallmark eosinophil gene CLC. A web of bidirectional eosinophil interactions with other esophageal populations was derived. Comparing esophageal eosinophils and mast cells revealed that esophageal eosinophils expressed genes involved in DNAX-activation protein-12 (also known as TYROBP) interactions, IgG receptor-triggered events, immunoregulation, and IL-10 signaling., Conclusions: In EoE, esophageal eosinophils exist as 2 populations, a minority population resembling blood eosinophils and the other population characterized by high de novo transcription of diverse sensing receptors and inflammatory mediators readying them to potentially intersect with diverse cell types., Competing Interests: Disclosure statement Funding for this study was derived from National Institutes of Health grant R01 AI045898 and the Campaign Urging Research for Eosinophilic Diseases (CURED) awarded to M.E.R. This project was supported in part by National Institutes of Health grant P30 DK078392 (Gene Expression Core, Research Flow Cytometry Core) of the Digestive Diseases Research Core Center in Cincinnati. Disclosure of potential conflict of interest: M.E. Rothenberg is a consultant for Pulm One, Spoon Guru, ClostraBio, Serpin Pharm, Allakos, Celldex, Nexstone One, Santa Ana Bio, Enzen Pharmaceuticals, Bristol Myers Squibb, Astra Zeneca, Pfizer, GlaxoSmith Kline, Regeneron/Sanofi, Revolo Biotherapeutics, and Guidepoint and has an equity interest in the first 9 listed and royalties from reslizumab (Teva Pharmaceuticals), PEESSv2 (Mapi Research Trust), and UpToDate. M.E. Rothenberg is an inventor of patents owned by Cincinnati Children’s Hospital. The rest of the authors declare that they have no relevant conflicts of interest., (Copyright © 2024 American Academy of Allergy, Asthma & Immunology. All rights reserved.)

Wheat ( Triticum aestivum L.) is one of the most important crops worldwide. Powdery mildew caused by Blumeria graminis f. sp. tritici is a destructive disease threatening wheat yield and quality. The utilization of resistant genes and cultivars is considered the most economical, environmentally friendly, and effective method to control powdery mildew. Wheat breeding line Jingzi 102 was highly resistant to powdery mildew at both seedling and adult plant stages. Genetic analysis of F 1 , F 2 , and F 2:3 populations of "Jingzi 102 × Shixin 828" showed that the resistance of Jingzi 102 against powdery mildew isolate E09 at the seedling stage was controlled by a single dominant gene, temporarily designated PmJZ . Using bulked segregant RNA sequencing combined with molecular markers analysis, PmJZ was located on the long arm of chromosome 2B and flanked by markers BJK695 - 1 and CIT02g - 20 with the genetic distances of 1.2 and 0.5 centimorgan, respectively, corresponding to the bread wheat genome of Chinese Spring (International Wheat Genome Sequencing Consortium RefSeq v2.1) 703.8 to 707.6 Mb. PmJZ is most likely different from the documented Pm genes on chromosome 2BL based on their physical positions, molecular markers analysis, and resistance spectrum. Based on the gene annotation information, five genes related to disease resistance could be considered as the candidate genes of PmJZ . To accelerate the application of PmJZ , the flanking markers BJK695-1 and CIT02g-20 can serve for marker-assisted selection of PmJZ in wheat disease-resistance breeding., Competing Interests: The author(s) declare no conflict of interest.

Objective: This study employed a comprehensive single-cell analysis approach to explore the role of cell apoptosis-related genes in muscle aging., Methods: The single-cell RNA sequencing data from the GSE143704 dataset were used to identify distinct cell clusters and assess gene expression patterns related to apoptosis activation. The "limma" package was used to identify hub genes, after which we performed Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to identify relevant pathways. Additionally, Gene Set Enrichment Analysis(GSEA) and Gene Set Variation Analysis (GSVA) were used to uncover relevant biological pathways. The Receiver Operating Characteristic Curve (ROC) was used to evaluate the diagnostic value of the hub genes. Single-sample Gene Set Enrichment Analysis (ssGSEA) was used to analyze the immune cell infiltration levels., Results: Single-cell sequencing data from muscle aging patients allowed the identification of various cell types, including epithelial cells, adipocytes, and tissue-resident macrophages. By conducting a differential expression analysis that intersected active and nonactive apoptosis, as well as comparing elderly and young samples, a total of 22 hub genes were identified (p < 0.05). The 22 hub genes have discriminative ability as potential biomarkers for diagnosing muscle aging. The enrichment analysis indicated that these genes were closely associated with diverse pathways, including "response to UV-B" and "extracellular matrix organization" (p < 0.05). Furthermore, GSEA and GSVA indicated that multiple pathways emerged-for example, the "complement and coagulation cascades", "proteasome", "insulin signaling pathway", and "MAPK signaling pathway". Additionally, the analysis of immune cell infiltration revealed positive correlations between most of the hub genes and immune cells., Conclusion: Our study identified 22 apoptosis-related genes involved in muscle aging and indicated their potential diagnostic value. These findings offer a novel perspective on the pathogenesis of muscle aging and present potential targets for therapeutic interventions., Competing Interests: Declaration of competing interest The authors declared no potential conflict of interest with respect to the research, authorship, and publication of this article., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Canid alphaherpesvirus 1 (CHV-1) infection can cause spontaneous abortions in pregnant dams, and in young puppies, fatal systemic infections are common. MicroRNAs (miRNAs) affect viral infection by binding to messenger RNAs, and inhibiting expression of host and/or viral genes. We conducted deep sequencing of small RNAs in CHV-1-infected and mock-infected Madin-Darby Canine Kidney (MDCK) epithelial cells, and detected sequences corresponding to 282 cellular miRNAs. Of these, 18 were significantly upregulated at 12 h post-infection, most of which were encoded on the X chromosome. We next quantified the mature forms of several of the miRNAs using stem loop RT-qPCR. Our results revealed a discordance between the levels of small RNAs corresponding to canine miRNAs, and levels of the corresponding mature miRNAs, which suggests a block in miRNA biogenesis in infected cells. Nevertheless, we identified several mature miRNAs that exhibited a statistically significant increase upon infection. These included cfa-miR-8908b, a miRNA of unknown function, and cfa-miR-146a, homologs of which target innate immune pathways and are known to play a role in other viral infections. Interestingly, ontology analysis predicted that cfa-miR-8908b targets factors involved in the ubiquitin-like protein conjugation pathway and peroxisome biogenesis among other cellular functions. This is the first study to evaluate changes in miRNA levels upon CHV-1 infection. Based on our findings, we developed a model whereby CHV-1 infection results in changes in levels of a limited number of cellular miRNAs that target elements of the host immune response, which may provide clues regarding novel therapeutic targets., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Background: Chronic rhinosinusitis with nasal polyps (CRSwNP) is characterized by type 2 inflammation in the United States, but the actual roles that eosinophils play in CRSwNP remain largely unclear., Objective: To reveal the roles and heterogeneity of eosinophils in nasal polyp (NP) tissue, we performed single cell RNA sequencing (scRNA-Seq) analysis of NP tissue., Methods: Sinonasal tissues (NP and control sinus tissue) and patient matched peripheral blood (PB) samples were obtained from 5 control patients and 5 patients with CRSwNP. Eosinophils were enriched before processing for scRNA-Seq. The gene expression profiles in eosinophils were determined by microwell-based scRNA-Seq technology (BD Rhapsody platform). We predicted the overall function of NP eosinophils by Gene Ontology (geneontology.org) enrichment and pathway analyses and confirmed expression of selected genes by flow cytometry., Results: After filtering out contaminating cells, we detected 5,542 eosinophils from control PB, 3,883 eosinophils from CRSwNP PB, 101 eosinophils from control sinus tissues (not included in further analyses), and 9,727 eosinophils from NPs by scRNA-Seq. We found that 204 genes were downregulated and 354 genes upregulated in NP eosinophils compared to all PB eosinophils (>1.5-fold, P adj  < .05). Upregulated genes in NP eosinophils were associated with activation, cytokine-mediated signaling, growth factor activity, NF-κB signaling, and antiapoptotic molecules. NP eosinophils displayed 4 clusters revealing potential heterogeneity of eosinophils in NP tissue., Conclusions: Elevated eosinophils in NP tissue appear to exist in several subtypes that may play important pathogenic roles in CRSwNP, in part by controlling inflammation and hyperproliferation of other cells., Competing Interests: Disclosure statement Supported in part by Regeneron Pharmaceuticals, the National Institutes of Health (grants P01AI145818, R01AI137174, and U19AI136443), and a grant from the Ernest S. Bazley Foundation. Disclosure of potential conflict of interest: W. W. Stevens has served on advisory boards for GlaxoSmithKline, Regeneron, and Melinta Therapeutics. A. T. Petershas served on advisory boards for Sanofi-Genzyme/Regeneron, Optinose, AstraZeneca, Novartis, and GSK; and has received research support from Optinose and Sanofi/Regeneron. L. C. Grammer reports personal fees from Astellas Pharmaceuticals. K. C. Welch reports consultant fees from Baxter, OptiNose, and Acclarent. D. B. Conley reports consulting fees from Medtronic and Sanofi/Regeneron. R. P. Schleimer reports personal fees from Intersect ENT, Merck, GlaxoSmithKline, Sanofi, AstraZeneca/Medimmune, Genentech, Actobio Therapeutics, Lyra Therapeutics, Astellas Pharma, and Otsuka; and has royalty rights to Siglec-8 and Siglec-8 ligand–related patents licensed by Johns Hopkins to Allakos. R. C. Kern reports consulting fees from Lyra Therapeutics, Medtronic, GSK, Genentech and Sanofi/Regeneron. B. S. Bochner has received remuneration for serving on scientific advisory board of Allakos and owns stock in Allakos; has served as consultant for GSK, Third Harmonic Bio, Lupagen, Acelyrin, and Sanofi/Regeneron; receives publication-related royalty payments from Elsevier and UpToDate; is coinventor of existing Siglec-8–related patents and thus may be entitled to a share of royalties received by Johns Hopkins University during development and potential sales of such products; and is a cofounder of Allakos, which makes him subject to certain restrictions under university policy (terms of this arrangement are being managed by Johns Hopkins University and Northwestern University in accordance with their conflict-of-interest policies). B. K. Tan reports personal fees from Sanofi/Regeneron/Genzyme and GSK. A. Kato has served on advisory board for AstraZeneca; reports gift for his research from Lyra Therapeutics; and has received research grants from Regeneron and AstraZeneca. The rest of the authors declare that they have no relevant conflicts of interest., (Copyright © 2024 American Academy of Allergy, Asthma & Immunology. Published by Elsevier Inc. All rights reserved.)

Verrucous carcinoma (VC) is a rare subtype of squamous cell carcinoma (SCC) characterized by its histological presentation as a low-grade tumor with no potential for metastasis, setting it apart from invasive SCC. However, distinguishing VC from its benign counterpart, verrucous hyperplasia (VH), is challenging due to their clinical and morphological similarities. Despite the importance of accurate diagnosis for determining treatment strategies, diagnosis of VH and VC relied only on lesion recurrence after resection. To address this challenge, we generated RNA profiling data from tissue samples of VH and VC patients to identify novel diagnostic markers. We analyzed differentially expressed (DE) mRNA and long non-coding RNA (lncRNA) in tissue samples from VH and VC patients. Additionally, ChIP-X Enrichment Analysis 3 (ChEA3) was conducted to identify the top five transcription factors potentially regulating the expression of DE mRNAs in VH and VC. Our analysis of mRNA and lncRNA expression profiles in VH and VC provides insights into the underlying molecular characteristics of these diseases and offers potential new diagnostic markers. The identification of specific DE genes and lncRNAs may enable clinicians to more accurately differentiate between VH and VC, leading to better treatment choices., Competing Interests: Declaration of competing interest No potential conflict of interest relevant to this article was reported., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Age-related hearing loss (ARHL) or presbycusis is associated with irreversible progressive damage in the inner ear, where the sound is transduced into electrical signal; but the detailed mechanism remains unclear. Here, we sought to determine the potential molecular mechanism involved in the pathogeneses of ARHL with bioinformatics methods. A single-cell transcriptome sequencing study was performed on the cochlear samples from young and aged mice. Detection of identified cell type marker allowed us to screen 18 transcriptional clusters, including myeloid cells, epithelial cells, B cells, endothelial cells, fibroblasts, T cells, inner pillar cells, neurons, inner phalangeal cells, and red blood cells. Cell-cell communications were analyzed between young and aged cochlear tissue samples by using the latest integration algorithms Cellchat. A total of 56 differentially expressed genes were screened between the two groups. Functional enrichment analysis showed these genes were mainly involved in immune, oxidative stress, apoptosis, and metabolic processes. The expression levels of crucial genes in cochlear tissues were further verified by immunohistochemistry. Overall, this study provides new theoretical support for the development of clinical therapeutic drugs., (© 2024. The Author(s).)

Objective: Mechanisms underpinning clinical evolution to systemic lupus erythematosus (SLE) from preceding antinuclear antibodies (ANA) positivity are poorly understood. This study aimed to understand blood immune cell transcriptional signatures associated with subclinical ANA positivity, and progression or non-progression to SLE., Methods: Bulk RNA-sequencing of peripheral blood mononuclear cells isolated at baseline from 35 ANA positive (ANA+) subjects with non-diagnostic symptoms was analysed using differential gene expression, weighted gene co-expression network analysis, deconvolution of cell subsets and functional enrichment analyses. ANA+ subjects, including those progressing to classifiable SLE at 12 months (n=15) and those with stable subclinical ANA positivity (n=20), were compared with 15 healthy subjects and 18 patients with SLE., Results: ANA+ subjects demonstrated extensive transcriptomic dysregulation compared with healthy controls with reduced CD4+naïve T-cells and resting NK cells, but higher activated dendritic cells. B-cell lymphopenia was evident in SLE but not ANA+ subjects. Two-thirds of dysregulated genes were common to ANA+ progressors and non-progressors. ANA+ progressors showed elevated modular interferon signature in which constituent genes were inducible by both type I interferon (IFN-I) and type II interferon (IFN-II) in vitro. Baseline downregulation of mitochondrial oxidative phosphorylation complex I components significantly associated with progression to SLE but did not directly correlate with IFN modular activity. Non-progressors demonstrated more diverse cytokine profiles., Conclusions: ANA positivity, irrespective of clinical trajectory, is profoundly dysregulated and transcriptomically closer to SLE than to healthy immune function. Metabolic derangements and IFN-I activation occur early in the ANA+ preclinical phase and associated with diverging transcriptomic profiles which distinguish subsequent clinical evolution., Competing Interests: Competing interests: LMC has received consultancy fees from UCB and Alumis. MYMY has received consultancy fees from Aurinia Pharmaceuticals and UCB and speaker fees from Alumis, Roche and Novartis. EMV has received consultancy fees from Roche, GSK, AstraZeneca, Aurinia Pharmaceuticals, Lilly and Novartis. He has also received research grants paid to his employer from Roche, AstraZeneca and Sandoz. All other authors have declared no competing interests., (© Author(s) (or their employer(s)) 2024. No commercial re-use. See rights and permissions. Published by BMJ on behalf of EULAR.)

Accumulation of alpha-synuclein (α-syn) is a key pathological hallmark of synucleinopathies and has been shown to negatively impact neuronal function and activity. α-syn is an important factor contributing to astrocyte overactivation, though the effect of astrocyte overactivation on neurons remains unclear. Single-cell RNA sequencing data of mouse brain frontal cortex and midbrain from Hua-Syn (A53T) and wild type mice were utilized from the GEO database. Enrichment analysis, protein-protein interaction networks, and cell-cell interaction networks all indicated enhanced communication between astrocytes and neurons, along with the involvement of TNF and inflammation-related signaling pathways. In vitro experiments were performed to further explore the mechanism of neurotoxicity in astrocyte-neuron crosstalk. Astrocytes were treated by α-syn, neuronal TNFR1 receptors were antagonized by R-7050, and the cells were co-cultured after 24 h treatment. ELISA results revealed that cytokines such as TNF-α and IL-6 were significantly upregulated in astrocytes following the endocytosis of α-syn. Immunofluorescence (IF) showed neuronal dendritic reduction, axon elongation and increased co-localisation of TNFR1 receptor expression. Western blot showed up-regulation of PKR, P-eIF2α and ATF4 protein expression. Conversely, after antagonizing neuronal TNFR1 receptors with the R-7050 chemical inhibitor, neuronal synaptic structure was significantly restored and the expression of PKR, P-eIF2α and ATF4 was down-regulated. In summary, TNF-α acts as a signaling molecule mediating the up-regulated astrocyte-neuron crosstalk, providing new insights into the pathogenesis of α-syn-related neurological disorders., 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.)

Localized scleroderma (LS) is an autoimmune disease characterized by inflammation and fibrosis, leading to severe cutaneous manifestations such as skin hardening, tightness, discoloration, and other textural changes that may result in disability. While LS shares similar histopathologic features and immune-fibroblast interactions with systemic sclerosis (SSc), its molecular mechanisms remain understudied. Endothelial cells (EC) are known to play a crucial role in SSc but have not been investigated in LS. Single-cell RNA sequencing (scRNA-seq) now allows for detailed examination of this cell type in the primary organ of interest for scleroderma, the skin. In this study, we analyzed skin-isolated cells from 27 LS patients (pediatric and adult) and 17 healthy controls using scRNA-seq. Given the known role of EC damage as an initial event in SSc and the histologic and clinical skin similarities to LS, we focused primarily on endothelial cells. Our analysis identified eight endothelial subclusters within the dataset, encompassing both disease and healthy samples. Interaction analysis revealed that signaling from diseased endothelial cells was predicted to promote fibrosis through SELE interaction with FGFBP1 and other target genes. We also observed high levels of JAG in arterial endothelial cells and NOTCH in capillary endothelial cells, indicating the activation of a signaling pathway potentially responsible for epidermal abnormalities and contributing to LS pathogenesis. In summary, our scRNA-seq analysis identified potential disease-propagating endothelial cell clusters with upregulated pathways in LS skin, highlighting their importance in disease progression.

Barkground: Circular RNAs (circRNAs) play important regulatory roles in a variety of biological processes in mammals. Multiple birth-traits in goats are affected by several factors, but the expression and function of circRNAs in follicular development of goats are not clear. In this study, we aimed to investigate the possible regulatory mechanisms of circRNA and collected five groups of large follicles (Follicle diameter > 6 mm) and small follicles (1 mm < Follicle diameter < 3 mm) from Leizhou goats in estrus for RNA sequencing., Results: RNA sequencing showed that 152 circRNAs were differentially expressed in small and large follicles. Among them, 101 circRNAs were up-regulated in large follicles and 51 circRNAs were up-regulated in small follicles. GO and KEGG enrichment analyses showed that parental genes of the differential circRNAs were significantly enriched in important pathways, such as ovarian steroidogenesis, GnRH signaling pathway, animal autophagy and oxytocin signalling pathway. BioSignal analysis revealed that 152 differentially expressed circRNAs could target 91 differential miRNAs including miR-101 family (chi-miR-101-3p, chi-miR-101-5p), miR-202 family (chi-miR-202-5p, chi-miR-202-3p),60 circRNAs with translation potential. Based on the predicted sequencing results, the ceRNA networks chicirc&#95;008762/chi-miR-338-3p/ARHGAP18 and chicirc&#95;040444/chi-miR-338-3p/STAR were constructed in this study. Importantly, the new gene circCFAP20DC was first discovered in goats. The EDU assay and flow cytometry results indicated that circCFAP20DC enhanced the proliferation of follicular granulosa cells(GCs). Real-time quantitative PCR and western blotting assays showed that circCFAP20DC activated the Retinoblastoma(RB) pathway and promoted the progression of granulosa cells from G1 to S phase., Conclusion: Differential circRNAs in goat size follicles may have important biological functions for follicular development. The novel gene circCFAP20DC activates the RB pathway, promoting the progression of GCs from G1 to S phase. This, in turn, enhances the proliferation of follicular GCs in goats., (© 2024. The Author(s).)

Normal veins could develop to varicose vein (VV) by some risk factors, and might further progress to shallow vein thrombosis (SVT). However, the molecular mechanism of key genes associated with the progression and regression of VV are still not thorough enough. In this study, the healthy control (HC), VV, and SVT vascular samples were collected for transcriptome sequencing. The differentially expressed genes (DEGs) were screened by "DESeq2", including DEGs1 (HC vs. VV), DEGs2 (HC vs. SVT) and DEGs3 (VV vs. SVT). And their functional enrichment analyses were conducted by "ClusterProfiler". The receiver operating characteristic (ROC) curve was used to obtain the key genes (KGs) of the pathogenesis of VV and SVT. The qRT-PCR assay was performed to validate the expressions of KGs. Immune cell infiltration analyses were conducted based on ssGSEA method. The competitive endogenous RNAs (ceRNAs) regulatory network was constructed. The target drugs of KGs were predicted using DrugBank database. The biofunctions of DACT3 were further investigated through a series of experiments in vitro. All of these DEGs were associated with inflammation and immunity related functions. Immune cell infiltration was significantly different between VV and SVT. Six key genes including PLP2, DACT3, LRRC25, PILRA, MSX1 and APOD that were associated with the progression and regression of VV were screened. The expression of LRRC25 and PILRA was significantly negatively associated with central memory T cell, and significantly positively associated with B cell. Besides, XIST was the critical regulator of multiple KGs. Cimetidine was potential drug for VV and SVT therapy. Overexpression of DACT3 significantly inhibited the proliferation and migration of vascular smooth muscle cells (VSMCs), and affected their cell cycle and phenotypic transition. This study identified six key genes associated with the progression and regression of VV. Among them, DACT3 was proved to hinder VV progression. These findings may help to deepen understanding its underlying mechanisms., (© 2024. The Author(s).)

PANoptosis is engaged in the program of immune response and carcinogenicity. Nonetheless, the actual impacts of PANoptosis on clinical management and oncology immunity in hepatocellular carcinoma (HCC) are not fully grasped. RNA-seq-derived computations were conducted to sort out the molecular subtypes and elucidate the disparities based on PANoptosis molecules. Single-cell sequencing (scRNA-seq) tools including Cytotrace and Addmodulescore were extracted to characterize diversification potency and quantify the PANoptosis motion. Transcriptional factors were inferred by the pySCENIC package and Cellchat program scrutinized the intercellular exchange across cell compartments. The PANoptosis score system originated by incorporating 10 machine learning algorithms and 101 compositions to project clinical results and deteriorate tendencies. Circulatory PANoptosis-associated protein HSP90AA1 was determined by enzyme-linked immunosorbent assay (ELISA). HCC individuals could be categorized into low- and high-PANoptosis groups with diverse biogenic and pharmacotherapy heterogeneity. Individuals in the elevated PANoptosis subtype were characterized as "hot tumor" conveying the increased presence of immunogenicity while reiterating an explicit negative connection with tumor stemness. Compared to immune and stromal cells, cancerous cells showcased decreased PANoptosis and heightened PANoptosis malignant cell subgroups might be tied to a substantial level of genomic expression of SREBF2, JUND, GATAD1, ZBTB20, SMAD5 and implied a more aggressive potential. The PANoptosis index, derived from machine learning, has been established to provide succinct frameworks for predicting outcomes and clarified the noteworthy utility of conventional regimens, as the differentiated power of HCC occurred together with vascular invasion and hepatocellular adenoma (HCA). The experiment confirmed that the circulating HSP90AA1 was aberrantly augmented in HCC patients, thus demonstrating its potential as a discriminatory biomarker. We systematically deciphered the molecular and immune ecosystem traits of PANoptosis in bulk and scRNA-seq degrees, which may deliver advantageous insights for customized treatment, awareness of the pathological process and prognosis scrutiny for HCC patients., (© 2024. The Author(s).)

Background: Several observational studies have indicated an association between interstitial cystitis and the composition of the gut microbiota; however, the causality and underlying mechanisms remain unclear. Understanding the link between gut microbiota and interstitial cystitis could inform strategies for prevention and treatment., Methods: A two-sample Mendelian randomization analysis was conducted using published genome-wide association study summary statistics. We employed inverse variance weighted, weighted mode, MR-Egger, weighted median, simple mode, and cML-MA methods to investigate the causal relationship between gut microbiota and interstitial cystitis. Sensitivity analysis was performed to validate the results. Relevant gut microbiota was examined through reverse MR. Single nucleotide polymorphisms were annotated using FUMA to identify genes associated with these genetic variants, thereby revealing potential host gene-microbiota associations in interstitial cystitis patients., Results: Eight bacterial taxa were identified in our analysis as associated with interstitial cystitis. Among these, Butyricimonas , Coprococcus , Lactobacillales , Lentisphaerae , and Bilophila wadsworthia were positively correlated with interstitial cystitis risk, while taxa such as Desulfovibrio piger , Oscillibacter unclassified and Ruminococcus lactaris exhibited protective effects against interstitial cystitis. The robustness of these associations was confirmed through sensitivity analyses. Reverse MR analysis did not reveal evidence of reverse causality. Single nucleotide polymorphisms were annotated using FUMA and subjected to biological analysis. Seven hub genes (SPTBN1, PSME4, CHAC2, ERLEC1, ASB3, STAT5A, and STAT3) were identified as differentially expressed between interstitial cystitis patients and healthy individuals, representing potential therapeutic targets., Conclusion: Our two-sample Mendelian randomization study established a causal relationship between gut microbiota and interstitial cystitis. Furthermore, our identification of a host gene-microbiota association offers a new avenue for investigating the potential pathogenesis of interstitial cystitis and suggests avenues for the development of personalized treatment 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, Zhao, Zhou, Lv, Jin, Zhou, Liu and Feng.)

Purpose: Proliferation of stromal connective tissue is a hallmark of pancreatic adenocarcinoma (PAAD). The engagement of activated cancer-associated fibroblasts (CAFs) contributes to the progression of PAAD through their involvement in tumor fibrogenesis. However, the prognostic significance of CAF-based risk signature in PAAD has not been explored., Methods: The single-cell RNA sequencing (scRNA-seq) data sourced from GSE155698 within the Gene Expression Omnibus (GEO) database was supplemented by bulk RNA sequencing data from The Cancer Genome Atlas (TCGA) and microarray data retrieved from the GEO database. The scRNA-seq data underwent processing via the Seurat package to identify distinct CAF clusters utilizing specific CAF markers. Differential gene expression analysis between normal and tumor samples was conducted within the TCGA-PAAD cohort. Univariate Cox regression analysis pinpointed genes associated with CAF clusters, identifying prognostic CAF-related genes. These genes were utilized in LASSO regression to craft a predictive risk signature. Subsequently, integrating clinicopathological traits and the risk signature, a nomogram model was constructed., Results: Our scRNA-seq analysis unveiled four distinct CAF clusters in PAAD, with two linked to PAAD prognosis. Among 207 identified DEGs, 148 exhibited significant correlation with these CAF clusters, forming the basis of a seven-gene risk signature. This signature emerged as an independent predictor in multivariate analysis for PAAD and demonstrated predictive efficacy in immunotherapeutic outcomes. Additionally, a novel nomogram, integrating age and the CAF-based risk signature, exhibited robust predictability and reliability in prognosticating PAAD. Moreover, the risk signature displayed substantial correlations with stromal and immune scores, as well as specific immune cell types., Conclusions: The prognosis of PAAD can be accurately predicted using the CAF-based risk signature, and a thorough analysis of the PAAD CAF signature may aid in deciphering the patient's immunotherapy response and presenting fresh cancer treatment options.

Background: Cervical cancer (CC) is a neoplasia with a high heterogeneity. We aimed to explore the characteristics of tumor microenvironment (TME) for CC treatment., Methods: HPV positive (+) and negative (-) samples from cervical cancer (CC) patients were sourced from the Gene Expression Omnibus (GEO) database. The single-cell RNA sequencing (scRNA-seq) data were processed and annotated for cell types utilizing the Seurat package. Following this, the expression levels and biological roles of the marker genes were analyzed applying real-time PCR (RT-PCR) and transwell assays. Furthermore, the enrichment of genes with significantly differential expressions and copy number variations was assessed by the ClusterProlifer and inferCNV software packages., Results: Seven main cell clusters were classified based on a total of 12,431 cells. The HPV- CC samples exhibited a higher immune cell infiltration level, while epithelial cells and myofibroblasts had higher proportion in the HPV+ CC samples with extensive heterogeneity. Immune pathways including antigen treatment and presentation, immunoglobulin production and T cell mediated immunity were significantly activated in the HPV- CC group with lower cell cycle and proliferation activity. However, the anti-tumor immunity of these cells was inhibited in HPV+ CC group with higher cell proliferation activity. Moreover, the amplification and loss of CNVs also supported that these cells in HPV- CC samples were prone to anti-tumor activation. Further cell validation results showed that except GZMA, the levels of APOC1, CEACAM6, FOXP3, SFRP4 and TFF3 were all higher in CC cells Hela, and that silencing TFF3 could inhibit the migration and invasion of CC cells in-vitro ., Conclusion: This study highlighted the critical role of HPV infection in CC progression, providing a novel molecular basis for optimizing the current preventive screening and personalized treatment for the cancer., Competing Interests: The authors declare that they have no competing interests., (© 2024 Zhang et al.)

Background: Severe acute pancreatitis (SAP) is characterized by inflammation, with inflammatory immune cells playing a pivotal role in disease progression. This study aims to understand variations in specific immune cell subtypes in SAP, uncover their mechanisms of action, and identify potential biological markers for predicting Acute Pancreatitis (AP) severity., Methods: We collected peripheral blood from 7 untreated SAP patients and employed single-cell RNA sequencing for the first time to construct a transcriptome atlas of peripheral blood mononuclear cells (PBMCs) in SAP. Integrating SAP transcriptomic data with 6 healthy controls from the GEO database facilitated the analysis of immune cell roles in SAP. We obtained comprehensive transcriptomic datasets from AP samples in the GEO database and identified potential biomarkers associated with AP severity using the "Scissor" tool in single-cell transcriptomic data., Results: This study presents the inaugural construction of a peripheral blood single-cell atlas for SAP patients, identifying 20 cell subtypes. Notably, there was a significant decrease in effector T cell subsets and a noteworthy increase in monocytes compared to healthy controls. Moreover, we identified a novel monocyte subpopulation expressing high levels of PPBP and PF4 which was significantly elevated in SAP. The proportion of monocyte subpopulations with high CCL3 expression was also markedly increased compared to healthy controls, as verified by flow cytometry. Additionally, cell communication analysis revealed insights into immune and inflammation-related signaling pathways in SAP patient monocytes. Finally, our findings suggest that the subpopulation with high CCL3 expression, along with upregulated pro-inflammatory genes such as S100A12 , IL1B , and CCL3 , holds promise as biomarkers for predicting AP severity., Conclusion: This study reveals monocytes' crucial role in SAP initiation and progression, characterized by distinct pro-inflammatory features intricately linked to AP severity. A monocyte subpopulation with elevated PPBP and CCL3 levels emerges as a potential biomarker and therapeutic target., 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, Wang, Wu, Tao, Li, Zheng, Xu, Du, Zhao, Liang, Xu and Wang.)

Aldehyde dehydrogenase 1, family member A2, is a retinoic acid-synthesizing enzyme encoded by Aldh1a2 in mice and ALDH1A2 in humans. This enzyme is indispensable for kidney development, but its role in kidney physiology and pathophysiology remains to be fully defined. In this review, we mined single-cell and single-nucleus RNA sequencing databases of mouse and human kidneys and found that glomerular parietal epithelial cells (PECs) express a full set of genes encoding proteins needed for cellular vitamin A uptake, intracellular transport, and metabolism into retinoic acid. In particular, Aldh1a2/ALDH1A2 mRNAs are selectively enriched in mouse and human PECs. Aldh1a2 expression in PECs is greatly increased in a mouse model of anti-glomerular basement membrane glomerulonephritis and moderately induced in a mouse model of ischemia-reperfusion acute kidney injury. Aldh1a2 expression in PECs is substantially repressed in a chronic kidney disease mouse model combining diabetes, hypertension, and partial nephrectomy and is moderately repressed in mouse models of focal segmental glomerulosclerosis and diabetic nephropathy. Single-nucleus RNA sequencing data show that ALDH1A2 mRNA expression in PECs is diminished in patients with chronic kidney disease associated with diabetes, hypertension and polycystic kidney disease. In addition to data mining, we also performed Spearman's rank correlation coefficient analyses and identified gene transcripts correlated with Aldh1a2/ALDH1A2 transcripts in mouse PECs and PEC subtypes, and in human PECs of healthy subjects and patients with AKI or CKD. Furthermore, we conducted Gene Ontology pathway analyses and identified the biological pathways enriched among these Aldh1a2/ALDH1A2 -correlated genes. Our data mining and analyses led us to hypothesize that ALDH1A2 - mediated retinoic acid synthesis in PECs plays a yet-undefined role in the kidney and that its dysregulation mediates injury. Conditional, PEC-selective Aldh1a2 knockout, RNA silencing and transgenic mouse models will be useful tools to test this hypothesis. Clinical studies on genetics, epigenetics, expression and functions of ALDH1A2 and other genes needed for retinoic acid biosynthesis and signaling are also warranted., 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 Liu, Fermin, Xu, Kopp and Xu.)

Fibrosis is characterized by the aberrant deposition of extracellular matrix (ECM) due to dysregulated tissue repair responses, imposing a significant global burden on fibrosis-related diseases. Although alpha-smooth muscle actin (α-SMA/ ACTA2 )-expressing myofibroblasts are considered as key player in fibrogenesis, the origin of ECM-producing cells remains controversial. To address this issue, we integrated and analyzed large-scale single-cell transcriptomic datasets from patients with distinct fibrotic diseases involving the heart, lung, liver, or kidney. Unexpectedly, not all ACTA2- expressing cells were ECM-producing cells identified by expressing collagen genes; instead, the majority of ECM-producing cells were myofibroblasts and fibroblasts derived from circulating bone marrow precursor, and to a lesser extent from local pericytes and vascular smooth cells in all fibrotic diseases. This was confirmed in sex-mismatched kidney transplants by the discovery that ECM-producing cells originated from recipient, not donor, bone marrow-derived progenitor cells (BMPCs). Moreover, these BMPCs-derived ECM-producing cells exhibited a proinflammatory phenotype. Thus, bone marrow-derived proinflammatory and profibrotic fibroblasts/myofibroblasts with stem cell properties serve as a major source of ECM-producing cells and may play a driving role in tissue fibrosis across a wide range of human fibrotic diseases. Targeting these ECM-producing cells may provide a novel therapy for diseases with fibrosis., Competing Interests: Competing Interests: The authors have declared that no competing interest exists., (© The author(s).)

Rationale: The chronic lung disease bronchopulmonary dysplasia (BPD) is the most severe complication of extreme prematurity. BPD results in impaired lung alveolar and vascular development and long-term respiratory morbidity, for which only supportive therapies exist. Umbilical cord-derived mesenchymal stromal cells (UC-MSCs) improve lung structure and function in experimental BPD. Results of clinical trials with MSCs for many disorders do not yet match the promising preclinical studies. A lack of specific criteria to define functionally distinct MSCs persists. Objectives: To determine and correlate single-cell UC-MSC transcriptomic profiles with therapeutic potential. Methods: UC-MSCs from five term donors and human neonatal dermal fibroblasts (HNDFs; control cells of mesenchymal origin) transcriptomes were investigated using single-cell RNA sequencing (scRNA-seq) analysis. The lung-protective effect of UC-MSCs with a distinct transcriptome and control HNDFs was tested in vivo in hyperoxia-induced neonatal lung injury in rats. Measurements and Main Results: UC-MSCs showed limited transcriptomic heterogeneity but were different from HNDFs. Gene Ontology enrichment analysis revealed distinct (progenitor-like and fibroblast-like) UC-MSC subpopulations. Only treatment with progenitor-like UC-MSCs improved lung function and structure and attenuated pulmonary hypertension in hyperoxia-exposed rat pups. Moreover, scRNA-seq identified major histocompatibility complex class I as a molecular marker of nontherapeutic cells and associated with decreased lung retention. Conclusions: UC-MSCs with a progenitor-like transcriptome, but not with a fibroblast-like transcriptome, provide lung protection in experimental BPD. High expression of major histocompatibility complex class I is associated with reduced therapeutic benefit. scRNA-seq may be useful to identify subsets of MSCs with superior repair capacity for clinical application.

Background: Milk is essential for mammalian nutrition because it provides vital nutrients for growth and development. Milk composition, which is influenced by genetic and environmental factors, supports lactation, a complex process crucial for milk production and quality. Recent research has focused on noncoding RNAs, particularly microRNAs (miRNAs), which are present in body fluids and regulate gene expression post-transcriptionally. This study comprehensively characterizes miRNAs in milk of four livestock species, namely Bubalus bubalis, Capra hircus, Equus asinus, and Ovis aries and identifies potential target genes., Results: High-throughput sequencing of milk RNA resulted in distinct read counts across species: B. bubalis (8,790,441 reads), C. hircus (12,976,275 reads), E. asinus (9,385,067 reads), and O. aries (7,295,297 reads). E. asinus had the highest RNA mapping rate (94.6%) and O. aries the lowest (84.8%). A substantially greater proportion of miRNAs over other small RNAs was observed for the donkey milk sample (7.74%) compared to buffalo (0.87%), goat (1.57%), and sheep (1.12%). Shared miRNAs, which included miR-200a, miR-200b, miR-200c, and miR-23a among others, showed varying expression levels across species, confirmed by qPCR analysis. Functional annotation of predicted miRNA target genes highlighted diverse roles, with an enrichment in functions linked to metabolism and immunity. Pathway analysis identified immune response pathways as significant, with several miRNAs targeting specific genes across species, suggesting their regulatory function in milk., Conclusions: Both conserved and species-specific miRNAs were detected in milk of the investigated species. The identified target genes of these miRNAs have important roles in neonatal development, adaptation, growth, and immune response. Furthermore, they influence milk and meat production traits in livestock., (© 2024. The Author(s).)

Background: Although the patient survival rate for many malignancies has been improved with immune checkpoint inhibitors (ICIs), some patients experience various immune-related adverse events (irAEs). IrAEs impact several organ systems, including the kidney. With anti-programmed cell death protein 1 (PD-1) therapy (pembrolizumab), kidney-related adverse events occur relatively rarely compared with other irAEs. However, the occurrence of AKI usually leads to anti-PD-1 therapy interruption or discontinuation. Therefore, there is an urgent need to clarify the mechanisms of renal irAEs (R-irAEs) to facilitate early management. This study aimed to analyse the characteristics of peripheral blood mononuclear cells (PBMCs) in R-irAEs., Methods: PBMCs were collected from three patients who developed R-irAEs after anti-PD-1 therapy and three patients who did not. The PBMCs were subjected to scRNA-seq to identify cell clusters and differentially expressed genes (DEGs). Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analyses were performed to investigate the most active biological processes in immune cells., Results: Fifteen cell clusters were identified across the two groups. FOS, RPS26, and JUN were the top three upregulated genes in CD4 + T cells. The DEGs in CD4 + T cells were enriched in Th17 differentiation, Th1 and Th2 cell differentiation, NF-kappa B, Nod-like receptor, TNF, IL-17, apoptosis, and NK cell-mediated cytotoxicity signaling pathways. RPS26, TRBV25-1, and JUN were the top three upregulated genes in CD8 + T cells. The DEGs in CD8 + T cells were enriched in Th17 cell differentiation, antigen processing and presentation, natural killer cell-mediated cytotoxicity, the intestinal immune network for IgA production, the T-cell receptor signalling pathway, Th1 and Th2 cell differentiation, the phagosome, and cell adhesion molecules., Conclusions: In conclusion, R-irAEs are associated with immune cell dysfunction. DEGs and their enriched pathways identified in CD4 + T cells and CD8 + T cells play important roles in the development of renal irAEs related to anti-PD-1 therapy. These findings offer fresh perspectives on the pathogenesis of renal damage caused by anti-PD-1 therapy., (© 2024. The Author(s).)

Background: Hepatocellular carcinoma (HCC) is a prevalent and heterogeneous tumor with limited treatment options and unfavorable prognosis. The crucial role of a disintegrin and metalloprotease (ADAM) gene family in the tumor microenvironment of HCC remains unclear., Methods: This study employed a novel multi-omics integration strategy to investigate the potential roles of ADAM family signals in HCC. A series of single-cell and spatial omics algorithms were utilized to uncover the molecular characteristics of ADAM family genes within HCC. The GSVA package was utilized to compute the scores for ADAM family signals, subsequently stratified into three categories: high, medium, and low ADAM signal levels through unsupervised clustering. Furthermore, we developed and rigorously validated an innovative and robust clinical prognosis assessment model by employing 99 mainstream machine learning algorithms in conjunction with co-expression feature spectra of ADAM family genes. To validate our findings, we conducted PCR and IHC experiments to confirm differential expression patterns within the ADAM family genes., Results: Gene signals from the ADAM family were notably abundant in endothelial cells, liver cells, and monocyte macrophages. Single-cell sequencing and spatial transcriptomics analyses have both revealed the molecular heterogeneity of the ADAM gene family, further emphasizing its significant impact on the development and progression of HCC. In HCC tissues, the expression levels of ADAM9, ADAM10, ADAM15, and ADAM17 were markedly elevated. Elevated ADAM family signal scores were linked to adverse clinical outcomes and disruptions in the immune microenvironment and metabolic reprogramming. An ADAM prognosis signal, developed through the utilization of 99 machine learning algorithms, could accurately forecast the survival duration of HCC, achieving an AUC value of approximately 0.9., Conclusions: This study represented the inaugural report on the deleterious impact and prognostic significance of ADAM family signals within the tumor microenvironment of HCC., 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 Chen, Yuan, Guan, Cui, Fu, Wei and Liu.)

Thyroid hormone receptor alpha (THRα) is a nuclear hormone receptor that binds triiodothyronine (T3) and acts as an important transcription factor in development, metabolism, and reproduction. In mammals, THRα has two major splicing isoforms, THRα1 and THRα2. The better-characterized isoform, THRα1, is a transcriptional stimulator of genes involved in cell metabolism and growth. The less-well-characterized isoform, THRα2, lacks the ligand-binding domain (LBD) and is thought to act as an inhibitor of THRα1 activity. The ratio of THRα1 to THRα2 splicing isoforms is therefore critical for transcriptional regulation in different tissues and during development. However, the expression patterns of both isoforms have not been studied in healthy human tissues or in the developing brain. Given the lack of commercially available isoform-specific antibodies, we addressed this question by analyzing four bulk RNA-sequencing datasets and two scRNA-sequencing datasets to determine the RNA expression levels of human THRA1 and THRA2 transcripts in healthy adult tissues and in the developing brain. We demonstrate how 10X Chromium scRNA-seq datasets can be used to perform splicing-sensitive analyses of isoforms that differ at the 3'-end. In all datasets, we found a strong predominance of THRA2 transcripts at all examined stages of human brain development and in the central nervous system of healthy human adults.

Low back pain (LBP) is largely attributed to intervertebral disc degeneration (IVDD), of which the endplate changes are an important component. However, the alterations in cell fate and properties within the endplates during degeneration remain unknown. Here, we firstly performed the single-cell RNA-sequencing analysis (scRNA-seq) of the cells focusing on degenerative human endplates. By unsupervised clustering of the 8,534 single-cell based on the gene expression, we identified nine distinct cell types. We employed Gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, and the single-cell regulatory network inference and clustering (SCENIC) to determine the enriched pathways and transcriptional activities across seven chondrocyte subpopulations. Furthermore, two cell fates of chondrocyte differentiation were found by trajectory analysis, one was enriched in inflammation-related genes, and the other was related to extracellular matrix (ECM). Additionally, the intercellular interactions of macrophages (MA) and chondrocytes, T cells/natural killer cells (T/NK) and chondrocytes were examined by ligand-receptor pairs analysis, showing the important regulative function of FN1 from MA and CD74 from T/NK during endplate degeneration. Overall, our findings provide novel perspectives on the endplate degeneration at the single-cell level and a whole-transcriptome size., (© 2024. The Author(s).)

Plant height (PH) is a critical agronomic trait in Brassica napus , significantly impacting yield. Consequently, identifying genes associated with plant height is a pivotal objective in oilseed rape breeding. This study employed a combination of bulk segregant analysis sequencing (BSA-seq) and RNA sequencing (RNA-seq) for analysis. A novel quantitative trait locus (QTL), qPH&#95;C02 , was identified between 63,989,634 and 64,945,122 bp on chromosome C02, from which eight candidate genes were screened. The Gene Ontology (GO) analysis revealed enrichment in peroxisomes, while the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated enrichment in the oxidative phosphorylation (OP) pathway. It is hypothesized that the observed differences in plant height and silique length may be attributed to the regulation of peroxidase activity in the OP pathway, which in turn alters plant energy metabolism and controls nutrient uptake. Subsequently, we will further test this hypothesis. The results of this study will contribute to our understanding of the genetic basis for differences in plant height and provide a foundation for the selection and breeding of Brassica napus varieties with desired plant shapes.

Oncostatin M (OSM) is a member of the interleukin-6 (IL-6) family of cytokines and has been found to have anti-inflammatory and pro-inflammatory properties in various cellular and disease contexts. OSM signals through two receptor complexes, one of which includes OSMRβ. Here, we investigated OSM-OSMRβ signaling in adult mouse hematopoietic stem cells (HSCs) using the conditional Osmr fl/fl mouse model B6;129-Osmr tm1.1Nat /J. We crossed Osmr fl/fl mice to interferon-inducible Mx1-Cre, which is robustly induced in adult HSCs. From these mice, we isolated HSCs by flow cytometry, stimulated with recombinant OSM or vehicle for 1 hour, and assessed gene expression changes in control versus Osmr knockout HSCs by RNA-seq. This data may be utilized to investigate OSMRβ -dependent and -independent OSM signaling as well as the transcriptional effects of an IL-6 family cytokine on mouse HSCs to further define its anti-inflammatory versus pro-inflammatory properties., (© 2024. The Author(s).)

Introduction: Osteosarcoma is a common type of bone cancer characterized by a poor prognosis due to its metastatic nature. The tumor microenvironment (TME) plays a critical role in tumor metastasis and therapy response. Therefore, our study aims to explore the metastatic mechanism of osteosarcoma, potentially opening new avenues for cancer treatment., Methods: In this study, we collected data from the GSE152048, GSE14359, and GSE49003 datasets. Differentially expressed genes (DEGs) were identified in osteosarcoma cases with primary and metastatic features using R software and the limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to investigate metastasis-related genes. A protein-protein interaction (PPI) network was established using the STRING database to further analyze these metastasis-associated genes. The abundances of different cell types with a mixed cell population were estimated using the CIBERSORT approach. The scRNA-seq data were analyzed by the Seurat package in R software, and intercellular communications were elucidated using the CellChat R package., Results: In this study, 92 DEGs related to metastasis were identified, including 41 upregulated and 51 downregulated genes in both the GSE14359 and GSE49003 datasets. Metastasis-associated pathways were identified, including those involving the cyclin-dependent protein kinase holoenzyme complex, transferase complex, transferring phosphorus-containing groups, SCF ubiquitin ligase complex, and the serine/threonine protein kinase complex. KEGG and PPI network analyses revealed 15 hub genes, including Skp2, KIF20A, CCNF, TROAP, PHB, CKS1B, MCM3, CCNA2, TRIP13, CENPM, Hsp90AB1, JUN, CKS2, TK1, and KIF4A. Skp2 has been known as an E3 ubiquitin ligase involved in osteosarcoma progression. The proportion of CD8+ T cells was found to be higher in metastatic osteosarcoma tissues, and high expression of PHB was associated with a favorable prognosis in osteosarcoma patients. Additionally, 23 cell clusters were classified into eight cell types, including chondrocytes, MSC, T cells, monocytes, tissue stem cells, neurons, endothelial cells, and macrophages. The 15 hub genes were expressed across various cell types, and interactions between different cell types were observed., Conclusion: Our study reveals the intricate communication between tumor microenvironment components and tumor metastasis in osteosarcoma., 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 Li, Bai, Zhang, Chen and Shang.)

Spermatogenesis is a complex biological process crucial for male reproduction and is characterized by intricate interactions between testicular somatic cells and germ cells. Due to the cellular heterogeneity of the testes, investigating different cell types across developmental stages has been challenging. Single-cell RNA sequencing (scRNA-seq) has emerged as a valuable approach for addressing this limitation. Here, we conducted an unbiased transcriptomic study of spermatogenesis in sexually mature 4-month-old Hezuo pigs using 10× Genomics-based scRNA-seq. A total of 16,082 cells were collected from Hezuo pig testes, including germ cells (spermatogonia (SPG), spermatocytes (SPCs), spermatids (SPTs), and sperm (SP)) and somatic cells (Sertoli cells (SCs), Leydig cells (LCs), myoid cells (MCs), endothelial cells (ECs), and natural killer (NK) cells/macrophages). Pseudo-time analysis revealed that LCs and MCs originated from common progenitors in the Hezuo pig. Functional enrichment analysis indicated that the differentially expressed genes (DEGs) in the different types of testicular germ cells were enriched in the PI3K-AKT, Wnt, HIF-1, and adherens junction signaling pathways, while the DEGs in testicular somatic cells were enriched in ECM-receptor interaction and antigen processing and presentation. Moreover, genes related to spermatogenesis, male gamete generation, sperm part, sperm flagellum, and peptide biosynthesis were expressed throughout spermatogenesis. Using immunohistochemistry, we verified several stage-specific marker genes (such as UCHL1 , WT1 , SOX9 , and ACTA2 ) for SPG, SCs, and MCs. By exploring the changes in the transcription patterns of various cell types during spermatogenesis, our study provided novel insights into spermatogenesis and testicular cells in the Hezuo pig, thereby laying the foundation for the breeding and preservation of this breed.

Both silicon (Si) and zinc (Zn) ions are essential elements to bone health and their mechanisms for promoting osteogenesis have aroused the extensive attention of researchers. Thereinto, the mechanism by which dual ions promote osteogenic differentiation remains to be elucidated. Herein, the effects of Si and Zn ions on the cytological behaviors of mBMSCs were firstly studied. Then, the molecular mechanism of Si-Zn dual ions regulating the osteogenic differentiation of mBMSCs was investigated via transcriptome sequencing technology. In the single-ion system, Si ion at the concentration of 1.5 mM (Si-1.5) had better comprehensive effects of cell proliferation, ALP activity and osteogenesis-related gene expression levels (ALP, Runx2, OCN, Col-I and BSP); Zn ion at the concentration of 50 μM (Zn-50) demonstrated better combining effects of cell proliferation, ALP activity and same osteogenic genes expression levels. In the dual-ion system, the Si (1.5 mM)-Zn (50 μM) group (Si1.5-Zn50) synthetically enhanced ALP activity and osteogenesis genes compared with single-ion groups. Analysis of the transcriptome sequencing results showed that Si ion had a certain effect on promoting the osteogenic differentiation of mBMSCs; Zn ion had a stronger effect of contributing to a better osteogenic differentiation of mBMSCs than that of Si ion; the Si-Zn dual ions had a synergistic enhancement on conducting to the osteogenic differentiation of mBMSCs compared to single ion (Si or Zn). This study offers a blueprint for exploring the regulation mechanism of osteogenic differentiation by dual ions., (© 2024. The Author(s).)

Objective: Osteoporosis (OP) is a disease characterized by decreased bone mass, deteriorated microstructure, and increased fragility and fracture risk. The diagnosis and prevention of OP and its complications have become major public health challenges. Therefore, exploring the complex ecological connections between the immune and skeletal systems may provide new insights for clinical prevention and treatment strategies., Methods: First, we performed single-cell RNA sequencing on human lumbar lamina tissue and conducted clustering and subgroup analysis of quality-controlled single-cell transcriptome data to identify target subgroups. Subsequently, enrichment analysis and pseudotime analysis were performed. In addition, we conducted in-depth studies on the gene regulatory network between different cell subgroups and the communication between bone immune cells., Results: In this study, we identified several cell subgroups that may be involved in the progression of OP. For example, the CCL4 + NKT and CXCL8 + neutrophils subgroups promote OP progression by mediating an inflammatory environment that disrupts bone homeostasis, and the MNDA + Mac subgroup promotes osteoclast differentiation to promote OP. Moreover, the TNFAIP6 + Obl, NR4A2 + B and HMGN2 + erythrocyte subgroups promoted the balance of bone metabolism and suppressed OP. In the cell communication network, Obl closely interacts with immune cell subgroups through the CXCR4-CXCL12, CTGF-ITGB2, and TNFSF14-TNFRSF14 axes., Conclusion: Our research revealed specific subgroups and intercellular interactions that play crucial roles in the pathogenesis of OP, providing potential new insights for more precise therapeutic interventions for OP., 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: Evidence from observational studies indicates that inflammatory proteins play a vital role in Guillain-Barre Syndrome (GBS). Nevertheless, it is unclear how circulating inflammatory proteins are causally associated with GBS. Herein, we conducted a two-sample Mendelian randomization (MR) analysis to systematically explore the causal links of genetically determined systemic inflammatory proteins on GBS., Methods: A total of 8,293 participants of European ancestry were included in a genome-wide association study of 41 inflammatory proteins as instrumental variables. Five MR approaches, encompassing inverse-variance weighted, weighted median, MR-Egger, simple model, and weighted model were employed to explore the causal links between inflammatory proteins and GBS. MR-Egger regression was utilized to explore the pleiotropy. Cochran's Q statistic was implemented to quantify the heterogeneity. Furthermore, we performed single-cell RNA sequencing analysis and predicted potential drug targets through molecular docking technology., Results: By applying MR analysis, four inflammatory proteins causally associated with GBS were identified, encompassing IFN-γ (OR:1.96, 95%CI: 1.02-3.78, P IVW =0.045), IL-7 (OR:1.86, 95%CI: 1.07-3.23, P IVW =0.029), SCGF-β (OR:1.56, 95%CI: 1.11-2.19, P IVW =0.011), and Eotaxin (OR:1.99, 95%CI: 1.01-3.90, P IVW =0.046). The sensitivity analysis revealed no evidence of pleiotropy or heterogeneity. Additionally, significant genes were found through single-cell RNA sequencing analysis and several anti-inflammatory or neuroprotective small molecular compounds were identified by utilizing molecular docking technology., Conclusions: Our MR analysis suggested that IFN-γ, IL-7, SCGF-β, and Eotaxin were causally linked to the occurrence and development of GBS. These findings elucidated potential causal associations and highlighted the significance of these inflammatory proteins in the pathogenesis and prospective therapeutic targets for GBS., 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 Liu and Pan.)

Recent advancements have elucidated the multifaceted roles of the Schlafen (SLFN) family, including SLFN5, SLFN11, SLFN12, SLFN13, and SLFN14, which are implicated in immunological responses. However, little is known about the roles of this gene family in relation to malignancy development. The current study aimed to explore the diagnostic and prognostic potential of Schlafen family genes in colorectal adenocarcinoma (COAD) through bioinformatics analysis. Leveraging advanced bioinformatics tools of bulk RNA-sequencing and single-cell sequencing, we conducted in-depth analyses of gene expressions, functional enrichment, and survival patterns of patients with colorectal cancer compared to normal tissue. Among Schlafen family genes, the transcription levels of SLFN5 in COAD tissues were significantly elevated and correlated with poor survival outcomes. Furthermore, SLFN5 regulated the immune response via Janus kinase (JAK)/signal transduction and activator of transcription (STAT)/interferon (IFN)-alpha/beta signaling. These chemokines in inflammation are associated with diabetes and metabolism, suggesting their involvement in altered cellular energetics for COAD progress. In addition, an immune cell deconvolution analysis indicated a correlation between SLFN5 expression and immune-related cell populations, such as regulatory T cells (Tregs). These findings highlighted the potential clinical significance of SLFN5 in COAD and provided insights into its involvement in the tumor microenvironment and immune regulation. Meanwhile, the drug discovery data of SFLN5 with potential targeted small molecules suggested its therapeutic potential for COAD. Collectively, the current research demonstrated that SFLN5 play crucial roles in tumor development and serve as a prospective biomarker for COAD., Competing Interests: Competing Interests: The authors have declared that no competing interest exists., (© The author(s).)

Objective and Methods: To ascertain the connection between cuproptosis-related genes (CRGs) and the prognosis of hepatocellular carcinoma (HCC) via single-cell RNA sequencing (scRNA-seq) and RNA sequencing (RNA-seq) data, relevant data were downloaded from the GEO and TCGA databases. The differentially expressed CRGs (DE-CRGs) were filtered by the overlaps in differentially expressed genes (DEGs) between HCC patients and normal controls (NCs) in the scRNA-seq database, DE-CRGs between high- and low-CRG-activity cells, and DEGs between HCC patients and NCs in the TCGA database., Results: Thirty-three DE-CRGs in HCC were identified. A prognostic model (PM) was created employing six survival-related genes (SRGs) (NDRG2, CYB5A, SOX4, MYC, TM4SF1, and IFI27) via univariate Cox regression analysis and LASSO. The predictive ability of the model was validated via a nomogram and receiver operating characteristic curves. Research has employed tumor immune dysfunction and exclusion as a means to examine the influence of PM on immunological heterogeneity. Macrophage M0 levels were significantly different between the high-risk group (HRG) and the low-risk group (LRG), and a greater macrophage level was linked to a more unfavorable prognosis. The drug sensitivity data indicated a substantial difference in the half-maximal drug-suppressive concentrations of idarubicin and rapamycin between the HRG and the LRG. The model was verified by employing public datasets and our cohort at both the protein and mRNA levels., Conclusion: A PM using 6 SRGs (NDRG2, CYB5A, SOX4, MYC, TM4SF1, and IFI27) was developed via bioinformatics research. This model might provide a fresh perspective for assessing and managing HCC., (© 2024. The Author(s).)

RNA sequencing (RNA-Seq) has become a widely adopted technique for studying gene expression. However, conventional RNA-Seq analyses rely on gene expression (GE) values that aggregate all the transcripts produced under a single gene identifier, overlooking the complexity of transcript variants arising from different transcription start sites or alternative splicing. Transcript variants may encode proteins with diverse functional domains, or noncoding RNAs. This study explored the implications of neglecting transcript variants in RNA-Seq analyses. Among the 1334 transcription factor (TF) genes expressed in mouse embryonic stem (ES) or trophoblast stem (TS) cells, 652 were differentially expressed in TS cells based on GE values (365 upregulated and 287 downregulated, ≥absolute 2-fold changes, false discovery rate (FDR) p -value ≤ 0.05). The 365 upregulated genes expressed 883 transcript variants. Further transcript expression (TE) based analyses identified only 174 (<20%) of the 883 transcripts to be upregulated. The remaining 709 transcripts were either downregulated or showed no significant changes. Meanwhile, the 287 downregulated genes expressed 856 transcript variants and only 153 (<20%) of the 856 transcripts were downregulated. The other 703 transcripts were either upregulated or showed no significant change. Additionally, the 682 insignificant TF genes (GE values < absolute 2-fold changes and/or FDR p-values > 0.05) between ES and TS cells expressed 2215 transcript variants. These included 477 (>21%) differentially expressed transcripts (276 upregulated and 201 downregulated, ≥absolute 2-fold changes, FDR p-value ≤ 0.05). Hence, GE based RNA-Seq analyses do not represent accurate expression levels due to divergent transcripts expression from the same gene. Our findings show that by including transcript variants in RNA-Seq analyses, we can generate a precise understanding of a gene's functional and regulatory landscape; ignoring the variants may result in an erroneous interpretation.

Granulosa cells (GCs) play crucial roles in oocyte maturation. Through gap junctions and extracellular vesicles, they mediate the exchange of molecules such as microRNAs and messenger RNAs. Different ovarian cell types exhibit unique gene expression profiles, reflecting their specialized functions and stages. By combining RNA-seq data from various cell types forming the follicle, we aimed at capturing a wide range of expression patterns, offering insights into the functional diversity and complexity of the transcriptome regulation across GCs. Herein, we performed an integrated bioinformatics analysis of RNA sequencing datasets present in public databases, with a unique and standardized workflow., By combining the data from different studies, we successfully increased the robustness and reliability of our findings and discovered novel genes, miRNAs, and signaling pathways associated with GCs function and oocyte maturation. Moreover, our results provide a valuable resource for further wet-lab research on GCs biology and their impact on oocyte development and competence., (© 2024. The Author(s).)

Background: With the widespread use of antibiotics, more attention has been paid to their side effects. We paid extra attention to the impact of antibiotics on children's bodies. Therefore, we analyzed the characteristic changes in the gut microbiota of children after antibiotic treatment to explore the pathogenesis of antibiotic-associated diseases in more depth and to provide a basis for diagnosis and treatment., Methods: We recruited 28 children with bronchopneumonia in the western district of Zhuhai, China, and divided them into three treatment groups based on antibiotic type. We took stool samples from children before and 3-5 days after antibiotic treatment. 16S rRNA gene sequencing was used to analyze the effects of antibiotic therapy on the gut microbiota of children. Continuous nonparametric data are represented as median values and analyzed using the Wilcoxon rank-sum test., Results: While alpha diversity analysis found no significant changes in the mean abundance of the gut microbiota of children after a short course of antibiotic treatment, beta diversity analysis demonstrated significant changes in the composition and diversity of the gut microbiota of children even after a short course of antibiotic therapy. We also found that meloxicillin sulbactam can inhibit the growth of Proteobacteria, Bacteroidetes, and Verrucomicrobia, ceftriaxone inhibits Verrucomicrobia and Bacteroides, and azithromycin inhibits Fusobacteria, Actinobacteria, Proteobacteria, and Verrucomicrobia. We further performed a comparative analysis at the genus level and found significantly different clusters in each group. Finally, we found that azithromycin had the greatest effect on the metabolic function of intestinal microbiota, followed by ceftriaxone, and no significant change in the metabolic process of intestinal microbiota after meloxicillin sulbactam treatment., Conclusions: Antibiotic treatment significantly affects the diversity of intestinal microbiota in children, even after a short course of antibiotic treatment. Different classes of antibiotics affect diverse microbiota primarily, leading to varying alterations in metabolic function. Meanwhile, we identified a series of intestinal microbiota that differed significantly after antibiotic treatment. These groups of microbiota could be used as biomarkers to provide an additional basis for diagnosing and treating antibiotic-associated diseases., (© 2024. The Author(s).)

Application of polyethylene glycol (PEG) to a peripheral nerve injury at the time of primary neurorrhaphy is thought to prevent Wallerian degeneration via direct axolemma fusion. The molecular mechanisms of nerve fusion and recovery are unclear. Our study tested the hypothesis that PEG alters gene expression in neural and muscular environments as part of its restorative properties. Lewis rats underwent unilateral sciatic nerve transection with immediate primary repair. Subjects were randomly assigned to receive either PEG treatment or standard repair at the time of neurorrhaphy. Samples of sciatic nerve distal to the injury and tibialis muscle at the site of innervation were harvested at 24 hours and 4 weeks postoperatively. Total RNA sequencing and subsequent bioinformatics analyses were used to identify significant differences in differentially expressed genes (DEGs) and their related biological pathways (p<0.05) in PEG-treated subjects compared to non-PEG controls. No significant DEGs were identified in PEG-treated sciatic nerve compared to controls after 24 hours, but 1,480 DEGs were identified in PEG-treated tibialis compared to controls. At 4 weeks, 918 DEGs were identified in PEG-treated sciatic nerve, whereas only 3 DEGs remained in PEG-treated tibialis compared to controls. DEGs in sciatic were mostly upregulated (79%) and enriched in pathways present during nervous system development and growth, whereas DEGs in muscle were mostly downregulated (77%) and related to inflammation and tissue repair. Our findings indicate that PEG application during primary neurorrhaphy leads to significant differential gene regulation in the neural and muscular environment that is associated with improved functional recovery in animals treated with PEG compared to sham non-PEG controls. A detailed understanding of key molecules underlying PEG function in recovery after peripheral nerve repair may facilitate amplification of PEG effects through systemic or focal treatments at the time of neurotmesis., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Weiss et al. 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 polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh&#95;A02G0645&#95;N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Purpose: Transdifferentiation of corneal fibroblasts to myofibroblasts in the stroma is a central mechanistic event in corneal wound healing. This study sought to characterize genes and pathways influencing transdifferentiation of human corneal fibroblasts (hCSFs) to human corneal myofibroblasts (hCMFs) using RNA sequencing (RNA-seq) to develop comprehensive mechanistic information and identify newer targets for corneal fibrosis management., Methods: Primary hCSFs were derived from donor human corneas. hCMFs were generated by treating primary hCSFs with transforming growth factor β1 (TGFβ1; 5 ng/mL) for 72 hours under serum-free conditions. RNA was extracted using the RNeasy Plus Mini Kit and subjected to RNA-seq analysis after quality control testing. Differential gene expression, pathway enrichment, and protein-protein network analyses were performed using DESeq2, GSEA/PANTHER/Reactome, and Cytoscape/cytoHubba, respectively., Results: RNA-seq analysis of hCMFs and hCSFs identified 3843 differentially expressed genes and transcripts (adjusted P < 0.05). The log(fold change) ≥ ±1.5 filter showed 816 upregulated and 739 downregulated genes between two cell types. Pathway enrichment analysis showed the highest normalized enrichment score for epithelial-to-mesenchymal transition (5.569), followed by mTORC1 signaling (2.949), angiogenesis (2.176), and TGFβ signaling (2.008). Protein-protein interaction network analysis identified the top 20 nodes influencing corneal myofibroblast development. The expression of a novel MXRA5 in corneal stroma and its association with corneal fibrosis was verified by real-time quantitative reverse transcription PCR and immunofluorescence. RNA-seq and gene count files were submitted to the NCBI Gene Expression Omnibus (GSE260476)., Conclusions: This study identified several novel genes involved in myofibroblast development, offering potential targets for developing newer therapeutic strategies for corneal fibrosis.

Background: The Testis is an important reproductive organ in male mammals and the site for spermatogenesis, androgen synthesis, and secretion. Non-coding RNAs (ncRNAs) play an important regulatory role in various biological processes. However, the regulatory role of ncRNAs in the development of yak testes and spermatogenesis remains largely unclear., Result: In this study, we compared the expression profiles of circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in yak testicular tissue samples collected at 6 months (Y6M), 18 months (Y18M), and 4 years (Y4Y). Using RNA sequencing (RNA-Seq), we observed a significant difference in the expression patterns of ncRNAs in the samples collected at different testicular development stages. Twenty-two differentially expressed (DE) circRNAs, 69 DE miRNAs, and 64 DE mRNAs were detected in Y6M, Y18M, and Y4Y testicular samples, respectively. The results of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the source genes of DE circRNAs, predicted target genes of DE miRNAs, and DE mRNAs were specifically associated with signaling pathways and GO terms that were related to sperm synthesis, sperm vitality, and testicular development, such as cell cycle, Wnt signaling pathway, MAPK signaling pathway, GnRH signaling pathway, and spermatogenesis. The analysis of the circRNA-miRNA-mRNA network revealed that some DE ncRNAs, including miR-574, miR-449a, CDC42, and CYP11A1, among others, may be involved in testicular spermatogenesis. Concurrently, various circRNA-miRNA interaction pairs were observed., Conclusion: Our findings provide a database of circRNAs, miRNAs, and mRNAs expression profiles in testicular tissue of yaks at different developmental stages and a detailed understanding of the regulatory network of ncRNAs in yak testicular development and provide data that can help elucidate the molecular mechanisms underlying yak testicular development., (© 2024. The Author(s).)

Hyperforin is the compound responsible for the effectiveness of St. John's wort (Hypericum perforatum) as an antidepressant, but its complete biosynthetic pathway remains unknown. Gene discovery based on co-expression analysis of bulk RNA-sequencing data or genome mining failed to discover the missing steps in hyperforin biosynthesis. In this study, we sequenced the 1.54-Gb tetraploid H. perforatum genome assembled into 32 chromosomes with the scaffold N50 value of 42.44 Mb. By single-cell RNA sequencing, we identified a type of cell, "Hyper cells", wherein hyperforin biosynthesis de novo takes place in both the leaves and flowers. Through pathway reconstitution in yeast and tobacco, we identified and characterized four transmembrane prenyltransferases (HpPT1-4) that are localized at the plastid envelope and complete the hyperforin biosynthetic pathway. The hyperforin polycyclic scaffold is created by a reaction cascade involving an irregular isoprenoid coupling and a tandem cyclization. Our findings reveal how and where hyperforin is biosynthesized, enabling synthetic-biology reconstitution of the complete pathway. Thus, this study not only deepens our comprehension of specialized metabolism at the cellular level but also provides strategic guidance for elucidation of the biosynthetic pathways of other specializied metabolites in plants., (Copyright © 2024 The Author. Published by Elsevier Inc. All rights reserved.)

Myelodysplastic syndromes (MDS) are myeloid malignancies with heterogeneous genotypes and phenotypes, characterized by ineffective haematopoiesis and a high risk of progression towards acute myeloid leukaemia (AML). Prognosis for patients treated with hypomethylating agents (HMAs), as is azacytidine, the main drug used as frontline therapy for MDS is mostly based on cytogenetics and next generation sequencing (NGS) of the initial myeloid clone. Although the critical influence of the epigenetic landscape upon cancer cells survival and development as well on tumour environment establishment is currently recognized and approached within current clinical practice in MDS, the heterogenous response of the patients to epigenetic therapy is suggesting a more complex mechanism of action, as is the case of RNA methylation. In this sense, the newly emerging field of epitranscriptomics could provide a more comprehensive perspective upon the modulation of gene expression in malignancies, as is the proof-of-concept of MDS. We initially did RNA methylation sequencing on MDS patients (n = 6) treated with azacytidine and compared responders with non-responders. Afterwards, the genes identified were assessed in vitro and afterwards validated on a larger cohort of MDS patients treated with azacytidine (n = 58). Our data show that a more accurate prognosis could be based on analysing the methylome and thus we used methylation sequencing to differentially split high-grade MDS patients with identical demographical and cytogenetic features, between azacytidine responders and non-responders., (© 2024 The Author(s). Journal of Cellular and Molecular Medicine published by Foundation for Cellular and Molecular Medicine and John Wiley & Sons Ltd.)

Camptotheca acuminata Decne., a significant natural source of the anticancer drug camptothecin (CPT), synthesizes CPT through the monoterpene indole alkaloid (MIA) pathway. In this study, we used single-cell RNA sequencing (scRNA-seq) to generate datasets encompassing over 60,000 cells from C. acuminata shoot apexes and leaves. After cell clustering and annotation, we identified five major cell types in shoot apexes and four in leaves. Analysis of MIA pathway gene expression revealed that most of them exhibited heightened expression in proliferating cells (PCs) and vascular cells (VCs). In contrast to MIA biosynthesis in Catharanthus roseus, CPT biosynthesis in C. acuminata did not exhibit multicellular compartmentalization. Some putative genes encoding enzymes and transcription factors (TFs) related to the biosynthesis of CPT and its derivatives were identified through co-expression analysis. These include 19 cytochrome P450 genes, 8 O-methyltransferase (OMT) genes, and 62 TFs. Additionally, these pathway genes exhibited dynamic expression patterns during VC and EC development. Furthermore, by integrating gene and transposable element (TE) expression data, we constructed novel single-cell transcriptome atlases for C. acuminata. This approach significantly facilitated the identification of rare cell types, including peripheral zone cells (PZs). Some TE families displayed cell type specific, tissue specific, or developmental stage-specific expression patterns, suggesting crucial roles for these TEs in cell differentiation and development. Overall, this study not only provides novel insights into CPT biosynthesis and spatial-temporal TE expression characteristics in C. acuminata, but also serves as a valuable resource for further comprehensive investigations into the development and physiology of this species., (© 2024 Scandinavian Plant Physiology Society.)

Colorectal cancer (CRC) exhibits considerable heterogeneity on tumour location. However, there is still a lack of comprehensive annotation regarding the characteristics and differences between the left-sided (L-CRC) and right-sided (R-CRC) CRC. Here, we performed single-cell RNA sequencing (scRNA-seq) on immune and stromal cells from 12 L-CRC and 10 R-CRC patients. We found that L-CRC exhibited stronger tumour invasion and poor prognosis compared with R-CRC. In addition, functional enrichment analysis of a normal cohort showed that fibroblasts of left colon are associated with tumour-related pathways. This suggested that the heterogeneity observed in both L-CRC and R-CRC may be influenced by the specific location within the colon itself. Further, we identified a potentially novel MYH11+ cancer-associated fibroblast (CAF) subset predominantly enriched in L-CRC. Moreover, we found that MYH11+ CAFs may promote tumour migration via interacting with macrophages, and was associated with poor prognosis in CRC. In summary, our study revealed the crucial role of MYH11+ CAFs in predicting a poor prognosis, thereby contributing valuable insights to the exploration of heterogeneity in L-CRC and R-CRC., (© 2024 The Author(s). Journal of Cellular and Molecular Medicine published by Foundation for Cellular and Molecular Medicine and John Wiley & Sons Ltd.)

Background & Aims: Cystic fibrosis (CF) is considered a multisystemic disorder in which CF-associated liver disease (CFLD) is the third most common cause of mortality. Currently, no effective treatment is available for CFLD because its pathophysiology is still unclear. Interestingly, CFLD exhibits identical vascular characteristics as non-cirrhotic portal hypertension, recently classified as porto-sinusoidal vascular disorders (PSVD)., Methods: Since endothelial cells (ECs) are an important component in PSVD, we performed single-cell RNA sequencing (scRNA-seq) on four explant livers from CFLD patients to identify differential endothelial characteristics which could contribute to the disease. We comprehensively characterized the endothelial compartment and compared it with publicly available scRNA-seq datasets from cirrhotic and healthy livers. Key gene signatures were validated ex vivo on patient tissues., Results: We found that ECs from CF liver explants are more closely related to healthy than cirrhotic patients. In CF patients we also discovered a distinct population of liver sinusoidal ECs-coined CF LSECs-upregulating genes involved in the complement cascade and coagulation. Finally, our immunostainings further validated the predominant periportal location of CF LSECs., Conclusions: Our work showed novel aspects of human liver ECs at the single-cell level thereby supporting endothelial involvement in CFLD, and reinforcing the hypothesis that ECs could be a driver of PSVD. Therefore, considering the vascular compartment in CF and CFLD may help developing new therapeutic approaches for these diseases., (© 2024 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Mesenchymal cells within theplacental villi play a crucial role in shaping the morphology of branching structures and driving the development of blood vessels. However, the markers and functions of placental villous pericytes (PVPs) as distinct subgroups of placental villous mesenchymal cells, remain unclear. Therefore, in this study, the markers and functions of PVPs were investigated. Single-cell sequencing data from the first-trimester placental villi was obtained and the Seurat tool was used to identify PVP markers. Gene Ontology (GO) analysis of specific genes was performed using the DAVID database. The Cellchat tool was employed to investigate the interaction signals between the PVPs and other cells. Expression of the PVP markers was confirmed using immunofluorescence. Presence of extracellular vesicles in the placental villous mesenchyme and PVPs was examined by transmission electron microscopy. Our findings revealed that renin (REN) and amphiregulin (AREG)-positive fibroblasts in the placental villi specifically expressed several classic pericyte markers. In the first trimester, certain conserved functions of pericytes were observed and they displayed tissue-specific functions such as in the integrin-mediated signaling pathway and extracellular exosomes. Moreover, the placental villous mesenchyme was found to be rich in extracellular vesicles. AREG is specifically transcribed in the first trimester PVPs, however, its protein was located in syncytiotrophoblasts. These insights contribute to a comprehensive understanding of early placental development and offer new therapeutic targets for placenta-derived pregnancy complications., Competing Interests: Declaration of Competing Interest The author declares that he has no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Society for Biology of Reproduction & the Institute of Animal Reproduction and Food Research of Polish Academy of Sciences in Olsztyn. Published by Elsevier B.V. All rights reserved.)

Purpose: We aimed to investigate the transcriptomic alterations that occur in the subacromial bursa (SAB) following degenerative or traumatic shoulder diseases., Materials and Methods: RNA sequencing was employed to evaluate the transcriptomic alterations of the SAB in individuals afflicted with degenerative rotator cuff tear (RCT), traumatic RCT and proximal humerus fracture (PHF). To gain insights into the biological significance of differentially expressed genes (DEGs), we conducted an enrichment analysis utilizing Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. We further utilized single-cell RNA sequencing datasets of SAB from a recently published study to explore the associated cellular dynamics and alterations., Results: We detected 1,790 up-regulated and 1,964 down-regulated DEGs between degenerative RCT and PHF, 2,085 up-regulated and 1,919 down-regulated DEGs between degenerative RCT and traumatic RCT, and 20 up-regulated and 12 down-regulated DEGs between traumatic RCT and PHF. Given the similar expression pattern between traumatic RCT and PHF, they were integrated as the traumatic group. In comparison with the traumatic group, 1,983 up-regulated and 2,205 down-regulated DEGs were detected in degenerative SAB. Enrichment analysis of up-regulated DEGs uncovered an elevated inflammatory and immunologic responses in degenerative SAB. Single-cell transcriptomic analysis revealed macrophage represented the immune cell with the most DEGs between the degenerative and traumatic RCT., Conclusion: Our results revealed that the SAB in degenerative RCT exhibited a different transcriptional signature compared to that in traumatic RCT, and enrichment analysis showed immunologic and inflammatory activations. Macrophages may play a fundamental role in this process.