Molecular Subtypes of Bladder Cancer: Component Signatures and Potential Value in Clinical Decision-making

Bladder cancer may be classified into "molecular subtypes" based on gene expression. These are associated with treatment response and patient outcomes. The gene expression signatures that define these subtypes are diverse, including signatures of epithelial differentiation, stromal involvement, cell cycle activity, and immune cell infiltration. Multiple different systems are described. While earlier studies considered molecular subtypes to be intrinsic properties of cancer, recent data have shown molecular subtypes change as tumors progress and evolve, and often differ between histologically distinct regions of a tumor. The data also indicate that some signatures that define molecular subtypes may be treated as independent continuous variables, rather than categorical subtypes, and these individual signatures may be more clinically informative. This review describes molecular subtypes of urothelial carcinoma, including histologic subtypes and tumors with divergent differentiation, and explores potential future uses in patient management.

Prediction of Responsiveness to PD-L1/PD-1 Inhibitors Using miRNA Profiles Associated With PD-L1 Expression in Lung Adenocarcinoma and Squamous Cell Carcinoma

BACKGROUND/AIM

MicroRNAs (miRNAs) regulate programmed cell death ligand 1 (PD-L1) and play a crucial role in tumor immune response. However, the relationship between miRNA expression patterns and PD-L1 remains unclear in both lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). We investigated PD-L1-related miRNAs that can predict treatment response in patients treated with PD-L1/PD-1 inhibitors.

PATIENTS AND METHODS

We selected miRNAs that were correlated with PD-L1 expression within the LUAD and LUSC datasets obtained from The Cancer Genome Atlas (TCGA) and Clinical Proteomic Tumor Analysis Consortium (CPTAC). We validated whether the miRNA profile could be used to predict the prognosis of patients treated with PD-L1/PD-1 inhibitors.

RESULTS

Based on four public datasets, we selected 66 and 23 miRNAs associated with PD-L1 expression in LUAD and LUSC, respectively. From the above miRNAs, we identified 5 miRNAs in LUSC and 1 miRNA in LUAD that could predict the response to PD-L1/PD-1 inhibitors in a validation set of patients treated with PD-L1/PD-1 inhibitors. In LUSC, the miRNA profile exhibited a high predictive capability for the response to PD-L1/PD-1 treatment [area under the curve (AUC)=0.963] and accurately predicted prognosis (p=0.031). In LUAD, the miRNA profile was relatively less predictive than in LUSC (AUC=0.691 and p=0.213). Additionally, we observed variations in the PD-L1-associated miRNA profiles, as well as in the associated pathways, between LUAD and LUSC.

CONCLUSION

The PD-L1-associated miRNA profile may predict treatment response in LUSC patients treated with PD-L1/PD-1 inhibitors and help select the PD-L1/PD-1 inhibitor treatment group.

The impact of exercise on gene regulation in association with complex trait genetics

Endurance exercise training is known to reduce risk for a range of complex diseases. However, the molecular basis of this effect has been challenging to study and largely restricted to analyses of either few or easily biopsied tissues. Extensive transcriptome data collected across 15 tissues during exercise training in rats as part of the Molecular Transducers of Physical Activity Consortium has provided a unique opportunity to clarify how exercise can affect tissue-specific gene expression and further suggest how exercise adaptation may impact complex disease-associated genes. To build this map, we integrate this multi-tissue atlas of gene expression changes with gene-disease targets, genetic regulation of expression, and trait relationship data in humans. Consensus from multiple approaches prioritizes specific tissues and genes where endurance exercise impacts disease-relevant gene expression. Specifically, we identify a total of 5523 trait-tissue-gene triplets to serve as a valuable starting point for future investigations [Exercise; Transcription; Human Phenotypic Variation].

Comparative transcriptome analysis of bone marrow resident versus culture-expanded mouse mesenchymal stem/stromal cells

BACKGROUND AIMS

Mesenchymal stem/stromal cells (MSCs) are defined as culture-expanded populations, and although these cells recapitulate many properties of bone marrow (BM) resident skeletal stem/progenitor cells, few studies have directly compared these populations to evaluate how culture adaptation and expansion impact critical quality attributes.

METHODS

We analyzed by RNA sequencing Lin-SCA1+ MSCs enriched from BM by immunodepletion (ID) and after subsequent culture expansion (Ex) and Lin-LEPR+ MSCs sorted (S) directly from BM. Pairwise comparisons were used to identify differentially expressed genes (DEGs) between populations, and gene set enrichment analysis was employed to identify biological pathways/processes unique to each population. K-means cluster analysis resolved isolation status-dependent changes in transcription in pseudotime.

RESULTS

Hierarchical clustering segregated populations by isolation process, and principal component analysis identified transcripts related to vasculature development, ossification and inflammatory/cytokine signaling as key drivers of population variance. Pairwise comparisons identified 3849 DEGs in ID versus S BM-MSCs mapping to Gene Ontology (GO) terms related to immune and metabolic processes and 334 DEGs in Ex versus ID BM-MSCs mapping to GO terms related to tissue development, cell growth and replication and organelle organization. K-means cluster analysis revealed significant differences in transcripts encoding stemness and differentiation markers, extracellular matrix structural constituents and remodeling enzymes and paracrine-acting factors between populations.

CONCLUSIONS

These comparative analyses reveal significant differences in gene expression signatures between BM resident and culture-expanded MSCs, thereby providing new insight into how culture adaptation/expansion endows the latter with unique quality attributes.

Single-cell profiling transcriptomic reveals cellular heterogeneity and cellular crosstalk in choroidal neovascularization model

Choroidal neovascularization (CNV) is a hallmark of neovascular age-related macular degeneration (nAMD) and a major contributor to vision loss in nAMD cases. However, the identification of specific cell types associated with nAMD remains challenging. Herein, we performed single-cell sequencing to comprehensively explore the cellular diversity and understand the foundational components of the retinal pigment epithelium (RPE)/choroid complex. We unveiled 10 distinct cell types within the RPE/choroid complex. Notably, we observed significant heterogeneity within endothelial cells (ECs), fibroblasts, and macrophages, underscoring the intricate nature of the cellular composition in the RPE/choroid complex. Within the EC category, four distinct clusters were identified and EC cluster 0 was tightly associated with choroidal neovascularization. We identified five clusters of fibroblasts actively involved in the pathogenesis of nAMD, influencing fibrotic responses, angiogenic effects, and photoreceptor function. Additionally, three clusters of macrophages were identified, suggesting their potential roles in regulating the progression of nAMD through immunomodulation and inflammation regulation. Through CellChat analysis, we constructed a complex cell-cell communication network, revealing the role of EC clusters in interacting with fibroblasts and macrophages in the context of nAMD. These interactions were found to govern angiogenic effects, fibrotic responses, and inflammatory processes. In summary, this study reveals noteworthy cellular heterogeneity in the RPE/choroid complex and provides valuable insights into the pathogenesis of CNV. These findings will open up potential avenues for deep understanding and targeted therapeutic interventions in nAMD.

Kidney toxicity and transcriptome analyses of male ICR mice acutely exposed to the mushroom toxin α-amanitin

Amatoxins are responsible for most fatal mushroom poisoning cases, as it causes both hepatotoxicity and nephrotoxicity. However, studies on amatoxin nephrotoxicity are limited. Here, we investigated nephrotoxicity over 4 days and nephrotoxicity/hepatotoxicity over 14 days in mice. The organ weight ratio, serological indices, and tissue histology results indicated that a nephrotoxicity mouse model was established with two stages: (1) no apparent effects within 24 h; and (2) the appearance of adverse effects, with gradual worsening within 2-14 days. For each stage, the kidney transcriptome revealed patterns of differential mRNA expression and significant pathway changes, and Western blot analysis verified the expression of key proteins. Amanitin-induced nephrotoxicity was directly related to RNA polymerase II because mRNA levels decreased, RNA polymerase II-related pathways were significantly enriched at the transcription level, and RNA polymerase II protein was degraded in the early poisoning stage. In the late stage, nephrotoxicity was more severe than hepatotoxicity. This is likely associated with inflammation because inflammation-related pathways were significantly enriched and NF-κB activation was increased in the kidney.

Phenotypic and transcriptome profiling of spikes reveals the regulation of light regimens on spike growth and fertile floret number in wheat

The spike growth phase is critical for the establishment of fertile floret (grain) numbers in wheat (Triticum aestivum L.). Then, how to shorten the spike growth phase and increase grain number synergistically? Here, we showed high-resolution analyses of floret primordia (FP) number, morphology and spike transcriptomes during the spike growth phase under three light regimens. The development of all FP in a spike could be divided into four distinct stages: differentiation (Stage I), differentiation and morphology development concurrently (Stage II), morphology development (Stage III), and polarization (Stage IV). Compared to the short photoperiod, the long photoperiod shortened spike growth and stimulated early flowering by shortening Stage III; however, this reduced assimilate accumulation, resulting in fertile floret loss. Interestingly, long photoperiod supplemented with red light shortened the time required to complete Stages I-II, then raised assimilates supply in the spike and promoted anther development before polarization initiation, thereby increasing fertile FP number during Stage III, and finally maintained fertile FP development during Stage IV until they became fertile florets via a predicted dynamic gene network. Our findings proposed a light regimen, critical stages and candidate regulators that achieved a shorter spike growth phase and a higher fertile floret number in wheat.

microRNA-9a-5p disrupts the ELAVL1/VEGF axis to alleviate traumatic brain injury

