Outlier detection using iterative adaptive mini-minimum spanning tree generation with applications on medical data.

As an important technique for data pre-processing, outlier detection plays a crucial role in various real applications and has gained substantial attention, especially in medical fields. Despite the importance of outlier detection, many existing methods are vulnerable to the distribution of outliers and require prior knowledge, such as the outlier proportion. To address this problem to some extent, this article proposes an adaptive mini-minimum spanning tree-based outlier detection (MMOD) method, which utilizes a novel distance measure by scaling the Euclidean distance. For datasets containing different densities and taking on different shapes, our method can identify outliers without prior knowledge of outlier percentages. The results on both real-world medical data corpora and intuitive synthetic datasets demonstrate the effectiveness of the proposed method compared to state-of-the-art methods.

Image Synthesis and Modified BlendMask Instance Segmentation for Automated Nanoparticle Phenotyping.

Automated nanoparticle phenotyping is a critical aspect of high-throughput drug research, which requires analyzing nanoparticle size, shape, and surface topography from microscopy images. To automate this process, we present an instance segmentation pipeline that partitions individual nanoparticles on microscopy images. Our pipeline makes two key contributions. Firstly, we synthesize diverse and approximately realistic nanoparticle images to improve robust learning. Secondly, we improve the BlendMask model to segment tiny, overlapping, or sparse particle images. Specifically, we propose a parameterized approach for generating novel pairs of single particles and their masks, encouraging greater diversity in the training data. To synthesize more realistic particle images, we explore three particle placement rules and an image selection criterion. The improved one-stage instance segmentation network extracts distinctive features of nanoparticles and their context at both local and global levels, which addresses the data challenges associated with tiny, overlapping, or sparse nanoparticles. Extensive experiments demonstrate the effectiveness of our pipeline for automating nanoparticle partitioning and phenotyping in drug research using microscopy images.

Caps Captioning: A Modern Image Captioning Approach Based on Improved Capsule Network.

In image captioning models, the main challenge in describing an image is identifying all the objects by precisely considering the relationships between the objects and producing various captions. Over the past few years, many methods have been proposed, from an attribute-to-attribute comparison approach to handling issues related to semantics and their relationships. Despite the improvements, the existing techniques suffer from inadequate positional and geometrical attributes concepts. The reason is that most of the abovementioned approaches depend on Convolutional Neural Networks (CNNs) for object detection. CNN is notorious for failing to detect equivariance and rotational invariance in objects. Moreover, the pooling layers in CNNs cause valuable information to be lost. Inspired by the recent successful approaches, this paper introduces a novel framework for extracting meaningful descriptions based on a parallelized capsule network that describes the content of images through a high level of understanding of the semantic contents of an image. The main contribution of this paper is proposing a new method that not only overrides the limitations of CNNs but also generates descriptions with a wide variety of words by using Wikipedia. In our framework, capsules focus on the generation of meaningful descriptions with more detailed spatial and geometrical attributes for a given set of images by considering the position of the entities as well as their relationships. Qualitative experiments on the benchmark dataset MS-COCO show that our framework outperforms state-of-the-art image captioning models when describing the semantic content of the images.

Establishing a consensus for the hallmarks of cancer based on gene ontology and pathway annotations.

BACKGROUND: The hallmarks of cancer provide a highly cited and well-used conceptual framework for describing the processes involved in cancer cell development and tumourigenesis. However, methods for translating these high-level concepts into data-level associations between hallmarks and genes (for high throughput analysis), vary widely between studies. The examination of different strategies to associate and map cancer hallmarks reveals significant differences, but also consensus. RESULTS: Here we present the results of a comparative analysis of cancer hallmark mapping strategies, based on Gene Ontology and biological pathway annotation, from different studies. By analysing the semantic similarity between annotations, and the resulting gene set overlap, we identify emerging consensus knowledge. In addition, we analyse the differences between hallmark and gene set associations using Weighted Gene Co-expression Network Analysis and enrichment analysis. CONCLUSIONS: Reaching a community-wide consensus on how to identify cancer hallmark activity from research data would enable more systematic data integration and comparison between studies. These results highlight the current state of the consensus and offer a starting point for further convergence. In addition, we show how a lack of consensus can lead to large differences in the biological interpretation of downstream analyses and discuss the challenges of annotating changing and accumulating biological data, using intermediate knowledge resources that are also changing over time.

