Cell surface RNAs control neutrophil recruitment.

RNAs localizing to the outer cell surface have been recently identified in mammalian cells, including RNAs with glycan modifications known as glycoRNAs. However, the functional significance of cell surface RNAs and their production are poorly known. We report that cell surface RNAs are critical for neutrophil recruitment and that the mammalian homologs of the sid-1 RNA transporter are required for glycoRNA expression. Cell surface RNAs can be readily detected in murine neutrophils, the elimination of which substantially impairs neutrophil recruitment to inflammatory sites in vivo and reduces neutrophils' adhesion to and migration through endothelial cells. Neutrophil glycoRNAs are predominantly on cell surface, important for neutrophil-endothelial interactions, and can be recognized by P-selectin (Selp). Knockdown of the murine Sidt genes abolishes neutrophil glycoRNAs and functionally mimics the loss of cell surface RNAs. Our data demonstrate the biological importance of cell surface glycoRNAs and highlight a noncanonical dimension of RNA-mediated cellular functions.

Colorectal Cancer Cells Enter a Diapause-like DTP State to Survive Chemotherapy.

Cancer cells enter a reversible drug-tolerant persister (DTP) state to evade death from chemotherapy and targeted agents. It is increasingly appreciated that DTPs are important drivers of therapy failure and tumor relapse. We combined cellular barcoding and mathematical modeling in patient-derived colorectal cancer models to identify and characterize DTPs in response to chemotherapy. Barcode analysis revealed no loss of clonal complexity of tumors that entered the DTP state and recurred following treatment cessation. Our data fit a mathematical model where all cancer cells, and not a small subpopulation, possess an equipotent capacity to become DTPs. Mechanistically, we determined that DTPs display remarkable transcriptional and functional similarities to diapause, a reversible state of suspended embryonic development triggered by unfavorable environmental conditions. Our study provides insight into how cancer cells use a developmentally conserved mechanism to drive the DTP state, pointing to novel therapeutic opportunities to target DTPs.

DNA-Demethylating Agents Target Colorectal Cancer Cells by Inducing Viral Mimicry by Endogenous Transcripts.

DNA-demethylating agents have shown clinical anti-tumor efficacy via an unknown mechanism of action. Using a combination of experimental and bioinformatics analyses in colorectal cancer cells, we demonstrate that low-dose 5-AZA-CdR targets colorectal cancer-initiating cells (CICs) by inducing viral mimicry. This is associated with induction of dsRNAs derived at least in part from endogenous retroviral elements, activation of the MDA5/MAVS RNA recognition pathway, and downstream activation of IRF7. Indeed, disruption of virus recognition pathways, by individually knocking down MDA5, MAVS, or IRF7, inhibits the ability of 5-AZA-CdR to target colorectal CICs and significantly decreases 5-AZA-CdR long-term growth effects. Moreover, transfection of dsRNA into CICs can mimic the effects of 5-AZA-CdR. Together, our results represent a major shift in understanding the anti-tumor mechanisms of DNA-demethylating agents and highlight the MDA5/MAVS/IRF7 pathway as a potentially druggable target against CICs.

Kled: an ultra-fast and sensitive structural variant detection tool for long-read sequencing data.

Structural Variants (SVs) are a crucial type of genetic variant that can significantly impact phenotypes. Therefore, the identification of SVs is an essential part of modern genomic analysis. In this article, we present kled, an ultra-fast and sensitive SV caller for long-read sequencing data given the specially designed approach with a novel signature-merging algorithm, custom refinement strategies and a high-performance program structure. The evaluation results demonstrate that kled can achieve optimal SV calling compared to several state-of-the-art methods on simulated and real long-read data for different platforms and sequencing depths. Furthermore, kled excels at rapid SV calling and can efficiently utilize multiple Central Processing Unit (CPU) cores while maintaining low memory usage. The source code for kled can be obtained from https://github.com/CoREse/kled.

Epigenetic therapy induces transcription of inverted SINEs and ADAR1 dependency.

