DeepFusion: A deep bimodal information fusion network for unraveling protein-RNA interactions using in vivo RNA structures.

RNA-binding proteins (RBPs) are key post-transcriptional regulators, and the malfunctions of RBP-RNA binding lead to diverse human diseases. However, prediction of RBP binding sites is largely based on RNA sequence features, whereas in vivo RNA structural features based on high-throughput sequencing are rarely incorporated. Here, we designed a deep bimodal information fusion network called DeepFusion for unraveling protein-RNA interactions by incorporating structural features derived from DMS-seq data. DeepFusion integrates two sub-models to extract local motif-like information and long-term context information. We show that DeepFusion performs best compared with other cutting-edge methods with only sequence inputs on two datasets. DeepFusion's performance is further improved with bimodal input after adding in vivo DMS-seq structural features. Furthermore, DeepFusion can be used for analyzing RNA degradation, demonstrating significantly different RBP-binding scores in genes with slow degradation rates versus those with rapid degradation rates. DeepFusion thus provides enhanced abilities for further analysis of functional RNAs. DeepFusion's code and data are available at http://bioinfo.org/deepfusion/.

Optical Genome Mapping Reveals the Landscape of Structural Variations and Their Clinical Significance in HBOC-Related Breast Cancer.

BACKGROUND: Structural variations (SVs) are common genetic alterations in the human genome. However, the profile and clinical relevance of SVs in patients with hereditary breast and ovarian cancer (HBOC) syndrome (germline BRCA1/2 mutations) remains to be fully elucidated. METHODS: Twenty HBOC-related cancer samples (5 breast and 15 ovarian cancers) were studied by optical genome mapping (OGM) and next-generation sequencing (NGS) assays. RESULTS: The SV landscape in the 5 HBOC-related breast cancer samples was comprehensively investigated to determine the impact of intratumor SV heterogeneity on clinicopathological features and on the pattern of genetic alteration. SVs and copy number variations (CNVs) were common genetic events in HBOC-related breast cancer, with a median of 212 SVs and 107 CNVs per sample. The most frequently detected type of SV was insertion, followed by deletion. The 5 HBOC-related breast cancer samples were divided into SVhigh and SVlow groups according to the intratumor heterogeneity of SVs. SVhigh tumors were associated with higher Ki-67 expression, higher homologous recombination deficiency (HRD) scores, more mutated genes, and altered signaling pathways. Moreover, 60% of the HBOC-related breast cancer samples displayed chromothripsis, and 8 novel gene fusion events were identified by OGM and validated by transcriptome data. CONCLUSIONS: These findings suggest that OGM is a promising tool for the detection of SVs and CNVs in HBOC-related breast cancer. Furthermore, OGM can efficiently characterize chromothripsis events and novel gene fusions. SVhigh HBOC-related breast cancers were associated with unfavorable clinicopathological features. SVs may therefore have predictive and therapeutic significance for HBOC-related breast cancers in the clinic.

Corrigendum to "TREAT: Therapeutic RNAs exploration inspired by artificial intelligence technology" [Comput Struct Biotechnol J 20 (2022) 5680-5689].

[This corrects the article DOI: 10.1016/j.scr.2023.103115.].

Dominant neoantigen verification in hepatocellular carcinoma by a single-plasmid system coexpressing patient HLA and antigen.

