A microengineered Brain-Chip to model neuroinflammation in humans.

Species differences in brain and blood-brain barrier (BBB) biology hamper the translation of findings from animal models to humans, impeding the development of therapeutics for brain diseases. Here, we present a human organotypic microphysiological system (MPS) that includes endothelial-like cells, pericytes, glia, and cortical neurons and maintains BBB permeability at in vivo relevant levels. This human Brain-Chip engineered to recapitulate critical aspects of the complex interactions that mediate neuroinflammation and demonstrates significant improvements in clinical mimicry compared to previously reported similar MPS. In comparison to Transwell culture, the transcriptomic profiling of the Brain-Chip displayed significantly advanced similarity to the human adult cortex and enrichment in key neurobiological pathways. Exposure to TNF-α recreated the anticipated inflammatory environment shown by glia activation, increased release of proinflammatory cytokines, and compromised barrier permeability. We report the development of a robust brain MPS for mechanistic understanding of cell-cell interactions and BBB function during neuroinflammation.

Analytics and visualization tools to characterize single-cell stochasticity using bacterial single-cell movie cytometry data.

BACKGROUND: Time-lapse microscopy live-cell imaging is essential for studying the evolution of bacterial communities at single-cell resolution. It allows capturing detailed information about the morphology, gene expression, and spatial characteristics of individual cells at every time instance of the imaging experiment. The image analysis of bacterial "single-cell movies" (videos) generates big data in the form of multidimensional time series of measured bacterial attributes. If properly analyzed, these datasets can help us decipher the bacterial communities' growth dynamics and identify the sources and potential functional role of intra- and inter-subpopulation heterogeneity. Recent research has highlighted the importance of investigating the role of biological "noise" in gene regulation, cell growth, cell division, etc. Single-cell analytics of complex single-cell movie datasets, capturing the interaction of multiple micro-colonies with thousands of cells, can shed light on essential phenomena for human health, such as the competition of pathogens and benign microbiome cells, the emergence of dormant cells ("persisters"), the formation of biofilms under different stress conditions, etc. However, highly accurate and automated bacterial bioimage analysis and single-cell analytics methods remain elusive, even though they are required before we can routinely exploit the plethora of data that single-cell movies generate. RESULTS: We present visualization and single-cell analytics using R (ViSCAR), a set of methods and corresponding functions, to visually explore and correlate single-cell attributes generated from the image processing of complex bacterial single-cell movies. They can be used to model and visualize the spatiotemporal evolution of attributes at different levels of the microbial community organization (i.e., cell population, colony, generation, etc.), to discover possible epigenetic information transfer across cell generations, infer mathematical and statistical models describing various stochastic phenomena (e.g., cell growth, cell division), and even identify and auto-correct errors introduced unavoidably during the bioimage analysis of a dense movie with thousands of overcrowded cells in the microscope's field of view. CONCLUSIONS: ViSCAR empowers researchers to capture and characterize the stochasticity, uncover the mechanisms leading to cellular phenotypes of interest, and decipher a large heterogeneous microbial communities' dynamic behavior. ViSCAR source code is available from GitLab at https://gitlab.com/ManolakosLab/viscar .

Modeling alpha-synuclein pathology in a human brain-chip to assess blood-brain barrier disruption.

Parkinson's disease and related synucleinopathies are characterized by the abnormal accumulation of alpha-synuclein aggregates, loss of dopaminergic neurons, and gliosis of the substantia nigra. Although clinical evidence and in vitro studies indicate disruption of the Blood-Brain Barrier in Parkinson's disease, the mechanisms mediating the endothelial dysfunction is not well understood. Here we leveraged the Organs-on-Chips technology to develop a human Brain-Chip representative of the substantia nigra area of the brain containing dopaminergic neurons, astrocytes, microglia, pericytes, and microvascular brain endothelial cells, cultured under fluid flow. Our αSyn fibril-induced model was capable of reproducing several key aspects of Parkinson's disease, including accumulation of phosphorylated αSyn (pSer129-αSyn), mitochondrial impairment, neuroinflammation, and compromised barrier function. This model may enable research into the dynamics of cell-cell interactions in human synucleinopathies and serve as a testing platform for target identification and validation of novel therapeutics.

A Novel Microphysiological Colon Platform to Decipher Mechanisms Driving Human Intestinal Permeability.

