JAK-STAT signaling maintains homeostasis in T cells and macrophages.

Immune cells need to sustain a state of constant alertness over a lifetime. Yet, little is known about the regulatory processes that control the fluent and fragile balance that is called homeostasis. Here we demonstrate that JAK-STAT signaling, beyond its role in immune responses, is a major regulator of immune cell homeostasis. We investigated JAK-STAT-mediated transcription and chromatin accessibility across 12 mouse models, including knockouts of all STAT transcription factors and of the TYK2 kinase. Baseline JAK-STAT signaling was detected in CD8+ T cells and macrophages of unperturbed mice-but abrogated in the knockouts and in unstimulated immune cells deprived of their normal tissue context. We observed diverse gene-regulatory programs, including effects of STAT2 and IRF9 that were independent of STAT1. In summary, our large-scale dataset and integrative analysis of JAK-STAT mutant and wild-type mice uncovered a crucial role of JAK-STAT signaling in unstimulated immune cells, where it contributes to a poised epigenetic and transcriptional state and helps prepare these cells for rapid response to immune stimuli.

Structural cells are key regulators of organ-specific immune responses.

The mammalian immune system implements a remarkably effective set of mechanisms for fighting pathogens1. Its main components are haematopoietic immune cells, including myeloid cells that control innate immunity, and lymphoid cells that constitute adaptive immunity2. However, immune functions are not unique to haematopoietic cells, and many other cell types display basic mechanisms of pathogen defence3-5. To advance our understanding of immunology outside the haematopoietic system, here we systematically investigate the regulation of immune genes in the three major types of structural cells: epithelium, endothelium and fibroblasts. We characterize these cell types across twelve organs in mice, using cellular phenotyping, transcriptome sequencing, chromatin accessibility profiling and epigenome mapping. This comprehensive dataset revealed complex immune gene activity and regulation in structural cells. The observed patterns were highly organ-specific and seem to modulate the extensive interactions between structural cells and haematopoietic immune cells. Moreover, we identified an epigenetically encoded immune potential in structural cells under tissue homeostasis, which was triggered in response to systemic viral infection. This study highlights the prevalence and organ-specific complexity of immune gene activity in non-haematopoietic structural cells, and it provides a high-resolution, multi-omics atlas of the epigenetic and transcriptional networks that regulate structural cells in the mouse.

Ultra-high-throughput single-cell RNA sequencing and perturbation screening with combinatorial fluidic indexing.

Cell atlas projects and high-throughput perturbation screens require single-cell sequencing at a scale that is challenging with current technology. To enable cost-effective single-cell sequencing for millions of individual cells, we developed 'single-cell combinatorial fluidic indexing' (scifi). The scifi-RNA-seq assay combines one-step combinatorial preindexing of entire transcriptomes inside permeabilized cells with subsequent single-cell RNA-seq using microfluidics. Preindexing allows us to load several cells per droplet and computationally demultiplex their individual expression profiles. Thereby, scifi-RNA-seq massively increases the throughput of droplet-based single-cell RNA-seq, and provides a straightforward way of multiplexing thousands of samples in a single experiment. Compared with multiround combinatorial indexing, scifi-RNA-seq provides an easy and efficient workflow. Compared to cell hashing methods, which flag and discard droplets containing more than one cell, scifi-RNA-seq resolves and retains individual transcriptomes from overloaded droplets. We benchmarked scifi-RNA-seq on various human and mouse cell lines, validated it for primary human T cells and applied it in a highly multiplexed CRISPR screen with single-cell transcriptome readout of T cell receptor activation.

Hematopoietic expression of a chimeric murine-human CALR oncoprotein allows the assessment of anti-CALR antibody immunotherapies in vivo.

