Make Interactive Complex Heatmaps in R

Exploring macrophage heterogeneity in IgA nephropathy: Mechanisms of renal impairment and current therapeutic targets

The lack of understanding of the mechanism of renal injury in IgA nephropathy (IgAN) hinders the development of personalized treatment plans and targeted therapies. Improved insight into the cause of renal dysfunction in IgAN is necessary to enhance the effectiveness of strategies for slowing the progression of the disease. This study examined single cell RNA sequencing (scRNA seq) and bulk-RNA seq data and found that the gene expression of renal intrinsic cells (RIC) was significantly changed in patients with renal impairment, with a primary focus on energy metabolism. We discovered a clear metabolic reprogramming of RIC during renal function impairment (RF) using the 'scMetabolism' package, which manifested as a weakening of oxidative phosphorylation, alterations in fatty acid metabolism, and changes in glycolysis. Cellular communication analysis revealed that communication between macrophages (Ma) and RIC became more active and impacted cell function through the ligand-receptor-transcription factor (L-R-TF) axis in patients with RF. Our studies showed a notable upsurge in the expression of gene CLU and the infiltration of CLU+ Ma in patients with RF. CLU is a multifunctional protein, extensively involved in processes such as cell apoptosis and immune responses. Data obtained from the Nephroseq V5 database and multiplex immunohistochemistry (mIHC) were used to validate the findings, which were found to be robustly correlated with estimated glomerular filtration rate (eGFR) of the IgAN patients, as demonstrated by linear regression (LR). This study provides new insights into the cellular and molecular changes that occur in IgAN during renal impairment, revealing that elevated expression of CLU and CLU+ Ma percolation are common features in patients with RF. These findings offer potential targets and strategies for personalized management and targeted therapy of IgAN.

Metabolite profiling and transcriptome analyses reveal defense regulatory network against pink tea mite invasion in tea plant

BACKGROUND

The tea plant Camellia sinensis (L.) O. Kuntze is a perennial crop, invaded by diversity of insect pest species, and pink tea mite is one of the most devastating pests for sustainable tea production. However, molecular mechanism of defense responses against pink tea mites in tea is still unknown. In this study, metabolomics and transcriptome profiles of susceptible and resistant tea varieties were compared before and after pink tea mite infestation.

RESULTS

Metabolomics analysis revealed that abundance levels of polyphenol-related compounds changed significantly before and after infestation. At the transcript level, nearly 8 GB of clean reads were obtained from each sequenced library, and a comparison of infested plants of resistant and susceptible tea varieties revealed 9402 genes with significant differential expression. An array of genes enriched in plant pathogen interaction and biosynthetic pathways of phenylpropanoids showed significant differential regulation in response to pink tea mite invasion. In particular, the functional network linkage of disease resistant proteins, phenylalanine ammonia lyase, flavanone -3-hydroxylase, hydroxycinnamoyl-CoA shikimate transferase, brassinosteroid-6-oxidase 1, and gibberellin 2 beta-dioxygenase induced dynamic defense signals to suppress prolonged pink tea mite attacks. Further integrated analyses identified a complex network of transcripts and metabolites interlinked with precursors of various flavonoids that are likely modulate resistance against to pink tea mite.

CONCLUSIONS

Our results characterized the profiles of insect induced metabolic and transcript reprogramming and identified a defense regulatory network that can potentially be used to fend off pink tea mites damage.

Distinct negative-sense RNA viruses induce a common set of transcripts encoding proteins forming an extensive network

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

Key genes and immune pathways in T-cell mediated rejection post-liver transplantation identified via integrated RNA-seq and machine learning

Liver transplantation is the definitive treatment for end-stage liver disease, yet T-cell mediated rejection (TCMR) remains a major challenge. This study aims to identify key genes associated with TCMR and their potential biological processes and mechanisms. The GSE145780 dataset was subjected to differential expression analysis, weighted gene co-expression network analysis (WGCNA), and machine learning algorithms to pinpoint key genes associated with TCMR. Gene Set Enrichment Analysis (GSEA), immune infiltration analysis, and regulatory networks were constructed to ascertain the biological relevance of these genes. Expression validation was performed using single-cell RNA-seq (scRNA-seq) data and liver biopsy tissues from patients. We identified 5 key genes (ITGB2, FCER1G, IL-18, GBP1, and CD53) that are associated with immunological functions, such as chemotactic activity, antigen processing, and T cell differentiation. GSEA highlighted enrichment in chemokine signaling and antigen presentation pathways. A lncRNA-miRNA-mRNA network was delineated, and drug target prediction yielded 26 potential drugs. Evaluation of expression levels in non-rejection (NR) and TCMR groups exhibited significant disparities in T cells and myeloid cells. Tissue analyses from patients corroborated the upregulation of GBP1, IL-18, CD53, and FCER1G in TCMR cases. Through comprehensive analysis, this research has identified 4 genes intimately connected with TCMR following liver transplantation, shedding light on the underlying immune activation pathways and suggesting putative targets for therapeutic intervention.

Molecular mechanisms and diagnostic model of glioma-related epilepsy

Epilepsy is one of the most common symptoms in patients with gliomas; however, the mechanisms underlying its interaction are not yet clear. Moreover, epidemiological studies have not accurately identified patients with glioma-related epilepsy (GRE), and there is an urgent need to identify the molecular mechanisms and markers of its occurrence. We analyzed the demographics, transcriptome, whole-genome, and methylation sequences of 997 patients with glioma, to determine the genetic differences between glioma and GRE patients and to determine the upregulated molecular function, cellular composition, biological processes involved, signaling pathways, and immune cell infiltration. Twelve machine learning algorithms were refined into 113 combinatorial algorithms for building diagnostic recognition models. A total of 342 patients with GRE were identified with WHO grade 2 (174), grade 3 (107), and grade 4 (61). The mean age of the patients with GREs, with IDH mutations (n = 217 [63%]) and 1p19q non-codeletion (n = 169 [49%]), was 38 years old. GRE molecular functions were mainly passive transmembrane transporter protein activity, ion channel activity, and gated channel activity. Cellular components were enriched in the cation-channel and transmembrane transporter complexes. Cerebral cortical development regulates the membrane potential and synaptic organization as major biological processes. The signaling pathways mainly focused on cholinergic, GABAergic, and glutamatergic synapses. LASSO, combined with Random Forest, was the best diagnostic model and identified nine diagnostic genes. This study provides new insights and future perspectives for resolving the molecular mechanisms of GRE.

Children with eosinophilic esophagitis non-responsive to combination therapy have distinct esophageal transcriptomic and microbiome profile

BACKGROUND

A combination of proton-pump inhibitors (PPI) and topical steroids (TS) is used to treat children with eosinophilic esophagitis (EoE). However, a subset of children do not respond to this combination therapy. We aimed to identify the esophageal transcriptional, cell composition, and microbial differences between the non-responders (EoE-PPI-TSnr; n = 7) and responders (EoE-PPI-TSr; n = 7) to the combination therapy for EoE and controls (n = 9) using metatranscriptomics.

METHODS

Differential gene expression analysis was used to identify transcriptional differences, validated using the EoE diagnostic panel (EDP). Deconvolution analysis was performed to identify differences in their cell type composition. Microbiome analysis was conducted from esophageal biopsies RNAseq data, and microbial abundance was correlated with esophageal gene expression.

RESULTS

In all, 3164 upregulated and 3154 downregulated genes distinguished EoE-PPI-TSnr from EoE-PPI-TSr. Eosinophilic inflammatory response, cytokine signaling, and collagen formation pathways were significantly upregulated in EoE-PPI-TSnr. There was a 56% overlap in dysregulated genes between EoE-PPI-TSnr and EDP, with a perfect agreement in the directionality of modulation. Eosinophils, dendritic cells (DCs), immature DCs, megakaryocytic-erythroid progenitors, and T helper type 1 cells were significantly higher in EoE-PPI-TSnr. There was no significant difference in microbiome diversity. The relative abundance of Fusobacterium sp. and Acinetobacter sp. notably differed in EoE-PPI-TSnr and correlated with the key pathways.

CONCLUSION

Our results provide critical insights into the molecular, cellular, and microbial factors associated with the lack of response to PPI and TS combination therapy in children with EoE. This study advances our understanding of the pathobiology of EoE while guiding personalized treatment strategies.

Time series analysis combined with transcriptome sequencing to explore characteristic genes and potential molecular mechanisms associated with ultrasound-guided microwave ablation of glioma

OBJECTIVE

This study aimed to explore marker genes and their potential molecular mechanisms involved in US-guided MWA for glioma in mice.

METHOD

The differentially expressed genes (DEGs1 and DEGs2) and lncRNAs (DELs1 and DELs2) were obtained between Non (glioma tissues without MWA) and T0 groups (0h after MWA), as well as between Non and T24 groups (24h after MWA). The down-regulation cluster genes (CONDOWNDEGs) and upregulation cluster genes (CONUPDEGs) were identified by time series analysis. Candidate genes were obtained by overlapping CONDOWNDEGs with downregulation DEGs (DOWNDEGs)1 and DOWNDEGs2, as well as CONUPDEGs with up-regulation DEGs (UPDEGs)1 and UPDEGs2. The expressions of immune checkpoints and inflammatory factors, gene set enrichment analysis (GSEA), and protein subcellular localization were performed. The eXpression2Kinases (X2K), GeneMANIA, transcription factor (TF), and competing endogenous (ce) RNA regulatory networks were conducted. The expression of marker genes was validated by quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

Five marker genes (IL32, VCAM1, IL34, NFKB1 and CXCL13) were identified, which were connected with immune-related functions. Two immune checkpoints (CD96 and TIGIT) and six inflammatory factors played key roles in US-guided MWA for glioma. ceRNA regulatory networks revealed that miR-625-5p, miR-625-3p, miR-31-5p and miR-671-5p were associated with target genes. qRT-PCR indicated both IL32, VCAM1, and NFKB1 were potential markers under US-guided MWA-related time series analysis.