Plasma microRNA (miR)-9 has been identified as a promising diagnostic biomarker for traumatic brain injury (TBI). This study aims to investigate the possible role and mechanisms of miR-9a-5p affecting TBI. Microarray-based gene expression profiling of TBI was used for screening differentially expressed miRNAs and genes. TBI rat models were established. miR-9a-5p, ELAVL1 and VEGF expression in the brain tissue of TBI rats was detected. The relationship among miR-9a-5p, ELAVL1 and VEGF was tested. TBI modeled rats were injected with miR-9a-5p-, ELAVL1 or VEGF-related sequences to identify their effects on TBI. miR-9a-5p was poorly expressed in the brain tissue of rats with TBI. ELAVL1 was a downstream target gene of miR-9a-5p, which could negatively regulate its expression. Enforced miR-9a-5p expression prevented brain tissue damage in TBI rats by targeting ELAVL1. Meanwhile, ELAVL1 could increase the expression of VEGF, which was highly expressed in the brain tissue of rats with TBI. In addition, ectopically expressed miR-9a-5p alleviated brain tissue damage in TBI rats by downregulating the ELAVL1/VEGF axis. Overall, miR-9a-5p can potentially reduce brain tissue damage in TBI rats by targeting ELAVL1 and down-regulating VEGF expression.

Transcriptome analysis provides new insights into the immune response of Ruditapes philippinarum infected with Vibrio alginolyticus

Manila clam (Ruditapes philippinarum) is a bivalve species with commercial value, but it is easily infected by pathogenic microorganisms in aquaculture, which restricts the shellfish industry. Notably, the impact of Vibrio alginolyticus on clam culture is obvious. In this study, RNA-seq was performed to analyze clam hepatopancreas tissue in 48 h (challenge group, G48h) and 96 h (challenge group, G96h) after infection with V. alginolyticus and 0 h after injection of PBS (control group, C). The results showed that a total of 1670 differentially expressed genes were detected in the G48h vs C group, and 1427 differentially expressed genes were detected in the G96h vs C group. In addition, KEGG analysis showed that DEGs were significantly enriched in pathways such as Lysosome and Mitophagy. Moreover, 15 immune related DEGs were selected for qRT-PCR analysis to verify the accuracy of RNA-seq, and the results showed that the expression level of DEGs was consistent with that of RNA-seq. Therefore, the results obtained in this study provides a preliminary understanding of the immune defense of R. philippinarum and molecular insights for genetic breeding of V. alginolyticus resistance in Manila clam.

Revealing the pathogenesis of keloids based on the status: Active vs inactive

Recently, the pathomechanisms of keloids have been extensively researched using transcriptomic analysis, but most studies did not consider the activity of keloids. We aimed to profile the transcriptomics of keloids according to their clinical activity and location within the keloid lesion, compared with normal and mature scars. Tissue samples were collected (keloid based on its activity (active and inactive), mature scar from keloid patients and normal scar (NS) from non-keloid patients). To reduce possible bias, all keloids assessed in this study had no treatment history and their location was limited to the upper chest or back. Multiomics assessment was performed by using single-cell RNA sequencing and multiplex immunofluorescence. Increased mesenchymal fibroblasts (FBs) was the main feature in keloid patients. Noticeably, the proportion of pro-inflammatory FBs was significantly increased in active keloids compared to inactive ones. To explore the nature of proinflammatory FBs, trajectory analysis was conducted and CCN family associated with mechanical stretch exhibited higher expression in active keloids. For vascular endothelial cells (VECs), the proportion of tip and immature cells increased in keloids compared to NS, especially at the periphery of active keloids. Also, keloid VECs highly expressed genes with characteristics of mesenchymal activation compared to NS, especially those from the active keloid center. Multiomics analysis demonstrated the distinct expression profile of active keloids. Clinically, these findings may provide the future appropriate directions for development of treatment modalities of keloids. Prevention of keloids could be possible by the suppression of mesenchymal activation between FBs and VECs and modulation of proinflammatory FBs may be the key to the control of active keloids.

Identification of key mRNAs and signaling pathways in obsessive compulsive disorder based on weighted gene co-expression network analysis and cytoHubba plugin

PURPOSE

Obsessive-compulsive disorder (OCD) is a complex psychiatric disorder. Genetic and broad environmental factors are common risk factors for OCD. The purpose of this study is to explore the molecular mechanism of OCD and to find new molecular targets for the diagnosis and management of OCD.

METHODS

All data were downloaded from public dataset. Key modules and candidate key mRNAs were identified based on weighted gene co-expression network analysis (WGCNA). The "limma" R package was used for differential expression analysis of mRNAs. Subsequently, functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) was also carried out. In addition, a diagnostic model was constructed. Finally, the infiltration level of immune cells in OCD and its correlation with multicentric key DEmRNAs were analyzed.

RESULTS

Green and red modules were selected as the hub modules. A total of 447 mRNAs were considered candidate key mRNAs according to GS > 0.2 and MM > 0.3. A total of 26 DEmRNAs in the same direction were identified in the GSE60190 and GSE78104 datasets. A total of 26 DEmRNAs were intersected with candidate key mRNAs in WGCNA to obtain 10 intersection DEmRNAs (HSPB1, ITPK1, CBX7, PPP1R10, TAOK1, PISD, MKNK2, RWDD1, PPA1, and RELN). However, only four DEmRNAs (HSPB1, TAOK1, MKNK2, and PPA1) predicted related drugs. Subsequently, receiver operating characteristic analysis shows that the diagnostic model has high diagnostic value. Moreover, six multicentric key DEmRNAs (SNRPF, SNRNP70, PRPF8, NOP56, EPRS, and CCT2) were screened by UpSet package. Finally, six multicentric key DEmRNAs were found to be associated with immune cells.

CONCLUSION

The key molecules obtained in this study lay a foundation for further research on the molecular mechanism of OCD.

Signature reversion of three disease-associated gene signatures prioritizes cancer drug repurposing candidates

Drug repurposing is promising because approving a drug for a new indication requires fewer resources than approving a new drug. Signature reversion detects drug perturbations most inversely related to the disease-associated gene signature to identify drugs that may reverse that signature. We assessed the performance and biological relevance of three approaches for constructing disease-associated gene signatures (i.e., limma, DESeq2, and MultiPLIER) and prioritized the resulting drug repurposing candidates for four low-survival human cancers. Our results were enriched for candidates that had been used in clinical trials or performed well in the PRISM drug screen. Additionally, we found that pamidronate and nimodipine, drugs predicted to be efficacious against the brain tumor glioblastoma (GBM), inhibited the growth of a GBM cell line and cells isolated from a patient-derived xenograft (PDX). Our results demonstrate that by applying multiple disease-associated gene signature methods, we prioritized several drug repurposing candidates for low-survival cancers.

Early Onset Intrahepatic Cholangiocarcinoma: Clinical Characteristics, Oncological Outcomes, and Genomic/Transcriptomic Features

INTRODUCTION

Data on clinical characteristics and disease-specific prognosis among patients with early onset intrahepatic cholangiocarcinoma (ICC) are currently limited.

METHODS

Patients undergoing hepatectomy for ICC between 2000 and 2020 were identified by using a multi-institutional database. The association of early (≤50 years) versus typical onset (>50 years) ICC with recurrence-free (RFS) and disease-specific survival (DSS) was assessed in the multi-institutional database and validated in an external cohort. The genomic and transcriptomic profiles of early versus late onset ICC were analyzed by using the Total Cancer Genome Atlas (TCGA) and Memorial Sloan Kettering Cancer Center databases.

RESULTS

Among 971 patients undergoing resection for ICC, 22.7% (n = 220) had early-onset ICC. Patients with early-onset ICC had worse 5-year RFS (24.1% vs. 29.7%, p < 0.05) and DSS (36.5% vs. 48.9%, p = 0.03) compared with patients with typical onset ICC despite having earlier T-stage tumors and lower rates of microvascular invasion. In the validation cohort, patients with early-onset ICC had worse 5-year RFS (7.4% vs. 20.5%, p = 0.002) compared with individuals with typical onset ICC. Using the TCGA cohort, 652 and 266 genes were found to be upregulated (including ATP8A2) and downregulated (including UTY and KDM5D) in early versus typical onset ICC, respectively. Genes frequently implicated as oncogenic drivers, including CDKN2A, IDH1, BRAF, and FGFR2 were infrequently mutated in the early-onset ICC patients.

CONCLUSIONS

Early-onset ICC has distinct clinical and genomic/transcriptomic features. Morphologic and clinicopathologic characteristics were unable to fully explain differences in outcomes among early versus typical onset ICC patients. The current study offers a preliminary landscape of the molecular features of early-onset ICC.

Prognostic value and potential mechanism of cellular senescence and tumor microenvironment in hepatocellular carcinoma: Insights from bulk transcriptomics and single-cell sequencing analysis

The high mortality rate and postoperative recurrence of hepatocellular carcinoma (HCC) contribute to the burden on society and healthcare. The prognostic value and underlying mechanisms of cellular senescence and tumor microenvironment (TME) in HCC remain unclear. Bulk transcriptomic data were obtained from 368 HCC samples in The Cancer Genome Atlas-liver hepatocellular carcinoma cohort and 64 samples from the GSE116174 dataset. Single-cell RNA sequencing (scRNA-seq) data of HCC were obtained from the GSE149614 dataset, including 18 tumor samples from 10 patients. Prognosis-related cellular senescence genes and immune cells were identified through univariate analysis. Least absolute shrinkage and selection operator regression analysis was performed to construct the CellAge score and TME score, both of which were identified as independent prognostic factors for HCC based on multivariate Cox analysis. The combined CellAge and TME scores showed improved prognostic stratification for HCC patients, as confirmed by multivariate Cox analysis (p < .001). The gene set enrichment analysis (GSEA) revealed enrichment of the extracellular matrix receptor interaction signaling pathway in the group with high CellAge scores and low TME scores, which exhibited a worse prognosis. Single-cell sequencing results revealed higher expression activity of the cAMP response element modulator (CREM) extended transcription factor in HCC cells and most immune cells, indicating its involvement in TME remodeling. Finally, the tumor immune dysfunction and exclusion (TIDE) analysis demonstrated that the combined scores could predict the outcomes of immune therapy in patients with HCC. In conclusion, cellular senescence contributes to TME remodeling in HCC, and the developed CellAge and TME scores serve as independent prognostic factors and predictors of immune therapy in HCC.

