Recent Advances in Bioimage Analysis Methods for Detecting Skeletal Deformities in Biomedical and Aquaculture Fish Species.

Detecting skeletal or bone-related deformities in model and aquaculture fish is vital for numerous biomedical studies. In biomedical research, model fish with bone-related disorders are potential indicators of various chemically induced toxins in their environment or poor dietary conditions. In aquaculture, skeletal deformities are affecting fish health, and economic losses are incurred by fish farmers. This survey paper focuses on showcasing the cutting-edge image analysis tools and techniques based on artificial intelligence that are currently applied in the analysis of bone-related deformities in aquaculture and model fish. These methods and tools play a significant role in improving research by automating various aspects of the analysis. This paper also sheds light on some of the hurdles faced when dealing with high-content bioimages and explores potential solutions to overcome these challenges.

Multi-Task Pre-Training of Deep Neural Networks for Digital Pathology.

In this work, we investigate multi-task learning as a way of pre-training models for classification tasks in digital pathology. It is motivated by the fact that many small and medium-size datasets have been released by the community over the years whereas there is no large scale dataset similar to ImageNet in the domain. We first assemble and transform many digital pathology datasets into a pool of 22 classification tasks and almost 900k images. Then, we propose a simple architecture and training scheme for creating a transferable model and a robust evaluation and selection protocol in order to evaluate our method. Depending on the target task, we show that our models used as feature extractors either improve significantly over ImageNet pre-trained models or provide comparable performance. Fine-tuning improves performance over feature extraction and is able to recover the lack of specificity of ImageNet features, as both pre-training sources yield comparable performance.

Error curves for evaluating the quality of feature rankings.

In this article, we propose a method for evaluating feature ranking algorithms. A feature ranking algorithm estimates the importance of descriptive features when predicting the target variable, and the proposed method evaluates the correctness of these importance values by computing the error measures of two chains of predictive models. The models in the first chain are built on nested sets of top-ranked features, while the models in the other chain are built on nested sets of bottom ranked features. We investigate which predictive models are appropriate for building these chains, showing empirically that the proposed method gives meaningful results and can detect differences in feature ranking quality. This is first demonstrated on synthetic data, and then on several real-world classification benchmark problems.

A Random Forests Framework for Modeling Haplotypes as Mosaics of Reference Haplotypes.

Many genomic data analyses such as phasing, genotype imputation, or local ancestry inference share a common core task: matching pairs of haplotypes at any position along the chromosome, thereby inferring a target haplotype as a succession of pieces from reference haplotypes, commonly called a mosaic of reference haplotypes. For that purpose, these analyses combine information provided by linkage disequilibrium, linkage and/or genealogy through a set of heuristic rules or, most often, by a hidden Markov model. Here, we develop an extremely randomized trees framework to address the issue of local haplotype matching. In our approach, a supervised classifier using extra-trees (a particular type of random forests) learns how to identify the best local matches between haplotypes using a collection of observed examples. For each example, various features related to the different sources of information are observed, such as the length of a segment shared between haplotypes, or estimates of relationships between individuals, gametes, and haplotypes. The random forests framework was fed with 30 relevant features for local haplotype matching. Repeated cross-validations allowed ranking these features in regard to their importance for local haplotype matching. The distance to the edge of a segment shared by both haplotypes being matched was found to be the most important feature. Similarity comparisons between predicted and true whole-genome sequence haplotypes showed that the random forests framework was more efficient than a hidden Markov model in reconstructing a target haplotype as a mosaic of reference haplotypes. To further evaluate its efficiency, the random forests framework was applied to imputation of whole-genome sequence from 50k genotypes and it yielded average reliabilities similar or slightly better than IMPUTE2. Through this exploratory study, we lay the foundations of a new framework to automatically learn local haplotype matching and we show that extra-trees are a promising approach for such purposes. The use of this new technique also reveals some useful lessons on the relevant features for the purpose of haplotype matching. We also discuss potential improvements for routine implementation.

Characterization of a temporoparietal junction subtype of Alzheimer's disease.

