The metabolomic window into hepatobiliary disease

The emergent discipline of metabolomics has attracted considerable research effort in hepatology. Here we review the metabolomic data for non-alcoholic fatty liver disease (NAFLD), non-alcoholic steatohepatitis (NASH), cirrhosis, hepatocellular carcinoma (HCC), cholangiocarcinoma (CCA), alcoholic liver disease (ALD), hepatitis B and C, cholecystitis, cholestasis, liver transplantation, and acute hepatotoxicity in animal models. A metabolomic window has permitted a view into the changing biochemistry occurring in the transitional phases between a healthy liver and hepatocellular carcinoma or cholangiocarcinoma. Whether provoked by obesity and diabetes, alcohol use or oncogenic viruses, the liver develops a core metabolomic phenotype (CMP) that involves dysregulation of bile acid and phospholipid homeostasis. The CMP commences at the transition between the healthy liver (Phase 0) and NAFLD/NASH, ALD or viral hepatitis (Phase 1). This CMP is maintained in the presence or absence of cirrhosis (Phase 2) and whether or not either HCC or CCA (Phase 3) develops. Inflammatory signalling in the liver triggers the appearance of the CMP. Many other metabolomic markers distinguish between Phases 0, 1, 2 and 3. A metabolic remodelling in HCC has been described but metabolomic data from all four Phases demonstrate that the Warburg shift from mitochondrial respiration to cytosolic glycolysis foreshadows HCC and may occur as early as Phase 1. The metabolic remodelling also involves an upregulation of fatty acid β-oxidation, also beginning in Phase 1. The storage of triglycerides in fatty liver provides high energy-yielding substrates for Phases 2 and 3 of liver pathology. The metabolomic window into hepatobiliary disease sheds new light on the systems pathology of the liver.

Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

Functional error modeling for uncertainty quantification in hydrogeology

Approximate models (proxies) can be employed to reduce the computational costs of estimating uncertainty. The price to pay is that the approximations introduced by the proxy model can lead to a biased estimation. To avoid this problem and ensure a reliable uncertainty quantification, we propose to combine functional data analysis and machine learning to build error models that allow us to obtain an accurate prediction of the exact response without solving the exact model for all realizations. We build the relationship between proxy and exact model on a learning set of geostatistical realizations for which both exact and approximate solvers are run. Functional principal components analysis (FPCA) is used to investigate the variability in the two sets of curves and reduce the dimensionality of the problem while maximizing the retained information. Once obtained, the error model can be used to predict the exact response of any realization on the basis of the sole proxy response. This methodology is purpose-oriented as the error model is constructed directly for the quantity of interest, rather than for the state of the system. Also, the dimensionality reduction performed by FPCA allows a diagnostic of the quality of the error model to assess the informativeness of the learning set and the fidelity of the proxy to the exact model. The possibility of obtaining a prediction of the exact response for any newly generated realization suggests that the methodology can be effectively used beyond the context of uncertainty quantification, in particular for Bayesian inference and optimization.

Regular exercise decreases liver tumors development in hepatocyte-specific PTEN-deficient mice independently of steatosis.

BACKGROUND & AIMS Unhealthy lifestyles predispose to non-alcoholic steatohepatitis (NASH), which may further result in the development of hepatocellular carcinoma (HCC). Although NASH patients benefit from physical activity, it is unknown whether regular exercise reduces the risk of developing HCC. Therefore, we studied the effect of regular exercise on the development of HCC in male hepatocyte-specific PTEN-deficient mice (AlbCrePten(flox/flox)), which develop steatohepatitis and HCC spontaneously. METHODS Mice were fed a standardized 10% fat diet and were randomly divided into exercise or sedentary groups. The exercise group ran on a motorized treadmill for 60 minutes/day, 5 days/week during 32 weeks. RESULTS After 32 weeks of regular exercise, 71% of exercised mice developed nodules larger than 15 mm(3) vs 100% of mice in the sedentary group. The mean number of tumors per liver was reduced by exercise, as well as the total tumoral volume per liver. Exercise did not affect steatosis and had no effect on the Non-alcoholic fatty liver disease (NAFLD) Activity Score (NAS). Exercise decreased tumor cell proliferation. Mechanistically, exercise stimulated the phosphorylation of AMPK and its substrate raptor, which decreased the kinase activity of mTOR. CONCLUSIONS These data show a benefit of regular exercise on the development of HCC in an experimental model of NASH and offer a rationale for encouraging predisposed patients to increase their physical activity for the prevention of HCC.

Non-rigid Free-Form 2D-3D Registration Using Statistical Deformation Model

This paper presents a non-rigid free-from 2D-3D registration approach using statistical deformation model (SDM). In our approach the SDM is first constructed from a set of training data using a non-rigid registration algorithm based on b-spline free-form deformation to encode a priori information about the underlying anatomy. A novel intensity-based non-rigid 2D-3D registration algorithm is then presented to iteratively fit the 3D b-spline-based SDM to the 2D X-ray images of an unseen subject, which requires a computationally expensive inversion of the instantiated deformation in each iteration. In this paper, we propose to solve this challenge with a fast B-spline pseudo-inversion algorithm that is implemented on graphics processing unit (GPU). Experiments conducted on C-arm and X-ray images of cadaveric femurs demonstrate the efficacy of the present approach.