BACKGROUND: Previous studies confirmed that most neoantigens predicted by algorithms do not work in clinical practice, and experimental validations remain indispensable for confirming immunogenic neoantigens. In this study, we identified the potential neoantigens with tetramer staining, and established the Co-HA system, a single-plasmid system coexpressing patient human leukocyte antigen (HLA) and antigen, to detect the immunogenicity of neoantigens and verify new dominant hepatocellular carcinoma (HCC) neoantigens. METHODS: First, we enrolled 14 patients with HCC for next-generation sequencing for variation calling and predicting potential neoantigens. Then, the Co-HA system was established. To test the feasibility of the system, we constructed target cells coexpressing HLA-A*11:01 and the reported KRAS G12D neoantigen as well as specific T-cell receptor (TCR)-T cells. The specific cytotoxicity generated by this neoantigen was shown using the Co-HA system. Moreover, potential HCC-dominant neoantigens were screened out by tetramer staining and validated by the Co-HA system using methods including flow cytometry, enzyme-linked immunospot assay and ELISA. Finally, antitumor test in mouse mode and TCR sequencing were performed to further evaluate the dominant neoantigen. RESULTS: First, 2875 somatic mutations in 14 patients with HCC were identified. The main base substitutions were C>T/G>A transitions, and the main mutational signatures were 4, 1 and 16. The high-frequency mutated genes included HMCN1, TTN and TP53. Then, 541 potential neoantigens were predicted. Importantly, 19 of the 23 potential neoantigens in tumor tissues also existed in portal vein tumor thrombi. Moreover, 37 predicted neoantigens restricted by HLA-A*11:01, HLA-A*24:02 or HLA-A*02:01 were performed by tetramer staining to screen out potential HCC-dominant neoantigens. HLA-A*24:02-restricted epitope 5'-FYAFSCYYDL-3' and HLA-A*02:01-restricted epitope 5'-WVWCMSPTI-3' demonstrated strong immunogenicity in HCC, as verified by the Co-HA system. Finally, the antitumor efficacy of 5'-FYAFSCYYDL-3'-specific T cells was verified in the B-NDG-B2mtm1Fcrntm1(mB2m) mouse and their specific TCRs were successfully identified. CONCLUSION: We found the dominant neoantigens with high immunogenicity in HCC, which were verified with the Co-HA system.

Drug screening and biomarker gene investigation in cancer therapy through the human transcriptional regulatory network.

A complex and vast biological network regulates all biological functions in the human body in a sophisticated manner, and abnormalities in this network can lead to disease and even cancer. The construction of a high-quality human molecular interaction network is possible with the development of experimental techniques that facilitate the interpretation of the mechanisms of drug treatment for cancer. We collected 11 molecular interaction databases based on experimental sources and constructed a human protein-protein interaction (PPI) network and a human transcriptional regulatory network (HTRN). A random walk-based graph embedding method was used to calculate the diffusion profiles of drugs and cancers, and a pipeline was constructed by using five similarity comparison metrics combined with a rank aggregation algorithm, which can be implemented for drug screening and biomarker gene prediction. Taking NSCLC as an example, curcumin was identified as a potentially promising anticancer drug from 5450 natural small molecules, and combined with differentially expressed genes, survival analysis, and topological ranking, we obtained BIRC5 (survivin), which is both a biomarker for NSCLC and a key target for curcumin. Finally, the binding mode of curcumin and survivin was explored using molecular docking. This work has a guiding significance for antitumor drug screening and the identification of tumor markers.

Biological knowledge graph-guided investigation of immune therapy response in cancer with graph neural network.

The determination of transcriptome profiles that mediate immune therapy in cancer remains a major clinical and biological challenge. Despite responses induced by immune-check points inhibitors (ICIs) in diverse tumor types and all the big breakthroughs in cancer immunotherapy, most patients with solid tumors do not respond to ICI therapies. It still remains a big challenge to predict the ICI treatment response. Here, we propose a framework with multiple prior knowledge networks guided for immune checkpoints inhibitors prediction-DeepOmix-ICI (or ICInet for short). ICInet can predict the immune therapy response by leveraging geometric deep learning and prior biological knowledge graphs of gene-gene interactions. Here, we demonstrate more than 600 ICI-treated patients with ICI response data and gene expression profile to apply on ICInet. ICInet was used for ICI therapy responses prediciton across different cancer types-melanoma, gastric cancer and bladder cancer, which includes 7 cohorts from different data sources. ICInet is able to robustly generalize into multiple cancer types. Moreover, the performance of ICInet in those cancer types can outperform other ICI biomarkers in the clinic. Our model [area under the curve (AUC = 0.85)] generally outperformed other measures, including tumor mutational burden (AUC = 0.62) and programmed cell death ligand-1 score (AUC = 0.74). Therefore, our study presents a prior-knowledge guided deep learning method to effectively select immunotherapy-response-associated biomarkers, thereby improving the prediction of immunotherapy response for precision oncology.

