Search

In the human and veterinary fields, oral vaccines generate considerable interest. In dogs, these vaccines are newly developed, and understanding their mechanisms is crucial. Mesenteric lymph nodes (MLNs) and Peyer's patches (PPs) are important sites for gastrointestinal mucosal induction, yet canine MLNs lack sufficient information. To address this, we collected MLN samples from healthy dogs, performed flow cytometry to characterize immune cells, and conducted single-cell RNA sequencing (scRNA-seq) to explore subpopulations, particularly B and T lymphocytes. This effort enabled the characterization of canine MLN's main cell populations and the construction of a predictive atlas, as well as the identification of particularities of this area., (© 2024. The Author(s).)

Objective: The aim of this research was to gain a thorough understanding of the processes involved in cell communication and discover potential indicators for treating multiple myeloma (MM) through the use of single-cell RNA sequencing (scRNA-seq). And explored the expression of multiple myeloma-related subgroups on metal ion-related pathways to explore the relationship between MM and metal ions., Methods: We performed a fair examination using single-cell RNA sequencing on 32 bone marrow specimens collected from 22 individuals at different points of MM advancement and 9 individuals without any health issues. To analyze the scRNA-seq data, we employed advanced computational algorithms, including Slingshot, Monocle2, and other methodologies. Specifically, Slingshot and Monocle2 enabled us to simulate the biological functionalities of different cell populations and map trajectories of cell developmental pathways. Additionally, we utilized the UMAP algorithm, a powerful dimension reduction technique, to cluster cells and identify genes that were differentially expressed across clusters., Results: Our study revealed distinct gene expression patterns and molecular pathways within each patient, which exhibited associations with disease progression. The analysis provided insights into the tumor microenvironment (TME), intra- and inter-patient heterogeneity, and cell-cell interactions mediated by ligand-receptor signaling. And found that multiple myeloma-related subgroups were expressed higher levels in MMP and TIMP pathways, there were some associations., Conclusion: Our study presents a fresh perspective for future research endeavors and clinical interventions in the field of MM. The identified gene expression patterns and molecular pathways hold immense potential as therapeutic targets for the treatment of multiple myeloma. The utilization of scRNA-seq technology has significantly contributed to a more precise understanding of the complex cellular processes and interactions within MM. Through these advancements, we are now better equipped to unravel the underlying mechanisms driving the development and progression of this complex disease., Competing Interests: The authors declare that the research is conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Li, Lin, Zhao, Huang, Fan, Cen and Ma.)

Background: Intracerebral hemorrhage (ICH) is a severe stroke subtype with high morbidity, disability, and mortality rates. Currently, no biomarkers for ICH are available for use in clinical practice. We aimed to explore the roles of RNAs in ICH pathogenesis and identify potential diagnostic biomarkers., Methods: We collected 233 individual blood samples from two independent cohorts, including 64 patients with ICH, 59 patients with ischemic stroke (IS), 60 patients with hypertension (HTN) and 50 healthy controls (CTRL) for RNA sequencing. Differentially expressed genes (DEGs) analysis, gene set enrichment analysis (GSEA), and weighted correlation network analysis (WGCNA) were performed to identify ICH-specific modules. The immune cell composition was evaluated with ImmuneCellAI. Multiple machine learning algorithms to select potential biomarkers for ICH diagnosis, and further validated by quantitative real-time polymerase chain reaction (RT-PCR). Receiver operating characteristic (ROC) curve analysis and decision curve analysis (DCA) were performed to evaluate the diagnostic value of the signature for ICH. Finally, we generated M1 and M2 macrophages to investigate the expression of candidate genes., Results: In both cohorts, 519 mRNAs and 131 lncRNAs were consistently significantly differentially expressed between ICH patients and HTN controls. Gene function analysis suggested that immune system processes may be involved in ICH pathology. ImmuneCellAI analysis revealed that the abundances of 11 immune cell types were altered after ICH in both cohorts. WGCNA and GSEA identified 18 immune-related DEGs. Multiple algorithms identified an RNA panel (CKAP4, BCL6, TLR8) with high diagnostic value for discriminating ICH patients from HTN controls, CTRLs and IS patients (AUCs: 0.93, 0.95 and 0.82; sensitivities: 81.3%, 84.4% and 75%; specificities: 100%, 96% and 79.7%, respectively). Additionally, CKAP4 and TLR8 mRNA and protein levels decreased in RAW264.7 M1 macrophages and increased in RAW264.7 M2 macrophages, while BCL6 expression increased in M1 macrophages but not in M2 macrophages, which may provide potential therapeutic targets for ICH., Conclusions: This study demonstrated that the expression levels of lncRNAs and mRNAs are associated with ICH, and an RNA panel (CKAP4, BCL6, TLR8) was developed as a potential diagnostic tool for distinguishing ICH from IS and controls, which could provide useful insight into ICH diagnosis and pathogenesis., 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 Bai, Liu, Wang, Sun, Wang, Liu, Hao, Zhou, Yuan and Liu.)

Colorectal cancer (CRC) is one of most common tumors worldwide, causing a prominent global health burden. Cell senescence is a complex physiological state, characterized by proliferation arrest. Here, we investigated the role of cellular senescence in the heterogeneity of CRC. Based on senescence-associated genes, CRC samples were classified into different senescence patterns with different survival, cancer-related biological processes and immune cell infiltrations. A senescence-related model was then developed to calculate the senescence-related score to comprehensively explore the heterogeneity of each CRC sample such as stromal activities, immunoreactivities and drug sensitivity. Single-cell analysis revealed there were different immune cell infiltrations between low and high senescence-related model genes enrichment groups, which was confirmed by multiplex immunofluorescence staining. Pseudotime analysis indicated model genes play a pivotal role in the evolution of B cells. Besides, intercellular communication modeled by NicheNet showed tumor cells with higher enrichment of senescence-related model genes highly expressed CXCL2/3 and CCL3/4, which attracted immunosuppressive cell infiltration and promoted tumor metastasis. Finally, top 6 hub genes were identified from senescence-related model genes by PPI analysis. And RT-qPCR revealed the expression differences of hub genes between normal and CRC cell lines, indicating to some extent the clinical practicability of senescence-related model. To sum up, our study explores the impact of cellular senescence on the prognosis, TME and treatment of CRC based on senescence patterns. This provides a new perspective for CRC treatment., (© 2024. The Author(s).)

Oral squamous cell carcinoma (OSCC) has an extremely poor prognosis. Recent studies have suggested that mitophagy-related genes (MRGs) are closely correlated with the development and occurrence of cancer, but the role they play in oral cancer has not yet been explained.We conducted a comprehensive analysis of integrated single-cell and bulk RNA sequencing (RNA-seq) data retrieved from Gene Expression Omnibus (GEO) datasets and The Cancer Genome Atlas (TCGA) database. Multiple methods were combined to provide a comprehensive understanding of the genetic expression patterns and biology of OSCC, such as analysis of pseudotime series, CellChat cell communication, immune infiltration, Gene Ontology (GO), LASSO Cox regression, gene set variation analysis (GSVA), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment analysis (GSEA), Tumor Mutation Burden (TMB) and drug sensitivity assessments. The findings of this study demonstrated significantly greater activity of MRGs in NK cells than in other cells in OSCC. A reliable prognostic model was developed using 12 candidate genes strongly associated with mitochondrial autophagy. T stage, N stage and risk score were revealed as independent prognostic factors. Distinctively enriched pathways and immune cells were observed in different risk groups. Notably, low-risk patients were more responsive to chemotherapy. In addition, a nomogram model with excellent predictive ability was established by combining the risk scores and clinical features. The activity of MRGs suggest the potential for the development of new targeted therapies. The construction of a robust prognostic model also provides reference value for individualized prediction and clinical decision-making in patients with OSCC., (© 2024. The Author(s).)

Rhizopus microsporus is a species in the order Mucorales that is known to cause mucormycosis, but it is poorly understood as a host of viruses. Here, we examined 25 clinical strains of R. microsporus for viral infection with a conventional double-stranded RNA (dsRNA) assay using agarose gel electrophoresis (AGE) and the recently established fragmented and primer-ligated dsRNA sequencing (FLDS) protocol. By AGE, five virus-infected strains were detected. Then, full-length genomic sequences of 12 novel RNA viruses were revealed by FLDS, which were related to the families Mitoviridae, Narnaviridae, and Endornaviridae, ill-defined groups of single-stranded RNA (ssRNA) viruses with similarity to the established families Virgaviridae and Phasmaviridae, and the proposed family "Ambiguiviridae." All the characterized viruses, except a potential phasmavirid with a negative-sense RNA genome, had positive-sense RNA genomes. One virus belonged to a previously established species within the family Mitoviridae , whereas the other 11 viruses represented new species or even new genera. These results show that the fungal pathogen R. microsporus harbors diverse RNA viruses and extend our understanding of the diversity of RNA viruses in the fungal order Mucorales, division Mucoromycota. Identifying RNA viruses from clinical isolates of R. microsporus may expand the repertoire of natural therapeutic agents for mucormycosis in the future.IMPORTANCEThe diversity of mycoviruses in fungal hosts in the division Mucoromycota has been underestimated, mainly within the species Rhizopus microsporus . Only five positive-sense RNA genomes had previously been discovered in this species. Because current sequencing methods poorly complete the termini of genomes, we used fragmented and primer-ligated double-stranded RNA sequencing to acquire the full-length genomes. Eleven novel mycoviruses were detected in this study, including the first negative-sense RNA genome reported in R. microsporus . Our findings extend the understanding of the viral diversity in clinical strains of Mucoromycota, may provide insights into the pathogenesis and ecology of this fungus, and may offer therapeutic options., Competing Interests: The authors declare no conflict of interest.