Single-cell, single-nucleus, and spatial transcriptomics characterization of the immunological landscape in the healthy and PSC human liver

BACKGROUND & AIMS

Primary sclerosing cholangitis (PSC) is an immune-mediated cholestatic liver disease for which there is an unmet need to understand the cellular composition of the affected liver and how it underlies disease pathogenesis. We aimed to generate a comprehensive atlas of the PSC liver using multi-omic modalities and protein-based functional validation.

METHODS

We employed single-cell and single-nucleus RNA sequencing (47,156 cells and 23,000 nuclei) and spatial transcriptomics (one sample by 10x Visium and five samples with Nanostring GeoMx DSP) to profile the cellular ecosystem in 10 PSC livers. Transcriptomic profiles were compared to 24 neurologically deceased donor livers (107,542 cells) and spatial transcriptomics controls, as well as 18,240 cells and 20,202 nuclei from three PBC livers. Flow cytometry was performed to validate PSC-specific differences in immune cell phenotype and function.

RESULTS

PSC explants with parenchymal cirrhosis and prominent periductal fibrosis contained a population of cholangiocyte-like hepatocytes that were surrounded by diverse immune cell populations. PSC-associated biliary, mesenchymal, and endothelial populations expressed chemokine and cytokine transcripts involved in immune cell recruitment. Additionally, expanded CD4+ T cells and recruited myeloid populations in the PSC liver expressed the corresponding receptors to these chemokines and cytokines, suggesting potential recruitment. Tissue-resident macrophages, by contrast, were reduced in number and exhibited a dysfunctional and downregulated inflammatory response to lipopolysaccharide and interferon-γ stimulation.

CONCLUSIONS

We present a comprehensive atlas of the PSC liver and demonstrate an exhaustion-like phenotype of myeloid cells and markers of chronic cytokine expression in late-stage PSC lesions. This atlas expands our understanding of the cellular complexity of PSC and has potential to guide the development of novel treatments.

IMPACT AND IMPLICATIONS

Primary sclerosing cholangitis (PSC) is a rare liver disease characterized by chronic inflammation and irreparable damage to the bile ducts, which eventually results in liver failure. Due to a limited understanding of the underlying pathogenesis of disease, treatment options are limited. To address this, we sequenced healthy and diseased livers to compare the activity, interactions, and localization of immune and non-immune cells. This revealed that hepatocytes lining PSC scar regions co-express cholangiocyte markers, whereas immune cells infiltrate the scar lesions. Of these cells, macrophages, which typically contribute to tissue repair, were enriched in immunoregulatory genes and demonstrated a lack of responsiveness to stimulation. These cells may be involved in maintaining hepatic inflammation and could be a target for novel therapies.

Cross-Talking Pathways of Rapidly Accelerated Fibrosarcoma-1 (RAF-1) in Alzheimer's Disease

Alzheimer's disease (AD) becomes one of the main global burden diseases with the aging population. This study was to investigate the potential molecular mechanisms of rapidly accelerated fibrosarcoma-1 (RAF-1) in AD through bioinformatics analysis. Differential gene expression analysis was performed in GSE132903 dataset. We used weight gene correlation network analysis (WGCNA) to evaluate the relations among co-expression modules and construct global regulatory network. Cross-talking pathways of RAF-1 in AD were identified by functional enrichment analysis. Totally, 2700 differentially expressed genes (DEGs) were selected between AD versus non-dementia control and RAF-1-high versus low group. Among them, DEGs in turquoise module strongly associated with AD and high expression of RAF-1 were enriched in vascular endothelial growth factor (VEGF), neurotrophin, mitogen-activated protein kinase (MAPK) signaling pathway, oxidative phosphorylation, GABAergic synapse, and axon guidance. Moreover, cross-talking pathways of RAF-1, including MAPK, VEGF, neurotrophin signaling pathways, and axon guidance, were identified by global regulatory network. The performance evaluation of AUC was 84.2%. The gene set enrichment analysis (GSEA) indicated that oxidative phosphorylation and synapse-related biological processes were enriched in RAF-1-high and AD group. Our findings strengthened the potential roles of high RAF-1 level in AD pathogenesis, which were mediated by MAPK, VEGF, neurotrophin signaling pathways, and axon guidance.

The landscape of immune dysregulation in pediatric sepsis at a single-cell resolution

Recognizing immune dysregulation as a hallmark of sepsis pathophysiology, leukocytes have attracted major attention of investigation. While adult and pediatric sepsis are clinically distinct, their immunological delineation remains limited. Single cell technologies facilitated the characterization of immune signatures. We tackled to delineate immunological profiles of pediatric sepsis at a single-cell level by analyzing blood samples from six septic children, at both acute and recovery phases, and four healthy children. 16 single-cell transcriptomic datasets were analyzed and compared to adult sepsis dataset. We showed a unique shift in neutrophil subpopulations and functions between acute and recovery phases, along with the regulatory role of resistin. Neutrophil signatures were comparable between adult and pediatric sepsis. Innate-like CD4 T cells were predominantly and uniquely observed in acute phase of pediatric sepsis. Our study serves as a rich source of information about the phenotypic diversity and trajectory of circulating immune cells during pediatric sepsis.

Hepatocellular Carcinoma Epigenetic Patterns Correspond to Differences in Ethnoracial Status and Treatment Response in a Single-Center Retrospective Study

PURPOSE

To correlate epigenetic patterns with ethnoracial status and locoregional therapy (LRT) response in patients with hepatocellular carcinoma (HCC).

MATERIALS AND METHODS

DNA and RNA were extracted from 47 distinct formalin-fixed paraffin-embedded tumor samples from 42 patients with HCC (n = 14 Black, n = 19 White, n = 9 Hispanic). LRT response was determined using computed tomography (CT) or magnetic resonance (MR) imaging 3 months posttreatment of 35 tumors (n = 22 complete response, n = 13 retreatment candidates). RNA expression and DNA methylation were used to stratify patients by ethnoracial status and treatment response using partial least-squares discriminant analysis (PLS-DA). Results were validated using hierarchical clustering. Ingenuity pathway analysis was performed to identify upstream regulators and pathways.

RESULTS

PLS-DA identified 100 genes and 12 methylated regions that differentiated tumors from Black from White/Hispanic patients. Hierarchical clustering clustered samples with the top 16 genes or the top 5 methylation regions. Dysregulated pathways included adrenomedullin pathway (P = .030), EIF2 signaling (P = .007), and several metabolic pathways. AGTR1 (log2fold = 1.59) and GSTM3 (log2fold = 2.53) represented potential differentially expressed therapeutic targets. PLS-DA identified 100 genes and 150 methylation regions that differentiated between complete responders and retreatment candidates. Hierarchical clustering clustered samples with the top 30 genes or the top 13 methylation regions. Dysregulated pathways included metabolic and DNA repair-related pathways. ASAP2 (log2fold = 0.29) and RAD50 (log2fold = 0.22) represented potential differentially expressed therapeutic targets.

CONCLUSIONS

Variation in gene expression and DNA methylation patterns in patients with HCC corresponded to ethnoracial status and LRT response. These initial results suggest tumor profiling has the potential to close ethnoracial disparities and improve treatment stratification.

Integrated Single-Cell Transcriptomic Atlas of Human Kidney Endothelial Cells

Key Points

Background

Kidney endothelial cells are exposed to different microenvironmental conditions that support specific physiologic processes. However, the heterogeneity of human kidney endothelial cells has not yet been systematically described.

Methods

We reprocessed and integrated seven human kidney control single-cell/single-nucleus RNA sequencing datasets of >200,000 kidney cells in the same process.

Results

We identified five major cell types, 29,992 of which were endothelial cells. Endothelial cell reclustering identified seven subgroups that differed in molecular characteristics and physiologic functions. Mapping new data to a normal kidney endothelial cell atlas allows rapid data annotation and analysis. We confirmed that endothelial cell types in the human kidney were also highly conserved in the mouse kidney and identified endothelial marker genes that were conserved in humans and mice, as well as differentially expressed genes between corresponding subpopulations. Furthermore, combined analysis of single-cell transcriptome data with public genome-wide association study data showed a significant enrichment of endothelial cells, especially arterial endothelial cells, in BP heritability. Finally, we identified M1 and M12 from coexpression networks in endothelial cells that may be deeply involved in BP regulation.

Conclusions

