Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Dark Energy Survey Year 6 Results: MagLim++ Lens Sample Selection and Measurements of Galaxy Clustering

Galaxy clustering is a sensitive probe of the expansion history and growth of structure of the universe, and key degeneracies can be broken by combining these data with measurements of cosmic shear and galaxy-galaxy lensing (a so-called 3$\times$2pt analysis). The largest and least biased statistical samples of galaxies for use in clustering analyses can be collected photometrically through large imaging surveys. However, selecting clean photometric subsamples for cosmology are crucial for avoiding contamination that can bias cosmological constraints. Here we present the MagLim++ galaxy sample, selected to optimize for cosmological constraining power and incorporating an array of novel quality cuts to identify and remove residual contamination. This sample comes from the full six years of observations from the Dark Energy Survey. We present measurements of the two-point angular clustering ($w(θ)$) of 9,186,205 galaxies distributed over 4031 sq. degrees and in six tomographic redshift bins centered at $\bar{z}\approx$ [0.31, 0.44, 0.62, 0.78, 0.90, 1.01]. These measurements are used as part of the 3$\times$2pt and other DES Y6 legacy cosmological analyses in companion works. We describe the battery of null tests and mitigation schemes implemented to address observational, astrophysical, and methodological systematics in the analysis. The resulting $w(θ)$ measurements have a S/N = 149 (90.2 for linear scales only), which we use to place galaxy-clustering-only constraints on the matter density of the Universe, $Ω_m=0.311^{+0.023}_{-0.035}$, and amplitude of galaxy clustering in each redshift bin, $b_iσ_8=[1.16^{+0.04}_{-0.06},\ 1.40^{+0.04}_{-0.06},\ 1.57^{+0.04}_{-0.06},\ 1.59^{+0.04}_{-0.05},\ 1.50^{+0.04}_{-0.05},\ 1.74^{+0.06}_{-0.08}]$.