CONCLUSION

The use of US-guided MWA might be a practical method for influencing the function of target genes, regulating time frames to decrease inflammation, and stimulating immune responses in glioma therapy.

Precision Oncology and Systemic Targeted Therapy in Pseudomyxoma Peritonei

PURPOSE

Pseudomyxoma peritonei (PMP) is a rare and poorly understood malignant condition characterized by the accumulation of intra-abdominal mucin produced from peritoneal metastases. Currently, cytoreductive surgery remains the mainstay of treatment but disease recurrence and death after relapse frequently occur in patients with PMP. New therapeutic strategies are therefore urgently needed for these patients.

EXPERIMENTAL DESIGN

A total of 120 PMP samples from 50 patients were processed to generate a collection of 50 patient-derived organoid (PDO) and xenograft (PDX) models. Whole exome sequencing, immunohistochemistry analyses, and in vitro and in vivo drug efficacy studies were performed.

RESULTS

In this study, we have generated a collection of PMP preclinical models and identified druggable targets, including BRAFV600E, KRASG12C, and KRASG12D, that could also be detected in intra-abdominal mucin biopsies of patients with PMP using droplet digital PCR. Preclinical models preserved the histopathological markers from the original patient sample. The BRAFV600E inhibitor encorafenib reduced cell viability of BRAFV600E PMP-PDO models. Proof-of-concept in vivo experiments showed that a systemic treatment with encorafenib significantly reduced tumor growth and prolonged survival in subcutaneous and orthotopic BRAFV600E-PMP-PDX mouse models.

CONCLUSIONS

Our study demonstrates for the first time that systemic targeted therapies can effectively control PMP tumors. BRAF signaling pathway inhibition represents a new therapeutic opportunity for patients with BRAFV600E PMP who have a poor prognosis. Importantly, our present data and collection of preclinical models pave the way for evaluating the efficacy of other systemic targeted therapies toward extending the promise of precision oncology to patients with PMP.

Seasonal Dynamics of Eukaryotic Microbial Communities in the Water-Receiving Reservoir of the Long-Distance Water Diversion Project, China

Inter-basin water transfer projects, such as the Yellow River to Qingdao Water Diversion Project (YQWD), are essential for addressing water scarcity, but impact local aquatic ecosystems. This study investigates the seasonal characteristics of eukaryotic microbial communities in the Jihongtan Reservoir, the main water-receiving body of YQWD, over a one-year period using 18S rDNA amplicon sequencing. The results showed that the eukaryotic microbial diversity did not exhibit significant seasonal variation (p > 0.05), but there was a notable variance in the community structure (p < 0.05). Arthropoda and Paracyclopina, representing the most dominant phylum and the most dominant genus, respectively, both exhibited the lowest abundance during the winter. The Chlorophyta, as the second-dominant phylum, demonstrates its higher abundance in the spring and winter. The Mantel test and PLS-PM (Partial Least Squares Path Modeling) revealed that water temperature (WT), dissolved oxygen (DO), and pH influenced the seasonal dynamic of eukaryotic microbial communities significantly, of which WT was the primary driving factor. In addition to environmental factors, water diversion is likely to be an important influencing factor. The results of the co-occurrence network and robustness suggested that the spring network is the most complex and exhibits the highest stability. Moreover, keystone taxa within networks have been identified, revealing that these key groups encompass both abundant and rare species, with specificity to different seasons. These insights are vital for understanding the seasonal variation of microbial communities in the Jihongtan Reservoir during ongoing water diversions.

System analysis based on T-cell exhaustion-related genes identifies PTPRT as a promising diagnostic and prognostic biomarker for gastric cancer

Multiple investigations have demonstrated the crucial involvement of T-cell exhaustion (TEX) in anti-tumor immune response and their strong correlation with prognosis. This study aimed at creating a strong signature using TEX for gastric cancer through bioinformatics analysis and experimental validation. We utilized data from The Cancer Genome Atlas (TCGA) databases to retrieve RNA-seq data from patients with stomach adenocarcinoma (STAD). Genes related to TEX were discovered using gene set variance analysis (GSVA) and weighted gene correlation network analysis (WGCNA). Subsequently, prognostic signature based on TEX was developed using LASSO-Cox analysis. Relationship between key genes and immune cells were examined. Finally, biological function of a key TEX-related gene PTPRT in gastric cancer was verified by in vivo experiment. A total of 29 TEX-related biomarkers were screened by WGCNA and random forest. Among them, five core signatures (PTPRT, CAV2, PPIH, PRDM2, and FGF1), further identified by LASSO-Cox, were considered as strong predictors of prognosis for gastric cancer and associated with immune infiltration. PTPRT gene had the largest number of SNPs, with the most mutation types. In vivo experiments revealed that PTPRT overexpression significantly inhibited tumor malignant progression and accelerated apoptosis through stimulating the secretion of killer cytokines such as TNF-α and IFN-γ. In addition, flow cytometry revealed that PTPRT overexpression alleviated TEX by increasing the abundance of CD8+ T cells, with inhibition of cell surface PD-1 and Tim-3. The predictive prognostic value of TEX gene expression levels was evaluated in patients with gastric cancer, providing a new perspective for precision immuno-oncology studies.

Identifying key inflammatory genes in psoriasis via weighted gene co-expression network analysis: Potential targets for therapy

Psoriasis is a globally prevalent chronic inflammatory skin disease. This study aimed to scrutinize the hub genes related to inflammation and potential molecular mechanisms in psoriasis. Utilizing mRNA expression profiles from public datasets GSE13355, GSE78097, and GSE14905, we set up a comprehensive analysis. Initially, we selected differentially expressed genes (DEGs) from psoriasis and control samples in GSE13355, followed by calculating inflammatory indices using genomic set variation analysis (GSVA). Weighted gene co-expression network analysis (WGCNA) was then applied to link significant modules with the inflammatory index. This process helped us identify differentially expressed inflammation-related genes (DE-IRGs). A protein-protein interaction (PPI) network was established, with the molecular complex detection (MCODE) plug-in pinpointing six chemokine genes (CCR7, CCL2, CCL19, CXCL8, CXCL1, and CXCL2) as central hub genes. These genes demonstrated pronounced immunohistochemical staining in psoriatic tissues compared to normal skin. Notably, the CCR7 gene exhibited the highest potential for m6A modification sites. Furthermore, we constructed transcription factor-microRNA-mRNA networks, identifying 139 microRNAs and 52 transcription factors associated with the hub genes. For the LASSO logistic regression model, the area under the curve (AUC) in the training set was 1, and in the two validation cohorts GSE78097 and GSE14905 were 1 and 0.872, respectively. In conclusion, our study highlights six chemokine genes (CCR7, CCL2, CCL19, CXCL8, CXCL1, and CXCL2) as potential biomarkers in psoriasis, providing insights into the immune and inflammatory responses as pivotal instances in disease pathogenesis. These findings pave the way for exploring new therapeutic targets, particularly focusing on chemokine-associated pathways in psoriasis treatment.

Patient-derived Organoid Pharmacotyping As A Predictive Tool for Therapeutic Selection in Pancreatic Ductal Adenocarcinoma

OBJECTIVE

We integrate a new approach to chemosensitivity data for clinically-relevant regimen matching, and demonstrate the relationship with clinical outcomes in a large PDO biobank.

SUMMARY BACKGROUND DATA

Pancreatic ductal adenocarcinoma (PDAC) usually recurs following potentially curative resection. Prior studies related patient-derived organoid (PDO) chemosensitivity with clinical responses.

METHODS

PDOs were established from pre-treatment biopsies in a multi-institution clinical trial (n=21) and clinical specimens at a high-volume pancreatectomy center (n=74, of which 48 were pre-treated). PDO in vitro chemosensitivities to standard-of-care chemotherapeutics (pharmacotypes) were matched to potential clinically-relevant regimens by a weighted nearest-neighbors analysis. Clinical outcomes were then compared for patients who had well-matched versus poorly-matched treatment according to this metric.

RESULTS

Our function matched 91% of PDOs to a standard-of-care regimen (9% pan-resistant). PDOs poorly-matched to the neoadjuvant regimen received would have matched to an alternative in 34% of cases. Patients receiving neoadjuvant chemotherapy well-matched to their pharmacotype experienced improved CA 19-9 response (60% decreased to normal when well-matched, 29% when poorly-matched, P<0.05) and lymph node down-staging (33% N0 after poorly-matched, 69% after well-matched, P<0.05). Patients receiving both well-matched neoadjuvant and adjuvant chemotherapy experienced improved recurrence-free- and overall survival (median RFS 8.5 mo poorly-matched, 15.9 mo well-matched, P<0.05; median OS 19.5 vs. 30.3 mo, P<0.05).

CONCLUSION

In vitro PDO pharmacotyping can inform PDAC therapy selection. We demonstrate improved outcomes including survival for patients treated with regimens well-matched to their PDO chemosensitivities. A subsequent prospective study using PDO pharmacotype matching could improve oncologic outcomes and improve quality of life by avoiding therapies not expected to be effective.

Bioinformatics Identification and Experimental Verification of Disulfidptosis-Related Genes in the Progression of Osteoarthritis

BACKGROUND

Osteoarthritis (OA) is a disabling and highly prevalent condition affecting millions worldwide. Recently discovered, disulfidptosis represents a novel form of cell death induced by the excessive accumulation of cystine. Despite its significance, a systematic exploration of disulfidptosis-related genes (DRGs) in OA is lacking.

METHODS

