Search NASA⌕ Search

SEARCH · Search NASA

Results for “differential expression”

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 487 records · Page 27

Non-coding structural variants disrupt FOXG1 transcriptional regulation in early neurodevelopment

The FOXG1 transcription factor is a crucial regulator of embryonic brain development. Pathogenic FOXG1 variants cause FOXG1 syndrome. Although structural variants in the non-coding region downstream of FOXG1 have been reported in 38 individuals with similar characteristics, the regulatory pathomechanisms remain unknown. Here, we identify two non-coding structural variants in individuals with FOXG1 syndrome-like features, allowing us to delineate a ~ 124 kb commonly affected regulatory region. Using epigenomic profiling and in vivo enhancer assays, we characterize and validate regulatory elements within the commonly affected regulatory region and wider FOXG1 TAD. We see strong activation of previously validated forebrain enhancers, and identify an enhancer cluster and progenitor-specific enhancer region that are strongly activated during forebrain-directed neural progenitor cell differentiation, a process in which FOXG1 is an important regulator. Perturbation of these elements results in varying degrees of reduced FOXG1 transcription in forebrain neural progenitor cells and in population shifts within these cells, while removal of the TAD boundary leads to aberrant expression of the neighbouring PRKD1 gene. Our findings characterize enhancer and architectural elements essential for proper FOXG1 transcription during neurodevelopment, therefore improving variant interpretation in this region.

Hamerlinck, Lisa↗

The chromatin remodeler ADNP regulates neurodevelopmental disorder risk genes and neocortical neurogenesis

Although chromatin remodelers are among the most important risk genes associated with neurodevelopmental disorders (NDDs), the roles of these complexes during brain development are in many cases unclear. Here, we focused on the recently discovered ChAHP chromatin remodeling complex. The zinc finger and homeodomain transcription factor ADNP is a core subunit of this complex, and de novoADNPmutations lead to intellectual disability and autism spectrum disorder. However, germlineAdnpknockout mice were previously shown to exhibit early embryonic lethality, obscuring subsequent roles for the ChAHP complex in neurogenesis. To circumvent this early developmental arrest, we generated a conditionalAdnpmutant allele. Using single-cell transcriptomics, cut&run-seq, and histological approaches, we show that during neocortical development, Adnp orchestrates the production of late-born, upper-layer neurons through a two-step process. First, Adnp is required to sustain progenitor proliferation specifically during the developmental window for upper-layer cortical neurogenesis. Accordingly, we found that Adnp recruits the ChAHP subunit Chd4 to genes associated with progenitor proliferation. Second, in postmitotic differentiated neurons, we define a network of risk genes linked to NDDs that are regulated by Adnp and Chd4. Taken together, these data demonstrate that ChAHP is critical for driving the expansion of upper-layer cortical neurons and for regulating neuronal gene expression programs, suggesting that these processes may potentially contribute to NDD etiology.

Science & Technology - Other Topics↗

Divergence of Velocity Fields in Electrochemical Systems

The passage of current through electrochemical systems results in the development of concentration gradients in the electrolytic phase that can be modeled using concentrated solution theory. Application of this theory requires knowledge of three concentration-dependent transport coefficients, which are often taken to be conductivity, diffusion coefficient, and the cation transference number with respect to the solvent velocity. The governing diffusion equation for molar concentration contains two additional terms - the thermodynamic factor which is related to activity of the electrolytic species and the solvent velocity. The main advance in this paper is the derivation of an expression for the divergence of the solvent velocity. Solving this equation requires knowledge of the partial molar volume of the electrolyte. Analogous expressions are derived for the mass average and molar average velocities. These velocities occur naturally in the diffusion equation if concentration is expressed as weight fraction and mole fraction of the electrolytic species, respectively.

diffusion↗

Sex- and age-specific aspects of human peripheral T-cell dynamics

