Measuring cell identity in noisy biological systems

Cortical arealization of interneurons defines shared and distinct molecular programs in developing human and macaque brains

Cortical interneurons generated from ganglionic eminence via a long-distance journey of tangential migration display evident cellular and molecular differences across brain regions, which seeds the heterogeneous cortical circuitry in primates. However, whether such regional specifications in interneurons are intrinsically encoded or gained through interactions with the local milieu remains elusive. Here, we recruit 685,692 interneurons from cerebral cortex and subcortex including ganglionic eminence within the developing human and macaque species. Our integrative and comparative analyses reveal that less transcriptomic alteration is accompanied by interneuron migration within the ganglionic eminence subdivisions, in contrast to the dramatic changes observed in cortical tangential migration, which mostly characterize the transcriptomic specification for different destinations and for species divergence. Moreover, the in-depth survey of temporal regulation illustrates species differences in the developmental dynamics of cell types, e.g., the employment of CRH in primate interneurons during late-fetal stage distinguishes from their postnatal emergence in mice, and our entropy quantifications manifest the interneuron diversities gradually increase along the developmental ages in human and macaque cerebral cortices. Overall, our analyses depict the spatiotemporal features appended to cortical interneurons, providing a new proxy for understanding the relationship between cellular diversity and functional progression.

Dynamics of alternative polyadenylation in single root cells of Arabidopsis thaliana

Introduction

Single-cell RNA-seq (scRNA-seq) technologies have been widely used to reveal the diversity and complexity of cells, and pioneering studies on scRNA-seq in plants began to emerge since 2019. However, existing studies on plants utilized scRNA-seq focused only on the gene expression regulation. As an essential post-transcriptional mechanism for regulating gene expression, alternative polyadenylation (APA) generates diverse mRNA isoforms with distinct 3' ends through the selective use of different polyadenylation sites in a gene. APA plays important roles in regulating multiple developmental processes in plants, such as flowering time and stress response.

Methods

In this study, we developed a pipeline to identify and integrate APA sites from different scRNA-seq data and analyze APA dynamics in single cells. First, high-confidence poly(A) sites in single root cells were identified and quantified. Second, three kinds of APA markers were identified for exploring APA dynamics in single cells, including differentially expressed poly(A) sites based on APA site expression, APA markers based on APA usages, and APA switching genes based on 3' UTR (untranslated region) length change. Moreover, cell type annotations of single root cells were refined by integrating both the APA information and the gene expression profile.

Results

We comprehensively compiled a single-cell APA atlas from five scRNA-seq studies, covering over 150,000 cells spanning four major tissue branches, twelve cell types, and three developmental stages. Moreover, we quantified the dynamic APA usages in single cells and identified APA markers across tissues and cell types. Further, we integrated complementary information of gene expression and APA profiles to annotate cell types and reveal subtle differences between cell types.

Discussion

This study reveals that APA provides an additional layer of information for determining cell identity and provides a landscape of APA dynamics during Arabidopsis root development.

A common regulatory switch controls a suite of C4 traits in multiple cell types

The C4 photosynthetic pathway provided a major advantage to plants growing in hot, dry environments, including the ancestors of our most productive crops. Two traits were essential for the evolution of this pathway: increased vein density and the functionalization of bundle sheath cells for photosynthesis. Although GRAS transcriptional regulators, including SHORT ROOT (SHR), have been implicated in mediating leaf patterning in both C3 and C4 species, little is known about what controls the specialized features of the cells that mediate C4 metabolism and physiology. We show in the model monocot, Setaria viridis, that SHR regulates components of multiple cell identities, including chloroplast biogenesis and photosynthetic gene expression in bundle sheath cells, a central feature of C4 plants. Furthermore, we found that it also contributes to the two-cell compartmentalization of the characteristic four-carbon shuttle pathway. Disruption of SHR function clearly reduced photosynthetic capacity and seed yield in mutant plants under heat stress. Together, these results show how cell identities are remodeled by SHR to host the suite of traits characteristic of C4 regulation, which are a main engineering target in non-C4 crops to improve climate resilience.

Spatiotemporal expression patterns of anxiety disorder-associated genes

