ACDC: Automated Cell Detection and Counting for Time-Lapse Fluorescence Microscopy.

Advances in microscopy imaging technologies have enabled the visualization of live-cell dynamic processes using time-lapse microscopy imaging. However, modern methods exhibit several limitations related to the training phases and to time constraints, hindering their application in the laboratory practice. In this work, we present a novel method, named Automated Cell Detection and Counting (ACDC), designed for activity detection of fluorescent labeled cell nuclei in time-lapse microscopy. ACDC overcomes the limitations of the literature methods, by first applying bilateral filtering on the original image to smooth the input cell images while preserving edge sharpness, and then by exploiting the watershed transform and morphological filtering. Moreover, ACDC represents a feasible solution for the laboratory practice, as it can leverage multi-core architectures in computer clusters to efficiently handle large-scale imaging datasets. Indeed, our Parent-Workers implementation of ACDC allows to obtain up to a 3.7× speed-up compared to the sequential counterpart. ACDC was tested on two distinct cell imaging datasets to assess its accuracy and effectiveness on images with different characteristics. We achieved an accurate cell-count and nuclei segmentation without relying on large-scale annotated datasets, a result confirmed by the average Dice Similarity Coefficients of 76.84 and 88.64 and the Pearson coefficients of 0.99 and 0.96, calculated against the manual cell counting, on the two tested datasets

Overcoming intratumoural heterogeneity for reproducible molecular risk stratification: a case study in advanced kidney cancer

Abstract Background Metastatic clear cell renal cell cancer (mccRCC) portends a poor prognosis and urgently requires better clinical tools for prognostication as well as for prediction of response to treatment. Considerable investment in molecular risk stratification has sought to overcome the performance ceiling encountered by methods restricted to traditional clinical parameters. However, replication of results has proven challenging, and intratumoural heterogeneity (ITH) may confound attempts at tissue-based stratification. Methods We investigated the influence of confounding ITH on the performance of a novel molecular prognostic model, enabled by pathologist-guided multiregion sampling (n = 183) of geographically separated mccRCC cohorts from the SuMR trial (development, n = 22) and the SCOTRRCC study (validation, n = 22). Tumour protein levels quantified by reverse phase protein array (RPPA) were investigated alongside clinical variables. Regularised wrapper selection identified features for Cox multivariate analysis with overall survival as the primary endpoint. Results The optimal subset of variables in the final stratification model consisted of N-cadherin, EPCAM, Age, mTOR (NEAT). Risk groups from NEAT had a markedly different prognosis in the validation cohort (log-rank p = 7.62 × 10−7; hazard ratio (HR) 37.9, 95% confidence interval 4.1–353.8) and 2-year survival rates (accuracy = 82%, Matthews correlation coefficient = 0.62). Comparisons with established clinico-pathological scores suggest favourable performance for NEAT (Net reclassification improvement 7.1% vs International Metastatic Database Consortium score, 25.4% vs Memorial Sloan Kettering Cancer Center score). Limitations include the relatively small cohorts and associated wide confidence intervals on predictive performance. Our multiregion sampling approach enabled investigation of NEAT validation when limiting the number of samples analysed per tumour, which significantly degraded performance. Indeed, sample selection could change risk group assignment for 64% of patients, and prognostication with one sample per patient performed only slightly better than random expectation (median logHR = 0.109). Low grade tissue was associated with 3.5-fold greater variation in predicted risk than high grade (p = 0.044). Conclusions This case study in mccRCC quantitatively demonstrates the critical importance of tumour sampling for the success of molecular biomarker studies research where ITH is a factor. The NEAT model shows promise for mccRCC prognostication and warrants follow-up in larger cohorts. Our work evidences actionable parameters to guide sample collection (tumour coverage, size, grade) to inform the development of reproducible molecular risk stratification methods

Functional Transcription Factor Target Networks Illuminate Control of Epithelial Remodelling

Cell identity is governed by gene expression, regulated by transcription factor (TF) binding at cis-regulatory modules. Decoding the relationship between TF binding patterns and gene regulation is nontrivial, remaining a fundamental limitation in understanding cell decision-making. We developed the NetNC software to predict functionally active regulation of TF targets; demonstrated on nine datasets for the TFs Snail, Twist, and modENCODE Highly Occupied Target (HOT) regions. Snail and Twist are canonical drivers of epithelial to mesenchymal transition (EMT), a cell programme important in development, tumour progression and fibrosis. Predicted “neutral” (non-functional) TF binding always accounted for the majority (50% to 95%) of candidate target genes from statistically significant peaks and HOT regions had higher functional binding than most of the Snail and Twist datasets examined. Our results illuminated conserved gene networks that control epithelial plasticity in development and disease. We identified new gene functions and network modules including crosstalk with notch signalling and regulation of chromatin organisation, evidencing networks that reshape Waddington’s epigenetic landscape during epithelial remodelling. Expression of orthologous functional TF targets discriminated breast cancer molecular subtypes and predicted novel tumour biology, with implications for precision medicine. Predicted invasion roles were validated using a tractable cell model, supporting our approach

Carbonic Anhydrase 9 Expression Increases with Vascular Endothelial Growth Factor-Targeted Therapy and Is Predictive of Outcome in Metastatic Clear Cell Renal Cancer