Juvenile localized and systemic scleroderma are rare autoimmune diseases which cause significant disability and morbidity in children. The mechanisms driving juvenile scleroderma remain unclear, necessitating further cellular and molecular level studies. The Visium CytAssist spatial transcriptomics (ST) platform, which preserves the spatial location of cells and simultaneously sequences the whole transcriptome, was employed to profile the histopathological slides from skin lesions of juvenile scleroderma patients. (1) Spatial domains were identified from ST data and exhibited strong concordance with the pathologist's annotations of anatomical structures. (2) The integration of paired ST data and single-cell RNA sequencing (scRNA-seq) from the same patients validated the comparable accuracy of the two platforms and facilitated the estimation of cell type composition in ST data. (3) The pathologist-annotated immune infiltrates, such as perivascular immune infiltrates, were clearly delineated by the ST analysis, underscoring the biological relevance of the findings. This is the first study utilizing spatial transcriptomics to investigate skin lesions in juvenile scleroderma patients. The validity of the ST data was corroborated by gene expression analyses and the pathologist's assessments. Integration with scRNA-seq data facilitated the cell type-level analysis and validation. Analyses of immune infiltrates through combined ST data and pathological review enhances our understanding of the pathogenesis of juvenile scleroderma.

Single-cell RNA sequencing (scRNA-seq) technology has emerged as a powerful tool for dissecting cellular heterogeneity and understanding the intricate biology of diseases, including cancer. Endometrial cancer (EC) stands out as the most prevalent gynecological malignancy in Europe and the second most diagnosed worldwide, yet its cellular complexity remains poorly understood. In this review, we explore the contributions of scRNA-seq studies to shed light on the tumor cells and cellular landscape of EC. We discuss the diverse tumoral and microenvironmental populations identified through scRNA-seq, highlighting the implications for understanding disease progression. Furthermore, we address potential limitations inherent in scRNA-seq studies, such as technical biases and sample size constraints, emphasizing the need for larger-scale research encompassing a broader spectrum of EC histological subtypes. Notably, a significant proportion of scRNA-seq analyses have focused on primary endometrioid carcinoma tumors, underscoring the need to incorporate additional histological and aggressive types to comprehensively capture the heterogeneity of EC. By critically evaluating the current state of scRNA-seq research in EC, this review underscores the importance of advancing towards more comprehensive studies to accelerate our understanding of this complex disease., Competing Interests: JC reports the following: Consulting/Advisor: Roche, AstraZeneca, Seattle Genetics, Daiichi Sankyo, Lilly, Merck Sharp&Dohme, Leuko, Bioasis, Clovis Oncology, Boehringer Ingelheim, Ellipses, Hibercell, BioInvent, Gemoab, Gilead, Menarini, Zymeworks, Reveal Genomics, Scorpion Therapeutics, Expres2ion Biotechnologies, Jazz Pharmaceuticals, Abbvie, BridgeBio, Biontech. Honoraria: Roche, Novartis, Eisai, Pfizer, Lilly, Merck Sharp&Dohme, Daiichi Sankyo, Astrazeneca, Gilead, Steamline Therapeutics. Research funding to the Institution: Roche, Ariad pharmaceuticals, AstraZeneca, Baxalta GMBH/Servier Affaires, Bayer healthcare, Eisai, F.Hoffman-La Roche, Guardanth health, Merck Sharp&Dohme, Pfizer, Piqur Therapeutics,Iqvia, Queen Mary University of London. Stock: MAJ3 Capital, Leuko relative. Travel, accommodation, expenses: Roche, Novartis, Eisai, Pfizer, Daiichi Sankyo, Astrazeneca, Gilead, Merck Sharp&Dhome, Steamline Therapeutics. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 González-Martínez, Pérez-Mies, Cortés and Palacios.)

Background: Obesity is a global health issue with detrimental effects on various human organs, including the reproductive system. Observational human data and several lines of animal experimental data suggest that maternal obesity impairs ovarian function and early embryo development, but the precise pathogenesis remains unclear., Methods: We established a high-fat diet (HFD)-induced obese female mouse model to assess systemic metabolism, ovarian morphology, and oocyte function in mice. For the first time, this study employed single-cell RNA sequencing to explore the altered transcriptomic landscape of preimplantation embryos at different stages in HFD-induced obese mice. Differential gene expression analysis, enrichment analysis and protein-protein interactions network analysis were performed., Results: HFD-induced obese female mice exhibited impaired glucolipid metabolism and insulin resistance. The ovaries of HFD mice had a reduced total follicle number, an increased proportion of atretic follicles, and irregular granulosa cell arrangement. Furthermore, the maturation rate of embryonic development by in vitro fertilization of oocytes was significantly decreased in HFD mice. Additionally, the transcriptional landscapes of preimplantation embryos at different stages in mice induced by different diets were significantly distinguished. The maternal-to-zygotic transition was also affected by the failure to remove maternal RNAs and to turn off zygotic genome expression., Conclusions: HFD-induced obesity impaired ovarian morphology and oocyte function in female mice and further led to alterations in the transcriptional landscape of preimplantation embryos at different stages of HFD mice., (© 2024. The Author(s).)

Lung adenocarcinoma (LUAD), the most common subtype of lung cancer globally, has seen improved prognosis with advancements in diagnostic, surgical, radiotherapy, and molecular therapy techniques, while its 5-year survival rate remains low. Molecular biomarkers provide prognostic value. Oxidative stress factors, such as reactive nitrogen species and ROS, are crucial in various stages of tumor progression, influencing cell transformation, proliferation, angiogenesis, and metastasis. ROS demonstrate dual roles, affecting tumor cells, hypoxia sensitivity, and the microenvironment. Comprehensive analysis of oxidative stress in LUAD has not been conducted to date. Therefore, we systematically investigated the regulatory patterns of oxidative stress in LUAD based on oxidative stress-related genes and correlated these patterns with cellular infiltration characteristics of the tumor immune microenvironment. The model utilizes single-factor Cox analysis to screen key differential genes with prognostic value and employs least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis to construct a prognostic-related prediction model. Ten candidate genes were selected based on this model. The risk score was constructed using the coefficients and expression levels of these ten genes. Furthermore, the impact of this risk score on overall survival (OS) was determined. Two genes with the most significant differential expression, SFTPB and S100P, were selected through qRT-PCR. Cell experiments including CCK-8, Edu, transwell assays confirmed their effects on lung cancer cells growth, consistent with the results of bioinformatics analysis. These findings suggested that this model held potential clinical value for evaluating the prognosis of lung adenocarcinoma., 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: The occurrence and progression of hepatocellular carcinoma (HCC) are significantly affected by DNA damage response (DDR). Exploring DDR-related biomarkers can help predict the prognosis and immune characteristics of HCC., Methods: First, the single-cell RNA sequencing (scRNA-seq) dataset GSE242889 was processed and performed manual annotation. Then we found the marker genes of DDR-active subgroups based on "AUCell" algorithm. The "Limma" R package was used to identify differentially expressed genes (DEGs) between tumor and normal samples of HCC. The risk prognostic model was constructed by filtering genes using univariate Cox and LASSO regression analyses. Finally, the signatures were analyzed for immune infiltration, gene mutation, and drug sensitivity. Last but not least, KPNA2, which had the largest coefficient in our model was validated by experiments including western blot, MTT, colony formation and γ-H2AX assays., Results: We constructed a prognostic model based on 5 DDR marker genes including KIF2C, CDC20, KPNA2, UBE2S and ADH1B for HCC. We also proved that the model had an excellent performance in both training and validation cohorts. Patients in the high-risk group had a poorer prognosis, different immune features, gene mutation frequency, immunotherapy response and drug sensitivity compared with the low-risk group. Besides, our experimental results proved that KPNA2 was up-regulated in liver cancer cells than in hepatocytes. More importantly, the knockdown of KPNA2 significantly inhibited cell variability, proliferation and promoted DNA damage., Conclusions: We innovatively integrated scRNA-seq and bulk RNA sequencing to construct the DDR-related prognostic model. Our model could effectively predict the prognosis, immune landscape and therapy response of HCC., 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. Published by Elsevier B.V.)