Genomic and transcriptomic profiling indicates the prognosis significance of mutational signature for TMB-high subtype in Chinese patients with gastric cancer.

INTRODUCTION: Gastric cancer (GC)is the third leading cause of cancer-related deaths in China and immunotherapy emerging as a revolutionary treatment for GC recently. Tumor mutational burden (TMB) is a predictive biomarker of immunotherapy in multiple cancers. However, the prognostic significance and subtype of TMB in GC is not fully understood. OBJECTIVES: This study aims to evaluate the prognostic value of TMB in Chinese GC and further classify TMB-high GC (GCTMB-H) patients combing with mutational signatures. METHODS: Genomic profiling of 435 cancer-gene panel was performed using 206 GC samples from Chinese people. Actionable genetic alterations were compared across all the samples to generate actionable subtyping. The prognostic value of TMB in Chinese GC was evaluated. Mutational signatures were analyzed on TMB-H subtype to stratify the prognosis of TMB. Transcriptomic analysis was applied to compare the distributed immunocytes among different subtypes. RESULTS: 88.3% (182/206) of GC samples had at least one mutation, while 45.1% (93/206) had at least one somatic copy number alteration (SCNA). 29.6% (61/206) of GC samples were TMB-H, including 13 MSI-H and 48 MSS tumors. According to distinct genetic alteration profiles of 69 actionable genes, we classified GC samples into eight molecular subtypes, including TMB-H, ERBB2 amplified, ATM mutated, BRCA2 mutated, CDKN2A/B deleted, PI3KCA mutated, KRAS mutated, and less-mutated subtype. TMB-H subtype presented a remarkable immune-activated phenotype as determined by transcriptomic analysis that was further validated in the TCGA GC cohort. GCTMB-H patients exhibited significantly better survival (P = 0.047). But Signature 1-high GCTMB-H patients had relatively worse prognosis (P = 0.0209, HR = 2.571) than Signature 1-low GCTMB-H patients from Chinese GC cohort, also validated in TCGA GC cohort, presenting highly activated carbohydrate, fatty acid or lipid metabolism. CONCLUSION: The Signature 1-high GCTMB-H could be a marker of poor prognosis and is associated with metabolism disorder.

TREAT: Therapeutic RNAs exploration inspired by artificial intelligence technology.

Recent advances in RNA engineering have enabled the development of RNA-based therapeutics for a broad spectrum of applications. Developing RNA therapeutics start with targeted RNA screening and move to the drug design and optimization. However, existing target screening tools ignore noncoding RNAs and their disease-relevant regulatory relationships. And designing therapeutic RNAs encounters high computational complexity of multi-objective optimization to overcome the immunogenicity, instability and inefficient translational production. To unlock the therapeutic potential of noncoding RNAs and enable one-stop screening and design of therapeutic RNAs, we have built the platform TREAT. It incorporates 43,087,953 regulatory relationships between coding and noncoding genes from 81 biological networks under different physiological conditions. TREAT introduces graph representation learning with Random Walk Diffusions to perform disease-relevant target screening, in addition to the commonly utilized Topological Degree and PageRank algorithms. Design and optimization of large RNAs or interfering RNAs are both available. To reduce the computational complexity of multi-objective optimization for large RNA, we stratified the features into local and global features. The local features are evaluated on the fixed-length or dynamic-length local bins, whereas the latter are inspired by AI language models of protein sequence. Then the global assessment is performed on refined candidates, thus reducing the enormous search space. Overall, TREAT is a one-stop platform for the screening and designing of therapeutic RNAs, with particular attention to noncoding RNAs and cutting-edge AI technology embedded, leading the progress of innovative therapeutics for challenging diseases. TREAT is freely accessible at https://rna.org.cn/treat.

Multi-modality artificial intelligence in digital pathology.

