Mapping the topography of spatial gene expression with interpretable deep learning.

Spatially resolved transcriptomics technologies provide high-throughput measurements of gene expression in a tissue slice, but the sparsity of this data complicates the analysis of spatial gene expression patterns such as gene expression gradients. We address these issues by deriving a topographic map of a tissue slice-analogous to a map of elevation in a landscape-using a novel quantity called the isodepth. Contours of constant isodepth enclose spatial domains with distinct cell type composition, while gradients of the isodepth indicate spatial directions of maximum change in gene expression. We develop GASTON, an unsupervised and interpretable deep learning algorithm that simultaneously learns the isodepth, spatial gene expression gradients, and piecewise linear functions of the isodepth that model both continuous gradients and discontinuous spatial variation in the expression of individual genes. We validate GASTON by showing that it accurately identifies spatial domains and marker genes across several biological systems. In SRT data from the brain, GASTON reveals gradients of neuronal differentiation and firing, and in SRT data from a tumor sample, GASTON infers gradients of metabolic activity and epithelial-mesenchymal transition (EMT)-related gene expression in the tumor microenvironment.

The landscape of alternative polyadenylation in single cells of the developing mouse embryo.

3' untranslated regions (3' UTRs) post-transcriptionally regulate mRNA stability, localization, and translation rate. While 3'-UTR isoforms have been globally quantified in limited cell types using bulk measurements, their differential usage among cell types during mammalian development remains poorly characterized. In this study, we examine a dataset comprising ~2 million nuclei spanning E9.5-E13.5 of mouse embryonic development to quantify transcriptome-wide changes in alternative polyadenylation (APA). We observe a global lengthening of 3' UTRs across embryonic stages in all cell types, although we detect shorter 3' UTRs in hematopoietic lineages and longer 3' UTRs in neuronal cell types within each stage. An analysis of RNA-binding protein (RBP) dynamics identifies ELAV-like family members, which are concomitantly induced in neuronal lineages and developmental stages experiencing 3'-UTR lengthening, as putative regulators of APA. By measuring 3'-UTR isoforms in an expansive single cell dataset, our work provides a transcriptome-wide and organism-wide map of the dynamic landscape of alternative polyadenylation during mammalian organogenesis.

Functional Genetic Screen to Identify Interneurons Governing Behaviorally Distinct Aspects of Drosophila Larval Motor Programs.

Drosophila larval crawling is an attractive system to study rhythmic motor output at the level of animal behavior. Larval crawling consists of waves of muscle contractions generating forward or reverse locomotion. In addition, larvae undergo additional behaviors, including head casts, turning, and feeding. It is likely that some neurons (e.g., motor neurons) are used in all these behaviors, but the identity (or even existence) of neurons dedicated to specific aspects of behavior is unclear. To identify neurons that regulate specific aspects of larval locomotion, we performed a genetic screen to identify neurons that, when activated, could elicit distinct motor programs. We used 165 Janelia CRM-Gal4 lines-chosen for sparse neuronal expression-to ectopically express the warmth-inducible neuronal activator TrpA1, and screened for locomotor defects. The primary screen measured forward locomotion velocity, and we identified 63 lines that had locomotion velocities significantly slower than controls following TrpA1 activation (28°). A secondary screen was performed on these lines, revealing multiple discrete behavioral phenotypes, including slow forward locomotion, excessive reverse locomotion, excessive turning, excessive feeding, immobile, rigid paralysis, and delayed paralysis. While many of the Gal4 lines had motor, sensory, or muscle expression that may account for some or all of the phenotype, some lines showed specific expression in a sparse pattern of interneurons. Our results show that distinct motor programs utilize distinct subsets of interneurons, and provide an entry point for characterizing interneurons governing different elements of the larval motor program.