Alzheimer's disease (AD) subtypes have been described according to genetics, neuropsychology, neuropathology, and neuroimaging. Thirty-one patients with clinically probable AD were selected based on perisylvian metabolic decrease on FDG-PET. They were compared to 25 patients with a typical pattern of decreased posterior metabolism. Tree-based machine learning was used on those 56 images to create a classifier that was subsequently applied to 207 Alzheimer's Disease Neuroimaging Initiative (ADNI) patients with AD. Machine learning was also used to discriminate between the two ADNI groups based on neuropsychological scores. Compared to AD patients with a typical precuneus metabolic decrease, the new subtype showed stronger hypometabolism in the temporoparietal junction. The classifier was able to distinguish the two groups in the ADNI population. Both groups could only be distinguished cognitively by Trail Making Test-A scores. This study further confirms that there is more than a typical metabolic pattern in probable AD with amnestic presentation.

Unsupervised Gene Network Inference with Decision Trees and Random Forests.

In this chapter, we introduce the reader to a popular family of machine learning algorithms, called decision trees. We then review several approaches based on decision trees that have been developed for the inference of gene regulatory networks (GRNs). Decision trees have indeed several nice properties that make them well-suited for tackling this problem: they are able to detect multivariate interacting effects between variables, are non-parametric, have good scalability, and have very few parameters. In particular, we describe in detail the GENIE3 algorithm, a state-of-the-art method for GRN inference.

Random Forests Based Group Importance Scores and Their Statistical Interpretation: Application for Alzheimer's Disease.

Machine learning approaches have been increasingly used in the neuroimaging field for the design of computer-aided diagnosis systems. In this paper, we focus on the ability of these methods to provide interpretable information about the brain regions that are the most informative about the disease or condition of interest. In particular, we investigate the benefit of group-based, instead of voxel-based, analyses in the context of Random Forests. Assuming a prior division of the voxels into non overlapping groups (defined by an atlas), we propose several procedures to derive group importances from individual voxel importances derived from Random Forests models. We then adapt several permutation schemes to turn group importance scores into more interpretable statistical scores that allow to determine the truly relevant groups in the importance rankings. The good behaviour of these methods is first assessed on artificial datasets. Then, they are applied on our own dataset of FDG-PET scans to identify the brain regions involved in the prognosis of Alzheimer's disease.

dynGENIE3: dynamical GENIE3 for the inference of gene networks from time series expression data.

The elucidation of gene regulatory networks is one of the major challenges of systems biology. Measurements about genes that are exploited by network inference methods are typically available either in the form of steady-state expression vectors or time series expression data. In our previous work, we proposed the GENIE3 method that exploits variable importance scores derived from Random forests to identify the regulators of each target gene. This method provided state-of-the-art performance on several benchmark datasets, but it could however not specifically be applied to time series expression data. We propose here an adaptation of the GENIE3 method, called dynamical GENIE3 (dynGENIE3), for handling both time series and steady-state expression data. The proposed method is evaluated extensively on the artificial DREAM4 benchmarks and on three real time series expression datasets. Although dynGENIE3 does not systematically yield the best performance on each and every network, it is competitive with diverse methods from the literature, while preserving the main advantages of GENIE3 in terms of scalability.

Landmark detection in 2D bioimages for geometric morphometrics: a multi-resolution tree-based approach.

The detection of anatomical landmarks in bioimages is a necessary but tedious step for geometric morphometrics studies in many research domains. We propose variants of a multi-resolution tree-based approach to speed-up the detection of landmarks in bioimages. We extensively evaluate our method variants on three different datasets (cephalometric, zebrafish, and drosophila images). We identify the key method parameters (notably the multi-resolution) and report results with respect to human ground truths and existing methods. Our method achieves recognition performances competitive with current existing approaches while being generic and fast. The algorithms are integrated in the open-source Cytomine software and we provide parameter configuration guidelines so that they can be easily exploited by end-users. Finally, datasets are readily available through a Cytomine server to foster future research.

SCENIC: single-cell regulatory network inference and clustering.

