Numerous studies indicate that non-coding RNAs (ncRNAs) have critical functions across biological processes, and single-nucleotide polymorphisms (SNPs) could contribute to diseases or traits through influencing ncRNA expression. However, the associations between SNPs and ncRNA expression are largely unknown. Therefore, genome-wide expression quantitative trait loci (eQTL) analysis to assess the effects of SNPs on ncRNA expression, especially in multiple cancer types, will help to understand how risk alleles contribute toward tumorigenesis and cancer development. Using genotype data and expression profiles of ncRNAs of >8700 samples from The Cancer Genome Atlas (TCGA), we developed a computational pipeline to systematically identify ncRNA-related eQTLs (ncRNA-eQTLs) across 33 cancer types. We identified a total of 6 133 278 and 721 122 eQTL-ncRNA pairs in cis-eQTL and trans-eQTL analyses, respectively. Further survival analyses identified 8312 eQTLs associated with patient survival times. Furthermore, we linked ncRNA-eQTLs to genome-wide association study (GWAS) data and found 262 332 ncRNA-eQTLs overlapping with known disease- and trait-associated loci. Finally, a user-friendly database, ncRNA-eQTL (http://ibi.hzau.edu.cn/ncRNA-eQTL), was developed for free searching, browsing and downloading of all ncRNA-eQTLs. We anticipate that such an integrative and comprehensive resource will improve our understanding of the mechanistic basis of human complex phenotypic variation, especially for ncRNA- and cancer-related studies.

Drug combination is now a hot research topic in the pharmaceutical industry, but experiment-based methodologies are extremely costly in time and money. Many computational methods have been proposed to address these problems by starting from existing drug combinations. However, in most cases, only molecular structure information is included, which covers too limited a set of drug characteristics to efficiently screen drug combinations. Here, we integrated similarity-based multifeature drug data to improve the prediction accuracy by using the neighbor recommender method combined with ensemble learning algorithms. By conducting feature assessment analysis, we selected the most useful drug features and achieved 0.964 AUC in the ensemble models. The comparison results showed that the ensemble models outperform traditional machine learning algorithms such as suppeort vector machine (SVM), naïve Bayes (NB), and logistic regression (GLM). Furthermore, we predicted 7 candidate drug combinations for a specific drug, paclitaxel, and successfully verified that the two of the predicted combinations have promising effects.

Developing safe and effective therapeutic modalities remains a critical challenge for improving the prognosis of patients with colorectal cancer (CRC). In this regard, targeting epigenetic regulation in cancers has recently emerged as a promising therapeutic approach. Since several natural compounds have recently been shown to be important epigenetic modulators, we hypothesized that Ginseng might exert its anti-cancer activity by regulating DNA methylation alterations in CRC. In this study, a series of cell culture studies were conducted, followed by their interrogation in patient-derived 3D organoid models to evaluate Ginseng’s anti-cancer activity in CRC. Genome-wide methylation alterations were interrogated by undertaking MethylationEpic BeadChip microarrays. First, 50% inhibitory concentrations (IC₅₀) were determined by cell viability assays, and subsequent Ginseng treatment demonstrated a significant anti-cancer effect on clonogenicity and cellular migration in CRC cells. Treatment with Ginseng potentiated cellular apoptosis through regulation of apoptosis-related genes in CRC cells. Furthermore, Ginseng treatment downregulated the expression of DNA methyltransferases (DNMTs) and decreased the global DNA methylation levels in CRC cells. The genome-wide methylation profiling identified Ginseng-induced hypomethylation of transcriptionally silenced tumor suppressor genes. Finally, cell culture-based findings were successfully validated in patient-derived 3D organoids. In conclusion, we demonstrate that Ginseng exerts its anti-tumorigenic potential by regulating cellular apoptosis via the downregulation of DNMTs and reversing the methylation status of transcriptionally silenced genes in CRC.

Genetic disease genes are considered a promising source of drug targets. Most diseases are caused by more than one pathogenic factor; thus, it is reasonable to consider that chemical agents targeting multiple disease genes are more likely to have desired activities. This is supported by a comprehensive analysis on the relationships between agent activity/druggability and target genetic characteristics. The therapeutic potential of agents increases steadily with increasing number of targeted disease genes, and can be further enhanced by strengthened genetic links between targets and diseases. By using the multi-label classification models for genetics-based drug activity prediction, we provide universal tools for prioritizing drug candidates. All of the documented data and the machine-learning prediction service are available at SCG-Drug (http://zhanglab.hzau.edu.cn/scgdrug).

Head and neck squamous cell carcinoma (HNSC) is one of most common malignancies with high mortality worldwide. Importantly, the molecular heterogeneity of HNSC complicates the clinical diagnosis and treatment, leading to poor overall survival outcomes. To dissect the complex heterogeneity, recent studies have reported multiple molecular subtyping systems. For instance, HNSC can be subdivided to four distinct molecular subtypes: atypical, basal, classical, and mesenchymal, of which the mesenchymal subtype is characterized by upregulated epithelial-mesenchymal transition (EMT) and associated with poorer survival outcomes. Despite a wealth of studies into the complex molecular heterogeneity, the regulatory mechanism specific to this aggressive subtype remain largely unclear. Herein, we developed a network-based bioinformatics framework that integrates lncRNA and mRNA expression profiles to elucidate the subtype-specific regulatory mechanisms. Applying the framework to HNSC, we identified a clinically relevant lncRNA LNCOG as a key master regulator mediating EMT underlying the mesenchymal subtype. Five genes with strong prognostic values, namely ANXA5, ITGA5, CCBE1, P4HA2, and EPHX3, were predicted to be the putative targets of LNCOG and subsequently validated in other independent datasets. By integrative analysis of the miRNA expression profiles, we found that LNCOG may act as a ceRNA to sponge miR-148a-3p thereby upregulating ITGA5 to promote HNSC progression. Furthermore, our drug sensitivity analysis demonstrated that the five putative targets of LNCOG were also predictive of the sensitivities of multiple FDA-approved drugs. In summary, our bioinformatics framework facilitates the dissection of cancer subtype-specific lncRNA regulatory mechanisms, providing potential novel biomarkers for more optimized treatment of HNSC.