This study utilized three OA-related datasets and DRGs. Differentially expressed (DE)-DRGs were derived by intersecting the differentially expressed genes (DEGs) from GSE114007 with DRGs. Feature genes underwent screening through three machine learning algorithms. High diagnostic value genes were identified using the receiver operating characteristic curve. Hub genes were confirmed through expression validation. These hub genes were then employed to construct a nomogram and conduct enrichment, immune, and correlation analyses. An additional validation of hub genes was performed through in vitro cell experiments.

RESULTS

SLC3A2 and PDLIM1 were designated as hub genes, displaying excellent diagnostic performance. PDLIM1 exhibited low expression in early chondrocyte differentiation, rising significantly in the late stage, while SLC3A2 showed high overall expression, declining in the late differentiation stage. Cellular experiments corroborated the correlation of SLC3A2 and PDLIM1 with chondrocyte inflammation.

CONCLUSIONS

Two hub genes, SLC3A2 and PDLIM1, were identified in relation to disulfidptosis, providing potential directions for diagnosing and treating OA.

Knockdown of SDC-1 Gene Alleviates the Metabolic Pathway for the Development of MODS

This study aims to reveal the metabolic differences between SDC-1 knockout mice and wild-type mice and the metabolic differences caused by shock in SDC-1 knockout mice by integrating transcriptomics and metabolomics. A total of 1009 differential metabolites were differentially expressed based on untargeted metabolomics and high-resolution mass spectrometry detection techniques. According to Kyoto Encyclopedia of Genes and Genomes enrichment, SDC-1 knockout significantly altered fat digestion and absorption, GnRH signaling pathway, fructose and mannose metabolism, and some other amino-related metabolic pathways and significantly modulated positively regulated longevity regulatory pathways, longevity regulatory pathways-worm, nicotinamide and niacinamide metabolism, and vitamin digestion and absorption pathways after its shock. Our findings indicate that SDC-1 knockout may have potential therapeutic effects in hemorrhagic shock by increasing nicotinamide metabolism.

The IRE1α/XBP1 signaling axis drives myoblast fusion in adult skeletal muscle

Skeletal muscle regeneration involves a signaling network that regulates the proliferation, differentiation, and fusion of muscle precursor cells to injured myofibers. IRE1α, one of the arms of the unfolded protein response, regulates cellular proteostasis in response to ER stress. Here, we demonstrate that inducible deletion of IRE1α in satellite cells of mice impairs skeletal muscle regeneration through inhibiting myoblast fusion. Knockdown of IRE1α or its downstream target, X-box protein 1 (XBP1), also inhibits myoblast fusion during myogenesis. Transcriptome analysis revealed that knockdown of IRE1α or XBP1 dysregulates the gene expression of molecules involved in myoblast fusion. The IRE1α-XBP1 axis mediates the gene expression of multiple profusion molecules, including myomaker (Mymk). Spliced XBP1 (sXBP1) transcription factor binds to the promoter of Mymk gene during myogenesis. Overexpression of myomaker in IRE1α-knockdown cultures rescues fusion defects. Inducible deletion of IRE1α in satellite cells also inhibits myoblast fusion and myofiber hypertrophy in response to functional overload. Collectively, our study demonstrates that IRE1α promotes myoblast fusion through sXBP1-mediated up-regulation of the gene expression of multiple profusion molecules, including myomaker.

Tissue gene expression profiles and communication networks inform candidate blood biomarker identification in psoriasis and atopic dermatitis

Overlapping clinical and pathomechanistic features can complicate the diagnosis and treatment of inflammatory skin diseases, including psoriasis and atopic dermatitis (AD). Spatial transcriptomics allows the identification of disease- and cell-specific molecular signatures that may advance biomarker development and future treatments. This study identified transcriptional signatures in keratinocytes and sub-basal CD4+ and CD8+ T lymphocytes from patients with psoriasis and AD. In silico prediction of ligand:receptor interactions delivered key signalling pathways (interferon, effector T cells, stroma cell and matrix biology, neuronal development, etc.). Targeted validation of selected transcripts, including CCL22, RELB, and JUND, in peripheral blood T cells suggests the chosen approach as a promising tool also in other inflammatory diseases. Psoriasis and AD are characterized by transcriptional dysregulation in T cells and keratinocytes that may be targeted therapeutically. Spatial transcriptomics is a valuable tool in the search for molecular signatures that can be used as biomarkers and/or therapeutic targets.

Integrative bioinformatics analysis to identify ferroptosis-related genes in non-obstructive azoospermia

PURPOSE

The objective of this study was to discern ferroptosis-related genes (FRGs) linked to non-obstructive azoospermia and investigate the associated molecular mechanisms.

METHOD

A dataset related to azoospermia was retrieved from the Gene Expression Omnibus database, and FRGs were sourced from GeneCards. Ferroptosis-related differentially expressed genes (FRDEGs) were discerned. Subsequently, these genes underwent analyses encompassing Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, as well as protein-protein interaction (PPI) networks and assessments of functional similarity. Following the identification of hub genes, an exploration of immune infiltration, single-cell expression, diagnostic utility, and interactions involving hub genes, RNA-binding proteins (RBPs), transcription factors (TFs), microRNAs (miRNAs), and drugs was conducted.

RESULTS

A total of 35 differentially expressed FRGs were discerned. These genes demonstrated enrichment in functions and pathways associated with ferroptosis. From the PPI network, eight hub genes were selected. Functional similarity analysis highlighted the potential pivotal roles of HMOX1 and GPX4 in azoospermia. Analysis of immune cell infiltration indicated a significant decrease in activated dendritic cells in the azoospermia group, with notable correlations between hub genes, particularly SAT1 and HMGCR, and immune cell infiltration. Unique expression patterns of hub genes across various cell types in the human testis were observed, with GPX4 prominently enriched in spermatid/sperm. Eight hub genes exhibited robust diagnostic value (AUC > 0.75). Lastly, a comprehensive hub gene-miRNA-TF-RBP-drug network was constructed.

CONCLUSION

In summary, our investigation unveiled eight FRDEGs associated with azoospermia, which hold potential as biomarkers for the diagnosis and treatment of azoospermia.

Detection of neoplastic-immune hybrid cells with metastatic properties in uveal melanoma

BACKGROUND

Uveal melanoma is the most common non-cutaneous melanoma and is an intraocular malignancy affecting nearly 7,000 individuals per year worldwide. Of these, approximately 50% will progress to metastatic disease for which there are currently no effective curative therapies. Despite advances in molecular profiling and metastatic stratification of uveal melanoma tumors, little is known regarding their underlying biology of metastasis. Our group has identified a disseminated neoplastic cell population characterized by co-expression of immune and melanoma proteins, circulating hybrid cells (hybrids), in patients with uveal melanoma. Compared to circulating tumor cells, which lack expression of immune proteins, hybrids are detected at an increased prevalence in peripheral blood and can be used as a non-invasive biomarker to predict metastatic progression.

METHODS

To ascertain mechanisms underlying enhanced hybrid cell dissemination we identified hybrid cells within primary uveal melanoma tumors using single cell RNA sequencing (n = 8) and evaluated their gene expression and predicted ligand-receptor interactions in relation to other melanoma and immune cells within the primary tumor. We then verified expression of upregulated hybrid pathways within patient-matched tumor and peripheral blood hybrids (n = 4) using cyclic immunofluorescence and quantified their protein expression relative to other non-hybrid tumor and disseminated tumor cells.

RESULTS

Among the top upregulated genes and pathways in hybrid cells were those involved in enhanced cell motility and cytoskeletal rearrangement, immune evasion, and altered cellular metabolism. In patient-matched tumor and peripheral blood, we verified gene expression by examining concordant protein expression for each pathway category: TMSB10 (cell motility), CD74 (immune evasion) and GPX1 (metabolism). Both TMSB10 and GPX1 were expressed on significantly higher numbers of disseminated hybrid cells compared to circulating tumor cells, and CD74 and GPX1 were expressed on more disseminated hybrids than tumor-resident hybrids. Lastly, we identified that hybrid cells express ligand-receptor signaling pathways implicated in promoting metastasis including GAS6-AXL, CXCL12-CXCR4, LGALS9-P4HB and IGF1-IGFR1.

CONCLUSION

These findings highlight the importance of TMSB10, GPX1 and CD74 for successful hybrid cell dissemination and survival in circulation. Our results contribute to the understanding of uveal melanoma tumor progression and interactions between tumor cells and immune cells in the tumor microenvironment that may promote metastasis.

Impact of Missense Mutations on Spike Protein Stability and Binding Affinity in the Omicron Variant

The global effort to combat the COVID-19 pandemic faces ongoing uncertainty with the emergence of Variants of Concern featuring numerous mutations on the Spike (S) protein. In particular, the Omicron Variant is distinguished by 32 mutations, including 10 within its receptor-binding domain (RBD). These mutations significantly impact viral infectivity and the efficacy of vaccines and antibodies currently in use for therapeutic purposes. In our study, we employed structure-based computational saturation mutagenesis approaches to predict the effects of Omicron missense mutations on RBD stability and binding affinity, comparing them to the original Wuhan-Hu-1 strain. Our results predict that mutations such as G431W and P507W induce the most substantial destabilizations in the Wuhan-Hu-1-S/Omicron-S RBD. Notably, we postulate that mutations in the Omicron-S exhibit a higher percentage of enhancing binding affinity compared to Wuhan-S. We found that the mutations at residue positions G447, Y449, F456, F486, and S496 led to significant changes in binding affinity. In summary, our findings may shed light on the widespread prevalence of Omicron mutations in human populations. The Omicron mutations that potentially enhance their affinity for human receptors may facilitate increased viral binding and internalization in infected cells, thereby enhancing infectivity. This informs the development of new neutralizing antibodies capable of targeting Omicron's immune-evading mutations, potentially aiding in the ongoing battle against the COVID-19 pandemic.