Weaverdyck, N. [BCCP, Berkeley; UC, Berkeley; LBL,↗
DOE OSTI2026

Multiresolution Quantum Chemistry: Nonlinear Response Properties at the Basis Set Limit

We benchmark the accuracy of Dunning correlation-consistent Gaussian basis sets for computing frequencydependent second-order hyperpolarizabilities relevant to second-harmonic generation (SHG), using multiresolution analysis (MRA) as a reference. Basis set errors are analyzed using a unit-sphere representation of the effective hyperpolarizability vector, enabling direct assessment of directional error structure. We introduce a relative RMS total error metric that integrates directional deviations over the unit sphere and complement it with signed projection errors that distinguish over- and underestimation. Unsupervised clustering based on these signed directional metrics reveals four distinct convergence behaviors across a set of 68 molecules. Unitsphere visualizations of representative systems show that basis set errors are often highly anisotropic and localized along specific bond directions, even when global error measures appear small. Doubly augmented basis sets consistently outperform singly augmented ones, and core-polarization functions are required for uniform convergence in second-row systems. Overall, this work demonstrates that directional analysis combined with clustering provides a robust framework for understanding basis set convergence in nonlinear optical response properties.

Basis sets↗
DOE OSTI2025

Predicting Dynamic-to-Static Correction Factor from Petrophysical Data and Chemostratigraphy using Unsupervised Machine Learning

Estimating static mechanical properties of stratigraphic layers is critical for optimizing subsurface engineering applications. To estimate dynamic-to-static correction factor F ds (static-to-dynamic Young’s modulus ratio) across the Caney shale interval in Oklahoma, USA, we integrated triaxial test measurements and petrophysical data, including well logs and X-ray fluorescence (XRF) using unsupervised machine learning (ML). We used a novel workflow that includes principal component analysis (PCA) to reduce data set dimensionality of well logs and XRF data sets—both separately and combined—creating three scenarios, and later applied inverse distance weighting (IDW) to derive F ds profiles for these scenarios. Furthermore, we applied K-means clustering on each scenario to predict depositional facies, and built a stiffness zonation profile through chemostratigraphic analysis of the terrigenous elements to validate the predicted F ds . The predicted F ds profile from each scenario using the PCA-IDW method was compared with the constant F ds approach from our previous study by calculating the root mean square error (RMSE). The combined data sets scenario yielded the lowest RMSE value of 0.113, while the RMSE values for the well logs and XRF scenarios were 0.131 and 0.129, respectively. In addition, the predicted F ds from the XRF scenario well-matched the stiffness zonation from the chemostratigraphic analysis that was built using the optimized K-means clustering of nine clusters for that scenario. These methods and findings offer a valuable tool for refining lithological classification and improving the F ds profile, potentially enhancing drilling and stimulation strategies for subsurface energy engineering applications.

clastic rock↗
DOE OSTI2024

Exploring Architectural-Aware Affinity Policies in Modern HPC Runtimes

Modern commodity and High-Performance Computing (HPC) systems are evolving with complex CPU architectures. These architectures now feature higher core and NUMA domain counts and implement features such as hyperthreading. When considering significant differences in hardware configurations, library availability, and hardware-tailored system/software stacks, which could substantially vary from one system to another, performance portability is hard to achieve. Throughout the years, this trend resulted in an increasingly high burden on application developers to fine-tune their workloads for each architecture. This work explores how hardware-dependent aspects such as locality/process/thread affinity affect performance in modern CPU architectures. We focus our study on the Global Memory and Threading (GMT) distributed runtime system as a representative of Partitioned Global Address Space (PGAS) software stacks commonly adopted for productivity. In particular, to appreciate performance implications, we evaluate GMT’s thread affinity policies, and, introduce two new ones which exploit architectural awareness. Finally, we explore alternative NUMA configurations via different process bindings and perform a scalability study on three HPC clusters with varying CPU architectures and NUMA layouts. Our analysis indicates that more complex architectures are more affected by affinity and binding policies and highlights the importance of setting proper runtime configurations to achieve superior performance.

Di Dio Lavore, Ian↗
DOE OSTI2025

Evaluating HPC Scheduling Strategies for Urgent Workloads

Scientific computing centers increasingly face workloads with diverse urgency requirements, driven by applications that demand rapid or even immediate execution. Appropriately configured scheduling policies can significantly improve both user satisfaction and overall cluster utilization. In this work, we present a systematic analysis of scheduler configurations under scenarios where a fraction of jobs have urgent computing needs. We evaluate multiple job scheduling simulators, develop a lightweight job-submission emulation framework, and create tools to analyze and visualize the resulting scheduling data. Our study identifies key trade-offs between responsiveness, fairness, and efficiency, and offers a set of practical scheduling configurations (particularly for Slurm) that can be tailored to HPC environments supporting mixed-urgency workloads.

Maheshwari, Ketan [ORNL] (ORCID:000000033800662X)↗
DOE OSTI2026

Probing the Coherence of Metal Absorbers toward Background Quasar Pairs

Abstract We investigate the spatial coherence of metal absorption lines in the circumgalactic medium using 4115 quasar pairs from Sloan Digital Sky Survey DR16. We identify 184 Mg ii (15 with dual detections) and 50 C iv (11 with dual detections) absorber pairs. The fraction of pairs with fractional equivalent width difference δW r ≤ 0.5 declines with projected separation: for Mg ii , it drops from ∼52% at 0–25 kpc to ∼7% at 100–200 kpc, while for C iv , it decreases from ∼100% to ∼21% over 0–200 kpc. In every separation bin, C iv absorbers with δW r ≤ 0.5 maintain a higher fraction than Mg ii absorbers. Similar results are obtained when the nondetections are treated as 2 σ W r upper limits in a survival analysis. The transverse autocorrelation function shows strong Mg ii clustering within 0–25 kpc, whereas the C iv excess remains relatively flat out to ∼100 kpc. These findings suggest that high-ionization C iv gas retains coherence over larger transverse scales than low-ionization Mg ii gas, in agreement with previous lensed quasar and high-resolution studies.

Li, Jing [Guangxi University] (ORCID:0009000148751↗
DOE OSTI2025

NiAl–MoO 2 S 2 Nanoparticles: Structural Evolution and Mechanistic Insights into High-Performance Selenium Oxyanion Removal across Diverse pH Conditions

Advancing sorbent materials for the selective removal of toxic oxyanions from water requires synthetic control, tunable chemistry, and an atomic-level understanding of structure–function relationships. Here, we report the synthesis and detailed characterization of NiAl–MoO 2 S 2 , a novel layered double hydroxide (LDH) nanomaterial designed for the efficient sequestration of selenium oxoanions (SeO 3 2– and SeO 4 2– ) from complex aqueous environments. The material is synthesized through a room-temperature ion-exchange process, wherein interlayer NO 3 – anions in NiAl–LDH are replaced with MoO 2 S 2 2– clusters, forming high-surface-area, flower-like nanoparticles. Comprehensive structural analysis using the synchrotron X-ray pair distribution function, X-ray absorption spectroscopy, and X-ray photoelectron spectroscopy reveals a distinct chemical transformation of intercalated [MoO 2 S 2 ] 2– into [Mo 2 O 2 S 6 ] 2– -like clusters, generating redox-active interlayers that drive selenium capture. This tailored interfacial chemistry underpins the material’s exceptional sorption performance, achieving distribution coefficients (K d ) ≥ 10 6 mL/g and maximum capacities of 343 mg/g for SeO 4 2– and 514 mg/g for SeO 3 2– , outperforming state-of-the-art inorganic sorbents. Importantly, NiAl–MoO 2 S 2 maintains high selectivity and capacity across acidic, neutral, and alkaline pH, efficiently removing selenium from ppm to sub-10 ppb trace levels, even in the presence of competing ions typical of natural and industrial waters. The selenium uptake proceeds via reductive precipitation coupled with the oxidation of molybdenum and sulfide within the LDH framework. This study highlights the power of strategic synthetic modification and interlayer functionalization in LDHs to unlock new structural motifs and redox chemistries, offering a scalable route to advanced materials for environmental remediation.

Adsorption↗
DOE OSTI2025

Cosmological constraints from the Planck cluster catalogue with DES shear profiles and Chandra observations

We present cosmological constraints from the Planck PSZ2 cosmological cluster sample, using weak-lensing shear profiles from Dark Energy Survey (DES) data and X-ray observations from the Chandra telescope for the mass calibration. We compute hydrostatic mass estimates for all clusters in the PSZ2 sample with a scaling relation between their Sunyaev-Zeldovich signal and X-ray derived hydrostatic mass, calibrated with the Chandra data. We introduce a method to correct these masses with a hydrostatic mass bias using shear profiles from wide-field galaxy surveys. We simultaneously fit the number counts of the PSZ2 sample and the mass calibration with the DES data, finding $Ω_\text{m}=0.312^{+0.018}_{-0.024}$, $σ_8=0.777\pm 0.024$, $S_8\equiv σ_8 \sqrt{Ω_\text{m} / 0.3}=0.791^{+0.023}_{-0.021}$, and $(1-b)=0.844^{+0.055}_{-0.062}$ for our baseline analysis when combined with BAO data. When considering a hydrostatic mass bias evolving with mass, we find $Ω_\text{m}=0.353^{+0.025}_{-0.031}$, $σ_8=0.751\pm 0.023$, and $S_8=0.814^{+0.019}_{-0.020}$. We verify the robustness of our results by exploring a variety of analysis settings, with a particular focus on the definition of the halo centre used for the extraction of shear profiles. We compare our results with a number of other analyses, in particular two recent analyses of cluster samples obtained from SPT and eROSITA data that share the same mass calibration data set. We find that our results are in overall agreement with most late-time probes, in very mild tension with CMB results (1.6$σ$), and in significant tension with results from eROSITA clusters (2.9$σ$). We confirm that our mass calibration is consistent with the eROSITA analysis by comparing masses for clusters present in both Planck and eROSITA samples, eliminating it as a potential cause of tension.

Aymerich, G. [Orsay, IAS; AIM, Saclay] (ORCID:0009↗
DOE OSTI2025

The Many-Body Expansion for Metals: II. Nonadditive Terms in Clusters Composed of Metals with n s 1 , n s 2 , and n s 2 p 1 Configurations

The many-body expansion (MBE) was applied to homometallic and heterometallic trimers of metals with ns 1 , ns 2 , and ns 2 p 1 configurations to investigate its convergence, the magnitude and nature (stabilizing/destabilizing) of the individual terms and seek an understanding of their variation across the different families of clusters. In particular, we examined the series of alkali metals (Li 3 , Na 3 , K 3 , Li 2 Na, LiNa 2 ), alkali metal borides (Li 2 B and LiB 2 ), and alkaline earth metals (Be 3 , Be 2 Mg, BeMg 3 , and Mg 3 ) trimers, as well as sodium clusters Na n , n = 2–5. Here, we found that there is no uniform contribution (stabilizing or destabilizing) across the series in the different families of trimers. For instance, the 2-B term stabilizes the ground states of the Na 3 (doublet), Na 4 (singlet), and Na 5 (doublet) clusters, and the 3-B term destabilizes them; however, the opposite holds for the quartet state of the Li 3 , Li 2 Na, LiNa 2 , and Na 3 clusters (destabilizing 2-B, stabilizing 3-B). Substituting Li with B in the quartet state of Li 3 results in a significant reduction of the 3-B term amounting to 16% (Li 2 B) and 5% (LiB 3 ) of the binding energy. On the contrary, the ground states of the alkaline earth metal clusters (Be 3 , Be 2 Mg, BeMg 3 , and Mg 3 ) are stabilized by the 3-B term, while the 2-B term destabilizes them. Overall, we find that the 3-B terms significantly stabilize the high-spin multiplicity states of the ns 1 configurations and the low-spin states of the ns 2 configurations. Finally, as the size of the metal increases, the contribution of the 3-B term to the binding energy decreases due to the longer metal–metal bond distances.

Alkali metals↗
DOE OSTI2025

Microstructurally informed synchrotron x-ray analysis revealing helium defect transitions in ultrafine grained tungsten

The formation of insoluble gaseous defects in materials due to nuclear transmutation or ion implantation involves the diffusion of impurity atoms to form atomic defect clusters that coalesce into bubbles or cavities and ultimately degrade the material properties. Transmission electron microscopy (TEM) is limited in its ability to resolve sub-nanometer gas clusters whereas X-ray diffraction (XRD) provides information pertaining to local atomic changes. Here, in this study, helium (He) implanted ultrafine grained tungsten is explored through a multimodal defect characterization campaign combining TEM-informed Small Angle X-ray Scattering (SAXS) analysis, XRD lattice parameter measurements, and nanoscale He cluster quantification from a region of reciprocal space accessible via Wide Angle X-ray Scattering (WAXS). Moderate elevated temperature implantations are shown to produce high concentrations of sub-nanoscale He clusters and small, homogeneously distributed cavities, which collectively are linked to lattice expansion and further substantiated through complementary atomistic simulations. Increased implantation temperatures encourage the diffusion of these defects to the grain boundaries (GBs), leading to lattice relaxation and the growth of large GB cavities manifesting as bimodal size distributions in the SAXS analysis. Overall, our results demonstrate the utility of multimodal synchrotron X-ray analysis in bridging the gap between microscale He cavity quantification and atomic-scale defect analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tea ( Camellia sinensis ) Extract-Mediated Green Synthesis of Co 3 O 4 and Co 3 O 4 @Graphene Nanocomposites for Multifunctional Applications in Pollutant Degradation, Sensing, and Energy Storage

A novel solution-mixing method was proposed to synthesize Co 3 O 4 /graphene nanocomposites (Co 3 O 4 @Gr) using a green tea leaf (Camellia sinensis) extract as the reducing agent. XRD analysis shows that the as-prepared Co 3 O 4 @Gr exhibits a cubic spinel crystal structure. From morphological analysis, the obtained Co 3 O 4 NS forms spherical clusters that are uniformly distributed on the graphene surface. FT-IR and Raman analyses confirmed the strong molecular and vibrational interactions between the Co 3 O 4 NS and Gr. The suppressed PL intensity peak of the Co 3 O 4 @Gr NCs indicated significant inhibition in the recombination of charge carriers between the hybrid orbitals within the composites. As a result, the catalytic efficiency of Co 3 O 4 @Gr NCs increased to 80% compared to pristine Co 3 O 4 , which exhibited only 45% efficiency against methylene blue (MB) dye. Moreover, the as-prepared NCs exhibited a detection limit of 0.01−224 μM, demonstrating a superior low-DPA detection with high sensitivity. The Co 3 O 4 @Gr/GCE exhibits admirable selectivity for various pesticides, fungicides, and metal ions, with outstanding reproducibility and stability. From electrochemical investigations, the highest specific capacitance values of the as-synthesized Co 3 O 4 @Gr were 349 F/g at a scan rate of 5 mV/s and 158 F/g at a current density of 1 A/g.

Capacitors↗
DOE OSTI2024

Formation and transformation of iron oxy-hydroxide precursor clusters to ferrihydrite

Iron (Fe) oxy-hydroxide minerals such as ferrihydrite (Fh) are ubiquitous in Earth-surface environments and important in biogeochemical element cycling. Recent research has suggested that their formation is preceded by the precipitation of ultrasmall (~1 nm) Keggin-like Fe oxy-hydroxide clusters. However, relatively little is understood about the structure of the precursor clusters and the impacts of pH and time on their growth and transformation to more stable phases. Here, we used a new method that involves mixed flow reactors (MFR) to synthesize these Fe oxy-hydroxide precursor clusters at pH 1.0, 1.5, 2.5, and 4.5. In situ and ex situ synchrotron scattering measurements and laboratory small-angle X-ray scattering (SAXS) were used to study the structure and size of Fe oxy-hydroxide clusters and their transformation products, respectively. Results show that with increasing pH, the particle size and structural order of samples increase, forming solids that resemble 2-line Fh at pH 4.5. The experimental data were compared with X-ray pair distribution functions (PDF) calculated for a range of Fe(III) oxyhydroxide clusters, including Fe 13 Keggin isomers computed previously using density functional theory (DFT), which yielded at best only partial agreement at short range (<5 Å). Aging of the clusters synthesized at pH 1.5 and 2.5 results in growth and transformation via Ostwald ripening to mixtures of goethite (Gt) and lepidocrocite (Lp). This process was inhibited by immediately reacting the early-formed clusters with phosphate (PO 4 3– ), suggesting that oxyanion surface complexes can stabilize the initial clusters by preventing growth and crystallization to more stable phases.

58 GEOSCIENCES↗
DOE OSTI2024

X-ray scattering based scanning tomography for imaging and structural characterization of cellulose in plants

X-ray and neutron scattering have long been used for structural characterization of cellulose in plants. Due to averaging over the illuminated sample volume, these measurements traditionally overlooked the compositional and morphological heterogeneity within the sample. Here, a scanning tomographic imaging method is described, using contrast derived from the X-ray scattering intensity, for virtually sectioning the sample to reveal its internal structure at a resolution of a few micrometres. This method provides a means for retrieving the local scattering signal that corresponds to any voxel within the virtual section, enabling characterization of the local structure using traditional data-analysis methods. This is accomplished through tomographic reconstruction of the spatial distribution of a handful of mathematical components identified by non-negative matrix factorization from the large dataset of X-ray scattering intensity. Joint analysis of multiple datasets, to find similarity between voxels by clustering of the decomposed data, could help elucidate systematic differences between samples, such as those expected from genetic modifications, chemical treatments or fungal decay. The spatial distribution of the microfibril angle can also be analyzed, based on the tomographically reconstructed scattering intensity as a function of the azimuthal angle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Advances in ArborX to support exascale applications

ArborX is a performance portable geometric search library developed as part of the Exascale Computing Project (ECP). In this paper, we explore a collaboration between ArborX and a cosmological simulation code HACC. Large cosmological simulations on exascale platforms encounter a bottleneck due to the in-situ analysis requirements of halo finding, a problem of identifying dense clusters of dark matter (halos). This problem is solved by using a density-based DBSCAN clustering algorithm. With each MPI rank handling hundreds of millions of particles, it is imperative for the DBSCAN implementation to be efficient. In addition, the requirement to support exascale supercomputers from different vendors necessitates performance portability of the algorithm. We describe how this challenge problem guided ArborX development, and enhanced the performance and the scope of the library. We explore the improvements in the basic algorithms for the underlying search index to improve the performance, and describe several implementations of DBSCAN in ArborX. Further, we report the history of the changes in ArborX and their effect on the time to solve a representative benchmark problem, as well as demonstrate the real world impact on production end-to-end cosmology simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multiprobe cosmology from the abundance of SPT clusters and DES galaxy clustering and weak lensing

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the Universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3 × 2 ⁢pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining Λ cold dark matter (Λ ⁢CDM) parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure Ω m = 0.300 ± 0.017 and 𝜎 8 = 0.797 ± 0.026. Compared to constraints from Planck primary cosmic microwave background (CMB) anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 (1.2⁢𝜎) for the two-parameter difference. We further obtain 𝑆 8 ≡𝜎 8 ⁢(Ω m /0.3) 0.5 = 0.796 ± 0.013 which is lower than the Planck measurement at the 1.6⁢𝜎 level. The combined SPT cluster, DES 3 ×2 ⁢pt, and Planck datasets mildly prefer a nonzero positive neutrino mass, with a 95% upper limit ∑ 𝑚 𝜈 < 0.25 eV on the sum of neutrino masses. Assuming a 𝑤⁢CDM model, we constrain the dark energy equation of state parameter 𝑤 = −1.1⁢5$^{+0.23}_{−0.17}$ and when combining with Planck primary CMB anisotropies, we recover 𝑤 = −1.2⁢0$^{+0.15}_{−0.09}$, a 1.7⁢𝜎 difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology and our analysis paves the way for upcoming joint analyses of next-generation datasets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Production of the light-activated elsinochrome phytotoxin in the soybean pathogen Coniothyrium glycines hints at virulence factor

The Dothideomycete pathogenConiothyrium glycinescauses red leaf blotch of soybean, a major disease in Africa. It is one of two fungal plant pathogens on the USDA PPQ Select Agents and Toxins list of pathogens important to the biosecurity of the United States, reflective of its potential to be highly destructive if introduced. Despite its importance, there are no published reports regarding the molecular basis of host infection. Examination of theC. glycinesgenome revealed a secondary metabolite gene cluster that is similar in gene content and organization to clusters that synthesize light-activated perylenequinone toxins, such as cercosporin. Perylenequinones are non-host specific toxins that, upon exposure to light, generate reactive oxygen species, which have near-universal toxicity to plant hosts.Coniothyrium glycinesisolates from eastern and southern Africa were cultured axenically under light and dark conditions. Light-grown cultures produced red-pink pigmentation typical of perylenequinones. Differential gene expression analysis showed that six of the eight genes in the biosynthetic gene cluster, including the polyketide synthase gene, were significantly upregulated in light. Liquid chromatography-mass spectrometry confirmed production of the perylenequinone elsinochrome A, a known virulence factor in other fungal pathogens. On leaves incubated in the dark, significantly fewer lesions formed and symptoms were delayed, compared to leaves incubated in the light. In addition, we identified orthologous gene clusters in more distantly related Dothideomycete plant pathogens where their presence was previously unknown, indicating a broader importance of these toxins to agriculture and fungal ecology. This work provides the first evidence that elsinochrome A may contribute to the virulence ofC. glycines.

Science & Technology - Other Topics↗
DOE OSTI2025

Single cell RNA sequencing reveals shifts in cell maturity and function of endogenous and infiltrating cell types in response to acute intervertebral disc injury

Intervertebral disc (IVD) degeneration contributes to disabling back pain. Degeneration can be initiated by injury and progressively leads to an irreversible loss of cells and function. IVD function restoration through cell replacement therapies have had limited success due to knowledge gaps in the critical cell populations important for repair. Here, in this study, we used single cell RNA sequencing to identify the transcriptional changes of IVD resident and infiltrating cell populations from Control and Injured coccygeal IVDs extracted from 12-week-old female C57BL/6J mice 7 days post injury. Clustering, gene ontology, and pseudotime trajectory analyses determined transcriptomic divergences with injury, flow cytometry identified they types of infiltrating immune cells, and immunofluorescence was utilized to define mesenchymal stem cell (MSC) localization. We identified 11 distinct clusters that included IVD, immune, vascular cells, and MSCs. Differential gene expression analysis determined that Outer Annulus Fibrosus, Neutrophils, Saa2-High MSCs, Macrophages, and Krt18 + Nucleus Pulposus (NP) cells were the major drivers of transcriptomic differences between Control and Injured cells. Gene ontology revealed that the most upregulated biological pathways were angiogenesis and T cell-related while wound healing and ECM regulation were downregulated. Pseudotime trajectory analyses revealed that IVD injury directed cells towards increased differentiation in all clusters, except for Krt18 + NP cells which remained in a less mature cell state. Saa2-High and Grem1-High MSCs populations shifted towards more differentiated IVD cells profiles with injury and localized distinctly within the IVD. This study revealed novel MSC populations with the potential to be leveraged for future IVD repair studies.

Cartilage↗