Understanding lung immunity requires an unbiased profiling of tissue-resident T cells at their precise anatomical locations within the lung, but such information has not been characterized in the immunized mouse model. In this pilot study, using 10x Genomics Chromium and Visium platform, we performed an integrative analysis of spatial transcriptome with single-cell RNA-seq and single-cell ATAC-seq on lung cells from mice after immunization using a well-established Klebsiella pneumoniae infection model. We built an optimized deconvolution pipeline to accurately decipher specific cell-type compositions by anatomic location. We discovered that combining scATAC-seq and scRNA-seq data may provide more robust cell-type identification, especially for lineage-specific T helper cells. Combining all three modalities, we observed a dynamic change in the location of T helper cells as well as their corresponding chemokines. In summary, our proof-of-principle study demonstrated the power and potential of single-cell multi-omics analysis to uncover spatial- and cell-type-dependent mechanisms of lung immunity.

Coxsackievirus B3 (CVB3)-induced myocarditis is commonly employed to study viral pathogenesis in mice. Chronically affected mice may develop dilated cardiomyopathy, which may involve the mediation of immune and nonimmune cells. To dissect this complexity, we performed single-cell RNA sequencing on heart cells from healthy and myocarditic mice, leading us to note significant proportions of myeloid cells, T cells, and fibroblasts. Although the transcriptomes of myeloid cells were mainly of M2 phenotype, the Th17 cells, CTLs, and Treg cells had signatures critical for cytotoxic functions. Fibroblasts were heterogeneous expressing genes important in fibrosis and regulation of inflammation and immune responses. The intercellular communication networks revealed unique interactions and signaling pathways in the cardiac cellulome, whereas myeloid cells and T cells had upregulated unique transcription factors modulating cardiac remodeling functions. Together, our data suggest that M2 cells, T cells, and fibroblasts may cooperatively or independently participate in the pathogenesis of viral myocarditis.

Although new genomics-based pipelines have potential to augment antibody discovery, these methods remain in their infancy due to an incomplete understanding of the selection process that governs B cell clonal selection, expansion, and antigen specificity. Furthermore, it remains unknown how factors such as aging and reduction of tolerance influence B cell selection. Here we perform single-cell sequencing of antibody repertoires and transcriptomes of murine B cells following immunizations with a model therapeutic antigen target. We determine the relationship between antibody repertoires, gene expression signatures, and antigen specificity across 100,000 B cells. Recombinant expression and characterization of 227 monoclonal antibodies revealed the existence of clonally expanded and class-switched antigen-specific B cells that were more frequent in young mice. Although integrating multiple repertoire features such as germline gene usage and transcriptional signatures failed to distinguish antigen-specific from nonspecific B cells, other features such as immunoglobulin G (IgG) subtype and sequence composition correlated with antigen specificity.

Thermal fluctuation of local magnetization intercoupled with charge carriers and phonons offers a path to enhance thermoelectric performance. Thermopower enhancement by spin fluctuations (SF) has been observed before. However, the crucial evidence for enhancing thermoelectric-figure-of-merit (zT) by SF has not been reported until now. Here we report that the SF leads to nearly 80% zT enhancement in ferromagnetic CrTe near and below T C ∼ 335 K. The ferromagnetism is originated from the collective electronic and localized magnetic moments. The field-dependent transport properties demonstrate the profound impact of SF on the electrons and phonons. Under an external magnetic field, the enhancement in thermopower is suppressed, and the thermal conductivity is enhanced, evidencing the existence of a strong SF. The anomalous thermoelectric transport properties are analyzed based on theoretical models, and a good agreement with experimental data is found. This study contributes to the fundamental understanding of SF for designing high-performance spin-driven thermoelectrics.

Two-dimensional (2D) van der Waals (vdW) ferromagnets have opened new avenues for manipulating spin at the limits of single or few atomic layers, and for creating unique magneto-exciton devices through the coupling of ferromagnetic (FM) orders and excitons. However, 2D vdW ferromagnets explored so far have rarely possessed exciton behaviors; to date, FM CrI3 have been revealed to show ligand-field photoluminescence correlated with FM ordering, but typically with a broad emission peak. Here, we report a straightforward approach to realize strong coupling of narrow helical emission and FM orders in CrI3 through microsphere cavity. The resonant whispering-gallery modes (WGM) of SiO2 microspheres cause strong oscillation helical emissions with a full width at half-maximum (FWHM) of ∼5 nm under continuous wave excitation. Reversible magnetic coding of helical luminescence is realized in the range of 950-1100 nm. This work enables numerous opportunities for creating magnetic encoding lasing for photonic integrated chips.