Background: The diversity of the antigenic T cell receptor (TCR) repertoire clonally expressed on T lymphocytes is a key element of the adaptive immune system protective functions. A decline in diversity in the older adults is associated with health deterioration. This diversity is generated by the rearrangement of TRB genes coding for TCR chains during lymphocyte differentiation in the thymus, but is essentially maintained by peripheral T lymphocytes proliferation for most of life. Deep sequencing of rearranged TRB genes from blood cells allows the monitoring of peripheral T cell repertoire dynamics. We analysed two aspects of rearranged TRB diversity, related to T lymphocyte proliferation and to the distribution of the T cell clone size, in a collection of repertoires obtained from 1 to 74 years-old donors. Results: Our results show that peripheral T lymphocytes expansion differs according to the recombination status of their TRB loci. Their proliferation rate changes with age, with different patterns in men and women. T cell clone size becomes more heterogeneous with time, and, in adults, is always more even in women. Importantly, a longitudinal analysis of TRB repertoires obtained at ten years intervals from individual men and women confirms the findings of this cross-sectional study. Conclusions: Peripheral T lymphocyte proliferation partially depends on their thymic developmental history. The rate of proliferation of T cells differing in their TRB rearrangement status is different in men and women before the age of 18 years old, but similar thereafter.

59 BASIC BIOLOGICAL SCIENCES↗

Mechanism of tethered agonist-mediated signaling by polycystin-1

Mutations of polycystin-1 (PC1) are the major cause (85% of cases) of autosomal dominant polycystic kidney disease (ADPKD), which is the fourth leading cause of kidney failure. PC1 is thought to function as an atypical G protein-coupled receptor, yet the mechanism by which PC1 regulates G-protein signaling remains poorly understood. A significant portion of ADPKD mutations of PC1 encode a protein with defects in maturation or reduced function that may be amenable to functional rescue. In this work, we have combined complementary biochemical and cellular assay experiments and accelerated molecular simulations, which revealed an allosteric transduction pathway in activation of the PC1 C-terminal fragment. Our findings will facilitate future rational drug design efforts targeting the PC1 signaling function. Polycystin-1 (PC1) is an important unusual G protein-coupled receptor (GPCR) with 11 transmembrane domains, and its mutations account for 85% of cases of autosomal dominant polycystic kidney disease (ADPKD). PC1 shares multiple characteristics with Adhesion GPCRs. These include a GPCR proteolysis site that autocatalytically divides these proteins into extracellular, N-terminal, and membrane-embedded, C-terminal fragments (CTF), and a tethered agonist (TA) within the N-terminal stalk of the CTF that is suggested to activate signaling. However, the mechanism by which a TA can activate PC1 is not known. Here, we have combined functional cellular signaling experiments of PC1 CTF expression constructs encoding wild type, stalkless, and three different ADPKD stalk variants with all-atom Gaussian accelerated molecular dynamics (GaMD) simulations to investigate TA-mediated signaling activation. Correlations of residue motions and free-energy profiles calculated from the GaMD simulations correlated with the differential signaling abilities of wild type and stalk variants of PC1 CTF. They suggested an allosteric mechanism involving residue interactions connecting the stalk, Tetragonal Opening for Polycystins (TOP) domain, and putative pore loop in TA-mediated activation of PC1 CTF. Key interacting residues such as N3074–S3585 and R3848–E4078 predicted from the GaMD simulations were validated by mutagenesis experiments. Together, these complementary analyses have provided insights into a TA-mediated activation mechanism of PC1 CTF signaling, which will be important for future rational drug design targeting PC1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conformation-specific synthetic intrabodies modulate mTOR signaling with subcellular spatial resolution

Subcellular compartmentalization is integral to the spatial regulation of mechanistic target of rapamycin (mTOR) signaling. However, the biological outputs associated with location-specific mTOR signaling events are poorly understood and challenging to decouple. Here, we engineered synthetic intracellular antibodies (intrabodies) that are capable of modulating mTOR signaling with genetically programmable spatial resolution. Epitope-directed phage display was exploited to generate high affinity synthetic antibody fragments (Fabs) against the FKBP12–Rapamycin binding site of mTOR (mTOR FRB ). We determined high-resolution crystal structures of two unique Fabs that discriminate distinct conformational states of mTOR FRB through recognition of its substrate recruitment interface. By leveraging these conformation-specific binders as intracellular probes, we uncovered the structural basis for an allosteric mechanism governing mTOR complex 1 (mTORC1) stability mediated by subtle structural adjustments within mTOR FRB . Furthermore, our results demonstrated that synthetic binders emulate natural substrates by employing divergent yet complementary hydrophobic residues at defined positions, underscoring the broad molecular recognition capability of mTOR FRB . Intracellular signaling studies showed differential time-dependent inhibition of S6 kinase 1 and Akt phosphorylation by genetically encoded intrabodies, thus supporting a mechanism of inhibition analogous to the natural product rapamycin. Finally, we implemented a feasible approach to selectively modulate mTOR signaling in the nucleus through spatially programmed intrabody expression. These findings establish intrabodies as versatile tools for dissecting the conformational regulation of mTORC1 and should be useful to explore how location-specific mTOR signaling influences disease progression.

