Search

Long non-coding RNAs (lncRNAs) act as versatile regulators of many biological processes and play vital roles in various diseases. lncRNASNP is dedicated to providing a comprehensive repository of single nucleotide polymorphisms (SNPs) and somatic mutations in lncRNAs and their impacts on lncRNA structure and function. Since the last release in 2018, there has been a huge increase in the number of variants and lncRNAs. Thus, we updated the lncRNASNP to version 3 by expanding the species to eight eukaryotic species (human, chimpanzee, pig, mouse, rat, chicken, zebrafish, and fruitfly), updating the data and adding several new features. SNPs in lncRNASNP have increased from 11 181 387 to 67 513 785. The human mutations have increased from 1 174 768 to 2 387 685, including 1 031 639 TCGA mutations and 1 356 046 CosmicNCVs. Compared with the last release, updated and new features in lncRNASNP v3 include (i) SNPs in lncRNAs and their impacts on lncRNAs for eight species, (ii) SNP effects on miRNA−lncRNA interactions for eight species, (iii) lncRNA expression profiles for six species, (iv) disease & GWAS-associated lncRNAs and variants, (v) experimental & predicted lncRNAs and drug target associations and (vi) SNP effects on lncRNA expression (eQTL) across tumor & normal tissues. The lncRNASNP v3 is freely available at http://gong_lab.hzau.edu.cn/lncRNASNP3/.

Immune checkpoint blockade (ICB) therapy provides remarkable clinical benefits for multiple cancer types. Much work is currently being conducted to investigate the mechanisms of ICB therapy at the transcriptional level. Integrating the data produced by these studies will help us give more insight into the transcriptomic features of ICB therapy. We collected the transcriptome and clinical data of ICB-treated patient samples from the Gene Expression Omnibus, ArrayExpress, The Cancer Genome Atlas, and dbGaP databases. On the basis of the clinical information, all samples are initially classified into response/nonresponse or pretreatment/on-treatment groups. Differential expression, pathway enrichment, and immune cell infiltration analyses are performed between the samples from different groups. We also introduce the Response Score (RS) calculated by integrating the variability degree and the frequency of the dysregulated genes in the responders to evaluate the impact of gene expression on the response. Finally, all the abovementioned contents are integrated into the ICBatlas database. ICBatlas provides the transcriptome features of ICB therapy through the analysis of 1,515 ICB-treated samples from 25 studies across nine cancer types. The data in ICBatlas include clinical outcomes, treatment-related genes, biological pathways, and immune cell infiltration. Users can investigate the abovementioned transcriptome features in the response (R vs. NR) or treatment (Pre vs. On) modules at the data set, cancer type, or immune checkpoint level and compare the degree of gene impact on the response in the RS module. ICBatlas is the first database to show the transcriptome features on ICB therapy in human cancers and freely available at http://bioinfo.life.hust.edu.cn/ICBatlas/.

Supplementary Table from ICBatlas: A Comprehensive Resource for Depicting Immune Checkpoint Blockade Therapy Characteristics from Transcriptome Profiles

Supplementary Data from ICBatlas: A Comprehensive Resource for Depicting Immune Checkpoint Blockade Therapy Characteristics from Transcriptome Profiles

Immune checkpoint blockade (ICB) therapy provides remarkable clinical benefits for multiple cancer types. Much work is currently being conducted to investigate the mechanisms of ICB therapy at the transcriptional level. Integrating the data produced by these studies will help us give more insight into the transcriptomic features of ICB therapy. We collected the transcriptome and clinical data of ICB-treated patient samples from the Gene Expression Omnibus, ArrayExpress, The Cancer Genome Atlas, and dbGaP databases. On the basis of the clinical information, all samples are initially classified into response/nonresponse or pretreatment/on-treatment groups. Differential expression, pathway enrichment, and immune cell infiltration analyses are performed between the samples from different groups. We also introduce the Response Score (RS) calculated by integrating the variability degree and the frequency of the dysregulated genes in the responders to evaluate the impact of gene expression on the response. Finally, all the abovementioned contents are integrated into the ICBatlas database. ICBatlas provides the transcriptome features of ICB therapy through the analysis of 1,515 ICB-treated samples from 25 studies across nine cancer types. The data in ICBatlas include clinical outcomes, treatment-related genes, biological pathways, and immune cell infiltration. Users can investigate the abovementioned transcriptome features in the response (R vs. NR) or treatment (Pre vs. On) modules at the data set, cancer type, or immune checkpoint level and compare the degree of gene impact on the response in the RS module. ICBatlas is the first database to show the transcriptome features on ICB therapy in human cancers and freely available at http://bioinfo.life.hust.edu.cn/ICBatlas/.

