Single-cell transcriptomics sheds light on the identity and metabolism of developing leaf cells.

As the main photosynthetic instruments of vascular plants, leaves are crucial and complex plant organs. A strict organization of leaf mesophyll and epidermal cell layers orchestrates photosynthesis and gas exchange. In addition, water and nutrients for leaf growth are transported through the vascular tissue. To establish the single-cell transcriptomic landscape of these different leaf tissues, we performed high-throughput transcriptome sequencing of individual cells isolated from young leaves of Arabidopsis (Arabidopsis thaliana) seedlings grown in two different environmental conditions. The detection of approximately 19,000 different transcripts in over 1,800 high-quality leaf cells revealed 14 cell populations composing the young, differentiating leaf. Besides the cell populations comprising the core leaf tissues, we identified subpopulations with a distinct identity or metabolic activity. In addition, we proposed cell-type-specific markers for each of these populations. Finally, an intuitive web tool allows for browsing the presented dataset. Our data present insights on how the different cell populations constituting a developing leaf are connected via developmental, metabolic, or stress-related trajectories.

A redundant transcription factor network steers spatiotemporal Arabidopsis triterpene synthesis.

Plant specialized metabolites modulate developmental and ecological functions and comprise many therapeutic and other high-value compounds. However, the mechanisms determining their cell-specific expression remain unknown. Here we describe the transcriptional regulatory network that underlies cell-specific biosynthesis of triterpenes in Arabidopsis thaliana root tips. Expression of thalianol and marneral biosynthesis pathway genes depends on the phytohormone jasmonate and is limited to outer tissues. We show that this is promoted by the activity of redundant bHLH-type transcription factors from two distinct clades and coactivated by homeodomain factors. Conversely, the DOF-type transcription factor DAG1 and other regulators prevent expression of the triterpene pathway genes in inner tissues. We thus show how precise expression of triterpene biosynthesis genes is determined by a robust network of transactivators, coactivators and counteracting repressors.

Comparative analysis of antibody- and lipid-based multiplexing methods for single-cell RNA-seq.

BACKGROUND: Multiplexing of samples in single-cell RNA-seq studies allows a significant reduction of the experimental costs, straightforward identification of doublets, increased cell throughput, and reduction of sample-specific batch effects. Recently published multiplexing techniques using oligo-conjugated antibodies or -lipids allow barcoding sample-specific cells, a process called "hashing." RESULTS: Here, we compare the hashing performance of TotalSeq-A and -C antibodies, custom synthesized lipids and MULTI-seq lipid hashes in four cell lines, both for single-cell RNA-seq and single-nucleus RNA-seq. We also compare TotalSeq-B antibodies with CellPlex reagents (10x Genomics) on human PBMCs and TotalSeq-B with different lipids on primary mouse tissues. Hashing efficiency was evaluated using the intrinsic genetic variation of the cell lines and mouse strains. Antibody hashing was further evaluated on clinical samples using PBMCs from healthy and SARS-CoV-2 infected patients, where we demonstrate a more affordable approach for large single-cell sequencing clinical studies, while simultaneously reducing batch effects. CONCLUSIONS: Benchmarking of different hashing strategies and computational pipelines indicates that correct demultiplexing can be achieved with both lipid- and antibody-hashed human cells and nuclei, with MULTISeqDemux as the preferred demultiplexing function and antibody-based hashing as the most efficient protocol on cells. On nuclei datasets, lipid hashing delivers the best results. Lipid hashing also outperforms antibodies on cells isolated from mouse brain. However, antibodies demonstrate better results on tissues like spleen or lung.

Distinct Transcriptional Programs in Ascitic and Solid Cancer Cells Induce Different Responses to Chemotherapy in High-Grade Serous Ovarian Cancer.

High-grade serous ovarian cancer (HGSOC) is responsible for the largest number of ovarian cancer deaths. The frequent therapy-resistant relapses necessitate a better understanding of mechanisms driving therapy resistance. Therefore, we mapped more than a hundred thousand cells of HGSOC patients in different phases of the disease, using single-cell RNA sequencing. Within patients, we compared chemonaive with chemotreated samples. As such, we were able to create a single-cell atlas of different HGSOC lesions and their treatment. This revealed a high intrapatient concordance between spatially distinct metastases. In addition, we found remarkable baseline differences in transcriptomics of ascitic and solid cancer cells, resulting in a different response to chemotherapy. Moreover, we discovered different robust subtypes of cancer-associated fibroblasts (CAF) in all patients. Besides inflammatory CAFs, vascular CAFs, and matrix CAFs, we identified a new CAF subtype that was characterized by high expression of STAR, TSPAN8, and ALDH1A1 and clearly enriched after chemotherapy. Together, tumor heterogeneity in both cancer and stromal cells contributes to therapy resistance in HGSOC and could form the basis of novel therapeutic strategies that differentiate between ascitic and solid disease. IMPLICATIONS: The newly characterized differences between ascitic and solid cancer cells before and after chemotherapy could inform novel treatment strategies for metastatic HGSOC.

TIM3+ TRBV11-2 T cells and IFNγ signature in patrolling monocytes and CD16+ NK cells delineate MIS-C.

