Decoding the hallmarks of allograft dysfunction with a comprehensive pan-organ transcriptomic atlas.

The pathogenesis of allograft (dys)function has been increasingly studied using 'omics'-based technologies, but the focus on individual organs has created knowledge gaps that neither unify nor distinguish pathological mechanisms across allografts. Here we present a comprehensive study of human pan-organ allograft dysfunction, analyzing 150 datasets with more than 12,000 samples across four commonly transplanted solid organs (heart, lung, liver and kidney, n = 1,160, 1,241, 1,216 and 8,853 samples, respectively) that we leveraged to explore transcriptomic differences among allograft dysfunction (delayed graft function, acute rejection and fibrosis), tolerance and stable graft function. We identified genes that correlated robustly with allograft dysfunction across heart, lung, liver and kidney transplantation. Furthermore, we developed a transfer learning omics prediction framework that, by borrowing information across organs, demonstrated superior classifications compared to models trained on single organs. These findings were validated using a single-center prospective kidney transplant cohort study (a collective 329 samples across two timepoints), providing insights supporting the potential clinical utility of our approach. Our study establishes the capacity for machine learning models to learn across organs and presents a transcriptomic transplant resource that can be employed to develop pan-organ biomarkers of allograft dysfunction.

Data integration and inference of gene regulation using single-cell temporal multimodal data with scTIE.

Single-cell technologies offer unprecedented opportunities to dissect gene regulatory mechanisms in context-specific ways. Although there are computational methods for extracting gene regulatory relationships from scRNA-seq and scATAC-seq data, the data integration problem, essential for accurate cell type identification, has been mostly treated as a standalone challenge. Here we present scTIE, a unified method that integrates temporal multimodal data and infers regulatory relationships predictive of cellular state changes. scTIE uses an autoencoder to embed cells from all time points into a common space by using iterative optimal transport, followed by extracting interpretable information to predict cell trajectories. Using a variety of synthetic and real temporal multimodal data sets, we show scTIE achieves effective data integration while preserving more biological signals than existing methods, particularly in the presence of batch effects and noise. Furthermore, on the exemplar multiome data set we generated from differentiating mouse embryonic stem cells over time, we show scTIE captures regulatory elements highly predictive of cell transition probabilities, providing new potentials to understand the regulatory landscape driving developmental processes.

BIDCell: Biologically-informed self-supervised learning for segmentation of subcellular spatial transcriptomics data.

Recent advances in subcellular imaging transcriptomics platforms have enabled high-resolution spatial mapping of gene expression, while also introducing significant analytical challenges in accurately identifying cells and assigning transcripts. Existing methods grapple with cell segmentation, frequently leading to fragmented cells or oversized cells that capture contaminated expression. To this end, we present BIDCell, a self-supervised deep learning-based framework with biologically-informed loss functions that learn relationships between spatially resolved gene expression and cell morphology. BIDCell incorporates cell-type data, including single-cell transcriptomics data from public repositories, with cell morphology information. Using a comprehensive evaluation framework consisting of metrics in five complementary categories for cell segmentation performance, we demonstrate that BIDCell outperforms other state-of-the-art methods according to many metrics across a variety of tissue types and technology platforms. Our findings underscore the potential of BIDCell to significantly enhance single-cell spatial expression analyses, enabling great potential in biological discovery.

Pcdh19 mediates olfactory sensory neuron coalescence during postnatal stages and regeneration.

The mouse olfactory system regenerates constantly throughout life. While genes critical for the initial projection of olfactory sensory neurons (OSNs) to the olfactory bulb have been identified, what genes are important for maintaining the olfactory map during regeneration are still unknown. Here we show a mutation in Protocadherin 19 (Pcdh19), a cell adhesion molecule and member of the cadherin superfamily, leads to defects in OSN coalescence during regeneration. Surprisingly, lateral glomeruli were more affected and males in particular showed a more severe phenotype. Single cell analysis unexpectedly showed OSNs expressing the MOR28 odorant receptor could be subdivided into two major clusters. We showed that at least one protocadherin is differentially expressed between OSNs coalescing on the medial and lateral glomeruli. Moreover, females expressed a slightly different complement of genes from males. These features may explain the differential effects of mutating Pcdh19 on medial and lateral glomeruli in males and females.

Atlas-scale single-cell multi-sample multi-condition data integration using scMerge2.