A Novel Function of TLR2 and MyD88 in the Regulation of Leukocyte Cell Migration Behavior During Wounding in Zebrafish Larvae.

Toll-like receptor (TLR) signaling via myeloid differentiation factor 88 protein (MyD88) has been indicated to be involved in the response to wounding. It remains unknown whether the putative role of MyD88 in wounding responses is due to a control of leukocyte cell migration. The aim of this study was to explore in vivo whether TLR2 and MyD88 are involved in modulating neutrophil and macrophage cell migration behavior upon zebrafish larval tail wounding. Live cell imaging of tail-wounded larvae was performed in tlr2 and myd88 mutants and their corresponding wild type siblings. In order to visualize cell migration following tissue damage, we constructed double transgenic lines with fluorescent markers for macrophages and neutrophils in all mutant and sibling zebrafish lines. Three days post fertilization (dpf), tail-wounded larvae were studied using confocal laser scanning microscopy (CLSM) to quantify the number of recruited cells at the wounding area. We found that in both tlr2-/- and myd88-/- groups the recruited neutrophil and macrophage numbers are decreased compared to their wild type sibling controls. Through analyses of neutrophil and macrophage migration patterns, we demonstrated that both tlr2 and myd88 control the migration direction of distant neutrophils upon wounding. Furthermore, in both the tlr2 and the myd88 mutants, macrophages migrated more slowly toward the wound edge. Taken together, our findings show that tlr2 and myd88 are involved in responses to tail wounding by regulating the behavior and speed of leukocyte migration in vivo.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions.

MOTIVATION: The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3' untranslated regions (3'UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies. RESULTS: We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3'UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3'UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses. AVAILABILITY AND IMPLEMENTATION: The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Anti-tuberculosis effect of isoniazid scales accurately from zebrafish to humans.

BACKGROUND AND PURPOSE: There is a clear need for innovation in anti-tuberculosis drug development. The zebrafish larva is an attractive disease model in tuberculosis research. To translate pharmacological findings to higher vertebrates, including humans, the internal exposure of drugs needs to be quantified and linked to observed response. EXPERIMENTAL APPROACH: In zebrafish studies, drugs are usually dissolved in the external water, posing a challenge to quantify internal exposure. We developed experimental methods to quantify internal exposure, including nanoscale blood sampling, and to quantify the bacterial burden, using automated fluorescence imaging analysis, with isoniazid as the test compound. We used pharmacokinetic-pharmacodynamic modelling to quantify the exposure-response relationship responsible for the antibiotic response. To translate isoniazid response to humans, quantitative exposure-response relationships in zebrafish were linked to simulated concentration-time profiles in humans, and two quantitative translational factors on sensitivity to isoniazid and stage of infection were included. KEY RESULTS: Blood concentration was only 20% of the external drug concentration. The bacterial burden increased exponentially, and an isoniazid dose corresponding to 15 mg·L-1 internal concentration (minimum inhibitory concentration) leads to bacteriostasis of the mycobacterial infection in the zebrafish. The concentration-effect relationship was quantified, and based on that relationship and the translational factors, the isoniazid response was translated to humans, which correlated well with observed data. CONCLUSIONS AND IMPLICATIONS: This proof of concept study confirmed the potential of zebrafish larvae as tuberculosis disease models in translational pharmacology and contributes to innovative anti-tuberculosis drug development, which is very clearly needed.

Automated image analysis system for studying cardiotoxicity in human pluripotent stem cell-Derived cardiomyocytes.