Cancer therapies that target epigenetic repressors can mediate their effects by activating retroelements within the human genome. Retroelement transcripts can form double-stranded RNA (dsRNA) that activates the MDA5 pattern recognition receptor1-6. This state of viral mimicry leads to loss of cancer cell fitness and stimulates innate and adaptive immune responses7,8. However, the clinical efficacy of epigenetic therapies has been limited. To find targets that would synergize with the viral mimicry response, we sought to identify the immunogenic retroelements that are activated by epigenetic therapies. Here we show that intronic and intergenic SINE elements, specifically inverted-repeat Alus, are the major source of drug-induced immunogenic dsRNA. These inverted-repeat Alus are frequently located downstream of 'orphan' CpG islands9. In mammals, the ADAR1 enzyme targets and destabilizes inverted-repeat Alu dsRNA10, which prevents activation of the MDA5 receptor11. We found that ADAR1 establishes a negative-feedback loop, restricting the viral mimicry response to epigenetic therapy. Depletion of ADAR1 in patient-derived cancer cells potentiates the efficacy of epigenetic therapy, restraining tumour growth and reducing cancer initiation. Therefore, epigenetic therapies trigger viral mimicry by inducing a subset of inverted-repeats Alus, leading to an ADAR1 dependency. Our findings suggest that combining epigenetic therapies with ADAR1 inhibitors represents a promising strategy for cancer treatment.

5-Fluorouracil efficacy requires anti-tumor immunity triggered by cancer-cell-intrinsic STING.

5-Fluorouracil (5-FU) is a widely used chemotherapeutic drug, but the mechanisms underlying 5-FU efficacy in immunocompetent hosts in vivo remain largely elusive. Through modeling 5-FU response of murine colon and melanoma tumors, we report that effective reduction of tumor burden by 5-FU is dependent on anti-tumor immunity triggered by the activation of cancer-cell-intrinsic STING. While the loss of STING does not induce 5-FU resistance in vitro, effective 5-FU responsiveness in vivo requires cancer-cell-intrinsic cGAS, STING, and subsequent type I interferon (IFN) production, as well as IFN-sensing by bone-marrow-derived cells. In the absence of cancer-cell-intrinsic STING, a much higher dose of 5-FU is needed to reduce tumor burden. 5-FU treatment leads to increased intratumoral T cells, and T-cell depletion significantly reduces the efficacy of 5-FU in vivo. In human colorectal specimens, higher STING expression is associated with better survival and responsiveness to chemotherapy. Our results support a model in which 5-FU triggers cancer-cell-initiated anti-tumor immunity to reduce tumor burden, and our findings could be harnessed to improve therapeutic effectiveness and toxicity for colon and other cancers.

THItoGene: a deep learning method for predicting spatial transcriptomics from histological images.

Spatial transcriptomics unveils the complex dynamics of cell regulation and transcriptomes, but it is typically cost-prohibitive. Predicting spatial gene expression from histological images via artificial intelligence offers a more affordable option, yet existing methods fall short in extracting deep-level information from pathological images. In this paper, we present THItoGene, a hybrid neural network that utilizes dynamic convolutional and capsule networks to adaptively sense potential molecular signals in histological images for exploring the relationship between high-resolution pathology image phenotypes and regulation of gene expression. A comprehensive benchmark evaluation using datasets from human breast cancer and cutaneous squamous cell carcinoma has demonstrated the superior performance of THItoGene in spatial gene expression prediction. Moreover, THItoGene has demonstrated its capacity to decipher both the spatial context and enrichment signals within specific tissue regions. THItoGene can be freely accessed at https://github.com/yrjia1015/THItoGene.

HNetGO: protein function prediction via heterogeneous network transformer.

Protein function annotation is one of the most important research topics for revealing the essence of life at molecular level in the post-genome era. Current research shows that integrating multisource data can effectively improve the performance of protein function prediction models. However, the heavy reliance on complex feature engineering and model integration methods limits the development of existing methods. Besides, models based on deep learning only use labeled data in a certain dataset to extract sequence features, thus ignoring a large amount of existing unlabeled sequence data. Here, we propose an end-to-end protein function annotation model named HNetGO, which innovatively uses heterogeneous network to integrate protein sequence similarity and protein-protein interaction network information and combines the pretraining model to extract the semantic features of the protein sequence. In addition, we design an attention-based graph neural network model, which can effectively extract node-level features from heterogeneous networks and predict protein function by measuring the similarity between protein nodes and gene ontology term nodes. Comparative experiments on the human dataset show that HNetGO achieves state-of-the-art performance on cellular component and molecular function branches.