The recent emergence of multi-sample multi-condition single-cell multi-cohort studies allows researchers to investigate different cell states. The effective integration of multiple large-cohort studies promises biological insights into cells under different conditions that individual studies cannot provide. Here, we present scMerge2, a scalable algorithm that allows data integration of atlas-scale multi-sample multi-condition single-cell studies. We have generalized scMerge2 to enable the merging of millions of cells from single-cell studies generated by various single-cell technologies. Using a large COVID-19 data collection with over five million cells from 1000+ individuals, we demonstrate that scMerge2 enables multi-sample multi-condition scRNA-seq data integration from multiple cohorts and reveals signatures derived from cell-type expression that are more accurate in discriminating disease progression. Further, we demonstrate that scMerge2 can remove dataset variability in CyTOF, imaging mass cytometry and CITE-seq experiments, demonstrating its applicability to a broad spectrum of single-cell profiling technologies.

Lipidomics Profiling and Risk of Coronary Artery Disease in the BioHEART-CT Discovery Cohort.

The current coronary artery disease (CAD) risk scores for predicting future cardiovascular events rely on well-recognized traditional cardiovascular risk factors derived from a population level but often fail individuals, with up to 25% of first-time heart attack patients having no risk factors. Non-invasive imaging technology can directly measure coronary artery plaque burden. With an advanced lipidomic measurement methodology, for the first time, we aim to identify lipidomic biomarkers to enable intervention before cardiovascular events. With 994 participants from BioHEART-CT Discovery Cohort, we collected clinical data and performed high-performance liquid chromatography with mass spectrometry to determine concentrations of 683 plasma lipid species. Statin-naive participants were selected based on subclinical CAD (sCAD) categories as the analytical cohort (n = 580), with sCAD+ (n = 243) compared to sCAD- (n = 337). Through a machine learning approach, we built a lipid risk score (LRS) and compared the performance of the existing Framingham Risk Score (FRS) in predicting sCAD+. We obtained individual classifiability scores and determined Body Mass Index (BMI) as the modifying variable. FRS and LRS models achieved similar areas under the receiver operating characteristic curve (AUC) in predicting the validation cohort. LRS enhanced the prediction of sCAD+ in the healthy-weight group (BMI < 25 kg/m2), where FRS performed poorly and identified individuals at risk that FRS missed. Lipid features have strong potential as biomarkers to predict CAD plaque burden and can identify residual risk not captured by traditional risk factors/scores. LRS compliments FRS in prediction and has the most significant benefit in healthy-weight individuals.

scTIE: data integration and inference of gene regulation using single-cell temporal multimodal data.

Single-cell technologies offer unprecedented opportunities to dissect gene regulatory mechanisms in context-specific ways. Although there are computational methods for extracting gene regulatory relationships from scRNA-seq and scATAC-seq data, the data integration problem, essential for accurate cell type identification, has been mostly treated as a standalone challenge. Here we present scTIE, a unified method that integrates temporal multimodal data and infers regulatory relationships predictive of cellular state changes. scTIE uses an autoencoder to embed cells from all time points into a common space using iterative optimal transport, followed by extracting interpretable information to predict cell trajectories. Using a variety of synthetic and real temporal multimodal datasets, we demonstrate scTIE achieves effective data integration while preserving more biological signals than existing methods, particularly in the presence of batch effects and noise. Furthermore, on the exemplar multiome dataset we generated from differentiating mouse embryonic stem cells over time, we demonstrate scTIE captures regulatory elements highly predictive of cell transition probabilities, providing new potentials to understand the regulatory landscape driving developmental processes.

Overcoming cohort heterogeneity for the prediction of subclinical cardiovascular disease risk.

Cardiovascular disease remains a leading cause of mortality with an estimated half a billion people affected in 2019. However, detecting signals between specific pathophysiology and coronary plaque phenotypes using complex multi-omic discovery datasets remains challenging due to the diversity of individuals and their risk factors. Given the complex cohort heterogeneity present in those with coronary artery disease (CAD), we illustrate several different methods, both knowledge-guided and data-driven approaches, for identifying subcohorts of individuals with subclinical CAD and distinct metabolomic signatures. We then demonstrate that utilizing these subcohorts can improve the prediction of subclinical CAD and can facilitate the discovery of novel biomarkers of subclinical disease. Analyses acknowledging cohort heterogeneity through identifying and utilizing these subcohorts may be able to advance our understanding of CVD and provide more effective preventative treatments to reduce the burden of this disease in individuals and in society as a whole.

NEMoE: a nutrition aware regularized mixture of experts model to identify heterogeneous diet-microbiome-host health interactions.