BACKGROUND: Cardiotoxicity, characterized by severe cardiac dysfunction, is a major problem in patients treated with different classes of anticancer drugs. Development of predictable human-based models and assays for drug screening are crucial for preventing potential drug-induced adverse effects. Current animal in vivo models and cell lines are not always adequate to represent human biology. Alternatively, human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) show great potential for disease modelling and drug-induced toxicity screenings. Fully automated high-throughput screening of drug toxicity on hiPSC-CMs by fluorescence image analysis is, however, very challenging, due to clustered cell growth patterns and strong intracellular and intercellular variation in the expression of fluorescent markers. RESULTS: In this paper, we report on the development of a fully automated image analysis system for quantification of cardiotoxic phenotypes from hiPSC-CMs that are treated with various concentrations of anticancer drugs doxorubicin or crizotinib. This high-throughput system relies on single-cell segmentation by nuclear signal extraction, fuzzy C-mean clustering of cardiac α-actinin signal, and finally nuclear signal propagation. When compared to manual segmentation, it generates precision and recall scores of 0.81 and 0.93, respectively. CONCLUSIONS: Our results show that our fully automated image analysis system can reliably segment cardiomyocytes even with heterogeneous α-actinin signals.

Fungal metabarcoding data integration framework for the MycoDiversity DataBase (MDDB).

Fungi have crucial roles in ecosystems, and are important associates for many organisms. They are adapted to a wide variety of habitats, however their global distribution and diversity remains poorly documented. The exponential growth of DNA barcode information retrieved from the environment is assisting considerably the traditional ways for unraveling fungal diversity and detection. The raw DNA data in association to environmental descriptors of metabarcoding studies are made available in public sequence read archives. While this is potentially a valuable source of information for the investigation of Fungi across diverse environmental conditions, the annotation used to describe environment is heterogenous. Moreover, a uniform processing pipeline still needs to be applied to the available raw DNA data. Hence, a comprehensive framework to analyses these data in a large context is still lacking. We introduce the MycoDiversity DataBase, a database which includes public fungal metabarcoding data of environmental samples for the study of biodiversity patterns of Fungi. The framework we propose will contribute to our understanding of fungal biodiversity and aims to become a valuable source for large-scale analyses of patterns in space and time, in addition to assisting evolutionary and ecological research on Fungi.

Predicting Metabolism from Gene Expression in an Improved Whole-Genome Metabolic Network Model of Danio rerio.

Zebrafish is a useful modeling organism for the study of vertebrate development, immune response, and metabolism. Metabolic studies can be aided by mathematical reconstructions of the metabolic network of zebrafish. These list the substrates and products of all biochemical reactions that occur in the zebrafish. Mathematical techniques such as flux-balance analysis then make it possible to predict the possible metabolic flux distributions that optimize, for example, the turnover of food into biomass. The only available genome-scale reconstruction of zebrafish metabolism is ZebraGEM. In this study, we present ZebraGEM 2.0, an updated and validated version of ZebraGEM. ZebraGEM 2.0 is extended with gene-protein-reaction associations (GPRs) that are required to integrate genetic data with the metabolic model. To demonstrate the use of these GPRs, we performed an in silico genetic screening for knockouts of metabolic genes and validated the results against published in vivo genetic knockout and knockdown screenings. Among the single knockout simulations, we identified 74 essential genes, whose knockout stopped growth completely. Among these, 11 genes are known have an abnormal knockout or knockdown phenotype in vivo (partial), and 41 have human homologs associated with metabolic diseases. We also added the oxidative phosphorylation pathway, which was unavailable in the published version of ZebraGEM. The updated model performs better than the original model on a predetermined list of metabolic functions. We also determined a minimal feed composition. The oxidative phosphorylation pathways were validated by comparing with published experiments in which key components of the oxidative phosphorylation pathway were pharmacologically inhibited. To test the utility of ZebraGEM2.0 for obtaining new results, we integrated gene expression data from control and Mycobacterium marinum-infected zebrafish larvae. The resulting model predicts impeded growth and altered histidine metabolism in the infected larvae.

An efficient and robust hybrid method for segmentation of zebrafish objects from bright-field microscope images.

