Pan-Cancer Identification of Prognostic-Associated Metabolic Pathways.

Metabolic dysregulation has been reported involving in the clinical outcomes of multiple cancers. However, systematical identification of the impact of metabolic pathways on cancer prognosis is still lacking. Here, we performed a pan-cancer analysis of popular metabolic checkpoint genes and pathways with cancer prognosis by integrating information of clinical survival with gene expression and pathway activity in multiple cancer patients. By discarding the effects of age and sex, we revealed extensive and significant associations between the survival of cancer patients and the expression of metabolic checkpoint genes, as well as the activities of three primary metabolic pathways: amino acid metabolism, carbohydrate metabolism, lipid metabolism, and eight nonprimary metabolic pathways. Among multiple cancers, we found the survival of kidney renal clear cell carcinoma and low-grade glioma exhibit high metabolic dependence. Our work systematically assesses the impact of metabolic checkpoint genes and pathways on cancer prognosis, providing clues for further study of cancer diagnosis and therapy.

Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

Numerical and experimental investigations on the propagation property of a mid-infrared 7 × 1 multimode fiber combiner.

Mid-infrared fiber combiners have great potential in power and spectral combination. However, studies on mid-infrared transmission optical field distributions using these combiners are limited. In this study, we designed and fabricated a 7 × 1 multimode fiber combiner based on sulfur-based glass fibers and observed approximately 80% per-port transmission efficiency at 4.778 µm wavelength. We investigated the propagation properties of the prepared combiners and explored the effects of transmission wavelength, output fiber length, and fusion deviation on the transmitted optical field and beam quality factor M2. Additionally, we assessed the effect of coupling on the excitation mode and spectral combination of the mid-infrared fiber combiner for multiple light sources. Our results provide an in-depth understanding of the propagation properties of the mid-infrared multimode fiber combiners, which may find applications in high-beam-quality laser devices.

Prototypical multiple instance learning for predicting lymph node metastasis of breast cancer from whole-slide pathological images.

Computerized identification of lymph node metastasis of breast cancer (BCLNM) from whole-slide pathological images (WSIs) can largely benefit therapy decision and prognosis analysis. Besides the general challenges of computational pathology, like extra-high resolution, very expensive fine-grained annotation, etc., two particular difficulties with this task lie in (1) modeling the significant inter-tumoral heterogeneity in BCLNM pathological images, and (2) identifying micro-metastases, i.e., metastasized tumors with tiny foci. Towards this end, this paper presents a novel weakly supervised method, termed as Prototypical Multiple Instance Learning (PMIL), to learn to predict BCLNM from WSIs with slide-level class labels only. PMIL introduces the well-established vocabulary-based multiple instance learning (MIL) paradigm into computational pathology, which is characterized by utilizing the so-called prototypes to model pathological data and construct WSI features. PMIL mainly consists of two innovatively designed modules, i.e., the prototype discovery module which acquires prototypes from training data by unsupervised clustering, and the prototype-based slide embedding module which builds WSI features by matching constitutive patches against the prototypes. Relative to existing MIL methods for WSI classification, PMIL has two substantial merits: (1) being more explicit and interpretable in modeling the inter-tumoral heterogeneity in BCLNM pathological images, and (2) being more effective in identifying micro-metastases. Evaluation is conducted on two datasets, i.e., the public Camelyon16 dataset and the Zbraln dataset created by ourselves. PMIL achieves an AUC of 88.2% on Camelyon16 and 98.4% on Zbraln (at 40x magnification factor), which consistently outperforms other compared methods. Comprehensive analysis will also be carried out to further reveal the effectiveness and merits of the proposed method.

circRIP: an accurate tool for identifying circRNA-RBP interactions.