Anxiety disorders (ADs) are the most common form of mental disorder that affects millions of individuals worldwide. Although physiological studies have revealed the neural circuits related to AD symptoms, how AD-associated genes are spatiotemporally expressed in the human brain still remains unclear. In this study, we integrated genome-wide association studies of four human AD subtypes-generalized anxiety disorder, social anxiety disorder, panic disorder, and obsessive-compulsive disorder-with spatial gene expression patterns. Our investigation uncovered a novel division among AD-associated genes, marked by significant and distinct expression enrichments in the cerebral nuclei, limbic, and midbrain regions. Each gene cluster was associated with specific anxiety-related behaviors, signaling pathways, region-specific gene networks, and cell types. Notably, we observed a significant negative correlation in the temporal expression patterns of these gene clusters during various developmental stages. Moreover, the specific brain regions enriched in each gene group aligned with neural circuits previously associated with negative decision-making and anxious temperament. These results suggest that the two distinct gene clusters may underlie separate neural systems involved in anxiety. As a result, our findings bridge the gap between genes and neural circuitry, shedding light on the mechanisms underlying AD-associated behaviors.

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.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

Crafting a blueprint for single-cell RNA sequencing

Droplet-based single-cell RNA sequencing (scRNA-Seq) has moved rapidly from a technology with great potential to a method applied to ever-broader questions. The detailed information that scRNA-Seq offers has proven incredibly powerful in resolving cell responses to developmental and environmental cues. However, to maximize the potential of this technology, a panoply of upstream, practical points require consideration. Principal among these are the optimization of cell-isolation procedures, accommodating biotic/abiotic stress responses, and discerning the number of cells and sequencing reads needed. To complement excellent reviews outlining applications and data analysis tools for scRNA-Seq, we here discuss these considerations and provide practical tips to tailor experimental design and ensure the best possible outcome.

Plant Metabolic Network 15: A resource of genome-wide metabolism databases for 126 plants and algae

