Search NASA⌕ Search

SEARCH · Search NASA

Results for “RNA-sequencing”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Single‐Cell Nanodroplet Processing Proteomics Pipeline for Analysis of Human‐Derived Microglia

Single-cell omics tools provide unique insights into heterogeneous cell populations and their responses to stimuli. For example, single-cell RNA sequencing has identified several transcriptionally distinct populations of microglia, which are resident immune cells of the central nervous system (CNS) that are responsive to CNS injury, infection, and neurodegeneration. To date, single-cell studies of microglia have focused on RNA-sequencing or cytometry by time of flight (CyTOF), which provide indirect readouts of protein abundance or quantification of a limited number of targets. Herein, we present a workflow based on FACS-assisted isolation, cryopreservation, and nanodroplet-based processing for single-cell mass spectrometry proteomics analysis of the postmortem human brain cortex-derived microglia. From a single microglial cell, 1039 proteins could be identified on average. As a proof-of-principle, we applied single-cell proteomics for exploring the heterogeneity of brain microglia at the cellular level. This pilot proteomics data partially recapitulates the prior microglia subtypes. Specifically, we determined that mitochondrial proteins, in particular members of NADH dehydrogenase (Complex I), cytochrome b-c1 (Complex III), cytochrome c oxidase (Complex IV), F1-ATPase (Complex V), and Na+/K+-ATPase complex, drive variation across microglia. This pipeline offers the potential for identifying functionally and analytically relevant protein targets for microglia in Alzheimer's disease and other neurological disorders.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Inference of cell type-specific gene regulatory networks on cell lineages from single cell omic datasets

Abstract Cell type-specific gene expression patterns are outputs of transcriptional gene regulatory networks (GRNs) that connect transcription factors and signaling proteins to target genes. Single-cell technologies such as single cell RNA-sequencing (scRNA-seq) and single cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq), can examine cell-type specific gene regulation at unprecedented detail. However, current approaches to infer cell type-specific GRNs are limited in their ability to integrate scRNA-seq and scATAC-seq measurements and to model network dynamics on a cell lineage. To address this challenge, we have developed single-cell Multi-Task Network Inference (scMTNI), a multi-task learning framework to infer the GRN for each cell type on a lineage from scRNA-seq and scATAC-seq data. Using simulated and real datasets, we show that scMTNI is a broadly applicable framework for linear and branching lineages that accurately infers GRN dynamics and identifies key regulators of fate transitions for diverse processes such as cellular reprogramming and differentiation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Wox3 -patterning module organizes planar growth in grass leaves and ligules

Grass leaves develop from a ring of primordial initial cells within the periphery of the shoot apical meristem, a pool of organogenic stem cells that generates all of the organs of the plant shoot. At maturity, the grass leaf is a flattened, strap-like organ comprising a proximal supportive sheath surrounding the stem and a distal photosynthetic blade. The sheath and blade are partitioned by a hinge-like auricle and the ligule, a fringe of epidermally derived tissue that grows from the adaxial (top) leaf surface. Together, the ligule and auricle comprise morphological novelties that are specific to grass leaves. Understanding how the planar outgrowth of grass leaves and their adjoining ligules is genetically controlled can yield insight into their evolutionary origins. Here we use single-cell RNA-sequencing analyses to identify a ‘rim’ cell type present at the margins of maize leaf primordia. Cells in the leaf rim have a distinctive identity and share transcriptional signatures with proliferating ligule cells, suggesting that a shared developmental genetic programme patterns both leaves and ligules. Moreover, we show that rim function is regulated by genetically redundant Wuschel-like homeobox3 (WOX3) transcription factors. Higher-order mutations in maize Wox3 genes greatly reduce leaf width and disrupt ligule outgrowth and patterning. Together, these findings illustrate the generalizable use of a rim domain during planar growth of maize leaves and ligules, and suggest a parsimonious model for the homology of the grass ligule as a distal extension of the leaf sheath margin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Discovery of FoTO1 and Taxol genes enables biosynthesis of baccatin III