In common medical procedures, the time-consuming and expensive nature of obtaining test results plagues doctors and patients. Digital pathology research allows using computational technologies to manage data, presenting an opportunity to improve the efficiency of diagnosis and treatment. Artificial intelligence (AI) has a great advantage in the data analytics phase. Extensive research has shown that AI algorithms can produce more up-to-date and standardized conclusions for whole slide images. In conjunction with the development of high-throughput sequencing technologies, algorithms can integrate and analyze data from multiple modalities to explore the correspondence between morphological features and gene expression. This review investigates using the most popular image data, hematoxylin-eosin stained tissue slide images, to find a strategic solution for the imbalance of healthcare resources. The article focuses on the role that the development of deep learning technology has in assisting doctors' work and discusses the opportunities and challenges of AI.

High Throughput Transcriptome Data Analysis and Computational Verification Reveal Immunotherapy Biomarkers of Compound Kushen Injection for Treating Triple-Negative Breast Cancer.

BACKGROUND: Although notable therapeutic and prognostic benefits of compound kushen injection (CKI) have been found when it was used alone or in combination with chemotherapy or radiotherapy for triple-negative breast cancer (TNBC) treatment, the effects of CKI on TNBC microenvironment remain largely unclear. This study aims to construct and validate a predictive immunotherapy signature of CKI on TNBC. METHODS: The UPLC-Q-TOF-MS technology was firstly used to investigate major constituents of CKI. RNA sequencing data of CKI-perturbed TNBC cells were analyzed to detect differential expression genes (DEGs), and the GSVA algorithm was applied to explore significantly changed pathways regulated by CKI. Additionally, the ssGSEA algorithm was used to quantify immune cell abundance in TNBC patients, and these patients were classified into distinct immune infiltration subgroups by unsupervised clustering. Then, prognosis-related genes were screened from DEGs among these subgroups and were further overlapped with the DEGs regulated by CKI. Finally, a predictive immunotherapy signature of CKI on TNBC was constructed based on the LASSO regression algorithm to predict mortality risks of TNBC patients, and the signature was also validated in another TNBC cohort. RESULTS: Twenty-three chemical components in CKI were identified by UPLC-Q-TOF-MS analysis. A total of 3692 DEGs were detected in CKI-treated versus control groups, and CKI significantly activated biological processes associated with activation of T, natural killer and natural killer T cells. Three immune cell infiltration subgroups with 1593 DEGs were identified in TNBC patients. Then, two genes that can be down-regulated by CKI with hazard ratio (HR) > 1 and 26 genes that can be up-regulated by CKI with HR < 1 were selected as key immune- and prognosis-related genes regulated by CKI. Lastly, a five-gene prognostic signature comprising two risky genes (MARVELD2 and DYNC2I2) that can be down-regulated by CKI and three protective genes (RASSF2, FERMT3 and RASSF5) that can be up-regulated by CKI was developed, and it showed a good performance in both training and test sets. CONCLUSIONS: This study proposes a predictive immunotherapy signature of CKI on TNBC, which would provide more evidence for survival prediction and treatment guidance in TNBC as well as a paradigm for exploring immunotherapy biomarkers in compound medicines.

Single-cell Long Non-coding RNA Landscape of T Cells in Human Cancer Immunity.

The development of new biomarkers or therapeutic targets for cancer immunotherapies requires deep understanding of T cells. To date, the complete landscape and systematic characterization of long noncoding RNAs (lncRNAs) in T cells in cancer immunity are lacking. Here, by systematically analyzing full-length single-cell RNA sequencing (scRNA-seq) data of more than 20,000 libraries of T cells across three cancer types, we provided the first comprehensive catalog and the functional repertoires of lncRNAs in human T cells. Specifically, we developed a custom pipeline for de novotranscriptome assembly and obtained a novel lncRNA catalog containing 9433 genes. This increased the number of current human lncRNA catalog by 16% and nearly doubled the number of lncRNAs expressed in T cells. We found that a portion of expressed genes in single T cells were lncRNAs which had been overlooked by the majority of previous studies. Based on metacell maps constructed by the MetaCell algorithm that partitions scRNA-seq datasets into disjointed and homogenous groups of cells (metacells), 154 signature lncRNA genes were identified. They were associated with effector, exhausted, and regulatory T cell states. Moreover, 84 of them were functionally annotated based on the co-expression networks, indicating that lncRNAs might broadly participate in the regulation of T cell functions. Our findings provide a new point of view and resource for investigating the mechanisms of T cell regulation in cancer immunity as well as for novel cancer-immune biomarker development and cancer immunotherapies.

KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis.

Gene set enrichment (GSE) analysis plays an essential role in extracting biological insight from genome-scale experiments. ORA (overrepresentation analysis), FCS (functional class scoring), and PT (pathway topology) approaches are three generations of GSE methods along the timeline of development. Previous versions of KOBAS provided services based on just the ORA method. Here we presented version 3.0 of KOBAS, which is named KOBAS-i (short for KOBAS intelligent version). It introduced a novel machine learning-based method we published earlier, CGPS, which incorporates seven FCS tools and two PT tools into a single ensemble score and intelligently prioritizes the relevant biological pathways. In addition, KOBAS has expanded the downstream exploratory visualization for selecting and understanding the enriched results. The tool constructs a novel view of cirFunMap, which presents different enriched terms and their correlations in a landscape. Finally, based on the previous version's framework, KOBAS increased the number of supported species from 1327 to 5944. For an easier local run, it also provides a prebuilt Docker image that requires no installation, as a supplementary to the source code version. KOBAS can be freely accessed at http://kobas.cbi.pku.edu.cn, and a mirror site is available at http://bioinfo.org/kobas.

DeepOmix: A scalable and interpretable multi-omics deep learning framework and application in cancer survival analysis.

Integrative analysis of multi-omics data can elucidate valuable insights into complex molecular mechanisms for various diseases. However, due to their different modalities and high dimension, utilizing and integrating different types of omics data suffers from great challenges. There is an urgent need to develop a powerful method to improve survival prediction and detect functional gene modules from multi-omics data. To deal with these problems, we present DeepOmix (a scalable and interpretable multi-Omics Deep learning framework and application in cancer survival analysis), a flexible, scalable, and interpretable method for extracting relationships between the clinical survival time and multi-omics data based on a deep learning framework. DeepOmix enables the non-linear combination of variables from different omics datasets and incorporates prior biological information defined by users (such as signaling pathways and tissue networks). Benchmark experiments demonstrate that DeepOmix outperforms the other five cutting-edge prediction methods. Besides, Lower Grade Glioma (LGG) is taken as the case study to perform the prognosis prediction and illustrate the functional module nodes which are associated with the prognostic result in the prediction model.

ncFANs v2.0: an integrative platform for functional annotation of non-coding RNAs.

Increasing evidence proves the essential regulatory roles of non-coding RNAs (ncRNAs) in biological processes. However, characterizing the specific functions of ncRNAs remains a challenging task, owing to the intensive consumption of the experimental approaches. Here, we present an online platform ncFANs v2.0 that is a significantly enhanced version of our previous ncFANs to provide multiple computational methods for ncRNA functional annotation. Specifically, ncFANs v2.0 was updated to embed three functional modules, including ncFANs-NET, ncFANs-eLnc and ncFANs-CHIP. ncFANs-NET is a new module designed for data-free functional annotation based on four kinds of pre-built networks, including the co-expression network, co-methylation network, long non-coding RNA (lncRNA)-centric regulatory network and random forest-based network. ncFANs-eLnc enables the one-stop identification of enhancer-derived lncRNAs from the de novo assembled transcriptome based on the user-defined or our pre-annotated enhancers. Moreover, ncFANs-CHIP inherits the original functions for microarray data-based functional annotation and supports more chip types. We believe that our ncFANs v2.0 carries sufficient convenience and practicability for biological researchers and facilitates unraveling the regulatory mechanisms of ncRNAs. The ncFANs v2.0 server is freely available at http://bioinfo.org/ncfans or http://ncfans.gene.ac.

Rejuvenation remodels transcriptional network to improve rhizogenesis in mature Juglans tree.