Accurate segmentation of zebrafish from bright-field microscope images is crucial to many applications in the life sciences. Early zebrafish stages are used, and in these stages the zebrafish is partially transparent. This transparency leads to edge ambiguity as is typically seen in the larval stages. Therefore, segmentation of zebrafish objects from images is a challenging task in computational bio-imaging. Popular computational methods fail to segment the relevant edges, which subsequently results in inaccurate measurements and evaluations. Here we present a hybrid method to accomplish accurate and efficient segmentation of zebrafish specimens from bright-field microscope images. We employ the mean shift algorithm to augment the colour representation in the images. This improves the discrimination of the specimen to the background and provides a segmentation candidate retaining the overall shape of the zebrafish. A distance-regularised level set function is initialised from this segmentation candidate and fed to an improved level set method, such that we can obtain another segmentation candidate which preserves the explicit contour of the object. The two candidates are fused using heuristics, and the hybrid result is refined to represent the contour of the zebrafish specimen. We have applied the proposed method on two typical datasets. From experiments, we conclude that the proposed hybrid method improves both efficiency and accuracy of the segmentation of the zebrafish specimen. The results are going to be used for high-throughput applications with zebrafish.

Three-dimensional reconstruction and measurements of zebrafish larvae from high-throughput axial-view in vivo imaging.

High-throughput imaging is applied to provide observations for accurate statements on phenomena in biology and this has been successfully applied in the domain of cells, i.e. cytomics. In the domain of whole organisms, we need to take the hurdles to ensure that the imaging can be accomplished with a sufficient throughput and reproducibility. For vertebrate biology, zebrafish is a popular model system for High-throughput applications. The development of the Vertebrate Automated Screening Technology (VAST BioImager), a microscope mounted system, enables the application of zebrafish high-throughput screening. The VAST BioImager contains a capillary that holds a zebrafish for imaging. Through the rotation of the capillary, multiple axial-views of a specimen can be acquired. For the VAST BioImager, fluorescence and/or confocal microscopes are used. Quantitation of a specific signal as derived from a label in one fluorescent channel requires insight in the zebrafish volume to be able to normalize quantitation to volume units. However, from the setup of the VAST BioImager, a specimen volume cannot be straightforwardly derived. We present a high-throughput axial-view imaging architecture based on the VAST BioImager. We propose profile-based 3D reconstruction to produce 3D volumetric representations for zebrafish larvae using the axial-views. Volume and surface area can then be derived from the 3D reconstruction to obtain the shape characteristics in high-throughput measurements. In addition, we develop a calibration and a validation of our methodology. From our measurements we show that with a limited amount of views, accurate measurements of volume and surface area for zebrafish larvae can be obtained. We have applied the proposed method on a range of developmental stages in zebrafish and produced metrical references for the volume and surface area for each stage.

Fast Post-Processing Pipeline for Optical Projection Tomography.

To improve the effectiveness and efficiency of optical projection tomography (OPT) 3-D reconstruction, we present a fast post-processing pipeline, including cropping, background subtraction, center of rotation (COR) correction, and 3-D reconstruction. Regarding to the COR correction, a novel algorithm based on interest point detection of sinogram is proposed by considering the principle of OPT imaging. Instead of locating the COR on single sinogram, we select equally spaced sinograms in the detected full range of specimen to make the located COR more convincing. The presented post-processing pipeline is implemented in a parallel manner and the experiments show that the average runtime for each image of size 1036 ×1360 ×400 pixels is less than 1 min. To quantify and compare the reconstructed results of different COR correction approaches, the coefficient of variation instead of variance is employed. The results indicate that the proposed COR correction outperforms the three traditional COR alignment approaches in terms of effectiveness and computational complexity.

Modeling biological gradient formation: combining partial differential equations and Petri nets.

Both Petri nets and differential equations are important modeling tools for biological processes. In this paper we demonstrate how these two modeling techniques can be combined to describe biological gradient formation. Parameters derived from partial differential equation describing the process of gradient formation are incorporated in an abstract Petri net model. The quantitative aspects of the resulting model are validated through a case study of gradient formation in the fruit fly.

Cellular adhesome screen identifies critical modulators of focal adhesion dynamics, cellular traction forces and cell migration behaviour.