BACKGROUND & AIMS: The limited availability of organoid systems that mimic the molecular signatures and architecture of human intestinal epithelium has been an impediment to allowing them to be harnessed for the development of therapeutics as well as physiological insights. We developed a microphysiological Organ-on-Chip (Emulate, Inc, Boston, MA) platform designed to mimic properties of human intestinal epithelium leading to insights into barrier integrity. METHODS: We combined the human biopsy-derived leucine-rich repeat-containing G-protein-coupled receptor 5-positive organoids and Organ-on-Chip technologies to establish a micro-engineered human Colon Intestine-Chip (Emulate, Inc, Boston, MA). We characterized the proximity of the model to human tissue and organoids maintained in suspension by RNA sequencing analysis, and their differentiation to intestinal epithelial cells on the Colon Intestine-Chip under variable conditions. Furthermore, organoids from different donors were evaluated to understand variability in the system. Our system was applied to understanding the epithelial barrier and characterizing mechanisms driving the cytokine-induced barrier disruption. RESULTS: Our data highlight the importance of the endothelium and the in vivo tissue-relevant dynamic microenvironment in the Colon Intestine-Chip in the establishment of a tight monolayer of differentiated, polarized, organoid-derived intestinal epithelial cells. We confirmed the effect of interferon-γ on the colonic barrier and identified reorganization of apical junctional complexes, and induction of apoptosis in the intestinal epithelial cells as mediating mechanisms. We show that in the human Colon Intestine-Chip exposure to interleukin 22 induces disruption of the barrier, unlike its described protective role in experimental colitis in mice. CONCLUSIONS: We developed a human Colon Intestine-Chip platform and showed its value in the characterization of the mechanism of action of interleukin 22 in the human epithelial barrier. This system can be used to elucidate, in a time- and challenge-dependent manner, the mechanism driving the development of leaky gut in human beings and to identify associated biomarkers.

An information-theoretic approach for measuring the distance of organ tissue samples using their transcriptomic signatures.

MOTIVATION: Recapitulating aspects of human organ functions using in vitro (e.g. plates, transwells, etc.), in vivo (e.g. mouse, rat, etc.), or ex vivo (e.g. organ chips, 3D systems, etc.) organ models is of paramount importance for drug discovery and precision medicine. It will allow us to identify potential side effects and test the effectiveness of new therapeutic approaches early in their design phase, and will inform the development of better disease models. Developing mathematical methods to reliably compare the 'distance/similarity' of organ models from/to the real human organ they represent is an understudied problem with important applications in biomedicine and tissue engineering. RESULTS: We introduce the Transcriptomic Signature Distance (TSD), an information-theoretic distance for assessing the transcriptomic similarity of two tissue samples, or two groups of tissue samples. In developing TSD, we are leveraging next-generation sequencing data as well as information retrieved from well-curated databases providing signature gene sets characteristic for human organs. We present the justification and mathematical development of the new distance and demonstrate its effectiveness and advantages in different scenarios of practical importance using several publicly available RNA-seq datasets. AVAILABILITY AND IMPLEMENTATION: The computation of both TSD versions (simple and weighted) has been implemented in R and can be downloaded from https://github.com/Cod3B3nd3R/Transcriptomic-Signature-Distance. CONTACT: dimitris.manatakis@emulatebio.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multi-criteria protein structure comparison and structural similarities analysis using pyMCPSC.

Protein Structure Comparison (PSC) is a well developed field of computational proteomics with active interest from the research community, since it is widely used in structural biology and drug discovery. With new PSC methods continuously emerging and no clear method of choice, Multi-Criteria Protein Structure Comparison (MCPSC) is commonly employed to combine methods and generate consensus structural similarity scores. We present pyMCPSC, a Python based utility we developed to allow users to perform MCPSC efficiently, by exploiting the parallelism afforded by the multi-core CPUs of today's desktop computers. We show how pyMCPSC facilitates the analysis of similarities in protein domain datasets and how it can be extended to incorporate new PSC methods as they are becoming available. We exemplify the power of pyMCPSC using a case study based on the Proteus_300 dataset. Results generated using pyMCPSC show that MCPSC scores form a reliable basis for identifying the true classification of a domain, as evidenced both by the ROC analysis as well as the Nearest-Neighbor analysis. Structure similarity based "Phylogenetic Trees" representation generated by pyMCPSC provide insight into functional grouping within the dataset of domains. Furthermore, scatter plots generated by pyMCPSC show the existence of strong correlation between protein domains belonging to SCOP Class C and loose correlation between those of SCOP Class D. Such analyses and corresponding visualizations help users quickly gain insights about their datasets. The source code of pyMCPSC is available under the GPLv3.0 license through a GitHub repository (https://github.com/xulesc/pymcpsc).

