ICAT: a novel algorithm to robustly identify cell states following perturbations in single-cell transcriptomes.

MOTIVATION: The detection of distinct cellular identities is central to the analysis of single-cell RNA sequencing (scRNA-seq) experiments. However, in perturbation experiments, current methods typically fail to correctly match cell states between conditions or erroneously remove population substructure. Here, we present the novel, unsupervised algorithm Identify Cell states Across Treatments (ICAT) that employs self-supervised feature weighting and control-guided clustering to accurately resolve cell states across heterogeneous conditions. RESULTS: Using simulated and real datasets, we show ICAT is superior in identifying and resolving cell states compared with current integration workflows. While requiring no a priori knowledge of extant cell states or discriminatory marker genes, ICAT is robust to low signal strength, high perturbation severity, and disparate cell type proportions. We empirically validate ICAT in a developmental model and find that only ICAT identifies a perturbation-unique cellular response. Taken together, our results demonstrate that ICAT offers a significant improvement in defining cellular responses to perturbation in scRNA-seq data. AVAILABILITY AND IMPLEMENTATION: https://github.com/BradhamLab/icat.

Polychrome labeling reveals skeletal triradiate and elongation dynamics and abnormalities in patterning cue-perturbed embryos.

The larval skeleton of the sea urchin Lytechinus variegatus is an ideal model system for studying skeletal patterning; however, our understanding of the etiology of skeletal patterning in sea urchin larvae is limited due to the lack of approaches to live-image skeleton formation. Calcium-binding fluorochromes have been used to study the temporal dynamics of bone growth and healing. To date, only calcein green has been used in sea urchin larvae to fluorescently label the larval skeleton. Here, we optimize labeling protocols for two additional calcium-binding fluorochromes: xylenol orange and calcein blue- and demonstrate that these fluorochromes can be used individually or in nested pulse-chase experiments to understand the temporal dynamics of skeletogenesis and patterning. Using a pulse-chase approach, we show that the initiation of skeletogenesis begins around 15 â€‹h post fertilization. We also assess the timing of triradiate formation in embryos treated with a range of patterning perturbagens and demonstrate that triradiate formation is delayed and asynchronous in embryos ventralized via treatment with either nickel or chlorate. Finally, we measure the extent of fluorochrome incorporation in triple-labeled embryos to determine the elongation rate of numerous skeletal elements throughout early skeletal patterning and compare this to the rate of skeletal growth in embryos treated with axitinib to inhibit VEGFR. We find that skeletal elements elongate much more slowly in axitinib-treated embryos, and that axitinib treatment is sufficient to induce abnormal orientation of the triradiates.

Cell Cycle Dynamics during Stomatal Development: Window of MUTE Action and Ramification of Its Loss-of-Function on an Uncommitted Precursor.

Plants develop in the absence of cell migration. As such, cell division and differentiation need to be coordinated for functional tissue formation. Cellular valves on the plant epidermis, stomata, are generated through a stereotypical sequence of cell division and differentiation events. In Arabidopsis, three master regulatory transcription factors, SPEECHLESS (SPCH), MUTE and FAMA, sequentially drive initiation, proliferation and differentiation of stomata. Among them, MUTE switches the cell cycle mode from proliferative asymmetric division to terminal symmetric division and orchestrates the execution of the single symmetric division event. However, it remains unclear to what extent MUTE regulates the expression of cell cycle genes through the symmetric division and whether MUTE accumulation itself is gated by the cell cycle. Here, we show that MUTE directly upregulates the expression of cell cycle components throughout the terminal cell cycle phases of a stomatal precursor, not only core cell cycle engines but also check-point regulators. Time-lapse live imaging using the multicolor Plant Cell Cycle Indicator revealed that MUTE accumulates up to the early G2 phase, whereas its successor and direct target, FAMA, accumulate at late G2 through terminal mitosis. In the absence of MUTE, meristemoids fail to differentiate and their G1 phase elongates as they reiterate asymmetric divisions. Together, our work provides the framework of cell cycle and master regulatory transcription factors to coordinate a single symmetric cell division and suggests a mechanism for the eventual cell cycle arrest of an uncommitted stem-cell-like precursor at the G1 phase.

Cell biology of the leaf epidermis: Fate specification, morphogenesis, and coordination.