Background: Increased oxidative stress (OS) activity following intracerebral hemorrhage (ICH) had significantly impacting patient prognosis. Identifying optimal genes associated with OS could enhance the understanding of OS after ICH., Methods: We employed single-cell RNA sequencing (scRNA-seq) to investigate the heterogeneity of OS across various cellular tiers following ICH, aiming to acquire biological insights into ICH. We utilized AUCell, Ucell, singscore, ssgsea, and AddModuleScore algorithms, along with correlation analysis, to identify hub genes influencing high OS post-ICH. Furthermore, we employed four machine learning algorithms, eXtreme Gradient Boosting, Boruta, Random Forest, and Least Absolute Shrinkage and Selection Operator, to identify the optimal feature genes. To validate the accuracy of our analysis, we conducted validation in ICH animal experiments., Results: After analyzing the scRNA-seq dataset using various algorithms, we found that OS activity exhibited heterogeneity across different cellular layers following ICH, with particularly heightened activity observed in monocytes. Further integration of bulk data and machine learning algorithms revealed that ANXA2 and COTL1 were closely associated with high OS after ICH. Our animal experiments demonstrated an increase in OS expression post-ICH. Additionally, the protein expression of ANXA2 and COTL1 was significantly elevated and co-localized with microglia. Pearson correlation coefficient analysis revealed a significant correlation between ANXA2 and OS, indicating strong consistency (r = 0.84, p < 0.05). Similar results were observed for COTL1 and OS (r = 0.69, p < 0.05)., Conclusions: Following ICH, ANXA2 and COTL1 might penetrate the brain via monocytes, localize within microglia, and enhance OS activity. This might help us better understand OS after ICH., 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: Cervical cancer (CC) is the fourth most common malignancy among women globally and serves as the main cause of cancer-related deaths among women in developing countries. The early symptoms of CC are often not apparent, with diagnoses typically made at advanced stages, which lead to poor clinical prognoses. In recent years, numerous studies have shown that there is a close relationship between mast cells (MCs) and tumor development. However, research on the role MCs played in CC is still very limited at that time. Thus, the study conducted a single-cell multi-omics analysis on human CC cells, aiming to explore the mechanisms by which MCs interact with the tumor microenvironment in CC. The goal was to provide a scientific basis for the prevention, diagnosis, and treatment of CC, with the hope of improving patients' prognoses and quality of life., Method: The present study acquired single-cell RNA sequencing data from ten CC tumor samples in the ArrayExpress database. Slingshot and AUCcell were utilized to infer and assess the differentiation trajectory and cell plasticity of MCs subpopulations. Differential expression analysis of MCs subpopulations in CC was performed, employing Gene Ontology, gene set enrichment analysis, and gene set variation analysis. CellChat software package was applied to predict cell communication between MCs subpopulations and CC cells. Cellular functional experiments validated the functionality of TNFRSF12A in HeLa and Caski cell lines. Additionally, a risk scoring model was constructed to evaluate the differences in clinical features, prognosis, immune infiltration, immune checkpoint, and functional enrichment across various risk scores. Copy number variation levels were computed using inference of copy number variations., Result: The obtained 93,524 high-quality cells were classified into ten cell types, including T&#95;NK cells, endothelial cells, fibroblasts, smooth muscle cells, epithelial cells, B cells, plasma cells, MCs, neutrophils, and myeloid cells. Furthermore, a total of 1,392 MCs were subdivided into seven subpopulations: C0 CTSG+ MCs, C1 CALR+ MCs, C2 ALOX5+ MCs, C3 ANXA2+ MCs, C4 MGP+ MCs, C5 IL32+ MCs, and C6 ADGRL4+ MCs. Notably, the C2 subpopulation showed close associations with tumor-related MCs, with Slingshot results indicating that C2 subpopulation resided at the intermediate-to-late stage of differentiation, potentially representing a crucial transition point in the benign-to-malignant transformation of CC. CNVscore and bulk analysis results further confirmed the transforming state of the C2 subpopulation. CellChat analysis revealed TNFRSF12A as a key receptor involved in the actions of C2 ALOX5+ MCs. Moreover, in vitro experiments indicated that downregulating the TNFRSF12A gene may partially inhibit the development of CC. Additionally, a prognosis model and immune infiltration analysis based on the marker genes of the C2 subpopulation provided valuable guidance for patient prognosis and clinical intervention strategies., Conclusions: We first identified the transformative tumor-associated MCs subpopulation C2 ALOX5+ MCs within CC, which was at a critical stage of tumor differentiation and impacted the progression of CC. In vitro experiments confirmed the inhibitory effect of knocking down the TNFRSF12A gene on the development of CC. The prognostic model constructed based on the C2 ALOX5+MCs subset demonstrated excellent predictive value. These findings offer a fresh perspective for clinical decision-making in CC., 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 Zhao, Hong, Zhou, Huang, Lin, Zhang, Liang and Tang.)