Adventitious rooting of walnut species (Juglans L.) is known to be rather difficult, especially for mature trees. The adventitious root formation (ARF) capacities of mature trees can be significantly improved by rejuvenation. However, the underlying gene regulatory networks (GRNs) of rejuvenation remain largely unknown. To characterize such regulatory networks, we carried out the transcriptomic study using RNA samples of the cambia and peripheral tissues on the bottom of rejuvenated and mature walnut (Juglans hindsii × J. regia) cuttings during the ARF. The RNA sequencing data suggested that zeatin biosynthesis, energy metabolism and substance metabolism were activated by rejuvenation, whereas photosynthesis, fatty acid biosynthesis and the synthesis pathways for secondary metabolites were inhibited. The inter- and intra-module GRNs were constructed using differentially expressed genes. We identified 35 hub genes involved in five modules associated with ARF. Among these hub genes, particularly, beta-glucosidase-like (BGLs) family members involved in auxin metabolism were overexpressed at the early stage of the ARF. Furthermore, BGL12 from the cuttings of Juglans was overexpressed in Populus alba × P. glandulosa. Accelerated ARF and increased number of ARs were observed in the transgenic poplars. These results provide a high-resolution atlas of gene activity during ARF and help to uncover the regulatory modules associated with the ARF promoted by rejuvenation.

Imbalance of TGF-ß1/BMP-7 pathways induced by M2-polarized macrophages promotes hepatocellular carcinoma aggressiveness.

The transforming growth factor-beta (TGF-ß) signaling pathway is the predominant cytokine signaling pathway in the development and progression of hepatocellular carcinoma (HCC). Bone morphogenetic protein (BMP), another member of the TGF-ß superfamily, has been frequently found to participate in crosstalk with the TGF-ß pathway. However, the complex interaction between the TGF-ß and BMP pathways has not been fully elucidated in HCC. We found that the imbalance of TGF-ß1/BMP-7 pathways was associated with aggressive pathological features and poor clinical outcomes in HCC. The induction of the imbalance of TGF-ß1/BMP-7 pathways in HCC cells could significantly promote HCC cell invasion and stemness by increasing inhibitor of differentiation 1 (ID1) expression. We also found that the microRNA (miR)-17-92 cluster, originating from the extracellular vesicles (EVs) of M2-polarized tumor-associated macrophages (M2-TAMs), stimulated the imbalance of TGF-ß1/BMP-7 pathways in HCC cells by inducing TGF-ß type II receptor (TGFBR2) post-transcriptional silencing and inhibiting activin A receptor type 1 (ACVR1) post-translational ubiquitylation by targeting Smad ubiquitylation regulatory factor 1 (Smurf1). In vivo, short hairpin (sh)-MIR17HG and ACVR1 inhibitors profoundly attenuated HCC cell growth and metastasis by rectifying the imbalance of TGF-ß1/BMP-7 pathways. Therefore, we proposed that the imbalance of TGF-ß1/BMP-7 pathways is a feasible prognostic biomarker and recovering the imbalance of TGF-ß1/BMP-7 pathways might be a potential therapeutic strategy for HCC.

FangNet: Mining herb hidden knowledge from TCM clinical effective formulas using structure network algorithm.

The use of herbs to treat various human diseases has been recorded for thousands of years. In Asia's current medical system, numerous herbal formulas have been repeatedly verified to confirm their effectiveness in different periods, which is a great resource for drug innovation and discovery. Through the mining of these clinical effective formulas by network pharmacology and bioinformatics analysis, important biologically active ingredients derived from these natural products might be discovered. As modern medicine requires a combination of multiple drugs for the treatment of complex diseases, previously clinical formulas are also combinations of various herbs according to the main causes and accompanying symptoms. However, the herbs that play a major role in the treatment of diseases are always unclear. Therefore, how to rank each herb's relative importance and determine the core herbs, is the first step to assisting herb selection for active ingredients discovery. To solve this problem, we built the platform FangNet, which ranks all herbs on their relative topological importance using the PageRank algorithm, based on the constructed symptom-herb network from a collection of clinical empirical prescriptions. Three types of herb hidden knowledge, including herb importance rank, herb-herb co-occurrence, and associations to symptoms, were provided in an interactive visualization. Moreover, FangNet has designed role-based permission for teams to store, analyze, and jointly interpret their clinical formulas, in an easy and secure collaboration environment, aiming at creating a central hub for massive symptom-herb connections. FangNet can be accessed at http://fangnet.org or http://fangnet.herb.ac.cn.