MCDHGN: heterogeneous network-based cancer driver gene prediction and interpretability analysis.

MOTIVATION: Accurately predicting the driver genes of cancer is of great significance for carcinogenesis progress research and cancer treatment. In recent years, more and more deep-learning-based methods have been used for predicting cancer driver genes. However, deep-learning algorithms often have black box properties and cannot interpret the output results. Here, we propose a novel cancer driver gene mining method based on heterogeneous network meta-paths (MCDHGN), which uses meta-path aggregation to enhance the interpretability of predictions. RESULTS: MCDHGN constructs a heterogeneous network by using several types of multi-omics data that are biologically linked to genes. And the differential probabilities of SNV, DNA methylation, and gene expression data between cancerous tissues and normal tissues are extracted as initial features of genes. Nine meta-paths are manually selected, and the representation vectors obtained by aggregating information within and across meta-path nodes are used as new features for subsequent classification and prediction tasks. By comparing with eight homogeneous and heterogeneous network models on two pan-cancer datasets, MCDHGN has better performance on AUC and AUPR values. Additionally, MCDHGN provides interpretability of predicted cancer driver genes through the varying weights of biologically meaningful meta-paths. AVAILABILITY AND IMPLEMENTATION: https://github.com/1160300611/MCDHGN.

DVA: predicting the functional impact of single nucleotide missense variants.

BACKGROUND: In the past decade, single nucleotide variants (SNVs) have been identified as having a significant relationship with the development and treatment of diseases. Among them, prioritizing missense variants for further functional impact investigation is an essential challenge in the study of common disease and cancer. Although several computational methods have been developed to predict the functional impacts of variants, the predictive ability of these methods is still insufficient in the Mendelian and cancer missense variants. RESULTS: We present a novel prediction method called the disease-related variant annotation (DVA) method that predicts the effect of missense variants based on a comprehensive feature set of variants, notably, the allele frequency and protein-protein interaction network feature based on graph embedding. Benchmarked against datasets of single nucleotide missense variants, the DVA method outperforms the state-of-the-art methods by up to 0.473 in the area under receiver operating characteristic curve. The results demonstrate that the proposed method can accurately predict the functional impact of single nucleotide missense variants and substantially outperforms existing methods. CONCLUSIONS: DVA is an effective framework for identifying the functional impact of disease missense variants based on a comprehensive feature set. Based on different datasets, DVA shows its generalization ability and robustness, and it also provides innovative ideas for the study of the functional mechanism and impact of SNVs.

Tumour-derived exosome SNHG17 induced by oestrogen contributes to ovarian cancer progression via the CCL13-CCR2-M2 macrophage axis.

Oestrogen is known to be strongly associated with ovarian cancer. There was much work to show the importance of lncRNA SNHG17 in ovarian cancer. However, no study has revealed the molecular regulatory mechanism and functional effects between oestrogen and SNHG17 in the development and metastasis of ovarian cancer. In this study, we found that SNHG17 expression was significantly increased in ovarian cancer and positively correlated with oestrogen treatment. Oestrogen could promote M2 macrophage polarization as well as ovarian cancer cells SKOV3 and ES2 cell exosomal SNHG17 expression. When exposure to oestrogen, exosomal SNHG17 promoted ovarian cancer cell proliferation, migration, invasion and epithelial-mesenchymal transition (EMT) in vitro, and tumour growth and lung metastasis in vivo by accelerating M2-like phenotype of macrophages. Mechanically, exosomal SNHG17 could facilitate the release of CCL13 from M2 macrophage via the PI3K-Akt signalling pathway. Moreover, CCL13-CCR2 axis was identified to be involved in ovarian cancer tumour behaviours driven by oestrogen. There results demonstrate a novel mechanism that exosomal SNHG17 exerts an oncogenic effect on ovarian cancer via the CCL13-CCR2-M2 macrophage axis upon oestrogen treatment, of which SNHG17 may be a potential biomarker and therapeutic target for ovarian cancer responded to oestrogen.

N6-Methyladenosine-Modified KREMEN2 Promotes Tumorigenesis and Malignant Progression of High-Grade Serous Ovarian Cancer.

