CiteFuse enables multi-modal analysis of CITE-seq data.

MOTIVATION: Multi-modal profiling of single cells represents one of the latest technological advancements in molecular biology. Among various single-cell multi-modal strategies, cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) allows simultaneous quantification of two distinct species: RNA and cell-surface proteins. Here, we introduce CiteFuse, a streamlined package consisting of a suite of tools for doublet detection, modality integration, clustering, differential RNA and protein expression analysis, antibody-derived tag evaluation, ligand-receptor interaction analysis and interactive web-based visualization of CITE-seq data. RESULTS: We demonstrate the capacity of CiteFuse to integrate the two data modalities and its relative advantage against data generated from single-modality profiling using both simulations and real-world CITE-seq data. Furthermore, we illustrate a novel doublet detection method based on a combined index of cell hashing and transcriptome data. Finally, we demonstrate CiteFuse for predicting ligand-receptor interactions by using multi-modal CITE-seq data. Collectively, we demonstrate the utility and effectiveness of CiteFuse for the integrative analysis of transcriptome and epitope profiles from CITE-seq data. AVAILABILITY AND IMPLEMENTATION: CiteFuse is freely available at http://shiny.maths.usyd.edu.au/CiteFuse/ as an online web service and at https://github.com/SydneyBioX/CiteFuse/ as an R package. CONTACT: pengyi.yang@sydney.edu.au. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Autoencoder-based cluster ensembles for single-cell RNA-seq data analysis.

BACKGROUND: Single-cell RNA-sequencing (scRNA-seq) is a transformative technology, allowing global transcriptomes of individual cells to be profiled with high accuracy. An essential task in scRNA-seq data analysis is the identification of cell types from complex samples or tissues profiled in an experiment. To this end, clustering has become a key computational technique for grouping cells based on their transcriptome profiles, enabling subsequent cell type identification from each cluster of cells. Due to the high feature-dimensionality of the transcriptome (i.e. the large number of measured genes in each cell) and because only a small fraction of genes are cell type-specific and therefore informative for generating cell type-specific clusters, clustering directly on the original feature/gene dimension may lead to uninformative clusters and hinder correct cell type identification. RESULTS: Here, we propose an autoencoder-based cluster ensemble framework in which we first take random subspace projections from the data, then compress each random projection to a low-dimensional space using an autoencoder artificial neural network, and finally apply ensemble clustering across all encoded datasets to generate clusters of cells. We employ four evaluation metrics to benchmark clustering performance and our experiments demonstrate that the proposed autoencoder-based cluster ensemble can lead to substantially improved cell type-specific clusters when applied with both the standard k-means clustering algorithm and a state-of-the-art kernel-based clustering algorithm (SIMLR) designed specifically for scRNA-seq data. Compared to directly using these clustering algorithms on the original datasets, the performance improvement in some cases is up to 100%, depending on the evaluation metric used. CONCLUSIONS: Our results suggest that the proposed framework can facilitate more accurate cell type identification as well as other downstream analyses. The code for creating the proposed autoencoder-based cluster ensemble framework is freely available from https://github.com/gedcom/scCCESS.

scReClassify: post hoc cell type classification of single-cell rNA-seq data.

BACKGROUND: Single-cell RNA-sequencing (scRNA-seq) is a fast emerging technology allowing global transcriptome profiling on the single cell level. Cell type identification from scRNA-seq data is a critical task in a variety of research such as developmental biology, cell reprogramming, and cancers. Typically, cell type identification relies on human inspection using a combination of prior biological knowledge (e.g. marker genes and morphology) and computational techniques (e.g. PCA and clustering). Due to the incompleteness of our current knowledge and the subjectivity involved in this process, a small amount of cells may be subject to mislabelling. RESULTS: Here, we propose a semi-supervised learning framework, named scReClassify, for 'post hoc' cell type identification from scRNA-seq datasets. Starting from an initial cell type annotation with potentially mislabelled cells, scReClassify first performs dimension reduction using PCA and next applies a semi-supervised learning method to learn and subsequently reclassify cells that are likely mislabelled initially to the most probable cell types. By using both simulated and real-world experimental datasets that profiled various tissues and biological systems, we demonstrate that scReClassify is able to accurately identify and reclassify misclassified cells to their correct cell types. CONCLUSIONS: scReClassify can be used for scRNA-seq data as a post hoc cell type classification tool to fine-tune cell type annotations generated by any cell type classification procedure. It is implemented as an R package and is freely available from https://github.com/SydneyBioX/scReClassify.

Multiplexed Temporal Quantification of the Exercise-regulated Plasma Peptidome.