VarBen: Generating in Silico Reference Data Sets for Clinical Next-Generation Sequencing Bioinformatics Pipeline Evaluation.

Next-generation sequencing is increasingly being adopted as a valuable method for the detection of somatic variants in clinical oncology. However, it is still challenging to reach a satisfactory level of robustness and standardization in clinical practice when using the currently available bioinformatics pipelines to detect variants from raw sequencing data. Moreover, appropriate reference data sets are lacking for clinical bioinformatics pipeline development, validation, and proficiency testing. Herein, we developed the Variant Benchmark tool (VarBen), an open-source software for variant simulation to generate customized reference data sets by directly editing the original sequencing reads. VarBen can introduce a variety of variants, including single-nucleotide variants, small insertions and deletions, and large structural variants, into targeted, exome, or whole-genome sequencing data, and can handle sequencing data from both the Illumina and Ion Torrent sequencing platforms. To demonstrate the feasibility and robustness of VarBen, we performed variant simulation on different sequencing data sets and compared the simulated variants with real-world data. The validation study showed that the simulated data are highly comparable to real-world data and that VarBen is a reliable tool for variant simulation. In addition, our collaborative study of somatic variant calling in 20 laboratories emphasizes the need for laboratories to evaluate their bioinformatics pipelines with customized reference data sets. VarBen may help users develop and validate their bioinformatics pipelines using locally generated sequencing data.

NONCODEV6: an updated database dedicated to long non-coding RNA annotation in both animals and plants.

NONCODE (http://www.noncode.org/) is a comprehensive database of collection and annotation of noncoding RNAs, especially long non-coding RNAs (lncRNAs) in animals. NONCODEV6 is dedicated to providing the full scope of lncRNAs across plants and animals. The number of lncRNAs in NONCODEV6 has increased from 548 640 to 644 510 since the last update in 2017. The number of human lncRNAs has increased from 172 216 to 173 112. The number of mouse lncRNAs increased from 131 697 to 131 974. The number of plant lncRNAs is 94 697. The relationship between lncRNAs in human and cancer were updated with transcriptome sequencing profiles. Three important new features were also introduced in NONCODEV6: (i) updated human lncRNA-disease relationships, especially cancer; (ii) lncRNA annotations with tissue expression profiles and predicted function in five common plants; iii) lncRNAs conservation annotation at transcript level for 23 plant species. NONCODEV6 is accessible through http://www.noncode.org/.

HERB: a high-throughput experiment- and reference-guided database of traditional Chinese medicine.

Pharmacotranscriptomics has become a powerful approach for evaluating the therapeutic efficacy of drugs and discovering new drug targets. Recently, studies of traditional Chinese medicine (TCM) have increasingly turned to high-throughput transcriptomic screens for molecular effects of herbs/ingredients. And numerous studies have examined gene targets for herbs/ingredients, and link herbs/ingredients to various modern diseases. However, there is currently no systematic database organizing these data for TCM. Therefore, we built HERB, a high-throughput experiment- and reference-guided database of TCM, with its Chinese name as BenCaoZuJian. We re-analyzed 6164 gene expression profiles from 1037 high-throughput experiments evaluating TCM herbs/ingredients, and generated connections between TCM herbs/ingredients and 2837 modern drugs by mapping the comprehensive pharmacotranscriptomics dataset in HERB to CMap, the largest such dataset for modern drugs. Moreover, we manually curated 1241 gene targets and 494 modern diseases for 473 herbs/ingredients from 1966 references published recently, and cross-referenced this novel information to databases containing such data for drugs. Together with database mining and statistical inference, we linked 12 933 targets and 28 212 diseases to 7263 herbs and 49 258 ingredients and provided six pairwise relationships among them in HERB. In summary, HERB will intensively support the modernization of TCM and guide rational modern drug discovery efforts. And it is accessible through http://herb.ac.cn/.