Time-resolved molecular measurements reveal changes in astronauts during spaceflight.

From the early days of spaceflight to current missions, astronauts continue to be exposed to multiple hazards that affect human health, including low gravity, high radiation, isolation during long-duration missions, a closed environment and distance from Earth. Their effects can lead to adverse physiological changes and necessitate countermeasure development and/or longitudinal monitoring. A time-resolved analysis of biological signals can detect and better characterize potential adverse events during spaceflight, ideally preventing them and maintaining astronauts' wellness. Here we provide a time-resolved assessment of the impact of spaceflight on multiple astronauts (n = 27) by studying multiple biochemical and immune measurements before, during, and after long-duration orbital spaceflight. We reveal space-associated changes of astronauts' physiology on both the individual level and across astronauts, including associations with bone resorption and kidney function, as well as immune-system dysregulation.

Time-resolved molecular measurements reveal changes in astronauts during spaceflight.

From the early days of spaceflight to current missions, astronauts continue to be exposed to multiple hazards that affect human health, including low gravity, high radiation, isolation during long-duration missions, a closed environment and distance from Earth. Their effects can lead to adverse physiological changes and necessitate countermeasure development and/or longitudinal monitoring. A time-resolved analysis of biological signals can detect and better characterize potential adverse events during spaceflight, ideally preventing them and maintaining astronauts' wellness. Here we provide a time-resolved assessment of the impact of spaceflight on multiple astronauts (n=27) by studying multiple biochemical and immune measurements before, during, and after long-duration orbital spaceflight. We reveal space-associated changes of astronauts' physiology on both the individual level and across astronauts, including associations with bone resorption and kidney function, as well as immune-system dysregulation.

Applying differential network analysis to longitudinal gene expression in response to perturbations.

Differential Network (DN) analysis is a method that has long been used to interpret changes in gene expression data and provide biological insights. The method identifies the rewiring of gene networks in response to external perturbations. Our study applies the DN method to the analysis of RNA-sequencing (RNA-seq) time series datasets. We focus on expression changes: (i) in saliva of a human subject after pneumococcal vaccination (PPSV23) and (ii) in primary B cells treated ex vivo with a monoclonal antibody drug (Rituximab). The DN method enabled us to identify the activation of biological pathways consistent with the mechanisms of action of the PPSV23 vaccine and target pathways of Rituximab. The community detection algorithm on the DN revealed clusters of genes characterized by collective temporal behavior. All saliva and some B cell DN communities showed characteristic time signatures, outlining a chronological order in pathway activation in response to the perturbation. Moreover, we identified early and delayed responses within network modules in the saliva dataset and three temporal patterns in the B cell data.

Temporal response characterization across individual multiomics profiles of prediabetic and diabetic subjects.

Longitudinal deep multiomics profiling, which combines biomolecular, physiological, environmental and clinical measures data, shows great promise for precision health. However, integrating and understanding the complexity of such data remains a big challenge. Here we utilize an individual-focused bottom-up approach aimed at first assessing single individuals' multiomics time series, and using the individual-level responses to assess multi-individual grouping based directly on similarity of their longitudinal deep multiomics profiles. We used this individual-focused approach to analyze profiles from a study profiling longitudinal responses in type 2 diabetes mellitus. After generating periodograms for individual subject omics signals, we constructed within-person omics networks and analyzed personal-level immune changes. The results identified both individual-level responses to immune perturbation, and the clusters of individuals that have similar behaviors in immune response and which were associated to measures of their diabetic status.

Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells.

VEGF inhibitor drugs are part of standard care in oncology and ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them. Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper. The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals.

LIF, a mitogen for choroidal endothelial cells, protects the choriocapillaris: implications for prevention of geographic atrophy.

In the course of our studies aiming to discover vascular bed-specific endothelial cell (EC) mitogens, we identified leukemia inhibitory factor (LIF) as a mitogen for bovine choroidal EC (BCE), although LIF has been mainly characterized as an EC growth inhibitor and an anti-angiogenic molecule. LIF stimulated growth of BCE while it inhibited, as previously reported, bovine aortic EC (BAE) growth. The JAK-STAT3 pathway mediated LIF actions in both BCE and BAE cells, but a caspase-independent proapoptotic signal mediated by cathepsins was triggered in BAE but not in BCE. LIF administration directly promoted activation of STAT3 and increased blood vessel density in mouse eyes. LIF also had protective effects on the choriocapillaris in a model of oxidative retinal injury. Analysis of available single-cell transcriptomic datasets shows strong expression of the specific LIF receptor in mouse and human choroidal EC. Our data suggest that LIF administration may be an innovative approach to prevent atrophy associated with AMD, through protection of the choriocapillaris.