As the outermost layer of plants, the epidermis serves as a critical interface between plants and the environment. During leaf development, the differentiation of specialized epidermal cell types, including stomatal guard cells, pavement cells, and trichomes, occurs simultaneously, each providing unique and pivotal functions for plant growth and survival. Decades of molecular-genetic and physiological studies have unraveled key players and hormone signaling specifying epidermal differentiation. However, most studies focus on only one cell type at a time, and how these distinct cell types coordinate as a unit is far from well-comprehended. Here we provide a review on the current knowledge of regulatory mechanisms underpinning the fate specification, differentiation, morphogenesis, and positioning of these specialized cell types. Emphasis is given to their shared developmental origins, fate flexibility, as well as cell cycle and hormonal controls. Furthermore, we discuss computational modeling approaches to integrate how mechanical properties of individual epidermal cell types and entire tissue/organ properties mutually influence each other. We hope to illuminate the underlying mechanisms coordinating the cell differentiation that ultimately generate a functional leaf epidermis.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes.

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

Spatially mapping gene expression in sea urchin primary mesenchyme cells.

During sea urchin embryogenesis, primary mesenchyme cells (PMCs) follow a stereotypical migratory program, arrange into a primary pattern, then begin to secrete a bilaterally symmetric calcium carbonate skeleton. Recently identified genes are expressed in spatially-restricted domains within the PMC population (Sun & Ettensohn, 2014). To better understand the molecular mechanisms orchestrating PMC positioning, we are characterizing the expression profiles of PMC subset-specific genes. To deconvolve the spatiotemporal expression patterns within PMCs, we detect cell-specific mRNA expression with combined RNA fluorescence in situ hybridization and immunolabeling of PMCs. Subsequent confocal microscopy provides 3D position and expression information for individual PMCs. We extract PMC positions and relative gene expression levels, then model these results using open-source 3D modeling software. This versatile protocol can be extended to other models and systems.

Zygotic LvBMP5-8 is required for skeletal patterning and for left-right but not dorsal-ventral specification in the sea urchin embryo.

Skeletal patterning in the sea urchin embryo requires coordinated signaling between the pattern-dictating ectoderm and the skeletogenic primary mesenchyme cells (PMCs); recent studies have begun to uncover the molecular basis for this process. Using an unbiased RNA-Seq-based screen, we have previously identified the TGF-ß superfamily ligand, LvBMP5-8, as a skeletal patterning gene in Lytechinus variegatus embryos. This result is surprising, since both BMP5-8 and BMP2/4 ligands have been implicated in sea urchin dorsal-ventral (DV) and left-right (LR) axis specification. Here, we demonstrate that zygotic LvBMP5-8 is required for normal skeletal patterning on the left side, as well as for normal PMC positioning during gastrulation. Zygotic LvBMP5-8 is required for expression of the left-side marker soxE, suggesting that LvBMP5-8 is required for left-side specification. Interestingly, we also find that LvBMP5-8 knockdown suppresses serotonergic neurogenesis on the left side. While LvBMP5-8 overexpression is sufficient to dorsalize embryos, we find that zygotic LvBMP5-8 is not required for normal DV specification or development. In addition, ectopic LvBMP5-8 does not dorsalize LvBMP2/4 morphant embryos, indicating that, in the absence of BMP2/4, BMP5-8 is insufficient to specify dorsal. Taken together, our data demonstrate that zygotic LvBMP5-8 signaling is essential for left-side specification, and for normal left-side skeletal and neural patterning, but not for DV specification. Thus, while both BMP2/4 and BMP5-8 regulate LR axis specification, BMP2/4 but not zygotic BMP5-8 regulates DV axis specification in sea urchin embryos.

RNA-Seq identifies SPGs as a ventral skeletal patterning cue in sea urchins.

The sea urchin larval skeleton offers a simple model for formation of developmental patterns. The calcium carbonate skeleton is secreted by primary mesenchyme cells (PMCs) in response to largely unknown patterning cues expressed by the ectoderm. To discover novel ectodermal cues, we performed an unbiased RNA-Seq-based screen and functionally tested candidates; we thereby identified several novel skeletal patterning cues. Among these, we show that SLC26a2/7 is a ventrally expressed sulfate transporter that promotes a ventral accumulation of sulfated proteoglycans, which is required for ventral PMC positioning and skeletal patterning. We show that the effects of SLC perturbation are mimicked by manipulation of either external sulfate levels or proteoglycan sulfation. These results identify novel skeletal patterning genes and demonstrate that ventral proteoglycan sulfation serves as a positional cue for sea urchin skeletal patterning.