Manifold Approximation Research Articles

A comprehensive human embryo reference tool using single-cell RNA-sequencing data.

Stem cell-based embryo models offer unprecedented experimental tools for studying early human development. The usefulness of embryo models hinges on their molecular, cellular and structural fidelities to their in vivo counterparts. To authenticate human embryo models, single-cell RNA sequencing has been utilized for unbiased transcriptional profiling. However, an organized and integrated human single-cell RNA-sequencing dataset, serving as a universal reference for benchmarking human embryo models, remains unavailable. Here we developed such a reference through the integration of six published human datasets covering development from the zygote to the gastrula. Lineage annotations are contrasted and validated with available human and nonhuman primate datasets. Using stabilized Uniform Manifold Approximation and Projection, we constructed an early embryogenesis prediction tool, where query datasets can be projected on the reference and annotated with predicted cell identities. Using this reference tool, we examined published human embryo models, highlighting the risk of misannotation when relevant references are not utilized for benchmarking and authentication.

A machine learning approach to stratify patients with hypermobile Ehlers-Danlos syndrome/hypermobility spectrum disorders according to disorders of gut brain interaction, comorbidities and quality of life.

A high prevalence of disorders of gut-brain interaction (DGBI) exist in patients with hypermobile Ehlers-Danlos Syndrome (hEDS) and hypermobility spectrum disorders (HSD). However, it is unknown if clusters of hEDS/HSD patients exist which overlap with different DGBIs and whether this overlap influences presence of comorbidities and quality of life. We aimed to study these knowledge gaps. A prospectively collected hEDS/HSD cohort of 1044 individuals were studied. We undertook Uniform Manifold Approximation and Projection-enabled (UMAP) dimension reduction to create a representation of nonlinear interactions between hEDS/HSD and DGBIs, from which individuals were stratified into clusters. Somatization, Postural Tachycardia Syndrome (PoTS), autonomic symptoms, psychological factors and quality of life were statistically compared between clusters. The mean age of patients was 40 ± 13.2 years; 87.8% were female. Patients segregated into three clusters: Cluster 0 (n = 466): hEDS/HSD+ functional foregut disorders (FFD) + irritable bowel syndrome (IBS); Cluster 1 (n = 180): hEDS/HSD+ IBS and Cluster 2 (n = 337): hEDS/HSD alone. In cluster 0, we demonstrated increased somatization (p <0.0001), anxiety (p <0.0001), depression (p <0.0001), PoTS prevalence (p = 0.003), autonomic symptoms (p <0.0001) and reduced quality of life (p <0.0001) compared to cluster 2. Cluster 0 had greater comorbidity burden than cluster 1. Within hEDS/HSD, subgroups exist with a high prevalence of FFD and IBS. These subgroups have a higher prevalence of psychological disorders, dysautonomia and poorer quality of life compared with hEDS/HSD alone. Further research should focus on healthcare utilization, management and prognosis in hEDS/HSD and DGBI overlap.

Abstract 4135951: Single Nuclei RNA Sequencing Reveals Distinct Cellular Composition of Pericoronary Adipose Tissue

Introduction: Adipose tissue is a metabolically active organ with local paracrine interactions on adjacent tissues. Pericoronary adipose tissue (PcAT) is an emerging cardiovascular disease (CVD) risk factor and likely has paracrine effects on coronary vasculature that may promote coronary artery disease (CAD) development. PcAT may be functionally and clinically distinct from epicardial adipose tissue (EAT) given its proximity to the coronary artery adventitia. The cellular composition of PcAT and EAT and whether they have unique inflammatory cell signatures and functional expression at a single cell level has not been explored. Hypothesis: PcAT will have diverse cellular composition with heightened inflammatory cell content compared to EAT. Methods: Visual dissection of paired EAT not adjacent to an epicardial coronary artery and PcAT in direct apposition to an epicardial artery (left anterior descending artery) was performed from a single human heart immediately after explant (prior to heart transplantation). Nuclei were isolated using the Chromium Nuclei Isolation Kit (10x Genomics) with RNAse inhibitor. Cells were counted with a Cellometer K2 and Chromium Next GEM Single Cell libraries were prepared, followed by reverse transcription, PCR and NovaSeq 6000 platform sequencing. Raw sequence reads from single nuclei RNA-seq (snRNA-seq) were processed and after QC filtering, cells were clustered using Seurat. Clusters were annotated using singleR. Results: In PcAT and EAT, &gt;7000 cells were successfully extracted with 73755 and 48750 mean snRNA-seq reads per cell, respectively. snRNA-seq cellular composition analysis delineated a heightened inflammatory milieu and distinct cellular composition of human PcAT vs. EAT. A UMAP (Uniform Manifold Approximation and Projection) plot displays higher inflammatory cell infiltrate with a neutrophil predominance and lesser adipocyte composition in PcAT compared with EAT ( Figure 1) . Conclusions: We have generated snRNA-seq data from human PcAT and shown distinct cellular composition compared to EAT. Given the proximity of PcAT to coronary vasculature, these results are supportive of paracrine inflammatory effects and highlight adipose tissue biology in cardiovascular disease. Understanding this biology may identify novel therapeutic targets for CVD.