Supplementary Figure from ICBatlas: A Comprehensive Resource for Depicting Immune Checkpoint Blockade Therapy Characteristics from Transcriptome Profiles

Motivation Immune cells are important components of the immune system and are crucial for disease initiation, progression, prognosis and survival. Although several computational methods have been designed for predicting the abundance of immune cells, very few tools are applicable to mouse. Given that, mouse is the most widely used animal model in biomedical research, there is an urgent need to develop a precise algorithm for predicting mouse immune cells. Results We developed a tool named Immune Cell Abundance Identifier for mouse (ImmuCellAI-mouse), for estimating the abundance of 36 immune cell (sub)types from gene expression data in a hierarchical strategy of three layers. Reference expression profiles and robust marker gene sets of immune cell types were curated. The abundance of cells in three layers was predicted separately by calculating the ssGSEA enrichment score of the expression deviation profile per cell type. Benchmark results showed high accuracy of ImmuCellAI-mouse in predicting most immune cell types, with correlation coefficients between predicted value and real cell proportion of most cell types being larger than 0.8. We applied ImmuCellAI-mouse to a mouse breast tumor dataset and revealed the dynamic change of immune cell infiltration during treatment, which is consistent with the findings of the original study but with more details. We also constructed an online server for ImmuCellAI-mouse, on which users can upload expression matrices for analysis. ImmuCellAI-mouse will be a useful tool for studying the immune microenvironment, cancer immunology and immunotherapy in mouse models, providing an indispensable supplement for human disease studies. Availability and implementation Software is available at http://bioinfo.life.hust.edu.cn/ImmuCellAI-mouse/. Supplementary information Supplementary data are available at Bioinformatics online.