To understand and engineer plant metabolism, we need a comprehensive and accurate annotation of all metabolic information across plant species. As a step towards this goal, we generated genome-scale metabolic pathway databases of 126 algal and plant genomes, ranging from model organisms to crops to medicinal plants (https://plantcyc.org). Of these, 104 have not been reported before. We systematically evaluated the quality of the databases, which revealed that our semi-automated validation pipeline dramatically improves the quality. We then compared the metabolic content across the 126 organisms using multiple correspondence analysis and found that Brassicaceae, Poaceae, and Chlorophyta appeared as metabolically distinct groups. To demonstrate the utility of this resource, we used recently published sorghum transcriptomics data to discover previously unreported trends of metabolism underlying drought tolerance. We also used single-cell transcriptomics data from the Arabidopsis root to infer cell type-specific metabolic pathways. This work shows the quality and quantity of our resource and demonstrates its wide-ranging utility in integrating metabolism with other areas of plant biology.

Molecular Mechanisms Driving Switch Behavior in Xylem Cell Differentiation

Plant xylem cells conduct water and mineral nutrients. Although most plant cells are totipotent, xylem cells are unusual and undergo terminal differentiation. Many genes regulating this process are well characterized, including the Vascular-related NAC Domain 7 (VND7), MYB46, and MYB83 transcription factors, which are proposed to act in interconnected feedforward loops (FFLs). Less is known regarding the molecular mechanisms underlying the terminal transition to xylem cell differentiation. Here, we generate whole-root and single-cell data, which demonstrate that VND7 initiates sharp switching of root cells to xylem cell identity. Based on these data, we identified 4 candidate VND7 downstream target genes capable of generating this switch. Although MYB46 responds to VND7 induction, it is not among these targets. This system provides an important model to study the emergent properties that may give rise to totipotency relative to terminal differentiation and reveals xylem cell subtypes.

Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data

We introduce quanTIseq, a method to quantify the fractions of ten immune cell types from bulk RNA-sequencing data. quanTIseq was extensively validated in blood and tumor samples using simulated, flow cytometry, and immunohistochemistry data.quanTIseq analysis of 8000 tumor samples revealed that cytotoxic T cell infiltration is more strongly associated with the activation of the CXCR3/CXCL9 axis than with mutational load and that deconvolution-based cell scores have prognostic value in several solid cancers. Finally, we used quanTIseq to show how kinase inhibitors modulate the immune contexture and to reveal immune-cell types that underlie differential patients' responses to checkpoint blockers.Availability: quanTIseq is available at http://icbi.at/quantiseq .

A Conceptual Framework for Cell Identity Transitions in Plants

Multicellular organisms develop from a single cell that proliferates to form different cell types with specialized functions. Sixty years ago, Waddington suggested the 'epigenetic landscape' as a useful metaphor for the process. According to this view, cells move through a rugged identity space along genetically encoded trajectories, until arriving at one of the possible final fates. In plants in particular, these trajectories have strong spatial correlates, as cell identity is intimately linked to its relative position within the plant. During regeneration, however, positional signals are severely disrupted and differentiated cells are able to undergo rapid non-canonical identity changes. Moreover, while pluripotent properties have long been ascribed to plant cells, the introduction of induced pluripotent stem cells in animal studies suggests such plasticity may not be unique to plants. As a result, current concepts of differentiation as a gradual and hierarchical process are being reformulated across biological fields. Traditional studies of plant regeneration have placed strong emphasis on the emergence of patterns and tissue organization, and information regarding the events occurring at the level of individual cells is only now beginning to emerge. Here, I review the historical and current concepts of cell identity and identity transitions, and discuss how new views and tools may instruct the future understanding of differentiation and plant regeneration.

Changes in Gene Expression in Space and Time Orchestrate Environmentally Mediated Shaping of Root Architecture

Shaping of root architecture is a quintessential developmental response that involves the concerted action of many different cell types, is highly dynamic, and underpins root plasticity. To determine to what extent the environmental regulation of lateral root development is a product of cell-type preferential activities, we tracked transcriptomic responses to two different treatments that both change root development in Arabidopsis thaliana at an unprecedented level of temporal detail. We found that individual transcripts are expressed with a very high degree of temporal and spatial specificity, yet biological processes are commonly regulated, in a mechanism we term response nonredundancy. Using causative gene network inference to compare the genes regulated in different cell types and during responses to nitrogen and a biotic interaction, we found that common transcriptional modules often regulate the same gene families but control different individual members of these families, specific to response and cell type. This reinforces that the activity of a gene cannot be defined simply as molecular function; rather, it is a consequence of spatial location, expression timing, and environmental responsiveness.

Root Regeneration Triggers an Embryo-like Sequence Guided by Hormonal Interactions

Plant roots can regenerate after excision of their tip, including the stem cell niche. To determine which developmental program mediates such repair, we applied a combination of lineage tracing, single-cell RNA sequencing, and marker analysis to test different models of tissue reassembly. We show that multiple cell types can reconstitute stem cells, demonstrating the latent potential of untreated plant cells. The transcriptome of regenerating cells prior to stem cell activation resembles that of an embryonic root progenitor. Regeneration defects are more severe in embryonic than in adult root mutants. Furthermore, the signaling domains of the hormones auxin and cytokinin mirror their embryonic dynamics and manipulation of both hormones alters the position of new tissues and stem cell niche markers. Our findings suggest that plant root regeneration follows, on a larger scale, the developmental stages of embryonic patterning and is guided by spatial information provided by complementary hormone domains.

The potential of single-cell profiling in plants

Single-cell transcriptomics has been employed in a growing number of animal studies, but the technique has yet to be widely used in plants. Nonetheless, early studies indicate that single-cell RNA-seq protocols developed for animal cells produce informative datasets in plants. We argue that single-cell transcriptomics has the potential to provide a new perspective on plant problems, such as the nature of the stem cells or initials, the plasticity of plant cells, and the extent of localized cellular responses to environmental inputs. Single-cell experimental outputs require different analytical approaches compared with pooled cell profiles and new tools tailored to single-cell assays are being developed. Here, we highlight promising new single-cell profiling approaches, their limitations as applied to plants, and their potential to address fundamental questions in plant biology.

Transcriptome dynamics of the stomatal lineage: birth, amplification, and termination of a self-renewing population

Developmental transitions can be described in terms of morphology and the roles of individual genes, but also in terms of global transcriptional and epigenetic changes. Temporal dissections of transcriptome changes, however, are rare for intact, developing tissues. We used RNA sequencing and microarray platforms to quantify gene expression from labeled cells isolated by fluorescence-activated cell sorting to generate cell-type-specific transcriptomes during development of an adult stem-cell lineage in the Arabidopsis leaf. We show that regulatory modules in this early lineage link cell types that had previously been considered to be under separate control and provide evidence for recruitment of individual members of gene families for different developmental decisions. Because stomata are physiologically important and because stomatal lineage cells exhibit exemplary division, cell fate, and cell signaling behaviors, this dataset serves as a valuable resource for further investigations of fundamental developmental processes.

Quantification of cell identity from single-cell gene expression profiles

The definition of cell identity is a central problem in biology. While single-cell RNA-seq provides a wealth of information regarding cell states, better methods are needed to map their identity, especially during developmental transitions. Here, we use repositories of cell type-specific transcriptomes to quantify identities from single-cell RNA-seq profiles, accurately classifying cells from Arabidopsis root tips and human glioblastoma tumors. We apply our approach to single cells captured from regenerating roots following tip excision. Our technique exposes a previously uncharacterized transient collapse of identity distant from the injury site, demonstrating the biological relevance of a quantitative cell identity index.

Peripheral blood cells inform on the presence of breast cancer: a population-based case-control study

Tumor–host interactions extend beyond the local microenvironment and cancer development largely depends on the ability of malignant cells to hijack and exploit the normal physiological processes of the host. Here, we established that many genes within peripheral blood cells show differential expression when an untreated breast cancer (BC) is present, and harnessed this fact to construct a 50-gene signature that distinguish BC patients from population-based controls. Our results were derived from a series of large datasets within our unique population-based Norwegian Women and Cancer cohort that allowed us to investigate the influence of medications and tumor characteristics on our blood-based test, and were further tested in two external datasets. Our 50-gene signature contained cytostatic signals including the specific suppression of the immune response and medications influencing transcription involved in those processes were identified as confounders. Through analysis of the biological processes differentially expressed in blood, we were able to provide a rationale as to why the systemic response of the host may be a reliable marker of BC, characterized by the underexpression of both immune-specific pathways and “universal” cell programs driven by MYC (i.e., metabolism, growth and cell cycle). In conclusion, gene expression of peripheral blood cells is markedly perturbed by the specific presence of carcinoma in the breast and these changes simultaneously engage a number of systemic cytostatic signals emerging connections with immune escape of BC.

Lessons from miR-143/145: the importance of cell-type localization of miRNAs

miR-143 and miR-145 are co-expressed microRNAs (miRNAs) that have been extensively studied as potential tumor suppressors. These miRNAs are highly expressed in the colon and are consistently reported as being downregulated in colorectal and other cancers. Through regulation of multiple targets, they elicit potent effects on cancer cell growth and tumorigenesis. Importantly, a recent discovery demonstrates that miR-143 and miR-145 are not expressed in colonic epithelial cells; rather, these two miRNAs are highly expressed in mesenchymal cells such as fibroblasts and smooth muscle cells. The expression patterns of miR-143 and miR-145 and other miRNAs were initially determined from tissue level data without consideration that multiple different cell types, each with their own unique miRNA expression patterns, make up each tissue. Herein, we discuss the early reports on the identification of dysregulated miR-143 and miR-145 expression in colorectal cancer and how lack of consideration of cellular composition of normal tissue led to the misconception that these miRNAs are downregulated in cancer. We evaluate mechanistic data from miR-143/145 studies in context of their cell type-restricted expression pattern and the potential of these miRNAs to be considered tumor suppressors. Further, we examine other examples of miRNAs being investigated in inappropriate cell types modulating pathways in a non-biological fashion. Our review highlights the importance of determining the cellular expression pattern of each miRNA, so that downstream studies are conducted in the appropriate cell type.

Transcriptional output in a prospective design conditionally on follow-up and exposure: the multistage model of cancer

Transcriptomics as the analysis of mRNA and microRNA could be implemented in prospective studies both in peripheral blood and tissues. Its application in cancer epidemiology could provide a new understanding of the functional changes underlying the multistage model of carcinogenesis, as well as the relationship between these changes and exposure to carcinogens. Transcriptomics is not merely another -omics technology for risk assessment in traditional prospective studies. Instead, this novel approach has the potential to estimate the distribution of gene expression conditionally on different exposures, and to study the length of the different stages of carcinogenesis. If it proves to be a valid approach, transcriptomics could be an opportunity to make meaningful advances in our understanding of the carcinogenic process.