Abstract 4114276: Identification of the heterogeneity in the risk of heart failure in those with left bundle branch block pattern using artificial intelligence

Introduction: The association between left bundle branch block (LBBB) and heart failure is well-established. However, the lack of a practical risk assessment tool for patients with LBBB poses a clinical challenge to the management despite the heterogeneity of the LBBB population. Hypothesis: Artificial intelligence (AI) can extract features to stratify the risk of heart failure (HF) within the heterogeneity of LBBB. Methods: All 12-lead ECGs diagnosed with LBBB recorded from 1/1/2015 to 12/31/2019 were identified in our institution. ECG data were analyzed as 2D matrices with a shape of 12×2500. The Uniform Manifold Approximation and Projection (UMAP) was used to visualize the heterogeneity of LBBB ECG. Additionally, a 2-dimensional convolutional neural network model was trained to detect LBBB ECGs associated with past HF diagnosis. The model was then applied to ECGs from an external cohort of patients with LBBB but without a history of HF. Cumulative incidences of HF admission were compared by stratifying according to tertiles of the model output (low, intermediate and high AI score groups) Results: We identified 15,124 LBBB ECGs in our institution (training dataset) and 2,463 individuals with LBBB ECGs in the external cohort (external test dataset). The UMAP projection of the ECGs revealed distinct heterogeneity in the data ( Fig A ). This heterogeneity was not fully explained by patient demographics, ECG features, or institutions. When the external cohort was divided into 3 groups according to the tertile of AI prediction, patients with high AI scores had a higher risk of HF admission (high vs low AI score group: hazard ratio (HR), 2.10; 95% confidence interval (CI), 1.66-2.65; intermediate vs low AI score group: HR, 1.39; 95%CI, 1.08-1.78: log-rank p &lt;0.0001; Fig 1B ). This finding was unchanged after adjusting for multiple clinical characteristics and ECG parameters, including QRS duration. Conclusion: We observed significant heterogeneity in ECGs presenting LBBB, and the AI model excellently stratified the risk of heart failure admission in patients with LBBB from a single recording of 12-lead ECG, suggesting great potential for clinical utility in managing patients with LBBB.

Abstract 4129223: Large Language Modeling-Enabled Analysis of Atrial Fibrillation Topics and Sentiments on Social Media

Purpose: We utilized large language modeling (LLM) and automated topic modeling techniques to describe public perceptions regarding AF on Reddit, a social media platform. Methods: In this qualitative study, we curated all AF-related discussions on Reddit from 2006 to 2023 by searching for all posts and comments containing AF-related terms (“atrial fibrillation”, “afib”, and “atrial flutter”). Discussions were embedded using a state-of-the-art embedding model (BAAI/bge-base-en-v1.5). Embeddings are reduced in dimensionality using Uniform Manifold Approximation and Project and clustered using K-Means Clustering to find topics and groups. Topic and group labels are enhanced using large language models (Llama2 and GPT4, respectively). Results: We identified a total of 86,323 discussions related to AF from 38,183 unique users between April 14, 2006 and November 20, 2023. Discussions around wearables and ablations increased over time, while discussions around anticoagulation had periodic spikes. The topic modeling pipeline automatically identified 64 topics and 9 groups (Figure 1). Groups emphasized the following themes: lifestyle factors and triggers associated with AF (groups 1, 5, 7); patient experiences with management strategies including anticoagulation, rate control, pharmacological rhythm control, and ablation (groups 2, 4, 5); experiences with wearable devices (group 8); psychological burden of living with AF (group 6); and online advertisements for AF-related medications (groups 3, 9). Conclusions: Our AI-enabled topic modeling analysis revealed public perceptions about AF including triggers, experiences with different management strategies, impact of wearables, and psychological burden of the diagnosis, which may help guide potential strategies for improving shared decision making around the AF patient journey.

Abstract 4138373: Machine learning identifies clinically distinct phenotypes in patients with aortic regurgitation.

Background: Aortic regurgitation (AR) is a prevalent valve disease with a long latent period to symptoms. Recent data has suggested the role of novel markers of myocardial overload in assessing onset of decompensation. Method: We sought to evaluate the role of unsupervised cluster analyses in identifying different clinical clusters, including clinical status, and a large number of echocardiographic variables including left ventricular (LV) volumes, and their association with mortality. Patients with ≥moderate-severe chronic AR identified using echocardiography at Mayo Clinic, Rochester were retrospectively analyzed. Primary outcome was all-cause mortality censored at aortic valve surgery/last follow-up. Uniform Manifold Approximation and Projection (UMAP) with K -means algorithm was used to cluster patients using clinical and, echocardiographic variables at the time of presentation. Missing data were imputed with the Multiple Imputation by Chained Equations (MICE) method. A supervised approach trained on the training set was used to find cluster membership in a hold-out validation set. Log-rank tests were used to assess differences in mortality rates between the clusters, both in the training and validation sets. Results: Three distinct clusters were identified among 1100 patients (log-rank P for survival &lt;0.001). Cluster 1 (n=337), which included younger males with severe AR but fewer symptoms, showed the best survival, 75.6% (69.5, 82.3). Cluster 2 (n=235), older and more females with elevated filling pressures, showed intermediate survival of 64.2 % (56.8, 72.5). Cluster 3 (n=253), characterized by severe symptomatic AR, demonstrated the lowest survival of 45.3 % (34.4, 59.8) at 5 years. Similar clusters were formed in the internal validation cohort. Conclusion: Distinct clusters with variable echocardiographic features and mortality differences exist within patients with chronic ≥moderate-severe AR. Recognizing these clusters can refine individual risk stratification and clinical decision-making after verification in future prospective studies.

