Thrombospondin-1 Drives Cardiac Remodeling in Chronic Kidney Disease.

Patients with chronic kidney disease (CKD) face a high risk of cardiovascular disease. Previous studies reported that endogenous thrombospondin 1 (TSP1) involves right ventricular remodeling and dysfunction. Here we show that a murine model of CKD increased myocardial TSP1 expression and produced left ventricular hypertrophy, fibrosis, and dysfunction. TSP1 knockout mice were protected from these features. In vitro, indoxyl sulfate is driving deleterious changes in cardiomyocyte through the TSP1. In patients with CKD, TSP1 and aryl hydrocarbon receptor were both differentially expressed in the myocardium. Our findings summon large clinical studies to confirm the translational role of TSP1 in patients with CKD.

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.

Construction and optimization of multi-platform precision pathways for precision medicine.

In the enduring challenge against disease, advancements in medical technology have empowered clinicians with novel diagnostic platforms. Whilst in some cases, a single test may provide a confident diagnosis, often additional tests are required. However, to strike a balance between diagnostic accuracy and cost-effectiveness, one must rigorously construct the clinical pathways. Here, we developed a framework to build multi-platform precision pathways in an automated, unbiased way, recommending the key steps a clinician would take to reach a diagnosis. We achieve this by developing a confidence score, used to simulate a clinical scenario, where at each stage, either a confident diagnosis is made, or another test is performed. Our framework provides a range of tools to interpret, visualize and compare the pathways, improving communication and enabling their evaluation on accuracy and cost, specific to different contexts. This framework will guide the development of novel diagnostic pathways for different diseases, accelerating the implementation of precision medicine into clinical practice.

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.

Mass cytometry analysis reveals altered immune profiles in patients with coronary artery disease.

Objective: The importance of inflammation in atherosclerosis is well accepted, but the role of the adaptive immune system is not yet fully understood. To further explore this, we assessed the circulating immune cell profile of patients with coronary artery disease (CAD) to identify discriminatory features by mass cytometry. Methods: Mass cytometry was performed on patient samples from the BioHEART-CT study, gated to detect 82 distinct cell subsets. CT coronary angiograms were analysed to categorise patients as having CAD (CAD+) or having normal coronary arteries (CAD-). Results: The discovery cohort included 117 patients (mean age 61 ± 12 years, 49% female); 79 patients (68%) were CAD+. Mass cytometry identified changes in 15 T-cell subsets, with higher numbers of proliferating, highly differentiated and cytotoxic cells and decreases in naïve T cells. Five T-regulatory subsets were related to an age and gender-independent increase in the odds of CAD incidence when expressing CCR2 (OR 1.12), CCR4 (OR 1.08), CD38 and CD45RO (OR 1.13), HLA-DR (OR 1.06) and Ki67 (OR 1.22). Markers of proliferation and differentiation were also increased within B cells, while plasmacytoid dendritic cells were decreased. This combination of changes was assessed using SVM models in discovery and validation cohorts (area under the curve = 0.74 for both), confirming the robust nature of the immune signature detected. Conclusion: We identified differences within immune subpopulations of CAD+ patients which are indicative of a systemic immune response to coronary atherosclerosis. This immune signature needs further study via incorporation into risk scoring tools for the precision diagnosis of CAD.

The impact of similarity metrics on cell-type clustering in highly multiplexed in situ imaging cytometry data.

Motivation: The advent of highly multiplexed in situ imaging cytometry assays has revolutionized the study of cellular systems, offering unparalleled detail in observing cellular activities and characteristics. These assays provide comprehensive insights by concurrently profiling the spatial distribution and molecular features of numerous cells. In navigating this complex data landscape, unsupervised machine learning techniques, particularly clustering algorithms, have become essential tools. They enable the identification and categorization of cell types and subsets based on their molecular characteristics. Despite their widespread adoption, most clustering algorithms in use were initially developed for cell suspension technologies, leading to a potential mismatch in application. There is a critical gap in the systematic evaluation of these methods, particularly in determining the properties that make them optimal for in situ imaging assays. Addressing this gap is vital for ensuring accurate, reliable analyses and fostering advancements in cellular biology research. Results: In our extensive investigation, we evaluated a range of similarity metrics, which are crucial in determining the relationships between cells during the clustering process. Our findings reveal substantial variations in clustering performance, contingent on the similarity metric employed. These variations underscore the importance of selecting appropriate metrics to ensure accurate cell type and subset identification. In response to these challenges, we introduce FuseSOM, a novel ensemble clustering algorithm that integrates hierarchical multiview learning of similarity metrics with self-organizing maps. Through a rigorous stratified subsampling analysis framework, we demonstrate that FuseSOM outperforms existing best-practice clustering methods specifically tailored for in situ imaging cytometry data. Our work not only provides critical insights into the performance of clustering algorithms in this novel context but also offers a robust solution, paving the way for more accurate and reliable in situ imaging cytometry data analysis. Availability and implementation: The FuseSOM R package is available on Bioconductor and is available under the GPL-3 license. All the codes for the analysis performed can be found at Github.