High-grade serous ovarian cancer (HGSOC) remains the most lethal female cancer by far. Herein, clinical HGSOC samples had higher N6-methyladenosine (m6A) modification than normal ovarian tissue, and its dysregulation had been reported to drive aberrant transcription and translation programs. However, Kringle-containing transmembrane protein 2 (KREMEN2) and its m6A modification have not been fully elucidated in HGSOC. In this study, the data from the high-throughput messenger RNA (mRNA) sequencing of clinical samples were processed using the weighted correlation network analysis and functional enrichment analysis. Results revealed that KREMEN2 was a driver gene in the tumorigenesis of HGSOC and a potential target of m6A demethylase fat-mass and obesity-associated protein (FTO). KREMEN2 and FTO levels were upregulated and downregulated, respectively, and correlation analysis showed a significant negative correlation in HGSOC samples. Importantly, upregulated KREMEN2 was remarkably associated with lymph node metastasis, distant metastasis, peritoneal metastasis, and high International Federation of Gynecology and Obstetrics stage (Ⅲ/Ⅳ), independent of the age of patients. KREMEN2 promoted the growth of HGSOC in vitro and in vivo, which was dependent on FTO. The methylated RNA immunoprecipitation qPCR and RNA immunoprecipitation assays were performed to verify the m6A level and sites of KREMEN2. FTO overexpression significantly decreased m6A modification in the 3' and 5' untranslated regions of KREMEN2 mRNA and downregulated its expression. In addition, we found that FTO-mediated m6A modification of KREMEN2 mRNA was recognized and stabilized by the m6A reader IGF2BP1 rather than by IGF2BP2 or IGF2BP3. This study highlights the m6A modification of KREMEN2 and extends the importance of RNA epigenetics in HGSOC.

A Microfluidic Chip-Based Automated System for Whole-Course Monitoring the Drug Responses of Organoids.

Tumor patients-derived organoids, as a promising preclinical prediction model, have been utilized to evaluate ex vivo drug responses for formulating optimal therapeutic strategies. Detecting adenosine triphosphate (ATP) has been widely used in existing organoid-based drug response tests. However, all commercial ATP detection kits containing the cell lysis procedure can only be applied for single time point ATP detection, resulting in the neglect of dynamic ATP variations in living cells. Meanwhile, due to the limited number of viable organoids from a single patient, it is impractical to exhaustively test all potential time points in search of optimal ones. In this work, a multifunctional microfluidic chip was developed to perform all procedures of organoid-based drug response tests, including establishment, culturing, drug treatment, and ATP monitoring of organoids. An ATP sensor was developed to facilitate the first successful attempt on whole-course monitoring the growth status of fragile organoids. To realize a clinically applicable automatic system for the drug testing of lung cancer, a microfluidic chip based automated system was developed to perform entire organoid-based drug response test, bridging the gap between laboratorial manipulation and clinical practices, as it outperformed previous methods by improving data repeatability, eliminating human error/sample loss, and more importantly, providing a more accurate and comprehensive evaluation of drug effects.

Clonal expansion of shared T cell receptors reveals the existence of immune commonality among different lesions of synchronous multiple primary lung cancer.

The increase in the detection rate of synchronous multiple primary lung cancer (MPLC) has posed remarkable clinical challenges due to the limited understanding of its pathogenesis and molecular features. Here, comprehensive comparisons of genomic and immunologic features between MPLC and solitary lung cancer nodule (SN), as well as different lesions of the same patient, were performed. Compared with SN, MPLC displayed a lower rate of EGFR mutation but higher rates of BRAF, MAP2K1, and MTOR mutation, which function exactly in the upstream and downstream of the same signaling pathway. Considerable heterogeneity in T cell receptor (TCR) repertoire exists among not only different patients but also among different lesions of the same patient. Invasive lesions of MPLC exhibited significantly higher TCR diversity and lower TCR expansion than those of SN. Intriguingly, different lesions of the same patient always shared a certain proportion of TCR clonotypes. Significant clonal expansion could be observed in shared TCR clonotypes, particularly in those existing in all lesions of the same patient. In conclusion, this study provided evidences of the distinctive mutational landscape, activation of oncogenic signaling pathways, and TCR repertoire in MPLC as compared with SN. The significant clonal expansion of shared TCR clonotypes demonstrated the existence of immune commonality among different lesions of the same patient and shed new light on the individually tailored precision therapy for MPLC.