Exercise is extremely beneficial to whole body health reducing the risk of a number of chronic human diseases. Some of these physiological benefits appear to be mediated via the secretion of peptide/protein hormones into the blood stream. The plasma peptidome contains the entire complement of low molecular weight endogenous peptides derived from secretion, protease activity and PTMs, and is a rich source of hormones. In the current study we have quantified the effects of intense exercise on the plasma peptidome to identify novel exercise regulated secretory factors in humans. We developed an optimized 2D-LC-MS/MS method and used multiple fragmentation methods including HCD and EThcD to analyze endogenous peptides. This resulted in quantification of 5,548 unique peptides during a time course of exercise and recovery. The plasma peptidome underwent dynamic and large changes during exercise on a time-scale of minutes with many rapidly reversible following exercise cessation. Among acutely regulated peptides, many were known hormones including insulin, glucagon, ghrelin, bradykinin, cholecystokinin and secretogranins validating the method. Prediction of bioactive peptides regulated with exercise identified C-terminal peptides from Transgelins, which were increased in plasma during exercise. In vitro experiments using synthetic peptides identified a role for transgelin peptides on the regulation of cell-cycle, extracellular matrix remodeling and cell migration. We investigated the effects of exercise on the regulation of PTMs and proteolytic processing by building a site-specific network of protease/substrate activity. Collectively, our deep peptidomic analysis of plasma revealed that exercise rapidly modulates the circulation of hundreds of bioactive peptides through a network of proteases and PTMs. These findings illustrate that peptidomics is an ideal method for quantifying changes in circulating factors on a global scale in response to physiological perturbations such as exercise. This will likely be a key method for pinpointing exercise regulated factors that generate health benefits.

Preoperative comprehensive geriatric assessment and multidisciplinary team input in older people undergoing elective orthopaedic surgery: A feasibility trial.

OBJECTIVE: To determine the feasibility of preoperative comprehensive geriatric assessment (CGA) and multidisciplinary team (MDT) input for older people undergoing elective orthopaedic surgery in a tertiary New Zealand setting. METHODS: This single-centre retrospective study included elective orthopaedic patients older than 65 years (and Maori/Pasifika aged greater than 55 years) with hyperpolypharmacy, frailty, neurocognitive disorders and poor functional status. Patients attended a preoperative clinic where they had a geriatrician-led CGA along with MDT input. The feasibility of this preoperative model was assessed using outcomes of acceptability, accessibility and adherence. A qualitative description of patient demographics along with clinic assessment and interventions further describes this pilot experience. RESULTS: Sixty patients met inclusion criteria. This group were vulnerable older people (median age 77 years), with a high incidence of hyperpolypharmacy (85%), frailty (80%) and neurocognitive disorders (30%). Acceptability was high (97%), along with CGA accessibility (100%); however, MDT accessibility varied (53-90%). Adherence to MDT intervention was low; with only 26% of patients completing physiotherapy sessions and only 29% adhering to dietary advice. Accurate recall was a significant factor contributing to poor adherence. Comprehensive geriatric assessment was demonstrated to be a broad and flexible intervention. CONCLUSIONS: CGA with MDT input is an acceptable and accessible intervention to be utilised as part of improved preoperative care for the older person undergoing elective orthopaedic surgery. Further consideration around methods to increase adherence in this patient group should be explored. Future research should focus on refining the intervention, and quantifying impact on patient outcomes.

SnapKin: a snapshot deep learning ensemble for kinase-substrate prediction from phosphoproteomics data.

A major challenge in mass spectrometry-based phosphoproteomics lies in identifying the substrates of kinases, as currently only a small fraction of substrates identified can be confidently linked with a known kinase. Machine learning techniques are promising approaches for leveraging large-scale phosphoproteomics data to computationally predict substrates of kinases. However, the small number of experimentally validated kinase substrates (true positive) and the high data noise in many phosphoproteomics datasets together limit their applicability and utility. Here, we aim to develop advanced kinase-substrate prediction methods to address these challenges. Using a collection of seven large phosphoproteomics datasets, and both traditional and deep learning models, we first demonstrate that a 'pseudo-positive' learning strategy for alleviating small sample size is effective at improving model predictive performance. We next show that a data resampling-based ensemble learning strategy is useful for improving model stability while further enhancing prediction. Lastly, we introduce an ensemble deep learning model ('SnapKin') by incorporating the above two learning strategies into a 'snapshot' ensemble learning algorithm. We propose SnapKin, an ensemble deep learning method, for predicting substrates of kinases from large-scale phosphoproteomics data. We demonstrate that SnapKin consistently outperforms existing methods in kinase-substrate prediction. SnapKin is freely available at https://github.com/PYangLab/SnapKin.

Genome-wide analysis in Drosophila reveals diet-by-gene interactions and uncovers diet-responsive genes.