IMMU-03. EXPLORING CANCER-SPECIFIC T CELLS AND ANTIGENS IN GLIOBLASTOMA BASED ON SINGLE-CELL SEQUENCING OF CD8+ TILS

Abstract BACKGROUND The effectiveness of cancer immunotherapy against glioblastoma (GBM) remains limited. This study aims to identify cancer-specific antigens to develop antigen-based cancer immunotherapies for GBM. We investigated candidate tumor antigen-specific T cells in GBM by single-cell RNA sequencing (scRNA-seq) and single-cell TCR sequencing (scTCR-seq), and we explored candidate antigens through genetic analysis of tumor tissues. METHODS Flow cytometry analysis was conducted on fresh tumor digest samples to evaluate tumor-infiltrating lymphocytes (TILs) of GBM. Single-cell RNA and TCR sequencing were performed on CD8+ T cells in TILs. Both data generated by Cell Ranger software were loaded into Seurat at R studio. The dimensional reduction for clustering was performed with uniform manifold approximation and projection (UMAP). Additionally, we performed bulk RNA sequencing (RNA-seq) and whole exome sequencing (WES) on tumor tissues. RESULTS Two out of 20 patients (10%) exhibited abundant TILs in GBM. From these two patients, approximately 20,000 CD8+ T cells in total were analyzed in single-cell analysis. Clustering based on UMAP revealed a predominance of clusters expressing exhaustion markers, with high TCR clonality centered on exhausted T cell (TEX) clusters, like our previous date of lung cancer TILs recognizing tumor antigens. Moreover, within the TEX clusters, the expression of CXCL13, often reported as an antigen-specific marker in recent years, was significantly higher in our data compared to two publicly available datasets. Subsequently, based on RNA-seq and WES, 61 candidate neoantigens were estimated using machine-learning prediction algorithms from NEC Corporation. Furthermore, 5 overexpressed cancer/testis antigens, and 3 bacterial species-derived microbial peptides using Kraken2 were selected as antigen candidates. CONCLUSIONS We identified prospective tumor antigen-specific T cell subsets with high CXCL13 expression from two GBM patients. We plan to identify cancer-specific T cells and antigens from these candidates by T cell activation assay.

PATH-22. SINGLE-CELL RNA (SCRNA) SEQUENCING OF PAIRED PRIMARY AND RECURRENT ATYPICAL TERATOID/RHABDOID TUMORS (AT/RT) REVEALED MOLECULAR DIFFERENCES AND THERAPY-RESISTANT PRIMARY TUMOR CELLS

Abstract Atypical teratoid/rhabdoid tumors (AT/RT) are the most common malignant brain tumors during infancy. They split into four molecular types. The major three (AT/RT-SHH, AT/RT-TYR, and AT/RT-MYC) all carry mutations in SMARCB1, the fourth smaller type is characterized by SMARCA4 mutations (AT/RT-SMARCA4). AT/RT are associated with a dismal outcome, and 50% of the patients suffer a progress or recurrence of the tumor. Molecular mechanisms of disease recurrence are currently unclear. Here, we identified molecular characteristics of AT/RT primaries and recurrences, which could help to improve the therapy for relapsed patients. Previously, we investigated tumor tissue from 26 patients revealing signatures of recurrences in comparison to matched primary tumor samples using histology, DNA methylation analysis, and bulk RNA sequencing. To further identify mechanisms regarding tumor progression and recurrence, we performed single-cell RNA (scRNA) sequencing for five of these primary-recurrence pairs so far using formalin-fixed paraffin-embedded (FFPE) tissue. Using Uniform Manifold Approximation and Projection (UMAP) for dimensionality reduction, we observed a distinct separation in the gene expression profiles of primary and recurrent tumor cells. In addition, we aimed to identify primary tumor cells which potentially persisted after therapy and developed to recurrences. Therefore, we identified primary tumor cells showing highest similarity to recurrent tumor cells using CIBERSORT and trajectory analysis. Differentially expressed genes of these specific primary tumor cells were investigated. They are involved in cell migration, -adhesion, -proliferation, and -stability, among others. Overall, molecular characteristics that occur in the course of the disease were identified, which may help to define new therapeutic targets for AT/RT recurrences.

PATH-21. SINGEL CELL TRANSCRIPTOMICS OF SPINAL EPENDYMOMA SUBTYPES IDENTIFIES INTRATUMORAL HETEROGENEITY