AbstractBackgroundThere is a lack of biomarkers to predict outcome with targeted therapy in metastatic clear cell renal cancer (mccRCC). This may be because dynamic molecular changes occur with therapy.ObjectiveTo explore if dynamic, targeted-therapy-driven molecular changes correlate with mccRCC outcome.Design, setting, and participantsMultiple frozen samples from primary tumours were taken from sunitinib-naïve (n=22) and sunitinib-treated mccRCC patients (n=23) for protein analysis. A cohort (n=86) of paired, untreated and sunitinib/pazopanib-treated mccRCC samples was used for validation. Array comparative genomic hybridisation (CGH) analysis and RNA interference (RNAi) was used to support the findings.InterventionThree cycles of sunitinib 50mg (4 wk on, 2 wk off).Outcome measurements and statistical analysisReverse phase protein arrays (training set) and immunofluorescence automated quantitative analysis (validation set) assessed protein expression.Results and limitationsDifferential expression between sunitinib-naïve and treated samples was seen in 30 of 55 proteins (p<0.05 for each). The proteins B-cell CLL/lymphoma 2 (BCL2), mutL homolog 1 (MLH1), carbonic anhydrase 9 (CA9), and mechanistic target of rapamycin (mTOR) (serine/threonine kinase) had both increased intratumoural variance and significant differential expression with therapy. The validation cohort confirmed increased CA9 expression with therapy. Multivariate analysis showed high CA9 expression after treatment was associated with longer survival (hazard ratio: 0.48; 95% confidence interval, 0.26–0.87; p=0.02). Array CGH profiles revealed sunitinib was associated with significant CA9 region loss. RNAi CA9 silencing in two cell lines inhibited the antiproliferative effects of sunitinib. Shortcomings of the study include selection of a specific protein for analysis, and the specific time points at which the treated tissue was analysed.ConclusionsCA9 levels increase with targeted therapy in mccRCC. Lower CA9 levels are associated with a poor prognosis and possible resistance, as indicated by the validation cohort.Patient summaryDrug treatment of advanced kidney cancer alters molecular markers of treatment resistance. Measuring carbonic anhydrase 9 levels may be helpful in determining which patients benefit from therapy

Network biology and machine learning approaches to metastasis and treatment response

Cancer causes 13% of human deaths worldwide, 90% of which involve metastasis. The reactivation of embryonic processes in epithelial cancers—and the epithelial-mesenchymal transition (EMT) in particular—results in increased cell motility and invasiveness, and is a known mechanism for initiating metastasis. The reverse process, the mesenchymal-epithelial transition (MET), is implicated in the process of cells colonising pre-metastatic niches. Understanding the relationships between EMT, MET and metastasis is therefore highly relevant to cancer research and treatment. Key challenges include deciphering the large, uncharted space of gene function, mapping the complex signalling networks involved and understanding how the EMT and MET programmes function in vivo within specific environments and disease contexts. Inference and analysis of small-scale networks from human tumour tissue samples, scored for protein expression, provides insight into pleiotropy, complex interactions and context-specific behaviour. Small sets of proteins (10–50, representative of key biological processes) are scored using quantitative antibody-based technologies (e.g. immunofluorescence) to give static expression values. A novel inference algorithm specifically for these data, Gabi, is presented, which produces signed, directed networks. On synthetic data, inferred networks often recapitulate the information flow between proteins in ground truth connectivity. Directionality predictions are highly accurate (90% correct) if the input network structure is itself accurate. The Gabi algorithm was applied to study multiple carcinomas (renal, breast, ovarian), providing novel insights into the relationships between EMT players and fundamental processes dysregulated in cancers (e.g. apoptosis, proliferation). Survival analysis on these cohorts shows further evidence for association of EMT with poor outcome. A patent-pending method is presented for stratifying response to sunitinib in metastatic renal cancer patients. The method is based on a proportional hazards model with predictive features selected automatically using regularisation (Bayesian information criterion). The final algorithm includes N-cadherin expression, a determinant of mesenchymal properties, and shows significant predictive power (p = 7.6x10-7, log-rank test). A separate method stratifies response to tamoxifen in estrogen-receptor positive, node-negative breast cancer patients using a cross-validated support vector machine (SVM). The algorithm was predictive on blind-test data (p = 4.92 x 10-6, log-rank test). Methods developed have been made available within a web application (TMA Navigator) and an R package (rTMA). TMA Navigator produces visual data summaries, networks and survival analysis for uploaded tissue microarray (TMA) scores. rTMA expands on TMA Navigator capabilities for advanced workflows within a programming environment

TMA Navigator: network inference, patient stratification and survival analysis with tissue microarray data

Tissue microarrays (TMAs) allow multiplexed analysis of tissue samples and are frequently used to estimate biomarker protein expression in tumour biopsies. TMA Navigator (www.tmanavigator.org) is an open access web application for analysis of TMA data and related information, accommodating categorical, semi-continuous and continuous expression scores. Non-biological variation, or batch effects, can hinder data analysis and may be mitigated using the ComBat algorithm, which is incorporated with enhancements for automated application to TMA data. Unsupervised grouping of samples (patients) is provided according to Gaussian mixture modelling of marker scores, with cardinality selected by Bayesian information criterion regularization. Kaplan-Meier survival analysis is available, including comparison of groups identified by mixture modelling using the Mantel-Cox log-rank test. TMA Navigator also supports network inference approaches useful for TMA datasets, which often constitute comparatively few markers. Tissue and cell-type specific networks derived from TMA expression data offer insights into the molecular logic underlying pathophenotypes, towards more effective and personalized medicine. Output is interactive, and results may be exported for use with external programs. Private anonymous access is available, and user accounts may be generated for easier data management.</p