Genetic and environmental factors play a major role in metabolic health. However, they do not act in isolation, as a change in an environmental factor such as diet may exert different effects based on an individual's genotype. Here, we sought to understand how such gene-diet interactions influenced nutrient storage and utilization, a major determinant of metabolic disease. We subjected 178 inbred strains from the Drosophila genetic reference panel (DGRP) to diets varying in sugar, fat, and protein. We assessed starvation resistance, a holistic phenotype of nutrient storage and utilization that can be robustly measured. Diet influenced the starvation resistance of most strains, but the effect varied markedly between strains such that some displayed better survival on a high carbohydrate diet (HCD) compared to a high-fat diet while others had opposing responses, illustrating a considerable gene × diet interaction. This demonstrates that genetics plays a major role in diet responses. Furthermore, heritability analysis revealed that the greatest genetic variability arose from diets either high in sugar or high in protein. To uncover the genetic variants that contribute to the heterogeneity in starvation resistance, we mapped 566 diet-responsive SNPs in 293 genes, 174 of which have human orthologs. Using whole-body knockdown, we identified two genes that were required for glucose tolerance, storage, and utilization. Strikingly, flies in which the expression of one of these genes, CG4607 a putative homolog of a mammalian glucose transporter, was reduced at the whole-body level, displayed lethality on a HCD. This study provides evidence that there is a strong interplay between diet and genetics in governing survival in response to starvation, a surrogate measure of nutrient storage efficiency and obesity. It is likely that a similar principle applies to higher organisms thus supporting the case for nutrigenomics as an important health strategy.

Signaling Heterogeneity is Defined by Pathway Architecture and Intercellular Variability in Protein Expression.

Insulin's activation of PI3K/Akt signaling, stimulates glucose uptake by enhancing delivery of GLUT4 to the cell surface. Here we examined the origins of intercellular heterogeneity in insulin signaling. Akt activation alone accounted for ~25% of the variance in GLUT4, indicating that additional sources of variance exist. The Akt and GLUT4 responses were highly reproducible within the same cell, suggesting the variance is between cells (extrinsic) and not within cells (intrinsic). Generalized mechanistic models (supported by experimental observations) demonstrated that the correlation between the steady-state levels of two measured signaling processes decreases with increasing distance from each other and that intercellular variation in protein expression (as an example of extrinsic variance) is sufficient to account for the variance in and between Akt and GLUT4. Thus, the response of a population to insulin signaling is underpinned by considerable single-cell heterogeneity that is largely driven by variance in gene/protein expression between cells.

Akt phosphorylates insulin receptor substrate to limit PI3K-mediated PIP3 synthesis.

The phosphoinositide 3-kinase (PI3K)-Akt network is tightly controlled by feedback mechanisms that regulate signal flow and ensure signal fidelity. A rapid overshoot in insulin-stimulated recruitment of Akt to the plasma membrane has previously been reported, which is indicative of negative feedback operating on acute timescales. Here, we show that Akt itself engages this negative feedback by phosphorylating insulin receptor substrate (IRS) 1 and 2 on a number of residues. Phosphorylation results in the depletion of plasma membrane-localised IRS1/2, reducing the pool available for interaction with the insulin receptor. Together these events limit plasma membrane-associated PI3K and phosphatidylinositol (3,4,5)-trisphosphate (PIP3) synthesis. We identified two Akt-dependent phosphorylation sites in IRS2 at S306 (S303 in mouse) and S577 (S573 in mouse) that are key drivers of this negative feedback. These findings establish a novel mechanism by which the kinase Akt acutely controls PIP3 abundance, through post-translational modification of the IRS scaffold.

For the body to work properly, cells must constantly 'talk' to each other using signalling molecules. Receiving a chemical signal triggers a series of molecular events in a cell, a so-called 'signal transduction pathway' that connects a signal with a precise outcome. Disturbing cell signalling can trigger disease, and strict control mechanisms are therefore in place to ensure that communication does not break down or become erratic. For instance, just as a thermostat turns off the heater once the right temperature is reached, negative feedback mechanisms in cells switch off signal transduction pathways when the desired outcome has been achieved. The hormone insulin is a signal for growth that increases in the body following a meal to promote the storage of excess blood glucose (sugar) in muscle and fat cells. The hormone binds to insulin receptors at the cell surface and switches on a signal transduction pathway that makes the cell take up glucose from the bloodstream. If the signal is not engaged diseases such as diabetes develop. Conversely, if the signal cannot be adequately switched of cancer can develop. Determining exactly how insulin works would help to understand these diseases better and to develop new treatments. Kearney et al. therefore set out to examine the biochemical 'fail-safes' that control insulin signalling. Experiments using computer simulations of the insulin signalling pathway revealed a potential new mechanism for negative feedback, which centred on a molecule known as Akt. The models predicted that if the negative feedback were removed, then Akt would become hyperactive and accumulate at the cell's surface after stimulation with insulin. Further manipulation of the 'virtual' insulin signalling pathway and studies of live cells in culture confirmed that this was indeed the case. The cell biology experiments also showed how Akt, once at the cell surface, was able to engage the negative feedback and shut down further insulin signalling. Akt did this by inactivating a protein required to pass the signal from the insulin receptor to the rest of the cell. Overall, this work helps to understand cell communication by revealing a previously unknown, and critical component of the insulin signalling pathway.