Abstract Spinal ependymomas comprise glial neoplasms that are classified by histological and molecular features. Two of the respective types, myxopapillary ependymomas (MPE) and spinal ependymomas (SP-EPN) are largely understudied. For each of these, our group identified clinically relevant subtypes with poorer (MPE-A; SP-EPN-A) and more favorable progression-free survival (MPE-B; SP-EPN-B) employing global methylation and transcriptomic profiling. As cellular origin and pathogenesis of these tumors are largely unknown, therapeutic options for patients with unfavorable subtypes are unavailable. We performed single cell transcriptomic sequencing of 12 formalin-fixed and paraffin-embedded spinal ependymomas including all MPE and SP-EPN subtypes and obtained a median of 10,091 cells per tumor after quality control. Uniform Manifold Approximation and Projection (UMAP) analysis revealed intertumoral heterogeneity as most samples formed separate clusters with tumors of the same subtype located in proximity. The expression of literature-defined neoplastic consensus cell states was examined for each cell of the data set and revealed intratumoral heterogeneity. Tumors expressed predominantly two programs, being a ciliated (ependymal-like) cell state (61% of all neoplastic cells) and the astrocytic state (31 % of all neoplastic cells). Interestingly, third and fourth most abundant states (mesenchymal and hypoxia) were almost exclusively expressed in MPE-A and SP-EPN-A. Furthermore, distinct clones of chromosomal copy number variations were identified across all tumors. These clones can be distinguished in UMAP representation of MPE single cell transcriptomes and to a lesser extend in SP-EPN. Next, we sought to characterize the resemblance of integrated tumor subtypes to normal spinal cord cell populations and observed a high similarity to human spinal cord ependymal cells. All spinal ependymoma subtypes display intratumoral heterogeneity with respect to the transcriptomic cell state resemblance of the neoplastic cells. Together with the transcriptomic similarity to ependymal and astrocytic cells, our data suggest the tumor origin within the astro-ependymal lineage.

Unsupervised learning analysis on the proteomes of Zika virus

Background The Zika virus (ZIKV), which is transmitted by mosquito vectors to nonhuman primates and humans, causes devastating outbreaks in the poorest tropical regions of the world. Molecular epidemiology, supported by clustering phylogenetic gold standard studies using sequence data, has provided valuable information for tracking and controlling the spread of ZIKV. Unsupervised learning (UL), a form of machine learning algorithm, can be applied on the datasets without the need of known information for training. Methods In this work, unsupervised Random Forest (URF), followed by the application of dimensional reduction algorithms such as principal component analysis (PCA), Uniform Manifold Approximation and Projection (UMAP), t-distributed stochastic neighbor embedding (t-SNE), and autoencoders were used to uncover hidden patterns from polymorphic amino acid sites extracted on the proteome ZIKV multi-alignments, without the need of an underlying evolutionary model. Results The four UL algorithms revealed specific host and geographical clustering patterns for ZIKV. Among the four dimensionality reduction (DR) algorithms, the performance was better for UMAP. The four algorithms allowed the identification of imported viruses for specific geographical clusters. The UL dimension coordinates showed a significant correlation with phylogenetic tree branch lengths and significant phylogenetic dependence in Abouheif’s Cmean and Pagel’s Lambda tests (p value &lt; 0.01) that showed comparable performance with the phylogenetic method. This analytical strategy was generalizable to an external large dengue type 2 dataset. Conclusion These UL algorithms could be practical evolutionary analytical techniques to track the dispersal of viral pathogens.

BIOM-04. SINGLE-CELL RNA SEQUENCING PROTOCOL OF THE HUMAN METASTATIC PROSTATE SPINE TUMOR MICROENVIRONMENT

Abstract INTRODUCTION 90% of spine tumor patients can develop metastases with debilitating complications, such as sensorimotor dysfunction, paralysis, and spinal instability. As spine tumor patients have limited eligibility to participate in clinical trials and have heterogenous tumor pathology, there is a distinct lack of therapeutic targets and prognostication markers of disease. Thus, there is a distinct need for molecular profiling of spine tumors. OBJECTIVE We describe an OR-to-bench-top processing of spine tumors for single-cell omic studies, with samples from a metastatic prostate cancer patient to demonstrate the ability to compare molecular markers and microenvironments of difficult-to-process spine tumors and adjacent normal tissue. METHODS Under IRB protocol #Pro00101198, patients undergoing spine surgery with the Department of Neurosurgery at Duke University were consented for tissue collection. Vertebral osseous metastases and adjacent normal tissue were confirmed histologically (H&amp;E, immunohistochemistry) and radiologically (MRI, PET, CT, and x-ray). Tissues were manually and enzymatically homogenized in serum-free media with collagenase A and DNAse I to obtain a single cell suspension. Cells were then filter sterilized, erythrolyzed, counted to determine live-death ratio, resuspended to 10-20 million cells per mL, then cryopreserved for downstream standard library preparation for single-cell sequencing using the 10X genomics platform. Downstream analysis was performed using Seurat pipeline, InferCNV analysis, and SingleR. RESULTS P404S and PSA expression confirmed tumor was prostate in origin. Radiological imaging and H&amp;E staining confirmed the presence of T7 osseous metastasis. Uniform manifold approximation and project (UMAP) plots identified 12 unique clusters of immune cell subpopulations between normal and tumorous vertebrae, with their own subsets of differentially expressed genes. CONCLUSION We have developed a protocol to generate granular single cell molecular profiling for spine tumors and metastases. Applications of this protocol can aid in identifying therapeutic targets, molecular predictors of disease, and potentiate increased eligibility of patients for clinical trials.