We created a comprehensive reference atlas of normal human kidney endothelial cells that provides the molecular foundation for understanding how the identity and function of kidney endothelial cells are altered in disease, aging, and between species. Finally, we provide a publicly accessible online tool to explore the datasets described in this work (https://vascularmap.jiangnan.edu.cn).

Unraveling Key Factors for Hypoxia Tolerance in Contrasting Varieties of Cotton Rose by Comparative Morpho-physiological and Transcriptome Analysis

The cotton rose (Hibiscus mutabilis) is a plant species commonly found in tropical and subtropical regions. It is remarkably resilient to waterlogging stress; however, the underlying mechanism behind this trait is yet unknown. This study used hypoxia-tolerant "Danbanhong" (DBH) and more hypoxia-sensitive "Yurui" (YR) genotypes and compared their morpho-physiological and transcriptional responses to hypoxic conditions. Notably, DBH had a higher number of adventitious roots (20.3) compared to YR (10.0), with longer adventitious roots in DBH (18.3 cm) than in YR (11.2 cm). Furthermore, the formation of aerenchyma was 3-fold greater in DBH compared to YR. Transcriptomic analysis revealed that DBH had more rapid transcriptional responses to hypoxia than YR. Identification of a greater number of differentially expressed genes (DEGs) for aerenchyma, adventitious root formation and development, and energy metabolism in DBH supported that DBH had better morphological and transcriptional adaptation than YR. DEG functional enrichment analysis indicated the involvement of variety-specific biological processes in adaption to hypoxia. Plant hormone signaling transduction, MAPK signaling pathway and carbon metabolism played more pronounced roles in DBH, whereas the ribosome genes were specifically induced in YR. These results show that effective multilevel coordination of adventitious root development and aerenchyma, in conjunction with plant hormone signaling and carbon metabolism, is required for increased hypoxia tolerance. This study provides new insights into the characterization of morpho-physiological and transcriptional responses to hypoxia in H. mutabilis, shedding light on the molecular mechanisms of its adaptation to hypoxic environments.

Comparative transcriptomic and metabolomic profiles reveal fruit peel color variation in two red pomegranate cultivars

Pomegranate (Punica granatum L.) which belongs to family Lythraceae, is one of the most important fruit crops of many tropical and subtropical regions. A high variability in fruit color is observed among different pomegranate accessions, which arises from the qualitative and quantitative differences in anthocyanins. However, the mechanism of fruit color variation is still not fully elucidated. In the present study, we investigated the red color mutation between a red-skinned pomegranate 'Hongbaoshi' and a purple-red-skinned cultivar 'Moshiliu', by using transcriptomic and metabolomic approaches. A total of 51 anthocyanins were identified from fruit peels, among which 3-glucoside and 3,5-diglucoside of cyanidin (Cy), delphinidin (Dp), and pelargonidin (Pg) were dominant. High proportion of Pg in early stages of 'Hongbaoshi' but high Dp in late stages of 'Moshiliu' were characterized. The unique high levels of Cy and Dp anthocyanins accumulating from early developmental stages accounted for the purple-red phenotype of 'Moshiliu'. Transcriptomic analysis revealed an early down-regulated and late up-regulated of anthocyanin-related structure genes in 'Moshiliu' compared with 'Hongbaoshi'. Alao, ANR was specially expressed in 'Hongbaoshi', with extremely low expression levels in 'Moshiliu'. For transcription factors R2R3-MYB, the profiles demonstrated a much higher transcription levels of three subgroup (SG) 5 MYBs and a sharp decrease in expression of SG6 MYB LOC116202527 in high-anthocyanin 'Moshiliu'. SG4 MYBs exhibited two entirely different patterns, LOC116203744 and LOC116212505 were down-regulated whereas LOC116205515 and LOC116212778 were up-regulated in 'Moshiliu' pomegranate. The results indicate that specific SG members of the MYB family might promote the peel coloration in different manners and play important roles in color mutation in pomegranate.

Transcriptomics of intracranial aneurysms: current state and opportunities in flow diversion

Over the last 2 decades, the field of transcriptomics has emerged as a major subdiscipline in biology. Transcriptomic techniques have been used by many groups over this time to better understand intracranial aneurysm development, rupture, and treatment. However, only a few studies have applied transcriptomics to understand the mechanisms behind flow diversion (FD) specifically, despite its increasing importance in the neurointerventional armamentarium. FD is an increasingly safe and effective treatment option for intracranial aneurysms. However, the clinical understanding and use of FD has far outpaced the understanding of the underlying mechanisms. To make FD more predictable, clinically efficacious, and safe, it is important to understand the biological mechanisms at play that lead to successful and unsuccessful FD. In this review, the authors focus on the current understanding of FD biology, the recent advances in transcriptomics, and what future studies could be performed to deepen the understanding of FD. They propose the new concept of the FD microenvironment to be studied, which may unlock a deeper biological understanding. This review provides the background for prospective studies into the development of targeted aneurysm therapy, whether by modified devices or by medical adjuncts.

Investigation of heat stress responses and adaptation mechanisms by integrative metabolome and transcriptome analysis in tea plants (Camellia sinensis)

Extreme high temperature has deleterious impact on the yield and quality of tea production, which has aroused the attention of growers and breeders. However, the mechanisms by which tea plant varieties respond to extreme environmental heat is not clear. In this study, we analyzed physiological indices, metabolites and transcriptome differences in three different heat-tolerant tea plant F1 hybrid progenies. Results showed that the antioxidant enzyme activity, proline, and malondialdehyde were significantly decreased in heat-sensitive 'FWS' variety, and the accumulation of reactive oxygen molecules such as H2O2 and O2- was remarkably increased during heat stress. Metabolomic analysis was used to investigate the metabolite accumulation pattern of different varieties in response to heat stress. The result showed that a total of 810 metabolites were identified and more than 300 metabolites were differentially accumulated. Transcriptional profiling of three tea varieties found that such genes encoding proteins with chaperon domains were preferentially expressed in heat-tolerant varieties under heat stress, including universal stress protein (USP32, USP-like), chaperonin-like protein 2 (CLP2), small heat shock protein (HSP18.1), and late embryogenesis abundant protein (LEA5). Combining metabolomic with transcriptomic analyses discovered that the flavonoids biosynthesis pathway was affected by heat stress and most flavonols were up-regulated in heat-tolerant varieties, which owe to the preferential expression of key FLS genes controlling flavonol biosynthesis. Take together, molecular chaperons, or chaperon-like proteins, flavonols accumulation collaboratively contributed to the heat stress adaptation in tea plant. The present study elucidated the differences in metabolite accumulation and gene expression patterns among three different heat-tolerant tea varieties under extreme ambient high temperatures, which helps to reveal the regulatory mechanisms of tea plant adaptation to heat stress, and provides a reference for the breeding of heat-tolerant tea plant varieties.

Metabolomics of Multiple Sclerosis Lesions Demonstrates Lipid Changes Linked to Alterations in Transcriptomics-Based Cellular Profiles

BACKGROUND AND OBJECTIVES

People with multiple sclerosis (MS) have a dysregulated circulating metabolome, but the metabolome of MS brain lesions has not been studied. The aims of this study were to identify differences in the brain tissue metabolome in MS compared with controls and to assess its association with the cellular profile of corresponding tissue.

METHODS

MS tissues included samples from the edge and core of chronic active or inactive lesions and periplaque white matter (WM). Control specimens were obtained from normal WM. Metabolomic analysis was performed using mass-spectrometry coupled with liquid/gas chromatography and subsequently integrated with single-nucleus RNA-sequencing data by correlating metabolite abundances with relative cell counts, as well as individual genes using Multiomics Factor Analysis (MOFA).

RESULTS

Seventeen samples from 5 people with secondary progressive MS and 8 samples from 6 controls underwent metabolomic profiling identifying 783 metabolites. MS lesions had higher levels of sphingosines (false discovery rate-adjusted p-value[q] = 2.88E-05) and sphingomyelins and ceramides (q = 2.15E-07), but lower nucleotide (q = 0.05), energy (q = 0.001), lysophospholipid (q = 1.86E-07), and monoacylglycerol (q = 0.04) metabolite levels compared with control WM. Periplaque WM had elevated sphingomyelins and ceramides (q = 0.05) and decreased energy metabolites (q = 0.01) and lysophospholipids (q = 0.05) compared with control WM. Sphingolipids and membrane lipid metabolites were positively correlated with astrocyte and immune cell abundances and negatively correlated with oligodendrocytes. On the other hand, long-chain fatty acid, endocannabinoid, and monoacylglycerol pathways were negatively correlated with astrocyte and immune cell populations and positively correlated with oligodendrocytes. MOFA demonstrated associations between differentially expressed metabolites and genes involved in myelination and lipid biosynthesis.

DISCUSSION

MS lesions and perilesional WM demonstrated a significantly altered metabolome compared with control WM. Many of the altered metabolites were associated with altered cellular composition and gene expression, indicating an important role of lipid metabolism in chronic neuroinflammation in MS.

Integrative bioinformatics and experimental validation of hub genetic markers in acne vulgaris: Toward personalized diagnostic and therapeutic strategies

BACKGROUND

Acne vulgaris is a widespread chronic inflammatory dermatological condition. The precise molecular and genetic mechanisms of its pathogenesis remain incompletely understood. This research synthesizes existing databases, targeting a comprehensive exploration of core genetic markers.

METHODS

Gene expression datasets (GSE6475, GSE108110, and GSE53795) were retrieved from the GEO. Differentially expressed genes (DEGs) were identified using the limma package. Enrichment analyses were conducted using GSVA for pathway assessment and clusterProfiler for GO and KEGG analyses. PPI networks and immune cell infiltration were analyzed using the STRING database and ssGSEA, respectively. We investigated the correlation between hub gene biomarkers and immune cell infiltration using Spearman's rank analysis. ROC curve analysis validated the hub genes' diagnostic accuracy. miRNet, TarBase v8.0, and ChEA3 identified miRNA/transcription factor-gene interactions, while DrugBank delineated drug-gene interactions. Experiments utilized HaCaT cells stimulated with Propionibacterium acnes, treated with retinoic acid and methotrexate, and evaluated using RT-qPCR, ELISA, western blot, lentiviral transduction, CCK-8, wound-healing, and transwell assays.

RESULTS

There were 104 genes with consistent differences across the three datasets of paired acne and normal skin. Functional analyses emphasized the significant enrichment of these DEGs in immune-related pathways. PPI network analysis pinpointed hub genes PTPRC, CXCL8, ITGB2, and MMP9 as central players in acne pathogenesis. Elevated levels of specific immune cell infiltration in acne lesions corroborated the inflammatory nature of the disease. ROC curve analysis identified the acne diagnostic potential of four hub genes. Key miRNAs, particularly hsa-mir-124-3p, and central transcription factors like TFEC were noted as significant regulators. In vitro validation using HaCaT cells confirmed the upregulation of hub genes following Propionibacterium acnes exposure, while CXCL8 knockdown reduced pro-inflammatory cytokines, cell proliferation, and migration. DrugBank insights led to the exploration of retinoic acid and methotrexate, both of which mitigated gene expression upsurge and inflammatory mediator secretion.

CONCLUSION

This comprehensive study elucidated pivotal genes associated with acne pathogenesis, notably PTPRC, CXCL8, ITGB2, and MMP9. The findings underscore potential biomarkers, therapeutic targets, and the therapeutic potential of agents like retinoic acid and methotrexate. The congruence between bioinformatics and experimental validations suggests promising avenues for personalized acne treatments.

Metatranscriptomics sheds light on the links between the functional traits of fungal guilds and ecological processes in forest soil ecosystems

Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. We showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.

Identification of Potential Neddylation-related Key Genes in Ischemic Stroke based on Machine Learning Methods

Ischemic stroke (IS) is a complex neurological disease that can lead to severe disability or even death. Understanding the molecular mechanisms involved in the occurrence and progression of IS is of great significance for developing effective treatment strategies. In this context, the role of neddylation refers to the potential impact of neddylation on various cellular processes, which may contribute to the pathogenesis and outcome of IS. First, differential analysis was conducted on the GSE16561 dataset from the GEO database to identify 350 differentially expressed genes (DEGs) between the IS and Control groups. By intersecting the differential genes with neddylation-related genes, 11 neddylation-related DEGs were obtained. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses showed that the DEGs were mainly enriched in hematopoietic cell lineage and neutrophil degranulation, while the neddylation-related DEGs were mainly enriched in apoptosis and post-translational protein modification. Further Lasso-Cox and random forest analyses were performed on the 11 neddylation-related DEGs, identifying key genes SRPK1, BIRC2, and KLHL3. Additionally, validation of the key genes was carried out using the GSE58294 dataset and clinical patients. Finally, the correlation between the key genes and ferroptosis and cuproptosis was analyzed, and a ceRNA network was constructed. Our study helps to elucidate the complex role of neddylation in the mechanism of ischemic stroke, providing potential opportunities for the development of therapeutic interventions.

Transcriptome meta-analysis-based identification of hub transcription factors and RNA-binding proteins potentially orchestrating gene regulatory cascades and crosstalk in response to abiotic stresses in Arabidopsis thaliana

Deteriorating climatic conditions and increasing human population necessitate the development of robust plant varieties resistant to harsh environments. Manipulation of regulatory proteins such as transcription factors (TFs) and RNA-binding proteins (RBPs) would be a beneficial strategy in this regard. Further, understanding the complex interconnections between different classes of regulatory molecules would be essential for the identification of candidate genes/proteins for trait improvement. Most studies to date have analysed the roles of TFs or RBPs individually, in conferring stress resilience. However, it would be important to identify dominant/upstream TFs and RBPs inducing widespread transcriptomic alterations through other regulators (i.e., other TFs/RBPs targeted by the upstream regulators). To this end, the present study employed a transcriptome meta-analysis and computational approaches to obtain a comprehensive overview of regulatory interactions. This work identified dominant TFs and RBPs potentially influencing stress-mediated differential expression of other regulators, which could in turn influence gene expression, and consequently, physiological responses. Twenty transcriptomic studies [related to (i) UV radiation, (ii) wounding, (iii) salinity, (iv) cold, and (v) drought stresses in Arabidopsis thaliana] were analysed for differential gene expression, followed by the identification of differentially expressed TFs and RBPs. Subsequently, other TFs and RBPs which could be influencing these regulators were identified, and their interaction networks and hub nodes were analysed. As a result, an interacting module of Basic Leucine Zipper (bZIP) family TFs as well as Heterogeneous nuclear ribonucleoproteins (hnRNP) and Glycine-rich protein (GRP) family RBPs (among other TFs and RBPs) were shown to potentially influence the stress-induced differential expression of other TFs and RBPs under all the considered stress conditions. Some of the identified hub TFs and RBPs are known to be of major importance in orchestrating stress-induced transcriptomic changes influencing a variety of physiological processes from seed germination to senescence. This study highlighted the gene/protein candidates that could be considered for multiplexed genetic manipulation - a promising approach to develop robust, multi-stress-resilient plant varieties.

Integrated transcriptomic and CGAs analysis revealed IbGLK1 is a key transcription factor for chlorogenic acid accumulation in sweetpotato (Ipomoea batatas [L.] Lam.) blades

Sweetpotato blades are rich in the functional secondary metabolite chlorogenic acid (CGA), which deepen potential for effective utilization of the blade in industry. In this study, we evaluated the type and content of CGA in the blades of 16 sweetpotato genotypes and analyzed the correlation between CGA content and antioxidant capacity. Then we isolated and characterized IbGLK1, a GARP-type transcription factor, by comparative transcriptome analysis. A subcellular localization assay indicated that IbGLK1 is located in the nucleus. Overexpression and silencing of IbGLK1 in sweetpotato blade resulted in a 0.90-fold increase and 1.84-fold decrease, respectively, in CGA content compared to the control. Yeast one-hybrid and dual-luciferase assays showed that IbGLK1 binds and activates the promoters of IbHCT, IbHQT, IbC4H, and IbUGCT, resulting in the promotion of CGA biosynthesis. In conclusion, our study provides insights into a high-quality gene for the regulation of CGA metabolism and germplasm resources for breeding sweetpotato.

Proteomic and Transcriptomic Analyses Provide New Insights into the Mechanism Underlying Lipid Deterioration in Pecan Kernels during Storage

Pecan nuts are rich in lipids that tend to deteriorate during storage. Tandem mass-tag-based quantitative proteomics and transcriptomics were used to investigate the changes in the protein and gene profiles of stored pecan kernels for the first time. Our previous lipidomic data were jointly analyzed to elucidate the coordinated changes in lipid molecules and related proteins/genes. The mechanism underlying lipid deterioration in pecan kernels during storage was revealed by multiomics analyses. Lipid metabolism-related pathways were activated during pecan storage. Phospholipases, triacylglycerol lipases, lipoxygenases, and oil body-related proteins/genes were highly expressed during storage, revealing their involvement in lipid deterioration. These data provide rich information and will be valuable for future genetic or chemical research to alleviate lipid deterioration in pecans.

Comparative transcriptome analysis provides insights into the resistance regulation mechanism and inhibitory effect of fungicide phenamacril in Fusarium asiaticum

Fusarium asiaticum is a destructive phytopathogenic fungus that causes Fusarium head blight of wheat (FHB), leading to serious yield and economic losses to cereal crops worldwide. Our previous studies indicated that target-site mutations (K216R/E, S217P/L, or E420K/G/D) of Type I myosin FaMyo5 conferred high resistance to phenamacril. Here, we first constructed one sensitive strain H1S and three point mutation resistant strains HA, HC and H1R. Then we conducted comparative transcriptome analysis of these F. asiaticum strains after 1 and 10 μg·mL-1 phenamacril treatment. Results indicated that 2135 genes were differentially expressed (DEGs) among the sensitive and resistant strains. The DEGs encoding ammonium transporter MEP1/MEP2, nitrate reductase, copper amine oxidase 1, 4-aminobutyrate aminotransferase, amino-acid permease inda1, succinate-semialdehyde dehydrogenase, 2, 3-dihydroxybenzoic acid decarboxylase, etc., were significantly up-regulated in all the phenamacril-resistant strains. Compared to the control group, a total of 1778 and 2097 DEGs were identified in these strains after 1 and 10 μg·mL-1 phenamacril treatment, respectively. These DEGs involved in 4-aminobutyrate aminotransferase, chitin synthase 1, multiprotein-bridging factor 1, transcriptional regulatory protein pro-1, amino-acid permease inda1, ATP-dependent RNA helicase DED1, acetyl-coenzyme A synthetase, sarcoplasmic/endoplasmic reticulum calcium ATPase 2, etc., showed significantly down-regulated expression in phenamacril-sensitive strain but not in resistant strains after phenamacril treatment. In addition, cyanide hydratase, mating-type protein MAT-1, putative purine nucleoside permease, plasma membrane protein yro2, etc., showed significantly co-down-regulated expression in all the strains after phenamacril treatment. Taken together, This study provides deep insights into the resistance regulation mechanism and the inhibitory effect of fungicide phenamacril and these new annotated proteins or enzymes are worth for the discovery of new fungicide targets.

Multi-omics analysis reveals the association between specific solute carrier proteins gene expression patterns and the immune suppressive microenvironment in glioma

Glioma is the most prevalent malignant brain tumour. Currently, reshaping its tumour microenvironment has emerged as an appealing strategy to enhance therapeutic efficacy. As the largest group of transmembrane transport proteins, solute carrier proteins (SLCs) are responsible for the transmembrane transport of various metabolites and ions. They play a crucial role in regulating the metabolism and functions of malignant cells and immune cells within the tumour microenvironment, making them a promising target in cancer therapy. Through multidimensional data analysis and experimental validation, we investigated the genetic landscape of SLCs in glioma. We established a classification system comprising 7-SLCs to predict the prognosis of glioma patients and their potential responses to immunotherapy and chemotherapy. Our findings unveiled specific SLC expression patterns and their correlation with the immune-suppressive microenvironment and metabolic status. The 7-SLC classification system was validated in distinguishing subgroups within the microenvironment, specifically identifying subsets involving malignant cells and tumour-associated macrophages. Furthermore, the orphan protein SLC43A3, a core member of the 7-SLC classification system, was identified as a key facilitator of tumour cell proliferation and migration, suggesting its potential as a novel target for cancer therapy.

Th17 Cell-Related Gene Biomarkers in Osteoporosis: Comprehensive Bioinformatics Analysis and In Vivo Validation

The interaction between the bone and immune systems has a major role in osteoporosis regulation. However, the infiltration of T helper 17 (Th17) cells and their associated genes in osteoporosis remains unclear. The GSE35959 dataset was obtained from the Gene Expression Omnibus (GEO) database, and the Immune Cell Abundance Identifier (ImmuCellAI) program was used to evaluate the abundance of 24 immune cell types, including Th17 cells. Differential analysis and relevance analysis were performed to identify differentially expressed Th17 cell-related genes (DETh17RGs) in osteoporosis. The potential functions of DETh17RGs were analyzed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment. Hub DETh17RGs were obtained through comprehensive analysis using Weighted Gene Co-Expression Network Analysis (WGCNA) and the CytoHubba plug-in algorithm. The expression levels of hub genes were validated using additional osteoporosis datasets. Additionally, the transcript levels of Hub genes in a mice model of osteoporosis were examined using quantitative PCR (qPCR). 464 DETh17RGs were identified in this study, with 421 genes showing positive associations and 43 genes showing negative associations. Among these, seven genes (CD44, TGFB1, ACTN4, ARHGDIA, ESR1, TLN1, FLNA) were considered as Hub DETh17RGs. The qPCR transcript levels of hub DETh17RGs in a mice model of osteoporosis exhibited consistent expression trends with the bioinformatics analysis. This research enhances our understanding of the molecular mechanisms involving Th17 cells in the development of osteoporosis and contributes to the discovery of potential biomarkers.

Genome-wide analysis of early vascular tunic repair and regeneration for Botrylloides digenesis reveals striking similarities to human wound healing

The early stages of regeneration after injury are similar to those of wound healing. The ascidian Botrylloides diegensis can regenerate an entire adult from a small fragment of vascular tunic following the removal of all zooids in an injury-induced regeneration model. We investigated the molecular and cellular changes following injury to determine the differences between the healing process and the initiation of whole-body regeneration (WBR). We conducted transcriptome analysis at specific time points during regeneration and wound healing to identify differentially expressed genes (DEGs) and the unique biological processes associated with each state. Our findings revealed 296 DEGs at 10 h post-injury (hpi), with 71 highly expressed in healed tissue and 225 expressed during the WBR process. These DEGs were predicted to play roles in tissue reorganization, integrin signaling, extracellular matrix organization, and the innate immune system. Pathway analysis of the upregulated genes in the healed tunic indicated functional enrichment related to tissue repair, as has been observed in other species. Additionally, we examined the cell types in the tunic and ampullae in both tissue states using histology and in situ hybridization for six genes identified by transcriptome analysis. We observed strong mRNA expression in cells within the WBR tunic, and in small RNA-positive granules near the tunic edge. We hypothesized that many of these genes function in the compaction of the ampullae tunic, which is a pivotal process for WBR and dormancy in B. diegensis, and in an immune response. These findings establish surprising similarities between ascidian regeneration and human wound healing, emphasizing the potential for future investigations into human regenerative and repair mechanisms. This study provides valuable insights into the gene sets specifically activated during regeneration compared to wound healing, shedding light on the divergent activities of these processes.

Combined transcriptomics and metabolomics to analyse the response of Cuminum cyminum L. under Pb stress

Lead (Pb) can disrupt plant gene expression, modify metabolite contents, and influence the growth of plants. Cuminum cyminum L. is highly adaptable to adversity, but molecular mechanism by which it responds to Pb stress is unknown. For this study, transcriptomic and metabolomic sequencing was performed on root tissues of C. cyminum under Pb stress. Our results showed that high Pb stress increased the activity of peroxidase (POD), the contents of malondialdehyde (MDA) and proline by 80.03 %, 174.46 % and 71.24 %, respectively. Meanwhile, Pb stress decreased the activities of superoxide dismutase (SOD) and catalase (CAT) as well as contents of soluble sugars and GSH, which thus affected the growth of C. cyminum. In addition, Pb stress influenced the accumulation and transport of Pb in C. cyminum. Metabolomic results showed that Pb stress affected eight metabolic pathways involving 108 differentially expressed metabolites, primarily amino acids, organic acids, and carbohydrates. The differentially expressed genes identified through transcriptome analysis were mainly involved the oxidation reductase activity, transmembrane transport, phytohormone signaling, and MAPK signaling pathway. The results of this study will help to understand the molecular mechanisms of C. cyminum response to Pb stress, and provide a basis for screening seeds with strong resistance to heavy metals.

Identifying transcriptomic indicators of tertiary treated municipal effluent in longnose dace (Rhinichthys cataractae) caged under semi-controlled conditions in experimental raceways

To evaluate effects of tertiary treated wastewater treatment plant effluent (MWWE) on transcriptomic responses in longnose dace (Rhinichthys cataractae; LND) we conducted a semi-controlled study in experimental raceways (Advancing Canadian Water Assets facility) imbedded in the Pine Creek treatment plant (Calgary, AB). LND collected from a reference site in the Bow River (REF) were caged in raceways containing either 5 % Pine Creek effluent (PC) or Bow River water (BR; control) over 28 d. Liver transcriptomes were analyzed in males and females sampled on days 7, 14 and 28 from BR and PC, and compared to REF fish on day 0. Concurrent with the caging, selected environmental substances of concern were analyzed in the BR and PC. Significantly different unigenes (SDUs) in females (vs males) within both BR and PC raceways increased over time and compared to REF fish. Moreover, SDUs in females and males within the same treatment (i.e., BR, PC) showed a temporal increase as well as compared to REF fish. Time was the dominant factor affecting SDUs, whereas sex and treatment had less of an impact on the transcriptome profiling. Gene Set Enrichment Analysis of BR vs PC over time revealed effects on genes involved in growth, metabolism of carbohydrates and lipids, and immune system on day 7; however, by day 28, 80-100 % of the transcripts localized to enriched biomarkers were associated with tissue immune responses in both sexes. Exposure to 5 % effluent had significant effects on female liver somatic index but no effects were observed on other phenotypic health indices in either sex. BR was used as the source of reference water, but analyses showed trace amounts of ESOCs. Analyses did not point towards definitive response patterns that could be used in field-based ecotoxicogenomic studies on the impacts of well-treated MWWE but suggested compromised adaptive immune responses.

Transcriptome analysis of primary adult B-cell lineage acute lymphoblastic leukemia identifies pathogenic variants and gene fusions, and predicts subtypes for in depth molecular diagnosis

BACKGROUND

B-cell acute lymphoblastic leukemia (B-ALL) is classified into subgroups based on known driver oncogenes and molecular lesions, including translocations and recurrent mutations. However, the current diagnostic tests do not identify subtypes or oncogenic lesions for all B-ALL samples, creating a heterogeneous B-ALL group of unknown subtypes.

METHODS

We sorted primary adult B-ALL cells and performed transcriptome analysis by bulk RNA sequencing (RNA-seq).

RESULTS

Transcriptomic analysis of an adult B-ALL cohort allowed the classification of four patient samples with subtypes that were not previously revealed by standard gene panels. The leukemia of two patients were of the DUX4 subtype and two were CRLF2+ Ph-like B-ALL. Furthermore, single nucleotide variant analysis detected the oncogenic NRAS-G12D, KRAS-G12D, and KRAS-G13D mutations in three of the patient samples, presenting targetable mutations. Additional oncogenic variants and gene fusions were uncovered, as well as multiple variants in the PDE4DIP gene across five of the patient samples.

CONCLUSION

We demonstrate that RNA-seq is an effective tool for precision medicine in B-ALL by providing comprehensive molecular profiling of leukemia cells, identifying subtype and oncogenic lesions, and stratifying patients for appropriate therapy.

The transcription factor OsbZIP48 governs rice responses to zinc deficiency

Zinc (Zn) deficiency is the most prevalent micronutrient disorder in rice and leads to delayed development and decreased yield. Nevertheless, despite its primary importance, how rice responds to Zn deficiency remains poorly understood. This study presents genetic evidence supporting the crucial role of OsbZIP48 in regulating rice's response to Zn deficiency, consistent with earlier findings in the model plant Arabidopsis. Genetic inactivation of OsbZIP48 in rice seedlings resulted in heightened sensitivity to Zn deficiency and reduced Zn translocation from roots to shoots. Consistently, OsbZIP48 was constitutively expressed in roots, slightly induced by Zn deficiency in shoots and localized into nuclei induced by Zn deficiency. Comparative transcriptome analysis of the wild-type plants and osbzip48 mutant grown under Zn deficiency enabled the identification of OsbZIP48 target genes, including key Zn transporter genes (OsZIP4 and OsZIP8). We demonstrated that OsbZIP48 controlled the expressions of these genes by directly binding to their promoters, specifically to the Zn deficiency response element motif. This study establishes OsbZIP48 as a critical transcription factor in rice's response to Zn deficiency, offering valuable insights for developing Zn-biofortified rice varieties to combat global Zn limitation.

Novel insights into drought-induced regulation of ribosomal genes through DNA methylation in chickpea

Modifications within the epigenome of an organism in response to external environmental conditions allow it to withstand the hostile stress factors. Drought in chickpea is a severely limiting abiotic stress factor which is known to cause huge yield loss. To analyse the methylome of chickpea in response to drought stress conditions and how it affects gene expression, we performed whole-genome bisulfite sequencing (WGBS) and RNA-seq of two chickpea genotypes which contrast for drought tolerance. It was observed that the mCHH was most variable under drought stress and the drought tolerant (DT) genotype exhibited substantial genome-wide hypomethylation as compared to the drought sensitive (DS) genotype. Specifically, there was substantial difference in gene expression and methylation for the ribosomal genes for the tolerant and sensitive genotypes. The differential expression of these genes was in complete agreement with earlier reported transcriptomes in chickpea. Many of these genes were hypomethylated (q < 0.01) and downregulated under drought stress (p < 0.01) in the sensitive genotype. The gene RPS6 (ribosomal protein small subunit) was found to be downregulated and hypomethylated in the drought sensitive genotype which could possibly lead to reduced ribosomal biosynthesis. This study provides novel insights into regulation of drought-responsive genes in chickpea.

Atrial proteomic profiling reveals a switch towards profibrotic gene expression program in CREM-IbΔC-X mice with persistent atrial fibrillation

BACKGROUND

Overexpression of the CREM (cAMP response element-binding modulator) isoform CREM-IbΔC-X in transgenic mice (CREM-Tg) causes the age-dependent development of spontaneous AF.

PURPOSE

To identify key proteome signatures and biological processes accompanying the development of persistent AF through integrated proteomics and bioinformatics analysis.

METHODS

Atrial tissue samples from three CREM-Tg mice and three wild-type littermates were subjected to unbiased mass spectrometry-based quantitative proteomics, differential expression and pathway enrichment analysis, and protein-protein interaction (PPI) network analysis.

RESULTS

A total of 98 differentially expressed proteins were identified. Gene ontology analysis revealed enrichment for biological processes regulating actin cytoskeleton organization and extracellular matrix (ECM) dynamics. Changes in ITGAV, FBLN5, and LCP1 were identified as being relevant to atrial fibrosis and structural based on expression changes, co-expression patterns, and PPI network analysis. Comparative analysis with previously published datasets revealed a shift in protein expression patterns from ion-channel and metabolic regulators in young CREM-Tg mice to profibrotic remodeling factors in older CREM-Tg mice. Furthermore, older CREM-Tg mice exhibited protein expression patterns reminiscent of those seen in humans with persistent AF.

CONCLUSIONS

This study uncovered distinct temporal changes in atrial protein expression patterns with age in CREM-Tg mice consistent with the progressive evolution of AF. Future studies into the role of the key differentially abundant proteins identified in this study in AF progression may open new therapeutic avenues to control atrial fibrosis and substrate development in AF.

Single cell sequencing data identify distinct B cell and fibroblast populations in stricturing Crohn's disease

Single cell RNA sequencing of human full thickness Crohn's disease (CD) small bowel resection specimens was used to identify potential therapeutic targets for stricturing (S) CD. Using an unbiased approach, 16 cell lineages were assigned within 14,539 sequenced cells from patient-matched SCD and non-stricturing (NSCD) preparations. SCD and NSCD contained identical cell types. Amongst immune cells, B cells and plasma cells were selectively increased in SCD samples. B cell subsets suggested formation of tertiary lymphoid tissue in SCD and compared with NSCD there was an increase in IgG, and a decrease in IgA plasma cells, consistent with their potential role in CD fibrosis. Two Lumican-positive fibroblast subtypes were identified and subclassified based on expression of selectively enriched genes as fibroblast clusters (C) 12 and C9. Cells within these clusters expressed the profibrotic genes Decorin (C12) and JUN (C9). C9 cells expressed ACTA2; ECM genes COL4A1, COL4A2, COL15A1, COL6A3, COL18A1 and ADAMDEC1; LAMB1 and GREM1. GO and KEGG Biological terms showed extracellular matrix and stricture organization associated with C12 and C9, and regulation of WNT pathway genes with C9. Trajectory and differential gene analysis of C12 and C9 identified four sub-clusters. Intra sub-cluster gene analysis detected 13 co-regulated gene modules that aligned along predicted pseudotime trajectories. CXCL14 and ADAMDEC1 were key markers in module 1. Our findings support further investigation of fibroblast heterogeneity and interactions with local and circulating immune cells at earlier time points in fibrosis progression. Breaking these interactions by targeting one or other population may improve therapeutic management for SCD.

Profiling the activity of the para-caspase MALT1 in B-cell acute lymphoblastic leukemia for potential targeted therapeutic application

B-cell acute lymphoblastic leukemia (B-ALL) remains a hard-to-treat disease with a poor prognosis in adults. Mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) is a para-caspase required for B-cell receptor (BCR)-mediated NF-κB activation. Inhibition of MALT1 in preclinical models has proven efficacious in many B-cell malignancies including chronic lymphocytic leukemia, mantle cell lymphoma and diffuse large B-cell lymphoma. We sought to examine the role of MALT1 in B-ALL and determine the biological consequences of its inhibition. Targeting MALT1 with both Z-VRPR-fmk and MI-2 efficiently kills B-ALL cells independent of the cell-of-origin (pro, pre, mature) or the presence of the Philadelphia chromosome, and spares normal B cells. The mechanism of cell death was through apoptotic induction, mostly in cycling cells. The proteolytic activity of MALT1 can be studied by measuring its ability to cleave its substrates. Surprisingly, with the exception of mature B-ALL, we did not detect cleavage of MALT1 substrates at baseline, nor after proteasomal inhibition or following activation of pre-BCR. To explore the possibility of a distinct role for MALT1 in B-ALL, independent of signaling through BCR, we studied the changes in gene expression profiling following a 24-hour treatment with MI-2 in 12 B-ALL cell lines. Our transcriptome analysis revealed a strong inhibitory effect on MYC-regulated gene signatures, further confirmed by Myc protein downregulation, concomitant with an increase in the Myc degrader FBXW7. In conclusion, our evidence suggests a novel role for MALT1 in B-ALL through Myc regulation and provides support for clinical testing of MALT1 inhibitors in B-ALL.