Main Conclusion: Transcriptome analysis in potato varieties revealed genes associated with tuber yield-related traits and developed gene expression markers. This study aimed to identify genes involved in high tuber yield and its component traits in test potato varieties (Kufri Frysona, Kufri Khyati, and Kufri Mohan) compared to control (Kufri Sutlej). The aeroponic evaluation showed significant differences in yield-related traits in the varieties. Total RNA sequencing was performed using tuber and leaf tissues on the Illumina platform. The high-quality reads (QV > 25) mapping with the reference potato genomes revealed statistically significant (P < 0.05) differentially expressed genes (DEGs) into two categories: up-regulated (> 2 Log 2 fold change) and down-regulated (< -2 Log 2 fold change). DEGs were characterized by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Collectively, we identified genes participating in sugar metabolism, stress response, transcription factors, phytohormones, kinase proteins, and other genes greatly affecting tuber yield and its related traits. A few selected genes were UDP-glucose glucosyltransferase, glutathion S-transferase, GDSL esterase/lipase, transcription factors (MYB, WRKY, bHLH63, and BURP), phytohormones (auxin-induced protein X10A, and GA20 oxidase), kinase proteins (Kunitz-type tuber invertase inhibitor, BRASSINOSTEROID INSENSITIVE 1-associated receptor kinase 1) and laccase. Based on the selected 17 peptide sequences representing 13 genes, a phylogeny tree and motifs were analyzed. Real time-quantitative polymerase chain reaction (RT-qPCR) analysis was used to validate the RNA-seq results. RT-qPCR based gene expression markers were developed for the genes such as 101 kDa heat shock protein, catechol oxidase B chloroplastic, cysteine protease inhibitor 1, Kunitz-type tuber invertase inhibitor, and laccase to identify high yielding potato genotypes. Thus, our study paved the path for potential genes associated with tuber yield traits in potato under aeroponics., (© 2024. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Next-generation sequencing (NGS) has revolutionized genomic research by enabling high-throughput, cost-effective genome and transcriptome sequencing accelerating personalized medicine for complex diseases, including cancer. Whole genome/transcriptome sequencing (WGS/WTS) provides comprehensive insights, while targeted sequencing is more cost-effective and sensitive. In comparison to short-read sequencing, which still dominates the field due to high speed and cost-effectiveness, long-read sequencing can overcome alignment limitations and better discriminate similar sequences from alternative transcripts or repetitive regions. Hybrid sequencing combines the best strengths of different technologies for a more comprehensive view of genomic/transcriptomic variations. Understanding each technology's strengths and limitations is critical for translating cutting-edge technologies into clinical applications. In this study, we sequenced DNA and RNA libraries of reference samples using various targeted DNA and RNA panels and the whole transcriptome on both short-read and long-read platforms. This study design enables a comprehensive analysis of sequencing technologies, targeting protocols, and library preparation methods. Our expanded profiling landscape establishes a reference point for assessing current sequencing technologies, facilitating informed decision-making in genomic research and precision medicine., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Background: Reproducibility is a major concern in biomedical studies, and existing publication guidelines do not solve the problem. Batch effects and quality imbalances between groups of biological samples are major factors hampering reproducibility. Yet, the latter is rarely considered in the scientific literature., Results: Our analysis uses 40 clinically relevant RNA-seq datasets to quantify the impact of quality imbalance between groups of samples on the reproducibility of gene expression studies. High-quality imbalance is frequent (14 datasets; 35%), and hundreds of quality markers are present in more than 50% of the datasets. Enrichment analysis suggests common stress-driven effects among the low-quality samples and highlights a complementary role of transcription factors and miRNAs to regulate stress response. Preliminary ChIP-seq results show similar trends. Quality imbalance has an impact on the number of differential genes derived by comparing control to disease samples (the higher the imbalance, the higher the number of genes), on the proportion of quality markers in top differential genes (the higher the imbalance, the higher the proportion; up to 22%) and on the proportion of known disease genes in top differential genes (the higher the imbalance, the lower the proportion). We show that removing outliers based on their quality score improves the resulting downstream analysis., Conclusions: Thanks to a stringent selection of well-designed datasets, we demonstrate that quality imbalance between groups of samples can significantly reduce the relevance of differential genes, consequently reducing reproducibility between studies. Appropriate experimental design and analysis methods can substantially reduce the problem., (© 2024. The Author(s).)

Objectives: The present study used single-cell RNA sequencing (scRNA-seq) to characterize the cellular composition of ovarian carcinosarcoma (OCS) and identify its molecular characteristics., Methods: scRNA-seq was performed in resected primary OCS for an in-depth analysis of tumor cells and the tumor microenvironment. Immunohistochemistry staining was used for validation. The scRNA-seq data of OCS were compared with those of high-grade serous ovarian carcinoma (HGSOC) tumors and other OCS tumors., Results: Both malignant epithelial and malignant mesenchymal cells were observed in the OCS patient of this study. We identified four epithelial cell subclusters with different biological roles. Among them, epithelial subcluster 4 presented high levels of breast cancer type 1 susceptibility protein homolog ( BRCA1 ) and DNA topoisomerase 2-α ( TOP2A ) expression and was related to drug resistance and cell cycle. We analyzed the interaction between epithelial and mesenchymal cells and found that fibroblast growth factor (FGF) and pleiotrophin (PTN) signalings were the main pathways contributing to communication between these cells. Moreover, we compared the malignant epithelial and mesenchymal cells of this OCS tumor with our previous published HGSOC scRNA-seq data and OCS data. All the epithelial subclusters in the OCS tumor could be found in the HGSOC samples. Notably, the mesenchymal subcluster C14 exhibited specific expression patterns in the OCS tumor, characterized by elevated expression of cytochrome P450 family 24 subfamily A member 1 ( CYP24A1 ), collagen type XXIII α1 chain ( COL23A1 ), cholecystokinin ( CCK ), bone morphogenetic protein 7 ( BMP7 ), PTN , Wnt inhibitory factor 1 ( WIF1 ), and insulin-like growth factor 2 ( IGF2 ). Moreover, this subcluster showed distinct characteristics when compared with both another previously published OCS tumor and normal ovarian tissue., Conclusions: This study provides the single-cell transcriptomics signature of human OCS, which constitutes a new resource for elucidating OCS diversity.

Rainbow trout ( Oncorhynchus mykiss , Walbaum , 1792 ) is an important economic cold-water fish that is susceptible to heat stress. To date, the heat stress response in rainbow trout is more widely understood at the transcriptional level, while little research has been conducted at the translational level. To reveal the translational regulation of heat stress in rainbow trout, in this study, we performed a ribosome profiling assay of rainbow trout liver under normal and heat stress conditions. Comparative analysis of the RNA-seq data with the ribosome profiling data showed that the folding changes in gene expression at the transcriptional level are moderately correlated with those at the translational level. In total, 1213 genes were significantly altered at the translational level. However, only 32.8% of the genes were common between both levels, demonstrating that heat stress is coordinated across both transcriptional and translational levels. Moreover, 809 genes exhibited significant differences in translational efficiency (TE), with the TE of these genes being considerably affected by factors such as the GC content, coding sequence length, and upstream open reading frame (uORF) presence. In addition, 3468 potential uORFs in 2676 genes were identified, which can potentially affect the TE of the main open reading frames. In this study, Ribo-seq and RNA-seq were used for the first time to elucidate the coordinated regulation of transcription and translation in rainbow trout under heat stress. These findings are expected to contribute novel data and theoretical insights to the international literature on the thermal stress response in fish.

Hematopoiesis is a tightly regulated process that gets skewed toward myelopoiesis. This restrains lymphopoiesis, but the role of lymphocytes in this process is not well defined. To unravel the intricacies of neutrophil responses in COVID-19, we performed bulk RNAseq on neutrophils from healthy controls and COVID-19 patients. Principal component analysis revealed distinguishing neutrophil gene expression alterations in COVID-19 patients. ICU and ward patients displayed substantial transcriptional changes, with ICU patients exhibiting a more pronounced response. Intriguingly, neutrophils from COVID-19 patients, notably ICU patients, exhibited an enrichment of immunoglobulin (Ig) and B cell lineage-associated genes, suggesting potential lineage plasticity. We validated our RNAseq findings in a larger cohort. Moreover, by reanalyzing single-cell RNA sequencing (scRNAseq) data on human bone marrow (BM) granulocytes, we identified the cluster of granulocyte-monocyte progenitors (GMP) enriched with Ig and B cell lineage-associated genes. These cells with lineage plasticity may serve as a resource depending on the host's needs during severe systemic infection. This distinct B cell subset may play a pivotal role in promoting myelopoiesis in response to infection. The scRNAseq analysis of BM neutrophils in infected mice further supported our observations in humans. Finally, our studies using an animal model of acute infection implicate IL-7/GM-CSF in influencing neutrophil and B cell dynamics. Elevated GM-CSF and reduced IL-7 receptor expression in COVID-19 patients imply altered hematopoiesis favoring myeloid cells over B cells. Our findings provide novel insights into the relationship between the B-neutrophil lineages during severe infection, hinting at potential implications for disease pathogenesis., Importance: This study investigates the dynamics of hematopoiesis in COVID-19, focusing on neutrophil responses. Through RNA sequencing of neutrophils from healthy controls and COVID-19 patients, distinct gene expression alterations are identified, particularly in ICU patients. Notably, neutrophils from COVID-19 patients, especially in the ICU, exhibit enrichment of immunoglobulin and B cell lineage-associated genes, suggesting potential lineage plasticity. Validation in a larger patient cohort and single-cell analysis of bone marrow granulocytes support the presence of granulocyte-monocyte progenitors with B cell lineage-associated genes. The findings propose a link between B-neutrophil lineages during severe infection, implicating a potential role for these cells in altered hematopoiesis favoring myeloid cells over B cells. Elevated GM-CSF and reduced IL-7 receptor expression in stress hematopoiesis suggest cytokine involvement in these dynamics, providing novel insights into disease pathogenesis., Competing Interests: The authors declare no conflict of interest.

Inferring gene regulatory networks from single-cell RNA-sequencing trajectories has been an active area of research yet methods are still needed to identify regulators governing cell transitions. We developed DREAMIT (Dynamic Regulation of Expression Across Modules in Inferred Trajectories) to annotate transcription-factor activity along single-cell trajectory branches, using ensembles of relations to target genes. Using a benchmark representing several different tissues, as well as external validation with ATAC-Seq and Perturb-Seq data on hematopoietic cells, the method was found to have higher tissue-specific sensitivity and specificity over competing approaches., (© 2024. The Author(s).)

Rationale : Immunosuppressive tumor microenvironment (iTME) plays an important role in carcinogenesis, and some macrophage subsets are associated with iTME generation. However, the sub-population characterization of macrophages in oral carcinogenesis remains largely unclear. Here, we investigated the immunosuppressive status with focus on function of a macrophage subset that expressed indoleamine 2,3 dioxygenase 1 (Macro-IDO1) in oral carcinogenesis. Methods : We built a single cell transcriptome atlas from 3 patients simultaneously containing oral squamous cell carcinoma (OSCC), precancerous oral leukoplakia (preca-OLK) and paracancerous tissue (PCA). Through single-cell RNA sequencing and further validation using multicolor immunofluorescence staining and the in vitro / in vivo experiments, the immunosuppressive cell profiles were built and the role of a macrophage subset that expressed indoleamine 2,3 dioxygenase 1 (Macro-IDO1) in the malignant transformation of oral leukoplakia was evaluated. Results : The iTME formed at preca-OLK stage, as evidenced by increased exhausted T cells, Tregs and some special subsets of macrophages and fibroblasts. Macro-IDO1 was predominantly enriched in preca-OLK and OSCC, distributed near exhausted T cells and possessed tumor associated macrophage transformation potentials. Functional analysis revealed the established immunosuppressive role of Macro-IDO1 in preca-OLK and OSCC: enriching the immunosuppression related genes; having an established level of immune checkpoint score; exerting strong immunosuppressive interaction with T cells; positively correlating with the CD8-exhausted. The immunosuppression related gene expression of macrophages also increased in preca-OLK/OSCC compared to PCA. The use of the IDO1 inhibitor reduced 4NQO induced oral carcinogenesis in mice. Mechanistically, IFN-γ-JAK-STAT pathway was associated with IDO1 upregulation in OLK and OSCC. Conclusions : These results highlight that Macro-IDO1-enriched in preca-OLK possesses a strong immunosuppressive role and contributes to oral carcinogenesis, providing a potential target for preventing precancerous legions from transformation into OSCC., Competing Interests: Competing Interests: The authors have declared that no competing interest exists., (© The author(s).)

Peanut (Arachis hypogaea L.) is an important leguminous oil and economic crop that produces flowers aboveground and fruits belowground. Subterranean fruit-pod development, which significantly affects peanut production, involves complex molecular mechanisms that likely require the coordinated regulation of multiple genes in different tissues. To investigate the molecular mechanisms that underlie peanut fruit-pod development, we characterized the anatomical features of early fruit-pod development and integrated single-nucleus RNA-sequencing (snRNA-seq) and single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) data at the single-cell level. We identified distinct cell types, such as meristem, embryo, vascular tissue, cuticular layer, and stele cells within the shell wall. These specific cell types were used to examine potential molecular changes unique to each cell type during pivotal stages of fruit-pod development. snRNA-seq analyses of differentially expressed genes revealed cell-type-specific insights that were not previously obtainable from transcriptome analyses of bulk RNA. For instance, we identified MADS-box genes that contributes to the formation of parenchyma cells and gravity-related genes that are present in the vascular cells, indicating an essential role for the vascular cells in peg gravitropism. Overall, our single-nucleus analysis provides comprehensive and novel information on specific cell types, gene expression, and chromatin accessibility during the early stages of fruit-pod development. This information will enhance our understanding of the mechanisms that underlie fruit-pod development in peanut and contribute to efforts aimed at improving peanut production., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Pancreatic cancer (PC) poses a substantial global health challenge, ranking as the fourth leading cause of cancer-related deaths due to its high mortality rate. Late-stage diagnoses are common due to the absence of specific symptoms. Pancreatic ductal adenocarcinoma (PDAC) accounts for the majority of PC cases. Recent research has suggested a potential link between elevated serum levels of bile acids (BAs) and tumorigenesis of PDAC. This study aims to understand how taurochenodeoxycholic acid (TCDCA), a secondary BA, influences PDAC using RNA sequencing techniques on the Capan-1 cell line. We identified 2,950 differentially expressed genes (DEGs) following TCDCA treatment, with 1,597 upregulated and 1,353 downregulated genes. These DEGs were associated with critical PDAC pathways, including coagulation, angiogenesis, cell migration, and signaling regulation. Furthermore, we reviewed relevant literature highlighting genes like DKK-1, KRT80, UPLA, and SerpinB2, known for their roles in PDAC tumorigenesis and metastasis. Our study sheds light on the complex relationship between BAs and PDAC, offering insights into potential diagnostic markers and therapeutic targets. Further research is needed to unravel these findings' precise mechanisms and clinical implications, potentially improving PDAC diagnosis and treatment., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Zoltan Vereb reports financial support was provided by Hungarian Academy of Sciences. Viktoria Venglovecz reports financial support was provided by National, Research, Development and Innovation Office. Lajos Kemeny reports financial support was provided by Hungarian Centre of Excellence for Molecular Medicine. Zoltan Vereb reports financial support was provided by National, Research, Development and Innovation Office. If there are other authors, they 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 The Authors. Published by Elsevier B.V. All rights reserved.)

During maturation oocytes at the germinal vesicle (GV) stage progress to metaphase II (MII). However, during in vitro maturation a proportion often fail to progress. To understand these processes, we employed RNA sequencing to examine the transcriptome profile of these three groups of oocytes from the pig. We compared our findings with similar public oocyte data from humans. The transcriptomes in oocytes that failed to progress was similar to those that did. We found in both species, the most upregulated genes in MII oocytes were associated with chromosome segregation and cell cycle processes, while the most down regulated genes were relevant to ribosomal and mitochondrial pathways. Moreover, those genes involved in chromosome segregation during GV to MII transition were conserved in pig and human. We also compared MII and GV oocyte transcriptomes at the isoform transcript level in both species. Several thousands of genes (including DTNBP1, MAPK1, RAB35, GOLGA7, ATP1A1 and ATP2B1) identified as not different in expression at a gene transcript level were found to have differences in isoform transcript levels. Many of these genes were involved in ATPase-dependent or GTPase-dependent intracellular transport in pig and human, respectively. In conclusion, our study suggests the failure to progress to MII in vitro may not be regulated at the level of the genome and that many genes are differentially regulated at the isoform level, particular those involved ATPase- or GTPase-dependent intracellular transport., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Tang 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.)

