Abstract Summary We present MetaMarker, a pipeline for discovering metagenomic biomarkers from whole-metagenome sequencing samples. Different from existing methods, MetaMarker is based on a de novo approach that does not require mapping raw reads to a reference database. We applied MetaMarker on whole-metagenome sequencing of colorectal cancer (CRC) stool samples from France to discover CRC specific metagenomic biomarkers. We showed robustness of the discovered biomarkers by validating in independent samples from Hong Kong, Austria, Germany and Denmark. We further demonstrated these biomarkers could be used to build a machine learning classifier for CRC prediction. Availability and implementation MetaMarker is freely available at https://bitbucket.org/mkoohim/metamarker under GPLv3 license. Supplementary information Supplementary data are available at Bioinformatics online.
Glioblastoma Recurrence and the Role of O6-Methylguanine–DNA Methyltransferase Promoter Methylation
Tumor recurrence in glioblastoma multiforme (GBM) is often attributed to acquired resistance to the standard chemotherapeutic agent, temozolomide (TMZ). Promoter methylation of the DNA repair gene MGMT (O6-methylguanine–DNA methyltransferase) has been associated with sensitivity to TMZ, whereas increased expression of MGMT has been associated with TMZ resistance. Clinical studies have observed a downward shift in MGMT methylation percentage from primary to recurrent stage tumors; however, the evolutionary processes that drive this shift and more generally the emergence and growth of TMZ-resistant tumor subpopulations are still poorly understood. Here, we develop a mathematical model, parameterized using clinical and experimental data, to investigate the role of MGMT methylation in TMZ resistance during the standard treatment regimen for GBM—surgery, chemotherapy, and radiation. We first found that the observed downward shift in MGMT promoter methylation status between detection and recurrence cannot be explained solely by evolutionary selection. Next, our model suggests that TMZ has an inhibitory effect on maintenance methylation of MGMT after cell division. Finally, incorporating this inhibitory effect, we study the optimal number of TMZ doses per adjuvant cycle for patients with GBM with high and low levels of MGMT methylation at diagnosis.
Quantifying Uncertainty and Robustness in a Biomathematical Model–Based Patient-Specific Response Metric for Glioblastoma
Purpose Glioblastomas, lethal primary brain tumors, are known for their heterogeneity and invasiveness. A growing body of literature has been developed demonstrating the clinical relevance of a biomathematical model, the proliferation-invasion model, of glioblastoma growth. Of interest here is the development of a treatment response metric, days gained (DG). This metric is based on individual tumor kinetics estimated through segmented volumes of hyperintense regions on T1-weighted gadolinium-enhanced and T2-weighted magnetic resonance images. This metric was shown to be prognostic of time to progression. Furthermore, it was shown to be more prognostic of outcome than standard response metrics. Although promising, the original article did not account for uncertainty in the calculation of the DG metric, leaving the robustness of this cutoff in question. Methods We harnessed the Bayesian framework to consider the impact of two sources of uncertainty: (1) image acquisition and (2) interobserver error in image segmentation. We first used synthetic data to characterize what nonerror variants are influencing the final uncertainty in the DG metric. We then considered the original patient cohort to investigate clinical patterns of uncertainty and to determine how robust this metric is for predicting time to progression and overall survival. Results Our results indicate that the key clinical variants are the time between pretreatment images and the underlying tumor growth kinetics, matching our observations in the clinical cohort. Finally, we demonstrated that for this cohort, there was a continuous range of cutoffs between 94 and 105 for which the prediction of the time to progression was over 80% reliable. Conclusion Although additional validation must be performed, this work represents a key step in ascertaining the clinical utility of this metric.
Sex differences in GBM revealed by analysis of patient imaging, transcriptome, and survival data