Science & Technology - Other Topics↗

A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds

Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

97 MATHEMATICS AND COMPUTING↗

Phosphorylation barcodes direct biased chemokine signaling at CXCR3

Chemokine receptors, a group of G protein-coupled receptors (GPCRs), interact with transducers such as G proteins, ß-arrestins, and GPCR kinases (GRKs). In the chemokine system, many chemokine agonists act as “biased agonists” that preferentially activate distinct signaling effectors when binding to the same receptor, resulting in distinct physiological effects. Although one third of FDA-approved drugs target GPCRs, there has been limited success in targeting the chemokine system. Currently, there is little evidence that differential receptor phosphorylation, or “phosphorylation barcodes,” direct these biased responses at chemokine receptors. To address this knowledge gap, we used mass spectrometry to demonstrate that chemokines of CXCR3 promote different ensembles of phosphorylation barcodes that are associated with differential activation of G proteins, ß-arrestins and GRKs. Chemokine stimulation also resulted in distinct changes throughout the kinome in global phosphoproteomic studies. Mutation of specific CXCR3 phosphosites altered ß-arrestin conformation and impacted ß-arrestin activation in molecular dynamics simulations. T-cells expressing phosphorylation-deficient CXCR3 mutants resulted in distinct agonist- and receptor-specific chemotactic and signaling profiles that were not completely explained by engagement of G proteins, ß-arrestins, and GRKs alone. In conclusion, our results directly link distinct GPCR phosphorylation patterns with non-redundant chemokine signaling (Figure 1).

CXCR3↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

Acarbose impairs gut Bacteroides growth by targeting intracellular glucosidases

ABSTRACT Acarbose is a type 2 diabetes medicine that prevents dietary starch breakdown into glucose by inhibiting host amylase and glucosidase enzymes. Numerous gut species in theBacteroidesgenus enzymatically break down starch and change in relative abundance within the gut microbiome in acarbose-treated individuals. To mechanistically explain this observation, we used two model starch-degradingBacteroides,Bacteroides ovatus(Bo), andBacteroides thetaiotaomicron(Bt). Bt growth on starch polysaccharides is severely impaired by acarbose, whereas Bo growth is much less affected by the drug. TheBacteroidesuse a starch utilization system (Sus) to grow on starch. We hypothesized that Bo and Bt Sus enzymes are differentially inhibited by acarbose. Instead, we discovered that although acarbose primarily targets the Sus periplasmic GH97 enzymes in both organisms, the drug affects starch processing at multiple other points. Acarbose competes for transport through the TonB-dependent SusC proteins and binds to the Sus transcriptional regulators. Furthermore, Bo expresses a non-Sus GH97 (BoGH97D) when grown in starch with acarbose. The Bt homolog, BtGH97H, is not expressed in the same conditions, nor can overexpression of BoGH97D complement the Bt growth inhibition in the presence of acarbose. This work informs us about unexpected complexities of Sus function and regulation inBacteroides, including variation between related species. Furthermore, this indicates that the gut microbiome may be a source of variable response to acarbose treatment for diabetes. IMPORTANCE Acarbose is a type 2 diabetes medication that works primarily by stopping starch breakdown into glucose in the small intestine. This is accomplished by the inhibition of host enzymes, leading to better blood sugar control via reduced ability to derive glucose from dietary starches. The drug and undigested starch travel to the large intestine where acarbose interferes with the ability of some bacteria to grow on starch. However, little is known about how gut bacteria interact with acarbose, including microbes that can use starch as a carbon source. Here, we show that two gut species,Bacteroides ovatus(Bo) andBacteroides thetaiotaomicron(Bt), respond differently to acarbose: Bt growth is inhibited by acarbose, while Bo growth is less affected. We reveal a complex set of mechanisms involving differences in starch import and sensing behind the different Bo and Bt responses. This indicates the gut microbiome may be a source of variable response to acarbose treatment for diabetes via complex mechanisms in common gut microbes.