Abstract Plants make complex and potent therapeutic molecules 1,2 , but sourcing these molecules from natural producers or through chemical synthesis is difficult, which limits their use in the clinic. A prominent example is the anti-cancer therapeutic paclitaxel (sold under the brand name Taxol), which is derived from yew trees (Taxusspecies) 3 . Identifying the full paclitaxel biosynthetic pathway would enable heterologous production of the drug, but this has yet to be achieved despite half a century of research 4 . WithinTaxus’ large, enzyme-rich genome 5 , we suspected that the paclitaxel pathway would be difficult to resolve using conventional RNA-sequencing and co-expression analyses. Here, to improve the resolution of transcriptional analysis for pathway identification, we developed a strategy we term multiplexed perturbation × single nuclei (mpXsn) to transcriptionally profile cell states spanning tissues, cell types, developmental stages and elicitation conditions. Our data show that paclitaxel biosynthetic genes segregate into distinct expression modules that suggest consecutive subpathways. These modules resolved seven new genes, allowing a de novo 17-gene biosynthesis and isolation of baccatin III, the industrial precursor to Taxol, inNicotiana benthamianaleaves, at levels comparable with the natural abundance inTaxusneedles. Notably, we found that a nuclear transport factor 2 (NTF2)-like protein, FoTO1, is crucial for promoting the formation of the desired product during the first oxidation, resolving a long-standing bottleneck in paclitaxel pathway reconstitution. Together with a new β-phenylalanine-CoA ligase, the eight genes discovered here enable the de novo biosynthesis of 3’-N-debenzoyl-2’-deoxypaclitaxel. More broadly, we establish a generalizable approach to efficiently scale the power of co-expression analysis to match the complexity of large, uncharacterized genomes, facilitating the discovery of high-value gene sets.

Science & Technology - Other Topics↗
DOE OSTI2023

Robust reconstruction of single-cell RNA-seq data with iterative gene weight updates

Single-cell RNA-sequencing technologies have greatly enhanced our understanding of heterogeneous cell populations and underlying regulatory processes. However, structural (spatial or temporal) relations between cells are lost during cell dissociation. These relations are crucial for identifying associated biological processes. Many existing tissue-reconstruction algorithms use prior information about subsets of genes that are informative with respect to the structure or process to be reconstructed. When such information is not available, and in the general case when the input genes code for multiple processes, including being susceptible to noise, biological reconstruction is often computationally challenging. We propose an algorithm that iteratively identifies manifold-informative genes using existing reconstruction algorithms for single-cell RNA-seq data as subroutine. We show that our algorithm improves the quality of tissue reconstruction for diverse synthetic and real scRNA-seq data, including data from the mammalian intestinal epithelium and liver lobules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Streamlined spatial and environmental expression signatures characterize the minimalist duckweed Wolffia australiana

Single-cell genomics permits a new resolution in the examination of molecular and cellular dynamics, allowing global, parallel assessments of cell types and cellular behaviors through development and in response to environmental circumstances, such as interaction with water and the light–dark cycle of the Earth. Here, we leverage the smallest, and possibly most structurally reduced, plant, the semiaquaticWolffia australiana, to understand dynamics of cell expression in these contexts at the whole-plant level. We examined single-cell-resolution RNA-sequencing data and foundWolffiacells divide into four principal clusters representing the above- and below-water-situated parenchyma and epidermis. Although these tissues share transcriptomic similarity with model plants, they display distinct adaptations thatWolffiahas made for the aquatic environment. Within this broad classification, discrete subspecializations are evident, with select cells showing unique transcriptomic signatures associated with developmental maturation and specialized physiologies. Assessing this simplified biological system temporally at two key time-of-day (TOD) transitions, we identify additional TOD-responsive genes previously overlooked in whole-plant transcriptomic approaches and demonstrate that the core circadian clock machinery and its downstream responses can vary in cell-specific manners, even in this simplified system. Distinctions between cell types and their responses to submergence and/or TOD are driven by expression changes of unexpectedly few genes, characterizingWolffiaas a highly streamlined organism with the majority of genes dedicated to fundamental cellular processes.Wolffiaprovides a unique opportunity to apply reductionist biology to elucidate signaling functions at the organismal level, for which this work provides a powerful resource.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Transcriptomic data sets for Novosphingobium aromaticivorans grown with the β-5-linked aromatic dimer dehydrodiconiferyl alcohol and the related G-aromatic monomers vanillin and ferulic acid