Image analysis driven single-cell analytics for systems microbiology.

BACKGROUND: Time-lapse microscopy is an essential tool for capturing and correlating bacterial morphology and gene expression dynamics at single-cell resolution. However state-of-the-art computational methods are limited in terms of the complexity of cell movies that they can analyze and lack of automation. The proposed Bacterial image analysis driven Single Cell Analytics (BaSCA) computational pipeline addresses these limitations thus enabling high throughput systems microbiology. RESULTS: BaSCA can segment and track multiple bacterial colonies and single-cells, as they grow and divide over time (cell segmentation and lineage tree construction) to give rise to dense communities with thousands of interacting cells in the field of view. It combines advanced image processing and machine learning methods to deliver very accurate bacterial cell segmentation and tracking (F-measure over 95%) even when processing images of imperfect quality with several overcrowded colonies in the field of view. In addition, BaSCA extracts on the fly a plethora of single-cell properties, which get organized into a database summarizing the analysis of the cell movie. We present alternative ways to analyze and visually explore the spatiotemporal evolution of single-cell properties in order to understand trends and epigenetic effects across cell generations. The robustness of BaSCA is demonstrated across different imaging modalities and microscopy types. CONCLUSIONS: BaSCA can be used to analyze accurately and efficiently cell movies both at a high resolution (single-cell level) and at a large scale (communities with many dense colonies) as needed to shed light on e.g. how bacterial community effects and epigenetic information transfer play a role on important phenomena for human health, such as biofilm formation, persisters' emergence etc. Moreover, it enables studying the role of single-cell stochasticity without losing sight of community effects that may drive it.

Efficient Multicriteria Protein Structure Comparison on Modern Processor Architectures.

Fast increasing computational demand for all-to-all protein structures comparison (PSC) is a result of three confounding factors: rapidly expanding structural proteomics databases, high computational complexity of pairwise protein comparison algorithms, and the trend in the domain towards using multiple criteria for protein structures comparison (MCPSC) and combining results. We have developed a software framework that exploits many-core and multicore CPUs to implement efficient parallel MCPSC in modern processors based on three popular PSC methods, namely, TMalign, CE, and USM. We evaluate and compare the performance and efficiency of the two parallel MCPSC implementations using Intel's experimental many-core Single-Chip Cloud Computer (SCC) as well as Intel's Core i7 multicore processor. We show that the 48-core SCC is more efficient than the latest generation Core i7, achieving a speedup factor of 42 (efficiency of 0.9), making many-core processors an exciting emerging technology for large-scale structural proteomics. We compare and contrast the performance of the two processors on several datasets and also show that MCPSC outperforms its component methods in grouping related domains, achieving a high F-measure of 0.91 on the benchmark CK34 dataset. The software implementation for protein structure comparison using the three methods and combined MCPSC, along with the developed underlying rckskel algorithmic skeletons library, is available via GitHub.

Tracking single-cells in overcrowded bacterial colonies.

Cell tracking enables data extraction from timelapse "cell movies" and promotes modeling biological processes at the single-cell level. We introduce a new fully automated computational strategy to track accurately cells across frames in time-lapse movies. Our method is based on a dynamic neighborhoods formation and matching approach, inspired by motion estimation algorithms for video compression. Moreover, it exploits "divide and conquer" opportunities to solve effectively the challenging cells tracking problem in overcrowded bacterial colonies. Using cell movies generated by different labs we demonstrate that the accuracy of the proposed method remains very high (exceeds 97%) even when analyzing large overcrowded microbial colonies.

In silico modeling of the effects of alpha-synuclein oligomerization on dopaminergic neuronal homeostasis.