Microbiology↗

Salmonella Typhoid Toxin PltB Subunit and Its Non-typhoidal Salmonella Ortholog Confer Differential Host Adaptation and Virulence

Typhoidal and non-typhoidal Salmonelleae (NTS) cause typhoid fever and gastroenteritis, respectively, in humans. Salmonella typhoid toxin contributes to typhoid disease progression and chronic infection, but little is known about the role of its NTS ortholog. We found that typhoid toxin and its NTS ortholog induce different clinical presentations. The PltB subunit of each toxin exhibits different glycan-binding preferences that correlate with glycan expression profiles of host cells targeted by each bacterium at the primary infection or intoxication sites. Through co-crystal structures of PltB subunits bound to specific glycan receptor moieties, we show that they induce markedly different glycan-binding preferences and virulence outcomes. Furthermore, immunization with the NTS S. Javiana or its toxin offers cross-reactive protection against lethal-dose typhoid toxin challenge. Cumulatively, these results offer insights into the evolution of host adaptations in Salmonella AB toxins, their cell and tissue tropisms, and the design for improved typhoid vaccines and therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Transporters and ascorbate–glutathione metabolism for differential cadmium accumulation and tolerance in two contrasting willow genotypes

Salix matsudana Koidz is a low cadmium (Cd)-accumulating willow, whereas its cultivated variety, Salix matsudana var. matsudana f. umbraculifera Rehd., is a high Cd-accumulating and tolerant willow (HCW). The physiological and molecular mechanisms underlying differential Cd accumulation and tolerance in the two Salix species are poorly understood. In this study, we confirmed that the differential Cd translocation capacity from roots to the shoots leads to the differential Cd accumulation in their aboveground parts between these two willow genotypes. Cadmium accumulation happens preferentially in the transport pathway, and Cd is mainly located in the vacuolar, cell wall and intercellular space in HCW bark by cadmium location analysis at tissue and subcellular levels. Comparative transcriptome analysis revealed that higher expressions of several metal transporter genes (ATP-binding cassette transporters, K + transporters/channels, yellow stripe-like proteins, zinc-regulated transporter/iron-regulated transporter-like proteins, etc.) are involved in root uptake and translocation capacity in HCW; meanwhile, ascorbate–glutathione metabolic pathways play essential roles in Cd detoxification and higher tolerance of the Cd-accumulator HCW. These results lay the foundation for further understanding the molecular mechanisms of Cd accumulation in woody plants and provide new insights into molecular-assisted-screening woody plant varieties for phytoremediation.

59 BASIC BIOLOGICAL SCIENCES↗

Transverse momentum dependent PDFs at N3LO

We compute the quark and gluon transverse momentum dependent parton distribution functions at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. Our calculation is based on an expansion of the differential Drell-Yan and gluon fusion Higgs production cross sections about their collinear limit. This method allows us to employ cutting edge multiloop techniques for the computation of cross sections to extract these universal building blocks of the collinear limit of QCD. The corresponding perturbative matching kernels for all channels are expressed in terms of simple harmonic polylogarithms up to weight five. As a byproduct, we confirm a previous computation of the soft function for transverse momentum factorization at N 3 LO. Our results are the last missing ingredient to extend the q T subtraction methods to N 3 LO and to obtain resummed q T spectra at N 3 LL' accuracy both for gluon as well as for quark initiated processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Learning in Modal Space: Solving Time-Dependent Stochastic PDEs Using Physics-Informed Neural Networks