Identification and clinicopathological analysis of potential p73-regulated biomarkers in colorectal cancer via integrative bioinformatics

This study aims to decipher crucial biomarkers regulated by p73 for the early detection of colorectal cancer (CRC) by employing a combination of integrative bioinformatics and expression profiling techniques. The transcriptome profile of HCT116 cell line p53- / - p73+ / + and p53- / - p73 knockdown was performed to identify differentially expressed genes (DEGs). This was corroborated with three CRC tissue expression datasets available in Gene Expression Omnibus. Further analysis involved KEGG and Gene ontology to elucidate the functional roles of DEGs. The protein-protein interaction (PPI) network was constructed using Cytoscape to identify hub genes. Kaplan-Meier (KM) plots along with GEPIA and UALCAN database analysis provided the insights into the prognostic and diagnostic significance of these hub genes. Machine/deep learning algorithms were employed to perform TNM-stage classification. Transcriptome profiling revealed 1289 upregulated and 1897 downregulated genes. When intersected with employed CRC datasets, 284 DEGs were obtained. Comprehensive analysis using gene ontology and KEGG revealed enrichment of the DEGs in metabolic process, fatty acid biosynthesis, etc. The PPI network constructed using these 284 genes assisted in identifying 20 hub genes. Kaplan-Meier, GEPIA, and UALCAN analyses uncovered the clinicopathological relevance of these hub genes. Conclusively, the deep learning model achieved TNM-stage classification accuracy of 0.78 and 0.75 using 284 DEGs and 20 hub genes, respectively. The study represents a pioneer endeavor amalgamating transcriptomics, publicly available tissue datasets, and machine learning to unveil key CRC-associated genes. These genes are found relevant regarding the patients' prognosis and diagnosis. The unveiled biomarkers exhibit robustness in TNM-stage prediction, thereby laying the foundation for future clinical applications and therapeutic interventions in CRC management.