Circular ribonucleic acids (RNAs) (circRNAs) are formed by covalently linking the downstream splice donor and the upstream splice acceptor. One of the most important functions of circRNAs is mainly exerted through binding RNA-binding proteins (RBPs). However, there is no efficient algorithm for identifying genome-wide circRNA-RBP interactions. Here, we developed a unique algorithm, circRIP, for identifying circRNA-RBP interactions from RNA immunoprecipitation sequencing (RIP-Seq) data. A simulation test demonstrated the sensitivity and specificity of circRIP. By applying circRIP, we identified 95 IGF2BP3-binding circRNAs based on the IGF2BP3 RIP-Seq dataset. We further identified 2823 and 1333 circRNAs binding to >100 RBPs in K562 and HepG2 cell lines, respectively, based on enhanced cross-linking immunoprecipitation (eCLIP) data, demonstrating the significance to survey the potential interactions between circRNAs and RBPs. In this study, we provide an accurate and sensitive tool, circRIP (https://github.com/bioinfolabwhu/circRIP), to systematically identify RBP and circRNA interactions from RIP-Seq and eCLIP data, which can significantly benefit the research community for the functional exploration of circRNAs.

CSCD2: an integrated interactional database of cancer-specific circular RNAs.

The significant function of circRNAs in cancer was recognized in recent work, so a well-organized resource is required for characterizing the interactions between circRNAs and other functional molecules (such as microRNA and RNA-binding protein) in cancer. We previously developed cancer-specific circRNA database (CSCD), a comprehensive database for cancer-specific circRNAs, which is widely used in circRNA research. Here, we updated CSCD to CSCD2 (http://geneyun.net/CSCD2 or http://gb.whu.edu.cn/CSCD2), which includes significantly more cancer-specific circRNAs identified from a large number of human cancer and normal tissues/cell lines. CSCD2 contains >1000 samples (825 tissues and 288 cell lines) and identifies a large number of circRNAs: 1 013 461 cancer-specific circRNAs, 1 533 704 circRNAs from only normal samples and 354 422 circRNAs from both cancer and normal samples. In addition, CSCD2 predicts potential miRNA-circRNA and RBP-circRNA interactions using binding motifs from >200 RBPs and 2000 microRNAs. Furthermore, the potential full-length and open reading frame sequence of these circRNAs were also predicted. Collectively, CSCD2 provides a significantly enhanced resource for exploring the function and regulation of circRNAs in cancer.

YTHDF1 promotes mRNA degradation via YTHDF1-AGO2 interaction and phase separation.

OBJECTIVES: YTHDF1 is known as a m6 A reader protein, and many researches of YTHDF1 focused on the regulation of mRNA translation efficiency. However, YTHDF1 is also related to RNA degradation, but how YTHDF1 regulates mRNA degradation is indefinite. Liquid-liquid phase separation (LLPS) underlies the formation of membraneless compartments in mammal cells, and there are few reports focused on the correlation of RNA degradation with LLPS. In this research, we focused on the mechanism of YTHDF1 degraded mRNA through LLPS. MATERIALS AND METHODS: The CRISPR/Cas9 knock out system was used to establish the YTHDF1 knock out (YTHDF1-KO) cell lines (HEK293 and HeLa) and METTL14 knock out (METTL14-KO) cell line (HEK293). 4SU-TT-seq was used to check the half-life changes of mRNAs. Actinomycin D and qPCR were used to test the half-life changes of individual mRNA. RNA was stained with SYTO RNA-select dye in wild type (WT) and YTHDF1-KO HeLa cell lines. Co-localization of YTHDF1 and AGO2 was identified by immunofluorescence. The interaction domain of YTHDF1 and AGO2 was identified by western blot. Phase separation of YTHDF1 was performed in vitro and in vivo. Fluorescence recovery after photobleaching (FRAP) was performed on droplets as an assessment of their liquidity. RESULTS: In this research, we found that deletion of YTHDF1 led to massive RNA patches deposited in cytoplasm. The results of 4SU-TT-seq showed that deletion of YTHDF1 would prolong the half-life of mRNAs. Immunofluorescence data showed that YTHDF1 and AGO2 could co-localize in P-body, and Co-IP results showed that YTHDF1 could interact with AGO2 through YT521-B homology (YTH) domain. We confirmed that YTHDF1 could undergo phase separation in vitro and in vivo, and compared with AGO2, YTHDF1 was more important in P-body formation. The FRAP results showed that liquid AGO2 droplets would convert to gel/solid when YTHDF1 was deleted. As AGO2 plays important roles in miRISCs, we also found that miRNA-mediate mRNA degradation is related to YTHDF1. CONCLUSIONS: YTHDF1 recruits AGO2 through the YTH domain. YTHDF1 degrades targeting mRNAs by promoting P-body formation through LLPS. The deletion of YTHDF1 causes the P-body to change from liquid droplets to gel/solid droplets, and form AGO2/RNA patches, resulting in a degradation delay of mRNAs. These findings reveal a previously unrecognized crosstalk between YTHDF1 and AGO2, raising a new sight of mRNA post-transcriptional regulation by YTHDF1.

ADEIP: an integrated platform of age-dependent expression and immune profiles across human tissues.

Gene expression and immune status in human tissues are changed with aging. There is a need to develop a comprehensive platform to explore the dynamics of age-related gene expression and immune profiles across tissues in genome-wide studies. Here, we collected RNA-Seq datasets from GTEx project, containing 16 704 samples from 30 major tissues in six age groups ranging from 20 to 79 years old. Dynamic gene expression along with aging were depicted and gene set enrichment analysis was performed among those age groups. Genes from 34 known immune function categories and immune cell compositions were investigated and compared among different age groups. Finally, we integrated all the results and developed a platform named ADEIP (http://gb.whu.edu.cn/ADEIP or http://geneyun.net/ADEIP), integrating the age-dependent gene expression and immune profiles across tissues. To demonstrate the usage of ADEIP, we applied two datasets: severe acute respiratory syndrome coronavirus 2 and human mesenchymal stem cells-assoicated genes. We also included the expression and immune dynamics of these genes in the platform. Collectively, ADEIP is a powerful platform for studying age-related immune regulation in organogenesis and other infectious or genetic diseases.

Silencing of circIgf1r plays a protective role in neuronal injury via regulating astrocyte polarization during epilepsy.

Epilepsy is a common brain disorder, repeated seizures of epilepsy may lead to a series of brain pathological changes such as neuronal or glial damage. However, whether circular RNAs are involved in neuronal injury during epilepsy is not fully understood. Here, we screened circIgf1r in the status epilepticus model through circRNA sequencing, and found that it was upregulated after the status epilepticus model through QPCR analysis. Astrocytes polarizing toward neurotoxic A1 phenotype and neurons loss were observed after status epilepticus. Through injecting circIgf1r siRNA into the lateral ventricle, it was found that knocking down circIgf1r in vivo would induce the polarization of astrocytes to phenotype A2 and reduce neuronal loss. The results in vitro further confirmed that inhibiting the expression of circIgf1r in astrocytes could protect neurons by converting reactive astrocytes from A1 to the protective A2. In addition, knocking down circIgf1r in astrocytes could functionally promote astrocyte autophagy and relieve the destruction of 4-AP-induced autophagy flux. In terms of mechanism, circIgf1r promoted the polarization of astrocytes to phenotype A1 by inhibiting autophagy. Taken together, our results reveal circIgf1r may serve as a potential target for the prevention and treatment of neuron damage after epilepsy.

Circular RNA circYPEL2: A Novel Biomarker in Cervical Cancer.

Cervical cancer (CC) is one of the most threatening diseases in women. Circular RNAs (circRNAs) have been reported to be cancer hallmarks, but typical circRNAs in CC were rarely indicated. Through high-throughput sequencing in CC and normal cervix tissues, circYPEL2 (hsa_circ_0005600) was proposed as a candidate circRNA. CircYPEL2 exhibited significantly high expression in CC tissue and strong stability in CC cell lines. Furthermore, knockdown and overexpression of circYPEL2 indicated the potential involvement in CC proliferation, migration and invasion. Finally, the downstream regulatory genes of circYPEL2 were investigated by knockdown experiment in CC cell lines with high-throughput sequencing. In summary, our work identified circYPEL2 as a potential biomarker for clinical research of cervical cancer.

Breath therapy for patients with chronic nonspecific low back pain: A protocol of systematic review and meta-analysis.

BACKGROUND: Chronic nonspecific low back pain (CNLBP) has become a major global public health problem. Its high incidence rate and high disability rate are so damaging both to individuals and communities. At present, many countries' clinical guidelines recommend exercise therapy. Breath therapy is one of the exercise therapies, playing an important role in exercise therapy. Some studies have shown that breath therapy has a considerable therapeutic effect on low back pain, but there is no specific conclusion. The aim of our study is to answer the question: if breath therapy is effective and safe for CNLBP? METHODS: The following databases will be searched: English databases (including Web of Science, the Cochrane Library (Central), EMBASE, MEDLINE, Allied and Alternative Medicine) and Chinese databases (including Chinese National Knowledge Infrastructure, Wanfang data and Chinese Scientific Journals Database [VIP]). The literature search will be constructed around search terms for breath therapy, search terms for chronic nonspecific low back pain and search terms for randomized controlled trials. The primary outcomes were related to duration, intensity, attack frequency of pain, and the secondary outcomes were related to physical function, quality of life, and adverse events related to interventions. Endnote software 9.1 will be applied in selecting study, Review Manager software 5.3 will be applied in analyzing and synthesizing. RESULTS: The results will provide evidence to judge whether breath therapy is effective and safe for CNLBP. CONCLUSION: Our research will provide reliable evidence of breath therapy for CNLBP. REGISTRATION: International Prospective Register of Systematic Reviews (PROSPERO) CRD42020156340.

CircRNA inhibits DNA damage repair by interacting with host gene.

BACKGROUND: Deregulated circular RNAs (circRNAs) are associated with the development of cancer and therapy resistance. However, functional research of circRNAs mostly focus on potential miRNA or protein binding and more potential regulation of circRNA on host gene DNA in cancers are yet to be inspected. METHOD: We performed total RNA sequencing on clinical breast cancer samples and identified the expression patterns of circRNAs and corresponding host genes in patient blood, tumor and adjacent normal tissues. qPCR, northern blot and in situ hybridization were used to validate the dysregulation of circRNA circSMARCA5. A series of procedures including R-loop dot-blotting, DNA-RNA immunoprecipitation and mass spectrum, etc. were conducted to explore the regulation of circSMARCA5 on the transcription of exon 15 of SMARCA5. Moreover, immunofluorescence and in vivo experiments were executed to investigate the overexpression of circSMARCA5 with drug sensitivities. RESULTS: We found that circRNAs has average higher expression over its host linear genes in peripheral blood. Compared to adjacent normal tissues, circSMARCA5 is decreased in breast cancer tissues, contrary to host gene SMARCA5. The enforced expression of circSMARCA5 induced drug sensitivity of breast cancer cell lines in vitro and in vivo. Furthermore, we demonstrated that circSMARCA5 can bind to its parent gene locus, forming an R-loop, which results in transcriptional pausing at exon 15 of SMARCA5. CircSMARCA5 expression resulted in the downregulation of SMARCA5 and the production of a truncated nonfunctional protein, and the overexpression of circSMARCA5 was sufficient to improve sensitivity to cytotoxic drugs. CONCLUSION: Our results revealed a new regulatory mechanism for circRNA on its host gene and provided evidence that circSMARCA5 may serve as a therapeutic target for drug-resistant breast cancer patients.

GR/HDAC2/TGFßR1 pathway contributes to prenatal caffeine induced-osteoarthritis susceptibility in male adult offspring rats.

Prenatal caffeine exposure (PCE) induces developmental toxicity of multi-organ and susceptibility to multi-disease in offspring. However, the effects of PCE on osteoarthritis susceptibility in adult offspring and its intrauterine programming mechanism remain to be further investigated. Here, we found that PCE induced susceptibility to osteoarthritis in male adult offspring rats, which was related to the inhibited function of cartilage matrix synthesis from fetuses to adults. Meanwhile, PCE consistently downregulated the H3K9ac and expression levels of transforming growth factor ß receptor 1 (TGFßR1), and then blocked TGFß signaling pathway, which contributed to the suppressed cartilage matrix synthesis. Moreover, the high level of corticosterone caused by PCE reduced the H3K9ac level on TGFßR1 promoter region through acting on glucocorticoids receptor (GR) and recruiting histone deacetylase 2 (HDAC2) into the nucleus of fetal chondrocytes. Taken together, PCE induced osteoarthritis susceptibility in male adult offspring rats, which was attributed to the low-functional programming of TGFßR1 induced by corticosterone via GR/HDAC2 signaling.

Silencing lncRNA SNHG6 suppresses proliferation and invasion of breast cancer cells through miR-26a/VASP axis.

The important role of LncRNA in the development of breast cancer is attracting more and more attention. In the previous study, we found that the expression level of LncRNA SNHG6 in breast cancer tissues and cells was significantly increased, but its mechanism in the development of breast cancer was still unclear. Our study found that knockdown of SNHG6 significantly inhibited the proliferation, migration and invasion of breast cancer cells MCF-7 and MDA-MB-231 cells. Further study showed that knockdown of SNHG6 significantly inhibited the expression level of VASP. More importantly, SNHG6 and VASP both can bind directly to miR-26a, suggesting that SNHG6 could act as a ceRNA to sponge miR-26a, thereby promoting the expression of VASP, which leading to activated proliferation, migration and invasion of breast cancer cells. Taken together, this study revealed the important role of the SNHG6/miR-26a/VASP regulatory network in the development of breast cancer, and provided a reference for exploring new pathogenesis and biomarkers of breast cancer.

The RNA N6-methyladenosine modification landscape of human fetal tissues.

A single genome gives rise to diverse tissues through complex epigenomic mechanisms, including N6-methyladenosine (m6A), a widespread RNA modification that is implicated in many biological processes. Here, to explore the global landscape of m6A in human tissues, we generated 21 whole-transcriptome m6A methylomes across major fetal tissues using m6A sequencing. These data reveal dynamic m6A methylation, identify large numbers of tissue differential m6A modifications and indicate that m6A is positively correlated with gene expression homeostasis. We also report m6A methylomes of long intergenic non-coding RNA (lincRNA), finding that enhancer lincRNAs are enriched for m6A. Tissue m6A regions are often enriched for single nucleotide polymorphisms that are associated with the expression of quantitative traits and complex traits including common diseases, which may potentially affect m6A modifications. Finally, we find that m6A modifications preferentially occupy genes with CpG-rich promoters, features of which regulate RNA transcript m6A. Our data indicate that m6A is widely regulated by human genetic variation and promoters, suggesting a broad involvement of m6A in human development and disease.

Genome-wide identification of cancer-specific alternative splicing in circRNA.

Circular RNA (circRNA) is a group of RNA families generated by RNA circularization, which was discovered ubiquitously across different cancers. However, the internal structure of circRNA is difficult to determine due to alternative splicing that occurs in its exons and introns. Furthermore, cancer-specific alternative splicing of circRNA is less likely to be identified. Here, we proposed a de novo algorithm, CircSplice, that could identify internal alternative splicing in circRNA and compare differential circRNA splicing events between different conditions ( http://gb.whu.edu.cn/CircSplice or https://github.com/GeneFeng/CircSplice ). By applying CircSplice in clear cell renal cell carcinoma and bladder cancer, we detected 4498 and 2977 circRNA alternative splicing (circ-AS) events in the two datasets respectively and confirmed the expression of circ-AS events by RT-PCR. We further inspected the distributions and patterns of circ-AS in cancer and adjacent normal tissues. To further understand the potential functions of cancer-specific circ-AS, we classified those events into tumor suppressors and oncogenes and performed pathway enrichment analysis. This study is the first comprehensive view of cancer-specific circRNA alternative splicing, which could contribute significantly to regulation and functional research of circRNAs in cancers.

The prognostic value of matrix metalloproteinase-7 and matrix metalloproteinase-15 in acute myeloid leukemia.

Matrix metalloproteinases (MMPs), a family of zinc-dependent endopeptidases, are involved in a variety of physiological and pathological processes. We analyzed 11 data sets from Gene Expression Omnibus Database and found that MMP7 and MMP15 were highly expressed in multiple carcinomas. GSE13204 showed that MMP7 and MMP15 were overexpressed in acute myeloid leukemia (AML) patients. The Cancer Genome Atlas data set exhibited that high expression of MMP7 or MMP15 in bone marrow (BM) of AML patients predicted poor overall survival. The χ 2 test results indicated that high expression level of MMP7 and MMP15 were correlated with high-risk stratification and high BM blast cell percentage in AML patients. To confirm these findings, we performed reverse-transcription quantitative polymerase chain reaction (RT-qPCR) and found that MMP7 and MMP15 were highly expressed in three AML cell lines. Further study showed that MMP7 and MMP15 were highly expressed both in BM and peripheral blood in collected AML samples compared with healthy individuals. Additionally, long noncoding RNA (lncRNA) microarray of BM samples of AML patients revealed that multiple lncRNAs were correlated with MMP7 and MMP15, suggesting that lncRNAs might be involved in the pathogenesis of AML via modulating MMPs. In conclusion, our study uncovers the potential roles of MMP7 and MMP15 in the prognosis of AML.

CVm6A: A Visualization and Exploration Database for m6As in Cell Lines.

N6-methyladenosine (m6A) has been identified in various biological processes and plays important regulatory functions in diverse cells. However, there is still no visualization database for exploring global m6A patterns across cell lines. Here we collected all available MeRIP-Seq and m6A-CLIP-Seq datasets from public databases and identified 340,950 and 179,201 m6A peaks dependent on 23 human and eight mouse cell lines respectively. Those m6A peaks were further classified into mRNA and lncRNA groups. To better understand the potential function of m6A, we then mapped m6A peaks in different subcellular components and gene regions. Among those human m6A modification, 190,050 and 150,900 peaks were identified in cancer and non-cancer cells, respectively. Finally, all results were integrated and imported into a visualized cell-dependent m6A database CVm6A. We believe the specificity of CVm6A could significantly contribute to the research for the function and regulation of cell-dependent m6A modification in disease and development.

Loss of miR-146b-5p promotes T cell acute lymphoblastic leukemia migration and invasion via the IL-17A pathway.

Metastatic disease remains the primary cause of death for individuals with T cell acute lymphoblastic leukemia (T-ALL). microRNAs (miRNAs) play important roles in the pathogenesis of T-ALL by inhibiting gene expression at posttranscriptional levels. The goal of the current project is to identify any significant miRNAs in T-ALL metastasis. We observed miR-146b-5p to be downregulated in T-ALL patients and cell lines, and bioinformatics analysis implicated miR-146b-5p in the hematopoietic system. miR-146b-5p inhibited the migration and invasion in T-ALL cells. Interleukin-17A (IL-17A) was predicted to be a target of miR-146b-5p; this was confirmed by luciferase assays. Interestingly, T-ALL patients and cell lines secreted IL-17A and expressed the IL-17A receptor (IL-17RA). IL-17A/IL-17RA interactions promoted strong T-ALL cell migration and invasion responses. Gene set enrichment analysis (GSEA) and quantitative polymerase chain reaction (qPCR) analysis indicated that matrix metallopeptidase-9 (MMP9), was a potential downstream effector of IL-17A activation, and nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling was also implicated in this process. Moreover, IL-17A activation promoted T-ALL cell metastasis to the liver in IL17A -/- mouse models. These results indicate that reduced miR-146b-5p expression in T-ALL may lead to the upregulation of IL-17A, which then promotes T-ALL cell migration and invasion by upregulating MMP9 via NF-κB signaling.