Benchmarking imputation methods for network inference using a novel method of synthetic scRNA-seq data generation

Background

Single cell RNA-sequencing (scRNA-seq) has very rapidly become the new workhorse of modern biology providing an unprecedented global view on cellular diversity and heterogeneity. In particular, the structure of gene-gene expression correlation contains information on the underlying gene regulatory networks. However, interpretation of scRNA-seq data is challenging due to specific experimental error and biases that are unique to this kind of data including drop-out (or technical zeros).

Methods

To deal with this problem several methods for imputation of zeros for scRNA-seq have been developed. However, it is not clear how these processing steps affect inference of genetic networks from single cell data. Here, we introduce Biomodelling.jl, a tool for generation of synthetic scRNA-seq data using multiscale modelling of stochastic gene regulatory networks in growing and dividing cells.

Results

Our tool produces realistic transcription data with a known ground truth network topology that can be used to benchmark different approaches for gene regulatory network inference. Using this tool we investigate the impact of different imputation methods on the performance of several network inference algorithms.

Conclusions

Biomodelling.jl provides a versatile and useful tool for future development and benchmarking of network inference approaches using scRNA-seq data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04778-9

The influence of methylation status on a stochastic model of MGMT dynamics in glioblastoma: Phenotypic selection can occur with and without a downshift in promoter methylation status

Glioblastoma originates in the brain and is one of the most aggressive cancer types. Glioblastoma represents 15% of all brain tumours, with a median survival of 15 months. Although the current standard of care for such a tumour (the Stupp protocol) has shown positive results for the prognosis of patients, O-6-methylguanine-DNA methyltransferase (MGMT) driven drug resistance has been an issue of increasing concern and hence requires innovative approaches. In addition to the well established drug resistance factors such as tumour location and blood brain barriers, it is also important to understand how the genetic and epigenetic dynamics of the glioblastoma cells can play a role. One important aspect of this is the study of methylation status of MGMT following administration of temozolomide. In this paper, we extend our previously published model that simulated MGMT expression in glioblastoma cells to incorporate the promoter methylation status of MGMT. This methylation status has clinical significance and is used as a marker for patient outcomes. Using this model, we investigate the causative relationship between temozolomide treatment and the methylation status of the MGMT promoter in a population of cells. In addition by constraining the model to relevant biological data using Approximate Bayesian Computation, we were able to identify parameter regimes that yield different possible modes of resistances, namely, phenotypic selection of MGMT, a downshift in the methylation status of the MGMT promoter or both simultaneously. We analysed each of the parameter sets associated with the different modes of resistance, presenting representative solutions as well as discovering some similarities between them as well as unique requirements for each of them. Finally, we used them to devise optimal strategies for inhibiting MGMT expression with the aim of minimising live glioblastoma cell numbers.

Phenotypic selection through cell death: stochastic modelling of O-6-methylguanine-DNA methyltransferase dynamics

Glioblastoma (GBM) is the most aggressive malignant primary brain tumour with a median overall survival of 15 months. To treat GBM, patients currently undergo a surgical resection followed by exposure to radiotherapy and concurrent and adjuvant temozolomide (TMZ) chemotherapy. However, this protocol often leads to treatment failure, with drug resistance being the main reason behind this. To date, many studies highlight the role of O-6-methylguanine-DNA methyltransferase (MGMT) in conferring drug resistance. The mechanism through which MGMT confers resistance is not well studied-particularly in terms of computational models. With only a few reasonable biological assumptions, we were able to show that even a minimal model of MGMT expression could robustly explain TMZ-mediated drug resistance. In particular, we showed that for a wide range of parameter values constrained by novel cell growth and viability assays, a model accounting for only stochastic gene expression of MGMT coupled with cell growth, division, partitioning and death was able to exhibit phenotypic selection of GBM cells expressing MGMT in response to TMZ. Furthermore, we found this selection allowed the cells to pass their acquired phenotypic resistance onto daughter cells in a stable manner (as long as TMZ is provided). This suggests that stochastic gene expression alone is enough to explain the development of chemotherapeutic resistance.