ABSTRACT The transcriptomes of a 2-pyrone-4,6-dicarboxylic acid-producing strain of Novosphingobium aromaticivorans DSM12444 were determined when grown in minimal medium containing glucose alone or glucose plus vanillin, ferulic acid, or the β-5-linked aromatic dimer dehydrodiconiferyl alcohol as carbon sources. Here, we present the RNA-sequencing data we obtained.

Metz, Fletcher↗
DOE OSTI2024

The TLK-ASF1 histone chaperone pathway plays a critical role in IL-1β–mediated AML progression

Identifying and targeting microenvironment-driven pathways that are active across acute myeloid leukemia (AML) genetic subtypes should allow the development of more broadly effective therapies. The proinflammatory cytokine interleukin-1β (IL-1β) is abundant in the AML microenvironment and promotes leukemic growth. Through RNA-sequencing analysis, we identify that IL-1β–upregulated ASF1B (antisilencing function-1B), a histone chaperone, in AML progenitors compared with healthy progenitors. ASF1B, along with its paralogous protein ASF1A, recruits H3-H4 histones onto the replication fork during S-phase, a process regulated by Tousled-like kinase 1 and 2 (TLKs). Although ASF1s and TLKs are known to be overexpressed in multiple solid tumors and associated with poor prognosis, their functional roles in hematopoiesis and inflammation-driven leukemia remain unexplored. In this study, we identify that ASF1s and TLKs are overexpressed in multiple genetic subtypes of AML. We demonstrate that depletion of ASF1s significantly reduces leukemic cell growth in both in vitro and in vivo models using human cells. Using a murine model, we show that overexpression of ASF1B accelerates leukemia progression. Moreover, Asf1b or Tlk2 deletion delayed leukemia progression, whereas these proteins are dispensable for normal hematopoiesis. Through proteomics and phosphoproteomics analyses, we uncover that the TLK-ASF1 pathway promotes leukemogenesis by affecting the cell cycle and DNA damage pathways. Collectively, our findings identify the TLK1-ASF1 pathway as a novel mediator of inflammatory signaling and a promising therapeutic target for AML treatment across diverse genetic subtypes. Finally, selective inhibition of this pathway offers potential opportunities to intervene effectively, address intratumoral heterogeneity, and ultimately improve clinical outcomes in AML.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Single-cell proteomics of Arabidopsis leaf mesophyll reveals dynamic protein responses to water-deficit stress

Background The application of single-cell omics tools to biological systems can provide unique insights into diverse cellular populations and their heterogeneous responses to internal and external perturbations. Thus far, most single-cell studies in plant systems have been limited to RNA-sequencing approaches, which only provide indirect readouts of cellular functions. Results Here, we present a single-cell proteomics workflow for plant cells that integrates tape-sandwich protoplasting, piezoelectric cell sorting, nanoPOTS sample preparation, and ion mobility-based MS data acquisition method for label-free single-cell proteomics analysis of Arabidopsis leaf mesophyll cells. From a single leaf protoplast, over 3,000 proteins were quantified with high precision. The workflow is demonstrated to identify stress associated changes in protein abundance by analyzing 117 protoplasts from well-watered and water-deficit stressed plants. Additionally, we describe a new approach for constructing covarying protein networks at the single-cell level and demonstrate how single-cell protein covariation analysis can reveal previously unrecognized protein functions while also capturing stress-induced changes in protein–protein dynamics. Conclusions The label-free scProteomic approach presented here represents a significant advance through the demonstration of a facile protoplast isolation method combined with deep and precise proteomic coverage of Arabidopsis leaf mesophyll cell types. We believe this study will serve as an informative reference to future plant scProteomic investigations.

Arabidopsis↗
DOE OSTI2022

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus

ABSTRACT Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials – a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Combined Imaging and RNA-Seq on a Microfluidic Platform for Viral Infection Studies

