ENTRAIN: integrating trajectory inference and gene regulatory networks with spatial data to co-localize the receptor-ligand interactions that specify cell fate.

MOTIVATION: Cell fate is commonly studied by profiling the gene expression of single cells to infer developmental trajectories based on expression similarity, RNA velocity, or statistical mechanical properties. However, current approaches do not recover microenvironmental signals from the cellular niche that drive a differentiation trajectory. RESULTS: We resolve this with environment-aware trajectory inference (ENTRAIN), a computational method that integrates trajectory inference methods with ligand-receptor pair gene regulatory networks to identify extracellular signals and evaluate their relative contribution towards a differentiation trajectory. The output from ENTRAIN can be superimposed on spatial data to co-localize cells and molecules in space and time to map cell fate potentials to cell-cell interactions. We validate and benchmark our approach on single-cell bone marrow and spatially resolved embryonic neurogenesis datasets to identify known and novel environmental drivers of cellular differentiation. AVAILABILITY AND IMPLEMENTATION: ENTRAIN is available as a public package at https://github.com/theimagelab/entrain and can be used on both single-cell and spatially resolved datasets.

Minimally invasive longitudinal intravital imaging of cellular dynamics in intact long bone.

Intravital two-photon microscopy enables deep-tissue imaging at high temporospatial resolution in live animals. However, the endosteal bone compartment and underlying bone marrow pose unique challenges to optical imaging as light is absorbed, scattered and dispersed by thick mineralized bone matrix and the adipose-rich bone marrow. Early bone intravital imaging methods exploited gaps in the cranial sutures to bypass the need to penetrate through cortical bone. More recently, investigators have developed invasive methods to thin the cortical bone or implant imaging windows to image cellular dynamics in weight-bearing long bones. Here, we provide a step-by-step procedure for the preparation of animals for minimally invasive, nondestructive, longitudinal intravital imaging of the murine tibia. This method involves the use of mixed bone marrow radiation chimeras to unambiguously double-label osteoclasts and osteomorphs. The tibia is exposed by a simple skin incision and an imaging chamber constructed using thermoconductive T-putty. Imaging sessions up to 12 h long can be repeated over multiple timepoints to provide a longitudinal time window into the endosteal and marrow niches. The approach can be used to investigate cellular dynamics in bone remodeling, cancer cell life cycle and hematopoiesis, as well as long-lived humoral and cellular immunity. The procedure requires an hour to complete and is suitable for users with minimal prior expertise in small animal surgery.

Apoptotic cell fragments locally activate tingible body macrophages in the germinal center.

Germinal centers (GCs) that form within lymphoid follicles during antibody responses are sites of massive cell death. Tingible body macrophages (TBMs) are tasked with apoptotic cell clearance to prevent secondary necrosis and autoimmune activation by intracellular self antigens. We show by multiple redundant and complementary methods that TBMs derive from a lymph node-resident, CD169-lineage, CSF1R-blockade-resistant precursor that is prepositioned in the follicle. Non-migratory TBMs use cytoplasmic processes to chase and capture migrating dead cell fragments using a "lazy" search strategy. Follicular macrophages activated by the presence of nearby apoptotic cells can mature into TBMs in the absence of GCs. Single-cell transcriptomics identified a TBM cell cluster in immunized lymph nodes which upregulated genes involved in apoptotic cell clearance. Thus, apoptotic B cells in early GCs trigger activation and maturation of follicular macrophages into classical TBMs to clear apoptotic debris and prevent antibody-mediated autoimmune diseases.

Tracking the clonal dynamics of SARS-CoV-2-specific T cells in children and adults with mild/asymptomatic COVID-19.

