Functional characteristics of DNA N6-methyladenine modification based on long-read sequencing in pancreatic cancer.

Abnormalities of DNA modifications are closely related to the pathogenesis and prognosis of pancreatic cancer. The development of third-generation sequencing technology has brought opportunities for the study of new epigenetic modification in cancer. Here, we screened the N6-methyladenine (6mA) and 5-methylcytosine (5mC) modification in pancreatic cancer based on Oxford Nanopore Technologies sequencing. The 6mA levels were lower compared with 5mC and upregulated in pancreatic cancer. We developed a novel method to define differentially methylated deficient region (DMDR), which overlapped 1319 protein-coding genes in pancreatic cancer. Genes screened by DMDRs were more significantly enriched in the cancer genes compared with the traditional differential methylation method (P < 0.001 versus P = 0.21, hypergeometric test). We then identified a survival-related signature based on DMDRs (DMDRSig) that stratified patients into high- and low-risk groups. Functional enrichment analysis indicated that 891 genes were closely related to alternative splicing. Multi-omics data from the cancer genome atlas showed that these genes were frequently altered in cancer samples. Survival analysis indicated that seven genes with high expression (ADAM9, ADAM10, EPS8, FAM83A, FAM111B, LAMA3 and TES) were significantly associated with poor prognosis. In addition, the distinction for pancreatic cancer subtypes was determined using 46 subtype-specific genes and unsupervised clustering. Overall, our study is the first to explore the molecular characteristics of 6mA modifications in pancreatic cancer, indicating that 6mA has the potential to be a target for future clinical treatment.

DiseaseMeth version 3.0: a major expansion and update of the human disease methylation database.

DNA methylation has a growing potential for use as a biomarker because of its involvement in disease. DNA methylation data have also substantially grown in volume during the past 5 years. To facilitate access to these fragmented data, we proposed DiseaseMeth version 3.0 based on DiseaseMeth version 2.0, in which the number of diseases including increased from 88 to 162 and High-throughput profiles samples increased from 32 701 to 49 949. Experimentally confirmed associations added 448 pairs obtained by manual literature mining from 1472 papers in PubMed. The search, analyze and tools sections were updated to increase performance. In particular, the FunctionSearch now provides for the functional enrichment of genes from localized GO and KEGG annotation. We have also developed a unified analysis pipeline for identifying differentially DNA methylated genes (DMGs) from the original data stored in the database. 22 718 DMGs were found in 99 diseases. These DMGs offer application in disease evaluation using two self-developed online tools, Methylation Disease Correlation and Cancer Prognosis & Co-Methylation. All query results can be downloaded and can also be displayed through a box plot, heatmap or network module according to whichever search section is used. DiseaseMeth version 3.0 is freely available at http://diseasemeth.edbc.org/.

Accurate Prediction of Prognosis by Integrating Clinical and Molecular Characteristics in Colon Cancer.

Various factors affect the prognosis of patients with colon cancer. Complicated factors are found to be conducive to accurate assessment of prognosis. In this study, we developed a series of prognostic prediction models for survival time of colon cancer patients after surgery. Analysis of nine clinical characteristics showed that the most important factor was the positive lymph node ratio (LNR). High LNR was the most important clinical factor affecting 1- and 3-year survival; M0&age < 70 was the most important feature for 5 years. The performance of the model was improved through the integration of clinical characteristics and four types of molecule features (mRNA, lncRNA, miRNA, DNA methylation). The model provides guidance for clinical practice. According to the high-risk molecular features combined with age ≥ 70&T3, poorly differentiated or undifferentiated, M0&well differentiated, M0&T2, LNR high, T4&poorly differentiated, or undifferentiated, the survival time may be less than 1 year; for patients with high risk of molecular features combined with M0&T2, M0&T4, LNR 0& M0, LNR median &T3, and LNR high, the survival is predicted less than 3 years; and the survival of patients with M1&T3, M0 and high risk molecular features is less than 5 years. Using multidimensional and complex patient information, this study establishes potential criteria for clinicians to evaluate the survival of patients for colon cancer.

M6ADD: a comprehensive database of m6A modifications in diseases.