Cancer cells migrate from the primary tumour into surrounding tissue in order to form metastasis. Cell migration is a highly complex process, which requires continuous remodelling and re-organization of the cytoskeleton and cell-matrix adhesions. Here, we aimed to identify genes controlling aspects of tumour cell migration, including the dynamic organization of cell-matrix adhesions and cellular traction forces. In a siRNA screen targeting most cell adhesion-related genes we identified 200+ genes that regulate size and/or dynamics of cell-matrix adhesions in MCF7 breast cancer cells. In a subsequent secondary screen, the 64 most effective genes were evaluated for growth factor-induced cell migration and validated by tertiary RNAi pool deconvolution experiments. Four validated hits showed significantly enlarged adhesions accompanied by reduced cell migration upon siRNA-mediated knockdown. Furthermore, loss of PPP1R12B, HIPK3 or RAC2 caused cells to exert higher traction forces, as determined by traction force microscopy with elastomeric micropillar post arrays, and led to considerably reduced force turnover. Altogether, we identified genes that co-regulate cell-matrix adhesion dynamics and traction force turnover, thereby modulating overall motility behaviour.

Hierarchical classification strategy for Phenotype extraction from epidermal growth factor receptor endocytosis screening.

BACKGROUND: Endocytosis is regarded as a mechanism of attenuating the epidermal growth factor receptor (EGFR) signaling and of receptor degradation. There is increasing evidence becoming available showing that breast cancer progression is associated with a defect in EGFR endocytosis. In order to find related Ribonucleic acid (RNA) regulators in this process, high-throughput imaging with fluorescent markers is used to visualize the complex EGFR endocytosis process. Subsequently a dedicated automatic image and data analysis system is developed and applied to extract the phenotype measurement and distinguish different developmental episodes from a huge amount of images acquired through high-throughput imaging. For the image analysis, a phenotype measurement quantifies the important image information into distinct features or measurements. Therefore, the manner in which prominent measurements are chosen to represent the dynamics of the EGFR process becomes a crucial step for the identification of the phenotype. In the subsequent data analysis, classification is used to categorize each observation by making use of all prominent measurements obtained from image analysis. Therefore, a better construction for a classification strategy will support to raise the performance level in our image and data analysis system. RESULTS: In this paper, we illustrate an integrated analysis method for EGFR signalling through image analysis of microscopy images. Sophisticated wavelet-based texture measurements are used to obtain a good description of the characteristic stages in the EGFR signalling. A hierarchical classification strategy is designed to improve the recognition of phenotypic episodes of EGFR during endocytosis. Different strategies for normalization, feature selection and classification are evaluated. CONCLUSIONS: The results of performance assessment clearly demonstrate that our hierarchical classification scheme combined with a selected set of features provides a notable improvement in the temporal analysis of EGFR endocytosis. Moreover, it is shown that the addition of the wavelet-based texture features contributes to this improvement. Our workflow can be applied to drug discovery to analyze defected EGFR endocytosis processes.

Machine Learning approach to discriminate Saccharomyces cerevisiae yeast cells using sophisticated image features.

In biological research, Saccharomyces cerevisiae yeast cells are used to study the behaviour of proteins. This is a time consuming and not completely objective process. Hence, Image analysis platforms are developed to address these problems and to offer analysis per cell as well. The robust segmentation algorithms implemented in such platforms enables us to apply a machine learning approach on the measured cells. Such approach is based on a set of relevant individual cell features extracted from the microscope images of the yeast cells. In this paper, we composed a set of features to represent the intensity and morphology characteristics in a more sophisticated way. These features are based on first and second order histograms and wavelet-based texture measurement. To show the discrimination power of these features, we built a classification model to discriminate between different groups. The building process involved evaluation of a set of classification systems, data sampling techniques, data normalization schemes and attribute selection algorithms. The results show a significant ability to discriminate different cell strains and conditions; subsequently it reveals the benefits of the classification model based on the introduced features. This model is promising in revealing subtle patterns in future high-throughput yeast studies.

Macrophage-expressed perforins mpeg1 and mpeg1.2 have an anti-bacterial function in zebrafish.