BACKGROUND: Alpha-synuclein (ASYN) is central in Parkinson's disease (PD) pathogenesis. Converging pieces of evidence suggest that the levels of ASYN expression play a critical role in both familial and sporadic Parkinson's disease. ASYN fibrils are the main component of inclusions called Lewy Bodies (LBs) which are found mainly in the surviving neurons of the substantia nigra. Despite the accumulated knowledge regarding the involvement of ASYN in molecular mechanisms underlying the development of PD, there is much information missing which prevents understanding the causes of the disease and how to stop its progression. RESULTS: Using a Systems Biology approach, we develop a biomolecular reactions model that describes the intracellular ASYN dynamics in relation to overexpression, post-translational modification, oligomerization and degradation of the protein. Especially for the proteolysis of ASYN, the model takes into account the biological knowledge regarding the contribution of Chaperone Mediated Autophagy (CMA), macro-autophagic and proteasome pathways in the protein's degradation. Importantly, inhibitory phenomena, caused by ASYN, concerning CMA (more specifically the lysosomal-associated membrane protein 2a, abbreviated as Lamp2a receptor, which is the rate limiting step of CMA) and the proteasome are carefully modeled. The model is validated by simulation studies of known experimental overexpression data from SH-SY5Y cells and the unknown model parameters are estimated either computationally or by experimental fitting. The calibrated model is then tested under three hypothetical intervention scenarios and in all cases predicts increased cell viability that agrees with experimental evidence. The biomodel has been annotated and is made available in SBML format. CONCLUSIONS: The mathematical model presented here successfully simulates the dynamic phenomena of ASYN overexpression and oligomerization and predicts the biological system's behavior in a number of scenarios not used for model calibration. It allows, for the first time, to qualitatively estimate the protein levels that are capable of deregulating proteolytic homeostasis. In addition, it can help form new hypotheses for intervention that could be tested experimentally.

A learning scheme for reach to grasp movements: on EMG-based interfaces using task specific motion decoding models.

A learning scheme based on random forests is used to discriminate between different reach to grasp movements in 3-D space, based on the myoelectric activity of human muscles of the upper-arm and the forearm. Task specificity for motion decoding is introduced in two different levels: Subspace to move toward and object to be grasped. The discrimination between the different reach to grasp strategies is accomplished with machine learning techniques for classification. The classification decision is then used in order to trigger an EMG-based task-specific motion decoding model. Task specific models manage to outperform "general" models providing better estimation accuracy. Thus, the proposed scheme takes advantage of a framework incorporating both a classifier and a regressor that cooperate advantageously in order to split the task space. The proposed learning scheme can be easily used to a series of EMG-based interfaces that must operate in real time, providing data-driven capabilities for multiclass problems, that occur in everyday life complex environments.

Visualizing meta-features in proteomic maps.

BACKGROUND: The steps of a high-throughput proteomics experiment include the separation, differential expression and mass spectrometry-based identification of proteins. However, the last and more challenging step is inferring the biological role of the identified proteins through their association with interaction networks, biological pathways, analysis of the effect of post-translational modifications, and other protein-related information. RESULTS: In this paper, we present an integrative visualization methodology that allows combining experimentally produced proteomic features with protein meta-features, typically coming from meta-analysis tools and databases, in synthetic Proteomic Feature Maps. Using three proteomics analysis scenarios, we show that the proposed visualization approach is effective in filtering, navigating and interacting with the proteomics data in order to address visually challenging biological questions. The novelty of our approach lies in the ease of integration of any user-defined proteomic features in easy-to-comprehend visual representations that resemble the familiar 2D-gel images, and can be adapted to the user's needs. The main capabilities of the developed VIP software, which implements the presented visualization methodology, are also highlighted and discussed. CONCLUSIONS: By using this visualization and the associated VIP software, researchers can explore a complex heterogeneous proteomics dataset from different perspectives in order to address visually important biological queries and formulate new hypotheses for further investigation. VIP is freely available at http://pelopas.uop.gr/~egian/VIP/index.html.

Protein spot detection and quantification in 2-DE gel images using machine-learning methods.

Two-dimensional gel electrophoresis (2-DE) is the most established protein separation method used in expression proteomics. Despite the existence of sophisticated software tools, 2-DE gel image analysis still remains a serious bottleneck. The low accuracies of commercial software packages and the extensive manual calibration that they often require for acceptable results show that we are far from achieving the goal of a fully automated and reliable, high-throughput gel processing system. We present a novel spot detection and quantification methodology which draws heavily from unsupervised machine-learning methods. Using the proposed hierarchical machine learning-based segmentation methodology reduces both the number of faint spots missed (improves sensitivity) and the number of extraneous spots introduced (improves precision). The detection and quantification performance has been thoroughly evaluated and is shown to compare favorably (higher F-measure) to a commercially available software package (PDQuest). The whole image analysis pipeline that we have developed is fully automated and can be used for high-throughput proteomics analysis since it does not require any manual intervention for recalibration every time a new 2-DE gel image is to be analyzed. Furthermore, it can be easily parallelized for high performance and also applied without any modification to prealigned group average gels.