Enhancing discoveries of molecular QTL studies with small sample size using summary statistic imputation.

Quantitative trait locus (QTL) analyses of multiomic molecular traits, such as gene transcription (eQTL), DNA methylation (mQTL) and histone modification (haQTL), have been widely used to infer the functional effects of genome variants. However, the QTL discovery is largely restricted by the limited study sample size, which demands higher threshold of minor allele frequency and then causes heavy missing molecular trait-variant associations. This happens prominently in single-cell level molecular QTL studies because of sample availability and cost. It is urgent to propose a method to solve this problem in order to enhance discoveries of current molecular QTL studies with small sample size. In this study, we presented an efficient computational framework called xQTLImp to impute missing molecular QTL associations. In the local-region imputation, xQTLImp uses multivariate Gaussian model to impute the missing associations by leveraging known association statistics of variants and the linkage disequilibrium (LD) around. In the genome-wide imputation, novel procedures are implemented to improve efficiency, including dynamically constructing a reused LD buffer, adopting multiple heuristic strategies and parallel computing. Experiments on various multiomic bulk and single-cell sequencing-based QTL datasets have demonstrated high imputation accuracy and novel QTL discovery ability of xQTLImp. Finally, a C++ software package is freely available at https://github.com/stormlovetao/QTLIMP.

Struct2GO: protein function prediction based on graph pooling algorithm and AlphaFold2 structure information.

MOTIVATION: In recent years, there has been a breakthrough in protein structure prediction, and the AlphaFold2 model of the DeepMind team has improved the accuracy of protein structure prediction to the atomic level. Currently, deep learning-based protein function prediction models usually extract features from protein sequences and combine them with protein-protein interaction networks to achieve good results. However, for newly sequenced proteins that are not in the protein-protein interaction network, such models cannot make effective predictions. To address this, this article proposes the Struct2GO model, which combines protein structure and sequence data to enhance the precision of protein function prediction and the generality of the model. RESULTS: We obtain amino acid residue embeddings in protein structure through graph representation learning, utilize the graph pooling algorithm based on a self-attention mechanism to obtain the whole graph structure features, and fuse them with sequence features obtained from the protein language model. The results demonstrate that compared with the traditional protein sequence-based function prediction model, the Struct2GO model achieves better results. AVAILABILITY AND IMPLEMENTATION: The data underlying this article are available at https://github.com/lyjps/Struct2GO.

Serum Response Factor Reduces Gene Expression Noise and Confers Cell State Stability.

The role of serum response factor (Srf), a central mediator of actin dynamics and mechanical signaling, in cell identity regulation is debated to be either a stabilizer or a destabilizer. We investigated the role of Srf in cell fate stability using mouse pluripotent stem cells. Despite the fact that serum-containing cultures yield heterogeneous gene expression, deletion of Srf in mouse pluripotent stem cells leads to further exacerbated cell state heterogeneity. The exaggerated heterogeneity is detectible not only as increased lineage priming but also as the developmentally earlier 2C-like cell state. Thus, pluripotent cells explore more variety of cellular states in both directions of development surrounding naïve pluripotency, a behavior that is constrained by Srf. These results support that Srf functions as a cell state stabilizer, providing rationale for its functional modulation in cell fate intervention and engineering.

miR-206 alleviates LPS-induced inflammatory injury in cardiomyocytes via directly targeting USP33 to inhibit the JAK2/STAT3 signaling pathway.