Background: Xenotransplantation of genetically engineered porcine organs has the potential to address the challenge of organ donor shortage. Two cases of porcine-to-human kidney xenotransplantation were performed, yet the physiological effects on the xenografts and the recipients' immune responses remain largely uncharacterized., Methods: We performed single-cell RNA sequencing (scRNA-seq) and longitudinal RNA-seq analyses of the porcine kidneys to dissect xenotransplantation-associated cellular dynamics and xenograft-recipient interactions. We additionally performed longitudinal scRNA-seq of the peripheral blood mononuclear cells (PBMCs) to detect recipient immune responses across time., Findings: Although no hyperacute rejection signals were detected, scRNA-seq analyses of the xenografts found evidence of endothelial cell and immune response activation, indicating early signs of antibody-mediated rejection. Tracing the cells' species origin, we found human immune cell infiltration in both xenografts. Human transcripts in the longitudinal bulk RNA-seq revealed that human immune cell infiltration and the activation of interferon-gamma-induced chemokine expression occurred by 12 and 48 h post-xenotransplantation, respectively. Concordantly, longitudinal scRNA-seq of PBMCs also revealed two phases of the recipients' immune responses at 12 and 48-53 h. Lastly, we observed global expression signatures of xenotransplantation-associated kidney tissue damage in the xenografts. Surprisingly, we detected a rapid increase of proliferative cells in both xenografts, indicating the activation of the porcine tissue repair program., Conclusions: Longitudinal and single-cell transcriptomic analyses of porcine kidneys and the recipient's PBMCs revealed time-resolved cellular dynamics of xenograft-recipient interactions during xenotransplantation. These cues can be leveraged for designing gene edits and immunosuppression regimens to optimize xenotransplantation outcomes., Funding: This work was supported by NIH RM1HG009491 and DP5OD033430., Competing Interests: Declaration of interests J.D.B. is a founder and director of CDI Labs, Inc.; a founder of and consultant to Opentrons LabWorks/Neochromosome, Inc.; and serves or served on the scientific advisory boards of the following: CZ Biohub New York, LLC; Logomix, Inc.; Modern Meadow, Inc.; Rome Therapeutics, Inc.; Sangamo, Inc.; Tessera Therapeutics, Inc.; and the Wyss Institute. R.A.M. is on scientific advisory boards for eGenesis, Sanofi, Regeneron, CareDx, and Hansa Biopharma; is a consultant to Recombinetics; reports consulting fees from Hansa Medical, Regeneron, Thermo Fisher Scientific, Genentech, CareDx, One Lambda, ITB Med, Sanofi, and PPD Development; and reports grant support from Hansa Biopharma, all unrelated to the present work. R.A.M. also reports grant support from United Therapeutics Corporation, PBC. All other authors have no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Gene misexpression is the aberrant transcription of a gene in a context where it is usually inactive. Despite its known pathological consequences in specific rare diseases, we have a limited understanding of its wider prevalence and mechanisms in humans. To address this, we analyzed gene misexpression in 4,568 whole-blood bulk RNA sequencing samples from INTERVAL study blood donors. We found that while individual misexpression events occur rarely, in aggregate they were found in almost all samples and a third of inactive protein-coding genes. Using 2,821 paired whole-genome and RNA sequencing samples, we identified that misexpression events are enriched in cis for rare structural variants. We established putative mechanisms through which a subset of SVs lead to gene misexpression, including transcriptional readthrough, transcript fusions, and gene inversion. Overall, we develop misexpression as a type of transcriptomic outlier analysis and extend our understanding of the variety of mechanisms by which genetic variants can influence gene expression., Competing Interests: Declaration of interests The authors declare the following interests: T.V. has received PhD studentship funding from AstraZeneca. J.M. completed this work while employed by the University of Cambridge but is now an employee of Genomics plc. S.P. is a current employee and stockholder of AstraZeneca. D.J.R. is an employee of NHS Blood and Transplant. A.B. is currently an employee of Bayer AG, Research and Early Development Precision Medicine, Research & Development, Pharmaceutical Division, Wuppertal, DE. A.P. is a current employee and stockholder of AstraZeneca. D.S.P. is a current employee and stockholder of AstraZeneca. K.K. is a current employee and stockholder of AstraZeneca. J.D. serves on scientific advisory boards for AstraZeneca, Novartis, and UK Biobank and has received multiple grants from academic, charitable, and industry sources outside of the submitted work. M.I. is a trustee of the Public Health Genomics (PHG) Foundation, is a member of the Scientific Advisory Board of Open Targets, and has a research collaboration with AstraZeneca, which is unrelated to this study., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Small non-coding RNAs (sncRNAs) derived from tRNAs are known as tRNA-derived small RNAs (tsRNAs). These tsRNAs are further categorized into tRNA-derived fragments (tRFs) and tRNA halves (tiRNAs), which play significant roles in the various molecular mechanisms underlying certain human diseases. However, the generation of tsRNAs and their potential roles during Dengue virus (DENV) infection is not yet known. Here, we performed small RNA sequencing to identify the generation and alterations in tsRNAs expression profiles of DENV-infected Huh7 cells. Upon DENV infection, tRNA fragmentation was found to be increased. We identified a significant number of differentially expressed tsRNAs during DENV infection. Interestingly, the 3'tRF population showed upregulation, while the i-tRF population exhibited downregulation. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed to analyze the impact of differentially expressed tsRNAs on DENV pathogenesis. Our results suggest that differentially expressed tsRNAs are involved in transcriptional regulation via RNA polymerase II promoter and metabolic pathways. Overall, our study contributes significantly to our understanding of the roles played by tsRNAs in the complex dynamics of DENV infection., (© 2024. The Author(s).)

Introduction: Gastric cancer (GC) remains a major global health threat ranking as the fifth most prevalent cancer. Hypoxia, a characteristic feature of solid tumors, significantly contributes to the malignant progression of GC. Mitochondria are the major target of hypoxic injury that promotes mitochondrial dysfunction during the development of cancers including GC. However, the gene signature and prognostic model based on hypoxia- and mitochondrial dysfunction-related genes (HMDRGs) in the prediction of GC prognosis have not yet been established., Methods: The gene expression profile datasets of stomach cancer patients were retrieved from The Cancer Genome Atlas and the Gene Expression Omnibus databases. Prognostic genes were selected using Least Absolute Shrinkage and Selection Operator Cox (LASSO-Cox) regression analysis to construct a prognostic model. Immune infiltration was evaluated through ESTIMATE, CIBERSORT, and ssGSEA analyses. Tumor immune dysfunction and exclusion (TIDE) and immunophenoscore (IPS) were utilized to explore implications for immunotherapy. Furthermore, in vitro experiments were conducted to validate the functional roles of HMDRGs in GC cell malignancy., Results: In this study, five HMDRGs (ZFP36, SERPINE1, DUSP1, CAV1, and AKAP12) were identified for developing a prognostic model in GC. This model stratifies GC patients into high- and low-risk groups based on median risk scores. A nomogram predicting overall survival (OS) was constructed and showed consistent results with observed OS. Immune infiltration analysis indicated that individuals in the high-risk group tend to exhibit increased immune cell infiltration. Additionally, analysis of cancer immunotherapy responses revealed that high-risk group patients exhibit poorer responses to cancer immunotherapy compared to the low-risk group. Immunohistochemistry (IHC) staining indicated that the expression levels of HMDRGs were remarkably correlated with GC, of which, SERPINE1 displayed the most pronounced up-regulation, while ZFP36 exhibited the most notable down-regulation in GC patients. Furthermore, in vitro investigation validated that SERPINE1 and ZFP36 contribute to the malignant processes of GC cells correlated with mitochondrial dysfunction., Conclusions: This study presents a novel and efficient approach to evaluate GC prognosis and immunotherapy efficacy, and also provides insights into understanding the pathogenesis of GC., 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, Cui, Wang, Wang, Yu and Xiong.)

Background: Hemorrhagic enteritis, caused by Turkey Hemorrhagic Enteritis Virus (THEV), is a disease affecting turkey poults characterized by immunosuppression and bloody diarrhea. An avirulent THEV strain that retains the immunosuppressive ability is used as a live vaccine. Characterizing the splice map of THEV is an essential step that would allow studies of individual genes mediating its immunosuppressive functions. We used RNA sequencing to characterize the splice map of THEV for the first time, providing key insights into the THEV gene expression and mRNA structures., Methods: After infecting a turkey B-cell line with the vaccine strain, samples in triplicates were collected at 4-, 12-, 24-, and 72-hours post-infection. Total RNA was extracted, and poly-A-tailed mRNA sequenced. Reads were mapped to the THEV genome after trimming and transcripts assembled with StringTie. We performed PCR of THEV cDNA, cloned the PCR products, and used Sanger sequencing to validate all identified splice junctions., Results: Researchers previously annotated the THEV genome as encoding 23 open reading frames (ORFs). We identified 29 spliced transcripts from our RNA sequencing data, all containing novel exons although some exons matched some previously annotated ORFs. The three annotated splice junctions were also corroborated by our data. During validation we identified five additional unique transcripts, a subset of which were further validated by 3' rapid amplification of cDNA ends (3' RACE). Thus, we report that the genome of THEV contains 34 transcripts with the coding capacity for all annotated ORFs. However, we found six of the previously annotated ORFs to be truncated ORFs on the basis of the identification of an in-frame upstream start codon or the detection of additional coding exons. We also identified three of the annotated ORFs with longer or shorter isoforms, and seven novel unannotated ORFs that could potentially be translated; although it is beyond the scope of this manuscript to investigate whether they are translated., Conclusions: Similar to human adenoviruses, all THEV transcripts are spliced and organized into five transcription units under the control of their cognate promoters. The genes are expressed under temporal regulation and THEV also produces multiple distinctly spliced transcripts that code for the same protein. Studies of the newly identified potential proteins should be urgently performed as these proteins may have roles in THEV-induced immunosuppression. Also, knowing the splicing of THEV genes should be invaluable to future research focusing on studying THEV genes, as this will allow accurate cloning of the mRNAs., (© 2024. The Author(s).)