Assessment of the Hong Kong Liver Cancer Staging System in Europe.

BACKGROUND & AIMS European and American guidelines have endorsed the Barcelona Clinic Liver Cancer (BCLC) staging system. The aim of this study was to assess the performance of the recently developed Hong Kong Liver Cancer (HKLC) classification as a staging system for hepatocellular carcinoma (HCC) in Europe. METHODS We used a pooled set of 1693 HCC patients combining three prospective European cohorts. Discrimination ability between the nine substages and five stages of the HKLC classification system was assessed. To evaluate the predictive power of the HKLC and BCLC staging systems on overall survival, Nagelkerke pseudo R2, Bayesian Information Criterion and Harrell's concordance index were calculated. The number of patients who would benefit from a curative therapy was assessed for both staging system. RESULTS The HKLC classification in nine substages shows suboptimal discrimination between the staging groups. The classification in five stages shows better discrimination between groups. However, the BCLC classification performs better than the HKLC classification in the ability to predict OS. The HKLC treatment algorithm tags significantly more patients to curative therapy than the BCLC. CONCLUSIONS The BCLC staging system performs better for European patients than the HKLC staging system in predicting OS. Twice more patients are eligible for a curative therapy with the HKLC algorithm, whether this translates in survival benefit remains to be investigated. This article is protected by copyright. All rights reserved.

Global Optimization with Sparse and Local Gaussian Process Models

We present a novel surrogate model-based global optimization framework allowing a large number of function evaluations. The method, called SpLEGO, is based on a multi-scale expected improvement (EI) framework relying on both sparse and local Gaussian process (GP) models. First, a bi-objective approach relying on a global sparse GP model is used to determine potential next sampling regions. Local GP models are then constructed within each selected region. The method subsequently employs the standard expected improvement criterion to deal with the exploration-exploitation trade-off within selected local models, leading to a decision on where to perform the next function evaluation(s). The potential of our approach is demonstrated using the so-called Sparse Pseudo-input GP as a global model. The algorithm is tested on four benchmark problems, whose number of starting points ranges from 102 to 104. Our results show that SpLEGO is effective and capable of solving problems with large number of starting points, and it even provides significant advantages when compared with state-of-the-art EI algorithms.

Differentiating the Multipoint Expected Improvement for Optimal Batch Design

This work deals with parallel optimization of expensive objective functions which are modelled as sample realizations of Gaussian processes. The study is formalized as a Bayesian optimization problem, or continuous multi-armed bandit problem, where a batch of q > 0 arms is pulled in parallel at each iteration. Several algorithms have been developed for choosing batches by trading off exploitation and exploration. As of today, the maximum Expected Improvement (EI) and Upper Confidence Bound (UCB) selection rules appear as the most prominent approaches for batch selection. Here, we build upon recent work on the multipoint Expected Improvement criterion, for which an analytic expansion relying on Tallis’ formula was recently established. The computational burden of this selection rule being still an issue in application, we derive a closed-form expression for the gradient of the multipoint Expected Improvement, which aims at facilitating its maximization using gradient-based ascent algorithms. Substantial computational savings are shown in application. In addition, our algorithms are tested numerically and compared to state-of-the-art UCB-based batchsequential algorithms. Combining starting designs relying on UCB with gradient-based EI local optimization finally appears as a sound option for batch design in distributed Gaussian Process optimization.

Automatic assessment of time-resolved OCT images for selective retina therapy

Purpose In recent years, selective retina laser treatment (SRT), a sub-threshold therapy method, avoids widespread damage to all retinal layers by targeting only a few. While these methods facilitate faster healing, their lack of visual feedback during treatment represents a considerable shortcoming as induced lesions remain invisible with conventional imaging and make clinical use challenging. To overcome this, we present a new strategy to provide location-specific and contact-free automatic feedback of SRT laser applications. Methods We leverage time-resolved optical coherence tomography (OCT) to provide informative feedback to clinicians on outcomes of location-specific treatment. By coupling an OCT system to SRT treatment laser, we visualize structural changes in the retinal layers as they occur via time-resolved depth images. We then propose a novel strategy for automatic assessment of such time-resolved OCT images. To achieve this, we introduce novel image features for this task that when combined with standard machine learning classifiers yield excellent treatment outcome classification capabilities. Results Our approach was evaluated on both ex vivo porcine eyes and human patients in a clinical setting, yielding performances above 95 % accuracy for predicting patient treatment outcomes. In addition, we show that accurate outcomes for human patients can be estimated even when our method is trained using only ex vivo porcine data. Conclusion The proposed technique presents a much needed strategy toward noninvasive, safe, reliable, and repeatable SRT applications. These results are encouraging for the broader use of new treatment options for neovascularization-based retinal pathologies.

Experimental models for p53, hepatitis B and aflatoxin hepatocarcinogenesis