Identification and characterization of PAL genes involved in the regulation of stem development in Saccharum spontaneum L

BACKGROUND

Saccharum spontaneum L. is a closely related species of sugarcane and has become an important genetic component of modern sugarcane cultivars. Stem development is one of the important factors for affecting the yield, while the molecular mechanism of stem development remains poorly understanding in S. spontaneum. Phenylalanine ammonia-lyase (PAL) is a vital component of both primary and secondary metabolism, contributing significantly to plant growth, development and stress defense. However, the current knowledge about PAL genes in S. spontaneum is still limited. Thus, identification and characterization of the PAL genes by transcriptome analysis will provide a theoretical basis for further investigation of the function of PAL gene in sugarcane.

RESULTS

In this study, 42 of PAL genes were identified, including 26 SsPAL genes from S. spontaneum, 8 ShPAL genes from sugarcane cultivar R570, and 8 SbPAL genes from sorghum. Phylogenetic analysis showed that SsPAL genes were divided into three groups, potentially influenced by long-term natural selection. Notably, 20 SsPAL genes were existed on chromosomes 4 and 5, indicating that they are highly conserved in S. spontaneum. This conservation is likely a result of the prevalence of whole-genome replications within this gene family. The upstream sequence of PAL genes were found to contain conserved cis-acting elements such as G-box and SP1, GT1-motif and CAT-box, which collectively regulate the growth and development of S. spontaneum. Furthermore, quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that SsPAL genes of stem had a significantly upregulated than that of leaves, suggesting that they may promote the stem growth and development, particularly in the + 6 stem (The sixth cane stalk from the top to down) during the growth stage.