CDK4/6 inhibition enhances T-cell immunotherapy on hepatocellular carcinoma cells by rejuvenating immunogenicity

Hepatocellular carcinoma (HCC) poses a significant clinical challenge, necessitating the integration of immunotherapeutic approaches. Palbociclib, a selective CDK4/6 inhibitor, has demonstrated promising efficacy in preclinical HCC models and is being evaluated as a novel therapeutic option in clinical trials. Additionally, CDK4/6 inhibition induces cellular senescence, potentially influencing the tumor microenvironment and immunogenicity of cancer cells. In this study, we conducted comprehensive bioinformatic analyses using diverse HCC transcriptome datasets, including bulk and single-cell RNA-sequencing data from public databases. We also utilized human and mouse HCC cells to investigate functional aspects. Primary T cells isolated from mouse blood were employed to assess T cell immunity against HCC cells. Results revealed that CD8+ T-cell infiltration correlates with improved outcomes in HCC patients with suppressed CDK4/6 expression. Moreover, CDK4/6 expression was associated with alterations in the immune landscape and immune checkpoint expression within the liver tumor microenvironment. Furthermore, we found that treatment with Palbociclib and Doxorubicin induces cellular senescence and a senescence-associated secretory phenotype in HCC cells. Notably, pretreatment with Palbociclib augmented T cell-mediated cytotoxicity against HCC cells, despite upregulation of PD-L1, surpassing the effects of Doxorubicin pretreatment. In conclusion, our study elucidates a novel mechanism by which CDK4/6 inhibition enhances T-cell-associated cancer elimination and proposes a potential therapeutic strategy to enhance T-cell immunotherapy on HCC.

ARID1A suppresses R-loop-mediated STING-type I interferon pathway activation of anti-tumor immunity

Clinical trials have identified ARID1A mutations as enriched among patients who respond favorably to immune checkpoint blockade (ICB) in several solid tumor types independent of microsatellite instability. We show that ARID1A loss in murine models is sufficient to induce anti-tumor immune phenotypes observed in ARID1A mutant human cancers, including increased CD8+ T cell infiltration and cytolytic activity. ARID1A-deficient cancers upregulated an interferon (IFN) gene expression signature, the ARID1A-IFN signature, associated with increased R-loops and cytosolic single-stranded DNA (ssDNA). Overexpression of the R-loop resolving enzyme, RNASEH2B, or cytosolic DNase, TREX1, in ARID1A-deficient cells prevented cytosolic ssDNA accumulation and ARID1A-IFN gene upregulation. Further, the ARID1A-IFN signature and anti-tumor immunity were driven by STING-dependent type I IFN signaling, which was required for improved responsiveness of ARID1A mutant tumors to ICB treatment. These findings define a molecular mechanism underlying anti-tumor immunity in ARID1A mutant cancers.

LRRK2 mediates haloperidol-induced changes in indirect pathway striatal projection neurons

Haloperidol is used to manage psychotic symptoms in several neurological disorders through mechanisms that involve antagonism of dopamine D2 receptors that are highly expressed in the striatum. Significant side effects of haloperidol, known as extrapyramidal symptoms, lead to motor deficits similar to those seen in Parkinson's disease and present a major challenge in clinical settings. The underlying molecular mechanisms responsible for these side effects remain poorly understood. Parkinson's disease-associated LRRK2 kinase has an important role in striatal physiology and a known link to dopamine D2 receptor signaling. Here, we systematically explore convergent signaling of haloperidol and LRRK2 through pharmacological or genetic inhibition of LRRK2 kinase, as well as knock-in mouse models expressing pathogenic mutant LRRK2 with increased kinase activity. Behavioral assays show that LRRK2 kinase inhibition ameliorates haloperidol-induced motor changes in mice. A combination of electrophysiological and anatomical approaches reveals that LRRK2 kinase inhibition interferes with haloperidol-induced changes, specifically in striatal neurons of the indirect pathway. Proteomic studies and targeted intracellular pathway analyses demonstrate that haloperidol induces a similar pattern of intracellular signaling as increased LRRK2 kinase activity. Our study suggests that LRRK2 kinase plays a key role in striatal dopamine D2 receptor signaling underlying the undesirable motor side effects of haloperidol. This work opens up new therapeutic avenues for dopamine-related disorders, such as psychosis, also furthering our understanding of Parkinson's disease pathophysiology.

Lipid- and glucose-lowering effects of Rhamnan sulphate from Monostroma nitidum with altered gut microbiota in mice

Rhamnan sulphate (RS) is a sulphated polysaccharide found in green algae such as Monostroma nitidum that exhibits various biological functions, including anticoagulant, antitumour, antiviral, and anti-obesity properties. In our previous clinical trial, we demonstrated that RS intake improves constipation. However, no specific bacteria showed a significant (p < .05) change. Notably, these results were obtained after a short RS inoculation period of only 2 weeks. In the present study, to evaluate the long-term effects of RS on the gut microbiota, we orally administered RS to BALB/c mice for 11 weeks, analyzed their blood biochemical data, and performed 16s rRNA-sequencing. Oral administration of RS increased body weight with increased food intake, whereas plasma total cholesterol and fasting plasma glucose levels decreased. RS-fed mice showed lower fasting insulin levels (p < .1) and decreased homeostatic model assessment for insulin resistance (HOMA-IR, p < .0001), suggesting that RS improved insulin resistance. In the feces of mice, the amounts of acetic and propionic acids increased. In the gut microbiota, predictive metagenomic profiling using the phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt2) revealed functional alterations in Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways in RS-fed mice. Corresponding to the blood glucose-lowering effect, the glycolysis and tricarboxylic acid (TCA) cycle pathways were activated. In addition, the Firmicutes/Bacteroides (F/B) ratio, which may be associated with various health outcomes, was also reduced. These results suggest that the blood glucose-lowering effect, improvement in insulin resistance, and lipid-lowering effect of RS may be due to changes in the intestinal microbiota.

RNA expression profiling of peritoneal metastasis from pancreatic cancer treated with Pressurized Intraperitoneal Aerosol Chemotherapy (PIPAC)

Objectives

Pressurized Intraperitoneal Aerosol Chemotherapy (PIPAC) is an experimental treatment option in peritoneal metastasis from pancreatic cancer (PM-PC). Aims were to examine mRNA profile of fibrosis due to response after systemic chemotherapy and PIPAC (Regression) compared to treatment-naïve PM-PC and chronic cholecystitis-related peritoneal fibrosis (Controls).

Methods

Peritoneal biopsies (PBs) from PM-PC patients who had undergone systemic chemotherapy and PIPAC were evaluated with Peritoneal Regression Grading Score (PRGS). We extracted RNA from PBs with Regression (PRGS 1, n=11), treatment-naïve PM-PC (n=10), and Controls (n=10). Profiling of 800 mRNAs was performed (NanoString nCounter, PanCancer Immuno-Oncology 360 (IO-360) and 30 additional stroma-related mRNAs).

Results

Regression vs. PM-PC identified six up-regulated and 197 down-regulated mRNAs (FDR≤0.05), linked to TNF-α signaling via NF-kB, G2M checkpoint, epithelial-mesenchymal transition, estrogen response, and coagulation. Regression vs. Controls identified 43 significantly up-regulated mRNAs, linked to interferon-α response, and down-regulation of 99 mRNAs, linked to TNF-α signaling via NF-kB, inflammatory response, epithelial-mesenchymal transition, KRAS signaling, and hypoxia (FDR≤0.05).

Conclusions

In regressive fibrosis of PM-PC after systemic chemotherapy and PIPAC (Regression), downregulation of mRNAs related to key tumor biological pathways was identified. Regression also showed transcriptional differences from unspecific, benign fibrosis (Controls). Future studies should explore whether mRNA profiling of PBs with PM from PC or other primaries holds prognostic or predictive value.

Identification of key mitochondria-related genes and their relevance to the immune system linking Parkinson's disease and primary Sjögren's syndrome through integrated bioinformatics analyses

BACKGROUND

Mitochondria are the metabolic hubs of cells, regulating energy production and antigen presentation, which are essential for activation, proliferation, and function of immune cells. Recent evidence indicates that mitochondrial antigen presentation may have an impact on diseases such as Parkinson's disease (PD) and autoimmune diseases. However, there is limited knowledge about the mechanisms that regulate the presentation of mitochondrial antigens in these diseases.

METHODS

In the current study, RNA sequencing was performed on labial minor salivary gland (LSG) from 25 patients with primary Sjögren's syndrome (pSS) and 14 non-pSS aged controls. Additionally, we obtained gene expression omnibus datasets associated with PD patients from NCBI Gene Expression Omnibus (GEO) databases. Single-sample Gene Set Enrichment Analysis (ssGSEA), ESTIMATE and Spearman correlations were conducted to explore the association between mitochondrial related genes and the immune system. Furthermore, we applied weighted Gene Co-expression Network Analysis (WGCNA) to identify hub mitochondria-related genes and investigate the correlated networks in both diseases. Single cell transcriptome analysis, immunohistochemical (IHC) staining and quantitative real-time PCR (qRT-PCR) were used to verify the activation of the hub mitochondria-related pathway. Pearson correlations and the CIBERSORT algorithms were employed to further reveal the correlation between hub mitochondria-related pathways and immune infiltration.

RESULTS