Visibility graph based temporal community detection with applications in biological time series.

Temporal behavior is an essential aspect of all biological systems. Time series have been previously represented as networks. Such representations must address two fundamental problems on how to: (1) Create appropriate networks to reflect the characteristics of biological time series. (2) Detect characteristic dynamic patterns or events as network temporal communities. General community detection methods use metrics comparing the connectivity within a community to random models, or are based on the betweenness centrality of edges or nodes. However, such methods were not designed for network representations of time series. We introduce a visibility-graph-based method to build networks from time series and detect temporal communities within these networks. To characterize unevenly sampled time series (typical of biological experiments), and simultaneously capture events associated to peaks and troughs, we introduce the Weighted Dual-Perspective Visibility Graph (WDPVG). To detect temporal communities in individual signals, we first find the shortest path of the network between start and end nodes, identifying high intensity nodes as the main stem of our community detection algorithm that act as hubs for each community. Then, we aggregate nodes outside the shortest path to the closest nodes found on the main stem based on the closest path length, thereby assigning every node to a temporal community based on proximity to the stem nodes/hubs. We demonstrate the validity and effectiveness of our method through simulation and biological applications.

Digital Cell Sorter (DCS): a cell type identification, anomaly detection, and Hopfield landscapes toolkit for single-cell transcriptomics.

MOTIVATION: Analysis of singe cell RNA sequencing (scRNA-seq) typically consists of different steps including quality control, batch correction, clustering, cell identification and characterization, and visualization. The amount of scRNA-seq data is growing extremely fast, and novel algorithmic approaches improving these steps are key to extract more biological information. Here, we introduce: (i) two methods for automatic cell type identification (i.e., without expert curator) based on a voting algorithm and a Hopfield classifier, (ii) a method for cell anomaly quantification based on isolation forest, and (iii) a tool for the visualization of cell phenotypic landscapes based on Hopfield energy-like functions. These new approaches are integrated in a software platform that includes many other state-of-the-art methodologies and provides a self-contained toolkit for scRNA-seq analysis. RESULTS: We present a suite of software elements for the analysis of scRNA-seq data. This Python-based open source software, Digital Cell Sorter (DCS), consists in an extensive toolkit of methods for scRNA-seq analysis. We illustrate the capability of the software using data from large datasets of peripheral blood mononuclear cells (PBMC), as well as plasma cells of bone marrow samples from healthy donors and multiple myeloma patients. We test the novel algorithms by evaluating their ability to deconvolve cell mixtures and detect small numbers of anomalous cells in PBMC data. AVAILABILITY: The DCS toolkit is available for download and installation through the Python Package Index (PyPI). The software can be deployed using the Python import function following installation. Source code is also available for download on Zenodo: DOI 10.5281/zenodo.2533377. SUPPLEMENTARY INFORMATION: Supplemental Materials are available at PeerJ online.

Longitudinal saliva omics responses to immune perturbation: a case study.

Saliva omics has immense potential for non-invasive diagnostics, including monitoring very young or elderly populations, or individuals in remote locations. In this study, multiple saliva omics from an individual were monitored over three periods (100 timepoints) involving: (1) hourly sampling over 24 h without intervention, (2) hourly sampling over 24 h including immune system activation using the standard 23-valent pneumococcal polysaccharide vaccine, (3) daily sampling for 33 days profiling the post-vaccination response. At each timepoint total saliva transcriptome and proteome, and small RNA from salivary extracellular vesicles were profiled, including mRNA, miRNA, piRNA and bacterial RNA. The two 24-h periods were used in a paired analysis to remove daily variation and reveal vaccination responses. Over 18,000 omics longitudinal series had statistically significant temporal trends compared to a healthy baseline. Various immune response and regulation pathways were activated following vaccination, including interferon and cytokine signaling, and MHC antigen presentation. Immune response timeframes were concordant with innate and adaptive immunity development, and coincided with vaccination and reported fever. Overall, mRNA results appeared more specific and sensitive (timewise) to vaccination compared to other omics. The results suggest saliva omics can be consistently assessed for non-invasive personalized monitoring and immune response diagnostics.

