Integrative computational epigenomics to build data-driven gene regulation hypotheses

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

Variational autoencoding of gene landscapes during mouse CNS development uncovers layered roles of Polycomb Repressor Complex 2

A prominent aspect of most, if not all, central nervous systems (CNSs) is that anterior regions (brain) are larger than posterior ones (spinal cord). Studies in Drosophila and mouse have revealed that Polycomb Repressor Complex 2 (PRC2), a protein complex responsible for applying key repressive histone modifications, acts by several mechanisms to promote anterior CNS expansion. However, it is unclear what the full spectrum of PRC2 action is during embryonic CNS development and how PRC2 intersects with the epigenetic landscape. We removed PRC2 function from the developing mouse CNS, by mutating the key gene Eed, and generated spatio-temporal transcriptomic data. To decode the role of PRC2, we developed a method that incorporates standard statistical analyses with probabilistic deep learning to integrate the transcriptomic response to PRC2 inactivation with epigenetic data. This multi-variate analysis corroborates the central involvement of PRC2 in anterior CNS expansion, and also identifies several unanticipated cohorts of genes, such as proliferation and immune response genes. Furthermore, the analysis reveals specific profiles of regulation via PRC2 upon these gene cohorts. These findings uncover a differential logic for the role of PRC2 upon functionally distinct gene cohorts that drive CNS anterior expansion. To support the analysis of emerging multi-modal datasets, we provide a novel bioinformatics package that integrates transcriptomic and epigenetic datasets to identify regulatory underpinnings of heterogeneous biological processes.

Machine Learning in Epigenomics: Insights into Cancer Biology and Medicine

The recent deluge of genome-wide technologies for the mapping of the epigenome and resulting data in cancer samples has provided the opportunity for gaining insights into and understanding the roles of epigenetic processes in cancer. However, the complexity, high-dimensionality, sparsity, and noise associated with these data pose challenges for extensive integrative analyses. Machine Learning (ML) algorithms are particularly suited for epigenomic data analyses due to their flexibility and ability to learn underlying hidden structures. We will discuss four overlapping but distinct major categories under ML: dimensionality reduction, unsupervised methods, supervised methods, and deep learning (DL). We review the preferred use cases of these algorithms in analyses of cancer epigenomics data with the hope to provide an overview of how ML approaches can be used to explore fundamental questions on the roles of epigenome in cancer biology and medicine.

Unveiling the Pathogenesis of Psychiatric Disorders Using Network Models

Psychiatric disorders are complex brain disorders with a high degree of genetic heterogeneity, affecting millions of people worldwide. Despite advances in psychiatric genetics, the underlying pathogenic mechanisms of psychiatric disorders are still largely elusive, which impedes the development of novel rational therapies. There has been accumulating evidence suggesting that the genetics of complex disorders can be viewed through an omnigenic lens, which involves contextualizing genes in highly interconnected networks. Thus, applying network-based multi-omics integration methods could cast new light on the pathophysiology of psychiatric disorders. In this review, we first provide an overview of the recent advances in psychiatric genetics and highlight gaps in translating molecular associations into mechanistic insights. We then present an overview of network methodologies and review previous applications of network methods in the study of psychiatric disorders. Lastly, we describe the potential of such methodologies within a multi-tissue, multi-omics approach, and summarize the future directions in adopting diverse network approaches.

A Survey of Current Resources to Study lncRNA-Protein Interactions

Phenotypes are driven by regulated gene expression, which in turn are mediated by complex interactions between diverse biological molecules. Protein-DNA interactions such as histone and transcription factor binding are well studied, along with RNA-RNA interactions in short RNA silencing of genes. In contrast, lncRNA-protein interaction (LPI) mechanisms are comparatively unknown, likely directed by the difficulties in studying LPI. However, LPI are emerging as key interactions in epigenetic mechanisms, playing a role in development and disease. Their importance is further highlighted by their conservation across kingdoms. Hence, interest in LPI research is increasing. We therefore review the current state of the art in lncRNA-protein interactions. We specifically surveyed recent computational methods and databases which researchers can exploit for LPI investigation. We discovered that algorithm development is heavily reliant on a few generic databases containing curated LPI information. Additionally, these databases house information at gene-level as opposed to transcript-level annotations. We show that early methods predict LPI using molecular docking, have limited scope and are slow, creating a data processing bottleneck. Recently, machine learning has become the strategy of choice in LPI prediction, likely due to the rapid growth in machine learning infrastructure and expertise. While many of these methods have notable limitations, machine learning is expected to be the basis of modern LPI prediction algorithms.

A multi-modal data harmonisation approach for discovery of COVID-19 drug targets

Despite the volume of experiments performed and data available, the complex biology of coronavirus SARS-COV-2 is not yet fully understood. Existing molecular profiling studies have focused on analysing functional omics data of a single type, which captures changes in a small subset of the molecular perturbations caused by the virus. As the logical next step, results from multiple such omics analysis may be aggregated to comprehensively interpret the molecular mechanisms of SARS-CoV-2. An alternative approach is to integrate data simultaneously in a parallel fashion to highlight the inter-relationships of disease-driving biomolecules, in contrast to comparing processed information from each omics level separately. We demonstrate that valuable information may be masked by using the former fragmented views in analysis, and biomarkers resulting from such an approach cannot provide a systematic understanding of the disease aetiology. Hence, we present a generic, reproducible and flexible open-access data harmonisation framework that can be scaled out to future multi-omics analysis to study a phenotype in a holistic manner. The pipeline source code, detailed documentation and automated version as a R package are accessible. To demonstrate the effectiveness of our pipeline, we applied it to a drug screening task. We integrated multi-omics data to find the lowest level of statistical associations between data features in two case studies. Strongly correlated features within each of these two datasets were used for drug-target analysis, resulting in a list of 84 drug-target candidates. Further computational docking and toxicity analyses revealed seven high-confidence targets, amsacrine, bosutinib, ceritinib, crizotinib, nintedanib and sunitinib as potential starting points for drug therapy and development.

Circular RNA circCUL3 Accelerates the Warburg Effect Progression of Gastric Cancer through Regulating the STAT3/HK2 Axis

The Warburg effect is a significant hallmark of gastric cancer (GC), and increasing evidence emphasizes the crucial role of circular RNAs (circRNAs) in GC tumorigenesis. However, the precise molecular mechanisms by which circRNAs drive the GC Warburg effect are still elusive. The present study was designed to unveil the roles of circRNAs and the corresponding potential mechanism. High-regulated expression of circCUL3 was observed in both GC tissues and cell lines. Clinically, the high expression of circCUL3 was closely correlated with advanced clinical stage and overall survival in GC patients. Functionally, cellular experimental investigations demonstrated that circCUL3 promoted the proliferation, glucose consumption, lactate production, ATP quantity, and extracellular acidification rate (ECAR) of GC cells. In vivo, circCUL3 knockdown repressed tumor growth. Mechanistic analysis demonstrated that circCUL3 promoted signal transducer and activator of transcription (STAT)3 expression through sponging miR-515-5p; moreover, transcription factor STAT3 accelerated the transcriptional level of hexokinase 2 (HK2). In summary, the present findings provide mechanistic insights into circCUL3/miR-515-5p/STAT3/HK2 axis regulation on the GC Warburg effect, providing a novel possibility for an understanding of GC pathogenesis.