The transcriptome analysis revealed the presence of overlapping mitochondria-related genes and mitochondrial DNA damage in patients with pSS and PD. Reactive oxygen species (ROS), the senescence marker p53, and the inflammatory markers CD45 and Bcl-2 were found to be regionally distributed in LSGs of pSS patients. WGCNA analysis identified the STING pathway as the central mitochondria-related pathway closely associated with the immune system. Single cell analysis, IHC staining, and qRT-PCR confirmed the activation of the STING pathway. Subsequent, bioinformatic analysis revealed the proportion of infiltrating immune cells in the STING-high and STING-low groups of pSS and PD. Furthermore, the study demonstrated the association of the STING pathway with innate and adaptive immune cells, as well as functional cells, in the immune microenvironment of PD and pSS.

CONCLUSION

Our study uncovered a central pathway that connects mitochondrial dysfunction and the immune microenvironment in PD and pSS, potentially offering valuable insights into therapeutic targets for these conditions.

Proteome-based molecular subtyping and therapeutic target prediction in gastric cancer

Different molecular classifications for gastric cancer (GC) have been proposed based on multi-omics platforms with the long-term goal of improved precision treatment. However, the GC (phospho)proteome remains incompletely characterized, particularly at the level of tyrosine phosphorylation. In addition, previous multiomics-based stratification of patient cohorts has lacked identification of corresponding cell line models and comprehensive validation of broad or subgroup-selective therapeutic targets. To address these knowledge gaps, we applied a reverse approach, undertaking the most comprehensive (phospho)proteomic analysis of GC cell lines to date and cross-validating this using publicly available data. Mass spectrometry (MS)-based (phospho)proteomic and tyrosine phosphorylation datasets were subjected to individual or integrated clustering to identify subgroups that were subsequently characterized in terms of enriched molecular processes and pathways. Significant congruence was detected between cell line proteomic and specific patient-derived transcriptomic subclassifications. Many protein kinases exhibiting 'outlier' expression or phosphorylation in the cell line dataset exhibited genomic aberrations in patient samples and association with poor prognosis, with casein kinase I isoform delta/epsilon (CSNK1D/E) being experimentally validated as potential therapeutic targets. Src family kinases were predicted to be commonly hyperactivated in GC cell lines, consistent with broad sensitivity to the next-generation Src inhibitor eCF506. In addition, phosphoproteomic and integrative clustering segregated the cell lines into two subtypes, with epithelial-mesenchyme transition (EMT) and proliferation-associated processes enriched in one, designated the EMT subtype, and metabolic pathways, cell-cell junctions, and the immune response dominating the features of the other, designated the metabolism subtype. Application of kinase activity prediction algorithms and interrogation of gene dependency and drug sensitivity databases predicted that the mechanistic target of rapamycin kinase (mTOR) and dual specificity mitogen-activated protein kinase kinase 2 (MAP2K2) represented potential therapeutic targets for the EMT and metabolism subtypes, respectively, and this was confirmed using selective inhibitors. Overall, our study provides novel, in-depth insights into GC proteomics, kinomics, and molecular taxonomy and reveals potential therapeutic targets that could provide the basis for precision treatments.

Risk assessment model based on nucleotide metabolism-related genes highlights SLC27A2 as a potential therapeutic target in breast cancer

PURPOSE

Breast cancer (BC) is the most prevalent malignant tumor worldwide among women, with the highest incidence rate. The mechanisms underlying nucleotide metabolism on biological functions in BC remain incompletely elucidated. MATERIALS AND METHODS: We harnessed differentially expressed nucleotide metabolism-related genes from The Cancer Genome Atlas-BRCA, constructing a prognostic risk model through univariate Cox regression and LASSO regression analyses. A validation set and the GSE7390 dataset were used to validate the risk model. Clinical relevance, survival and prognosis, immune infiltration, functional enrichment, and drug sensitivity analyses were conducted.

RESULTS

Our findings identified four signature genes (DCTPP1, IFNG, SLC27A2, and MYH3) as nucleotide metabolism-related prognostic genes. Subsequently, patients were stratified into high- and low-risk groups, revealing the risk model's independence as a prognostic factor. Nomogram calibration underscored superior prediction accuracy. Gene Set Variation Analysis (GSVA) uncovered activated pathways in low-risk cohorts and mobilized pathways in high-risk cohorts. Distinctions in immune cells were noted between risk cohorts. Subsequent experiments validated that reducing SLC27A2 expression in BC cell lines or using the SLC27A2 inhibitor, Lipofermata, effectively inhibited tumor growth.

CONCLUSIONS

We pinpointed four nucleotide metabolism-related prognostic genes, demonstrating promising accuracy as a risk prediction tool for patients with BC. SLC27A2 appears to be a potential therapeutic target for BC among these genes.

Fungal diversity in sediments of the eastern tropical Pacific oxygen minimum zone revealed by metabarcoding

Oxygen minimum zones (OMZ) represent ~8% of the ocean, with the Pacific as the largest and top expanding area. These regions influence marine ecosystems, promoting anaerobic microbial communities. Nevertheless, only a fraction of microbial diversity has been studied, with fungi being the less explored component. So, herein we analyzed fungal diversity patterns in surface and subsurface sediments along a bathymetric transect using metabarcoding of the ITS1 region in the OMZ of the Mexican Pacific off Mazatlán. We identified 353 amplicon sequence variants (ASV), within the Ascomycota, Basidiomycota, and Rozellomycota. Spatial patterns evidenced higher alpha diversity in nearshore and subsurface subsamples, probably due to temporal fluctuations in organic matter inputs. Small-scale heterogeneity characterized the community with the majority of ASV (269 ASV) occurring in a single subsample, hinting at the influence of local biogeochemical conditions. This baseline data evidenced a remarkable fungal diversity presenting high variation along a bathymetric and vertical transects.

Deciphering infected cell types, hub gene networks and cell-cell communication in infectious bronchitis virus via single-cell RNA sequencing

Infectious bronchitis virus (IBV) is a coronavirus that infects chickens, which exhibits a broad tropism for epithelial cells, infecting the tracheal mucosal epithelium, intestinal mucosal epithelium, and renal tubular epithelial cells. Utilizing single-cell RNA sequencing (scRNA-seq), we systematically examined cells in renal, bursal, and tracheal tissues following IBV infection and identified tissue-specific molecular markers expressed in distinct cell types. We evaluated the expression of viral RNA in diverse cellular populations and subsequently ascertained that distal tubules and collecting ducts within the kidney, bursal mucosal epithelial cells, and follicle-associated epithelial cells exhibit susceptibility to IBV infection through immunofluorescence. Furthermore, our findings revealed an upregulation in the transcription of proinflammatory cytokines IL18 and IL1B in renal macrophages as well as increased expression of apoptosis-related gene STAT in distal tubules and collecting duct cells upon IBV infection leading to renal damage. Cell-to-cell communication unveiled potential interactions between diverse cell types, as well as upregulated signaling pathways and key sender-receiver cell populations after IBV infection. Integrating single-cell data from all tissues, we applied weighted gene co-expression network analysis (WGCNA) to identify gene modules that are specifically expressed in different cell populations. Based on the WGCNA results, we identified seven immune-related gene modules and determined the differential expression pattern of module genes, as well as the hub genes within these modules. Our comprehensive data provides valuable insights into the pathogenesis of IBV as well as avian antiviral immunology.

Mutational profile evaluates metastatic capacity of Chinese colorectal cancer patients, revealed by whole-exome sequencing

Colorectal cancer (CRC) is the third most common cancer and the prevalence rate of CRC is increasing in the China. In this study, whole-exome sequencing (WES) was performed on primary tissues of 47 CRC Chinese patients including 22 metastatic and 25 non-metastatic patients. By comparison with data from western colorectal cancer patients in the Cancer Genome Atlas (TCGA), we identified a number of genes that are more likely to be mutated in Chinese colorectal cancer patients, such as MUC12, MUC12, MUC2, MUC4, HYDIN and KMT2C. Interestingly, MUC family genes including MUC12, MUC2 and MUC4, have mutation rates of >20%, while the mutation frequency was extremely low in western colorectal cancer patients, which were <3% in TCGA and 0% in Memorial Sloan Kettering Cancer Center (MSKCC). We detected metastasis-specific mutated genes including TCF7L2, MST1L, HRNR and SMAD4, while MUC4, NEB, FLG and RFPL4A alteration is more prevalent in the non-metastasis group. Further analysis reveals mutation genes in metastasis group are more focus in the Wnt and Hippo signaling pathway. APC, SMAD4 and TCF7L2 accounted for the major genetic abnormalities in this pathway. In conclusion, this study identified the unique characteristics of gene mutations in Chinese patients with colorectal cancer, and is a valuable reference for personalized treatment in Chinese CRC patients.

RNA sequencing of formalin fixed paraffin-embedded heart tissue provides transcriptomic information about chemotherapy-induced cardiotoxicity

Gene expression of formalin-fixed paraffin-embedded (FFPE) tissue may serve for molecular studies on cardiovascular diseases. Chemotherapeutics, such as doxorubicin (DOX) may cause heart injury, but the mechanisms of these side effects of DOX are not well understood. This study aimed to investigate whether DOX-induced gene expression in archival FFPE heart tissue in experimental rats would correlate with the gene expression in fresh-frozen heart tissue by applying RNA sequencing technology. The results showed RNA from FFPE samples was degraded, resulting in a lower number of uniquely mapped reads. However, DOX-induced differentially expressed genes in FFPE were related to molecular mechanisms of DOX-induced cardiotoxicity, such as inflammation, calcium binding, endothelial dysfunction, senescence, and cardiac hypertrophy signaling. Our data suggest that, despite the limitations, RNA sequencing of archival FFPE heart tissue supports utilizing FFPE tissues from retrospective studies on cardiovascular disorders, including DOX-induced cardiotoxicity.

Exploring natural killer cell-related biomarkers in multiple myeloma: a novel nature killer cell-related model predicting prognosis and immunotherapy response using single-cell study

BACKGROUND