CONCLUSIONS

The results of this study revealed the molecular characteristics of SsPAL genes and indicated that they may play a vital role in stem growth and development of S. spontaneum. Altogether, our findings will promote the understanding of the molecular mechanism of S. spontaneum stem development, and also contribute to the sugarcane genetic improving.

Identification of core genes involved in the response of Apocynum venetum to salt stress based on transcriptome sequencing and WGCNA

Apocynum venetum L. belongs to the Apocynaceae family and is a plant that is highly resistant to stress. It is important in the fields of ecology, feeding, industry and medicine. The molecular mechanism underlying salt tolerance has not been elucidated. In this study, RNA-seq based transcriptome sequencing of A. venetum leaves after 0, 2, 6, 12, 24 and 48 h of treatment with 300 mM NaCl was performed. We conducted a comprehensive analysis of the transcriptome expression profiles of A. venetum under salt stress using the WGCNA method and identified red, black, and brown as the core modules regulating the salt tolerance of A. venetum. A co-expression regulatory network was constructed to identify the core genes in the module according to the correlations between genes. The genes TRINITY_DN102_c0_g1 (serine carboxypeptidase), TRINITY_DN3073_c0_g1 (SOS signaling pathway) and TRINITY_DN6732_c0_g1 (heat shock transcription factor) in the red module were determined to be the core genes. Two core genes in the black module, TRINITY_DN9926_c0_g1 and TRINITY_DN7962_c0_g1, are pioneer candidate salt tolerance-associated genes in A. venetum. The genes in the brown module were mainly enriched in two pathways, namely photosynthesis and osmotic balance. Among them, the TRINITY_DN6321_c0_g2 and TRINITY_DN244_c0_g1 genes encode aquaporin, which is helpful for maintaining the cell water balance and plays a protective role in defending A. venetum under abiotic stress. Our findings contribute to the identification of core genes involved in the response of A. venetum to salt stress.