Myeloproliferative neoplasms (MPNs) are characterized by a pathologic expansion of myeloid lineages. Mutations in JAK2, CALR and MPL genes are known to be three prominent MPN disease drivers. Mutant CALR (mutCALR) is an oncoprotein that interacts with and activates the thrombopoietin receptor (MPL) and represents an attractive target for targeted therapy of CALR mutated MPN. We generated a transgenic murine model with conditional expression of the human mutant exon 9 (del52) from the murine endogenous Calr locus. These mice develop essential thrombocythemia like phenotype with marked thrombocytosis and megakaryocytosis. The disease exacerbates with age showing prominent signs of splenomegaly and anemia. The disease is transplantable and mutCALR stem cells show proliferative advantage when compared to wild type stem cells. Transcriptome profiling of hematopoietic stem cells revealed oncogenic and inflammatory gene expression signatures. To demonstrate the applicability of the transgenic animals for immunotherapy, we treated mice with monoclonal antibody raised against the human mutCALR. The antibody treatment lowered platelet and stem cell counts in mutant mice. Secretion of mutCALR did not constitute a significant antibody sink. This animal model not only recapitulates human MPN but also serves as a relevant model for testing immunotherapeutic strategies targeting epitopes of the human mutCALR.

Viral variant-resolved wastewater surveillance of SARS-CoV-2 at national scale.

SARS-CoV-2 surveillance by wastewater-based epidemiology is poised to provide a complementary approach to sequencing individual cases. However, robust quantification of variants and de novo detection of emerging variants remains challenging for existing strategies. We deep sequenced 3,413 wastewater samples representing 94 municipal catchments, covering >59% of the population of Austria, from December 2020 to February 2022. Our system of variant quantification in sewage pipeline designed for robustness (termed VaQuERo) enabled us to deduce the spatiotemporal abundance of predefined variants from complex wastewater samples. These results were validated against epidemiological records of >311,000 individual cases. Furthermore, we describe elevated viral genetic diversity during the Delta variant period, provide a framework to predict emerging variants and measure the reproductive advantage of variants of concern by calculating variant-specific reproduction numbers from wastewater. Together, this study demonstrates the power of national-scale WBE to support public health and promises particular value for countries without extensive individual monitoring.

Genomic epidemiology of superspreading events in Austria reveals mutational dynamics and transmission properties of SARS-CoV-2.

Superspreading events shaped the coronavirus disease 2019 (COVID-19) pandemic, and their rapid identification and containment are essential for disease control. Here, we provide a national-scale analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) superspreading during the first wave of infections in Austria, a country that played a major role in initial virus transmissions in Europe. Capitalizing on Austria's well-developed epidemiological surveillance system, we identified major SARS-CoV-2 clusters during the first wave of infections and performed deep whole-genome sequencing of more than 500 virus samples. Phylogenetic-epidemiological analysis enabled the reconstruction of superspreading events and charts a map of tourism-related viral spread originating from Austria in spring 2020. Moreover, we exploited epidemiologically well-defined clusters to quantify SARS-CoV-2 mutational dynamics, including the observation of low-frequency mutations that progressed to fixation within the infection chain. Time-resolved virus sequencing unveiled viral mutation dynamics within individuals with COVID-19, and epidemiologically validated infector-infectee pairs enabled us to determine an average transmission bottleneck size of 103 SARS-CoV-2 particles. In conclusion, this study illustrates the power of combining epidemiological analysis with deep viral genome sequencing to unravel the spread of SARS-CoV-2 and to gain fundamental insights into mutational dynamics and transmission properties.

The DNA methylation landscape of glioblastoma disease progression shows extensive heterogeneity in time and space.

Glioblastoma is characterized by widespread genetic and transcriptional heterogeneity, yet little is known about the role of the epigenome in glioblastoma disease progression. Here, we present genome-scale maps of DNA methylation in matched primary and recurring glioblastoma tumors, using data from a highly annotated clinical cohort that was selected through a national patient registry. We demonstrate the feasibility of DNA methylation mapping in a large set of routinely collected FFPE samples, and we validate bisulfite sequencing as a multipurpose assay that allowed us to infer a range of different genetic, epigenetic, and transcriptional characteristics of the profiled tumor samples. On the basis of these data, we identified subtle differences between primary and recurring tumors, links between DNA methylation and the tumor microenvironment, and an association of epigenetic tumor heterogeneity with patient survival. In summary, this study establishes an open resource for dissecting DNA methylation heterogeneity in a genetically diverse and heterogeneous cancer, and it demonstrates the feasibility of integrating epigenomics, radiology, and digital pathology for a national cohort, thereby leveraging existing samples and data collected as part of routine clinical practice.