ABSTRACT
The ability to quantify structural changes of the endoplasmic reticulum (ER) is crucial for understanding the structure and function of this organelle. However, the rapid movement and complex topology of ER networks make this challenging. Here, we construct a state-of-the-art semantic segmentation method that we call ERnet for the automatic classification of sheet and tubular ER domains inside individual cells. Data are skeletonized and represented by connectivity graphs, enabling precise and efficient quantification of network connectivity. ERnet generates metrics on topology and integrity of ER structures and quantifies structural change in response to genetic or metabolic manipulation. We validate ERnet using data obtained by various ER-imaging methods from different cell types as well as ground truth images of synthetic ER structures. ERnet can be deployed in an automatic high-throughput and unbiased fashion and identifies subtle changes in ER phenotypes that may inform on disease progression and response to therapy.
Subject(s)
Endoplasmic Reticulum , Semantics , Endoplasmic Reticulum/metabolismABSTRACT
Protein design and directed evolution have separately contributed enormously to protein engineering. Without being mutually exclusive, the former relies on computation from first principles, while the latter is a combinatorial approach based on chance. Advances in ultrahigh throughput (uHT) screening, next generation sequencing and machine learning may create alternative routes to engineered proteins, where functional information linked to specific sequences is interpreted and extrapolated in silico. In particular, the miniaturisation of functional tests in water-in-oil emulsion droplets with picoliter volumes and their rapid generation and analysis (>1 kHz) allows screening of >107-membered libraries in a day. Subsequently, decoding the selected clones by short or long-read sequencing methods leads to large sequence-function datasets that may allow extrapolation from experimental directed evolution to further improved mutants beyond the observed hits. In this work, we explore experimental strategies for how to draw up 'fitness landscapes' in sequence space with uHT droplet microfluidics, review the current state of AI/ML in enzyme engineering and discuss how uHT datasets may be combined with AI/ML to make meaningful predictions and accelerate biocatalyst engineering.
ABSTRACT
The HIV-1 matrix protein p17 (p17) is a pleiotropic molecule impacting on different cell types. Its interaction with many cellular proteins underlines the importance of the viral protein as a major determinant of human specific adaptation. We previously showed the proangiogenic capability of p17. Here, by integrating functional analysis and receptor binding, we identify a functional epitope that displays molecular mimicry with human erythropoietin (EPO) and promotes angiogenesis through common beta chain receptor (ßCR) activation. The functional EPO-like epitope was found to be present in the matrix protein of HIV-1 ancestors SIV originated in chimpanzees (SIVcpz) and gorillas (SIVgor) but not in that of HIV-2 and its ancestor SIVsmm from sooty mangabeys. According to biological data, evolution of the EPO-like epitope showed a clear differentiation between HIV-1/SIVcpz-gor and HIV-2/SIVsmm branches, thus highlighting this epitope on p17 as a divergent signature discriminating HIV-1 and HIV-2 ancestors. P17 is known to enhance HIV-1 replication. Similarly to other ßCR ligands, p17 is capable of attracting and activating HIV-1 target cells and promoting a proinflammatory microenvironment. Thus, it is tempting to speculate that acquisition of an epitope on the matrix proteins of HIV-1 ancestors capable of triggering ßCR may have represented a critical step to enhance viral aggressiveness and early human-to-human SIVcpz/gor dissemination. The hypothesis that the p17/ßCR interaction and ßCR abnormal stimulation may also play a role in sustaining chronic activation and inflammation, thus marking the difference between HIV-1 and HIV-2 in term of pathogenicity, needs further investigation.
Subject(s)
Erythropoietin/genetics , HIV Antigens/metabolism , HIV-1/metabolism , gag Gene Products, Human Immunodeficiency Virus/metabolism , Cells, Cultured , Epitopes/immunology , Erythropoietin/metabolism , Evolution, Molecular , HIV Antigens/genetics , HIV Seropositivity , HIV-1/genetics , HIV-2 , Humans , Molecular Mimicry , Simian Immunodeficiency Virus , gag Gene Products, Human Immunodeficiency Virus/geneticsABSTRACT
Although chemotherapy is the first-line treatment for ovarian cancer (OCa) patients, chemoresistance (CR) decreases their progression-free survival. This paper investigates the genetic interaction (GI) related to OCa-CR. To decrease the complexity of establishing gene networks, individual signature genes related to OCa-CR are identified using a gradient boosting decision tree algorithm. Additionally, the genetic interaction coefficient (GIC) is proposed to measure the correlation of two signature genes quantitatively and explain their joint influence on OCa-CR. Gene pair that possesses high GIC is identified as signature pair. A total of 24 signature gene pairs are selected that include 10 individual signature genes and the influence of signature gene pairs on OCa-CR is explored. Finally, a signature gene pair-based prediction of OCa-CR is identified. The area under curve (AUC) is a widely used performance measure for machine learning prediction. The AUC of signature gene pair reaches 0.9658, whereas the AUC of individual signature gene-based prediction is 0.6823 only. The identified signature gene pairs not only build an efficient GI network of OCa-CR but also provide an interesting way for OCa-CR prediction. This improvement shows that our proposed method is a useful tool to investigate GI related to OCa-CR.
Subject(s)
Databases, Nucleic Acid , Drug Resistance, Neoplasm/genetics , Gene Expression Profiling , Gene Expression Regulation, Neoplastic , Machine Learning , Ovarian Neoplasms , Female , Gene Regulatory Networks , Humans , Ovarian Neoplasms/genetics , Ovarian Neoplasms/metabolismABSTRACT
The novel coronavirus (2019-nCoV) has recently emerged, causing COVID-19 outbreaks and significant societal/global disruption. Importantly, COVID-19 infection resembles SARS-like complications. However, the lack of knowledge about the underlying genetic mechanisms of COVID-19 warrants the development of prospective control measures. In this study, we employed whole-genome alignment and digital DNA-DNA hybridization analyses to assess genomic linkage between 2019-nCoV and other coronaviruses. To understand the pathogenetic behavior of 2019-nCoV, we compared gene expression datasets of viral infections closest to 2019-nCoV with four COVID-19 clinical presentations followed by functional enrichment of shared dysregulated genes. Potential chemical antagonists were also identified using protein-chemical interaction analysis. Based on phylogram analysis, the 2019-nCoV was found genetically closest to SARS-CoVs. In addition, we identified 562 upregulated and 738 downregulated genes (adj. P ≤ 0.05) with SARS-CoV infection. Among the dysregulated genes, SARS-CoV shared ≤19 upregulated and ≤22 downregulated genes with each of different COVID-19 complications. Notably, upregulation of BCL6 and PFKFB3 genes was common to SARS-CoV, pneumonia and severe acute respiratory syndrome, while they shared CRIP2, NSG1 and TNFRSF21 genes in downregulation. Besides, 14 genes were common to different SARS-CoV comorbidities that might influence COVID-19 disease. We also observed similarities in pathways that can lead to COVID-19 and SARS-CoV diseases. Finally, protein-chemical interactions suggest cyclosporine, resveratrol and quercetin as promising drug candidates against COVID-19 as well as other SARS-like viral infections. The pathogenetic analyses, along with identified biomarkers, signaling pathways and chemical antagonists, could prove useful for novel drug development in the fight against the current global 2019-nCoV pandemic.
Subject(s)
COVID-19/virology , SARS-CoV-2/pathogenicity , Severe acute respiratory syndrome-related coronavirus/pathogenicity , Antiviral Agents/therapeutic use , COVID-19/complications , Case-Control Studies , Comorbidity , Genome, Viral , Humans , MicroRNAs/metabolism , Severe acute respiratory syndrome-related coronavirus/genetics , Transcription Factors/metabolism , COVID-19 Drug TreatmentABSTRACT
MOTIVATION: High-throughput gene expression can be used to address a wide range of fundamental biological problems, but datasets of an appropriate size are often unavailable. Moreover, existing transcriptomics simulators have been criticized because they fail to emulate key properties of gene expression data. In this article, we develop a method based on a conditional generative adversarial network to generate realistic transcriptomics data for Escherichia coli and humans. We assess the performance of our approach across several tissues and cancer-types. RESULTS: We show that our model preserves several gene expression properties significantly better than widely used simulators, such as SynTReN or GeneNetWeaver. The synthetic data preserve tissue- and cancer-specific properties of transcriptomics data. Moreover, it exhibits real gene clusters and ontologies both at local and global scales, suggesting that the model learns to approximate the gene expression manifold in a biologically meaningful way. AVAILABILITY AND IMPLEMENTATION: Code is available at: https://github.com/rvinas/adversarial-gene-expression. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Escherichia coli , Gene Expression Profiling , Humans , Gene Expression Profiling/methods , Gene ExpressionABSTRACT
MOTIVATION: Single-cell RNA sequencing allows high-resolution views of individual cells for libraries of up to millions of samples, thus motivating the use of deep learning for analysis. In this study, we introduce the use of graph neural networks for the unsupervised exploration of scRNA-seq data by developing a variational graph autoencoder architecture with graph attention layers that operates directly on the connectivity between cells, focusing on dimensionality reduction and clustering. With the help of several case studies, we show that our model, named CellVGAE, can be effectively used for exploratory analysis even on challenging datasets, by extracting meaningful features from the data and providing the means to visualize and interpret different aspects of the model. RESULTS: We show that CellVGAE is more interpretable than existing scRNA-seq variational architectures by analysing the graph attention coefficients. By drawing parallels with other scRNA-seq studies on interpretability, we assess the validity of the relationships modelled by attention, and furthermore, we show that CellVGAE can intrinsically capture information such as pseudotime and NF-ĸB activation dynamics, the latter being a property that is not generally shared by existing neural alternatives. We then evaluate the dimensionality reduction and clustering performance on 9 difficult and well-annotated datasets by comparing with three leading neural and non-neural techniques, concluding that CellVGAE outperforms competing methods. Finally, we report a decrease in training times of up to × 20 on a dataset of 1.3 million cells compared to existing deep learning architectures. AVAILABILITYAND IMPLEMENTATION: The CellVGAE code is available at https://github.com/davidbuterez/CellVGAE. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Gene Expression Profiling , Single-Cell Gene Expression Analysis , Gene Expression Profiling/methods , Sequence Analysis, RNA/methods , Workflow , Single-Cell Analysis/methods , Cluster AnalysisABSTRACT
MOTIVATION: Gene expression data are commonly used at the intersection of cancer research and machine learning for better understanding of the molecular status of tumour tissue. Deep learning predictive models have been employed for gene expression data due to their ability to scale and remove the need for manual feature engineering. However, gene expression data are often very high dimensional, noisy and presented with a low number of samples. This poses significant problems for learning algorithms: models often overfit, learn noise and struggle to capture biologically relevant information. In this article, we utilize external biological knowledge embedded within structures of gene interaction graphs such as protein-protein interaction (PPI) networks to guide the construction of predictive models. RESULTS: We present Gene Interaction Network Constrained Construction (GINCCo), an unsupervised method for automated construction of computational graph models for gene expression data that are structurally constrained by prior knowledge of gene interaction networks. We employ this methodology in a case study on incorporating a PPI network in cancer phenotype prediction tasks. Our computational graphs are structurally constructed using topological clustering algorithms on the PPI networks which incorporate inductive biases stemming from network biology research on protein complex discovery. Each of the entities in the GINCCo computational graph represents biological entities such as genes, candidate protein complexes and phenotypes instead of arbitrary hidden nodes of a neural network. This provides a biologically relevant mechanism for model regularization yielding strong predictive performance while drastically reducing the number of model parameters and enabling guided post-hoc enrichment analyses of influential gene sets with respect to target phenotypes. Our experiments analysing a variety of cancer phenotypes show that GINCCo often outperforms support vector machine, Fully Connected Multi-layer Perceptrons (MLP) and Randomly Connected MLPs despite greatly reduced model complexity. AVAILABILITY AND IMPLEMENTATION: https://github.com/paulmorio/gincco contains the source code for our approach. We also release a library with algorithms for protein complex discovery within PPI networks at https://github.com/paulmorio/protclus. This repository contains implementations of the clustering algorithms used in this article. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Algorithms , Neoplasms , Humans , Neural Networks, Computer , Software , Neoplasms/genetics , Bias , Gene Expression , Computational Biology/methodsABSTRACT
High-throughput screening (HTS), as one of the key techniques in drug discovery, is frequently used to identify promising drug candidates in a largely automated and cost-effective way. One of the necessary conditions for successful HTS campaigns is a large and diverse compound library, enabling hundreds of thousands of activity measurements per project. Such collections of data hold great promise for computational and experimental drug discovery efforts, especially when leveraged in combination with modern deep learning techniques, and can potentially lead to improved drug activity predictions and cheaper and more effective experimental design. However, existing collections of machine-learning-ready public datasets do not exploit the multiple data modalities present in real-world HTS projects. Thus, the largest fraction of experimental measurements, corresponding to hundreds of thousands of "noisy" activity values from primary screening, are effectively ignored in the majority of machine learning models of HTS data. To address these limitations, we introduce Multifidelity PubChem BioAssay (MF-PCBA), a curated collection of 60 datasets that includes two data modalities for each dataset, corresponding to primary and confirmatory screening, an aspect that we call multifidelity. Multifidelity data accurately reflect real-world HTS conventions and present a new, challenging task for machine learning: the integration of low- and high-fidelity measurements through molecular representation learning, taking into account the orders-of-magnitude difference in size between the primary and confirmatory screens. Here we detail the steps taken to assemble MF-PCBA in terms of data acquisition from PubChem and the filtering steps required to curate the raw data. We also provide an evaluation of a recent deep-learning-based method for multifidelity integration across the introduced datasets, demonstrating the benefit of leveraging all HTS modalities, and a discussion in terms of the roughness of the molecular activity landscape. In total, MF-PCBA contains over 16.6 million unique molecule-protein interactions. The datasets can be easily assembled by using the source code available at https://github.com/davidbuterez/mf-pcba.
Subject(s)
Benchmarking , High-Throughput Screening Assays , High-Throughput Screening Assays/methods , Drug Discovery/methods , Machine Learning , Biological AssayABSTRACT
The development of efficient and effective bioinformatics tools and pipelines for identifying peptides with dipeptidyl peptidase IV (DPP-IV) inhibitory activities from large-scale protein datasets is of great importance for the discovery and development of potential and promising antidiabetic drugs. In this study, we present a novel stacking-based ensemble learning predictor (termed StackDPPIV) designed for identification of DPP-IV inhibitory peptides. Unlike the existing method, which is based on single-feature-based methods, we combined five popular machine learning algorithms in conjunction with ten different feature encodings from multiple perspectives to generate a pool of various baseline models. Subsequently, the probabilistic features derived from these baseline models were systematically integrated and deemed as new feature representations. Finally, in order to improve the predictive performance, the genetic algorithm based on the self-assessment-report was utilized to determine a set of informative probabilistic features and then used the optimal one for developing the final meta-predictor (StackDPPIV). Experiment results demonstrated that StackDPPIV could outperform its constituent baseline models on both the training and independent datasets. Furthermore, StackDPPIV achieved an accuracy of 0.891, MCC of 0.784 and AUC of 0.961, which were 9.4%, 19.0% and 11.4%, respectively, higher than that of the existing method on the independent test. Feature analysis demonstrated that our feature representations had more discriminative ability as compared to conventional feature descriptors, which highlights the combination of different features was essential for the performance improvement. In order to implement the proposed predictor, we had built a user-friendly online web server at http://pmlabstack.pythonanywhere.com/StackDPPIV.
Subject(s)
Dipeptidyl Peptidase 4 , Peptides , Computational Biology , Dipeptidyl Peptidase 4/metabolism , Machine Learning , Peptides/pharmacology , ProteinsABSTRACT
Bone marrow hematopoietic stem cells (HSCs) are crucial to maintain lifelong production of all blood cells. Although HSCs divide infrequently, it is thought that the entire HSC pool turns over every few weeks, suggesting that HSCs regularly enter and exit cell cycle. Here, we combine flow cytometry with label-retaining assays (BrdU and histone H2B-GFP) to identify a population of dormant mouse HSCs (d-HSCs) within the lin(-)Sca1+cKit+CD150+CD48(-)CD34(-) population. Computational modeling suggests that d-HSCs divide about every 145 days, or five times per lifetime. d-HSCs harbor the vast majority of multilineage long-term self-renewal activity. While they form a silent reservoir of the most potent HSCs during homeostasis, they are efficiently activated to self-renew in response to bone marrow injury or G-CSF stimulation. After re-establishment of homeostasis, activated HSCs return to dormancy, suggesting that HSCs are not stochastically entering the cell cycle but reversibly switch from dormancy to self-renewal under conditions of hematopoietic stress.
Subject(s)
Adult Stem Cells/cytology , Hematopoietic Stem Cells/cytology , Adult Stem Cells/physiology , Animals , Antigens, Differentiation/metabolism , Bone Marrow/physiology , Bromouracil/analogs & derivatives , Fluorouracil/metabolism , Green Fluorescent Proteins , Hematopoietic Stem Cells/physiology , Homeostasis , Mice , Mice, Transgenic , Uridine/analogs & derivatives , Uridine/metabolismABSTRACT
Although still in its infancy, artificial intelligence (AI) analysis of kidney biopsy images is anticipated to become an integral aspect of renal histopathology. As these systems are developed, the focus will understandably be on developing ever more accurate models, but successful translation to the clinic will also depend upon other characteristics of the system.In the extreme, deployment of highly performant but "black box" AI is fraught with risk, and high-profile errors could damage future trust in the technology. Furthermore, a major factor determining whether new systems are adopted in clinical settings is whether they are "trusted" by clinicians. Key to unlocking trust will be designing platforms optimized for intuitive human-AI interactions and ensuring that, where judgment is required to resolve ambiguous areas of assessment, the workings of the AI image classifier are understandable to the human observer. Therefore, determining the optimal design for AI systems depends on factors beyond performance, with considerations of goals, interpretability, and safety constraining many design and engineering choices.In this article, we explore challenges that arise in the application of AI to renal histopathology, and consider areas where choices around model architecture, training strategy, and workflow design may be influenced by factors beyond the final performance metrics of the system.
Subject(s)
Artificial Intelligence , Trust , Humans , KidneyABSTRACT
Coeliac disease (CD) is a complex, multifactorial pathology caused by different factors, such as nutrition, immunological response and genetic factors. Many autoimmune diseases are comorbidities for CD, and a comprehensive and integrated analysis with bioinformatics approaches can help in evaluating the interconnections among all the selected pathologies. We first performed a detailed survey of gene expression data available in public repositories on CD and less commonly considered comorbidities. Then we developed an innovative pipeline that integrates gene expression, cell-type data and online resources (e.g. a list of comorbidities from the literature), using bioinformatics methods such as gene set enrichment analysis and semantic similarity. Our pipeline is written in R language, available at the following link: http://bioinformatica.isa.cnr.it/COELIAC_DISEASE/SCRIPTS/. We found a list of common differential expressed genes, gene ontology terms and pathways among CD and comorbidities and the closeness among the selected pathologies by means of disease ontology terms. Physicians and other researchers, such as molecular biologists, systems biologists and pharmacologists can use it to analyze pathology in detail, from differential expressed genes to ontologies, performing a comparison with the pathology comorbidities or with other diseases.
ABSTRACT
MOTIVATION: One of the branches of Systems Biology is focused on a deep understanding of underlying regulatory networks through the analysis of the biomolecules oscillations and their interplay. Synthetic Biology exploits gene or/and protein regulatory networks towards the design of oscillatory networks for producing useful compounds. Therefore, at different levels of application and for different purposes, the study of biomolecular oscillations can lead to different clues about the mechanisms underlying living cells. It is known that network-level interactions involve more than one type of biomolecule as well as biological processes operating at multiple omic levels. Combining network/pathway-level information with genetic information it is possible to describe well-understood or unknown bacterial mechanisms and organism-specific dynamics. RESULTS: Following the methodologies used in signal processing and communication engineering, a methodology is introduced to identify and quantify the extent of multi-omic oscillations. These are due to the process of multi-omic integration and depend on the gene positions on the chromosome. Ad hoc signal metrics are designed to allow further biotechnological explanations and provide important clues about the oscillatory nature of the pathways and their regulatory circuits. Our algorithms designed for the analysis of multi-omic signals are tested and validated on 11 different bacteria for thousands of multi-omic signals perturbed at the network level by different experimental conditions. Information on the order of genes, codon usage, gene expression and protein molecular weight is integrated at three different functional levels. Oscillations show interesting evidence that network-level multi-omic signals present a synchronized response to perturbations and evolutionary relations along taxa. AVAILABILITY AND IMPLEMENTATION: The algorithms, the code (in language R), the tool, the pipeline and the whole dataset of multi-omic signal metrics are available at: https://github.com/lodeguns/Multi-omicSignals. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Algorithms , Benchmarking , Bacteria/genetics , Gene Regulatory Networks , Systems BiologyABSTRACT
The blood-brain barrier (BBB) is the primary barrier with a highly selective semipermeable border between blood vascular endothelial cells and the central nervous system. Since BBB can prevent drugs circulating in the blood from crossing into the interstitial fluid of the brain where neurons reside, many researchers are working hard on developing drug delivery systems to penetrate the BBB which currently poses a challenge. Thus, blood-brain barrier penetrating peptides (B3PPs) are an alternative neurotherapeutic for brain-related disorder since they can facilitate drug delivery into the brain. In the meanwhile, developing computational methods that are effective for both the identification and characterization of B3PPs in a cost-effective manner plays an important role for basic reach and in the pharmaceutical industry. Even though few computational methods for B3PP identification have been developed, their performance might fail in terms of generalization ability and interpretability. In this study, a novel and efficient scoring card method-based predictor (termed SCMB3PP) is presented for improving B3PP identification and characterization. To overcome the limitation of black-box computational approaches, the SCMB3PP predictor can automatically estimate amino acid and dipeptide propensities to be B3PPs. Both cross-validation and independent tests indicate that SCMB3PP can achieve impressive performance and outperform various popular machine learning-based methods and the existing methods on multiple independent test datasets. Furthermore, SCMB3PP-derived amino acid propensities were utilized to identify informative biophysical and biochemical properties for characterizing B3PPs. Finally, an online user-friendly web server ( http://pmlabstack.pythonanywhere.com/SCMB3PP ) is established to identify novel and potential B3PP cost-effectively. This novel computational approach is anticipated to facilitate the large-scale identification of high potential B3PP candidates for follow-up experimental validation.
Subject(s)
Blood-Brain Barrier , Dipeptides , Dipeptides/chemistry , Dipeptides/metabolism , Propensity Score , Endothelial Cells , Peptides/metabolism , Amino Acids/chemistryABSTRACT
BACKGROUND: Traumatic brain injury (TBI) is an extremely heterogeneous and complex pathology that requires the integration of different physiological measurements for the optimal understanding and clinical management of patients. Information derived from intracranial pressure (ICP) monitoring can be coupled with information obtained from heart rate (HR) monitoring to assess the interplay between brain and heart. The goal of our study is to investigate events of simultaneous increases in HR and ICP and their relationship with patient mortality.. METHODS: In our previous work, we introduced a novel measure of brain-heart interaction termed brain-heart crosstalks (ctnp), as well as two additional brain-heart crosstalks indicators [mutual information ([Formula: see text]) and average edge overlap (ωct)] obtained through a complex network modeling of the brain-heart system. These measures are based on identification of simultaneous increase of HR and ICP. In this article, we investigated the relationship of these novel indicators with respect to mortality in a multicenter TBI cohort, as part of the Collaborative European Neurotrauma Effectiveness Research in TBI high-resolution work package. RESULTS: A total of 226 patients with TBI were included in this cohort. The data set included monitored parameters (ICP and HR), as well as laboratory, demographics, and clinical information. The number of detected brain-heart crosstalks varied (mean 58, standard deviation 57). The Kruskal-Wallis test comparing brain-heart crosstalks measures of survivors and nonsurvivors showed statistically significant differences between the two distributions (p values: 0.02 for [Formula: see text], 0.005 for ctnp and 0.006 for ωct). An inverse correlation was found, computed using the point biserial correlation technique, between the three new measures and mortality: - 0.13 for ctnp (p value 0.04), - 0.19 for ωct (p value 0.002969) and - 0.09 for [Formula: see text] (p value 0.1396). The measures were then introduced into the logistic regression framework, along with a set of input predictors made of clinical, demographic, computed tomography (CT), and lab variables. The prediction models were obtained by dividing the original cohort into four age groups (16-29, 30-49, 50-65, and 65-85 years of age) to properly treat with the age confounding factor. The best performing models were for age groups 16-29, 50-65, and 65-85, with the deviance of ratio explaining more than 80% in all the three cases. The presence of an inverse relationship between brain-heart crosstalks and mortality was also confirmed. CONCLUSIONS: The presence of a negative relationship between mortality and brain-heart crosstalks indicators suggests that a healthy brain-cardiovascular interaction plays a role in TBI.
Subject(s)
Brain Injuries, Traumatic/physiopathology , Brain/physiopathology , Heart Rate/physiology , Heart/physiology , Intracranial Pressure/physiology , Adolescent , Adult , Aged , Aged, 80 and over , Brain/diagnostic imaging , Brain Injuries, Traumatic/mortality , Cohort Studies , Humans , Middle Aged , Monitoring, Physiologic , Young AdultABSTRACT
BACKGROUND: High-throughput sequencing Chromosome Conformation Capture (Hi-C) allows the study of DNA interactions and 3D chromosome folding at the genome-wide scale. Usually, these data are represented as matrices describing the binary contacts among the different chromosome regions. On the other hand, a graph-based representation can be advantageous to describe the complex topology achieved by the DNA in the nucleus of eukaryotic cells. METHODS: Here we discuss the use of a graph database for storing and analysing data achieved by performing Hi-C experiments. The main issue is the size of the produced data and, working with a graph-based representation, the consequent necessity of adequately managing a large number of edges (contacts) connecting nodes (genes), which represents the sources of information. For this, currently available graph visualisation tools and libraries fall short with Hi-C data. The use of graph databases, instead, supports both the analysis and the visualisation of the spatial pattern present in Hi-C data, in particular for comparing different experiments or for re-mapping omics data in a space-aware context efficiently. In particular, the possibility of describing graphs through statistical indicators and, even more, the capability of correlating them through statistical distributions allows highlighting similarities and differences among different Hi-C experiments, in different cell conditions or different cell types. RESULTS: These concepts have been implemented in NeoHiC, an open-source and user-friendly web application for the progressive visualisation and analysis of Hi-C networks based on the use of the Neo4j graph database (version 3.5). CONCLUSION: With the accumulation of more experiments, the tool will provide invaluable support to compare neighbours of genes across experiments and conditions, helping in highlighting changes in functional domains and identifying new co-organised genomic compartments.
Subject(s)
Chromatin , Chromosomes , Chromatin/genetics , Genome , Genomics , Molecular ConformationABSTRACT
BACKGROUND: Single-cell RNA sequencing (scRNA-Seq) experiments are gaining ground to study the molecular processes that drive normal development as well as the onset of different pathologies. Finding an effective and efficient low-dimensional representation of the data is one of the most important steps in the downstream analysis of scRNA-Seq data, as it could provide a better identification of known or putatively novel cell-types. Another step that still poses a challenge is the integration of different scRNA-Seq datasets. Though standard computational pipelines to gain knowledge from scRNA-Seq data exist, a further improvement could be achieved by means of machine learning approaches. RESULTS: Autoencoders (AEs) have been effectively used to capture the non-linearities among gene interactions of scRNA-Seq data, so that the deployment of AE-based tools might represent the way forward in this context. We introduce here scAEspy, a unifying tool that embodies: (1) four of the most advanced AEs, (2) two novel AEs that we developed on purpose, (3) different loss functions. We show that scAEspy can be coupled with various batch-effect removal tools to integrate data by different scRNA-Seq platforms, in order to better identify the cell-types. We benchmarked scAEspy against the most used batch-effect removal tools, showing that our AE-based strategies outperform the existing solutions. CONCLUSIONS: scAEspy is a user-friendly tool that enables using the most recent and promising AEs to analyse scRNA-Seq data by only setting up two user-defined parameters. Thanks to its modularity, scAEspy can be easily extended to accommodate new AEs to further improve the downstream analysis of scRNA-Seq data. Considering the relevant results we achieved, scAEspy can be considered as a starting point to build a more comprehensive toolkit designed to integrate multi single-cell omics.
Subject(s)
RNA , Single-Cell Analysis , Machine Learning , RNA/genetics , Sequence Analysis, RNA , Exome SequencingABSTRACT
OBJECTIVE: In a previous study, we observed the presence of simultaneous increases in intracranial pressure (ICP) and the heart rate (HR), which we denominated cardio-cerebral crosstalk (CC), and we related the number of such events to patient outcomes in a paediatric cohort. In this chapter, we present an extension of this work to an adult cohort from the Collaborative European NeuroTrauma Effectiveness Research in TBI (CENTER-TBI) study. METHODS: We implemented a sliding window algorithm to detect CC events. We considered subwindows of 10-min observations. If simultaneous increases of at least 20% in ICP and HR occurred with respect to the minimum ICP and HR values in the time windows, a CC event was detected. Correlation between the number of CC events and mortality was then obtained. RESULTS: The cohort consisted of 226 adults (aged 16-85 years). The number of CC events that were detected varied (mean 50, standard deviation 58). A point biserial correlation coefficient of -0.13 between mortality and CC was found. Although the correlation was weaker than that seen in the paediatric cohort (-0.30), the negative direction was replicated. CONCLUSION: In this work, we first extracted CC events from ICP and HR observations of adult patients with traumatic brain injury and related the number of CC events to patient outcomes. Consistency with the previous results in the paediatric cohort was observed. The more crosstalk events occurred, the better the patient outcome was.
Subject(s)
Brain Injuries, Traumatic , Adolescent , Adult , Aged , Aged, 80 and over , Algorithms , Brain Injuries, Traumatic/complications , Child , Cohort Studies , Heart Rate , Humans , Intracranial Pressure , Middle Aged , Young AdultABSTRACT
Accurate identification of bitter peptides is of great importance for better understanding their biochemical and biophysical properties. To date, machine learning-based methods have become effective approaches for providing a good avenue for identifying potential bitter peptides from large-scale protein datasets. Although few machine learning-based predictors have been developed for identifying the bitterness of peptides, their prediction performances could be improved. In this study, we developed a new predictor (named iBitter-Fuse) for achieving more accurate identification of bitter peptides. In the proposed iBitter-Fuse, we have integrated a variety of feature encoding schemes for providing sufficient information from different aspects, namely consisting of compositional information and physicochemical properties. To enhance the predictive performance, the customized genetic algorithm utilizing self-assessment-report (GA-SAR) was employed for identifying informative features followed by inputting optimal ones into a support vector machine (SVM)-based classifier for developing the final model (iBitter-Fuse). Benchmarking experiments based on both 10-fold cross-validation and independent tests indicated that the iBitter-Fuse was able to achieve more accurate performance as compared to state-of-the-art methods. To facilitate the high-throughput identification of bitter peptides, the iBitter-Fuse web server was established and made freely available online. It is anticipated that the iBitter-Fuse will be a useful tool for aiding the discovery and de novo design of bitter peptides.