BIOM-67. DIFFERENTIAL DNA METHYLATION PATTERNS IN THE CSF OF PATIENTS WITH DIFFUSE GLIOMAS

Abstract Studies have shown that different central nervous system (CNS) tumor types exhibit unique DNA methylation profiles. DNA methylation-based machine learning classification of CNS tumor tissue has been demonstrated to be highly-accurate in differentiating between tumor types, including between diffuse glioma subtypes. While similar analyses have been performed using cell-free DNA (cfDNA) from blood, few studies have looked at the potential of DNA methylation as a biomarker for diffuse gliomas using CSF-derived cfDNA. In this study, we analyzed CSF cfDNA methylation patterns in diffuse glioma patients. DNA methylation profiling of CSF of patients with glioblastoma, IDH-wildtype (GBM) (n=37), IDH-mutant diffuse gliomas (n=12), and non-neoplastic conditions (n=5) was performed using cell-free methylated DNA immunoprecipitation and high-throughput sequencing (cfMeDIP-seq). Methylation quantification was estimated for each genomic region. Significance testing was performed between diagnostic groups to discover differentially-methylated regions (DMRs). The 300 most differentially-methylated regions (DMRs) for each comparison (IDH-wildtype vs. IDH-mutant, IDH-mutant vs. control, IDH-wildtype vs. control), for a total of 900 DMRs, were utilized as features in dimensionality reduction analysis and machine learning. Out of ~1×107 defined genomic regions, there were ~7×105 DMRs between GBM and control samples and ~1.5×106 DMRs between GBM and IDH-mutant samples, accounting for ~7% and ~15% of the genome, respectively. Uniform manifold approximation and projection (UMAP) and hierarchical clustering analysis revealed that samples were mostly clustered by diagnosis. Random forest analysis of the 900 DMRs with an 80/20 train/test split resulted in an overall accuracy of 0.815. In conclusion, CSF cfDNA methylation profiles of diffuse glioma patients exhibit differential methylation patterns between IDH-mutant and IDH-wildtype tumors, as well as between tumor and control, and could be informative biomarkers for the diagnosis of patients with diffuse gliomas.

CNSC-52. OLIGODENDROGLIOMA SPATIAL TRANSCRIPTOMIC SEQUENCING REVEALS THERAPY-DEPENDENT EVOLUTION OVER TIME

Abstract Oligodendroglioma evolution is incompletely understood. Here we analyze serial oligodendroglioma samples to define genomic and cellular mechanisms underlying treatment resistance. RNA sequencing was performed on 53,406 spatial transcriptomes from 31 oligodendrogliomas (IDH mutant, 1p/19q-codeleted) that were resected from 16 patients. There were 15 CNS WHO grade 2 and 16 grade 3 oligodendrogliomas, which were stratified for temporal progression across grades (n=5 grade 2 to 2, n=5 grade 2 to 3, n=6 grade 3 to 3) and treatments between serial resections (n=2 TMZ/RT, n=2 TMZ, n=1 mTOR inhibition, n=1 IDH inhibition, n=10 surgical monotherapy). Spatial transcriptomic clusters were defined using cell signature gene sets, cell cycle analysis, differentially expressed marker genes, immunohistochemistry, spatial copy number alterations (CNAs), and deconvolution of 18,628 oligodendroglioma single-cell transcriptomes from 16 independent patients. Specificity was tested using immunofluorescence on independent oligodendrogliomas that did not recur after surgical monotherapy (n=8 samples, median follow-up 8.4 years) and spatial transcriptomic sequencing on 8 matched-paired astrocytomas from 4 patients. Integrated uniform manifold approximation and projection revealed a proneural/NPC-like cluster that was enriched in grade 3 oligodendroglioma samples and represented the starting point in oligodendroglioma lineage differentiation using RNA velocity, Monocle, and subclonal spatiotemporal CNA evolution. Comparison across treatment regimens showed (1) enrichment of the proneural/NPC-like cluster in recurrent versus primary tumors that were treated with surgical monotherapy, (2) enrichment of a proneural/OPC-like cluster and depletion of an oligodendrocyte-like cluster after TMZ/RT, (3) enrichment of proneural/NPC-like and immediate-early response clusters after TMZ/RT, (4) enrichment of the oligodendrocyte-like cluster following mTOR inhibition, and (5) enrichment of cancer associated fibroblasts (CAFs) following IDH inhibition. These data show the spatiotemporal landscape of oligodendroglioma is associated with branched lineage differentiation that is modified by postoperative therapy, and suggest that therapy-dependent druggable dependencies may exist to treat oligodendrogliomas that are resistant to standard interventions.

