PLOS Computational Biology

scMSI: Accurately inferring the sub-clonal Micro-Satellite status by an integrated deconvolution model on length spectrum

Microsatellite instability (MSI) is an important genomic biomarker for cancer diagnosis and treatment, and sequencing-based approaches are often applied to identify MSI because of its fastness and efficiency. These approaches, however, may fail to identify MSI on one or more sub-clones for certain cancers with a high degree of heterogeneity, leading to erroneous diagnoses and unsuitable treatments. Besides, the computational cost of identifying sub-clonal MSI can be exponentially increased when multiple sub-clones with different length distributions share MSI status. Herein, this paper proposes “scMSI”, an accurate and efficient estimation of sub-clonal MSI to identify the microsatellite status. scMSI is an integrative Bayesian method to deconvolute the mixed-length distribution of sub-clones by a novel alternating iterative optimization procedure based on a subtle generative model. During the process of deconvolution, the optimized division of each sub-clone is attained by a heuristic algorithm, aligning with clone proportions that adhere optimally to the sample’s clonal structure. To evaluate the performance, 16 patients diagnosed with endometrial cancer, exhibiting positive responses to the treatment despite having negative MSI status based on sequencing-based approaches, were considered. Excitingly, scMSI reported MSI on sub-clones successfully, and the findings matched the conclusions on immunohistochemistry. In addition, testing results on a series of experiments with simulation datasets concerning a variety of impact factors demonstrated the effectiveness and superiority of scMSI in detecting MSI on sub-clones over existing approaches. scMSI provides a new way of detecting MSI for cancers with a high degree of heterogeneity.

HPVsim: An agent-based model of HPV transmission and cervical disease

In 2020, the WHO launched its first global strategy to accelerate the elimination of cervical cancer, outlining an ambitious set of targets for countries to achieve over the next decade. At the same time, new tools, technologies, and strategies are in the pipeline that may improve screening performance, expand the reach of prophylactic vaccines, and prevent the acquisition, persistence and progression of oncogenic HPV. Detailed mechanistic modelling can help identify the combinations of current and future strategies to combat cervical cancer. Open-source modelling tools are needed to shift the capacity for such evaluations in-country. Here, we introduce the Human papillomavirus simulator (HPVsim), a new open-source software package for creating flexible agent-based models parameterised with country-specific vital dynamics, structured sexual networks, and co-transmitting HPV genotypes. HPVsim includes a novel methodology for modelling cervical disease progression, designed to be readily adaptable to new forms of screening. The software itself is implemented in Python, has built-in tools for simulating commonly-used interventions, includes a comprehensive set of tests and documentation, and runs quickly (seconds to minutes) on a laptop. Performance is greatly enhanced by HPVsim’s multiscale modelling functionality. HPVsim is open source under the MIT License and available via both the Python Package Index (via pip install) and GitHub (hpvsim.org).

Mathematical modeling of the evolution of resistance and aggressiveness of high-grade serous ovarian cancer from patient CA-125 time series

A time-series analysis of serum Cancer Antigen 125 (CA-125) levels was performed in 791 patients with high-grade serous ovarian cancer (HGSOC) from the Australian Ovarian Cancer Study to evaluate the development of chemoresistance and response to therapy. To investigate chemoresistance and better predict the treatment effectiveness, we examined two traits: resistance (defined as the rate of CA-125 change when patients were treated with therapy) and aggressiveness (defined as the rate of CA-125 change when patients were not treated). We found that as the number of treatment lines increases, the data-based resistance increases (a decreased rate of CA-125 decay). We use mathematical models of two distinct cancer cell types, treatment-sensitive cells and treatment-resistant cells, to estimate the values and evolution of the two traits in individual patients. By fitting to individual patient HGSOC data, our models successfully capture the dynamics of the CA-125 level. The parameters estimated from the mathematical models show that patients with inferred low growth rates of treatment-sensitive cells and treatment-resistant cells (low model-estimated aggressiveness) and a high death rate of treatment-resistant cells (low model-estimated resistance) have longer survival time after completing their second-line of therapy. These findings show that mathematical models can characterize the degree of resistance and aggressiveness in individual patients, which improves our understanding of chemoresistance development and could predict treatment effectiveness in HGSOC patients.

Single-cell analysis reveals the pan-cancer invasiveness-associated transition of adipose-derived stromal cells into COL11A1-expressing cancer-associated fibroblasts

During the last ten years, many research results have been referring to a particular type of cancer-associated fibroblasts associated with poor prognosis, invasiveness, metastasis and resistance to therapy in multiple cancer types, characterized by a gene expression signature with prominent presence of genes COL11A1 , THBS2 and INHBA . Identifying the underlying biological mechanisms responsible for their creation may facilitate the discovery of targets for potential pan-cancer therapeutics. Using a novel computational approach for single-cell gene expression data analysis identifying the dominant cell populations in a sequence of samples from patients at various stages, we conclude that these fibroblasts are produced by a pan-cancer cellular transition originating from a particular type of adipose-derived stromal cells naturally present in the stromal vascular fraction of normal adipose tissue, having a characteristic gene expression signature. Focusing on a rich pancreatic cancer dataset, we provide a detailed description of the continuous modification of the gene expression profiles of cells as they transition from APOD -expressing adipose-derived stromal cells to COL11A1 -expressing cancer-associated fibroblasts, identifying the key genes that participate in this transition. These results also provide an explanation to the well-known fact that the adipose microenvironment contributes to cancer progression.

Ten simple rules for empowering women in STEM