Inner ear organoids recapitulate development and are intended to generate cell types of the otic lineage for applications such as basic science research and cell replacement strategies. Here, we use single-cell sequencing to study the cellular heterogeneity of late-stage mouse inner ear organoid sensory epithelia, which we validated by comparison with datasets of the mouse cochlea and vestibular epithelia. We resolved supporting cell sub-types, cochlear-like hair cells, and vestibular type I and type II-like hair cells. While cochlear-like hair cells aligned best with an outer hair cell trajectory, vestibular-like hair cells followed developmental trajectories similar to in vivo programs branching into type II and then type I extrastriolar hair cells. These results highlight the transcriptional accuracy of the organoid developmental program but will also inform future strategies to improve synaptic connectivity and regional specification.

Bladder cancer (BLCA) is more common in men but more aggressive in women. Sex-based differences in cancer biology are commonly studied using a murine model with BLCA generated by N-butyl-N-(4-hydroxybutyl)-nitrosamine (BBN). While tumors in the BBN model have been profiled, these profiles provide limited information on the tumor microenvironment. Here, we applied single-cell RNA sequencing to characterize cell-type specific transcriptional differences between male and female BBN-induced tumors. We found proportional and gene expression differences in epithelial and non-epithelial subpopulations between male and female tumors. Expression of several genes predicted sex-specific survival in several human BLCA datasets. We identified novel and clinically relevant sex-specific transcriptional signatures including immune cells in the tumor microenvironment and it validated the relevance of the BBN model for studying sex differences in human BLCA. This work highlights the importance of considering sex as a biological variable in the development of new and accurate cancer markers.

The analysis and interpretation of single-cell RNA sequencing (scRNA-seq) experiments are compromised by the presence of poor-quality cells. For meaningful analyses, such poor-quality cells should be excluded as they introduce noise in the data. We introduce SkewC, a quality-assessment tool, to identify skewed cells in scRNA-seq experiments. The tool's methodology is based on the assessment of gene coverage for each cell, and its skewness as a quality measure; the gene body coverage is a unique characteristic for each protocol, and different protocols yield highly different coverage profiles. This tool is designed to avoid misclustering or false clusters by identifying, isolating, and removing cells with skewed gene body coverage profiles. SkewC is capable of processing any type of scRNA-seq dataset, regardless of the protocol. We envision SkewC as a distinctive QC method to be incorporated into scRNA-seq QC processing to preclude the possibility of scRNA-seq data misinterpretation.

Human plasmacytoid dendritic cells (pDCs) play a central role in initiating and activating host immune responses during infection. To understand how the transcriptome of pDCs is impacted by HIV-1 infection and exogenous stimulation, we isolated pDCs from healthy controls, people with HIV-1 (PWH) before and during toll-like receptor 9 (TLR9) agonist treatment and performed single-cell (sc)-RNA sequencing. Our cluster analysis revealed four pDC clusters: pDC1, pDC2, cytotoxic-like pDC and an exhausted pDC cluster. The inducible cytotoxic-like pDC cluster is characterized by high expression of both antiviral and cytotoxic genes. Further analyses confirmed that cytotoxic-like pDCs are distinct from NK and T cells. Cell-cell communication analysis also demonstrated that cytotoxic-like pDCs exhibit similar incoming and outgoing cellular communicating signals as other pDCs. Thus, our study presents a detailed transcriptomic atlas of pDCs and provides new perspectives on the mechanisms of regulation and function of cytotoxic-like pDCs.

To gain deeper insights into transcriptomes and epigenomes of organoids, liver organoids from two states (expandable and more differentiated) were subjected to single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) analyses. Mitochondrial gene expression was higher in differentiated than in non-differentiated hepatocytes, with ATAC-seq peaks increasing near the mitochondrial control region. Differentiation of liver organoids resulted in the expression of transcription factors that act as enhancers and repressors. In addition, epigenetic mechanisms regulating the expression of alpha-fetoprotein (AFP) and albumin (ALB) differed in liver organoids and adult liver. Knockdown of PDX1, an essential transcription factor for pancreas development, led to the hepatic maturation of liver organoids through regulation of AFP and ALB expression. This integrative analysis of the transcriptomes and epigenomes of liver organoids at the single-cell level may contribute to a better understanding of the regulatory networks during liver development and the further development of mature in vitro human liver models.