BACKGROUND: Unrevealing the interplay between diet, the microbiome, and the health state could enable the design of personalized intervention strategies and improve the health and well-being of individuals. A common approach to this is to divide the study population into smaller cohorts based on dietary preferences in the hope of identifying specific microbial signatures. However, classification of patients based solely on diet is unlikely to reflect the microbiome-host health relationship or the taxonomic microbiome makeup. RESULTS: We present a novel approach, the Nutrition-Ecotype Mixture of Experts (NEMoE) model, for establishing associations between gut microbiota and health state that accounts for diet-specific cohort variability using a regularized mixture of experts model framework with an integrated parameter sharing strategy to ensure data-driven diet-cohort identification consistency across taxonomic levels. The success of our approach was demonstrated through a series of simulation studies, in which NEMoE showed robustness with regard to parameter selection and varying degrees of data heterogeneity. Further application to real-world microbiome data from a Parkinson's disease cohort revealed that NEMoE is capable of not only improving predictive performance for Parkinson's Disease but also for identifying diet-specific microbial signatures of disease. CONCLUSION: In summary, NEMoE can be used to uncover diet-specific relationships between nutritional-ecotype and patient health and to contextualize precision nutrition for different diseases. Video Abstract.

Deep multimodal graph-based network for survival prediction from highly multiplexed images and patient variables.

The spatial architecture of the tumour microenvironment and phenotypic heterogeneity of tumour cells have been shown to be associated with cancer prognosis and clinical outcomes, including survival. Recent advances in highly multiplexed imaging, including imaging mass cytometry (IMC), capture spatially resolved, high-dimensional maps that quantify dozens of disease-relevant biomarkers at single-cell resolution, that contain potential to inform patient-specific prognosis. Existing automated methods for predicting survival, on the other hand, typically do not leverage spatial phenotype information captured at the single-cell level. Furthermore, there is no end-to-end method designed to leverage the rich information in whole IMC images and all marker channels, and aggregate this information with clinical data in a complementary manner to predict survival with enhanced accuracy. To that end, we present a deep multimodal graph-based network (DMGN) with two modules: (1) a multimodal graph-based module that considers relationships between spatial phenotype information in all image regions and all clinical variables adaptively, and (2) a clinical embedding module that automatically generates embeddings specialised for each clinical variable to enhance multimodal aggregation. We demonstrate that our modules are consistently effective at improving survival prediction performance using two public breast cancer datasets, and that our new approach can outperform state-of-the-art methods in survival prediction.

scSTAR reveals hidden heterogeneity with a real-virtual cell pair structure across conditions in single-cell RNA sequencing data.

Cell-state transition can reveal additional information from single-cell ribonucleic acid (RNA)-sequencing data in time-resolved biological phenomena. However, most of the current methods are based on the time derivative of the gene expression state, which restricts them to the short-term evolution of cell states. Here, we present single-cell State Transition Across-samples of RNA-seq data (scSTAR), which overcomes this limitation by constructing a paired-cell projection between biological conditions with an arbitrary time span by maximizing the covariance between two feature spaces using partial least square and minimum squared error methods. In mouse ageing data, the response to stress in CD4+ memory T cell subtypes was found to be associated with ageing. A novel Treg subtype characterized by mTORC activation was identified to be associated with antitumour immune suppression, which was confirmed by immunofluorescence microscopy and survival analysis in 11 cancers from The Cancer Genome Atlas Program. On melanoma data, scSTAR improved immunotherapy-response prediction accuracy from 0.8 to 0.96.

Viral Integration Plays a Minor Role in the Development and Prognostication of Oral Squamous Cell Carcinoma.

Viruses are well known drivers of several human malignancies. A causative factor for oral cavity squamous cell carcinoma (OSCC) in patients with limited exposure to traditional risk factors, including tobacco use, is yet to be identified. Our study aimed to comprehensively evaluate the role of viral drivers in OSCC patients with low cumulative exposure to traditional risk factors. Patients under 50 years of age with OSCC, defined using strict anatomic criteria were selected for WGS. The WGS data was interrogated using viral detection tools (Kraken 2 and BLASTN), together examining >700,000 viruses. The findings were further verified using tissue microarrays of OSCC samples using both immunohistochemistry and RNA in situ hybridisation (ISH). 28 patients underwent WGS and comprehensive viral profiling. One 49-year-old male patient with OSCC of the hard palate demonstrated HPV35 integration. 657 cases of OSCC were then evaluated for the presence of HPV integration through immunohistochemistry for p16 and HPV RNA ISH. HPV integration was seen in 8 (1.2%) patients, all middle-aged men with predominant floor of mouth involvement. In summary, a wide-ranging interrogation of >700,000 viruses using OSCC WGS data showed HPV integration in a minority of male OSCC patients and did not carry any prognostic significance.

