Characterizing the complex composition of solid tumors is fundamental for understanding tumor initiation, progression and metastasis. While patient-derived samples provide valuable insight, they are heterogeneous on multiple molecular levels, and often originate from advanced tumor stages. Here, we use single-cell transcriptome and epitope profiling together with pathway and lineage analyses to study tumorigenesis from a developmental perspective in a mouse model of salivary gland squamous cell carcinoma. We provide a comprehensive cell atlas and characterize tumor-specific cells. We find that these cells are connected along a reproducible developmental trajectory: initiated in basal cells exhibiting an epithelial-to-mesenchymal transition signature, tumorigenesis proceeds through Wnt-differential cancer stem cell-like subpopulations before differentiating into luminal-like cells. Our work provides unbiased insights into tumor-specific cellular identities in a whole tissue environment, and emphasizes the power of using defined genetic model systems.

Vaccines play a critical role in combating the COVID-19 pandemic. Future control of the pandemic requires improved vaccines with high efficacy against newly emerging SARS-CoV-2 variants and the ability to reduce virus transmission. Here we compare immune responses and preclinical efficacy of the mRNA vaccine BNT162b2, the adenovirus-vectored spike vaccine Ad2-spike and the live-attenuated virus vaccine candidate sCPD9 in Syrian hamsters, using both homogeneous and heterologous vaccination regimens. Comparative vaccine efficacy was assessed by employing readouts from virus titrations to single-cell RNA sequencing. Our results show that sCPD9 vaccination elicited the most robust immunity, including rapid viral clearance, reduced tissue damage, fast differentiation of pre-plasmablasts, strong systemic and mucosal humoral responses, and rapid recall of memory T cells from lung tissue after challenge with heterologous SARS-CoV-2. Overall, our results demonstrate that live-attenuated vaccines offer advantages over currently available COVID-19 vaccines.

We have previously reported Lvm1 as a quantitative trait locus (QTL) on chromosome 13 that links to cardiac left ventricular mass (LVM) in a panel of AxB/BxA mouse recombinant inbred strains (RIS). When performing a gene expression QTL (eQTL) analysis, we detected 33 cis-eQTLs that correlated with LVM. Among the latter, a group of eight cis-eQTLs clustered in a genomic region smaller than 6 Mb and surrounding the Lvm1 peak on chr13. Co-variant analysis indicated that all eight genes correlated with the phenotype in a causal rather than a reactive fashion, a finding that (despite its functional interest) did not provide grounds to prioritize any of these candidate genes. As a complementary approach, we performed weighted gene co-expression network analysis, which allowed us to detect 49 modules of highly connected genes. The module that correlated best with LVM: (1) showed linkage to a module QTL whose boundaries matched closely those of the phenotypic Lvm1 QTL on chr13; (2) harbored a disproportionately high proportion of genes originating from a small genomic region on chromosome 13 (including the 8 previously detected cis-eQTL genes); (3) contained genes that, beyond their individual level of expression, correlated with LVM as a function of their inter-connectivity; and (4) showed increased abundance of polymorphic insertion-deletion elements in the same region. Taken together, these data suggest that a domain on chromosome 13 constitutes the biologic principle responsible for the organization and linkage of the gene co-expression module, and indicate a mechanism whereby genetic variants within chromosome domains may associate to phenotypic changes via coordinate changes in the expression of several genes. One other possible implication of these findings is that candidate genes to consider as contributors to a particular phenotype should extend further than those that are closest to the QTL peak.

Little is known about the functions of chromosome Y (chrY) genes beyond their effects on sex and reproduction. In hearts, postpubertal testosterone affects the size of cells and the expression of genes differently in male C57BL/6J than in their C57.Y(A) counterparts, where the original chrY has been substituted with that from A/J mice. We further compared the 2 strains to better understand how chrY polymorphisms may affect cardiac properties, the latter being sexually dimorphic but unrelated to sex and reproduction. Genomic regions showing occupancy with androgen receptors (ARs) were identified in adult male hearts from both strains by chromatin immunoprecipitation. AR chromatin immunoprecipitation peaks (showing significant enrichment for consensus AR binding sites) were mostly strain specific. Measurements of anogenital distances in male pups showed that the biologic effects of perinatal androgens were greater in C57BL/6J than in C57.Y(A). Although perinatal endocrine manipulations showed that these differences contributed to the strain-specific differences in the response of adult cardiac cells to testosterone, the amounts of androgens produced by fetal testes were not different in each strain. Nonetheless, chrY polymorphisms associated in newborn pups' hearts with strain-specific differences in genomic regions showing either AR occupancy, accessible chromatin sites, or trimethylation of histone H3 Lysine 4 marks, as well as with differential expression of 2 chrY-encoded histone demethylases. In conclusion, the effects of chrY on adult cardiac phenotypes appeared to result from an interaction of this chromosome with the organizational programming effects exerted by the neonatal testosterone surge and show several characteristics of being mediated by an epigenetic remodeling of chromatin.