Circulating cytokine and chemokine patterns associated with cytomegalovirus reactivation after stem cell transplantation.

Objectives: Human cytomegalovirus (HCMV) reactivation is the leading viral complication after allogeneic haematopoietic stem cell transplantation (allo-HSCT). Understanding of circulating cytokine/chemokine patterns which accompany HCMV reactivation and correlate with HCMV DNAemia magnitude is limited. We aimed to characterise plasma cytokine/chemokine profiles in 36 allo-HSCT patients (21 with HCMV reactivation and 15 without HCMV reactivation) at four time-points in the first 100-day post-transplant. Methods: The concentrations of 31 cytokines/chemokines in plasma samples were analysed using a multiplex bead-based immunoassay. Cytokine/chemokine concentrations were compared in patients with high-level HCMV DNAemia, low-level HCMV DNAemia or no HCMV reactivation, and correlated with immune cell frequencies measured using mass cytometry. Results: Increased plasma levels of T helper 1-type cytokines/chemokines (TNF, IL-18, IP-10, MIG) were detected in patients with HCMV reactivation at the peak of HCMV DNAemia, relative to non-reactivators. Stem cell factor (SCF) levels were significantly higher before the detection of HCMV reactivation in patients who went on to develop high-level HCMV DNAemia (810-52 740 copies/mL) vs. low-level HCMV DNAemia (< 250 copies/mL). High-level HCMV reactivators, but not low-level reactivators, developed an elevated inflammatory cytokine/chemokine profile (MIP-1α, MIP-1ß, TNF, LT-α, IL-13, IL-9, SCF, HGF) at the peak of reactivation. Plasma cytokine concentrations displayed unique correlations with circulating immune cell frequencies in patients with HCMV reactivation. Conclusion: This study identifies distinct circulating cytokine/chemokine signatures associated with the magnitude of HCMV DNAemia and the progression of HCMV reactivation after allo-HSCT, providing important insight into immune recovery patterns associated with HCMV reactivation and viral control.

Classification of the tumor immune microenvironment and associations with outcomes in patients with metastatic melanoma treated with immunotherapies.

BACKGROUND: Tumor microenvironment (TME) characteristics are potential biomarkers of response to immune checkpoint inhibitors in metastatic melanoma. This study developed a method to perform unsupervised classification of TME of metastatic melanoma. METHODS: We used multiplex immunohistochemical and quantitative pathology-derived assessment of immune cell compositions of intratumoral and peritumoral regions of metastatic melanoma baseline biopsies to classify TME in relation to response to anti-programmed cell death protein 1 (PD-1) monotherapy or in combination with anti-cytotoxic T-cell lymphocyte-4 (ipilimumab (IPI)+PD-1). RESULTS: Spatial profiling of CD8+T cells, macrophages, and melanoma cells, as well as phenotypic PD-1 receptor ligand (PD-L1) and CD16 proportions, were used to identify and classify patients into one of three mutually exclusive TME classes: immune-scarce, immune-intermediate, and immune-rich tumors. Patients with immune-rich tumors were characterized by a lower proportion of melanoma cells and higher proportions of immune cells, including higher PD-L1 expression. These patients had higher response rates and longer progression-free survival (PFS) than those with immune-intermediate and immune-scarce tumors. At a median follow-up of 18 months (95% CI: 6.7 to 49 months), the 1-year PFS was 76% (95% CI: 64% to 90%) for patients with an immune-rich tumor, 56% (95% CI: 44% to 72%) for those with an immune-intermediate tumor, and 33% (95% CI: 23% to 47%) for patients with an immune-scarce tumor. A higher response rate was observed in patients with an immune-scarce or immune-intermediate tumor when treated with IPI+PD-1 compared with those treated with PD-1 alone. CONCLUSIONS: Our study provides an automatic TME classification method that may predict the clinical efficacy of immunotherapy for patients with metastatic melanoma.

Non-invasive diagnosis of early cutaneous squamous cell carcinoma.