Children infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) develop less severe coronavirus disease 2019 (COVID-19) than adults. The mechanisms for the age-specific differences and the implications for infection-induced immunity are beginning to be uncovered. We show by longitudinal multimodal analysis that SARS-CoV-2 leaves a small footprint in the circulating T cell compartment in children with mild/asymptomatic COVID-19 compared to adult household contacts with the same disease severity who had more evidence of systemic T cell interferon activation, cytotoxicity and exhaustion. Children harbored diverse polyclonal SARS-CoV-2-specific naïve T cells whereas adults harbored clonally expanded SARS-CoV-2-specific memory T cells. A novel population of naïve interferon-activated T cells is expanded in acute COVID-19 and is recruited into the memory compartment during convalescence in adults but not children. This was associated with the development of robust CD4+ memory T cell responses in adults but not children. These data suggest that rapid clearance of SARS-CoV-2 in children may compromise their cellular immunity and ability to resist reinfection.

High titre neutralizing antibodies in response to SARS-CoV-2 infection require RBD-specific CD4 T cells that include proliferative memory cells.

Background: Long-term immunity to SARS-CoV-2 infection, including neutralizing antibodies and T cell-mediated immunity, is required in a very large majority of the population in order to reduce ongoing disease burden. Methods: We have investigated the association between memory CD4 and CD8 T cells and levels of neutralizing antibodies in convalescent COVID-19 subjects. Findings: Higher titres of convalescent neutralizing antibodies were associated with significantly higher levels of RBD-specific CD4 T cells, including specific memory cells that proliferated vigorously in vitro. Conversely, up to half of convalescent individuals had low neutralizing antibody titres together with a lack of receptor binding domain (RBD)-specific memory CD4 T cells. These low antibody subjects had other, non-RBD, spike-specific CD4 T cells, but with more of an inhibitory Foxp3+ and CTLA-4+ cell phenotype, in contrast to the effector T-bet+, cytotoxic granzymes+ and perforin+ cells seen in RBD-specific memory CD4 T cells from high antibody subjects. Single cell transcriptomics of antigen-specific CD4+ T cells from high antibody subjects similarly revealed heterogenous RBD-specific CD4+ T cells that comprised central memory, transitional memory and Tregs, as well as cytotoxic clusters containing diverse TCR repertoires, in individuals with high antibody levels. However, vaccination of low antibody convalescent individuals led to a slight but significant improvement in RBD-specific memory CD4 T cells and increased neutralizing antibody titres. Interpretation: Our results suggest that targeting CD4 T cell epitopes proximal to and within the RBD-region should be prioritized in booster vaccines.

Gene expression predicts dormant metastatic breast cancer cell phenotype.

BACKGROUND: Breast cancer can recur months to decades after an initial diagnosis and treatment. The mechanisms that control tumor cell dormancy remain poorly understood, making it difficult to predict which patients will recur and thus benefit from more rigorous screening and treatments. Unfortunately, the extreme rarity of dormant DTCs has been a major obstacle to their study. METHODS: To overcome this challenge, we developed an efficient system to isolate and study rare dormant breast cancer cells from metastatic organs including bones, which represent a major site of metastasis. After isolation of cells from the long bones, we used single cell RNA-sequencing (scRNA-seq) to profile proliferative and dormant PyMT-Bo1 breast cancer cells. We also compared this signature to dormant versus proliferative tumor cells isolated from the lungs. Finally, we compared our dormant signature to human datasets. RESULTS: We identified a group of genes including Cfh, Gas6, Mme and Ogn that were highly expressed in dormant breast cancer cells present in the bone and lung. Expression of these genes had no impact on dormancy in murine models, but their expression correlated with disease-free survival in primary human breast cancer tumors, suggesting that these genes have predictive value in determining which patients are likely to recur. CONCLUSIONS: Dormant breast cancer cells exhibit a distinct gene expression signature regardless of metastatic site. Genes enriched in dormant breast cancer cells correlate with recurrence-free survival in breast cancer patients.

Osteoclasts recycle via osteomorphs during RANKL-stimulated bone resorption.

Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases.

Prostate cancer cell-intrinsic interferon signaling regulates dormancy and metastatic outgrowth in bone.