Differential expression and global analysis of miR156/SQUAMOSA promoter binding-like proteins (SPL) module in oat

SQUAMOSA promoter binding-like proteins (SPLs) are important transcription factors that influence growth phase transition and reproduction in plants. SPLs are targeted by miR156 but the SPL/miR156 module is completely unknown in oat. We identified 28 oat SPL genes (AsSPLs) distributed across all 21 oat chromosomes except for 4C and 6D. The oat- SPL gene family represented six of eight SPL phylogenetic groups, with no AsSPLs in groups 3 and 7. A novel oat miR156 (AsmiR156) family with 21 precursors divided into 7 groups was characterized. A total of 16 AsSPLs were found to be targeted by AsmiR156. Intriguingly, AsSPL3s showed high transcript abundance during early inflorescence (GS-54), as compared to the lower abundance of AsmiR156, indicating their role in reproductive development. Unravelling the SPL/miR156 regulatory hub and alterations in expression patterns of AsSPLs could provide an essential toolbox for genetic improvement in the cultivated oat.

Generation of selenium-rich wheat mutants and exploration of responsive genes for selenium accumulation

KEY MESSAGE

Salt tolerance, selenium accumulation and expression of the responsive genes were analyzed in the wheat high selenium mutants. Selenium is an essential trace element for the human body, and its deficiency can lead to various diseases such as Keshan disease and large bone disease. Wheat, being a major staple crop, plays a crucial role in providing dietary selenium supplementation to combat this deficiency. Despite progress in understanding the molecular regulation of selenium accumulation in certain crops, the molecular mechanisms governing selenium accumulation-related gene expression in wheat plants remain poorly understood. In this study, three mutant wheat lines with elevated selenium content were identified. Under the treatment of Na2SeO3 or NaCl, the selenium-rich wheat mutants exhibited decreased sensitivity to both selenium and NaCl compared to the wild type. Additionally, there was an increase in the activities of SOD and POD, while the content of MDA decreased. Through qRT-PCR analysis, the expression of selenium-related genes was affected, revealing that some of these genes not only regulate the response of wheat to salt stress, but also play a role in the process of selenium accumulation. The transcriptome results revealed that the important genes encoding glutathione S-transferases, peroxidases, superoxide dismutases, and UDP-glucosyltransferases may function in the regulation of salt tolerance and selenium accumulation in wheat. These findings significantly contribute to the current understanding of the molecular regulation of selenium accumulation in wheat crops, while also offering novel germplasm resources for cultivating selenium-rich and salt-tolerant wheat lines.

Effects of trimer repeats on Psidium guajava L. gene expression and prospection of functional microsatellite markers

Most research on trinucleotide repeats (TRs) focuses on human diseases, with few on the impact of TR expansions on plant gene expression. This work investigates TRs' effect on global gene expression in Psidium guajava L., a plant species with widespread distribution and significant relevance in the food, pharmacology, and economics sectors. We analyzed TR-containing coding sequences in 1,107 transcripts from 2,256 genes across root, shoot, young leaf, old leaf, and flower bud tissues of the Brazilian guava cultivars Cortibel RM and Paluma. Structural analysis revealed TR sequences with small repeat numbers (5-9) starting with cytosine or guanine or containing these bases. Functional annotation indicated TR-containing genes' involvement in cellular structures and processes (especially cell membranes and signal recognition), stress response, and resistance. Gene expression analysis showed significant variation, with a subset of highly expressed genes in both cultivars. Differential expression highlighted numerous down-regulated genes in Cortibel RM tissues, but not in Paluma, suggesting interplay between tissues and cultivars. Among 72 differentially expressed genes with TRs, 24 form miRNAs, 13 encode transcription factors, and 11 are associated with transposable elements. In addition, a set of 20 SSR-annotated, transcribed, and differentially expressed genes with TRs was selected as phenotypic markers for Psidium guajava and, potentially for closely related species as well.

Importance of CD8 Tex cell-associated gene signatures in the prognosis and immunology of osteosarcoma

As a highly aggressive bone malignancy, osteosarcoma poses a significant therapeutic challenge, especially in the setting of metastasis or recurrence. This study aimed to investigate the potential of CD8-Tex cell-associated genes as prognostic biomarkers to reveal the immunogenomic profile of osteosarcoma and guide therapeutic decisions. mRNA expression data and clinical details of osteosarcoma patients were obtained from the TCGA database (TARGET-OS dataset). The GSE21257 dataset (from the GEO database) was used as an external validation set to provide additional information on osteosarcoma specimens. 84 samples from the TARGET-OS dataset were used as the training set, and 53 samples from the GSE21257 dataset served as the external validation cohort. Univariate Cox regression analysis was utilized to identify CD8 Tex cell genes associated with prognosis. The LASSO algorithm was performed for 1000 iterations to select the best subset to form the CD8 Tex cell gene signature (TRS). Final genes were identified using the multivariate Cox regression model of the LASSO algorithm. Risk scores were calculated to categorize patients into high- and low-risk groups, and clinical differences were explored by Kaplan-Meier survival analysis to assess model performance. Prediction maps were constructed to estimate 1-, 3-, and 5 year survival rates for osteosarcoma patients, including risk scores for CD8 Texcell gene markers and clinicopathologic factors. The ssGSEA algorithm was used to assess the differences in immune function between TRS-defined high- and low-risk groups. TME and immune cell infiltration were further assessed using the ESTIMATE and CIBERSORT algorithms. To explore the relationship between immune checkpoint gene expression levels and the two risk-defined groups. A CD8 Tex cell-associated gene signature was extracted from the TISCH database and prognostic markers including two genes were developed. The high-risk group showed lower survival, and model performance was validated by ROC curves and C-index. Predictive plots were constructed to demonstrate survival estimates, combining CD8 Tex cell gene markers and clinical factors. This study provides valuable insights into the molecular and immune characteristics of osteosarcoma and offers potential avenues for advances in therapeutic approaches.

A novel defined programmed cell death related gene signature for predicting the prognosis of serous ovarian cancer

PURPOSE

This study aims to explore the contribution of differentially expressed programmed cell death genes (DEPCDGs) to the heterogeneity of serous ovarian cancer (SOC) through single-cell RNA sequencing (scRNA-seq) and assess their potential as predictors for clinical prognosis.

METHODS

SOC scRNA-seq data were extracted from the Gene Expression Omnibus database, and the principal component analysis was used for cell clustering. Bulk RNA-seq data were employed to analyze SOC-associated immune cell subsets key genes. CIBERSORT and single-sample gene set enrichment analysis (ssGSEA) were utilized to calculate immune cell scores. Prognostic models and nomograms were developed through univariate and multivariate Cox analyses.

RESULTS

Our analysis revealed that 48 DEPCDGs are significantly correlated with apoptotic signaling and oxidative stress pathways and identified seven key DEPCDGs (CASP3, GADD45B, GNA15, GZMB, IL1B, ISG20, and RHOB) through survival analysis. Furthermore, eight distinct cell subtypes were characterized using scRNA-seq. It was found that G protein subunit alpha 15 (GNA15) exhibited low expression across these subtypes and a strong association with immune cells. Based on the DEGs identified by the GNA15 high- and low-expression groups, a prognostic model comprising eight genes with significant prognostic value was constructed, effectively predicting patient overall survival. Additionally, a nomogram incorporating the RS signature, age, grade, and stage was developed and validated using two large SOC datasets.

CONCLUSION

GNA15 emerged as an independent and excellent prognostic marker for SOC patients. This study provides valuable insights into the prognostic potential of DEPCDGs in SOC, presenting new avenues for personalized treatment strategies.