The major risk factors for liver cancer in Southeast Asia: HBV infection, aflatoxin exposure and p53 expression/mutation, were examined in experimental models. Four groups were examined for development of hepatocellular carcinoma (HCC) with and without neonatal exposure to aflatoxin (AFB$\sb1)$: (Group I.) Transgenic HBsAg mice with one p53 allele. (Group II) Transgenic HBsAg mice with two p53 alleles. (Group III) Non-transgenic litter mates with one p53 allele. (Group IV) Non-transgenic litter mates with two p53 alleles. HCC developed in Group I animals exposed to aflatoxin at an earlier time and were of a higher grade than those seen later in other groups. These results provide an explanation for as to why p53 is a target for deletion and/or mutation in human HCC especially when found in high risk areas where HBV infection and Aflatoxin B1 food contamination is high, and nicely illustrates a synergistic interaction among these three factors. None of the tumors analyzed had loss or mutation in the p53 gene.^ To determine the significance of the specific p53ser249 mutation found in HBV/aflatoxin associated human hepatomas in an in-vivo experimental model using transgenic mice, a two-nucleotide change in the mouse p53 gene at amino acid position 246, which is equivalent to that of 249 in human p53, was introduced. Transgenic mice with mutant p53 controlled by the albumin promoter were generated and shown to express the p53ser246 mutant RNA and protein specifically in liver. Three groups were examined for development of HCC with and without neonatal exposure to aflatoxin: (Group V) Transgenic p53ser246 mice with two p53 alleles. (Group VI) Transgenic p53ser246 mice with one p53 allele. (Group VII) Double transgenic for p53ser246 and HBsAg with two p53 alleles. One hundred percent of male mice with the three risk factors injected with aflatoxin developed high grade liver tumors, compared to 66.6% from group VI and only 14.2% of group V suggesting synergistic interaction between HBsAg and this particular ser246 p53 mutation.^ In order to examine the growth properties of hepatocytes and correlation with p53 loss and/or mutation, cell proliferation and ploidy analysis of liver from normal heterozyous, homozygous null mice and from transgenic mutant p53ser246, mice were studied. Loss of wild-type p53 increased G1/G0 ratios of cells as well as proliferation and decreased cell ploidy. The mutant p53ser246 did not show a significant effect on cell ploidy or proliferation. However a striking 5-10X increase in G1/G0 ratio suggests that this specific mutation specifically induces G0 to G1 transition, which in turn further predisposes hepatocytes to the damaging effect of Aflatoxin. (Abstract shortened by UMI.) ^

On Software Engineering Repositories and Their Open Problems

In the last decade, a large number of software repositories have been created for different purposes. In this paper we present a survey of the publicly available repositories and classify the most common ones as well as discussing the problems faced by researchers when applying machine learning or statistical techniques to them.

Methods and Tools for Objective Assessment of Psychomotor Skills in Laparoscopic Surgery

Training and assessment paradigms for laparoscopic surgical skills are evolving from traditional mentor–trainee tutorship towards structured, more objective and safer programs. Accreditation of surgeons requires reaching a consensus on metrics and tasks used to assess surgeons’ psychomotor skills. Ongoing development of tracking systems and software solutions has allowed for the expansion of novel training and assessment means in laparoscopy. The current challenge is to adapt and include these systems within training programs, and to exploit their possibilities for evaluation purposes. This paper describes the state of the art in research on measuring and assessing psychomotor laparoscopic skills. It gives an overview on tracking systems as well as on metrics and advanced statistical and machine learning techniques employed for evaluation purposes. The later ones have a potential to be used as an aid in deciding on the surgical competence level, which is an important aspect when accreditation of the surgeons in particular, and patient safety in general, are considered. The prospective of these methods and tools make them complementary means for surgical assessment of motor skills, especially in the early stages of training. Successful examples such as the Fundamentals of Laparoscopic Surgery should help drive a paradigm change to structured curricula based on objective parameters. These may improve the accreditation of new surgeons, as well as optimize their already overloaded training schedules.

Stress Detection by means of Stress Physiological Template

This paper describes a stress detection system based on fuzzy logic and two physiological signals: Galvanic Skin Response and Heart Rate. Instead of providing a global stress classification, this approach creates an individual stress templates, gathering the behaviour of individuals under situations with different degrees of stress. The proposed method is able to detect stress properly with a rate of 99.5%, being evaluated with a database of 80 individuals. This result improves former approaches in the literature and well-known machine learning techniques like SVM, k-NN, GMM and Linear Discriminant Analysis. Finally, the proposed method is highly suitable for real-time applications

Two way clustering of Microarray Data using a Hybrid Approach

The Microarray technique is rather powerful, as it allows to test up thousands of genes at a time, but this produces an overwhelming set of data files containing huge amounts of data, which is quite difficult to pre-process, separate, classify and correlate for interesting conclusions to be extracted. Modern machine learning, data mining and clustering techniques based on information theory, are needed to read and interpret the information contents buried in those large data sets. Independent Component Analysis method can be used to correct the data affected by corruption processes or to filter the uncorrectable one and then clustering methods can group similar genes or classify samples. In this paper a hybrid approach is used to obtain a two way unsupervised clustering for a corrected microarray data.