PyIOmica: longitudinal omics analysis and trend identification.

SUMMARY: PyIOmica is an open-source Python package focusing on integrating longitudinal multiple omics datasets, characterizing and categorizing temporal trends. The package includes multiple bioinformatics tools including data normalization, annotation, categorization, visualization and enrichment analysis for gene ontology terms and pathways. Additionally, the package includes an implementation of visibility graphs to visualize time series as networks. AVAILABILITY AND IMPLEMENTATION: PyIOmica is implemented as a Python package (pyiomica), available for download and installation through the Python Package Index (https://pypi.python.org/pypi/pyiomica), and can be deployed using the Python import function following installation. PyIOmica has been tested on Mac OS X, Unix/Linux and Microsoft Windows. The application is distributed under an MIT license. Source code for each release is also available for download on Zenodo (https://doi.org/10.5281/zenodo.3548040). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics.

Polled Digital Cell Sorter (p-DCS): Automatic identification of hematological cell types from single cell RNA-sequencing clusters.

BACKGROUND: Single cell RNA sequencing (scRNA-seq) brings unprecedented opportunities for mapping the heterogeneity of complex cellular environments such as bone marrow, and provides insight into many cellular processes. Single cell RNA-seq has a far larger fraction of missing data reported as zeros (dropouts) than traditional bulk RNA-seq, and unsupervised clustering combined with Principal Component Analysis (PCA) can be used to overcome this limitation. After clustering, however, one has to interpret the average expression of markers on each cluster to identify the corresponding cell types, and this is normally done by hand by an expert curator. RESULTS: We present a computational tool for processing single cell RNA-seq data that uses a voting algorithm to automatically identify cells based on approval votes received by known molecular markers. Using a stochastic procedure that accounts for imbalances in the number of known molecular signatures for different cell types, the method computes the statistical significance of the final approval score and automatically assigns a cell type to clusters without an expert curator. We demonstrate the utility of the tool in the analysis of eight samples of bone marrow from the Human Cell Atlas. The tool provides a systematic identification of cell types in bone marrow based on a list of markers of immune cell types, and incorporates a suite of visualization tools that can be overlaid on a t-SNE representation. The software is freely available as a Python package at https://github.com/sdomanskyi/DigitalCellSorter . CONCLUSIONS: This methodology assures that extensive marker to cell type matching information is taken into account in a systematic way when assigning cell clusters to cell types. Moreover, the method allows for a high throughput processing of multiple scRNA-seq datasets, since it does not involve an expert curator, and it can be applied recursively to obtain cell sub-types. The software is designed to allow the user to substitute the marker to cell type matching information and apply the methodology to different cellular environments.

Modeling disease progression in Multiple Myeloma with Hopfield networks and single-cell RNA-seq.

Associative memories in Hopfield's neural networks are mapped to gene expression pattern to model different paths of disease progression towards Multiple Myeloma (MM). The model is built using single cell RNA-seq data from bone marrow aspirates of MM patients as well as patients diagnosed with Monoclonal Gammopathy of Undetermined Significance (MGUS) and Smoldering Multiple Myeloma (SMM), two medical conditions that often progress to full MM. Results: We identify different clusters of MGUS, SMM, and MM cells, map them to Hopfield associative memory patterns, and model the dynamics of transition between the different patterns. The model is then used to identify genes that are differentialy expressed across different MM stages and whose simultaneous inhibition is associated to a delayed disease progression.

Cell cycle time series gene expression data encoded as cyclic attractors in Hopfield systems.

Modern time series gene expression and other omics data sets have enabled unprecedented resolution of the dynamics of cellular processes such as cell cycle and response to pharmaceutical compounds. In anticipation of the proliferation of time series data sets in the near future, we use the Hopfield model, a recurrent neural network based on spin glasses, to model the dynamics of cell cycle in HeLa (human cervical cancer) and S. cerevisiae cells. We study some of the rich dynamical properties of these cyclic Hopfield systems, including the ability of populations of simulated cells to recreate experimental expression data and the effects of noise on the dynamics. Next, we use a genetic algorithm to identify sets of genes which, when selectively inhibited by local external fields representing gene silencing compounds such as kinase inhibitors, disrupt the encoded cell cycle. We find, for example, that inhibiting the set of four kinases AURKB, NEK1, TTK, and WEE1 causes simulated HeLa cells to accumulate in the M phase. Finally, we suggest possible improvements and extensions to our model.

High-throughput screening and bioinformatic analysis to ascertain compounds that prevent saturated fatty acid-induced ß-cell apoptosis.

Pancreatic ß-cell lipotoxicity is a central feature of the pathogenesis of type 2 diabetes. To study the mechanism by which fatty acids cause ß-cell death and develop novel approaches to prevent it, a high-throughput screen on the ß-cell line INS1 was carried out. The cells were exposed to palmitate to induce cell death and compounds that reversed palmitate-induced cytotoxicity were ascertained. Hits from the screen were analyzed by an increasingly more stringent testing funnel, ending with studies on primary human islets treated with palmitate. MAP4K4 inhibitors, which were not part of the screening libraries but were ascertained by a bioinformatics analysis, and the endocannabinoid anandamide were effective at inhibiting palmitate-induced apoptosis in INS1 cells as well as primary rat and human islets. These targets could serve as the starting point for the development of therapeutics for type 2 diabetes.

Evolutionary and Topological Properties of Genes and Community Structures in Human Gene Regulatory Networks.

The diverse, specialized genes present in today's lifeforms evolved from a common core of ancient, elementary genes. However, these genes did not evolve individually: gene expression is controlled by a complex network of interactions, and alterations in one gene may drive reciprocal changes in its proteins' binding partners. Like many complex networks, these gene regulatory networks (GRNs) are composed of communities, or clusters of genes with relatively high connectivity. A deep understanding of the relationship between the evolutionary history of single genes and the topological properties of the underlying GRN is integral to evolutionary genetics. Here, we show that the topological properties of an acute myeloid leukemia GRN and a general human GRN are strongly coupled with its genes' evolutionary properties. Slowly evolving ("cold"), old genes tend to interact with each other, as do rapidly evolving ("hot"), young genes. This naturally causes genes to segregate into community structures with relatively homogeneous evolutionary histories. We argue that gene duplication placed old, cold genes and communities at the center of the networks, and young, hot genes and communities at the periphery. We demonstrate this with single-node centrality measures and two new measures of efficiency, the set efficiency and the interset efficiency. We conclude that these methods for studying the relationships between a GRN's community structures and its genes' evolutionary properties provide new perspectives for understanding evolutionary genetics.

Computational dynamics of acoustically driven microsphere systems.

We propose a computational framework for the self-consistent dynamics of a microsphere system driven by a pulsed acoustic field in an ideal fluid. Our framework combines a molecular dynamics integrator describing the dynamics of the microsphere system with a time-dependent integral equation solver for the acoustic field that makes use of fields represented as surface expansions in spherical harmonic basis functions. The presented approach allows us to describe the interparticle interaction induced by the field as well as the dynamics of trapping in counter-propagating acoustic pulses. The integral equation formulation leads to equations of motion for the microspheres describing the effect of nondissipative drag forces. We show (1) that the field-induced interactions between the microspheres give rise to effective dipolar interactions, with effective dipoles defined by their velocities and (2) that the dominant effect of an ultrasound pulse through a cloud of microspheres gives rise mainly to a translation of the system, though we also observe both expansion and contraction of the cloud determined by the initial system geometry.

Combinations of Kinase Inhibitors Protecting Myoblasts against Hypoxia.

Cell-based therapies to treat skeletal muscle disease are limited by the poor survival of donor myoblasts, due in part to acute hypoxic stress. After confirming that the microenvironment of transplanted myoblasts is hypoxic, we screened a kinase inhibitor library in vitro and identified five kinase inhibitors that protected myoblasts from cell death or growth arrest in hypoxic conditions. A systematic, combinatorial study of these compounds further improved myoblast viability, showing both synergistic and additive effects. Pathway and target analysis revealed CDK5, CDK2, CDC2, WEE1, and GSK3ß as the main target kinases. In particular, CDK5 was the center of the target kinase network. Using our recently developed statistical method based on elastic net regression we computationally validated the key role of CDK5 in cell protection against hypoxia. This method provided a list of potential kinase targets with a quantitative measure of their optimal amount of relative inhibition. A modified version of the method was also able to predict the effect of combinations using single-drug response data. This work is the first step towards a broadly applicable system-level strategy for the pharmacology of hypoxic damage.

A scalable method for molecular network reconstruction identifies properties of targets and mutations in acute myeloid leukemia.

A key aim of systems biology is the reconstruction of molecular networks. We do not yet, however, have networks that integrate information from all datasets available for a particular clinical condition. This is in part due to the limited scalability, in terms of required computational time and power, of existing algorithms. Network reconstruction methods should also be scalable in the sense of allowing scientists from different backgrounds to efficiently integrate additional data. We present a network model of acute myeloid leukemia (AML). In the current version (AML 2.1), we have used gene expression data (both microarray and RNA-seq) from 5 different studies comprising a total of 771 AML samples and a protein-protein interactions dataset. Our scalable network reconstruction method is in part based on the well-known property of gene expression correlation among interacting molecules. The difficulty of distinguishing between direct and indirect interactions is addressed by optimizing the coefficient of variation of gene expression, using a validated gold-standard dataset of direct interactions. Computational time is much reduced compared to other network reconstruction methods. A key feature is the study of the reproducibility of interactions found in independent clinical datasets. An analysis of the most significant clusters, and of the network properties (intraset efficiency, degree, betweenness centrality, and PageRank) of common AML mutations demonstrated the biological significance of the network. A statistical analysis of the response of blast cells from 11 AML patients to a library of kinase inhibitors provided an experimental validation of the network. A combination of network and experimental data identified CDK1, CDK2, CDK4, and CDK6 and other kinases as potential therapeutic targets in AML.

Control of asymmetric Hopfield networks and application to cancer attractors.

The asymmetric Hopfield model is used to simulate signaling dynamics in gene regulatory networks. The model allows for a direct mapping of a gene expression pattern into attractor states. We analyze different control strategies aimed at disrupting attractor patterns using selective local fields representing therapeutic interventions. The control strategies are based on the identification of signaling bottlenecks, which are single nodes or strongly connected clusters of nodes that have a large impact on the signaling. We provide a theorem with bounds on the minimum number of nodes that guarantee control of bottlenecks consisting of strongly connected components. The control strategies are applied to the identification of sets of proteins that, when inhibited, selectively disrupt the signaling of cancer cells while preserving the signaling of normal cells. We use an experimentally validated non-specific and an algorithmically-assembled specific B cell gene regulatory network reconstructed from gene expression data to model cancer signaling in lung and B cells, respectively. Among the potential targets identified here are TP53, FOXM1, BCL6 and SRC. This model could help in the rational design of novel robust therapeutic interventions based on our increasing knowledge of complex gene signaling networks.

Identification of drug combinations containing imatinib for treatment of BCR-ABL+ leukemias.

The BCR-ABL translocation is found in chronic myeloid leukemia (CML) and in Ph+ acute lymphoblastic leukemia (ALL) patients. Although imatinib and its analogues have been used as front-line therapy to target this mutation and control the disease for over a decade, resistance to the therapy is still observed and most patients are not cured but need to continue the therapy indefinitely. It is therefore of great importance to find new therapies, possibly as drug combinations, which can overcome drug resistance. In this study, we identified eleven candidate anti-leukemic drugs that might be combined with imatinib, using three approaches: a kinase inhibitor library screen, a gene expression correlation analysis, and literature analysis. We then used an experimental search algorithm to efficiently explore the large space of possible drug and dose combinations and identified drug combinations that selectively kill a BCR-ABL+ leukemic cell line (K562) over a normal fibroblast cell line (IMR-90). Only six iterations of the algorithm were needed to identify very selective drug combinations. The efficacy of the top forty-nine combinations was further confirmed using Ph+ and Ph- ALL patient cells, including imatinib-resistant cells. Collectively, the drug combinations and methods we describe might be a first step towards more effective interventions for leukemia patients, especially those with the BCR-ABL translocation.