Multisystem inflammatory syndrome in children (MIS-C) is a life-threatening disease occurring several weeks after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Deep immune profiling showed acute MIS-C patients had highly activated neutrophils, classical monocytes and memory CD8+ T-cells, with increased frequencies of B-cell plasmablasts and double-negative B-cells. Post treatment samples from the same patients, taken during symptom resolution, identified recovery-associated immune features including increased monocyte CD163 levels, emergence of a new population of immature neutrophils and, in some patients, transiently increased plasma arginase. Plasma profiling identified multiple features shared by MIS-C, Kawasaki Disease and COVID-19 and that therapeutic inhibition of IL-6 may be preferable to IL-1 or TNF-α. We identified several potential mechanisms of action for IVIG, the most commonly used drug to treat MIS-C. Finally, we showed systemic complement activation with high plasma C5b-9 levels is common in MIS-C suggesting complement inhibitors could be used to treat the disease.

The retrosplenial cortex (RSC) is an area interconnected with regions of the brain that display spatial correlates. Neurons in connected regions may encode an animal's position in the environment and location or proximity to objects or boundaries. RSC has also been shown to be important for spatial memory, such as tracking distance from and between landmarks, contextual information, and orientation within an environment. For these reasons, it is important to determine how neurons in RSC represent cues such as objects or boundaries and their relationship to the environment. In the current work, we performed electrophysiological recordings in RSC, whereas rats foraged in arenas that could contain an object or in which the environment was altered. We report RSC neurons display changes in mean firing rate responding to alterations of the environment. These alterations include the arena rotating, changing size or shape, or an object being introduced into the arena.

SARS-CoV-2 infection and vaccination elicit potent immune responses. Our study presents a comprehensive multimodal single-cell analysis of blood from COVID-19 patients and healthy volunteers receiving the SARS-CoV-2 vaccine and booster. We profiled immune responses via transcriptional analysis and lymphocyte repertoire reconstruction. COVID-19 patients displayed an enhanced interferon signature and cytotoxic gene upregulation, absent in vaccine recipients. B and T cell repertoire analysis revealed clonal expansion among effector cells in COVID-19 patients and memory cells in vaccine recipients. Furthermore, while clonal αβ T cell responses were observed in both COVID-19 patients and vaccine recipients, expansion of clonal γδ T cells was found only in infected individuals. Our dataset enables side-by-side comparison of immune responses to infection versus vaccination, including clonal B and T cell responses. Our comparative analysis shows that vaccination induces a robust, durable clonal B and T cell responses, without the severe inflammation associated with infection.

Mutations in the gene encoding DNA methyltransferase 3A (DNMT3A) are the most common cause of clonal hematopoiesis and are among the most common initiating events of acute myeloid leukemia (AML). Studies in germline and somatic Dnmt3a knockout mice have identified focal, canonical hypomethylation phenotypes in hematopoietic cells; however, the kinetics of methylation loss following acquired DNMT3A inactivation in hematopoietic cells is essentially unknown. Therefore, we evaluated a somatic, inducible model of hematopoietic Dnmt3a loss, and show that inactivation of Dnmt3a in murine hematopoietic cells results in a relatively slow loss of methylation at canonical sites throughout the genome; in contrast, remethylation of Dnmt3a deficient genomes in hematopoietic cells occurs much more quickly. This data suggests that slow methylation loss may contribute, at least in part, to the long latent period that characterizes clonal expansion and leukemia development in individuals with acquired DNMT3A mutations in hematopoietic stem cells.

Somatic stem cells are advantageous research targets for understanding the properties required to maintain stemness. Human bone marrow-mesenchymal stromal cells (BM-MSCs) were separated into pluripotent-like SSEA-3(+) Muse cells (Muse-MSCs) and multipotent SSEA-3(-) MSCs (MSCs) and were subjected to single-cell RNA sequencing analysis. Compared with MSCs, Muse-MSCs exhibited higher expression levels of the p53 repressor MDM2; signal acceptance-related genes EGF, VEGF, PDGF, WNT, TGFB, INHB, and CSF; ribosomal protein; and glycolysis and oxidative phosphorylation. Conversely, MSCs had higher expression levels of FGF and ANGPT; Rho family and caveola-related genes; amino acid and cofactor metabolism; MHC class I/II, and lysosomal enzyme genes than Muse-MSCs. Unsupervised clustering further divided Muse-MSCs into two clusters stratified by the expression of cell cycle-related genes, and MSCs into three clusters stratified by the expression of cell cycle-, cytoskeleton-, and extracellular matrix-related genes. This study evaluating the differentiation ability of BM-MSC subpopulations provides intriguing insights for understanding stemness.