Natural killer cells (NKs) may be involved in multiple myeloma (MM) progression. The present study elucidated the correlation between NKs and the progression of MM using single-cell binding transcriptome probes to identify NK cell-related biomarkers.

METHODS

Single-cell analysis was performed including cell and subtype annotation, cell communication, and pseudotime analysis. Hallmark pathway enrichment analysis of NKs and NKs-related differentially expressed genes (DEGs) were conducted using Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and protein-protein interaction (PPI) networks. Then, a risk model was structured based on biomarkers identified through univariate Cox regression analysis and least absolute shrinkage and selection operator regression analysis and subsequently validated. Additionally, correlation of clinical characteristics, gene set enrichment analysis, immune analysis, regulatory network, and drug forecasting were explored.

RESULTS

A total of 13 cell clusters were obtained and annotated, including 8 cell populations that consisted of NKs. Utilizing 123 PPI network node genes, 8 NK-related DEGs were selected to construct a prognostic model. Immune cell infiltration results suggested that 11 immune cells exhibited marked differences in the high and low-risk groups. Finally, the model was used to screen potential drug targets to enhance immunotherapy efficacy.

CONCLUSION

A new prognostic model for MM associated with NKs was constructed and validated. This model provides a fresh perspective for predicting patient outcomes, immunotherapeutic response, and candidate drugs.

Notum enhances gastric cancer stem-like cell properties through upregulation of Sox2 by PI3K/AKT signaling pathway

PURPOSE

Considerable evidence suggests that tumor cells with stemness features contribute to initiation, progression, recurrence of gastric cancer (GC) and resistance to therapy, but involvement of underlying regulators and mechanisms remain largely unclear. However, the clinical significance and biological function of Notum in GC tumor sphere formation and tumorigenesis remain unclear.

METHODS

Bioinformatics analysis, RT-qPCR, western blot and imunohistochemistry staining were applied to characterize Notum expression in GC specimens. The early diagnostic value of Notum was analyzed by logistic regression analysis method. Cancer stemness assays were used in Notum knockdown and overexpressing cells in vitro and in vivo. RNA-seq was employed to reveal the downstream effectors of Notum.

RESULTS

Notum is highly expressed in early stage of GC patients and stem-like GC cells. For discriminating the early-stage and advanced GC patients, the joint analysis had a better diagnostic value. Overexpression of Notum markedly increased stemness features of GC cells to promote tumor sphere formation and tumorigenesis. Conversely, Notum knockdown attenuated the stem-like cell properties in vitro and in vivo. Mechanically, Notum upregulates Sox2 through activating the PI3K/AKT signaling pathway. Notum inhibitor Caffeine exhibited a potent inhibitory effect on stemness features by impairing the PI3K/AKT signaling pathway activity and targeting Sox2.

CONCLUSION

Our findings confer a comprehensive and mechanistic function of Notum in GC tumor sphere formation and tumorigenesis that may provide a novel and promising target for early diagnosis and clinical therapy of GC.

Necroptosis-related KLRB1 was a potent tumor suppressor and immunotherapy determinant in breast cancer

Breast cancer is a multifaceted and diverse illness that impacts millions of people globally. Identifying the underlying causes of BRCA and creating efficient treatment plans are urgent. Necroptosis is widely involved in cancer development. However, the specific roles of necroptosis in cancer immunotherapy of breast cancer have not been explored. In this study, we aim to establish the connection between necroptosis and immunotherapy in BRCA. TCGA, METABRIC, GSE103091, GSE159956, and GSE96058 were included for bioinformatics analysis. NMF and CoxBoost algorithms were used to develop the necroptosis-related patterns and model, respectively. A necroptosis-related model was developed and determined KLRB1 as a critical tumor suppressor by in vitro validation. The mutation characteristics, immune characteristics, and molecular functions of KLRB1 were explored. We further examined how necroptosis-related KLRB1 functions in BRCA as a powerful tumor suppressor and regulates the activity of macrophages by in vitro validation, including CCK8, EdU, and Transwell assays. KLRB1 was also revealed to be an immunotherapy determinant.

Tropomyosin1 isoforms underlie epithelial to mesenchymal plasticity, metastatic dissemination, and resistance to chemotherapy in high-grade serous ovarian cancer

Phenotypic plasticity, defined as the ability of individual cells with stable genotypes to exert different phenotypes upon exposure to specific environmental cues, represent the quintessential hallmark of the cancer cell en route from the primary lesion to distant organ sites where metastatic colonization will occur. Phenotypic plasticity is driven by a broad spectrum of epigenetic mechanisms that allow for the reversibility of epithelial-to-mesenchymal and mesenchymal-to-epithelial transitions (EMT/MET). By taking advantage of the co-existence of epithelial and quasi-mesenchymal cells within immortalized cancer cell lines, we have analyzed the role of EMT-related gene isoforms in the regulation of epithelial mesenchymal plasticity (EMP) in high grade serous ovarian cancer. When compared with colon cancer, a distinct spectrum of downstream targets characterizes quasi-mesenchymal ovarian cancer cells, likely to reflect the different modalities of metastasis formation between these two types of malignancy, i.e. hematogenous in colon and transcoelomic in ovarian cancer. Moreover, upstream RNA-binding proteins differentially expressed between epithelial and quasi-mesenchymal subpopulations of ovarian cancer cells were identified that underlie differential regulation of EMT-related isoforms. In particular, the up- and down-regulation of RBM24 and ESRP1, respectively, represent a main regulator of EMT in ovarian cancer cells. To validate the functional and clinical relevance of our approach, we selected and functionally analyzed the Tropomyosin 1 gene (TPM1), encoding for a protein that specifies the functional characteristics of individual actin filaments in contractile cells, among the ovarian-specific downstream AS targets. The low-molecular weight Tpm1.8/9 isoforms are specifically expressed in patient-derived ascites and promote invasion through activation of EMT and Wnt signaling, together with a broad spectrum of inflammation-related pathways. Moreover, Tpm1.8/9 expression confers resistance to taxane- and platinum-based chemotherapy. Small molecule inhibitors that target the Tpm1 isoforms support targeting Tpm1.8/9 as therapeutic targets for the development of future tailor-made clinical interventions.

INX-315, a Selective CDK2 Inhibitor, Induces Cell Cycle Arrest and Senescence in Solid Tumors

Cyclin-dependent kinase 2 (CDK2) is thought to play an important role in driving proliferation of certain cancers, including those harboring CCNE1 amplification and breast cancers that have acquired resistance to CDK4/6 inhibitors (CDK4/6i). The precise impact of pharmacologic inhibition of CDK2 is not known due to the lack of selective CDK2 inhibitors. Here we describe INX-315, a novel and potent CDK2 inhibitor with high selectivity over other CDK family members. Using cell-based assays, patient-derived xenografts (PDX), and transgenic mouse models, we show that INX-315 (i) promotes retinoblastoma protein hypophosphorylation and therapy-induced senescence (TIS) in CCNE1-amplified tumors, leading to durable control of tumor growth; (ii) overcomes breast cancer resistance to CDK4/6i, restoring cell cycle control while reinstating the chromatin architecture of CDK4/6i-induced TIS; and (iii) delays the onset of CDK4/6i resistance in breast cancer by driving deeper suppression of E2F targets. Our results support the clinical development of selective CDK2 inhibitors.

SIGNIFICANCE

INX-315 is a novel, selective inhibitor of CDK2. Our preclinical studies demonstrate activity for INX-315 in both CCNE1-amplified cancers and CDK4/6i-resistant breast cancer. In each case, CDK2 inhibition induces cell cycle arrest and a phenotype resembling cellular senescence. Our data support the development of selective CDK2 inhibitors in clinical trials. See related commentary by Watts and Spencer, p. 386. This article is featured in Selected Articles from This Issue, p. 384.

A machine learning model identifies M3-like subtype in AML based on PML/RARα targets

The typical genomic feature of acute myeloid leukemia (AML) M3 subtype is the fusion event of PML/RARα, and ATRA/ATO-based combination therapy is current standard treatment regimen for M3 subtype. Here, a machine-learning model based on expressions of PML/RARα targets was developed to identify M3 patients by analyzing 1228 AML patients. Our model exhibited high accuracy. To enable more non-M3 AML patients to potentially benefit from ATRA/ATO therapy, M3-like patients were further identified. We found that M3-like patients had strong GMP features, including the expression patterns of M3 subtype marker genes, the proportion of myeloid progenitor cells, and deconvolution of AML constituent cell populations. M3-like patients exhibited distinct genomic features, low immune activity and better clinical survival. The initiative identification of patients similar to M3 subtype may help to identify more patients that would benefit from ATO/ATRA treatment and deepen our understanding of the molecular mechanism of AML pathogenesis.

Small RNAs from mitochondrial genome recombination sites are incorporated into T. gondii mitoribosomes

The mitochondrial genomes of apicomplexans comprise merely three protein-coding genes, alongside a set of thirty to forty genes encoding small RNAs (sRNAs), many of which exhibit homologies to rRNA from E. coli. The expression status and integration of these short RNAs into ribosomes remains unclear and direct evidence for active ribosomes within apicomplexan mitochondria is still lacking. In this study, we conducted small RNA sequencing on the apicomplexan Toxoplasma gondii to investigate the occurrence and function of mitochondrial sRNAs. To enhance the analysis of sRNA sequencing outcomes, we also re-sequenced the T. gondii mitochondrial genome using an improved organelle enrichment protocol and Nanopore sequencing. It has been established previously that the T. gondii genome comprises 21 sequence blocks that undergo recombination among themselves but that their order is not entirely random. The enhanced coverage of the mitochondrial genome allowed us to characterize block combinations at increased resolution. Employing this refined genome for sRNA mapping, we find that many small RNAs originated from the junction sites between protein-coding blocks and rRNA sequence blocks. Surprisingly, such block border sRNAs were incorporated into polysomes together with canonical rRNA fragments and mRNAs. In conclusion, apicomplexan ribosomes are active within polysomes and are indeed assembled through the integration of sRNAs, including previously undetected sRNAs with merged mRNA-rRNA sequences. Our findings lead to the hypothesis that T. gondii's block-based genome organization enables the dual utilization of mitochondrial sequences as both messenger RNAs and ribosomal RNAs, potentially establishing a link between the regulation of rRNA and mRNA expression.