Cancer initiation and progression are likely caused by the dysregulation of biological pathways. Gene set analysis (GSA) could improve the signal-to-noise ratio and identify potential biological insights on the gene set level. However, platforms exploring cancer multi-omics data using GSA methods are lacking. In this study, we upgraded our GSCALite to GSCA (gene set cancer analysis, http://bioinfo.life.hust.edu.cn/GSCA) for cancer GSA at genomic, pharmacogenomic and immunogenomic levels. In this improved GSCA, we integrated expression, mutation, drug sensitivity and clinical data from four public data sources for 33 cancer types. We introduced useful features to GSCA, including associations between immune infiltration with gene expression and genomic variations, and associations between gene set expression/mutation and clinical outcomes. GSCA has four main functional modules for cancer GSA to explore, analyze and visualize expression, genomic variations, tumor immune infiltration, drug sensitivity and their associations with clinical outcomes. We used case studies of three gene sets: (i) seven cell cycle genes, (ii) tumor suppressor genes of PI3K pathway and (iii) oncogenes of PI3K pathway to prove the advantage of GSCA over single gene analysis. We found novel associations of gene set expression and mutation with clinical outcomes in different cancer types on gene set level, while on single gene analysis level, they are not significant associations. In conclusion, GSCA is a user-friendly web server and a useful resource for conducting hypothesis tests by using GSA methods at genomic, pharmacogenomic and immunogenomic levels.

Although many studies have investigated functional molecules in extracellular vesicles (EVs), the exact number of ribonucleic acid molecules in a single-EV is unknown. Therefore, it is critical to explore the transcriptomic features and heterogeneity at the level of a single-EV. Here, using the 10x Genomics platform, the RNA cargos are profiled in single EVs derived from human K562 and mesenchymal stem cells. The key steps are labeling intact EVs using calcein-AM, detecting the EV concentration via flow cytometry, and using the CB2 algorithm with adaptive thresholds to effectively distinguish real EVs from background. The gene number in a single-EV varied from 6 to 148, with a mean of 52. Ribosomal genes, mitochondrial genes, and eukaryotic translation elongation factor 1 alpha has a high EV percentage in all EV samples. Hemoglobin genes are uniquely highly expressed in K562-EVs, and cytoskeleton genes are enriched in MSC-EVs. Ten or more clusters with different marker genes in each single-EV dataset demonstrated EV heterogeneity. Moreover, integrating EVs and their parental cells reveal both EVs and cells in each cluster, indicating different cell origins of various EVs. To the best of the author's knowledge, this study provides the first high-throughput transcriptome at the single-EV level and improves the understanding of EVs.

Transcription factors (TFs) as key regulators play crucial roles in biological processes. The identification of TF–target regulatory relationships is a key step for revealing functions of TFs and their regulations on gene expression. The accumulated data of chromatin immunoprecipitation sequencing (ChIP-seq) provide great opportunities to discover the TF–target regulations across different conditions. In this study, we constructed a database named hTFtarget, which integrated huge human TF target resources (7190 ChIP-seq samples of 659 TFs and high-confidence binding sites of 699 TFs) and epigenetic modification information to predict accurate TF–target regulations. hTFtarget offers the following functions for users to explore TF–target regulations: (1) browse or search general targets of a query TF across datasets; (2) browse TF–target regulations for a query TF in a specific dataset or tissue; (3) search potential TFs for a given target gene or non-coding RNA; (4) investigate co-association between TFs in cell lines; (5) explore potential co-regulations for given target genes or TFs; (6) predict candidate TF binding sites on given DNA sequences; (7) visualize ChIP-seq peaks for different TFs and conditions in a genome browser. hTFtarget provides a comprehensive, reliable and user-friendly resource for exploring human TF–target regulations, which will be very useful for a wide range of users in the TF and gene expression regulation community. hTFtarget is available at http://bioinfo.life.hust.edu.cn/hTFtarget.

Background: Cytogenetically normal acute myeloid leukemia (CN-AML) is a large proportion of AMLs with diverse prognostic outcomes. Identifying membrane protein genes as prognostic factors to stratify CN-AML patients will be critical to improve their outcomes. Purpose: This study aims to identify prognostic factors to stratify CN-AML patients to choose better treatments and improve their outcomes. Methods: CN-AML data were from TCGA cohort (n = 79) and four GEO datasets. We identified independent prognostic genes by Cox regression and Kaplan-Meier methods, and constructed linear regression model using LASSO algorithm. The prediction error curve was calculated using R package "pec". Results: Based on independent prognostic membrane genes, we constructed a regression model for CN-AML prognosis prediction: score = (0.0492 * CD52) - (0.0018 * CD96) + (0.0131 * EMP1) + (0.2058 * TSPAN2) + (0.0234 * STAB1) - (0.3658 * MBTPS1), which was named as MPG6 (6-Membrane Protein Gene) score. Tested in multiple CN-AML datasets, consistent results showed that CN-AML patients with high MPG6 score had poor survival, higher WBC count and shorter EFS. Comparing with other reported scoring models, the benchmark result of MPG6 achieved better association with survival in multiple cohorts. Moreover, by combining with other clinical indicators in CN-AML, MPG6 could improve the performance of survival prediction and serve as a robust prognostic factor. Conclusions: We identified the MPG6 score as a stable indicator with great potential for clinical application in risk stratification and outcome prediction in CN-AML.

Most relapsed chronic myeloid leukemia (CML) patients after tyrosine kinase inhibitor (TKI) discontinuation are in a chronic phase and could achieve remission through restarting the TKI treatment. Here we reported a case of sudden lymphoid blast crisis after 67 months of TKI discontinuation and depicted the patient by DNA and RNA sequencing to investigate intrinsic molecular features. The mutations of TGFBR2 and PCNT and the dysregulations of TGF-β and other pathways might accelerate the B cell transformation, which may serve as a blast crisis risk indicator of CML. Single-cell transcriptome data revealed that several clusters of immature B cells and late pro-B cells presented clone evolution during the treatment. After failing multiple lines of TKIs, conditioning chemotherapies and chimeric antigen receptor T cells (CAR-T) targeting CD19 and CD22 were performed to achieve remission. In conclusion, we report the first case of a CML patient with sudden lymphoid blast crisis after a long treatment-free remission and additional gene abnormalities other than BCR-ABL1 might participate in the progression, which need to be closely monitored, and CAR-T could be a solution to the chemoresistant progression.

Background Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide, but its regulatory mechanism remains unclear. Although many TFs and miRNAs are reported to be important in HCC, their co-regulation and FFL modules in HCC development are needed to be investigated. Methods The feed-forward loop (FFL) regulatory module was identified by analyzing the miRNA and transcription factor co-regulatory network for differentially expressed genes between tumors and matched adjacent tissue samples. Gene expression and regulatory role of HCC development by key FFL in vitro and in vivo were validated by qPCR, Western blotting, cell proliferation assay, migration and invasion assay and experiments in nude mice with hepatoma xenografts. Results Here, by bioinformatics analysis, we identified FFL regulatory module miR-9-5p/FOXO1/CPEB3 may play critical roles in HCC progression. Gain- and loss-of-function studies demonstrated that miR-9-5p promote hepatocarcinogenesis, while FOXO1 and CPEB3 inhibit hepatocarcinoma growth. Furthermore, CPEB3 was firstly identified as a direct downstream target of miR-9-5p and FOXO1 by luciferase reporter assay and ChIP-Seq data, which was negatively regulated by miR-9-5p and activated by FOXO1. Following, the miR-9-5p/FOXO1/CPEB3 FFL was associated with poor prognosis and promoted cell growth and tumorigenesis of HCC in both in vitro and in vivo experiments. Conclusion Our study newly identified the existence of miR-9-5p/FOXO1/CPEB3 FFL and revealed its regulatory role in HCC progression, which may represent a new potential therapeutic target for cancer treatment.

Extracellular vesicles (EVs) packing various molecules play vital roles in intercellular communication. Non-coding RNAs (ncRNAs) are important functional molecules and biomarkers in EVs. A comprehensive investigation of ncRNAs expression in EVs under different conditions is a fundamental step for functional discovery and application of EVs. Here, we curated 2030 small RNA-seq datasets for human EVs (1506 sEV and 524 lEV) in 24 conditions and over 40 diseases. We performed a unified reads dynamic assignment algorithm (RDAA) considering mismatch and multi-mapping reads to quantify the expression profiles of seven ncRNA types (miRNA, snoRNA, piRNA, snRNA, rRNA, tRNA and Y RNA). We constructed EVAtlas (http://bioinfo.life.hust.edu.cn/EVAtlas), a comprehensive database for ncRNA expression in EVs with four functional modules: (i) browse and compare the distribution of ncRNAs in EVs from 24 conditions and eight sources (plasma, serum, saliva, urine, sperm, breast milk, primary cell and cell line); (ii) prioritize candidate ncRNAs in condition related tissues based on their expression; (iii) explore the specifically expressed ncRNAs in EVs from 24 conditions; (iv) investigate ncRNA functions, related drugs, target genes and EVs isolation methods. EVAtlas contains the most comprehensive ncRNA expression in EVs and will be a key resource in this field.

The BIG Data Center at Beijing Institute of Genomics (BIG) of the Chinese Academy of Sciences provides a suite of database resources in support of worldwide research activities in both academia and industry. With the vast amounts of multi-omics data generated at unprecedented scales and rates, the BIG Data Center is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. Resources with significant updates in the past year include BioProject (a biological project library), BioSample (a biological sample library), Genome Sequence Archive (GSA, a data repository for archiving raw sequence reads), Genome Warehouse (GWH, a centralized resource housing genome-scale data), Genome Variation Map (GVM, a public repository of genome variations), Science Wikis (a catalog of biological knowledge wikis for community annotations) and IC4R (Information Commons for Rice). Newly released resources include EWAS Atlas (a knowledgebase of epigenome-wide association studies), iDog (an integrated omics data resource for dog) and RNA editing resources (for editome-disease associations and plant RNA editosome, respectively). To promote biodiversity and health big data sharing around the world, the Open Biodiversity and Health Big Data (BHBD) initiative is introduced. All of these resources are publicly accessible at http://bigd.big.ac.cn.

Extracellular vesicles (EVs), such as exosomes and microvesicles, acted as cell-to-cell communication vectors and potential biomarkers for diseases. microRNAs (miRNAs) are the most well studied molecules in EVs, thus a comprehensive investigation of miRNA expression profiles in EVs will be helpful to explore their functions and biomarkers. We curated 462 small RNA sequencing samples of EVs from 17 sources/diseases and constructed the EVmiRNA database (http://bioinfo.life.hust.edu.cn/EVmiRNA) to show the miRNA expression profiles. We found >1000 miRNAs expressed in these EVs and detected specific miRNAs for EVs of each source/disease. EVmiRNA provides three functional modules: (i) the miRNA expression profiles and the sample information of EVs from different sources (such as blood, breast milk etc.); (ii) the specifically expressed miRNAs in different EVs that would be helpful for biomarker identification; (iii) the miRNA annotations including the miRNA expression in EVs and TCGA cancer types, miRNA pathway regulations as well as miRNA function and publications. EVmiRNA has a user-friendly web interface with powerful browse and search functions, as well as data downloading. It is the first database focusing on miRNA expression profiles in EVs and will be useful for the research and application community of EV biomarker, miRNA function and liquid biopsy.

Summary Transcription factors (TFs) and microRNAs (miRNAs) are two kinds of important regulators for transcriptional and post-transcriptional regulations. Understanding cross-talks between the two regulators and their targets is critical to reveal complex molecular regulatory mechanisms. Here, we developed FFLtool, a web server for detecting potential feed forward loop (FFL) of TF-miRNA-target regulation in human. In FFLtool, we integrated comprehensive regulations of TF-target and miRNA-target, and developed two functional modules: (i) The ‘FFL Analysis’ module can detect potential FFLs and internal regulatory networks in a user-defined gene set. FFLtool also provides three levels of evidence to illustrate the reliability for each FFL and enrichment functions for co-target genes of the same TF and miRNA; (ii) The ‘Browse FFLs’ module displays FFLs comprised of differentially or specifically expressed TFs and miRNAs and their target genes in cancers. FFLtool is a valuable resource for investigating gene expression regulation and mechanism study in biological processes and diseases. Availability and implementation FFLtool is available on http://bioinfo.life.hust.edu.cn/FFLtool/. Supplementary information Supplementary data are available at Bioinformatics online.

Human specifically expressed genes (SEGs) usually serve as potential biomarkers for disease diagnosis and treatment. However, the regulation underlying their specific expression remains to be revealed. In this study, we constructed SEG regulation database (SEGreg; available at http://bioinfo.life.hust.edu.cn/SEGreg) for showing SEGs and their transcription factors (TFs) and microRNA (miRNA) regulations under different physiological conditions, which include normal tissue, cancer tissue and cell line. In total, SEGreg collected 6387, 1451, 4506 and 5320 SEGs from expression profiles of 34 cancer types and 55 tissues of The Cancer Genome Atlas, Cancer Cell Line Encyclopedia, Human Body Map and Genotype-Tissue Expression databases/projects, respectively. The cancer or tissue corresponding expressed miRNAs and TFs were identified from miRNA and gene expression profiles, and their targets were collected from several public resources. Then the regulatory networks of all SEGs were constructed and integrated into SEGreg. Through a user-friendly interface, users can browse and search SEGreg by gene name, data source, tissue, cancer type and regulators. In summary, SEGreg is a specialized resource to explore SEGs and their regulations, which provides clues to reveal the mechanisms of carcinogenesis and biological processes.

Long non-coding RNAs (lncRNAs) are emerging as important regulators in different biological processes through various ways. Because the related data, especially mutations in cancers, increased sharply, we updated the lncRNASNP to version 2 (http://bioinfo.life.hust.edu.cn/lncRNASNP2). lncRNASNP2 provides comprehensive information of SNPs and mutations in lncRNAs, as well as their impacts on lncRNA structure and function. lncRNASNP2 contains 7260238 SNPs on 141353 human lncRNA transcripts and 3921448 SNPs on 117405 mouse lncRNA transcripts. Besides the SNP information in the first version, the following new features were developed to improve the lncRNASNP2. (i) noncoding variants from COSMIC cancer data (859534) in lncRNAs and their effects on lncRNA structure and function; (ii) TCGA cancer mutations (315234) in lncRNAs and their impacts; (iii) lncRNA expression profiling of 20 cancer types in both tumor and its adjacent samples; (iv) expanded lncRNA-associated diseases; (v) optimized the results about lncRNAs structure change induced by variants; (vi) reduced false positives in miRNA and lncRNA interaction results. Furthermore, we developed online tools for users to analyze new variants in lncRNA. We aim to maintain the lncRNASNP as a useful resource for lncRNAs and their variants.

The distribution and abundance of immune cells, particularly T‐cell subsets, play pivotal roles in cancer immunology and therapy. T cells have many subsets with specific function and current methods are limited in estimating them, thus, a method for predicting comprehensive T‐cell subsets is urgently needed in cancer immunology research. Here, Immune Cell Abundance Identifier (ImmuCellAI), a gene set signature‐based method, is introduced for precisely estimating the abundance of 24 immune cell types including 18 T‐cell subsets, from gene expression data. Performance evaluation on both the sequencing data with flow cytometry results and public expression data indicate that ImmuCellAI can estimate the abundance of immune cells with superior accuracy to other methods especially on many T‐cell subsets. Application of ImmuCellAI to immunotherapy datasets reveals that the abundance of dendritic cells, cytotoxic T, and gamma delta T cells is significantly higher both in comparisons of on‐treatment versus pre‐treatment and responders versus non‐responders. Meanwhile, an ImmuCellAI result‐based model is built for predicting the immunotherapy response with high accuracy (area under curve 0.80–0.91). These results demonstrate the powerful and unique function of ImmuCellAI in tumor immune infiltration estimation and immunotherapy response prediction., Immune Cell Abundance Identifier (ImmuCellAI) is a gene set signature‐based method for precisely estimating the abundance of 24 immune cell types including 18 T‐cell subsets. Application of ImmuCellAI to immunotherapy datasets reveals the dynamic change of immune cell abundance. An ImmuCellAI result‐based model for predicting the immunotherapy response achieves high accuracy with area under curve 0.80–0.91.

The National Genomics Data Center (NGDC) provides a suite of database resources to support worldwide research activities in both academia and industry. With the rapid advancements in higher-throughput and lower-cost sequencing technologies and accordingly the huge volume of multi-omics data generated at exponential scales and rates, NGDC is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. In the past year, efforts for update have been mainly devoted to BioProject, BioSample, GSA, GWH, GVM, NONCODE, LncBook, EWAS Atlas and IC4R. Newly released resources include three human genome databases (PGG.SNV, PGG.Han and CGVD), eLMSG, EWAS Data Hub, GWAS Atlas, iSheep and PADS Arsenal. In addition, four web services, namely, eGPS Cloud, BIG Search, BIG Submission and BIG SSO, have been significantly improved and enhanced. All of these resources along with their services are publicly accessible at https://bigd.big.ac.cn.

Alternative polyadenylation (APA) is an important post-transcriptional regulation that recognizes different polyadenylation signals (PASs), resulting in transcripts with different 3′ untranslated regions, thereby influencing a series of biological processes and functions. Recent studies have revealed that some single nucleotide polymorphisms (SNPs) could contribute to tumorigenesis and development through dysregulating APA. However, the associations between SNPs and APA in human cancers remain largely unknown. Here, using genotype and APA data of 9082 samples from The Cancer Genome Atlas (TCGA) and The Cancer 3′UTR Altas (TC3A), we systematically identified SNPs affecting APA events across 32 cancer types and defined them as APA quantitative trait loci (apaQTLs). As a result, a total of 467 942 cis-apaQTLs and 30 721 trans-apaQTLs were identified. By integrating apaQTLs with survival and genome-wide association studies (GWAS) data, we further identified 2154 apaQTLs associated with patient survival time and 151 342 apaQTLs located in GWAS loci. In addition, we designed an online tool to predict the effects of SNPs on PASs by utilizing PAS motif prediction tool. Finally, we developed SNP2APA, a user-friendly and intuitive database (http://gong_lab.hzau.edu.cn/SNP2APA/) for data browsing, searching, and downloading. SNP2APA will significantly improve our understanding of genetic variants and APA in human cancers.