Early cutaneous squamous cell carcinoma (cSCC) can be challenging to diagnose using clinical criteria as it could present similar to actinic keratosis (AK) or Bowen's disease (BD), precursors of cSCC. Currently, histopathological assessment of an invasive biopsy is the gold standard for diagnosis. A non-invasive diagnostic approach would reduce patient and health system burden. Therefore, this study used non-invasive sampling by tape-stripping coupled with data-independent acquisition mass spectrometry (DIA-MS) proteomics to profile the proteome of histopathologically diagnosed AK, BD and cSCC, as well as matched normal samples. Proteomic data were analysed to identify proteins and biological functions that are significantly different between lesions. Additionally, a support vector machine (SVM) machine learning algorithm was used to assess the usefulness of proteomic data for the early diagnosis of cSCC. A total of 696 proteins were identified across the samples studied. A machine learning model constructed using the proteomic data classified premalignant (AK + BD) and malignant (cSCC) lesions at 77.5% accuracy. Differential abundance analysis identified 144 and 21 protein groups that were significantly changed in the cSCC, and BD samples compared to the normal skin, respectively (adj. p < 0.05). Changes in pivotal carcinogenic pathways such as LXR/RXR activation, production of reactive oxygen species, and Hippo signalling were observed that may explain the progression of cSCC from premalignant lesions. In summary, this study demonstrates that DIA-MS analysis of tape-stripped samples can identify non-invasive protein biomarkers with the potential to be developed into a complementary diagnostic tool for early cSCC.

MoleculeExperiment enables consistent infrastructure for molecule-resolved spatial omics data in bioconductor.

MOTIVATION: Imaging-based spatial transcriptomics (ST) technologies have achieved subcellular resolution, enabling detection of individual molecules in their native tissue context. Data associated with these technologies promise unprecedented opportunity toward understanding cellular and subcellular biology. However, in R/Bioconductor, there is a scarcity of existing computational infrastructure to represent such data, and particularly to summarize and transform it for existing widely adopted computational tools in single-cell transcriptomics analysis, including SingleCellExperiment and SpatialExperiment (SPE) classes. With the emergence of several commercial offerings of imaging-based ST, there is a pressing need to develop consistent data structure standards for these technologies at the individual molecule-level. RESULTS: To this end, we have developed MoleculeExperiment, an R/Bioconductor package, which (i) stores molecule and cell segmentation boundary information at the molecule-level, (ii) standardizes this molecule-level information across different imaging-based ST technologies, including 10× Genomics' Xenium, and (iii) streamlines transition from a MoleculeExperiment object to a SpatialExperiment object. Overall, MoleculeExperiment is generally applicable as a data infrastructure class for consistent analysis of molecule-resolved spatial omics data. AVAILABILITY AND IMPLEMENTATION: The MoleculeExperiment package is publicly available on Bioconductor at https://bioconductor.org/packages/release/bioc/html/MoleculeExperiment.html. Source code is available on Github at: https://github.com/SydneyBioX/MoleculeExperiment. The vignette for MoleculeExperiment can be found at https://bioconductor.org/packages/release/bioc/html/MoleculeExperiment.html.

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.

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.

Tolerogenic dendritic cells protect against acute kidney injury.

Ischemia reperfusion injury is a common precipitant of acute kidney injury that occurs following disrupted perfusion to the kidney. This includes blood loss and hemodynamic shock, as well as during retrieval for deceased donor kidney transplantation. Acute kidney injury is associated with adverse long-term clinical outcomes and requires effective interventions that can modify the disease process. Immunomodulatory cell therapies such as tolerogenic dendritic cells remain a promising tool, and here we tested the hypothesis that adoptively transferred tolerogenic dendritic cells can limit kidney injury. The phenotypic and genomic signatures of bone marrow-derived syngeneic or allogeneic, Vitamin-D3/IL-10-conditioned tolerogenic dendritic cells were assessed. These cells were characterized by high PD-L1:CD86, elevated IL-10, restricted IL-12p70 secretion and a suppressed transcriptomic inflammatory profile. When infused systemically, these cells successfully abrogated kidney injury without modifying infiltrating inflammatory cell populations. They also provided protection against ischemia reperfusion injury in mice pre-treated with liposomal clodronate, suggesting the process was regulated by live, rather than reprocessed cells. Co-culture experiments and spatial transcriptomic analysis confirmed reduced kidney tubular epithelial cell injury. Thus, our data provide strong evidence that peri-operatively administered tolerogenic dendritic cells have the ability to protect against acute kidney injury and warrants further exploration as a therapeutic option. This technology may provide a clinical advantage for bench-to-bedside translation to affect patient outcomes.

Proteomic profiling of cutaneous melanoma explains the aggressiveness of distant organ metastasis.

