A machine learning one-class logistic regression model to predict stemness for single cell transcriptomics and spatial omics.

Cell annotation is a crucial methodological component to interpreting single cell and spatial omics data. These approaches were developed for single cell analysis but are often biased, manually curated and yet unproven in spatial omics. Here we apply a stemness model for assessing oncogenic states to single cell and spatial omic cancer datasets. This one-class logistic regression machine learning algorithm is used to extract transcriptomic features from non-transformed stem cells to identify dedifferentiated cell states in tumors. We found this method identifies single cell states in metastatic tumor cell populations without the requirement of cell annotation. This machine learning model identified stem-like cell populations not identified in single cell or spatial transcriptomic analysis using existing methods. For the first time, we demonstrate the application of a ML tool across five emerging spatial transcriptomic and proteomic technologies to identify oncogenic stem-like cell types in the tumor microenvironment.

High-plex imaging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging.

Resolving the spatial distribution of RNA and protein in tissues at subcellular resolution is a challenge in the field of spatial biology. We describe spatial molecular imaging, a system that measures RNAs and proteins in intact biological samples at subcellular resolution by performing multiple cycles of nucleic acid hybridization of fluorescent molecular barcodes. We demonstrate that spatial molecular imaging has high sensitivity (one or two copies per cell) and very low error rate (0.0092 false calls per cell) and background (~0.04 counts per cell). The imaging system generates three-dimensional, super-resolution localization of analytes at ~2 million cells per sample. Cell segmentation is morphology based using antibodies, compatible with formalin-fixed, paraffin-embedded samples. We measured multiomic data (980 RNAs and 108 proteins) at subcellular resolution in formalin-fixed, paraffin-embedded tissues (nonsmall cell lung and breast cancer) and identified >18 distinct cell types, ten unique tumor microenvironments and 100 pairwise ligand-receptor interactions. Data on >800,000 single cells and ~260 million transcripts can be accessed at http://nanostring.com/CosMx-dataset .

Developmental basis of evolutionary lung loss in plethodontid salamanders.

One or more members of four living amphibian clades have independently dispensed with pulmonary respiration and lack lungs, but little is known of the developmental basis of lung loss in any taxon. We use morphological, molecular, and experimental approaches to examine the Plethodontidae, a dominant family of salamanders, all of which are lungless as adults. We confirm an early anecdotal report that plethodontids complete early stages of lung morphogenesis: Transient embryonic lung primordia form but regress by apoptosis before hatching. Initiation of pulmonary development coincides with expression of the lung-specification gene Wnt2b in adjacent mesoderm, and the lung rudiment expresses pulmonary markers Nkx2.1 and Sox9. Lung developmental-genetic pathways are at least partially conserved despite the absence of functional adult lungs for at least 25 and possibly exceeding 60 million years. Adult lung loss appears associated with altered expression of signaling molecules that mediate later stages of tracheal and pulmonary development.

Gene co-expression reveals the modularity and integration of C4 and CAM in Portulaca.

C4 photosynthesis and Crassulacean acid metabolism (CAM) have been considered as largely independent adaptations despite sharing key biochemical modules. Portulaca is a geographically widespread clade of over 100 annual and perennial angiosperm species that primarily use C4 but facultatively exhibit CAM when drought stressed, a photosynthetic system known as C4 + CAM. It has been hypothesized that C4 + CAM is rare because of pleiotropic constraints, but these have not been deeply explored. We generated a chromosome-level genome assembly of Portulaca amilis and sampled mRNA from P. amilis and Portulaca oleracea during CAM induction. Gene co-expression network analyses identified C4 and CAM gene modules shared and unique to both Portulaca species. A conserved CAM module linked phosphoenolpyruvate carboxylase to starch turnover during the day-night transition and was enriched in circadian clock regulatory motifs in the P. amilis genome. Preservation of this co-expression module regardless of water status suggests that Portulaca constitutively operate a weak CAM cycle that is transcriptionally and posttranscriptionally upregulated during drought. C4 and CAM mostly used mutually exclusive genes for primary carbon fixation, and it is likely that nocturnal CAM malate stores are shuttled into diurnal C4 decarboxylation pathways, but we found evidence that metabolite cycling may occur at low levels. C4 likely evolved in Portulaca through co-option of redundant genes and integration of the diurnal portion of CAM. Thus, the ancestral CAM system did not strongly constrain C4 evolution because photosynthetic gene networks are not co-regulated for both daytime and nighttime functions.