N6-methyladenosine (m6A) modification is an important regulatory factor affecting diseases, including multiple cancers and it is a developing direction for targeted disease therapy. Here, we present the M6ADD (m6A-diseases database) database, a public data resource containing manually curated data on potential m6A-disease associations for which some experimental evidence is available; the related high-throughput sequencing data are also provided and analysed by using different computational methods. To give researchers a tool to query the m6A modification data, the M6ADD was designed as a web-based comprehensive resource focusing on the collection, storage and online analysis of m6A modifications, aimed at exploring the associations between m6A modification and gene disorders and diseases. The M6ADD includes 222 experimentally confirmed m6A-disease associations, involving 59 diseases from a review of more than 2000 published papers. The M6ADD also includes 409,229 m6A-disease associations obtained by computational and statistical methods from 30 high-throughput sequencing datasets. In addition, we provide data on 5239 potential m6A regulatory proteins related to 24 cancers based on network analysis prediction methods. In addition, we have developed a tool to explore the function of m6A-modified genes through the protein-protein interaction networks. The M6ADD can be accessed at http://m6add.edbc.org/.

Prediction of transcription factors binding events based on epigenetic modifications in different human cells.

Aim: We aim to predict transcription factor (TF) binding events from knowledge of gene expression and epigenetic modifications. Materials & methods: TF-binding events based on the Encode project and The Cancer Genome Atlas data were analyzed by the random forest method. Results: We showed the high performance of TF-binding predictive models in GM12878, HeLa, HepG2 and K562 cell lines and applied them to other cell lines and tissues. The genes bound by the top TFs (MAX and MAZ) were significantly associated with cancer-related processes such as cell proliferation and DNA repair. Conclusion: We successfully constructed TF-binding predictive models in cell lines and applied them in tissues.

CLING: Candidate Cancer-Related lncRNA Prioritization via Integrating Multiple Biological Networks.

Identification and characterization of lncRNAs in cancer with a view to their application in improving diagnosis and therapy remains a major challenge that requires new and innovative approaches. We have developed an integrative framework termed "CLING", aimed to prioritize candidate cancer-related lncRNAs based on their associations with known cancer lncRNAs. CLING focuses on joint optimization and prioritization of all candidates for each cancer type by integrating lncRNA topological properties and multiple lncRNA-centric networks. Validation analyses revealed that CLING is more effective than prioritization based on a single lncRNA network. Reliable AUC (Area Under Curve) scores were obtained across 10 cancer types, ranging from 0.85 to 0.94. Several novel lncRNAs predicted in the top 10 candidates for various cancer types have been confirmed by recent biological experiments. Furthermore, using a case study on liver hepatocellular carcinoma as an example, CLING facilitated the successful identification of novel cancer lncRNAs overlooked by differential expression analyses (DEA). This time- and cost-effective computational model may provide a valuable complement to experimental studies and assist in future investigations on lncRNA involvement in the pathogenesis of cancers. We have developed a web-based server for users to rapidly implement CLING and visualize data, which is freely accessible at http://bio-bigdata.hrbmu.edu.cn/cling/. CLING has been successfully applied to predict a few potential lncRNAs from thousands of candidates for many cancer types.

SEA version 3.0: a comprehensive extension and update of the Super-Enhancer archive.