Multimodal measurement of single cells provides deep insights into the intricate relationships between individual molecular layers and the regulatory mechanisms underlying intercellular variations. Here, we reported DMF-DM-seq, a highly integrated, sensitive, and automated platform for single-cell mRNA and microRNA (miRNA) co-sequencing based on digital microfluidics. This platform first integrates the processes of single-cell isolation, lysis, component separation, and simultaneous sequencing library preparation of mRNA and miRNA within a single DMF device. Compared with the current half-cell measuring strategy, DMF-DM-seq enables complete separation of single-cell mRNA and miRNA via a magnetic field application, resulting in a higher miRNA detection ability. DMF-DM-seq revealed differential expression patterns of single cells of noncancerous breast cells and noninvasive and aggressive breast cancer cells at both mRNA and miRNA levels. The results demonstrated the anticorrelated relationship between miRNA and their mRNA targets. Further, we unravel the tumor growth and metastasis-associated biological processes enriched by miRNA-targeted genes, along with important miRNA-interaction networks involved in significant signaling pathways. We also deconstruct the miRNA regulatory mechanisms underlying different signaling pathways across different breast cell types. In summary, DMF-DM-seq offers a powerful tool for a comprehensive study of the expression heterogeneity of single-cell mRNA and miRNA, which will be widely applied in basic and clinical research.

Background: Extensive research has been conducted on embryonic developmental disorders linked to Polycystic Ovary Syndrome (PCOS), a pathological condition that affects 5-10% of women and is characterized by irregularities in the menstrual cycle and infertility. By employing RNA sequencing (RNA-seq), we performed an in-depth investigation of PCOS-related changes in gene expression patterns at the mouse blastocyst stage., Methods: The zygotes of female B6D2 mice were obtained and then differentiated into blastocysts in K + Simplex Optimised Medium (KSOM) cultures containing exo-NC (negative control for exosomes) or exo-LIPE-AS1 (a novel exosomal marker of PCOS). Subsequently, blastocysts were collected for RNA-seq. The bioinformatics was performed to analyze and compare the differences of gene expression profile between blastocysts of control and PCOS group., Results: There were 1150 differentially expressed genes (DEGs) between the two groups of mouse blastocysts; 243 genes were upregulated and 907 downregulated in the blastocysts of the exo-LIPE-AS1 group compared to those of the exo-NC group. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the genes involved in amino acid synthesis and glutathione metabolic pathways were down-regulated in exo-LIPE-AS1 group., Conclusion: This study has revealed that blastocyst developmental retardation may be associated with the downregulation of amino acid synthesis and glutathione metabolism, which may affect energy metabolism, biosynthesis, cellular osmotic pressure, antioxidant synthesis, ROS clearance or mitochondrial function, and ultimately cause blastocyst cell development abnormalities. Our research offers encouraging data on the mechanisms underlying aberrant embryonic development in patients with PCOS as well as potential treatment strategies., (© 2024. The Author(s).)

This study investigates the cellular origin and tissue heterogeneity in bipolar disorder (BD) by integrating multiomics data. Four distinct datasets were employed, including single-cell RNA sequencing (scRNA-seq) data (embryonic and fetal brain, n = 8, 1,266 cells), BD Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data (adult brain, n = 210), BD bulk RNA-seq data (adult brain, n = 314), and BD genome-wide association study (GWAS) summary data (n = 413,466). The integration of scRNA-seq data with multiomics data relevant to BD was accomplished using the single-cell disease relevance score (scDRS) algorithm. We have identified a novel brain cell cluster named ADCY1, which exhibits distinct genetic characteristics. From a high-resolution genetic perspective, glial cells emerge as the primary cytopathology associated with BD. Specifically, astrocytes were significantly related to BD at the RNA-seq level, while microglia showed a strong association with BD across multiple panels, including the transcriptome-wide association study (TWAS), ATAC-seq, and RNA-seq. Additionally, oligodendrocyte precursor cells displayed a significant association with BD in both ATAC-seq and RNA-seq panel. Notably, our investigation of brain regions affected by BD revealed significant associations between BD and all three types of glial cells in the dorsolateral prefrontal cortex (DLPFC). Through comprehensive analyses, we identified several BD-associated genes, including CRMP1, SYT4, UCHL1, and ZBTB18. In conclusion, our findings suggest that glial cells, particularly in specific brain regions such as the DLPFC, may play a significant role in the pathogenesis of BD. The integration of multiomics data has provided valuable insights into the etiology of BD, shedding light on potential mechanisms underlying this complex psychiatric disorder., (© 2024. The Author(s).)

Platelet-rich fibrin (PRF) is prepared by spontaneous coagulation of fractionated blood. When squeezed between two plates, PRF is separated into solid PRF membranes and a liquid exudate, the PRF serum. The question arises regarding how much the overall activity remains in the PRF membranes and what is discarded into the PRF serum. To this end, we have exposed gingival fibroblasts to lysates prepared from PRF membranes and PRF serum, followed by bulk RNA sequencing. A total of 268 up- and 136 down-regulated genes in gingival fibroblasts exposed to PRF membrane lysates were significantly regulated under the premise of a minimum log2 with 2.5-fold change and a minus log10 significance level of two, respectively. PRF serum only caused 62 up- and 32 down-regulated genes under these conditions. Among the 46 commonly up-regulated genes were CXCL1, CXCL5, CXCL6, CXCL8, IL33, IL6, and PTGS2/COX2, stanniocalcin-1-all linked to an inflammatory response. PRF membrane lysates further increased chemokines CCL2, CCL7, CXCL2, CXCL3, and IL1R1, IL1RL1, and IL1RN, as well as the paracrine factors IL11, LIF, IGF1, BMP2, BMP6, FGF2, and CCN2/CTGF, and all hyaluronan synthases. On the other hand, PRF serum increased DKK1. The genes commonly down-regulated by PRF membrane lysates and PRF serum included interferon-induced protein with tetratricopeptide repeats (IFIT1, IFIT2, IFIT3) and odd-skipped-related transcription factors (OSR1 and OSR2), as well as FGF18 and GDF15, respectively. Taken together, PRF membrane lysates, compared to PRF serum, cause a more complex response in gingival fibroblasts, but each increased chemokine expression in gingival fibroblasts.

The role of monocytes in postmenopausal osteoporosis is widely recognized; however, the mechanisms underlying monocyte reprogramming remain unknown. In this study, single-cell RNA sequencing (scRNA-seq) was conducted on CD14+ bone marrow monocytes obtained from 3 postmenopausal women with normal BMD and 3 women with postmenopausal osteoporosis (PMOP). Monocle2 was used to classify the monocytes into 7 distinct clusters. The proportion of cluster 1 significantly decreased in PMOP patients, while the proportion of cluster 7 increased. Further analysis via GSEA, transcription factor activity analysis, and sc-metabolic analysis revealed significant differences between clusters 1 and 7. Cluster 7 exhibited upregulated pathways associated with inflammation, immunity, and osteoclast differentiation, whereas cluster 1 demonstrated the opposite results. Monocle2, TSCAN, VECTOR, and scVelo data indicated that cluster 1 represented the initial subset and that cluster 7 represents one of the terminal subsets. BayesPrism and ssGSEA were employed to analyze the bulk transcriptome data obtained from the GEO database. The observed alterations in the proportions of 1 and 7 were validated and found to have diagnostic significance. CD16 serves as the marker gene for cluster 7, thus leading to an increased proportion of CD16+ monocytes in women with PMOP. Flow cytometry was used to assess the consistency of outcomes with those of the bioinformatic analysis. Subsequently, an additional scRNA-seq analysis was conducted on bone marrow mononuclear cells obtained from 3 patients with PMOP and 3 postmenopausal women with normal BMD. The differential proportions of cluster 1 and cluster 7 were once again confirmed, with the pathological effect of cluster 7 may attribute to cell-cell communication. The scRNA-seq findings suggest that an imbalance in monocyte subsets is a characteristic feature of PMOP. These findings elucidate the limitations of utilizing bulk transcriptome data for detecting alterations in monocytes, which may influence novel research inquiries., (© The Author(s) 2024. Published by Oxford University Press on behalf of the American Society for Bone and Mineral Research.)

The placenta is crucial for fetal development, yet the impact of environmental stressors such as arsenic exposure remains poorly understood. We apply single-cell RNA sequencing to analyze the response of the mouse placenta to arsenic, revealing cell-type-specific gene expression, function, and pathological changes. Notably, the Prap1 gene, which encodes proline-rich acidic protein 1 (PRAP1), is significantly upregulated in 26 placental cell types including various trophoblast cells. Our study shows a female-biased increase in PRAP1 in response to arsenic and localizes it in the placenta. In vitro and ex vivo experiments confirm PRAP1 upregulation following arsenic treatment and demonstrate that recombinant PRAP1 protein reduces arsenic-induced cytotoxicity and downregulates cell cycle pathways in human trophoblast cells. Moreover, PRAP1 knockdown differentially affects cell cycle processes, proliferation, and cell death depending on the presence of arsenic. Our findings provide insights into the placental response to environmental stress, offering potential preventative and therapeutic approaches for environment-related adverse outcomes in mothers and children., (© 2024. The Author(s).)