Scalable workflow for characterization of cell-cell communication in COVID-19 patients.

COVID-19 patients display a wide range of disease severity, ranging from asymptomatic to critical symptoms with high mortality risk. Our ability to understand the interaction of SARS-CoV-2 infected cells within the lung, and of protective or dysfunctional immune responses to the virus, is critical to effectively treat these patients. Currently, our understanding of cell-cell interactions across different disease states, and how such interactions may drive pathogenic outcomes, is incomplete. Here, we developed a generalizable and scalable workflow for identifying cells that are differentially interacting across COVID-19 patients with distinct disease outcomes and use this to examine eight public single-cell RNA-seq datasets (six from peripheral blood mononuclear cells, one from bronchoalveolar lavage and one from nasopharyngeal), with a total of 211 individual samples. By characterizing the cell-cell interaction patterns across epithelial and immune cells in lung tissues for patients with varying disease severity, we illustrate diverse communication patterns across individuals, and discover heterogeneous communication patterns among moderate and severe patients. We further illustrate patterns derived from cell-cell interactions are potential signatures for discriminating between moderate and severe patients. Overall, this workflow can be generalized and scaled to combine multiple scRNA-seq datasets to uncover cell-cell interactions.

SurvBenchmark: comprehensive benchmarking study of survival analysis methods using both omics data and clinical data.

Survival analysis is a branch of statistics that deals with both the tracking of time and the survival status simultaneously as the dependent response. Current comparisons of survival model performance mostly center on clinical data with classic statistical survival models, with prediction accuracy often serving as the sole metric of model performance. Moreover, survival analysis approaches for censored omics data have not been thoroughly investigated. The common approach is to binarize the survival time and perform a classification analysis. Here, we develop a benchmarking design, SurvBenchmark, that evaluates a diverse collection of survival models for both clinical and omics data sets. SurvBenchmark not only focuses on classical approaches such as the Cox model but also evaluates state-of-the-art machine learning survival models. All approaches were assessed using multiple performance metrics; these include model predictability, stability, flexibility, and computational issues. Our systematic comparison design with 320 comparisons (20 methods over 16 data sets) shows that the performances of survival models vary in practice over real-world data sets and over the choice of the evaluation metric. In particular, we highlight that using multiple performance metrics is critical in providing a balanced assessment of various models. The results in our study will provide practical guidelines for translational scientists and clinicians, as well as define possible areas of investigation in both survival technique and benchmarking strategies.

Cross-Platform Omics Prediction procedure: a statistical machine learning framework for wider implementation of precision medicine.

In this modern era of precision medicine, molecular signatures identified from advanced omics technologies hold great promise to better guide clinical decisions. However, current approaches are often location-specific due to the inherent differences between platforms and across multiple centres, thus limiting the transferability of molecular signatures. We present Cross-Platform Omics Prediction (CPOP), a penalised regression model that can use omics data to predict patient outcomes in a platform-independent manner and across time and experiments. CPOP improves on the traditional prediction framework of using gene-based features by selecting ratio-based features with similar estimated effect sizes. These components gave CPOP the ability to have a stable performance across datasets of similar biology, minimising the effect of technical noise often generated by omics platforms. We present a comprehensive evaluation using melanoma transcriptomics data to demonstrate its potential to be used as a critical part of a clinical screening framework for precision medicine. Additional assessment of generalisation was demonstrated with ovarian cancer and inflammatory bowel disease studies.

The Gut Microbiome in Parkinson's Disease: A Longitudinal Study of the Impacts on Disease Progression and the Use of Device-Assisted Therapies.

Background: Altered gut microbiome (GM) composition has been established in Parkinson's disease (PD). However, few studies have longitudinally investigated the GM in PD, or the impact of device-assisted therapies. Objectives: To investigate the temporal stability of GM profiles from PD patients on standard therapies and those initiating device-assisted therapies (DAT) and define multivariate models of disease and progression. Methods: We evaluated validated clinical questionnaires and stool samples from 74 PD patients and 74 household controls (HCs) at 0, 6, and 12 months. Faster or slower disease progression was defined from levodopa equivalence dose and motor severity measures. 19 PD patients initiating Deep Brain Stimulation or Levodopa-Carbidopa Intestinal Gel were separately evaluated at 0, 6, and 12 months post-therapy initiation. Results: Persistent underrepresentation of short-chain fatty-acid-producing bacteria, Butyricicoccus, Fusicatenibacter, Lachnospiraceae ND3007 group, and Erysipelotrichaceae UCG-003, were apparent in PD patients relative to controls. A sustained effect of DAT initiation on GM associations with PD was not observed. PD progression analysis indicated that the genus Barnesiella was underrepresented in faster progressing PD patients at t = 0 and t = 12 months. Two-stage predictive modeling, integrating microbiota abundances and nutritional profiles, improved predictive capacity (change in Area Under the Curve from 0.58 to 0.64) when assessed at Amplicon Sequence Variant taxonomic resolution. Conclusion: We present longitudinal GM studies in PD patients, showing persistently altered GM profiles suggestive of a reduced butyrogenic production potential. DATs exerted variable GM influences across the short and longer-term. We found that specific GM profiles combined with dietary factors improved prediction of disease progression in PD patients.