The latency associated with bone metastasis emergence in castrate-resistant prostate cancer is attributed to dormancy, a state in which cancer cells persist prior to overt lesion formation. Using single-cell transcriptomics and ex vivo profiling, we have uncovered the critical role of tumor-intrinsic immune signaling in the retention of cancer cell dormancy. We demonstrate that loss of tumor-intrinsic type I IFN occurs in proliferating prostate cancer cells in bone. This loss suppresses tumor immunogenicity and therapeutic response and promotes bone cell activation to drive cancer progression. Restoration of tumor-intrinsic IFN signaling by HDAC inhibition increased tumor cell visibility, promoted long-term antitumor immunity, and blocked cancer growth in bone. Key findings were validated in patients, including loss of tumor-intrinsic IFN signaling and immunogenicity in bone metastases compared to primary tumors. Data herein provide a rationale as to why current immunotherapeutics fail in bone-metastatic prostate cancer, and provide a new therapeutic strategy to overcome the inefficacy of immune-based therapies in solid cancers.

A niche-dependent myeloid transcriptome signature defines dormant myeloma cells.

The era of targeted therapies has seen significant improvements in depth of response, progression-free survival, and overall survival for patients with multiple myeloma. Despite these improvements in clinical outcome, patients inevitably relapse and require further treatment. Drug-resistant dormant myeloma cells that reside in specific niches within the skeleton are considered a basis of disease relapse but remain elusive and difficult to study. Here, we developed a method to sequence the transcriptome of individual dormant myeloma cells from the bones of tumor-bearing mice. Our analyses show that dormant myeloma cells express a distinct transcriptome signature enriched for immune genes and, unexpectedly, genes associated with myeloid cell differentiation. These genes were switched on by coculture with osteoblastic cells. Targeting AXL, a gene highly expressed by dormant cells, using small-molecule inhibitors released cells from dormancy and promoted their proliferation. Analysis of the expression of AXL and coregulated genes in human cohorts showed that healthy human controls and patients with monoclonal gammopathy of uncertain significance expressed higher levels of the dormancy signature genes than patients with multiple myeloma. Furthermore, in patients with multiple myeloma, the expression of this myeloid transcriptome signature translated into a twofold increase in overall survival, indicating that this dormancy signature may be a marker of disease progression. Thus, engagement of myeloma cells with the osteoblastic niche induces expression of a suite of myeloid genes that predicts disease progression and that comprises potential drug targets to eradicate dormant myeloma cells.

Memory B cells are reactivated in subcapsular proliferative foci of lymph nodes.

Vaccine-induced immunity depends on the generation of memory B cells (MBC). However, where and how MBCs are reactivated to make neutralising antibodies remain unknown. Here we show that MBCs are prepositioned in a subcapsular niche in lymph nodes where, upon reactivation by antigen, they rapidly proliferate and differentiate into antibody-secreting plasma cells in the subcapsular proliferative foci (SPF). This novel structure is enriched for signals provided by T follicular helper cells and antigen-presenting subcapsular sinus macrophages. Compared with contemporaneous secondary germinal centres, SPF have distinct single-cell molecular signature, cell migration pattern and plasma cell output. Moreover, SPF are found both in human and mouse lymph nodes, suggesting that they are conserved throughout mammalian evolution. Our data thus reveal that SPF is a seat of immunological memory that may be exploited to rapidly mobilise secondary antibody responses and improve vaccine efficacy.

Single Cell RNA Sequencing of Rare Immune Cell Populations.

Single-cell RNA sequencing (scRNA-Seq) is transforming our ability to characterize cells, particularly rare cells that are often overlooked in bulk population analytical approaches. This has lead to the discovery of new cell types and cellular states that echo the underlying heterogeneity and plasticity in the immune system. Technologies for the capture, sequencing, and bioinformatic analysis of single cells are rapidly improving, and scRNA-Seq is now becoming much more accessible to non-specialized laboratories. Here, we describe our experiences in adopting scRNA-Seq to the study of rare immune cells in their microanatomical niches.