In rare instances, pediatric SARS-CoV-2 infection results in a novel immunodysregulation syndrome termed multisystem inflammatory syndrome in children (MIS-C). We compared MIS-C immunopathology with severe COVID-19 in adults. MIS-C does not result in pneumocyte damage but is associated with vascular endotheliitis and gastrointestinal epithelial injury. In MIS-C, the cytokine release syndrome is characterized by IFNγ and not type I interferon. Persistence of patrolling monocytes differentiates MIS-C from severe COVID-19, which is dominated by HLA-DRlo classical monocytes. IFNγ levels correlate with granzyme B production in CD16+ NK cells and TIM3 expression on CD38+/HLA-DR+ T cells. Single-cell TCR profiling reveals a skewed TCRß repertoire enriched for TRBV11-2 and a superantigenic signature in TIM3+/CD38+/HLA-DR+ T cells. Using NicheNet, we confirm IFNγ as a central cytokine in the communication between TIM3+/CD38+/HLA-DR+ T cells, CD16+ NK cells, and patrolling monocytes. Normalization of IFNγ, loss of TIM3, quiescence of CD16+ NK cells, and contraction of patrolling monocytes upon clinical resolution highlight their potential role in MIS-C immunopathogenesis.

Maize specialized metabolome networks reveal organ-preferential mixed glycosides.

Despite the scientific and economic importance of maize, little is known about its specialized metabolism. Here, five maize organs were profiled using different reversed-phase liquid chromatography-mass spectrometry methods. The resulting spectral metadata, combined with candidate substrate-product pair (CSPP) networks, allowed the structural characterization of 427 of the 5,420 profiled compounds, including phenylpropanoids, flavonoids, benzoxazinoids, and auxin-related compounds, among others. Only 75 of the 427 compounds were already described in maize. Analysis of the CSPP networks showed that phenylpropanoids are present in all organs, whereas other metabolic classes are rather organ-enriched. Frequently occurring CSPP mass differences often corresponded with glycosyl- and acyltransferase reactions. The interplay of glycosylations and acylations yields a wide variety of mixed glycosides, bearing substructures corresponding to the different biochemical classes. For example, in the tassel, many phenylpropanoid and flavonoid-bearing glycosides also contain auxin-derived moieties. The characterized compounds and mass differences are an important step forward in metabolic pathway discovery and systems biology research. The spectral metadata of the 5,420 compounds is publicly available (DynLib spectral database, https://bioit3.irc.ugent.be/dynlib/).

OTULIN maintains skin homeostasis by controlling keratinocyte death and stem cell identity.

OTULIN is a deubiquitinase that specifically cleaves linear ubiquitin chains. Here we demonstrate that the ablation of Otulin selectively in keratinocytes causes inflammatory skin lesions that develop into verrucous carcinomas. Genetic deletion of Tnfr1, knockin expression of kinase-inactive Ripk1 or keratinocyte-specific deletion of Fadd and Mlkl completely rescues mice with OTULIN deficiency from dermatitis and tumorigenesis, thereby identifying keratinocyte cell death as the driving force for inflammation. Single-cell RNA-sequencing comparing non-lesional and lesional skin reveals changes in epidermal stem cell identity in OTULIN-deficient keratinocytes prior to substantial immune cell infiltration. Keratinocytes lacking OTULIN display a type-1 interferon and IL-1ß response signature, and genetic or pharmacologic inhibition of these cytokines partially inhibits skin inflammation. Finally, expression of a hypomorphic mutant Otulin allele, previously shown to cause OTULIN-related autoinflammatory syndrome in humans, induces a similar inflammatory phenotype, thus supporting the importance of OTULIN for restraining skin inflammation and maintaining immune homeostasis.

Non-cell autonomous and spatiotemporal signalling from a tissue organizer orchestrates root vascular development.

During plant development, a precise balance of cytokinin is crucial for correct growth and patterning, but it remains unclear how this is achieved across different cell types and in the context of a growing organ. Here we show that in the root apical meristem, the TMO5/LHW complex increases active cytokinin levels via two cooperatively acting enzymes. By profiling the transcriptomic changes of increased cytokinin at single-cell level, we further show that this effect is counteracted by a tissue-specific increase in CYTOKININ OXIDASE 3 expression via direct activation of the mobile transcription factor SHORTROOT. In summary, we show that within the root meristem, xylem cells act as a local organizer of vascular development by non-autonomously regulating cytokinin levels in neighbouring procambium cells via sequential induction and repression modules.

Cell-by-cell dissection of phloem development links a maturation gradient to cell specialization.

In the plant meristem, tissue-wide maturation gradients are coordinated with specialized cell networks to establish various developmental phases required for indeterminate growth. Here, we used single-cell transcriptomics to reconstruct the protophloem developmental trajectory from the birth of cell progenitors to terminal differentiation in the Arabidopsis thaliana root. PHLOEM EARLY DNA-BINDING-WITH-ONE-FINGER (PEAR) transcription factors mediate lineage bifurcation by activating guanosine triphosphatase signaling and prime a transcriptional differentiation program. This program is initially repressed by a meristem-wide gradient of PLETHORA transcription factors. Only the dissipation of PLETHORA gradient permits activation of the differentiation program that involves mutual inhibition of early versus late meristem regulators. Thus, for phloem development, broad maturation gradients interface with cell-type-specific transcriptional regulators to stage cellular differentiation.

Vascular transcription factors guide plant epidermal responses to limiting phosphate conditions.

Optimal plant growth is hampered by deficiency of the essential macronutrient phosphate in most soils. Plant roots can, however, increase their root hair density to efficiently forage the soil for this immobile nutrient. By generating and exploiting a high-resolution single-cell gene expression atlas of Arabidopsis roots, we show an enrichment of TARGET OF MONOPTEROS 5/LONESOME HIGHWAY (TMO5/LHW) target gene responses in root hair cells. The TMO5/LHW heterodimer triggers biosynthesis of mobile cytokinin in vascular cells and increases root hair density during low-phosphate conditions by modifying both the length and cell fate of epidermal cells. Moreover, root hair responses in phosphate-deprived conditions are TMO5- and cytokinin-dependent. Cytokinin signaling links root hair responses in the epidermis to perception of phosphate depletion in vascular cells.