Expression of a novel surfactant protein gene is associated with sites of extrapulmonary respiration in a lungless salamander.

Numerous physiological and morphological adaptations were achieved during the transition to lungless respiration that accompanied evolutionary lung loss in plethodontid salamanders, including those that enable efficient gas exchange across extrapulmonary tissue. However, the molecular basis of these adaptations is unknown. Here, we show that lungless salamanders express in the larval integument and the adult buccopharynx-principal sites of respiratory gas exchange in these species-a novel paralogue of the gene surfactant-associated protein C (SFTPC), which is a critical component of pulmonary surfactant expressed exclusively in the lung in other vertebrates. The paralogous gene appears to be found only in salamanders, but, similar to SFTPC, in lunged salamanders it is expressed only in the lung. This heterotopic gene expression, combined with predictions from structural modelling and respiratory tissue ultrastructure, suggests that lungless salamanders may produce pulmonary surfactant-like secretions outside the lungs and that the novel paralogue of SFTPC might facilitate extrapulmonary respiration in the absence of lungs. Heterotopic expression of the SFTPC paralogue may have contributed to the remarkable evolutionary radiation of lungless salamanders, which account for more than two thirds of urodele species alive today.

We are not so special.

New sequence data from choanoflagellates improves our understanding of the genetic changes that occurred along the branch of the evolutionary tree that gave rise to animals.

Convergent evolutionary reduction of atrial septation in lungless salamanders.

Nearly two thirds of the approximately 700 species of living salamanders are lungless. These species respire entirely through the skin and buccopharyngeal mucosa. Lung loss dramatically impacts the configuration of the circulatory system but the effects of evolutionary lung loss on cardiac morphology have long been controversial. For example, there is presumably little need for an atrial septum in lungless salamanders due to the absence of pulmonary veins and the presence of a single source of mixed blood flowing into the heart, but whether lungless salamanders possess an atrial septum and whether the sinoatrial aperture is located in the left or right atrium are unresolved; authors have stated opposing claims since the late 1800s. Here, we use micro-computed tomography (µ-CT) imaging, gross dissection and histological reconstruction to compare cardiac morphology among lungless plethodontid salamanders (Plethodontidae), salamanders with lungs, and the convergently lungless species Onychodactylus japonicus (Hynobiidae). Plethodontid salamanders have partial atrial septa and incomplete separation of the atrium into left and right halves. Partial septation is also seen in O. japonicus. Hence, lungless salamanders from two lineages convergently evolved similar morphology of the atrial septum. The partial septum in lungless salamanders can make it appear that the sinoatrial aperture is in the left atrium, but this interpretation is incorrect. Outgroup comparisons demonstrate that the aperture is located in a posterodorsal extension of the right atrium into the left side of the heart. Independent evolutionary losses of the atrial septum may have a similar developmental basis. In mammals, the lungs induce formation of the atrial septum by secreting morphogens to neighboring mesenchyme. We hypothesize that the lungs induce atrial septum development in amphibians in a similar fashion to mammals, and that atrial septum reduction in lungless salamanders is a direct result of lunglessness.

Graded levels of Pax2a and Pax8 regulate cell differentiation during sensory placode formation.