Whole genome duplication in oral squamous cell carcinoma in patients younger than 50 years: implications for prognosis and adverse clinicopathological factors.

INTRODUCTION: Oral squamous cell carcinoma (OSCC) in the young (<50 years), without known carcinogenic risk factors, is on the rise globally. Whole genome duplication (WGD) has been shown to occur at higher rates in cancers without an identifiable carcinogenic agent. We aimed to evaluate the prevalence of WGD in a cohort of OSCC patients under the age of 50 years. METHODS: Whole genome sequencing (WGS) was performed on 28 OSCC patients from the Sydney Head and Neck Cancer Institute (SHNCI) biobank. An additional nine cases were obtained from The Cancer Genome Atlas (TCGA). RESULTS: WGD was seen in 27 of 37 (73%) cases. Non-synonymous, somatic TP53 mutations occurred in 25 of 27 (93%) cases of WGD and were predicted to precede WGD in 21 (77%). WGD was significantly associated with larger tumor size (p = 0.01) and was frequent in patients with recurrences (87%, p = 0.36). Overall survival was significantly worse in those with WGD (p = 0.05). CONCLUSIONS: Our data, based on one of the largest WGS datasets of young patients with OSCC, demonstrates a high frequency of WGD and its association with adverse pathologic characteristics and clinical outcomes. TP53 mutations also preceded WGD, as has been described in other tumors without a clear mutagenic driver.

Nutritional Intake and Gut Microbiome Composition Predict Parkinson's Disease.

Background: Models to predict Parkinson's disease (PD) incorporating alterations of gut microbiome (GM) composition have been reported with varying success. Objective: To assess the utility of GM compositional changes combined with macronutrient intake to develop a predictive model of PD. Methods: We performed a cross-sectional analysis of the GM and nutritional intake in 103 PD patients and 81 household controls (HCs). GM composition was determined by 16S amplicon sequencing of the V3-V4 region of bacterial ribosomal DNA isolated from stool. To determine multivariate disease-discriminant associations, we developed two models using Random Forest and support-vector machine (SVM) methodologies. Results: Using updated taxonomic reference, we identified significant compositional differences in the GM profiles of PD patients in association with a variety of clinical PD characteristics. Six genera were overrepresented and eight underrepresented in PD patients relative to HCs, with the largest difference being overrepresentation of Lactobacillaceae at family taxonomic level. Correlation analyses highlighted multiple associations between clinical characteristics and select taxa, whilst constipation severity, physical activity and pharmacological therapies associated with changes in beta diversity. The random forest model of PD, incorporating taxonomic data at the genus level and carbohydrate contribution to total energy demonstrated the best predictive capacity [Area under the ROC Curve (AUC) of 0.74]. Conclusion: The notable differences in GM diversity and composition when combined with clinical measures and nutritional data enabled the development of a predictive model to identify PD. These findings support the combination of GM and nutritional data as a potentially useful biomarker of PD to improve diagnosis and guide clinical management.

Computational approaches for direct cell reprogramming: from the bulk omics era to the single cell era.

Recent advances in direct cell reprogramming have made possible the conversion of one cell type to another cell type, offering a potential cell-based treatment to many major diseases. Despite much attention, substantial roadblocks remain including the inefficiency in the proportion of reprogrammed cells of current experiments, and the requirement of a significant amount of time and resources. To this end, several computational algorithms have been developed with the goal of guiding the hypotheses to be experimentally validated. These approaches can be broadly categorized into two main types: transcription factor identification methods which aim to identify candidate transcription factors for a desired cell conversion, and transcription factor perturbation methods which aim to simulate the effect of a transcription factor perturbation on a cell state. The transcription factor perturbation methods can be broken down into Boolean networks, dynamical systems and regression models. We summarize the contributions and limitations of each method and discuss the innovation that single cell technologies are bringing to these approaches and we provide a perspective on the future direction of this field.