The goal of this work was to pioneer a novel, low-overhead protocol for simultaneously assaying cell-surface markers and intracellular gene expression in a single mammalian cell. The purpose of developing such a method is to be able to understand the mechanisms by which pathogens engage with individual mammalian cells, depending on their cell surface proteins, and how both host and pathogen gene expression changes are reflective of these mechanisms. The knowledge gained from such analyses of single cells will ultimately lead to more robust pathogen detection and countermeasures. Our method was aimed at streamlining both the upstream cell sample preparation using microfluidic methods, as well as the actual library making protocol. Specifically, we wanted to implement a random hexamer-based reverse transcription of all RNA within a single cell (as opposed to oligo dT-based which would only capture polyadenylated transcripts), and then use a CRISPR-based method called scDash to deplete ribosomal DNAs (since ribosomal RNAs make up the majority of the RNA in a mammalian cell). After significant troubleshooting, we demonstrate that we are able to prepare cDNA from RNA using the random hexamer primer, and perform the rDNA depletion. We also show that we can visualize individually stained cells, setting up the pipeline for connecting surface markers to RNA-sequencing profiles. Finally, we test a number of devices for various parts of the pipeline, including bead generation, optical barcoding and cell dispensing, and demonstrate that while some of these have potential, more work is needed to optimize this part of the pipeline.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

3′ RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

3’ RNA-seq↗
DOE OSTI2024

O-GlcNAc transferase regulates collagen deposition and fibrosis resolution in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic pulmonary disease that is characterized by an excessive accumulation of extracellular matrix (ECM) proteins (e.g. collagens) in the parenchyma, which ultimately leads to respiratory failure and death. While current therapies exist to slow the progression, no therapies are available to resolve fibrosis. We characterized the O-linked N-Acetylglucosamine (O-GlcNAc) transferase (OGT)/O-GlcNAc axis in IPF using single-cell RNA-sequencing (scRNA-seq) data and human lung sections and isolated fibroblasts from IPF and non-IPF donors. The underlying mechanism(s) of IPF were further investigated using multiple experimental models to modulate collagen expression and accumulation by genetically and pharmacologically targeting OGT. Furthermore, we hone in on the transforming growth factor-beta (TGF-β) effector molecule, Smad3, by co-expressing it with OGT to determine if it is modified and its subsequent effect on Smad3 activation. We found that OGT and O-GlcNAc levels are upregulated in patients with IPF compared to non-IPF. We report that the OGT regulates collagen deposition and fibrosis resolution, which is an evolutionarily conserved process demonstrated across multiple species. Co-expression of OGT and Smad3 showed that Smad3 is O-GlcNAc modified. Blocking OGT activity resulted in decreased phosphorylation at Ser-423/425 of Smad3 attenuating the effects of TGF-β1 induced collagen expression/deposition. OGT inhibition or knockdown successfully blocked and reversed collagen expression and accumulation, respectively. Smad3 is discovered to be a substrate of OGT and its O-GlcNAc modification(s) directly affects its phosphorylation state. These data identify OGT as a potential target in pulmonary fibrosis resolution, as well as other diseases that might have aberrant ECM/collagen accumulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Genomic and phenotypic comparison of two variants of multidrug-resistant Salmonella enterica serovar Heidelberg isolated during the 2015–2017 multi-state outbreak in cattle

Salmonella enterica subspecies enterica serovar Heidelberg (Salmonella Heidelberg) has caused several multistate foodborne outbreaks in the United States, largely associated with the consumption of poultry. However, a 2015–2017 multidrug-resistant (MDR) Salmonella Heidelberg outbreak was linked to contact with dairy beef calves. Traceback investigations revealed calves infected with outbreak strains of Salmonella Heidelberg exhibited symptoms of disease frequently followed by death from septicemia. To investigate virulence characteristics of Salmonella Heidelberg as a pathogen in bovine, two variants with distinct pulse-field gel electrophoresis (PFGE) patterns that differed in morbidity and mortality during the multistate outbreak were genotypically and phenotypically characterized and compared. Strain SX 245 with PFGE pattern JF6X01.0523 was identified as a dominant and highly pathogenic variant causing high morbidity and mortality in affected calves, whereas strain SX 244 with PFGE pattern JF6X01.0590 was classified as a low pathogenic variant causing less morbidity and mortality. Comparison of whole-genome sequences determined that SX 245 lacked ~200 genes present in SX 244, including genes associated with the IncI1 plasmid and phages; SX 244 lacked eight genes present in SX 245 including a second YdiV Anti-FlhC(2)FlhD(4) factor, a lysin motif domain containing protein, and a pentapeptide repeat protein. RNA-sequencing revealed fimbriae-related, flagella-related, and chemotaxis genes had increased expression in SX 245 compared to SX 244. Furthermore, SX 245 displayed higher invasion of human and bovine epithelial cells than SX 244. These data suggest that the presence and up-regulation of genes involved in type 1 fimbriae production, flagellar regulation and biogenesis, and chemotaxis may play a role in the increased pathogenicity and host range expansion of the Salmonella Heidelberg isolates involved in the bovine-related outbreak.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Diverse PFAS produce unique transcriptomic changes linked to developmental toxicity in zebrafish