Super-enhancers (SEs) are critical for the transcriptional regulation of gene expression. We developed the super-enhancer archive version 3.0 (SEA v. 3.0, http://sea.edbc.org) to extend SE research. SEA v. 3.0 provides the most comprehensive archive to date, consisting of 164 545 super-enhancers. Of these, 80 549 are newly identified from 266 cell types/tissues/diseases using an optimized computational strategy, and 52 have been experimentally confirmed with manually curated references. We now support super-enhancers in 11 species including 7 new species (zebrafish, chicken, chimp, rhesus, sheep, Xenopus tropicalis and stickleback). To facilitate super-enhancer functional analysis, we added several new regulatory datasets including 3 361 785 typical enhancers, chromatin interactions, SNPs, transcription factor binding sites and SpCas9 target sites. We also updated or developed new criteria query, genome visualization and analysis tools for the archive. This includes a tool based on Shannon Entropy to evaluate SE cell type specificity, a new genome browser that enables the visualization of SE spatial interactions based on Hi-C data, and an enhanced enrichment analysis interface that provides online enrichment analyses of SE related genes. SEA v. 3.0 provides a comprehensive database of all available SE information across multiple species, and will facilitate super-enhancer research, especially as related to development and disease.

LncACTdb 2.0: an updated database of experimentally supported ceRNA interactions curated from low- and high-throughput experiments.

We describe LncACTdb 2.0 (http://www.bio-bigdata.net/LncACTdb/), an updated and significantly expanded database which provides comprehensive information of competing endogenous RNAs (ceRNAs) in different species and diseases. We have updated LncACTdb 2.0 with more data and several new features, including (i) manually curating 2663 experimentally supported ceRNA interactions from >5000 published literatures; (ii) expanding the scope of the database up to 23 species and 213 diseases/phenotypes; (iii) curating more ceRNA types such as circular RNAs and pseudogenes; (iv) identifying and scoring candidate lncRNA-associated ceRNA interactions across 33 cancer types from TCGA data; (v) providing illustration of survival, network and cancer hallmark information for ceRNAs. Furthermore, several flexible online tools including LncACT-Get, LncACT-Function, LncACT-Survival, LncACT-Network and LncACTBrowser have been developed to perform customized analysis, functional analysis, survival analysis, network illustration and genomic visualization. LncACTdb 2.0 also provides newly designed, user-friendly web interfaces to search, browse and download all the data. The BLAST interface is convenient for users to query dataset by inputting custom sequences. The Hot points interface provides users the most studied items by others. LncACTdb 2.0 is a continually updated database and will serve as an important resource to explore ceRNAs in physiological and pathological processes.

Lnc2Cancer v2.0: updated database of experimentally supported long non-coding RNAs in human cancers.

Lnc2Cancer 2.0 (http://www.bio-bigdata.net/lnc2cancer) is an updated database that provides comprehensive experimentally supported associations between lncRNAs and human cancers. In Lnc2Cancer 2.0, we have updated the database with more data and several new features, including (i) exceeding a 4-fold increase over the previous version, recruiting 4989 lncRNA-cancer associations between 1614 lncRNAs and 165 cancer subtypes. (ii) newly adding about 800 experimentally supported circulating, drug-resistant and prognostic-related lncRNAs in various cancers. (iii) appending the regulatory mechanism of lncRNA in cancer, including microRNA (miRNA), transcription factor (TF), variant and methylation regulation. (iv) increasing more than 70 high-throughput experiments (microarray and next-generation sequencing) of lncRNAs in cancers. (v) Scoring the associations between lncRNA and cancer to evaluate the correlations. (vi) updating the annotation information of lncRNAs (version 28) and containing more detailed descriptions for lncRNAs and cancers. Moreover, a newly designed, user-friendly interface was also developed to provide a convenient platform for users. In particular, the functions of browsing data by cancer primary organ, biomarker type and regulatory mechanism, advanced search following several features and filtering the data by LncRNA-Cancer score were enhanced. Lnc2Cancer 2.0 will be a useful resource platform for further understanding the associations between lncRNA and human cancer.

Inferences of individual drug responses across diverse cancer types using a novel competing endogenous RNA network.

Differences in individual drug responses are an obstacle to progression in cancer treatment, and predicting responses would help to plan treatment. The accumulation of cancer molecular profiling and drug response data provides opportunities and challenges to identify novel molecular signatures and mechanisms of tumor responsiveness to drugs. This study evaluated drug responses with a competing endogenous RNA (ceRNA) system that depended on competition between diverse RNA species. We identified drug response-related ceRNA (DRCEs) by combining the sequence and expression data of long noncoding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA), and the survival data of cancer patients treated with drugs. We constructed a patient-drug two-layer integrated network and used a linear weighting method to predict individual drug responses. DRCEs were found to be significantly enriched in known cancer and drug-associated data resources, involved in biological processes known to mediate drug responses, and correlated to drug activity in cancer cell lines. The dysregulation of DRCE expression influenced drug response-associated functions and pathways, suggesting DRCEs as potential therapeutic targets affecting drug responses. A further case study in breast invasive carcinoma (BRCA) found that DRCE expression was consistent with the drug response pattern and the aberrant expression of the two NEAT1-related DRCEs may lead to poor response to tamoxifen therapy for patients with TP53 mutations. In summary, this study provides a framework for ceRNA-based evaluation of clinical drug responses across multiple cancer types. Understanding the underlying molecular mechanisms of drug responses will allow improved response to chemotherapy and outcomes of cancer treatment.

Lnc2Meth: a manually curated database of regulatory relationships between long non-coding RNAs and DNA methylation associated with human disease.

Lnc2Meth (http://www.bio-bigdata.com/Lnc2Meth/), an interactive resource to identify regulatory relationships between human long non-coding RNAs (lncRNAs) and DNA methylation, is not only a manually curated collection and annotation of experimentally supported lncRNAs-DNA methylation associations but also a platform that effectively integrates tools for calculating and identifying the differentially methylated lncRNAs and protein-coding genes (PCGs) in diverse human diseases. The resource provides: (i) advanced search possibilities, e.g. retrieval of the database by searching the lncRNA symbol of interest, DNA methylation patterns, regulatory mechanisms and disease types; (ii) abundant computationally calculated DNA methylation array profiles for the lncRNAs and PCGs; (iii) the prognostic values for each hit transcript calculated from the patients clinical data; (iv) a genome browser to display the DNA methylation landscape of the lncRNA transcripts for a specific type of disease; (v) tools to re-annotate probes to lncRNA loci and identify the differential methylation patterns for lncRNAs and PCGs with user-supplied external datasets; (vi) an R package (LncDM) to complete the differentially methylated lncRNAs identification and visualization with local computers. Lnc2Meth provides a timely and valuable resource that can be applied to significantly expand our understanding of the regulatory relationships between lncRNAs and DNA methylation in various human diseases.

MSDD: a manually curated database of experimentally supported associations among miRNAs, SNPs and human diseases.

The MiRNA SNP Disease Database (MSDD, http://www.bio-bigdata.com/msdd/) is a manually curated database that provides comprehensive experimentally supported associations among microRNAs (miRNAs), single nucleotide polymorphisms (SNPs) and human diseases. SNPs in miRNA-related functional regions such as mature miRNAs, promoter regions, pri-miRNAs, pre-miRNAs and target gene 3'-UTRs, collectively called 'miRSNPs', represent a novel category of functional molecules. miRSNPs can lead to miRNA and its target gene dysregulation, and resulting in susceptibility to or onset of human diseases. A curated collection and summary of miRSNP-associated diseases is essential for a thorough understanding of the mechanisms and functions of miRSNPs. Here, we describe MSDD, which currently documents 525 associations among 182 human miRNAs, 197 SNPs, 153 genes and 164 human diseases through a review of more than 2000 published papers. Each association incorporates information on the miRNAs, SNPs, miRNA target genes and disease names, SNP locations and alleles, the miRNA dysfunctional pattern, experimental techniques, a brief functional description, the original reference and additional annotation. MSDD provides a user-friendly interface to conveniently browse, retrieve, download and submit novel data. MSDD will significantly improve our understanding of miRNA dysfunction in disease, and thus, MSDD has the potential to serve as a timely and valuable resource.

LincSNP 2.0: an updated database for linking disease-associated SNPs to human long non-coding RNAs and their TFBSs.

We describe LincSNP 2.0 (http://bioinfo.hrbmu.edu.cn/LincSNP), an updated database that is used specifically to store and annotate disease-associated single nucleotide polymorphisms (SNPs) in human long non-coding RNAs (lncRNAs) and their transcription factor binding sites (TFBSs). In LincSNP 2.0, we have updated the database with more data and several new features, including (i) expanding disease-associated SNPs in human lncRNAs; (ii) identifying disease-associated SNPs in lncRNA TFBSs; (iii) updating LD-SNPs from the 1000 Genomes Project; and (iv) collecting more experimentally supported SNP-lncRNA-disease associations. Furthermore, we developed three flexible online tools to retrieve and analyze the data. Linc-Mart is a convenient way for users to customize their own data. Linc-Browse is a tool for all data visualization. Linc-Score predicts the associations between lncRNA and disease. In addition, we provided users a newly designed, user-friendly interface to search and download all the data in LincSNP 2.0 and we also provided an interface to submit novel data into the database. LincSNP 2.0 is a continually updated database and will serve as an important resource for investigating the functions and mechanisms of lncRNAs in human diseases.