Given the association of disturbances in non-esterified fatty acid (NEFA) metabolism with the development of Type 2 Diabetes and Non-Alcoholic Fatty Liver Disease, computational models of glucose-insulin dynamics have been extended to account for the interplay with NEFA. In this study, we use arteriovenous measurement across the subcutaneous adipose tissue during a mixed meal challenge test to evaluate the performance and underlying assumptions of three existing models of adipose tissue metabolism and construct a new, refined model of adipose tissue metabolism. Our model introduces new terms, explicitly accounting for the conversion of glucose to glyceraldehye-3-phosphate, the postprandial influx of glycerol into the adipose tissue, and several physiologically relevant delays in insulin signalling in order to better describe the measured adipose tissues fluxes. We then applied our refined model to human adipose tissue flux data collected before and after a diet intervention as part of the Yoyo study, to quantify the effects of caloric restriction on postprandial adipose tissue metabolism. Significant increases were observed in the model parameters describing the rate of uptake and release of both glycerol and NEFA. Additionally, decreases in the model's delay in insulin signalling parameters indicates there is an improvement in adipose tissue insulin sensitivity following caloric restriction.

MOTIVATION: Microscopy imaging is an essential tool for medical diagnosis and molecular biology. It is particularly useful for extracting information about disease states, tissue heterogeneity and cell specific parameters such as cell type or cell size from biological specimens. However, the information obtained from the images is likely to be subjected to sampling and observational bias with respect to the underlying cell size/type distributions. RESULTS: We present an algorithm, Estimate Tissue Cell Size/Type Distribution (EstiTiCS), for the adjustment of the underestimation of the number of small cells and the size of measured cells while accounting for the section thickness independent of the tissue type. We introduce the sources of bias under different tissue distributions and their effect on the measured values with simulation experiments. Furthermore, we demonstrate our method on histological sections of paraffin-embedded adipose tissue sample images from 57 people from a dietary intervention study. This data consists of measured cell size and its distribution over the dietary intervention period at four time points. Adjusting for the bias with EstiTiCS results in a closer fit to the true/expected adipocyte size distribution with earlier studies. Therefore, we conclude that our method is suitable as the final step in estimating the tissue wide cell type/size distribution from microscopy imaging pipeline. AVAILABILITY AND IMPLEMENTATION: Source code and its documentation are available at https://github.com/michaelLenz/EstiTiCS The whole pipeline of our method is implemented in R and makes use of the 'nloptr' package. Adipose tissue data used for this study are available on request. CONTACT: [email protected], [email protected].

MOTIVATION: Comparing time courses of gene expression with time courses of phenotypic data may provide new insights in cellular mechanisms. In this study, we compared the performance of five pattern recognition methods with respect to their ability to relate genes and phenotypic data: one classical method (k-means) and four methods especially developed for time series [Short Time-series Expression Miner (STEM), Linear Mixed Model mixtures, Dynamic Time Warping for -Omics and linear modeling with R/Bioconductor limma package]. The methods were evaluated using data available from toxicological studies that had the aim to relate gene expression with phenotypic endpoints (i.e. to develop biomarkers for adverse outcomes). Additionally, technical aspects (influence of noise, number of time points and number of replicates) were evaluated on simulated data. RESULTS: None of the methods outperforms the others in terms of biology. Linear modeling with limma is mostly influenced by noise. STEM is mostly influenced by the number of biological replicates in the dataset, whereas k-means and linear modeling with limma are mostly influenced by the number of time points. In most cases, the results of the methods complement each other. We therefore provide recommendations to integrate the five methods. AVAILABILITY: The Matlab code for the simulations performed in this research is available in the Supplementary Data (Word file). The microarray data analysed in this paper are available at ArrayExpress (E-TOXM-22 and E-TOXM-23) and Gene Expression Omnibus (GSE39291). The phenotypic data are available in the Supplementary Data (Excel file). Links to the pattern recognition tools compared in this paper are provided in the main text. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.