Dermatofibrosarcoma protuberans (DFSP) is a rare cutaneous sarcoma characterized by the COL1A1-PDGFB fusion gene. This study utilized single-cell RNA sequencing to dissect the cellular and molecular landscape of primary DFSP. Distinct DFSP cell clusters, exhibiting fibroblast-like traits, revealed variations in pathways associated with proliferation, inflammation and metabolism. Differential gene expression analysis during the differentiation from tumour stem cells to DFSP cells unveiled SMOC2, DCN and TGFBR3 as potential regulators of tumour invasion and immune infiltration through VEGF/TGF-β signalling modulation. Cellular communication analysis highlighted interactions within DFSP cell clusters and with endothelial cells, implicating molecules such as NAMPT, ANGPT2 and PTN in pathogenesis and treatment resistance. These findings offer insights into DFSP intratumour heterogeneity, elucidate molecular mechanisms underlying tumour behaviour, and suggest potential therapeutic targets., (© 2024 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Background: Neoadjuvant immunotherapy has evolved as an effective option to treat non-small cell lung cancer (NSCLC). B cells play essential roles in the immune system as well as cancer progression. However, the repertoire of B cells and its association with clinical outcomes remains unclear in NSCLC patients receiving neoadjuvant immunotherapy., Methods: Single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing data for LUAD samples were accessed from the TCGA and GEO databases. LUAD-related B cell marker genes were confirmed based on comprehensive analysis of scRNA-seq data. We then constructed the B cell marker gene signature (BCMGS) and validated it. In addition, we evaluated the association of BCGMS with tumor immune microenvironment (TIME) characteristics. Furthermore, we validated the efficacy of BCGMS in a cohort of NSCLC patients receiving neoadjuvant immunotherapy., Results: A BCMGS was constructed based on the TCGA cohort and further validated in three independent GSE cohorts. In addition, the BCMGS was proven to be significantly associated with TIME characteristics. Moreover, a relatively higher risk score indicated poor clinical outcomes and a worse immune response among NSCLC patients receiving neoadjuvant immunotherapy., Conclusions: We constructed an 18-gene prognostic signature derived from B cell marker genes based on scRNA-seq data, which had the potential to predict the prognosis and immune response of NSCLC patients receiving neoadjuvant immunotherapy., (© 2024. The Author(s), under exclusive licence to Federación de Sociedades Españolas de Oncología (FESEO).)

Background: Inflammatory heart disease can be triggered by a variety of causes, both infectious and noninfectious in nature. We hypothesized that inflammatory cardiomyopathy is potentially related to microbial infection., Methods: In this retrospective study, we used deep RNA sequencing on formalin-fixed paraffin-embedded heart tissue specimens to detect pathogenic agents. We first investigated 4 single-sample cases to test the feasibility of this diagnostic protocol and further 3 control-sample paired cases to improve the protocol with differential metatranscriptomics next-generation sequencing (mtNGS) analysis., Results: We demonstrate that differential mtNGS allows identification of various microbials as potentially pathogenic, for example, Cutibacterium acnes , Corynebacterium aurimucosum , and Pseudomonas denitrificans , which are usually commensal in healthy individuals. Differential mtNGS also allows characterization of human host response in each individual by profiling alterations of gene expression, networked pathways, and inferred immune cell compositions, information of which is beneficial for us to understand different etiologies and immunity roles in each case. Additionally, differential mtNGS allows the identification of genetic variants in patients that may contribute to their susceptibility to particular microbial infections., Conclusions: The demonstrated power of differential mtNGS in simultaneous capture of both the infectious microbial(s) and the status of human host immune response could help us better understand the pathogenesis of complex inflammatory cardiomyopathy, if conducted on a larger scale of the population. The developed differential mtNGS method could also shed light on its translation and adoption of such a laboratory test in clinic practice, allowing for a more effective diagnosis to guide therapeutic treatment of the disease., Competing Interests: Dr Dimitrova was an employee of Philips Healthcare North America. Dr Fallon received research funding from Philips Healthcare North America. The other authors report no conflicts.

Immune cell infiltration is a significant pathological process in abdominal aortic aneurysms (AAA). T cells, particularly CD4 + T cells, are essential immune cells responsible for substantial infiltration of the aorta. Regulatory T cells (Tregs) in AAA have been identified as tissue-specific; however, the time, location, and mechanism of acquiring the tissue-specific phenotype are still unknown. Using single-cell RNA sequencing (scRNA-seq) on CD4 + T cells from the AAA aorta and spleen, we discovered heterogeneity among CD4 + T cells and identified activated, proliferating and developed aorta Tregs. These Tregs originate in the peripheral tissues and acquire the tissue-specific phenotype in the aorta. The identification of precursors for Tregs in AAA provides new insight into the pathogenesis of AAA., (© 2024 The Authors. Immunology published by John Wiley & Sons Ltd.)

Colorectal cancer (CRC) is a prevalent cancer with high morbidity and mortality rates worldwide. Late diagnosis is a significant contributor to low survival rates in a minority of cases. The study aimed to perform a robust pipeline using integrated bioinformatics tools that will enable us to identify potential diagnostic and prognostic biomarkers for early detection of CRC by exploring differentially expressed genes (DEGs). In addition to, testing the capability of replacing chemotherapy with plant extract in CRC treatment by validating it using real-time PCR. RNA-seq data from cancerous and adjacent normal tissues were pre-processed and analyzed using various tools such as FastQC, Kallisto, DESeq@ R package, g:Profiler, GNEMANIA-CytoScape and CytoHubba, resulting in the identification of 1641 DEGs enriched in various signaling routes. MMP7, TCF21, and VEGFD were found to be promising diagnostic biomarkers for CRC. An in vitro experiment was conducted to examine the potential anticancer properties of 5-fluorouracile, Withania somnifera extract, and their combination. The extract was found to exhibit a positive trend in gene expression and potential therapeutic value by targeting the three genes; however, further trials are required to regulate the methylation promoter. Molecular docking tests supported the findings by revealing a stable ligand-receptor complex. In conclusion, the study's analysis workflow is precise and robust in identifying DEGs in CRC that may serve as biomarkers for diagnosis and treatment. Additionally, the identified DEGs can be used in future research with larger sample sizes to analyze CRC survival., (© 2023. The Author(s).)

Objectives: RNA sequencing (RNA-seq) and bioinformatic analysis were combined and used to explore the anti-inflammatory effects and mechanisms of naringenin (Nar) in lipopolysaccharide (LPS)-stimulated human periodontal ligament stem cells (hPDLSCs)., Methods: Cell counting kit-8, quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR), and enzyme-linked immunosorbent assay (ELISA) were adopted to detect the effects of Nar on the proliferation and expression of inflammatory factors in LPS-stimulated hPDLSCs, screening for the optimal anti-inflammatory concentration of Nar. Differentially expressed genes (DEGs) were screened using |log2FC|≥1 and P ≤0.05 as criteria. Volcano plot analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, the String database, and the MCODE module of Cytoscape were utilized to select core genes and enriched pathways. The effects on the nuclear factor κB (NF-κB) signaling pathway were verified using ELISA, qRT-PCR, and Western blot., Results: Appropriate concentrations of Nar could alleviate the expression of inflammatory factors and promote the proliferation of hPDLSCs stimulated by LPS. The best anti-inflammatory effect was achieved with 20 μmol/L Nar. RNA-seq showed significant enrichment of inflammation-related signaling pathways. The anti-inflammatory effect of Nar was mediated by inhibiting the NF-κB signaling pathway, similar to the effect of the NF-κB inhibitor BAY 11-7802., Conclusions: Nar could exert its anti-inflammatory effects by inhibiting the NF-κB signaling pathway, making it a potential therapeutic option for the adjuvant treatment of periodontitis.

Background: CD8+ T cells have been recognized as crucial factors in the prognosis of melanoma. However, there is currently a lack of gene markers that accurately describe their characteristics and functions in acral melanoma (AM), which hinders the development of personalized medicine., Methods: Firstly, we explored the composition differences of immune cells in AM using single-cell RNA sequencing (scRNA-seq) data and comprehensively characterized the immune microenvironment of AM in terms of composition, developmental differentiation, function, and cell communication. Subsequently, we constructed and validated a prognostic risk scoring model based on differentially expressed genes (DEGs) of CD8+ T cells using the TCGA-SKCM cohort through Lasso-Cox method. Lastly, immunofluorescence staining was performed to validate the expression of four genes (ISG20, CCL4, LPAR6, DDIT3) in AM and healthy skin tissues as included in the prognostic model., Results: The scRNA-seq data revealed that memory CD8+ T cells accounted for the highest proportion in the immune microenvironment of AM, reaching 70.5%. Cell-cell communication analysis showed extensive communication relationships among effector CD8+ T cells. Subsequently, we constructed a prognostic scoring model based on DEGs derived from CD8+ T cell sources. Four CD8+ T cell-related genes were included in the construction and validation of the prognostic model. Additionally, immunofluorescence results demonstrated that ISG20 and CCL4 were downregulated, while LPAR6 and DDIT3 were upregulated in AM tissues compared to normal skin tissues., Conclusion: Identifying biomarkers based on the expression levels of CD8+ T cell-related genes may be an effective approach for establishing prognostic models in AM patients. The independently prognostic risk evaluation model we constructed provides new insights and theoretical support for immunotherapy in AM., (© 2024 The Author(s). Skin Research and Technology published by John Wiley & Sons Ltd.)

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating (nascent) transcripts across diverse crops and other plants using capped small (cs)RNA sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many 'distal' elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Purpose: The purpose of this study was to elucidate the relationship between the tumor microenvironment (TME) and cellular diversity in bladder cancer (BLCA) progression, leveraging single-cell RNA sequencing (scRNA-seq) data to identify potential prognostic biomarkers and construct a prognostic model for BLCA., Methods: We analyzed scRNA-seq data of normal and tumor bladder cells from the Gene Expression Omnibus (GEO) database to uncover crucial markers within the bladder TME. The study compared gene expression in normal versus tumor bladder cells, identifying differentially expressed genes. These genes were subsequently assessed for their prognostic significance using patient follow-up data from The Cancer Genome Atlas. Prognostic models were constructed using Least Absolute Shrinkage and Selection Operator and multivariate Cox regression analyses, focusing on eight genes of interest. The predictive performance of the model was also tested against additional GEO data sets (GSE31684, GSE13507, and GSE32894)., Results: The prognostic model demonstrated reliable prediction of patient outcomes. Validation through gene set enrichment analysis and immune cell infiltration assessment supported the model's efficacy. The results from both the univariate and multivariate analyses suggest that the risk score is an independent prognostic factor with a hazard ratio of 2.97 (95% CI, 2.28 to 3.9, P < .001). In the validation cohort, the AUC at 1, 2, and 3 years is 0.74, 0.74, and 0.72, respectively., Conclusion: Our findings proposed biomarkers with prognostic potential, laying the groundwork for future in vitro validation and therapeutic exploration. This contributes to a deeper understanding of the genes associated with bladder TME and may improve prognostic precision in BLCA management.