Despite recent developments in managing metastatic melanomas, patients' overall survival remains low. Therefore, the current study aims to understand better the proteome-wide changes associated with melanoma metastasis that will assist with identifying targeted therapies. The latest development in mass spectrometry-based proteomics, together with extensive bioinformatics analysis, was used to investigate the molecular changes in 60 formalin-fixed and paraffin-embedded samples of primary and lymph nodes (LN) and distant organ metastatic melanomas. A total of 4631 proteins were identified, of which 72 and 453 were significantly changed between the LN and distant organ metastatic melanomas compared to the primary lesions (adj. p-value <0.05). An increase in proteins such as SLC9A3R1, CD20 and GRB2 and a decrease in CST6, SERPINB5 and ARG1 were associated with regional LN metastasis. By contrast, increased metastatic activities in distant organ metastatic melanomas were related to higher levels of CEACAM1, MC1R, AKT1 and MMP3-9 and decreased levels of CDKN2A, SDC1 and SDC4 proteins. Furthermore, machine learning analysis classified the lesions with up to 92% accuracy based on their metastatic status. The findings from this study provide up to date proteome-level information about the progression of melanomas to regional LN and distant organs, leading to the identification of protein signatures with potential for clinical translation.

Cross-platform comparison of immune signatures in immunotherapy-treated patients with advanced melanoma using a rank-based scoring approach.

BACKGROUND: Gene expression profiling is increasingly being utilised as a diagnostic, prognostic and predictive tool for managing cancer patients. Single-sample scoring approach has been developed to alleviate instability of signature scores due to variations from sample composition. However, it is a challenge to achieve comparable signature scores across different expressional platforms. METHODS: The pre-treatment biopsies from a total of 158 patients, who have received single-agent anti-PD-1 (n = 84) or anti-PD-1 + anti-CTLA-4 therapy (n = 74), were performed using NanoString PanCancer IO360 Panel. Multiple immune-related signature scores were measured from a single-sample rank-based scoring approach, singscore. We assessed the reproducibility and the performance in reporting immune profile of singscore based on NanoString assay in advance melanoma. To conduct cross-platform analyses, singscores between the immune profiles of NanoString assay and the previous orthogonal whole transcriptome sequencing (WTS) data were compared through linear regression and cross-platform prediction. RESULTS: singscore-derived signature scores reported significantly high scores in responders in multiple PD-1, MHC-1-, CD8 T-cell-, antigen presentation-, cytokine- and chemokine-related signatures. We found that singscore provided stable and reproducible signature scores among the repeats in different batches and cross-sample normalisations. The cross-platform comparisons confirmed that singscores derived via NanoString and WTS were comparable. When singscore of WTS generated by the overlapping genes to the NanoString gene set, the signatures generated highly correlated cross-platform scores (Spearman correlation interquartile range (IQR) [0.88, 0.92] and r2 IQR [0.77, 0.81]) and better prediction on cross-platform response (AUC = 86.3%). The model suggested that Tumour Inflammation Signature (TIS) and Personalised Immunotherapy Platform (PIP) PD-1 are informative signatures for predicting immunotherapy-response outcomes in advanced melanoma patients treated with anti-PD-1-based therapies. CONCLUSIONS: Overall, the outcome of this study confirms that singscore based on NanoString data is a feasible approach to produce reliable signature scores for determining patients' immune profiles and the potential clinical utility in biomarker implementation, as well as to conduct cross-platform comparisons, such as WTS.

Spatial analysis for highly multiplexed imaging data to identify tissue microenvironments.

Highly multiplexed in situ imaging cytometry assays have made it possible to study the spatial organization of numerous cell types simultaneously. We have addressed the challenge of quantifying complex multi-cellular relationships by proposing a statistical method which clusters local indicators of spatial association. Our approach successfully identifies distinct tissue architectures in datasets generated from three state-of-the-art high-parameter assays demonstrating its value in summarizing the information-rich data generated from these technologies.

Spatial mapping reveals granuloma diversity and histopathological superstructure in human tuberculosis.

The hallmark of tuberculosis (TB) is the formation of immune cell-enriched aggregates called granulomas. While granulomas are pathologically diverse, their tissue-wide heterogeneity has not been spatially resolved at the single-cell level in human tissues. By spatially mapping individual immune cells in every lesion across entire tissue sections, we report that in addition to necrotizing granulomas, the human TB lung contains abundant non-necrotizing leukocyte aggregates surrounding areas of necrotizing tissue. These cellular lesions were more diverse in composition than necrotizing lesions and could be stratified into four general classes based on cellular composition and spatial distribution of B cells and macrophages. The cellular composition of non-necrotizing structures also correlates with their proximity to necrotizing lesions, indicating these are foci of distinct immune reactions adjacent to necrotizing granulomas. Together, we show that during TB, diseased lung tissue develops a histopathological superstructure comprising at least four different types of non-necrotizing cellular aggregates organized as satellites of necrotizing granulomas.