Pick-up single-cell proteomic analysis for quantifying up to 3000 proteins in a Mammalian cell

The shotgun proteomic analysis is currently the most promising single-cell protein sequencing technology, however its identification level of ~1000 proteins per cell is still insufficient for practical applications. Here, we develop a pick-up single-cell proteomic analysis (PiSPA) workflow to achieve a deep identification capable of quantifying up to 3000 protein groups in a mammalian cell using the label-free quantitative method. The PiSPA workflow is specially established for single-cell samples mainly based on a nanoliter-scale microfluidic liquid handling robot, capable of achieving single-cell capture, pretreatment and injection under the pick-up operation strategy. Using this customized workflow with remarkable improvement in protein identification, 2449-3500, 2278-3257 and 1621-2904 protein groups are quantified in single A549 cells (n = 37), HeLa cells (n = 44) and U2OS cells (n = 27) under the DIA (MBR) mode, respectively. Benefiting from the flexible cell picking-up ability, we study HeLa cell migration at the single cell proteome level, demonstrating the potential in practical biological research from single-cell insight.

Key candidate genes for male sterility in peppers unveiled via transcriptomic and proteomic analyses

This study aimed to enhance the use of male sterility in pepper to select superior hybrid generations. Transcriptomic and proteomic analyses of fertile line 1933A and nucleic male sterility line 1933B of Capsicum annuum L. were performed to identify male sterility-related proteins and genes. The phylogenetic tree, physical and chemical characteristics, gene structure characteristics, collinearity and expression characteristics of candidate genes were analyzed. The study identified 2,357 differentially expressed genes, of which 1,145 and 229 were enriched in the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases, respectively. A total of 7,628 quantifiable proteins were identified and 29 important proteins and genes were identified. It is worth noting that the existence of CaPRX genes has been found in both proteomics and transcriptomics, and 3 CaPRX genes have been identified through association analysis. A total of 66 CaPRX genes have been identified at the genome level, which are divided into 13 subfamilies, all containing typical CaPRX gene conformal domains. It is unevenly distributed across 12 chromosomes (including the virtual chromosome Chr00). Salt stress and co-expression analysis show that male sterility genes are expressed to varying degrees, and multiple transcription factors are co-expressed with CaPRXs, suggesting that they are involved in the induction of pepper salt stress. The study findings provide a theoretical foundation for genetic breeding by identifying genes, metabolic pathways, and molecular mechanisms involved in male sterility in pepper.

The association between atherosclerosis and nonalcoholic fatty liver disease

Atherosclerosis (AS) is closely related to nonalcoholic fatty liver disease (NAFLD), which promotes and exacerbates the development of AS. However, it is uncertain how the precise underlying mechanism occurs. Here, we attempted to further explore the association underlying atherosclerosis and nonalcoholic fatty liver disease through integrated bioinformatics analysis. Microarray data for atherosclerosis and nonalcoholic fatty liver disease were retrieved from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was used to identify the genes related to atherosclerosis and nonalcoholic fatty liver disease showing co-expression. Additionally, the common gene targets associated with atherosclerosis and nonalcoholic fatty liver disease were also analyzed and screened using data from 3 public databases [comparative toxicogenomics database (CTD), DISEASES, and GeneCards]. The Gene Ontology (GO) enrichment analysis and the Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis were performed using Metascape R, respectively. The protein-protein interaction networks (PPI) network was constructed using Cytoscape. According to the results of an analysis of common genes, matrix metalloproteinase 9 (MMP9) is co-expressed up-regulated in AS and NAFLD and is enriched in inflammatory and immune-related collaterals. Consequently, MMP9 may work together through immunity and inflammation to treat AS and NAFLD and may be a potential therapeutic target in the future. The findings of this study provide new insights into the shared association between AS and NAFLD. MMP9 is co-expressed up-regulated in AS and NAFLD, which be able to reveal the presence of co-expressed genes in atherosclerosis and NAFLD.

Exploring the Biomarkers and Potential Mechanisms of Botulinum Toxin Type A in the Treatment of Microglial Inflammatory Activation through P2X7 Receptors based on Transcriptome Sequencing

AIMS

This study aims to explore the potential mechanism by which Botulinum toxin type A (BoNT/ A) inhibits microglial inflammatory activation through P2X7 receptors (P2X7R).

BACKGROUND

BoNT/A is a promising analgesic drug, and previous studies have established that it alleviates Neuropathic Pain (NP) by inhibiting microglial inflammatory activation. This study examined the biomarkers and potential mechanisms by which BoNT/A relieves neuropathic pain by mediating microglial P2X7R and analyzing transcriptome sequencing data from mouse BV-2 microglial cells.

OBJECTIVE

The P2X7R agonist Bz-ATP was used to induce microglial inflammatory activation, whilst RNAseq technology was used to explore the biomarkers and potential mechanisms through which BoNT/A suppresses microglial inflammation.

METHODS

RNA sequencing was performed on three BV-2 cell samples treated with a P2X7R specific activator (Bz-ATP) and three BV-2 cell samples pre-treated with BoNT/A. Only data that successfully passed quality control measures were included in subsequent analysis. Initially, Differentially Expressed Genes (DEGs) were identified from BoNT/A and control samples, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Biomarkers were then identified by constructing a Protein- Protein Interaction (PPI) network and utilizing the CytoHubba plug-in in Cytoscape software. Lastly, enrichment analysis and regulatory network analysis were performed to elucidate the potential mechanism of BoNT/A in the treatment of NP.

RESULTS

93 DEGs related to the "cell component size regulation" GO term and enriched in the "axon guidance" KEGG pathway were identified. Subsequently, 6 biomarkers were identified, namely PTPRF, CHDH, CKM, Ky, Sema3b, and Sema3f, which were enriched in pathways related to biosynthesis and metabolism, disease progression, signal transduction, and organelle function, including the "ribosome" and "Wnt signaling pathway." Finally, a competing endogenous RNA (ceRNAs) network was constructed from 6 mRNAs, 66 miRNAs, and 31 lncRNAs, forming a complex relationship network.

CONCLUSION

Six genes (PTPRF, Sema3b, Sema3f, CHDH, CKM, and Ky) were identified as biomarkers of microglial inflammatory activation following BoNT/A treatment. This finding may provide a valuable reference for the relief and treatment of neuropathic pain.

Detection of neoplastic-immune hybrid cells with metastatic properties in uveal melanoma

Background

Uveal melanoma is the most common non-cutaneous melanoma and is an intraocular malignancy affecting nearly 7,000 individuals per year worldwide. Of these, approximately 50% will progress to metastatic disease for which there are currently no effective therapies. Despite advances in molecular profiling and metastatic stratification of uveal melanoma tumors, little is known regarding their underlying biology of metastasis. Our group has identified a disseminated neoplastic cell population characterized by co-expression of immune and melanoma proteins, circulating hybrid cells (hybrids), in patients with uveal melanoma. Compared to circulating tumor cells, which lack expression of immune proteins, hybrids are detected at an increased prevalence in peripheral blood and can be used as a non-invasive biomarker to predict metastatic progression.

Methods

To ascertain mechanisms underlying enhanced hybrid cell dissemination we identified hybrid cells within primary uveal melanoma tumors using single cell RNA sequencing and evaluated their gene expression and predicted ligand-receptor interactions in relation to other melanoma and immune cells within the primary tumor. We then verified expression of upregulated hybrid pathways within patient-matched tumor and peripheral blood hybrids using cyclic immunofluorescence and quantified their protein expression relative to other non-hybrid tumor and disseminated tumor cells.

Results

Among the top upregulated genes and pathways in hybrid cells were those involved in enhanced cell motility and cytoskeletal rearrangement, immune evasion, and altered cellular metabolism. In patient-matched tumor and peripheral blood, we verified gene expression by examining concordant protein expression for each pathway category: TMSB10 (cell motility), CD74 (immune evasion) and GPX1 (metabolism). Both TMSB10 and GPX1 were expressed on significantly higher numbers of disseminated hybrid cells compared to circulating tumor cells, and CD74 and GPX1 were expressed on more disseminated hybrids than tumor-resident hybrids. Lastly, we identified that hybrid cells express ligand-receptor signaling pathways implicated in promoting metastasis including GAS6-AXL, CXCL12-CXCR4, LGALS9-P4HB and IGF1-IGFR1.

Conclusion

These findings highlight the importance of TMSB10, GPX1 and CD74 for successful hybrid cell dissemination and survival in circulation. Our results contribute to the understanding of uveal melanoma tumor progression and interactions between tumor cells and immune cells in the tumor microenvironment that may promote metastasis.

Fn14 promotes myoblast fusion during regenerative myogenesis

Skeletal muscle regeneration involves coordinated activation of an array of signaling pathways. Fibroblast growth factor-inducible 14 (Fn14) is a bona fide receptor for the TWEAK cytokine. Levels of Fn14 are increased in the skeletal muscle of mice after injury. However, the cell-autonomous role of Fn14 in muscle regeneration remains unknown. Here, we demonstrate that global deletion of the Fn14 receptor in mice attenuates muscle regeneration. Conditional ablation of Fn14 in myoblasts but not in differentiated myofibers of mice inhibits skeletal muscle regeneration. Fn14 promotes myoblast fusion without affecting the levels of myogenic regulatory factors in the regenerating muscle. Fn14 deletion in myoblasts hastens initial differentiation but impairs their fusion. The overexpression of Fn14 in myoblasts results in the formation of myotubes having an increased diameter after induction of differentiation. Ablation of Fn14 also reduces the levels of various components of canonical Wnt and calcium signaling both in vitro and in vivo. Forced activation of Wnt signaling rescues fusion defects in Fn14-deficient myoblast cultures. Collectively, our results demonstrate that Fn14-mediated signaling positively regulates myoblast fusion and skeletal muscle regeneration.