Deconstructing tissue-specific effects of genes and variants on proliferation is critical to understanding cellular transformation and systematically selecting cancer therapeutics. This requires scalable methods for multiplexed genetic screens tracking fitness across time, across lineages, and in a suitable niche, since physiological cues influence functional differences. Towards this, we present an approach, coupling single-cell cancer driver screens in teratomas with hit enrichment by serial teratoma reinjection, to simultaneously screen drivers across multiple lineages in vivo. Using this system, we analyzed population shifts and lineage-specific enrichment for 51 cancer associated genes and variants, profiling over 100,000 cells spanning over 20 lineages, across two rounds of serial reinjection. We confirmed that c-MYC alone or combined with myristoylated AKT1 potently drives proliferation in progenitor neural lineages, demonstrating signatures of malignancy. Additionally, mutant MEK1 S218D/S222D provides a proliferative advantage in mesenchymal lineages like fibroblasts. Our method provides a powerful platform for multi-lineage longitudinal study of oncogenesis.

Clinically important broadly reactive B cells evolve during multiple infections, with B cells re-activated after secondary infection differing from B cells activated after a primary infection. Here we studied CD27highCD38high plasmablasts from patients with a primary or secondary dengue virus infection. Three transcriptionally and functionally distinct clusters were identified. The largest cluster 0/1 was plasma cell-related, with cells coding for serotype cross-reactive antibodies of the IgG1 isotype, consistent with memory B cell activation during an extrafollicular response. Cells in clusters 2 and 3 expressed low levels of antibody genes and high levels of genes associated with oxidative phosphorylation, EIF2 pathway, and mitochondrial dysfunction. Clusters 2 and 3 showed a transcriptional footprint of T cell help, in line with activation from naive B cells or memory B cells. Our results contribute to the understanding of the parallel B cell activation events that occur in humans after natural primary and secondary infection.

Atypical B cells are a population of activated B cells that are commonly enriched in individuals with chronic immune activation but are also part of a normal immune response to infection or vaccination. To better define the role of atypical B cells in the human adaptive immune response, we performed single-cell sequencing of transcriptomes, cell surface markers, and B cell receptors in individuals with chronic exposure to the malaria parasite Plasmodium falciparum, a condition known to lead to accumulation of circulating atypical B cells. We identified three previously uncharacterized populations of atypical B cells with distinct transcriptional and functional profiles and observed marked differences among these three subsets in their ability to produce immunoglobulin G upon T-cell-dependent activation. Our findings help explain the conflicting observations in prior studies regarding the function of atypical B cells and highlight their different roles in the adaptive immune response in chronic inflammatory conditions.

Skeletal muscle is composed of post-mitotic myofibers that form a syncytium containing hundreds of myonuclei. Using a progressive exercise training model in the mouse and single nucleus RNA-sequencing (snRNA-seq) for high-resolution characterization of myonuclear transcription, we show myonuclear functional specialization in muscle. After 4 weeks of exercise training, snRNA-seq reveals that resident muscle stem cells, or satellite cells, are activated with acute exercise but demonstrate limited lineage progression while contributing to muscle adaptation. In the absence of satellite cells, a portion of nuclei demonstrates divergent transcriptional dynamics associated with mixed-fate identities compared with satellite cell replete muscles. These data provide a compendium of information about how satellite cells influence myonuclear transcription in response to exercise.

Little is known about the effects of high-fat diet (HFD)-induced obesity on resident colonic lamina propria (LP) macrophages (LPMs) function and metabolism. Here, we report that obesity and diabetes resulted in increased macrophage infiltration in the colon. These macrophages exhibited the residency phenotype CX3CR1hiMHCIIhi and were CD4-TIM4-. During HFD, resident colonic LPM exhibited a lipid metabolism gene expression signature that overlapped that used to define lipid-associated macrophages (LAMs). Via single-cell RNA sequencing, we identified a sub-cluster of macrophages, increased in HFD, that were responsible for the LAM signature. Compared to other macrophages in the colon, these cells were characterized by elevated glycolysis, phagocytosis, and efferocytosis signatures. CX3CR1hiMHCIIhi colonic resident LPMs had fewer lipid droplets (LDs) and decreased triacylglycerol (TG) content compared to equivalent cells in lean mice and exhibited increased phagocytic capacity, suggesting that HFD induces adaptive responses in LPMs to limit bacterial translocation.