Previous reports have confirmed that miR-206 participates in inflammatory cardiomyopathy, but its definite mechanism remains elusive. This study aims to elucidate the potential mechanism of miR-206 in septic cardiomyopathy (SCM). The primary mouse cardiomyocytes were isolated and exposed to lipopolysaccharides (LPS) to construct a septic injury model in vitro. Then, the gene transcripts and protein levels were detected by RT-qPCR and/or Western blot assay. Cell proliferation, apoptosis, and inflammatory responses were evaluated by CCK-8/EdU, flow cytometry, and ELISA assays, respectively. Dual luciferase assay, Co-IP, and ubiquitination experiments were carried out to validate the molecular interactions among miR-206, USP33, and JAK2/STAT3 signaling. miR-206 was significantly downregulated, but USP33 was upregulated in LPS-induced cardiomyocytes. Gain-of-function of miR-206 elevated the proliferation but suppressed the inflammatory responses and apoptosis in LPS-induced cardiomyocytes. USP33, as a member of the USP protein family, was confirmed to be a direct target of miR-206 and could catalyze deubiquitination of JAK2 to activate JAK2/STAT3 signaling. Rescue experiments presented that neither upregulation of USP33 nor JAK2/STAT3 signaling activation considerably reversed the protective effects of miR-206 upregulation in LPS-induced cardiomyocytes. The above data showed that miR-206 protected cardiomyocytes from LPS-induced inflammatory injuries by targeting the USP33/JAK2/STAT3 signaling pathway, which might be a novel target for SCM treatment.

Risk factors for short-term prognosis of end-stage liver disease complicated by invasive pulmonary aspergillosis.

BACKGROUND AND AIM: Patients with end-stage liver disease (ESLD) are susceptible to invasive pulmonary aspergillosis (IPA). This study aimed to investigate the risk factors affecting the occurrence and short-term prognosis of ESLD complicated by IPA. METHODS: This retrospective case-control study included 110 patients with ESLD. Of them, 27 ESLD-IPA received antifungal therapy with amphotericin B (AmB); 27 AmB-free-treated ESLD-IPA patients were enrolled through 1:1 propensity score matching. Fifty-six ESLD patients with other comorbid pulmonary infections were enrolled as controls. The basic features of groups were compared, while the possible risk factors affecting the occurrence and short-term outcomes of IPA were analyzed. RESULTS: Data analysis revealed invasive procedures, glucocorticoid exposure, and broad-spectrum antibiotic use were independent risk factors for IPA. The 54 patients with ESLD-IPA exhibited an overall treatment effectiveness and 28-d mortality rate of 50.00% and 20.37%, respectively, in whom patients treated with AmB-containing showed higher treatment efficacy than patients treated with AmB-free antifungal regimens (66.7% vs. 33.3%, respectively, χ2 = 6.000, P = 0.014). Multivariate logistic regression analysis revealed that the treatment regimen was the only predictor affecting patient outcomes, with AmB-containing regimens were 4.893 times more effective than AmB-free regimens (95% CI, 1.367-17.515; P = 0.015). The only independent predictors affecting the 28-d mortality rate were neutrophil-to-lymphocyte ratio and IPA diagnosis (OR = 1.140 and 10.037, P = 0.046 and 0.025, respectively). CONCLUSIONS: Glucocorticoid exposure, invasive procedures, and broad-spectrum antibiotic exposure increased the risk of IPA in ESLD patients. AmB alone or combined with other antifungals may serve as an economical, safe, and effective treatment option for ESLD-IPA.

GCNT3 Promotes Hepatocellular Carcinoma Progression and EMT by Activating the PI3K/AKT Pathway.

Primary liver cancer, specifically hepatocellular carcinoma (HCC), is a major global health concern. GCNT3 has been identified as an oncogene in various human malignancies. This investigation aimed to discover the GCNT3 function in HCC. The present study employed integrated bioinformatics analyses to assess the expression pattern, prognostic implications, and putative function of GCNT3 in HCC. Transwell flow cytometry, CCK-8, and wound healing assays were performed to examine HCC cell growth, cell cycle, apoptosis, invasion, and migration. In addition, the epithelial-mesenchymal transition (EMT) markers and PI3K/AKT mechanism markers were examined via western blot analysis to elucidate the underlying mechanisms. In HCC, GCNT3 was significantly overexpressed, which was connected with enhanced tumor aggressiveness and an unfavorable prognosis of individuals. In vitro experiments demonstrated that elevated levels of GCNT3 promoted cell growth, migration, cell cycle development, and invasion, in addition to EMT, while suppressing apoptosis. Conversely, knockdown of GCNT3 exerted the opposite effects. GCNT3 overexpression increased PI3K/AKT phosphorylation in HCC cells, and LY294002 counteracted the impacts of upregulated GCNT3 on cell cycle, migration, invasion, proliferation, and EMT in HCC. The investigation showed that GCNT3 may enhance HCC progression and EMT by stimulating PI3K/AKT mechanism.