The C57BL/6J.YA/J mouse strain is a chromosome-substituted line where the original male-specific portion of chromosome Y (MSY) from C57BL/6J mice was substituted for that from A/J mice. In hearts from male C57BL/6J.YA/J and C57BL/6J mice, orchidectomy (ORX) affected in a strictly strain-specific fashion the expression a subset of genes showing enrichment for functional categories, including that of circadian rhythms and cardiac contractility. We further tested whether: (1) there were strain-specific differences in cardiac circadian rhythms; (2) strain-dependent differences in the effects of ORX on contractility genes translated into differences in cardiac functions; and (3) differential contractility responses occurred preferentially at times when circadian rhythms also showed strain-specific differences.

Natural sequence variation within mitochondrial DNA (mtDNA) contributes to human phenotypes and may serve as natural genetic markers in human cells for clonal and lineage tracing. We recently developed a single-cell multi-omic approach, called 'mitochondrial single-cell assay for transposase-accessible chromatin with sequencing' (mtscATAC-seq), enabling concomitant high-throughput mtDNA genotyping and accessible chromatin profiling. Specifically, our technique allows the mitochondrial genome-wide inference of mtDNA variant heteroplasmy along with information on cell state and accessible chromatin variation in individual cells. Leveraging somatic mtDNA mutations, our method further enables inference of clonal relationships among native ex vivo-derived human cells not amenable to genetic engineering-based clonal tracing approaches. Here, we provide a step-by-step protocol for the use of mtscATAC-seq, including various cell-processing and flow cytometry workflows, by using primary hematopoietic cells, subsequent single-cell genomic library preparation and sequencing that collectively take ~3-4 days to complete. We discuss experimental and computational data quality control metrics and considerations for the extension to other mammalian tissues. Overall, mtscATAC-seq provides a broadly applicable platform to map clonal relationships between cells in human tissues, investigate fundamental aspects of mitochondrial genetics and enable additional modes of multi-omic discovery.

Hedgehog (Hh) signaling plays a key role in the development and maintenance of animal tissues. This signaling is mediated by the atypical G protein-coupled receptor (GPCR) Smoothened (Smo). Smo activation leads to signaling through several well-characterized effectors to activate Hh target gene expression. Recent studies have implicated activation of the heterotrimeric G protein subunit Gαi and the subsequent decrease in cellular cAMP levels in promoting the Hh response in flies and mammals. Although Hh stimulation decreases cAMP levels in some insect cell lines, here using a bioluminescence resonance energy transfer (BRET)-based assay we found that this stimulation had no detectable effect in Drosophila S2-R+ cells. However, we observed an unexpected and significant Gαs-dependent increase in cAMP levels in response to strong Smo activation in Smo-transfected cells. This effect was mediated by Smo's broadly conserved core, and was specifically activated in response to phosphorylation of the Smo C-terminus by GPCR kinase 2 (Gprk2). Genetic analysis of heterotrimeric G protein function in the developing Drosophila wing revealed a positive role for cAMP in the endogenous Hh response. Specifically, we found that mutation or depletion of Gαs diminished low-threshold Hh responses in Drosophila, whereas depletion of Gαi potentiated them (in contrast to previous findings). Our analysis suggested that regulated cAMP production is important for controlling the sensitivity of cellular responses to Hh in Drosophila.

Recent developments in droplet-based microfluidics allow the transcriptional profiling of thousands of individual cells in a quantitative, highly parallel and cost-effective way. A critical, often limiting step is the preparation of cells in an unperturbed state, not altered by stress or ageing. Other challenges are rare cells that need to be collected over several days or samples prepared at different times or locations.