Macrophage-expressed gene 1 (MPEG1) encodes an evolutionarily conserved protein with a predicted membrane attack complex/perforin domain associated with host defence against invading pathogens. In vertebrates, MPEG1/perforin-2 is an integral membrane protein of macrophages, suspected to be involved in the killing of intracellular bacteria by pore-forming activity. Zebrafish have 3 copies of MPEG1; 2 are expressed in macrophages, whereas the third could be a pseudogene. The mpeg1 and mpeg1.2 genes show differential regulation during infection of zebrafish embryos with the bacterial pathogens Mycobacterium marinum and Salmonella typhimurium. While mpeg1 is downregulated during infection with both pathogens, mpeg1.2 is infection inducible. Upregulation of mpeg1.2 is partially dependent on the presence of functional Mpeg1 and requires the Toll-like receptor adaptor molecule MyD88 and the transcription factor NFκB. Knockdown of mpeg1 alters the immune response to M. marinum infection and results in an increased bacterial burden. In Salmonella typhimurium infection, both mpeg1 and mpeg1.2 knockdown increase the bacterial burdens, but mpeg1 morphants show increased survival times. The combined results of these two in vivo infection models support the anti-bacterial function of the MPEG1/perforin-2 family and indicate that the intricate cross-regulation of the two mpeg1 copies aids the zebrafish host in combatting infection of various pathogens.

Underlying molecular mechanisms of DIO2 susceptibility in symptomatic osteoarthritis.

OBJECTIVES: To investigate how the genetic susceptibility gene DIO2 confers risk to osteoarthritis (OA) onset in humans and to explore whether counteracting the deleterious effect could contribute to novel therapeutic approaches. METHODS: Epigenetically regulated expression of DIO2 was explored by assessing methylation of positional CpG-dinucleotides and the respective DIO2 expression in OA-affected and macroscopically preserved articular cartilage from end-stage OA patients. In a human in vitro chondrogenesis model, we measured the effects when thyroid signalling during culturing was either enhanced (excess T3 or lentiviral induced DIO2 overexpression) or decreased (iopanoic acid). RESULTS: OA-related changes in methylation at a specific CpG dinucleotide upstream of DIO2 caused significant upregulation of its expression (ß=4.96; p=0.0016). This effect was enhanced and appeared driven specifically by DIO2 rs225014 risk allele carriers (ß=5.58, p=0.0006). During in vitro chondrogenesis, DIO2 overexpression resulted in a significant reduced capacity of chondrocytes to deposit extracellular matrix (ECM) components, concurrent with significant induction of ECM degrading enzymes (ADAMTS5, MMP13) and markers of mineralisation (ALPL, COL1A1). Given their concurrent and significant upregulation of expression, this process is likely mediated via HIF-2α/RUNX2 signalling. In contrast, we showed that inhibiting deiodinases during in vitro chondrogenesis contributed to prolonged cartilage homeostasis as reflected by significant increased deposition of ECM components and attenuated upregulation of matrix degrading enzymes. CONCLUSIONS: Our findings show how genetic variation at DIO2 could confer risk to OA and raised the possibility that counteracting thyroid signalling may be a novel therapeutic approach.

Ultra high content image analysis and phenotype profiling of 3D cultured micro-tissues.

In many situations, 3D cell cultures mimic the natural organization of tissues more closely than 2D cultures. Conventional methods for phenotyping such 3D cultures use either single or multiple simple parameters based on morphology and fluorescence staining intensity. However, due to their simplicity many details are not taken into account which limits system-level study of phenotype characteristics. Here, we have developed a new image analysis platform to automatically profile 3D cell phenotypes with 598 parameters including morphology, topology, and texture parameters such as wavelet and image moments. As proof of concept, we analyzed mouse breast cancer cells (4T1 cells) in a 384-well plate format following exposure to a diverse set of compounds at different concentrations. The result showed concentration dependent phenotypic trajectories for different biologically active compounds that could be used to classify compounds based on their biological target. To demonstrate the wider applicability of our method, we analyzed the phenotypes of a collection of 44 human breast cancer cell lines cultured in 3D and showed that our method correctly distinguished basal-A, basal-B, luminal and ERBB2+ cell lines in a supervised nearest neighbor classification method.