Per- and polyfluoroalkyl substances (PFAS) are a widespread and persistent class of contaminants posing significant environmental and human health concerns. Comprehensive understanding of the modes of action underlying toxicity among structurally diverse PFAS is mostly lacking. To address this need, we recently reported on our application of developing zebrafish to evaluate a large library of PFAS for developmental toxicity. In the present study, we prioritized 15 bioactive PFAS that induced significant morphological effects and performed RNA-sequencing to characterize early transcriptional responses at a single timepoint (48 h post fertilization) after early developmental exposures (8 h post fertilization). Internal concentrations of 5 of the 15 PFAS were measured from pooled whole fish samples across multiple timepoints between 24–120 h post fertilization, and additional temporal transcriptomics at several timepoints (48–96 h post fertilization) were conducted for Nafion byproduct 2. A broad range of differentially expressed gene counts were identified across the PFAS exposures. Most PFAS that elicited robust transcriptomic changes affected biological processes of the brain and nervous system development. While PFAS disrupted unique processes, we also found that similarities in some functional head groups of PFAS were associated with the disruption in expression of similar gene sets. Body burdens after early developmental exposures to select sulfonic acid PFAS, including Nafion byproduct 2, increased from the 24–96 h post fertilization sampling timepoints and were greater than those of sulfonamide PFAS of similar chain lengths. In parallel, the Nafion byproduct 2-induced transcriptional responses increased between 48 and 96 h post fertilization. PFAS characteristics based on toxicity, transcriptomic effects, and modes of action will contribute to further prioritization of PFAS structures for testing and informed hazard assessment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

OxyR regulates the oxidative stress response in Zymomonas mobilis during oxic growth and anoxic biofuel fermentation

The bacterium Zymomonas mobilis is widely studied for its potential as an industrial biofuel producer. Anoxic fermentation by Z. mobilis in lignocellulosic hydrolysate can generate bioethanol from renewable plant biomass. In this study, we deleted a gene from the Z. mobilis genome encoding a homolog of OxyR, a transcription factor that activates an oxidative stress response in bacteria to reduce reactive oxygen species (ROS). Deletion of this transcription factor inhibited growth of Z. mobilis in oxic, but not anoxic, conditions in laboratory media. RNA-Sequencing was perfromed on wild-type Z. mobilis (ZM4) and ∆oxyR in both oxic and anoxic conditions in rich media (ZRMG). This study reveals the indirect regulon of OxyR in Z. mobilis, which is important for both oxic growth and anoxic biofuel fermentation.

bioenergy↗
DOE OSTI2024

Investigating biological nitrogen fixation via single-cell transcriptomics

The extensive use of nitrogen fertilizers has detrimental environmental consequences, and it is essential for society to explore sustainable alternatives. One promising avenue is engineering root nodule symbiosis, a naturally occurring process in certain plant species within the nitrogen-fixing clade, into non-leguminous crops. Advancements in single-cell transcriptomics provide unprecedented opportunities to dissect the molecular mechanisms underlying root nodule symbiosis at the cellular level. This review summarizes key findings from single-cell studies in Medicago truncatula, Lotus japonicus, and Glycine max. We highlight how these studies address fundamental questions about the development of root nodule symbiosis, including the following findings: (i) single-cell transcriptomics has revealed a conserved transcriptional program in root hair and cortical cells during rhizobial infection, suggesting a common infection pathway across legume species; (ii) characterization of determinate and indeterminate nodules using single-cell technologies supports the compartmentalization of nitrogen fixation, assimilation, and transport into distinct cell populations; (iii) single-cell transcriptomics data have enabled the identification of novel root nodule symbiosis genes and provided new approaches for prioritizing candidate genes for functional characterization; and (iv) trajectory inference and RNA velocity analyses of single-cell transcriptomics data have allowed the reconstruction of cellular lineages and dynamic transcriptional states during root nodule symbiosis.

Lotus japonicus↗