We present SCENIC, a computational method for simultaneous gene regulatory network reconstruction and cell-state identification from single-cell RNA-seq data (http://scenic.aertslab.org). On a compendium of single-cell data from tumors and brain, we demonstrate that cis-regulatory analysis can be exploited to guide the identification of transcription factors and cell states. SCENIC provides critical biological insights into the mechanisms driving cellular heterogeneity.

Sepsis prediction in critically ill patients by platelet activation markers on ICU admission: a prospective pilot study.

BACKGROUND: Platelets have been involved in both immune surveillance and host defense against severe infection. To date, whether platelet phenotype or other hemostasis components could be associated with predisposition to sepsis in critical illness remains unknown. The aim of this work was to identify platelet markers that could predict sepsis occurrence in critically ill injured patients. METHODS: This single-center, prospective, observational, 7-month study was based on a cohort of 99 non-infected adult patients admitted to ICUs for elective cardiac surgery, trauma, acute brain injury, and post-operative prolonged ventilation and followed up during ICU stay. Clinical characteristics and severity score (SOFA) were recorded on admission. Platelet activation markers, including fibrinogen binding to platelets, platelet membrane P-selectin expression, plasma soluble CD40L, and platelet-leukocytes aggregates were assayed by flow cytometry at admission and 48 h later, and then at the time of sepsis diagnosis (Sepsis-3 criteria) and 7 days later for sepsis patients. Hospitalization data and outcomes were also recorded. METHODS: Of the 99 patients, 19 developed sepsis after a median time of 5 days. These patients had a higher SOFA score at admission; levels of fibrinogen binding to platelets (platelet-Fg) and of D-dimers were also significantly increased compared to the other patients. Levels 48 h after ICU admission no longer differed between the two patient groups. Platelet-Fg % was an independent predictor of sepsis (P = 0.0031). By ROC curve analysis, cutoff point for Platelet-Fg (AUC = 0.75) was 50%. In patients with a SOFA cutoff of 8, the risk of sepsis reached 87% when Platelet-Fg levels were above 50%. Patients with sepsis had longer ICU and hospital stays and higher death rate. CONCLUSIONS: Platelet-bound fibrinogen levels assayed by flow cytometry within 24 h of ICU admission help identifying critically ill patients at risk of developing sepsis.

Circulating microRNA-based screening tool for breast cancer.

Circulating microRNAs (miRNAs) are increasingly recognized as powerful biomarkers in several pathologies, including breast cancer. Here, their plasmatic levels were measured to be used as an alternative screening procedure to mammography for breast cancer diagnosis.A plasma miRNA profile was determined by RT-qPCR in a cohort of 378 women. A diagnostic model was designed based on the expression of 8 miRNAs measured first in a profiling cohort composed of 41 primary breast cancers and 45 controls, and further validated in diverse cohorts composed of 108 primary breast cancers, 88 controls, 35 breast cancers in remission, 31 metastatic breast cancers and 30 gynecologic tumors.A receiver operating characteristic curve derived from the 8-miRNA random forest based diagnostic tool exhibited an area under the curve of 0.81. The accuracy of the diagnostic tool remained unchanged considering age and tumor stage. The miRNA signature correctly identified patients with metastatic breast cancer. The use of the classification model on cohorts of patients with breast cancers in remission and with gynecologic cancers yielded prediction distributions similar to that of the control group.Using a multivariate supervised learning method and a set of 8 circulating miRNAs, we designed an accurate, minimally invasive screening tool for breast cancer.

Collaborative analysis of multi-gigapixel imaging data using Cytomine.

MOTIVATION: Collaborative analysis of massive imaging datasets is essential to enable scientific discoveries. RESULTS: We developed Cytomine to foster active and distributed collaboration of multidisciplinary teams for large-scale image-based studies. It uses web development methodologies and machine learning in order to readily organize, explore, share and analyze (semantically and quantitatively) multi-gigapixel imaging data over the internet. We illustrate how it has been used in several biomedical applications. AVAILABILITY AND IMPLEMENTATION: Cytomine (http://www.cytomine.be/) is freely available under an open-source license from http://github.com/cytomine/ A documentation wiki (http://doc.cytomine.be) and a demo server (http://demo.cytomine.be) are also available. CONTACT: info@cytomine.be SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cerebral functional connectivity periodically (de)synchronizes with anatomical constraints.

This paper studies the link between resting-state functional connectivity (FC), measured by the correlations of fMRI BOLD time courses, and structural connectivity (SC), estimated through fiber tractography. Instead of a static analysis based on the correlation between SC and FC averaged over the entire fMRI time series, we propose a dynamic analysis, based on the time evolution of the correlation between SC and a suitably windowed FC. Assessing the statistical significance of the time series against random phase permutations, our data show a pronounced peak of significance for time window widths around 20-30 TR (40-60 s). Using the appropriate window width, we show that FC patterns oscillate between phases of high modularity, primarily shaped by anatomy, and phases of low modularity, primarily shaped by inter-network connectivity. Building upon recent results in dynamic FC, this emphasizes the potential role of SC as a transitory architecture between different highly connected resting-state FC patterns. Finally, we show that the regions contributing the most to these whole-brain level fluctuations of FC on the supporting anatomical architecture belong to the default mode and the executive control networks suggesting that they could be capturing consciousness-related processes such as mind wandering.

Classifying pairs with trees for supervised biological network inference.

Networks are ubiquitous in biology, and computational approaches have been largely investigated for their inference. In particular, supervised machine learning methods can be used to complete a partially known network by integrating various measurements. Two main supervised frameworks have been proposed: the local approach, which trains a separate model for each network node, and the global approach, which trains a single model over pairs of nodes. Here, we systematically investigate, theoretically and empirically, the exploitation of tree-based ensemble methods in the context of these two approaches for biological network inference. We first formalize the problem of network inference as a classification of pairs, unifying in the process homogeneous and bipartite graphs and discussing two main sampling schemes. We then present the global and the local approaches, extending the latter for the prediction of interactions between two unseen network nodes, and discuss their specializations to tree-based ensemble methods, highlighting their interpretability and drawing links with clustering techniques. Extensive computational experiments are carried out with these methods on various biological networks that clearly highlight that these methods are competitive with existing methods.

Evaluation and Comparison of Anatomical Landmark Detection Methods for Cephalometric X-Ray Images: A Grand Challenge.

Cephalometric analysis is an essential clinical and research tool in orthodontics for the orthodontic analysis and treatment planning. This paper presents the evaluation of the methods submitted to the Automatic Cephalometric X-Ray Landmark Detection Challenge, held at the IEEE International Symposium on Biomedical Imaging 2014 with an on-site competition. The challenge was set to explore and compare automatic landmark detection methods in application to cephalometric X-ray images. Methods were evaluated on a common database including cephalograms of 300 patients aged six to 60 years, collected from the Dental Department, Tri-Service General Hospital, Taiwan, and manually marked anatomical landmarks as the ground truth data, generated by two experienced medical doctors. Quantitative evaluation was performed to compare the results of a representative selection of current methods submitted to the challenge. Experimental results show that three methods are able to achieve detection rates greater than 80% using the 4 mm precision range, but only one method achieves a detection rate greater than 70% using the 2 mm precision range, which is the acceptable precision range in clinical practice. The study provides insights into the performance of different landmark detection approaches under real-world conditions and highlights achievements and limitations of current image analysis techniques.

A membrane-type-1 matrix metalloproteinase (MT1-MMP)-discoidin domain receptor 1 axis regulates collagen-induced apoptosis in breast cancer cells.

During tumour dissemination, invading breast carcinoma cells become confronted with a reactive stroma, a type I collagen-rich environment endowed with anti-proliferative and pro-apoptotic properties. To develop metastatic capabilities, tumour cells must acquire the capacity to cope with this novel microenvironment. How cells interact with and respond to their microenvironment during cancer dissemination remains poorly understood. To address the impact of type I collagen on the fate of tumour cells, human breast carcinoma MCF-7 cells were cultured within three-dimensional type I collagen gels (3D COL1). Using this experimental model, we have previously demonstrated that membrane type-1 matrix metalloproteinase (MT1-MMP), a proteinase overexpressed in many aggressive tumours, promotes tumour progression by circumventing the collagen-induced up-regulation of BIK, a pro-apoptotic tumour suppressor, and hence apoptosis. Here we performed a transcriptomic analysis to decipher the molecular mechanisms regulating 3D COL1-induced apoptosis in human breast cancer cells. Control and MT1-MMP expressing MCF-7 cells were cultured on two-dimensional plastic plates or within 3D COL1 and a global transcriptional time-course analysis was performed. Shifting the cells from plastic plates to 3D COL1 activated a complex reprogramming of genes implicated in various biological processes. Bioinformatic analysis revealed a 3D COL1-mediated alteration of key cellular functions including apoptosis, cell proliferation, RNA processing and cytoskeleton remodelling. By using a panel of pharmacological inhibitors, we identified discoidin domain receptor 1 (DDR1), a receptor tyrosine kinase specifically activated by collagen, as the initiator of 3D COL1-induced apoptosis. Our data support the concept that MT1-MMP contributes to the inactivation of the DDR1-BIK signalling axis through the cleavage of collagen fibres and/or the alteration of DDR1 receptor signalling unit, without triggering a drastic remodelling of the transcriptome of MCF-7 cells.

Phenotype classification of zebrafish embryos by supervised learning.

Zebrafish is increasingly used to assess biological properties of chemical substances and thus is becoming a specific tool for toxicological and pharmacological studies. The effects of chemical substances on embryo survival and development are generally evaluated manually through microscopic observation by an expert and documented by several typical photographs. Here, we present a methodology to automatically classify brightfield images of wildtype zebrafish embryos according to their defects by using an image analysis approach based on supervised machine learning. We show that, compared to manual classification, automatic classification results in 90 to 100% agreement with consensus voting of biological experts in nine out of eleven considered defects in 3 days old zebrafish larvae. Automation of the analysis and classification of zebrafish embryo pictures reduces the workload and time required for the biological expert and increases the reproducibility and objectivity of this classification.

Mapping gene regulatory networks in Drosophila eye development by large-scale transcriptome perturbations and motif inference.

Genome control is operated by transcription factors (TFs) controlling their target genes by binding to promoters and enhancers. Conceptually, the interactions between TFs, their binding sites, and their functional targets are represented by gene regulatory networks (GRNs). Deciphering in vivo GRNs underlying organ development in an unbiased genome-wide setting involves identifying both functional TF-gene interactions and physical TF-DNA interactions. To reverse engineer the GRNs of eye development in Drosophila, we performed RNA-seq across 72 genetic perturbations and sorted cell types and inferred a coexpression network. Next, we derived direct TF-DNA interactions using computational motif inference, ultimately connecting 241 TFs to 5,632 direct target genes through 24,926 enhancers. Using this network, we found network motifs, cis-regulatory codes, and regulators of eye development. We validate the predicted target regions of Grainyhead by ChIP-seq and identify this factor as a general cofactor in the eye network, being bound to thousands of nucleosome-free regions.

Bridging physiological and evolutionary time-scales in a gene regulatory network.

Gene regulatory networks (GRNs) govern phenotypic adaptations and reflect the trade-offs between physiological responses and evolutionary adaptation that act at different time-scales. To identify patterns of molecular function and genetic diversity in GRNs, we studied the drought response of the common sunflower, Helianthus annuus, and how the underlying GRN is related to its evolution. We examined the responses of 32,423 expressed sequences to drought and to abscisic acid (ABA) and selected 145 co-expressed transcripts. We characterized their regulatory relationships in nine kinetic studies based on different hormones. From this, we inferred a GRN by meta-analyses of a Gaussian graphical model and a random forest algorithm and studied the genetic differentiation among populations (FST ) at nodes. We identified two main hubs in the network that transport nitrate in guard cells. This suggests that nitrate transport is a critical aspect of the sunflower physiological response to drought. We observed that differentiation of the network genes in elite sunflower cultivars is correlated with their position and connectivity. This systems biology approach combined molecular data at different time-scales and identified important physiological processes. At the evolutionary level, we propose that network topology could influence responses to human selection and possibly adaptation to dry environments.