Pax gene haploinsufficiency causes a variety of congenital defects. Renal-coloboma syndrome, resulting from mutations in Pax2, is characterized by kidney hypoplasia, optic nerve malformation, and hearing loss. Although this underscores the importance of Pax gene dosage in normal development, how differential levels of these transcriptional regulators affect cell differentiation and tissue morphogenesis is still poorly understood. We show that differential levels of zebrafish Pax2a and Pax8 modulate commitment and behavior in cells that eventually contribute to the otic vesicle and epibranchial placodes. Initially, a subset of epibranchial placode precursors lie lateral to otic precursors within a single Pax2a/8-positive domain; these cells subsequently move to segregate into distinct placodes. Using lineage-tracing and ablation analyses, we show that cells in the Pax2a/8+ domain become biased towards certain fates at the beginning of somitogenesis. Experiments involving either Pax2a overexpression or partial, combinatorial Pax2a and Pax8 loss of function reveal that high levels of Pax favor otic differentiation whereas low levels increase cell numbers in epibranchial ganglia. In addition, the Fgf and Wnt signaling pathways control Pax2a expression: Fgf is necessary to induce Pax2a, whereas Wnt instructs the high levels of Pax2a that favor otic differentiation. Our studies reveal the importance of Pax levels during sensory placode formation and provide a mechanism by which these levels are controlled.

Chondrogenic and gliogenic subpopulations of neural crest play distinct roles during the assembly of epibranchial ganglia.

In vertebrates, the sensory neurons of the epibranchial (EB) ganglia transmit somatosensory signals from the periphery to the CNS. These ganglia are formed during embryogenesis by the convergence and condensation of two distinct populations of precursors: placode-derived neuroblasts and neural crest- (NC) derived glial precursors. In addition to the gliogenic crest, chondrogenic NC migrates into the pharyngeal arches, which lie in close proximity to the EB placodes and ganglia. Here, we examine the respective roles of these two distinct NC-derived populations during development of the EB ganglia using zebrafish morphant and mutants that lack one or both of these NC populations. Our analyses of mutant and morphant zebrafish that exhibit deficiencies in chondrogenic NC at early stages reveal a distinct requirement for this NC subpopulation during early EB ganglion assembly and segmentation. Furthermore, restoration of wildtype chondrogenic NC in one of these mutants, prdm1a, is sufficient to restore ganglion formation, indicating a specific requirement of the chondrogenic NC for EB ganglia assembly. By contrast, analysis of the sox10 mutant, which lacks gliogenic NC, reveals that the initial assembly of ganglia is not affected. However, during later stages of development, EB ganglia are dispersed in the sox10 mutant, suggesting that glia are required to maintain normal EB ganglion morphology. These results highlight novel roles for two subpopulations of NC cells in the formation and maintenance of EB ganglia: chondrogenic NC promotes the early-stage formation of the developing EB ganglia while glial NC is required for the late-stage maintenance of ganglion morphology.

Female-specific increase in primordial germ cells marks sex differentiation in threespine stickleback (Gasterosteus aculeatus).

Gonadal sex differentiation is increasingly recognized as a remarkably plastic process driven by species-specific genetic or environmental determinants. Among aquatic vertebrates, gonadal sex differentiation is a frequent endpoint in studies of endocrine disruption with little appreciation of underlying developmental mechanisms. Work in model organisms has highlighted the diversity of master sex-determining genes rather than uncovering any broad similarities prompting the highly conserved developmental decision of testes versus ovaries. Here we use molecular genetic markers of chromosomal sex combined with traditional histology to examine the transition of the bipotential gonads to ovaries or testes in threespine stickleback (Gasterosteus aculeatus). Serially-sectioned threespine stickleback fry were analyzed for qualitative and quantitative indications of sexual differentiation, including changes in gonadal morphology, number of germ cells and the incidence of gonadal apoptosis. We show that threespine stickleback sampled from anadromous and lacustrine populations are differentiated gonochorists. The earliest sex-specific event is a premeiotic increase in primordial germ cell number followed by a female-specific spike in apoptosis in the undifferentiated gonad of genetic females. The data suggest that an increase in PGC number may direct the undifferentiated gonad toward ovarian differentiation.