Prediction of biological integrity based on environmental similarity--revealing the scale-dependent link between study area and top environmental predictors.

Indices of Biological integrity (IBI) are considered valid indicators of the overall health of a water body because the biological community is an endpoint within natural systems. However, prediction of biological integrity using information from multi-parameter environmental observations is a challenging problem due to the hierarchical organization of the natural environment, the existence of nonlinear inter-dependencies among variables as well as natural stochasticity and measurement noise. We present a method for predicting the Fish Index of Biological Integrity (IBI) using multiple environmental observations at the state-scale in Ohio. Instream (chemical and physical quality) and offstream parameters (regional and local upstream land uses, stream fragmentation, and point source density and intensity) are used for this purpose. The IBI predictions are obtained using the environmental site-similarity concept and following a simple to implement leave-one-out cross validation approach. An IBI prediction for a sampling site is calculated by averaging the observed IBI scores of observations clustered in the most similar branch of a dendrogram--a hierarchical clustering tree of environmental observations--built using the rest of the observations. The standardized Euclidean distance is used to assess dissimilarity between observations. The constructed predictive model was able to explain 61% of the IBI variability statewide. Stream fragmentation and regional land use explained 60% of the variability; the remaining 1% was explained by instream habitat quality. Metrics related to local land use, water quality, and point source density and intensity did not improve the predictive model at the state-scale. The impact of local environmental conditions was evaluated by comparing local characteristics between well- and mispredicted sites. Significant differences in local land use patterns and upstream fragmentation density explained some of the model's over-predictions. Local land use conditions explained some of the model's IBI under-predictions at the state-scale since none of the variables within this group were included in the best final predictive model. Under-predicted sites also had higher levels of downstream fragmentation. The proposed variables ranking and predictive modeling methodology is very well suited for the analysis of hierarchical environments, such as natural fresh water systems, with many cross-correlated environmental variables. It is computationally efficient, can be fully automated, does not make any pre-conceived assumptions on the variables interdependency structure (such as linearity), and it is able to rank variables in a database and generate IBI predictions using only non-parametric easy to implement hierarchical clustering.

Improving 2-DE gel image denoising using contourlets.

One of the most commonly used methods for protein separation is 2-DE. After 2-DE gel scanning, images with a plethora of spot features emerge that are usually contaminated by inherent noise. The objective of the denoising process is to remove noise to the extent that the true spots are recovered correctly and accurately i.e. without introducing distortions leading to the detection of false-spot features. In this paper we propose and justify the use of the contourlet transform as a tool for 2-DE gel images denoising. We compare its effectiveness with state-of-the-art methods such as wavelets-based multiresolution image analysis and spatial filtering. We show that contourlets not only achieve better average S/N performance than wavelets and spatial filters, but also preserve better spot boundaries and faint spots and alter less the intensities of informative spot features, leading to more accurate spot volume estimation and more reliable spot detection, operations that are essential to differential expression proteomics for biomarkers discovery.

Proteomic feature maps: a new visualization approach in proteomics analysis.

The different steps of a proteomics analysis workflow generate a plethora of features for each extracted proteomic object (a protein spot in 2D gel electrophoresis (2-DE), or a peptide peak in liquid chromatography-mass spectrometry (LC-MS) analysis). Yet, the joint visualization of multiple object features on 2D gel-like maps is rather limited in currently available proteomics software packages. We introduce a new, simple, and intuitive visualization method that utilizes spheres to represent proteomic objects on proteomic feature maps, and exploits the spheres size and color to provide simultaneous visualization of user-selected feature pairs. Our contribution, a unified and flexible visualization mechanism that can be easily applied at any stage of a 2-DE or a LC-MS based differential proteomics study, is demonstrated and discussed using five representative scenarios. The joint visualization of proteomic object features and their spatial distribution is a powerful tool for inspecting and comparing the proteomics analysis results, attracting the users attention to useful information, such as differential expression trends and patterns, and even assisting in the evaluation and refinement of a proteomics experiment.