Machine Learning-enhanced Signature of Metastasis-related T Cell Marker Genes for Predicting Overall Survival in Malignant Melanoma.

In this study, we aimed to investigate disparities in the tumor immune microenvironment (TME) between primary and metastatic malignant melanoma (MM) using single-cell RNA sequencing (scRNA-seq) and to identify metastasis-related T cell marker genes (MRTMGs) for predicting patient survival using machine learning techniques. We identified 6 distinct T cell clusters in 10×scRNA-seq data utilizing the Uniform Manifold Approximation and Projection (UMAP) algorithm. Four machine learning algorithms highlighted SRGN, PMEL, GPR143, EIF4A2, and DSP as pivotal MRTMGs, forming the foundation of the MRTMGs signature. A high MRTMGs signature was found to be correlated with poorer overall survival (OS) and suppression of antitumor immunity in MM patients. We developed a nomogram that combines the MRTMGs signature with the T stage and N stage, which accurately predicts 1-year, 3-year, and 5-year OS probabilities. Furthermore, in an immunotherapy cohort, a high MRTMG signature was associated with an unfavorable response to anti-programmed death 1 (PD-1) therapy. In conclusion, primary and metastatic MM display distinct TME landscapes with different T cell subsets playing crucial roles in metastasis. The MRTMGs signature, established through machine learning, holds potential as a valuable biomarker for predicting the survival of MM patients and their response to anti-PD-1 therapy.

Explainable machine learning models for predicting the acute toxicity of pesticides to sheepshead minnow (Cyprinodon variegatus).

A quantitative structure-activity relationship (QSAR) study was conducted on 313 pesticides to predict their acute toxicity to Sheepshead minnow (Cyprinodon variegatus) by using DRAGON descriptors. Essentials accounting for a reliable model were all considered carefully, giving full consideration to the OECD (Organization for Economic Co-operation and Development) principles for QSAR acceptability in regulation during the model construction and assessment process. Nine variables were selected through the forward stepwise regression method and used as inputs to construct both linear and nonlinear models. The obtained models were validated internally and externally. Generally, machine learning-based methods, namely support vector machine (SVM), random forest (RF), and projection pursuit regression (PPR), perform better than the multiple linear regression (MLR) model. The statistical results (R2 = 0.682-0.933, Q2LOO = 0.604-0.659, Q2F1 = 0.740-0.796, CCC = 0.861-0.882) of the developed models show that they are robust, reliable, reproducible, accurate and predictive. Comparatively, the RF model performs best, giving predictive correlation coefficient Q2 of 0.814, root mean squared error (RMSE) of 0.658 and mean absolute error (MAE) of 0.534 for the test set, respectively. The RF model (as well as SVM and PPR models) was visualized and explained by using the SHapley Additive explanation (SHAP) analysis to enhance its transparency and credibility. In addition, the applicability domain (AD) range of the RF model was characterized by the Williams plot and the tree manifold approximation and projection (TMAP) technology was utilized to illustrate similarity and diversity of the entire data space, to assist in the analysis of the outliers. Activity cliff detection was investigated by using Arithmetic Residuals in K-groups Analysis (ARKA) descriptors. It was found that none of the pesticides was identified as an activity cliff in the training set or a potential prediction cliff in the test set. Therefore, the RF model fulfills each OECD principle in regulation for QSAR models. The research in this work will aid in the in silico QSAR prediction of the acute toxicity to Sheepshead minnow (Cyprinodon variegatus) for untested and new toxic pesticides and can also be extended to other studies.

Single cell RNA-seq reveals cellular and transcriptional heterogeneity in the splenic CD11b+Ly6Chigh monocyte population expanded in sepsis-surviving mice