One of the open problems in scientific computing is the long-time integration of nonlinear stochastic partial differential equations (SPDEs), especially with arbitrary initial data. We address this problem by taking advantage of recent advances in scientific machine learning and the spectral dynamically orthogonal (DO) and borthogonal (BO) methods for representing stochastic processes. The recently introduced DO/BO methods reduce the SPDE to solving a system of deterministic PDEs and a system of stochastic ordinary differential equations. Specifically, we propose two new physics-informed neural networks (PINNs) for solving time-dependent SPDEs, namely the neural network (NN)-DO/BO methods. The proposed methods incorporate the DO/BO constraints into the loss function (along with the modal decomposition of the SPDE) with an implicit form instead of generating explicit expressions for the temporal derivatives of the DO/BO modes. Hence, the NN-DO/BO methods can overcome some of the drawbacks of the original DO/BO methods. For example, we do not need the assumption that the covariance matrix of the random coefficients is invertible as in the original DO method, and we can remove the assumption of no eigenvalue crossing as in the original BO method. Moreover, the NN-DO/BO methods can be used to solve time-dependent stochastic inverse problems with the same formulation and same computational complexity as for forward problems. Furthermore, we demonstrate the capability of the proposed methods via several numerical examples, namely: (1) A linear stochastic advection equation with deterministic initial condition: we obtain good results with the proposed methods, while the original DO/BO methods cannot be applied directly in this case. (2) Long-time integration of the stochastic Burgers' equation: we show the good performance of NN-DO/BO methods, especially the effectiveness of the NN-BO approach for such problems with many eigenvalue crossings during the whole time evolution, while the original BO method fails. (3) Nonlinear reaction diffusion equation: we consider both the forward problem and the inverse problems, including very noisy initial point values, to investigate the flexibility of the NN-DO/BO methods in handling inverse and mixed type problems. Taken together, these simulation results demonstrate that the NN-DO/BO methods can be employed to effectively quantify uncertainty propagation in a wide range of physical problems, but future work should address the efficiency issue of PINNs for forward problems.

97 MATHEMATICS AND COMPUTING↗

Intestinal single-cell atlas reveals novel lymphocytes in pigs with similarities to human cells

Lymphocytes can heavily influence intestinal health, but resolving intestinal lymphocyte function is challenging as the intestine contains a vastly heterogeneous mixture of cells. Pigs are an advantageous biomedical model, but deeper understanding of intestinal lymphocytes is warranted to improve model utility. Twenty-six cell types were identified in the porcine ileum by single-cell RNA sequencing and further compared with cells in human and murine ileum. Though general consensus of cell subsets across species was revealed, some porcine-specific lymphocyte subsets were identified. Differential tissue dissection and in situ analyses conferred spatial context, revealing similar locations of lymphocyte subsets in Peyer’s patches and epithelium in pig-to-human comparisons. Like humans, activated and effector lymphocytes were abundant in the ileum but not periphery of pigs, suggesting tissue-specific and/or activation-associated gene expression. Gene signatures for peripheral and ileal innate lymphoid cells newly discovered in pigs were defined and highlighted similarities to human innate lymphoid cells. Overall, we reveal novel lymphocyte subsets in pigs and highlight utility of pigs for intestinal research applications.

Life Sciences & Biomedicine - Other Topics↗

Single-cell profiling of Arabidopsis leaves to Pseudomonas syringae infection

Plant response to pathogen infection varies within a leaf, yet this heterogeneity is not well resolved. We expose Arabidopsis to Pseudomonas syringae or mock treatment and profile >11,000 individual cells using single-cell RNA sequencing. Integrative analysis of cell populations from both treatments identifies distinct pathogen-responsive cell clusters exhibiting transcriptional responses ranging from immunity to susceptibility. Pseudotime analyses through pathogen infection reveals a continuum of disease progression from an immune to a susceptible state. Confocal imaging of promoter-reporter lines for transcripts enriched in immune cell clusters shows expression surrounding substomatal cavities colonized or in close proximity to bacterial colonies, suggesting that cells within immune clusters represent sites of early pathogen invasion. Susceptibility clusters exhibit more general localization and are highly induced at later stages of infection. Overall, our work shows cellular heterogeneity within an infected leaf and provides insight into plant differential response to infection at a single-cell level.

59 BASIC BIOLOGICAL SCIENCES↗