The discovery of novel anticancer drugs is critical for the pharmaceutical research and development, and patient treatment. Repurposing existing drugs that may have unanticipated effects as potential candidates is one way to meet this important goal. Systematic investigation of efficient anticancer drugs could provide valuable insights into trends in the discovery of anticancer drugs, which may contribute to the systematic discovery of new anticancer drugs. In this study, we collected and analyzed 150 anticancer drugs approved by the US Food and Drug Administration (FDA). Based on drug mechanism of action, these agents are divided into two groups: 61 cytotoxic-based drugs and 89 target-based drugs. We found that in the recent years, the proportion of targeted agents tended to be increasing, and the targeted drugs tended to be delivered as signal drugs. For 89 target-based drugs, we collected 102 effect-mediating drug targets in the human genome and found that most targets located on the plasma membrane and most of them belonged to the enzyme, especially tyrosine kinase. From above 150 drugs, we built a drug-cancer network, which contained 183 nodes (150 drugs and 33 cancer types) and 248 drug-cancer associations. The network indicated that the cytotoxic drugs tended to be used to treat more cancer types than targeted drugs. From 89 targeted drugs, we built a cancer-drug-target network, which contained 214 nodes (23 cancer types, 89 drugs, and 102 targets) and 313 edges (118 drug-cancer associations and 195 drug-target associations). Starting from the network, we discovered 133 novel drug-cancer associations among 52 drugs and 16 cancer types by applying the common target-based approach. Most novel drug-cancer associations (116, 87%) are supported by at least one clinical trial study. In this study, we provided a comprehensive data source, including anticancer drugs and their targets and performed a detailed analysis in term of historical tendency and networks. Its application to identify novel drug-cancer associations demonstrated that the data collected in this study is promising to serve as a fundamental for anticancer drug repurposing and development.

MicroRNAs (miRNAs) are small endogenous non-coding functional RNAs that post-transcriptionally regulate gene expression. They play essential roles in nearly all biological processes including cell development and differentiation, DNA damage repair, cell death as well as intercellular communication. They are highly involved in cancer, acting as tumor suppressors and/or promoters to modulate cell proliferation, epithelial-mesenchymal transition and tumor invasion and metastasis. Recent studies have shown that more than half of miRNAs are located within protein-coding or non-coding genes. Intragenic miRNAs and their host genes either share the promoter or have independent transcription. Meanwhile, miRNAs work as partners or antagonists of their host genes by fine-tuning their target genes functionally associated with host genes. This review outlined the complicated relationship between intragenic miRNAs and host genes. Focusing on miRNAs known as oncogenes or tumor suppressors in specific cancer types, it studied co-expression relationships between these miRNAs and host genes in the cancer types using TCGA data sets, which validated previous findings and revealed common, tumor-specific and even subtype-specific patterns. These observations will help understand the function of intragenic miRNAs and further develop miRNA therapeutics in cancer.

The wide application of next-generation sequencing (NGS), mainly through whole genome, exome and transcriptome sequencing, provides a high-resolution and global view of the cancer genome. Coupled with powerful bioinformatics tools, NGS promises to revolutionize cancer research, diagnosis and therapy. In this paper, we review the recent advances in NGS-based cancer genomic research as well as clinical application, summarize the current integrative oncogenomic projects, resources and computational algorithms, and discuss the challenge and future directions in the research and clinical application of cancer genomic sequencing.

Recent studies have shown that disease-susceptibility variants frequently lie in cell-type-specific enhancer elements. To identify, interpret, and prioritize such risk variants, we must identify the enhancers active in disease-relevant cell types, their upstream transcription factor (TF) binding, and their downstream target genes. To address this need, we built HACER (http://bioinfo.vanderbilt.edu/AE/HACER/), an atlas of Human ACtive Enhancers to interpret Regulatory variants. The HACER atlas catalogues and annotates in-vivo transcribed cell-type-specific enhancers, as well as placing enhancers within transcriptional regulatory networks by integrating ENCODE TF ChIP-Seq and predicted/validated chromatin interaction data. We demonstrate the utility of HACER in (i) offering a mechanistic hypothesis to explain the association of SNP rs614367 with ER-positive breast cancer risk, (ii) exploring tumor-specific enhancers in selective MYC dysregulation and (iii) prioritizing/annotating non-coding regulatory regions targeting CCND1. HACER provides a valuable resource for studies of GWAS, non-coding variants, and enhancer-mediated regulation.

Adopting a systems approach, we devise a general workflow to define actionable subtypes in human cancers. Applied to small cell lung cancer (SCLC), the workflow identifies four subtypes based on global gene expression patterns and ontologies. Three correspond to known subtypes (SCLC-A, SCLC-N, and SCLC-Y), while the fourth is a previously undescribed ASCL1+ neuroendocrine variant (NEv2, or SCLC-A2). Tumor deconvolution with subtype gene signatures shows that all of the subtypes are detectable in varying proportions in human and mouse tumors. To understand how multiple stable subtypes can arise within a tumor, we infer a network of transcription factors and develop BooleaBayes, a minimally-constrained Boolean rule-fitting approach. In silico perturbations of the network identify master regulators and destabilizers of its attractors. Specific to NEv2, BooleaBayes predicts ELF3 and NR0B1 as master regulators of the subtype, and TCF3 as a master destabilizer. Since the four subtypes exhibit differential drug sensitivity, with NEv2 consistently least sensitive, these findings may lead to actionable therapeutic strategies that consider SCLC intratumoral heterogeneity. Our systems-level approach should generalize to other cancer types. Small-cell lung cancer (SCLC) is an extremely aggressive disease with poor prognosis. Despite significant advances in treatments of other cancer types, therapeutic strategies for SCLC have remained unchanged for decades. We hypothesize that distinct SCLC subtypes with differential drug sensitivities may be responsible for poor treatment outcomes. To this end, we applied a computational pipeline to identify and characterize SCLC subtypes. We found four subtypes, including one (termed “NEv2”) that had not previously been reported. Across a broad panel of drugs, we show that NEv2 is more resistant than other SCLC subtypes, suggesting that this subtype may be partly responsible for poor treatment outcomes. Importantly, we validate the existence of NEv2 cells in both human and mouse tumors. Reprogramming the identity of NEv2 cells into other subtypes may sensitize these cells to existing treatments. However, deciphering global mechanisms that regulate different subtypes is generally unfeasible. To circumvent this, we developed BooleaBayes, a modeling approach that only infers local regulatory mechanisms near stable cell subtypes. Using BooleaBayes, we found master regulators and master destabilizers for each subtype. These findings predict targets that may destabilize a particular subtype, including NEv2, and lead to successful therapy, by either knocking out master regulators or turning on master destabilizers.

The sarcomere is the contractile unit within cardiomyocytes driving heart muscle contraction. We sought to test the mechanisms regulating actin and myosin filament assembly during sarcomere formation. Therefore, we developed an assay using human cardiomyocytes to monitor sarcomere assembly. We report a population of muscle stress fibers, similar to actin arcs in non-muscle cells, which are essential sarcomere precursors. We show sarcomeric actin filaments arise directly from muscle stress fibers. This requires formins (e.g., FHOD3), non-muscle myosin IIA and non-muscle myosin IIB. Furthermore, we show short cardiac myosin II filaments grow to form ~1.5 μm long filaments that then ‘stitch’ together to form the stack of filaments at the core of the sarcomere (i.e., the A-band). A-band assembly is dependent on the proper organization of actin filaments and, as such, is also dependent on FHOD3 and myosin IIB. We use this experimental paradigm to present evidence for a unifying model of sarcomere assembly.