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.

MEHunter: transformer-based mobile element variant detection from long reads.

SUMMARY: Mobile genetic elements (MEs) are heritable mutagens that significantly contribute to genetic diseases. The advent of long-read sequencing technologies, capable of resolving large DNA fragments, offers promising prospects for the comprehensive detection of ME variants (MEVs). However, achieving high precision while maintaining recall performance remains challenging mainly brought by the variable length and similar content of MEV signatures, which are often obscured by the noise in long reads. Here, we propose MEHunter, a high-performance MEV detection approach utilizing a fine-tuned transformer model adept at identifying potential MEVs with fragmented features. Benchmark experiments on both simulated and real datasets demonstrate that MEHunter consistently achieves higher accuracy and sensitivity than the state-of-the-art tools. Furthermore, it is capable of detecting novel potentially individual-specific MEVs that have been overlooked in published population projects. AVAILABILITY AND IMPLEMENTATION: MEHunter is available from https://github.com/120L021101/MEHunter.

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.

Quantitatively Evaluating Interactions between Patient-Derived Organoids and Autologous Immune Cells by Microfluidic Chip.

The coculture of patient-derived tumor organoids (PDOs) and autologous immune cells has been considered as a useful ex vivo surrogate of in vivo tumor-immune environment. However, the immune interactions between PDOs and autologous immune cells, including immune-mediated killing behaviors and immune-related cytokine variations, have yet to be quantitatively evaluated. This study presents a microfluidic chip for quantifying interactions between PDOs and autologous immune cells (IOI-Chip). A baffle-well structure is designed to ensure efficient trapping, long-term coculturing, and in situ fluorescent observation of a limited amount of precious PDOS and autologous immune cells, while a microbeads-based immunofluorescence assay is designed to simultaneously quantify multiple kinds of immune-related cytokines in situ. The PDO apoptosis and 2 main immune-related cytokines, TNF-α and IFN-γ, are simultaneously quantified using samples from a lung cancer patient. This study provides, for the first time, a capability to quantify interactions between PDOs and autologous immune cells at 2 levels, the immune-mediated killing behavior, and multiple immune-related cytokines, laying the technical foundation of ex vivo assessment of patient immune response.

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.

Hollow Structure Derived Phosphide Nanosheets for Water Oxidation.

Avoiding the stacking of active sites in catalyst structural design is a promising route for realizing active oxygen evolution reaction (OER). Herein, using a CoFe Prussian blue analoge cube with hollow structure (C-CoFe PBA) as a derived support, a highly effective Ni2P-FeP4-Co2P catalyst with a larger specific surface area is reported. Benefiting from the abundant active sites and fast charge transfer capability of the phosphide nanosheets, the Ni2P-FeP4-Co2P catalyst in 1 m KOH requires only overpotentials of 248 and 277 mV to reach current density of 10 and 50 mA cm-2 and outperforms the commercial catalyst RuO2 and most reported non-noble metal OER catalysts. In addition, the two-electrode system consisting of Ni2P-FeP4-Co2P and Pt/C is able to achieve a current density of 10 and 50 mA cm-2 at 1.529 and 1.65 V. This work provides more ideas and directions for synthesizing transition metal catalysts for efficient OER performance.

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.