Range of chromatin accessibility configurations are permissive of GABAergic fate acquisition in developing mouse brain

In recent single-cell -omics studies, both the differential activity of transcription factors regulating cell fate determination and differential genome activation have been tested for utility as descriptors of cell types. Naturally, genome accessibility and gene expression are interlinked. To understand the variability in genomic feature activation in the GABAergic neurons of different spatial origins, we have mapped accessible chromatin regions and mRNA expression in single cells derived from the developing mouse central nervous system (CNS). We first defined a reference set of open chromatin regions for scATAC-seq read quantitation across samples, allowing comparison of chromatin accessibility between brain regions and cell types directly. Second, we integrated the scATAC-seq and scRNA-seq data to form a unified resource of transcriptome and chromatin accessibility landscape for the cell types in di- and telencephalon, midbrain and anterior hindbrain of E14.5 mouse embryo. Importantly, we implemented resolution optimization at the clustering, and automatized the cell typing step. We show high level of concordance between the cell clustering based on the chromatin accessibility and the transcriptome in analyzed neuronal lineages, indicating that both genome and transcriptome features can be used for cell type definition. Hierarchical clustering by the similarity in accessible chromatin reveals that the genomic feature activation correlates with neurotransmitter phenotype, selector gene expression, cell differentiation stage and neuromere origins.

Identifying individualized prognostic signature and unraveling the molecular mechanism of recurrence in early-onset colorectal cancer

BACKGROUND

The incidence and mortality of early-onset colorectal cancer (EOCRC; < 50 years old) is increasing worldwide, with a high recurrence rate. The inherent heterogeneity of EOCRC makes its treatment challenging. Hence, to further understand the biology and reveal the molecular mechanisms of EOCRC, a recurrence risk signature is needed to guide clinical management.

METHODS

Based on the relative expression orderings (REOs) of genes in each sample, a prognostic signature was developed and validated utilizing multiple independent datasets. The underlying molecular mechanisms between distinct prognostic groups were explored via integrative analysis of multi-omics data.

RESULTS

The prognostic signature consisting of 6 gene pairs (6-GPS) could predict the recurrence risk for EOCRC at the individual level. High-risk EOCRC classified by 6-GPS showed a poor prognosis but a good response to adjuvant chemotherapy. Moreover, high-risk EOCRC was characterized by epithelial-mesenchymal transition (EMT) and enriched angiogenesis, and had higher mutation burden, immune cell infiltration, and PD-1/PD-L1 expression. Furthermore, we identified four genes associated with relapse-free survival in EOCRC, including SERPINE1, PECAM1, CDH1, and ANXA1. They were consistently differentially expressed at the transcriptome and proteome levels between high-risk and low-risk EOCRCs. They were also involved in regulating cancer progression and immune microenvironment in EOCRC. Notably, the expression of SERPINE1 and ANXA1 positively correlated with M2-like macrophage infiltration.

CONCLUSION

Our results indicate that 6-GPS can robustly predict the recurrence risk of EOCRC, and that SERPINE1, PECAM1, CDH1, and ANXA1 may serve as potential therapeutic targets. This study provides valuable information for the precision treatment of EOCRC.

Crowdsourcing temporal transcriptomic coronavirus host infection data: Resources, guide, and novel insights

The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) reawakened the need to rapidly understand the molecular etiologies, pandemic potential, and prospective treatments of infectious agents. The lack of existing data on SARS-CoV-2 hampered early attempts to treat severe forms of coronavirus disease-2019 (COVID-19) during the pandemic. This study coupled existing transcriptomic data from severe acute respiratory syndrome-related coronavirus 1 (SARS-CoV-1) lung infection animal studies with crowdsourcing statistical approaches to derive temporal meta-signatures of host responses during early viral accumulation and subsequent clearance stages. Unsupervised and supervised machine learning approaches identified top dysregulated genes and potential biomarkers (e.g. CXCL10, BEX2, and ADM). Temporal meta-signatures revealed distinct gene expression programs with biological implications to a series of host responses underlying sustained Cxcl10 expression and Stat signaling. Cell cycle switched from G1/G0 phase genes, early in infection, to a G2/M gene signature during late infection that correlated with the enrichment of DNA damage response and repair genes. The SARS-CoV-1 meta-signatures were shown to closely emulate human SARS-CoV-2 host responses from emerging RNAseq, single cell, and proteomics data with early monocyte-macrophage activation followed by lymphocyte proliferation. The circulatory hormone adrenomedullin was observed as maximally elevated in elderly patients who died from COVID-19. Stage-specific correlations to compounds with potential to treat COVID-19 and future coronavirus infections were in part validated by a subset of twenty-four that are in clinical trials to treat COVID-19. This study represents a roadmap to leverage existing data in the public domain to derive novel molecular and biological insights and potential treatments to emerging human pathogens.

Comparative Genomic Analysis of Asian Cultivated Rice and Its Wild Progenitor (Oryza rufipogon) Has Revealed Evolutionary Innovation of the Pentatricopeptide Repeat Gene Family through Gene Duplication

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in land plants. However, current knowledge about the evolution of the PPR gene family remains largely limited. In this study, we performed a comparative genomic analysis of the PPR gene family in O. sativa and its wild progenitor, O. rufipogon, and outlined a comprehensive landscape of gene duplications. Our findings suggest that the majority of PPR genes originated from dispersed duplications. Although segmental duplications have only expanded approximately 11.30% and 13.57% of the PPR gene families in the O. sativa and O. rufipogon genomes, we interestingly obtained evidence that segmental duplication promotes the structural diversity of PPR genes through incomplete gene duplications. In the O. sativa and O. rufipogon genomes, 10 (~33.33%) and 22 pairs of gene duplications (~45.83%) had non-PPR paralogous genes through incomplete gene duplication. Segmental duplications leading to incomplete gene duplications might result in the acquisition of domains, thus promoting functional innovation and structural diversification of PPR genes. This study offers a unique perspective on the evolution of PPR gene structures and underscores the potential role of segmental duplications in PPR gene structural diversity.

Analysis of uveal melanoma scRNA sequencing data identifies neoplastic-immune hybrid cells that exhibit metastatic potential

Uveal melanoma (UM) is the most common non-cutaneous melanoma and is an intraocular malignancy that affects nearly 7,000 individuals per year worldwide. Of these, nearly 50% will progress to metastatic disease for which there are currently no effective therapies. Despite advances in the molecular profiling and metastatic stratification of class 1 and 2 UM tumors, little is known regarding the underlying biology of UM metastasis. Our group has identified a disseminated tumor cell population characterized by co-expression of immune and melanoma proteins, (circulating hybrid cells (CHCs), in patients with UM. Compared to circulating tumor cells, CHCs are detected at an increased prevalence in peripheral blood and can be used as a non-invasive biomarker to predict metastatic progression. To identify mechanisms underlying enhanced hybrid cell dissemination we sought to identify hybrid cells within a primary UM single cell RNA-seq dataset. Using rigorous doublet discrimination approaches, we identified UM hybrids and evaluated their gene expression, predicted ligand-receptor status, and cell-cell communication state in relation to other melanoma and immune cells within the primary tumor. We identified several genes and pathways upregulated in hybrid cells, including those involved in enhancing cell motility and cytoskeleton rearrangement, evading immune detection, and altering cellular metabolism. In addition, we identified that hybrid cells express ligand-receptor signaling pathways implicated in promoting cancer metastasis including IGF1-IGFR1, GAS6-AXL, LGALS9-P4HB, APP-CD74 and CXCL12-CXCR4. These results contribute to our understanding of tumor progression and interactions between tumor cells and immune cells in the UM microenvironment that may promote metastasis.

Proton Pump Inhibitors Modulate Gene Expression Profile in Esophageal Mucosa and Microbiome

OBJECTIVE

Proton pump inhibitors (PPIs) are commonly used to manage children with upper gastrointestinal symptoms and without a formal diagnosis. We investigated the effect of PPIs on esophageal mucosal transcriptome and active microbiota in children with normal esophagi. Furthermore, we examined whether the differences in host esophageal mucosal gene expression were driven by an underlying esophageal epithelial cell type composition.

METHODS

Using metatranscriptomics, the host transcriptional and active microbial profiles were captured from 17 esophageal biopsy samples (PPI naïve [PPI-], n = 7; PPI exposed [PPI+], n = 10) collected from children without any endoscopic and histologic abnormalities in their esophagus (normal esophagus). Deconvolution computational analysis was performed with xCell to assess if the observed epithelial gene expression changes were related to the cell type composition in the esophageal samples.

RESULTS

The median (IQR) age of our cohort was 14 years (12-16) with female (63%) preponderance. Both groups were similar in terms of their demographics and clinical features. Compared with PPI-, the PPI+ had upregulation of 27 genes including the MUC genes. The cell type composition was similar between the PPI- and PPI+ groups. Prevotella sp and Streptococcus sp were abundant in PPI+ group.

CONCLUSIONS

In children with normal esophagus, PPI exposure can be associated with upregulation of esophageal mucosal homeostasis and epithelial cell function genes in a cell-type independent manner, and an altered esophageal microbiome. Additional studies are warranted to validate our findings and to investigate the causal effect of PPIs on the normal esophageal epithelium and microbial communities.