BackgroundSepsis survivors exhibit immune dysregulation that contributes to poor long-term outcomes. Phenotypic and functional alterations within the myeloid compartment are believed to be a contributing factor. Here we dissect the cellular and transcriptional heterogeneity of splenic CD11b+Ly6Chigh myeloid cells that are expanded in mice that survive the cecal ligation and puncture (CLP) murine model of polymicrobial sepsis to better understand the basis of immune dysregulation in sepsis survivors.MethodsSham or CLP surgeries were performed on C57BL/6J and BALB/c mice. Four weeks later splenic CD11b+Ly6Chigh cells from both groups were isolated for phenotypic (flow cytometry) and functional (phagocytosis and glycolysis) characterization and RNA was obtained for single-cell RNA-seq (scRNA-seq) and subsequent analysis.ResultsCD11b+Ly6Chigh cells from sham and CLP surviving mice exhibit phenotypic and functional differences that relate to immune function, some of which are observed in both C57BL/6J and BALB/c strains and others that are not. To dissect disease-specific and strain-specific distinctions within the myeloid compartment, scRNA-seq analysis was performed on CD11b+Ly6Chigh cells from C57BL/6J and BALB/c sham and CLP mice. Uniform Manifold Approximation and Projection from both strains identified 13 distinct clusters of sorted CD11b+Ly6Chigh cells demonstrating significant transcriptional heterogeneity and expressing gene signatures corresponding to classical-monocytes, non-classical monocytes, M1- or M2-like macrophages, dendritic-like cells, monocyte-derived dendritic-like cells, and proliferating monocytic myeloid-derived suppressor cells (M-MDSCs). Frequency plots showed that the percentages of proliferating M-MDSCs (clusters 8, 11 and 12) were increased in CLP mice compared to sham mice in both strains. Pathway and UCell score analysis in CLP mice revealed that cell cycle and glycolytic pathways were upregulated in proliferating M-MDSCs in both strains. Notably, granule protease genes were upregulated in M-MDSCs from CLP mice. ScRNA-seq analyses also showed that phagocytic pathways were upregulated in multiple clusters including the classical monocyte cluster, confirming the increased phagocytic capacity in CD11b+Ly6Chigh cells from CLP mice observed in ex vivo functional assays in C57BL/6J mice.ConclusionThe splenic CD11b+Ly6Chigh myeloid populations expanded in survivors of CLP sepsis correspond to proliferating cells that have an increased metabolic demand and gene signatures consistent with M-MDSCs, a population known to have immunosuppressive capacity.

1H NMR-Based Metabolomic Profiling and Comparison of Human Milk Across Different Lactation Stages in Secretors and Non-Secretors: A Study of Chinese Lactating Women

BackgroundHuman milk oligosaccharides (HMOs) and other milk-derived metabolites are crucial for infant health, influencing gut microbiota and overall development. ObjectiveThis study aims to uncover insights into the variations of HMOs and non-HMO metabolites based on secretor (Se) status, lactation time, mode of delivery, and infant sex. MethodsAn exploratory cross-sectional study was designed to compare the levels of HMOs and non-HMOs metabolites in milk samples from 129 lactating Chinese women within 1 year postpartum. Nuclear magnetic resonance analysis was employed for the identification and quantification of the metabolites. The metabolites measured were grouped into sugars, free amino acids, fatty acids, and metabolites related to energy metabolism. The influence of delivery mode and infant sex on milk metabolite composition were explored. ResultsUniform Manifold Approximation and Projection (UMAP) analysis of HMOs profiles revealed distinct clustering based on Se status, with significant differences in 2'-FL and 3-FL levels observed between Se+ and Se- groups. A decreasing trend for 2'-FL and 6-’SL levels, along with an increase in 3-FL levels, was observed with increasing lactating period within 12 months postpartum. Non-HMOs metabolite analysis indicated that Se status only affected glutamate levels. An increase in glutamine levels was observed 3–9 months postpartum. A continuous increase in o-phosphocholine levels was noted in 12 months postpartum, along with reductions in citrate and sn-glycero-phosphocholine levels. Delivery mode and infant sex did not affect both HMOs and non-HMOs levels. ConclusionsMetabolomic analysis of human milk reveals significant variation of HMOs, but not in non-HMOs, based on Se status. Changes in certain HMOs and non-HMOs levels were also observed over the one year of lactation. Understanding how these metabolites change over time may influence recommendations for maternal diet, supplementation, and the timing of breastfeeding to ensure optimal nutrient delivery to the infant.

Signature analysis of high-throughput transcriptomics screening data for mechanistic inference and chemical grouping.

High-throughput transcriptomics (HTTr) uses gene expression profiling to characterize the biological activity of chemicals in in vitro cell-based test systems. As an extension of a previous study testing 44 chemicals, HTTr was used to screen an additional 1,751 unique chemicals from the EPA's ToxCast collection in MCF7 cells using 8 concentrations and an exposure duration of 6 h. We hypothesized that concentration-response modeling of signature scores could be used to identify putative molecular targets and cluster chemicals with similar bioactivity. Clustering and enrichment analyses were conducted based on signature catalog annotations and ToxPrint chemotypes to facilitate molecular target prediction and grouping of chemicals with similar bioactivity profiles. Enrichment analysis based on signature catalog annotation identified known mechanisms of action (MeOAs) associated with well-studied chemicals and generated putative MeOAs for other active chemicals. Chemicals with predicted MeOAs included those targeting estrogen receptor (ER), glucocorticoid receptor (GR), retinoic acid receptor (RAR), the NRF2/KEAP/ARE pathway, AP-1 activation, and others. Using reference chemicals for ER modulation, the study demonstrated that HTTr in MCF7 cells was able to stratify chemicals in terms of agonist potency, distinguish ER agonists from antagonists, and cluster chemicals with similar activities as predicted by the ToxCast ER Pathway model. Uniform manifold approximation and projection (UMAP) embedding of signature-level results identified novel ER modulators with no ToxCast ER Pathway model predictions. Finally, UMAP combined with ToxPrint chemotype enrichment was used to explore the biological activity of structurally related chemicals. The study demonstrates that HTTr can be used to inform chemical risk assessment by determining in vitro points of departure, predicting chemicals' MeOA and grouping chemicals with similar bioactivity profiles.