Background: Association testing between molecular phenotypes and genomic variants can help to understand how genotype affects phenotype. RNA sequencing provides access to molecular phenotypes such as gene expression and alternative splicing while DNA sequencing or microarray genotyping are the prevailing options to obtain genomic variants., Results: We genotype variants for 74 male Braunvieh cattle from both DNA (~ 13-fold coverage) and deep total RNA sequencing from testis, vas deferens, and epididymis tissue (~ 250 million reads per tissue). We show that RNA sequencing can be used to identify approximately 40% of variants (7-10 million) called from DNA sequencing, with over 80% precision. Within highly expressed coding regions, over 92% of expected variants were called with nearly 98% precision. Allele-specific expression and putative post-transcriptional modifications negatively impact variant genotyping accuracy from RNA sequencing and contribute to RNA-DNA differences. Variants called from RNA sequencing detect roughly 75% of eGenes identified using variants called from DNA sequencing, demonstrating a nearly 2-fold enrichment of eQTL variants. We observe a moderate-to-strong correlation in nominal association p-values (Spearman ρ 2  ~ 0.6), although only 9% of eGenes have the same top associated variant., Conclusions: We find hundreds of thousands of RNA-DNA differences in variants called from RNA and DNA sequencing on the same individuals. We identify several highly significant eQTL when using RNA sequencing variant genotypes which are not found with DNA sequencing variant genotypes, suggesting that using RNA sequencing variant genotypes for association testing results in an increased number of false positives. Our findings demonstrate that caution must be exercised beyond filtering for variant quality or imputation accuracy when analysing or imputing variants called from RNA sequencing., (© 2024. The Author(s).)

A 26-year-old man presented with a tumor in the left soleus muscle. The tumor was diagnosed as a locally advanced leiomyosarcoma. The patient was treated with irradiation followed by wide resection. One year after surgery, the patient presented with multiple lung metastases. Despite aggressive sequential chemotherapy, systemic metastatic tumors continued to develop. To explore therapeutic options for the patient, we performed DNA-based CGP with FoundationOne ® CDx (F1). F1 identified an out-of-strand rearrangement of the NOS1AP::NTRK1 gene, which has not been previously reported. In contrast, RNA sequencing revealed an in-frame LMNA::NTRK1 gene, which is an oncogenic fusion gene.

As an important warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) flourishes in warm areas around the world due to the existence of the C4 photosynthetic pathway. However, how C4 photosynthesis operates in bermudagrass leaves is still poorly understood. In this study, we performed single-cell RNA-sequencing on 5296 cells from bermudagrass leaf blades. Eight cell clusters corresponding to mesophyll, bundle sheath, epidermis and vascular bundle cells were successfully identified using known cell marker genes. Expression profiling indicated that genes encoding NADP-dependent malic enzymes (NADP-MEs) were highly expressed in bundle sheath cells, whereas NAD-ME genes were weakly expressed in all cell types, suggesting C4 photosynthesis of bermudagrass leaf blades might be NADP-ME type rather than NAD-ME type. The results also indicated that starch synthesis-related genes showed preferential expression in bundle sheath cells, whereas starch degradation-related genes were highly expressed in mesophyll cells, which agrees with the observed accumulation of starch-filled chloroplasts in bundle sheath cells. Gene co-expression analysis further revealed that different families of transcription factors were co-expressed with multiple C4 photosynthesis-related genes, suggesting a complex transcription regulatory network of C4 photosynthesis might exist in bermudagrass leaf blades. These findings collectively provided new insights into the cell-specific expression patterns and transcriptional regulation of photosynthetic genes in bermudagrass., 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 Masson SAS. All rights reserved.)

Peas (Pisum sativum L.) are widely cultivated in temperate regions and are susceptible hosts for various viruses across different families. The discovery and identification of new viruses in peas has significant implications for field disease management. Here, we identified a mixed infection of two viruses from field-collected peas exhibiting virus-like symptoms using metatranscriptome and small RNA sequencing techniques. Upon identification, one of the viruses was determined to be a newly isolated and discovered bymovirus from peas, named "pea bymovirus 1 (PBV1)". The other was identified as a novel variant of bean yellow mosaic virus (BYMV-HZ1). Subsequently, mechanical inoculation and RT-PCR assays confirmed that both viruses could be inoculated back onto peas and tobaccos, showing mixed infection by PBV1 and BYMV-HZ1. To our knowledge, this is the first isolation of a bymovirus from pea and the first documented case of mixed infection of peas by PBV1 and BYMV-HZ1 in China., 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 Inc. All rights reserved.)

Background: Cancer-associated fibroblasts (CAFs) are found in primary and advanced tumours. They are primarily involved in tumour progression through complex mechanisms with other types of cells in the tumour microenvironment. However, essential fibroblasts-related genes (FRG) in bladder cancer still need to be explored, and there is a shortage of an ideal predictive model or molecular subtype for the progression and immune therapeutic assessment for bladder cancer, especially muscular-invasive bladder cancer based on the FRG., Materials and Methods: CAF-related genes of bladder cancer were identified by analysing single-cell RNA sequence datasets, and bulk transcriptome datasets and gene signatures were used to characterize them. Then, 10 types of machine learning algorithms were utilised to determine the hallmark FRG and construct the FRG index (FRGI) and subtypes. Further molecular subtypes combined with CD8+ T-cells were established to predict the prognosis and immune therapy response., Results: Fifty-four BLCA-related FRG were screened by large-scale scRNA-sequence datasets. The machine learning algorithm established a 3-genes FRGI. High FRGI represented a worse outcome. Then, FRGI combined clinical variables to construct a nomogram, which shows high predictive performance for the prognosis of bladder cancer. Furthermore, the BLCA datasets were separated into two subtypes - fibroblast hot and cold types. In five independent BLCA cohorts, the fibroblast hot type showed worse outcomes than the cold type. Multiple cancer-related hallmark pathways are distinctively enriched in these two types. In addition, high FRGI or fibroblast hot type shows a worse immune therapeutic response. Then, four subtypes called CD8-FRG subtypes were established under the combination of FRG signature and activity of CD8+ T-cells, which turned out to be effective in predicting the prognosis and immune therapeutic response of bladder cancer in multiple independent datasets. Pathway enrichment analysis, multiple gene signatures, and epigenetic alteration characterize the CD8-FRG subtypes and provide a potential combination strategy method against bladder cancer., Conclusions: In summary, the authors established a novel FRGI and CD8-FRG subtype by large-scale datasets and organised analyses, which could accurately predict clinical outcomes and immune therapeutic response of BLCA after surgery., (Copyright © 2024 The Author(s). Published by Wolters Kluwer Health, Inc.)

Background: Abdominal obesity increases the risk for non-alcoholic fatty liver disease (NAFLD), now known as metabolic dysfunction-associated steatotic liver disease (MASLD)., Methods: To elucidate the directional cell-type level biological mechanisms underlying the association between abdominal obesity and MASLD, we integrated adipose and liver single nucleus RNA-sequencing and bulk cis-expression quantitative trait locus (eQTL) data with the UK Biobank genome-wide association study (GWAS) data using colocalization. Then we used colocalized cis-eQTL variants as instrumental variables in Mendelian randomization (MR) analyses, followed by functional validation experiments on the target genes of the cis-eQTL variants., Findings: We identified 17 colocalized abdominal obesity GWAS variants, regulating 17 adipose cell-type marker genes. Incorporating these 17 variants into MR discovers a putative tissue-of-origin, cell-type-aware causal effect of abdominal obesity on MASLD consistently with multiple MR methods without significant evidence for pleiotropy or heterogeneity. Single cell data confirm the adipocyte-enriched mean expression of the 17 genes. Our cellular experiments across human adipogenesis identify risk variant -specific epigenetic and transcriptional mechanisms. Knocking down two of the 17 genes, PPP2R5A and SH3PXD2B, shows a marked decrease in adipocyte lipidation and significantly alters adipocyte function and adipogenesis regulator genes, including DGAT2, LPL, ADIPOQ, PPARG, and SREBF1. Furthermore, the 17 genes capture a characteristic MASLD expression signature in subcutaneous adipose tissue., Interpretation: Overall, we discover a significant cell-type level effect of abdominal obesity on MASLD and trace its biological effect to adipogenesis., Funding: NIH grants R01HG010505, R01DK132775, and R01HL170604; the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant No. 802825), Academy of Finland (Grants Nos. 333021), the Finnish Foundation for Cardiovascular Research the Sigrid Jusélius Foundation and the Jane and Aatos Erkko Foundation; American Association for the Study of Liver Diseases (AASLD) Advanced Transplant Hepatology award and NIH/NIDDK (P30DK41301) Pilot and Feasibility award; NIH/NIEHS F32 award (F32ES034668); Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2021), the Academy of Finland grant (Contract no. 138006); Academy of Finland (Grant Nos 335443, 314383, 272376 and 266286), Sigrid Jusélius Foundation, Finnish Medical Foundation, Finnish Diabetes Research Foundation, Novo Nordisk Foundation (#NNF20OC0060547, NNF17OC0027232, NNF10OC1013354) and Government Research Funds to Helsinki University Hospital; Orion Research Foundation, Maud Kuistila Foundation, Finish Medical Foundation, and University of Helsinki., Competing Interests: Declaration of interests M.L. was funded by the Sigrid Jusélius Foundation during the last 36 months. The other authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier B.V. All rights reserved.)