Advanced classification of hot subdwarf binaries using artificial intelligence techniques and Gaia DR3 data

Context. Hot subdwarf stars are compact blue evolved objects, burning helium in their cores surrounded by a tiny hydrogen envelope. In the Hertzsprung-Russell Diagram they are located by the blue end of the Horizontal Branch. Most models agree on a quite probable common envelope binary evolution scenario in the Red Giant phase. However, the current binarity rate for these objects is yet unsolved, but key, question in this field. Aims. This study aims to develop a novel classification method for identifying hot subdwarf binaries within large datasets using Artificial Intelligence techniques and data from the third Gaia data release (GDR3). The results will be compared with those obtained previously using Virtual Observatory techniques on coincident samples. Methods. The methods used for hot subdwarf binary classification include supervised and unsupervised machine learning techniques. Specifically, we have used Support Vector Machines (SVM) to classify 3084 hot subdwarf stars based on their colour-magnitude properties. Among these, 2815 objects have Gaia DR3 BP/RP spectra, which were classified using Self-Organizing Maps (SOM) and Convolutional Neural Networks (CNN). In order to ensure spectral quality, previously to SOM and CNN classification, our 2815 BP/RP set were pre-analysed with two different approaches: the cosine similarity technique and the Uniform Manifold Approximation and Projection (UMAP) technique. Additional analysis onto a golden sample of 88 well-defined objects, is also presented. Results. The findings demonstrate a high agreement level (∼70–90%) with the classifications from the Virtual Observatory Sed Analyzer (VOSA) tool. This shows that the SVM, SOM, and CNN methods effectively classify sources with an accuracy comparable to human inspection or non-AI techniques. Notably, SVM in a radial basis function achieves 70.97% reproducibility for binary targets using photometry, and CNN reaches 84.94% for binary detection using spectroscopy. We also found that the single–binary differences are especially observable on the infrared flux in our Gaia DR3 BP/BR spectra, at wavelengths larger than ∼700 nm. Conclusions. We find that all the methods used are in fairly good agreement and are particularly effective to discern between single and binary systems. The agreement is also consistent with the results previously obtained with VOSA. In global terms, considering all quality metrics, CNN is the method that provides the best accuracy. The methods also appear effective for detecting peculiarities in the spectra. While promising, challenges in dealing with uncertain compositions highlight the need for caution, suggesting further research is needed to refine techniques and enhance automated classification reliability, particularly for large-scale surveys.

Pan-Cancer Single-Cell Transcriptomic Analysis Reveals Divergent Expression of Embryonic Proangiogenesis Gene Modules in Tumorigenesis.

Angiogenesis is indispensable for the sustained survival and progression of both embryonic development and tumorigenesis. This intricate process is tightly regulated by a multitude of pro-angiogenic genes. The presence of gene modules facilitating angiogenesis has been substantiated in both embryonic development and the context of tumor proliferation. However, it remains unresolved whether the pro-angiogenic gene modules expressed during embryonic development also exist in tumors. This study performed a pan-cancer single-cell RNA sequencing (scRNA-seq) analysis on samples from 332 patients across seven cancer types: thyroid carcinoma, lung cancer, breast cancer, hepatocellular carcinoma, colorectal cancer, ovarian carcinoma, and prostate adenocarcinoma. Data processing was carried out using the Seurat R package, with rigorous quality control to filter high-quality cells and mitigate batch effects across datasets. We used principal component analysis (PCA), shared nearest neighbor graph-based clustering, and Uniform Manifold Approximation and Projection (UMAP) to visualize cell types and identify distinct cell clusters. Myeloid cell subpopulations were further analyzed for the expression of embryonic pro-angiogenic gene modules (EPGM) and tumor pro-angiogenic gene modules (TPGM). The analysis identified nine major cell types within the tumor microenvironment, with myeloid cells consistently exhibiting elevated expression of both tumor pro-angiogenic gene modules (TPGM) and EPGM across all tumor types. In particular, myeloid cells, including macrophages and monocytes, showed high EPGM expression, indicating an active role of embryonic pro-angiogenesis pathways in tumors. A subset analysis revealed 20 distinct myeloid subtypes with varying EPGM and TPGM expression across different cancers. Treatment and disease stage influenced these gene expressions, with certain subtypes, such as HSPAhi/STAT1+ macrophages in breast cancer, displaying reduced pro-angiogenic gene activity post-treatment. This study provides evidence that tumors may exploit EPGM to enhance vascularization and support sustained growth, as evidenced by the elevated EPGM expression in tumor-associated myeloid cells. The consistent presence of EPGM in TAMs across multiple cancer types suggests a conserved mechanism wherein tumors harness embryonic angiogenic pathways to facilitate their progression. Distinct EPGM expression patterns in specific myeloid cell subsets indicate potential therapeutic targets, particularly in cases where EPGM activation contributes to resistance against anti-angiogenic therapies. These findings shed new light on the molecular mechanisms underlying tumor angiogenesis and highlight the prognostic relevance of EPGM expression in cancer, underscoring its potential as a biomarker for clinical applications.