Immunological Genome Project and systems immunology

X-linked Deficiency in ELF4 in Females with Skewed X Chromosome Inactivation

Deficiency in ELF4, X-linked (DEX) is a newly identified monogenic autoinflammatory disease. Most reported cases are male, leading to the recognition of DEX being primarily limited to male patients. Here we described 3 pediatric female patients with DEX from 3 unrelated families, who are all heterozygous for ELF4 mutations (c.320_c.321insA, c.329delA and c.685 A > G). Similar to reported male DEX patients, the main clinical features include recurring oral ulcers, abdominal pain and diarrhea with colonoscopy showing ulcers in the colon. Meanwhile, novel and effective treatment strategies, such as the use of the biologic vedolizumab and exclusive enteral nutrition (EEN), have provided additional options for the treatment of DEX. Finally, we observed skewed X chromosome inactivation patterns in all three female patients, with over-inactivation of the X chromosome carrying the wild-type allele confirmed in two of them.

A twin analysis to estimate genetic and environmental factors contributing to variation in weighted gene co-expression network module eigengenes

Multivariate network-based analytic methods such as weighted gene co-expression network analysis are frequently applied to human and animal gene-expression data to estimate the first principal component of a module, or module eigengene (ME). MEs are interpreted as multivariate summaries of correlated gene-expression patterns and network connectivity across genes within a module. As such, they have the potential to elucidate the mechanisms by which molecular genomic variation contributes to individual differences in complex traits. Although increasingly used to test for associations between modules and complex traits, the genetic and environmental etiology of MEs has not been empirically established. It is unclear if, and to what degree, individual differences in blood-derived MEs reflect random variation versus familial aggregation arising from heritable or shared environmental influences. We used biometrical genetic analyses to estimate the contribution of genetic and environmental influences on MEs derived from blood lymphocytes collected on a sample of N = 661 older male twins from the Vietnam Era Twin Study of Aging (VETSA) whose mean age at assessment was 67.7 years (SD = 2.6 years, range = 62-74 years). Of the 26 detected MEs, 14 (56%) had statistically significant additive genetic variation with an average heritability of 44% (SD = 0.08, range = 35%-64%). Despite the relatively small sample size, this demonstration of significant family aggregation including estimates of heritability in 14 of the 26 MEs suggests that blood-based MEs are reliable and merit further exploration in terms of their associations with complex traits and diseases.

A specific gene expression program underlies antigen archiving by lymphatic endothelial cells in mammalian lymph nodes

Lymph node (LN) lymphatic endothelial cells (LEC) actively acquire and archive foreign antigens. Here, we address questions of how LECs achieve durable antigen archiving and whether LECs with high levels of antigen express unique transcriptional programs. We used single cell sequencing in dissociated LN tissue and spatial transcriptomics to quantify antigen levels in LEC subsets and dendritic cell populations at multiple time points after immunization and determined that ceiling and floor LECs archive antigen for the longest duration. We identify, using spatial transcriptomics, antigen positive LEC-dendritic cell interactions. Using a prime-boost strategy we find increased antigen levels within LECs after a second immunization demonstrating that LEC antigen acquisition and archiving capacity can be improved over multiple exposures. Using machine learning we defined a unique transcriptional program within archiving LECs that predicted LEC archiving capacity in mouse and human independent data sets. We validated this modeling, showing we could predict lower levels of LEC antigen archiving in chikungunya virus-infected mice and demonstrated in vivo the accuracy of our prediction. Collectively, our findings establish unique properties of LECs and a defining transcriptional program for antigen archiving that can predict antigen archiving capacity in different disease states and organisms.

Suture-anchored cutaneous tension induces persistent hypertrophic scarring in a novel murine model

Background

Hypertrophic scars cause impaired skin appearance and function, seriously affecting physical and mental health. Due to medical ethics and clinical accessibility, the collection of human scar specimens is frequently restricted, and the establishment of scar experimental animal models for scientific research is urgently needed. The four most commonly used animal models of hypertrophic scars have the following drawbacks: the rabbit ear model takes a long time to construct; the immunodeficient mouse hypertrophic scar model necessitates careful feeding and experimental operations; female Duroc pigs are expensive to purchase and maintain, and their large size makes it difficult to produce a significant number of models; and mouse scar models that rely on tension require special skin stretch devices, which are often damaged and shed, resulting in unstable model establishment. Our group overcame the shortcomings of previous scar animal models and created a new mouse model of hypertrophic scarring induced by suture anchoring at the wound edge.

Methods

We utilized suture anchoring of incisional wounds to impose directional tension throughout the healing process, restrain wound contraction, and generate granulation tissue, thus inducing scar formation. Dorsal paired incisions were generated in mice, with wound edges on the upper back sutured to the rib cage and the wound edges on the lower back relaxed as a control. Macroscopic manifestation, microscopic histological analysis, mRNA sequencing, bioinformatics, and in vitro cell assays were also conducted to verify the reliability of this method.

Results

Compared with those in relaxed controls, the fibrotic changes in stretched wounds were more profound. Histologically, the stretched scars were hypercellular, hypervascular, and hyperproliferative with disorganized extracellular matrix deposition, and displayed molecular hallmarks of hypertrophic fibrosis. In addition, the stretched scars exhibited transcriptional overlap with mechanically stretched scars, and human hypertrophic and keloid scars. Phosphatidylinositol 3-kinase-serine/threonine-protein kinase B signaling was implicated as a profibrotic mediator of apoptosis resistance under suture-induced tension.

Conclusions

This straightforward murine model successfully induces cardinal molecular and histological features of pathological hypertrophic scarring through localized suture tension to inhibit wound contraction. The model enables us to interrogate the mechanisms of tension-induced fibrosis and evaluate anti-scarring therapies.

A specific and portable gene expression program underlies antigen archiving by lymphatic endothelial cells

Antigens from protein subunit vaccination traffic from the tissue to the draining lymph node, either passively via the lymph or carried by dendritic cells at the local injection site. Lymph node (LN) lymphatic endothelial cells (LEC) actively acquire and archive foreign antigens, and archived antigen can be released during subsequent inflammatory stimulus to improve immune responses. Here, we answer questions about how LECs achieve durable antigen archiving and whether there are transcriptional signatures associated with LECs containing high levels of antigen. We used single cell sequencing in dissociated LN tissue to quantify antigen levels in LEC and dendritic cell populations at multiple timepoints after immunization, and used machine learning to define a unique transcriptional program within archiving LECs that can predict LEC archiving capacity in independent data sets. Finally, we validated this modeling, showing we could predict antigen archiving from a transcriptional dataset of CHIKV infected mice and demonstrated in vivo the accuracy of our prediction. Collectively, our findings establish a unique transcriptional program in LECs that promotes antigen archiving that can be translated to other systems.

Integrative temporal multi-omics reveals uncoupling of transcriptome and proteome during human T cell activation

Engagement of the T cell receptor (TCR) triggers molecular reprogramming leading to the acquisition of specialized effector functions by CD4 helper and CD8 cytotoxic T cells. While transcription factors, chemokines, and cytokines are known drivers in this process, the temporal proteomic and transcriptomic changes that regulate different stages of human primary T cell activation remain to be elucidated. Here, we report an integrative temporal proteomic and transcriptomic analysis of primary human CD4 and CD8 T cells following ex vivo stimulation with anti-CD3/CD28 beads, which revealed major transcriptome-proteome uncoupling. The early activation phase in both CD4 and CD8 T cells was associated with transient downregulation of the mRNA transcripts and protein of the central glucose transport GLUT1. In the proliferation phase, CD4 and CD8 T cells became transcriptionally more divergent while their proteome became more similar. In addition to the kinetics of proteome-transcriptome correlation, this study unveils selective transcriptional and translational metabolic reprogramming governing CD4 and CD8 T cell responses to TCR stimulation. This temporal transcriptome/proteome map of human T cell activation provides a reference map exploitable for future discovery of biomarkers and candidates targeting T cell responses.

INHBA(+) cancer-associated fibroblasts generate an immunosuppressive tumor microenvironment in ovarian cancer

Effective targeting of cancer-associated fibroblasts (CAFs) is hindered by the lack of specific biomarkers and a poor understanding of the mechanisms by which different populations of CAFs contribute to cancer progression. While the role of TGFβ in CAFs is well-studied, less attention has been focused on a structurally and functionally similar protein, Activin A (encoded by INHBA). Here, we identified INHBA(+) CAFs as key players in tumor promotion and immunosuppression. Spatiotemporal analyses of patient-matched primary, metastatic, and recurrent ovarian carcinomas revealed that aggressive metastatic tumors enriched in INHBA(+) CAFs were also enriched in regulatory T cells (Tregs). In ovarian cancer mouse models, intraperitoneal injection of the Activin A neutralizing antibody attenuated tumor progression and infiltration with pro-tumorigenic subsets of myofibroblasts and macrophages. Downregulation of INHBA in human ovarian CAFs inhibited pro-tumorigenic CAF functions. Co-culture of human ovarian CAFs and T cells revealed the dependence of Treg differentiation on direct contact with INHBA(+) CAFs. Mechanistically, INHBA/recombinant Activin A in CAFs induced the autocrine expression of PD-L1 through SMAD2-dependent signaling, which promoted Treg differentiation. Collectively, our study identified an INHBA(+) subset of immunomodulatory pro-tumoral CAFs as a potential therapeutic target in advanced ovarian cancers which typically show a poor response to immunotherapy.

B cell receptor repertoire abnormalities in autoimmune disease

B cells play a crucial role in the immune response and contribute to various autoimmune diseases. Recent studies have revealed abnormalities in the B cell receptor (BCR) repertoire of patients with autoimmune diseases, with distinct features observed among different diseases and B cell subsets. Classically, BCR repertoire was used as an identifier of distinct antigen-specific clonotypes, but the recent advancement of analyzing large-scale repertoire has enabled us to use it as a tool for characterizing cellular biology. In this review, we provide an overview of the BCR repertoire in autoimmune diseases incorporating insights from our latest research findings. In systemic lupus erythematosus (SLE), we observed a significant skew in the usage of VDJ genes, particularly in CD27+IgD+ unswitched memory B cells and plasmablasts. Notably, autoreactive clones within unswitched memory B cells were found to be increased and strongly associated with disease activity, underscoring the clinical significance of this subset. Similarly, various abnormalities in the BCR repertoire have been reported in other autoimmune diseases such as rheumatoid arthritis. Thus, BCR repertoire analysis holds potential for enhancing our understanding of the underlying mechanisms involved in autoimmune diseases. Moreover, it has the potential to predict treatment effects and identify therapeutic targets in autoimmune diseases.

Lung dendritic-cell metabolism underlies susceptibility to viral infection in diabetes

People with diabetes feature a life-risking susceptibility to respiratory viral infection, including influenza and SARS-CoV-2 (ref. 1), whose mechanism remains unknown. In acquired and genetic mouse models of diabetes, induced with an acute pulmonary viral infection, we demonstrate that hyperglycaemia leads to impaired costimulatory molecule expression, antigen transport and T cell priming in distinct lung dendritic cell (DC) subsets, driving a defective antiviral adaptive immune response, delayed viral clearance and enhanced mortality. Mechanistically, hyperglycaemia induces an altered metabolic DC circuitry characterized by increased glucose-to-acetyl-CoA shunting and downstream histone acetylation, leading to global chromatin alterations. These, in turn, drive impaired expression of key DC effectors including central antigen presentation-related genes. Either glucose-lowering treatment or pharmacological modulation of histone acetylation rescues DC function and antiviral immunity. Collectively, we highlight a hyperglycaemia-driven metabolic-immune axis orchestrating DC dysfunction during pulmonary viral infection and identify metabolic checkpoints that may be therapeutically exploited in mitigating exacerbated disease in infected diabetics.

Endothelial cells are a key target of IFN-g during response to combined PD-1/CTLA-4 ICB treatment in a mouse model of bladder cancer

To explore mechanisms of response to combined PD-1/CTLA-4 immune checkpoint blockade (ICB) treatment in individual cell types, we generated scRNA-seq using a mouse model of invasive urothelial carcinoma with three conditions: untreated tumor, treated tumor, and tumor treated after CD4+ T cell depletion. After classifying tumor cells based on detection of somatic variants and assigning non-tumor cell types using SingleR, we performed differential expression analysis, overrepresentation analysis, and gene set enrichment analysis (GSEA) within each cell type. GSEA revealed that endothelial cells were enriched for upregulated IFN-g response genes when comparing treated cells to both untreated cells and cells treated after CD4+ T cell depletion. Functional analysis showed that knocking out IFNgR1 in endothelial cells inhibited treatment response. Together, these results indicated that IFN-g signaling in endothelial cells is a key mediator of ICB induced anti-tumor activity.

Pre-pregnancy gene expression signatures are associated with subsequent improvement/worsening of rheumatoid arthritis during pregnancy

BACKGROUND

While many women with rheumatoid arthritis (RA) improve during pregnancy and others worsen, there are no biomarkers to predict this improvement or worsening. In our unique RA pregnancy cohort that includes a pre-pregnancy baseline, we have examined pre-pregnancy gene co-expression networks to identify differences between women with RA who subsequently improve during pregnancy and those who worsen.

METHODS

Blood samples were collected before pregnancy (T0) from 19 women with RA and 13 healthy women enrolled in our prospective pregnancy cohort. RA improvement/worsening between T0 and 3rd trimester was assessed by changes in the Clinical Disease Activity Index (CDAI). Pre-pregnancy expression profiles were examined by RNA sequencing and differential gene expression analysis. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules correlated with the improvement/worsening of RA during pregnancy and to assess their functional relevance.

RESULTS

Of the 19 women with RA, 14 improved during pregnancy (RAimproved) while 5 worsened (RAworsened). At the T0 baseline, however, the mean CDAI was similar between the two groups. WGCNA identified one co-expression module related to B cell function that was significantly correlated with the worsening of RA during pregnancy and was significantly enriched in genes differentially expressed between the RAimproved and RAworsened groups. A neutrophil-related expression signature was also identified in the RAimproved group at the T0 baseline.

CONCLUSION

The pre-pregnancy gene expression signatures identified represent potential biomarkers to predict the subsequent improvement/worsening of RA during pregnancy, which has important implications for the personalized treatment of RA during pregnancy.

Systems Immunology Approaches to Metabolism

Over the last decade, immunometabolism has emerged as a novel interdisciplinary field of research and yielded significant fundamental insights into the regulation of immune responses. Multiple classical approaches to interrogate immunometabolism, including bulk metabolic profiling and analysis of metabolic regulator expression, paved the way to appreciating the physiological complexity of immunometabolic regulation in vivo. Studying immunometabolism at the systems level raised the need to transition towards the next-generation technology for metabolic profiling and analysis. Spatially resolved metabolic imaging and computational algorithms for multi-modal data integration are new approaches to connecting metabolism and immunity. In this review, we discuss recent studies that highlight the complex physiological interplay between immune responses and metabolism and give an overview of technological developments that bear the promise of capturing this complexity most directly and comprehensively.

Endothelial cells are a key target of IFN-g during response to combined PD-1/CTLA-4 ICB treatment in a mouse model of bladder cancer

To explore mechanisms of response to combined PD-1/CTLA-4 immune checkpoint blockade (ICB) treatment in individual cell types, we generated scRNA-seq using a mouse model of invasive urothelial carcinoma with three conditions: untreated tumor, treated tumor, and tumor treated after CD4+ T cell depletion. After classifying tumor cells based on detection of somatic variants and assigning non-tumor cell types using SingleR, we performed differential expression analysis, overrepresentation analysis, and gene set enrichment analysis (GSEA) within each cell type. GSEA revealed that endothelial cells were enriched for upregulated IFN-g response genes when comparing treated cells to both untreated cells and cells treated after CD4+ T cell depletion. Functional analysis showed that knocking out IFNgR1 in endothelial cells inhibited treatment response. Together, these results indicated that IFN-g signaling in endothelial cells is a key mediator of ICB induced anti-tumor activity.

A network view of human immune system and virus-human interaction

The immune system is highly networked and complex, which is continuously changing as encountering old and new pathogens. However, reductionism-based researches do not give a systematic understanding of the molecular mechanism of the immune response and viral pathogenesis. Here, we present HUMPPI-2022, a high-quality human protein-protein interaction (PPI) network, containing > 11,000 protein-coding genes with > 78,000 interactions. The network topology and functional characteristics analyses of the immune-related genes (IRGs) reveal that IRGs are mostly located in the center of the network and link genes of diverse biological processes, which may reflect the gene pleiotropy phenomenon. Moreover, the virus-human interactions reveal that pan-viral targets are mostly hubs, located in the center of the network and enriched in fundamental biological processes, but not for coronavirus. Finally, gene age effect was analyzed from the view of the host network for IRGs and virally-targeted genes (VTGs) during evolution, with IRGs gradually became hubs and integrated into host network through bridging functionally differentiated modules. Briefly, HUMPPI-2022 serves as a valuable resource for gaining a better understanding of the composition and evolution of human immune system, as well as the pathogenesis of viruses.

Identification and characterization of in vitro expanded hematopoietic stem cells

Hematopoietic stem cells (HSCs) cultured outside the body are the fundamental component of a wide range of cellular and gene therapies. Recent efforts have achieved > 200-fold expansion of functional HSCs, but their molecular characterization has not been possible since the majority of cells are non-HSCs and single cell-initiated cultures have substantial clone-to-clone variability. Using the Fgd5 reporter mouse in combination with the EPCR surface marker, we report exclusive identification of HSCs from non-HSCs in expansion cultures. By directly linking single-clone functional transplantation data with single-clone gene expression profiling, we show that the molecular profile of expanded HSCs is similar to proliferating fetal HSCs and reveals a gene expression signature, including Esam, Prdm16, Fstl1, and Palld, that can identify functional HSCs from multiple cellular states. This "repopulation signature" (RepopSig) also enriches for HSCs in human datasets. Together, these findings demonstrate the power of integrating functional and molecular datasets to better derive meaningful gene signatures and opens the opportunity for a wide range of functional screening and molecular experiments previously not possible due to limited HSC numbers.

Current understanding of the human microbiome in glioma

There is mounting evidence that the human microbiome is highly associated with a wide variety of central nervous system diseases. However, the link between the human microbiome and glioma is rarely noticed. The exact mechanism of microbiota to affect glioma remains unclear. Recent studies have demonstrated that the microbiome may affect the development, progress, and therapy of gliomas, including the direct impacts of the intratumoral microbiome and its metabolites, and the indirect effects of the gut microbiome and its metabolites. Glioma-related microbiome (gut microbiome and intratumoral microbiome) is associated with both tumor microenvironment and tumor immune microenvironment, which ultimately influence tumorigenesis, progression, and responses to treatment. In this review, we briefly summarize current knowledge regarding the role of the glioma-related microbiome, focusing on its gut microbiome fraction and a brief description of the intratumoral microbiome, and put forward the prospects in which microbiome can be applied in the future and some challenges still need to be solved.

Primary Infection by E. multilocularis Induces Distinct Patterns of Cross Talk between Hepatic Natural Killer T Cells and Regulatory T Cells in Mice

The larval stage of the helminthic cestode Echinococcus multilocularis can inflict tumor-like hepatic lesions that cause the parasitic disease alveolar echinococcosis in humans, with high mortality in untreated patients. Opportunistic properties of the disease have been established based on the increased incidence in immunocompromised patients and mouse models, indicating that an appropriate adaptive immune response is required for the control of the disease. However, cellular interactions and the kinetics of the local hepatic immune responses during the different stages of infection with E. multilocularis remain unknown. In a mouse model of oral infection that mimics the normal infection route in human patients, the networks of the hepatic immune response were assessed using single-cell RNA sequencing (scRNA-seq) of isolated hepatic CD3⁺ T cells at different infection stages. We observed an early and sustained significant increase in natural killer T (NKT) cells and regulatory T cells (Tregs). Early tumor necrosis factor (TNF)- and integrin-dependent interactions between these two cell types promote the formation of hepatic lesions. At late time points, downregulation of programmed cell death protein 1 (PD-1) and ectonucleoside triphosphate diphosphohydrolase 1 (ENTPD1)-dependent signaling suppress the resolution of parasite-induced pathology. The obtained data provide fresh insight into the adaptive immune responses and local regulatory pathways at different infection stages of E. multilocularis in mice.

Single-cell dual-omics reveals the transcriptomic and epigenomic diversity of cardiac non-myocytes

AIMS

The precise cellular identity and molecular features of non-myocytes (non-CMs) in a mammalian heart at a single-cell level remain elusive. Depiction of epigenetic landscape with transcriptomic signatures using the latest single-cell multi-omics has the potential to unravel the molecular programs underlying the cellular diversity of cardiac non-myocytes. Here, we characterized the molecular and cellular features of cardiac non-CM populations in the adult murine heart at the single-cell level.

METHODS AND RESULTS

Through single-cell dual omics analysis, we mapped the epigenetic landscapes, characterized the transcriptomic profiles and delineated the molecular signatures of cardiac non-CMs in the adult murine heart. Distinct cis-regulatory elements and trans-acting factors for the individual major non-CM cell types (endothelial cells, fibroblast, pericytes, and immune cells) were identified. In particular, unbiased sub-clustering and functional annotation of cardiac fibroblasts (FBs) revealed extensive FB heterogeneity and identified FB sub-types with functional states related to the cellular response to stimuli, cytoskeleton organization, and immune regulation, respectively. We further explored the function of marker genes Hsd11b1 and Gfpt2 that label major FB subpopulations and determined the distribution of Hsd11b1+ and Gfp2+ FBs in murine healthy and diseased hearts.

CONCLUSIONS

In summary, we characterized the non-CM cellular identity at the transcriptome and epigenome levels using single-cell omics approaches and discovered previously unrecognized cardiac fibroblast subpopulations with unique functional states.

Dissecting the molecular control of immune cell accumulation in the inflamed joint

Mechanisms governing entry and exit of immune cells into and out of inflamed joints remain poorly understood. We sought herein to identify the key molecular pathways regulating such migration. Using murine models of inflammation in conjunction with mice expressing a photoconvertible fluorescent protein, we characterized the migration of cells from joints to draining lymph nodes and performed RNA-Seq analysis on isolated cells, identifying genes associated with migration and retention. We further refined the gene list to those specific for joint inflammation. RNA-Seq data revealed pathways and genes previously highlighted as characteristic of rheumatoid arthritis in patient studies, validating the methodology. Focusing on pathways associated with cell migration, adhesion, and movement, we identified genes involved in the retention of immune cells in the inflamed joint, namely junctional adhesion molecule A (JAM-A), and identified a role for such molecules in T cell differentiation in vivo. Thus, using a combination of cell-tracking approaches and murine models of inflammatory arthritis, we identified genes, pathways, and anatomically specific tissue signatures regulating cell migration in a variety of inflamed sites. This skin- and joint-specific data set will be an invaluable resource for the identification of therapeutic targets for arthritis and other inflammatory disorders.

Bioinformatic Challenges Detecting Genetic Variation in Precision Medicine Programs

Precision medicine programs to identify clinically relevant genetic variation have been revolutionized by access to increasingly affordable high-throughput sequencing technologies. A decade of continual drops in per-base sequencing costs means it is now feasible to sequence an individual patient genome and interrogate all classes of genetic variation for < $1,000 USD. However, while advances in these technologies have greatly simplified the ability to obtain patient sequence information, the timely analysis and interpretation of variant information remains a challenge for the rollout of large-scale precision medicine programs. This review will examine the challenges and potential solutions that exist in identifying predictive genetic biomarkers and pharmacogenetic variants in a patient and discuss the larger bioinformatic challenges likely to emerge in the future. It will examine how both software and hardware development are aiming to overcome issues in short read mapping, variant detection and variant interpretation. It will discuss the current state of the art for genetic disease and the remaining challenges to overcome for complex disease. Success across all types of disease will require novel statistical models and software in order to ensure precision medicine programs realize their full potential now and into the future.

ImmuneData: an integrated data discovery system for immunology data repositories

To meet the increasing demand for data sharing, data reuse and meta-analysis in the immunology research community, we have developed the data discovery system ImmuneData. The system provides integrated access to five immunology data repositories funded by the National Institute of Allergy and Infectious Diseases, Division of Allergy, Immunology and Transplantation, including ImmPort, ImmuneSpace, ITN TrialShare, ImmGen and IEDB. ImmuneData restructures the data repositories' metadata into a uniform schema using domain experts' knowledge and state-of-the-art Natural Language Processing (NLP) technologies. It comes with a user-friendly web interface, accessible at http://www.immunedata.org/, and a Google-like search engine for biological researchers to find and access data easily. The vast quantity of synonyms used in biomedical research increase the likelihood of incomplete search results. Thus, our search engine converts queries submitted by users into ontology terms, which are then expended by NLP technologies to ensure that the search results will include all synonyms for a particular concept. The system also includes an advanced search function to build customized queries to meet higher-level users' needs. ImmuneData ensures the FAIR principle (Findability, Accessibility, Interoperability and Reusability) of the five data repositories to benefit data reuse in the immunology research community. The data pipeline constructing our system can be extended to other data repositories to build a more comprehensive biological data discovery system.

DATABASE URL

http://www.immunedata.org/.

Single-cell analysis of skin immune cells reveals an Angptl4-ifi20b axis that regulates monocyte differentiation during wound healing

The persistent inflammatory response at the wound site is a cardinal feature of nonhealing wounds. Prolonged neutrophil presence in the wound site due to failed clearance by reduced monocyte-derived macrophages delays the transition from the inflammatory to the proliferative phase of wound healing. Angiopoietin-like 4 protein (Angptl4) is a matricellular protein that has been implicated in many inflammatory diseases. However, its precise role in the immune cell response during wound healing remains unclear. Therefore, we performed flow cytometry and single-cell RNA sequencing to examine the immune cell landscape of excisional wounds from Angptl4+/+ and Angptl4−/− mice. Chemotactic immune cell recruitment and infiltration were not compromised due to Angptl4 deficiency. However, as wound healing progresses, Angptl4−/− wounds have a prolonged neutrophil presence and fewer monocyte-derived macrophages than Angptl4+/+ and Angptl4LysM−/− wounds. The underlying mechanism involves a novel Angptl4-interferon activated gene 202B (ifi202b) axis that regulates monocyte differentiation to macrophages, coordinating neutrophil removal and inflammation resolution. An unbiased kinase inhibitor screen revealed an Angptl4-mediated kinome signaling network involving S6K, JAK, and CDK, among others, that modulates the expression of ifi202b. Silencing ifi202b in Angptl4−/− monocytes, whose endogenous expression was elevated, rescued the impaired monocyte-to-macrophage transition in the in vitro reconstituted wound microenvironment using wound exudate. GSEA and IPA functional analyses revealed that ifi202b-associated canonical pathways and functions involved in the inflammatory response and monocyte cell fate were enriched. Together, we identified ifi202b as a key gatekeeper of monocyte differentiation. By modulating ifi202b expression, Angptl4 orchestrates the inflammatory state, innate immune landscape, and wound healing process.

Imaging dendritic cell functions

Dendritic cells (DCs) are crucial for the appropriate initiation of adaptive immune responses. During inflammation, DCs capture antigens, mature, and migrate to lymphoid tissues to present foreign material to naïve T cells. These cells get activated and differentiate either into pathogen-specific cytotoxic CD8⁺ T cells that destroy infected cells or into CD4⁺ T helper cells that, among other effector functions, orchestrate antibody production by B cells. DC-mediated antigen presentation is equally important in non-inflammatory conditions. Here, DCs mediate induction of tolerance by presenting self-antigens or harmless environmental antigens and induce differentiation of regulatory T cells or inactivation of self-reactive immune cells. Detailed insights into the biology of DCs are, therefore, crucial for the development of novel vaccines as well as the prevention of autoimmune diseases. As in many other life science areas, our understanding of DC biology would be extremely restricted without bioimaging, a compilation of methods that visualize biological processes. Spatiotemporal tracking of DCs relies on various imaging tools, which not only enable insights into their positioning and migration within tissues or entire organs but also allow visualization of subcellular and molecular processes. This review aims to provide an overview of the imaging toolbox and to provide examples of diverse imaging techniques used to obtain fundamental insights into DC biology.

Dynamic prostate cancer transcriptome analysis delineates the trajectory to disease progression

Comprehensive genomic studies have delineated key driver mutations linked to disease progression for most cancers. However, corresponding transcriptional changes remain largely elusive because of the bias associated with cross-study analysis. Here, we overcome these hurdles and generate a comprehensive prostate cancer transcriptome atlas that describes the roadmap to tumor progression in a qualitative and quantitative manner. Most cancers follow a uniform trajectory characterized by upregulation of polycomb-repressive-complex-2, G2-M checkpoints, and M2 macrophage polarization. Using patient-derived xenograft models, we functionally validate our observations and add single-cell resolution. Thereby, we show that tumor progression occurs through transcriptional adaption rather than a selection of pre-existing cancer cell clusters. Moreover, we determine at the single-cell level how inhibition of EZH2 - the top upregulated gene along the trajectory - reverts tumor progression and macrophage polarization. Finally, a user-friendly web-resource is provided enabling the investigation of dynamic transcriptional perturbations linked to disease progression.

Epstein-Barr virus nuclear antigen 2 extensively rewires the human chromatin landscape at autoimmune risk loci

The interplay between environmental and genetic factors plays a key role in the development of many autoimmune diseases. In particular, the Epstein-Barr virus (EBV) is an established contributor to multiple sclerosis, lupus, and other disorders. Previously, we showed that the EBV nuclear antigen 2 (EBNA2) transactivating protein occupies up to half of the risk loci for a set of seven autoimmune disorders. To further examine the mechanistic roles played by EBNA2 at these loci on a genome-wide scale, we globally examined gene expression, chromatin accessibility, chromatin looping, and EBNA2 binding in a B cell line that was (1) uninfected, (2) infected with a strain of EBV lacking EBNA2, or (3) infected with a strain that expresses EBNA2. We identified more than 400 EBNA2-dependent differentially expressed human genes and more than 5000 EBNA2 binding events in the human genome. ATAC-seq analysis revealed more than 2000 regions in the human genome with EBNA2-dependent chromatin accessibility, and HiChIP data revealed more than 1700 regions where EBNA2 altered chromatin looping interactions. Autoimmune genetic risk loci were highly enriched at the sites of these EBNA2-dependent chromatin-altering events. We present examples of autoimmune risk genotype-dependent EBNA2 events, nominating genetic risk mechanisms for autoimmune risk loci such as ZMIZ1 Taken together, our results reveal important interactions between host genetic variation and EBNA2-driven disease mechanisms. Further, our study highlights a critical role for EBNA2 in rewiring human gene regulatory programs through rearrangement of the chromatin landscape and nominates these interactions as components of genetic mechanisms that influence the risk of multiple autoimmune diseases.

Tree-Based Co-Clustering Identifies Chromatin Accessibility Patterns Associated With Hematopoietic Lineage Structure

Chromatin accessibility, as measured by ATACseq, varies between hematopoietic cell types in different lineages of the hematopoietic differentiation tree, e.g. T cells vs. B cells, but methods that associate variation in chromatin accessibility to the lineage structure of the differentiation tree are lacking. Using an ATACseq dataset recently published by the ImmGen consortium, we construct associations between chromatin accessibility and hematopoietic cell types using a novel co-clustering approach that accounts for the structure of the hematopoietic, differentiation tree. Under a model in which all loci and cell types within a co-cluster have a shared accessibility state, we show that roughly 80% of cell type associated accessibility variation can be captured through 12 cell type clusters and 20 genomic locus clusters, with the cell type clusters reflecting coherent components of the differentiation tree. Using publicly available ChIPseq datasets, we show that our clustering reflects transcription factor binding patterns with implications for regulation across cell types. We show that traditional methods such as hierarchical and kmeans clusterings lead to cell type clusters that are more dispersed on the tree than our tree-based algorithm. We provide a python package, chromcocluster, that implements the algorithms presented.

Regulation of dendritic cell function by Fc-γ-receptors and the neonatal Fc receptor

Dendritic cells (DCs) express receptors to sense pathogens and/or tissue damage and to communicate with other immune cells. Among those receptors, Fc receptors (FcRs) are triggered by the Fc region of antibodies produced during adaptive immunity. In this review, the role of FcγR and neonatal Fc receptor (FcRn) in DC immunity will be discussed. Their expression in DC subsets and impact on antigen uptake and presentation, DC maturation and polarisation of T cell responses will be described. Lastly, we will discuss the importance of FcR-mediated DC function in the context of immunity during viral infection, inflammatory disease, cancer and immunotherapy.

Electroimmunology and cardiac arrhythmia

Conduction disorders and arrhythmias remain difficult to treat and are increasingly prevalent owing to the increasing age and body mass of the general population, because both are risk factors for arrhythmia. Many of the underlying conditions that give rise to arrhythmia - including atrial fibrillation and ventricular arrhythmia, which frequently occur in patients with acute myocardial ischaemia or heart failure - can have an inflammatory component. In the past, inflammation was viewed mostly as an epiphenomenon associated with arrhythmia; however, the recently discovered inflammatory and non-canonical functions of cardiac immune cells indicate that leukocytes can be arrhythmogenic either by altering tissue composition or by interacting with cardiomyocytes; for example, by changing their phenotype or perhaps even by directly interfering with conduction. In this Review, we discuss the electrophysiological properties of leukocytes and how these cells relate to conduction in the heart. Given the thematic parallels, we also summarize the interactions between immune cells and neural systems that influence information transfer, extrapolating findings from the field of neuroscience to the heart and defining common themes. We aim to bridge the knowledge gap between electrophysiology and immunology, to promote conceptual connections between these two fields and to explore promising opportunities for future research.

Defective Epstein-Barr Virus Genomes and Atypical Viral Gene Expression in B-Cell Lines Derived from Multiple Myeloma Patients

Epstein-Barr virus (EBV) is a human gammaherpesvirus that is causally associated with various lymphomas and carcinomas. Although EBV is not typically associated with multiple myeloma (MM), it can be found in some B-cell lines derived from MM patients. Here, we analyzed two EBV-positive MM-patient-derived cell lines, IM9 and ARH77, and found defective viral genomes and atypical viral gene expression patterns. We performed transcriptome sequencing to characterize the viral and cellular properties of the two EBV-positive cell lines, compared to the canonical MM cell line 8226. Principal-component analyses indicated that IM9 and ARH77 clustered together and distinct from 8226. Immunological Genome Project analysis designated these cells as stem cell and bone marrow derived. IM9 and ARH77 displayed atypical viral gene expression, including leaky lytic cycle gene expression with an absence of lytic DNA amplification. Genome sequencing revealed that the EBV genomes in ARH77 contain large deletions, while IM9 has copy number losses in multiple EBV loci. Both IM9 and ARH77 showed EBV genome heterogeneity, suggesting cells harboring multiple and variant viral genomes. We identified atypical high-level expression of lytic genes BLRF1 and BLRF2. We demonstrated that short hairpin RNA (shRNA) depletion of BLRF2 altered viral and host gene expression, including a reduction in lytic gene activation and DNA amplification. These findings demonstrate that aberrant viral genomes and lytic gene expression persist in rare B cells derived from MM tumors, and they suggest that EBV may contribute to the etiology of MM. IMPORTANCE EBV is an oncogenic herpesvirus, but its mechanisms of oncogenesis are not fully understood. A role for EBV in MM has not yet been established. We analyzed EBV-positive B-cell lines derived from MM patients and found that the cells harbored defective viral genomes with aberrant viral gene expression patterns and cell gene signatures for bone marrow-derived lymphoid stem cells. These findings suggest that aberrant EBV latent infection may contribute to the etiology of MM.

Tryptophan metabolism drives dynamic immunosuppressive myeloid states in IDH-mutant gliomas

The dynamics and phenotypes of intratumoral myeloid cells during tumor progression are poorly understood. Here we define myeloid cellular states in gliomas by longitudinal single-cell profiling and demonstrate their strict control by the tumor genotype: in isocitrate dehydrogenase (IDH)-mutant tumors, differentiation of infiltrating myeloid cells is blocked, resulting in an immature phenotype. In late-stage gliomas, monocyte-derived macrophages drive tolerogenic alignment of the microenvironment, thus preventing T cell response. We define the IDH-dependent tumor education of infiltrating macrophages to be causally related to a complex re-orchestration of tryptophan metabolism, resulting in activation of the aryl hydrocarbon receptor. We further show that the altered metabolism of IDH-mutant gliomas maintains this axis in bystander cells and that pharmacological inhibition of tryptophan metabolism can reverse immunosuppression. In conclusion, we provide evidence of a glioma genotype-dependent intratumoral network of resident and recruited myeloid cells and identify tryptophan metabolism as a target for immunotherapy of IDH-mutant tumors.

Leukocyte dynamics after intracerebral hemorrhage in a living patient reveal rapid adaptations to tissue milieu

Intracerebral hemorrhage (ICH) is a devastating form of stroke with a high mortality rate and few treatment options. Discovery of therapeutic interventions has been slow given the challenges associated with studying acute injury in the human brain. Inflammation induced by exposure of brain tissue to blood appears to be a major part of brain tissue injury. Here, we longitudinally profiled blood and cerebral hematoma effluent from a patient enrolled in the Minimally Invasive Surgery with Thrombolysis in Intracerebral Hemorrhage Evacuation trial, offering a rare window into the local and systemic immune responses to acute brain injury. Using single-cell RNA-Seq (scRNA-Seq), this is the first report to our knowledge that characterized the local cellular response during ICH in the brain of a living patient at single-cell resolution. Our analysis revealed shifts in the activation states of myeloid and T cells in the brain over time, suggesting that leukocyte responses are dynamically reshaped by the hematoma microenvironment. Interestingly, the patient had an asymptomatic rebleed that our transcriptional data indicated occurred prior to detection by CT scan. This case highlights the rapid immune dynamics in the brain after ICH and suggests that sensitive methods such as scRNA-Seq would enable greater understanding of complex intracerebral events.

Immune repertoire: Revealing the "real-time" adaptive immune response in autoimmune diseases

The diversity of the immune repertoire (IR) enables the human immune system to distinguish multifarious antigens (Ags) that humans may encounter throughout life. At the same time, bias or abnormalities in the IR also pay a contribution to the pathogenesis of autoimmune diseases. Rapid advancements in high-throughput sequencing (HTS) technology have ushered in a new era of immune studies, revealing novel molecules and pathways that might result in autoimmunity. In the field of IR, HTS can monitor the immune response status and identify disease-specific immune repertoires. In this review, we summarize updated progress on the mechanisms of the IR and current related studies on four autoimmune diseases, particularly focusing on systemic lupus erythematosus (SLE). These autoimmune diseases can exhibit slightly or significantly skewed IRs and provide novel insights that inform our comprehending of disease pathogenesis and provide potential targets for diagnosis and treatment.

Steroid receptor coactivator 3 (SRC-3/AIB1) is enriched and functional in mouse and human Tregs

A subset of CD4 + lymphocytes, regulatory T cells (Tregs), are necessary for central tolerance and function as suppressors of autoimmunity against self-antigens. The SRC-3 coactivator is an oncogene in multiple cancers and is capable of potentiating numerous transcription factors in a wide variety of cell types. Src-3 knockout mice display broad lymphoproliferation and hypersensitivity to systemic inflammation. Using publicly available bioinformatics data and directed cellular approaches, we show that SRC-3 also is highly enriched in Tregs in mice and humans. Human Tregs lose phenotypic characteristics when SRC-3 is depleted or pharmacologically inhibited, including failure of induction from resting T cells and loss of the ability to suppress proliferation of stimulated T cells. These data support a model for SRC-3 as a coactivator that actively participates in protection from autoimmunity and may support immune evasion of cancers by contributing to the biology of Tregs.

The Role of APOBECs in Viral Replication

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) proteins are a diverse and evolutionarily conserved family of cytidine deaminases that provide a variety of functions from tissue-specific gene expression and immunoglobulin diversity to control of viruses and retrotransposons. APOBEC family expansion has been documented among mammalian species, suggesting a powerful selection for their activity. Enzymes with a duplicated zinc-binding domain often have catalytically active and inactive domains, yet both have antiviral function. Although APOBEC antiviral function was discovered through hypermutation of HIV-1 genomes lacking an active Vif protein, much evidence indicates that APOBECs also inhibit virus replication through mechanisms other than mutagenesis. Multiple steps of the viral replication cycle may be affected, although nucleic acid replication is a primary target. Packaging of APOBECs into virions was first noted with HIV-1, yet is not a prerequisite for viral inhibition. APOBEC antagonism may occur in viral producer and recipient cells. Signatures of APOBEC activity include G-to-A and C-to-T mutations in a particular sequence context. The importance of APOBEC activity for viral inhibition is reflected in the identification of numerous viral factors, including HIV-1 Vif, which are dedicated to antagonism of these deaminases. Such viral antagonists often are only partially successful, leading to APOBEC selection for viral variants that enhance replication or avoid immune elimination.

Omental macrophages secrete chemokine ligands that promote ovarian cancer colonization of the omentum via CCR1

The omentum is the most common site of ovarian cancer metastasis. Immune cell clusters called milky spots are found throughout the omentum. It is however unknown if these immune cells contribute to ovarian cancer metastasis. Here we report that omental macrophages promote the migration and colonization of ovarian cancer cells to the omentum through the secretion of chemokine ligands that interact with chemokine receptor 1 (CCR1). We found that depletion of macrophages reduces ovarian cancer colonization of the omentum. RNA-sequencing of macrophages isolated from mouse omentum and mesenteric adipose tissue revealed a specific enrichment of chemokine ligand CCL6 in omental macrophages. CCL6 and the human homolog CCL23 were both necessary and sufficient to promote ovarian cancer migration by activating ERK1/2 and PI3K pathways. Importantly, inhibition of CCR1 reduced ovarian cancer colonization. These findings demonstrate a critical mechanism of omental macrophage induced colonization by ovarian cancer cells via CCR1 signaling.

Krishnan et al. find that CCR1 ligands CCL6 and CCL23 secreted by murine and human macrophages, respectively, enhance metastatic colonization of ovarian cancer cells to the omentum in manner dependent on chemokine receptor 1 (CCR1). This study suggests that targeting CCR1 or CCL23 in ovarian cancer may be a therapeutic strategy.

Orphan Nuclear Receptor NR2F6 Suppresses T Follicular Helper Cell Accumulation through Regulation of IL-21

CD4 T follicular helper (Tfh) cells are specialized in helping B cells during the germinal center (GC) reaction and ultimately promote long-term humoral immunity. Here we report that loss of the nuclear orphan receptor NR2F6 causes enhanced survival and accumulation of Tfh cells, GC B cells, and plasma cells (PCs) following T cell-dependent immunization. Nr2f6-deficient CD4 T cell dysfunction is the primary cause of cell accumulation. Cytokine expression in Nr2f6-deficient Tfh cells is dysregulated, and Il21 expression is enhanced. Mechanistically, NR2F6 binds directly to the interleukin 21 (IL-21) promoter and a conserved noncoding sequence (CNS) near the Il21 gene in resting CD4⁺ T cells. During Tfh cell differentiation, this direct NR2F6 DNA interaction is abolished. Enhanced Tfh cell accumulation in Nr2f6-deficient mice can be reverted by blocking IL-21R signaling. Thus, NR2F6 is a critical negative regulator of IL-21 cytokine production in Tfh cells and prevents excessive Tfh cell accumulation.

•

Loss of NR2F6 results in accumulation of Tfh, GC B, and plasma cells after immunization

•

Increased GC populations depend on Nr2f6 loss within the CD4 compartment

•

NR2F6 directly binds to several sites within the Il21 promoter and CNS −36

•

NR2F6 restrains Il21 expression in CD4 cells; IL-21R blockade reduces Tfh accumulation

Network Approaches for Dissecting the Immune System

The immune system is a complex biological network composed of hierarchically organized genes, proteins, and cellular components that combat external pathogens and monitor the onset of internal disease. To meet and ultimately defeat these challenges, the immune system orchestrates an exquisitely complex interplay of numerous cells, often with highly specialized functions, in a tissue-specific manner. One of the major methodologies of systems immunology is to measure quantitatively the components and interaction levels in the immunologic networks to construct a computational network and predict the response of the components to perturbations. The recent advances in high-throughput sequencing techniques have provided us with a powerful approach to dissecting the complexity of the immune system. Here we summarize the latest progress in integrating omics data and network approaches to construct networks and to infer the underlying signaling and transcriptional landscape, as well as cell-cell communication, in the immune system, with a focus on hematopoiesis, adaptive immunity, and tumor immunology. Understanding the network regulation of immune cells has provided new insights into immune homeostasis and disease, with important therapeutic implications for inflammation, cancer, and other immune-mediated disorders.

Genetic and genomic alterations differentially dictate low-grade glioma growth through cancer stem cell-specific chemokine recruitment of T cells and microglia

BACKGROUND

One of the clinical hallmarks of low-grade gliomas (LGGs) arising in children with the neurofibromatosis type 1 (NF1) cancer predisposition syndrome is significant clinical variability with respect to tumor growth, associated neurologic deficits, and response to therapy. Numerous factors could contribute to this clinical heterogeneity, including the tumor cell of origin, the specific germline NF1 gene mutation, and the coexistence of additional genomic alterations. Since human specimens are rarely acquired, and have proven difficult to maintain in vitro or as xenografts in vivo, we have developed a series of Nf1 mutant optic glioma mouse strains representing each of these contributing factors.

METHODS

Optic glioma stem cells (o-GSCs) were generated from this collection of Nf1 genetically engineered mice, and analyzed for their intrinsic growth properties, as well as the production of chemokines that could differentially attract T cells and microglia.

RESULTS

The observed differences in Nf1 optic glioma growth are not the result of cell autonomous growth properties of o-GSCs, but rather the unique patterns of o-GSC chemokine expression, which differentially attract T cells and microglia. This immune profile collectively dictates the levels of chemokine C-C ligand 5 (Ccl5) expression, the key stromal factor that drives murine Nf1 optic glioma growth.

CONCLUSIONS

These findings reveal that genetic and genomic alterations create murine LGG biological heterogeneity through the differential recruitment of T cells and microglia by o-GSC-produced chemokines, which ultimately determine the expression of stromal factors that drive tumor growth.

Systems Biology Approaches to Understanding the Human Immune System

Systems biology is an approach to interrogate complex biological systems through large-scale quantification of numerous biomolecules. The immune system involves >1,500 genes/proteins in many interconnected pathways and processes, and a systems-level approach is critical in broadening our understanding of the immune response to vaccination. Changes in molecular pathways can be detected using high-throughput omics datasets (e.g., transcriptomics, proteomics, and metabolomics) by using methods such as pathway enrichment, network analysis, machine learning, etc. Importantly, integration of multiple omic datasets is becoming key to revealing novel biological insights. In this perspective article, we highlight the use of protein-protein interaction (PPI) networks as a multi-omics integration approach to unravel information flow and mechanisms during complex biological events, with a focus on the immune system. This involves a combination of tools, including: InnateDB, a database of curated interactions between genes and protein products involved in the innate immunity; NetworkAnalyst, a visualization and analysis platform for InnateDB interactions; and MetaBridge, a tool to integrate metabolite data into PPI networks. The application of these systems techniques is demonstrated for a variety of biological questions, including: the developmental trajectory of neonates during the first week of life, mechanisms in host-pathogen interaction, disease prognosis, biomarker discovery, and drug discovery and repurposing. Overall, systems biology analyses of omics data have been applied to a variety of immunology-related questions, and here we demonstrate the numerous ways in which PPI network analysis can be a powerful tool in contributing to our understanding of the immune system and the study of vaccines.

Detecting pathogenic variants in autoimmune diseases using high-throughput sequencing

Sequencing the first human genome in 2003 took 15 years and cost $2.7 billion. Advances in sequencing technologies have since decreased costs to the point where it is now feasible to resequence a whole human genome for $1000 in a single day. These advances have allowed the generation of huge volumes of high‐quality human sequence data used to construct increasingly large catalogs of both population‐level and disease‐causing variation. The existence of such databases, coupled with a high‐quality human reference genome, means we are able to interrogate and annotate all types of genetic variation and identify pathogenic variants for many diseases. Increasingly, sequencing‐based approaches are being used to elucidate the underlying genetic cause of autoimmune diseases, a group of roughly 80 polygenic diseases characterized by abnormal immune responses where healthy tissue is attacked. Although sequence data generation has become routine and affordable, significant challenges remain with no gold‐standard methodology to identify pathogenic variants currently available. This review examines the latest methodologies used to identify pathogenic variants in autoimmune diseases and considers available sequencing options and subsequent bioinformatic methodologies and strategies. The development of reliable and robust sequencing and analytic workflows to detect pathogenic variants is critical to realize the potential of precision medicine programs where patient variant information is used to inform clinical practice.

Dissecting the Tumor-Immune Landscape in Chimeric Antigen Receptor T-cell Therapy: Key Challenges and Opportunities for a Systems Immunology Approach

The adoptive transfer of genetically engineered chimeric antigen receptor (CAR) T cells has opened a new frontier in cancer therapy. Unlike the paradigm of targeted therapies, the efficacy of CAR T-cell therapy depends not only on the choice of target but also on a complex interplay of tumor, immune, and stromal cell communication. This presents both challenges and opportunities from a discovery standpoint. Whereas cancer consortia have traditionally focused on the genomic, transcriptomic, epigenomic, and proteomic landscape of cancer cells, there is an increasing need to expand studies to analyze the interactions between tumor, immune, and stromal cell populations in their relevant anatomical and functional compartments. Here, we focus on the promising application of systems biology to address key challenges in CAR T-cell therapy, from understanding the mechanisms of therapeutic resistance in hematologic and solid tumors to addressing important clinical challenges in biomarker discovery and therapeutic toxicity. We propose a systems biology view of key clinical objectives in CAR T-cell therapy and suggest a path forward for a biomedical discovery process that leverages modern technological approaches in systems biology.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

A transcriptome-based approach to identify functional modules within and across primary human immune cells

Genome-wide transcriptomic analyses have provided valuable insight into fundamental biology and disease pathophysiology. Many studies have taken advantage of the correlation in the expression patterns of the transcriptome to infer a potential biologic function of uncharacterized genes, and multiple groups have examined the relationship between co-expression, co-regulation, and gene function on a broader scale. Given the unique characteristics of immune cells circulating in the blood, we were interested in determining whether it was possible to identify functional co-expression modules in human immune cells. Specifically, we sequenced the transcriptome of nine immune cell types from peripheral blood cells of healthy donors and, using a combination of global and targeted analyses of genes within co-expression modules, we were able to determine functions for these modules that were cell lineage-specific or shared among multiple cell lineages. In addition, our analyses identified transcription factors likely important for immune cell lineage commitment and/or maintenance.

Activin-mediated alterations of the fibroblast transcriptome and matrisome control the biomechanical properties of skin wounds

Matrix deposition is essential for wound repair, but when excessive, leads to hypertrophic scars and fibrosis. The factors that control matrix deposition in skin wounds have only partially been identified and the consequences of matrix alterations for the mechanical properties of wounds are largely unknown. Here, we report how a single diffusible factor, activin A, affects the healing process across scales. Bioinformatics analysis of wound fibroblast transcriptome data combined with biochemical and histopathological analyses of wounds and functional in vitro studies identify that activin promotes pro-fibrotic gene expression signatures and processes, including glycoprotein and proteoglycan biosynthesis, collagen deposition, and altered collagen cross-linking. As a consequence, activin strongly reduces the wound and scar deformability, as identified by a non-invasive in vivo method for biomechanical analysis. These results provide mechanistic insight into the roles of activin in wound repair and fibrosis and identify the functional consequences of alterations in the wound matrisome at the biomechanical level.

The relationship between histopathology, gene expression, and biochemical and mechanical properties of wounds is largely unknown. Here, the authors show that activin A alters wound healing at multiple levels by promoting pro-fibrotic gene expression and matrix deposition, thereby affecting biomechanical properties of skin wounds.

SciBet as a portable and fast single cell type identifier

Fast, robust and technology-independent computational methods are needed for supervised cell type annotation of single-cell RNA sequencing data. We present SciBet, a supervised cell type identifier that accurately predicts cell identity for newly sequenced cells with order-of-magnitude speed advantage. We enable web client deployment of SciBet for rapid local computation without uploading local data to the server. Facing the exponential growth in the size of single cell RNA datasets, this user-friendly and cross-platform tool can be widely useful for single cell type identification.

Quantitative In-Depth Analysis of the Mouse Mast Cell Transcriptome Reveals Organ-Specific Mast Cell Heterogeneity

Mast cells (MCs) are primarily resident hematopoietic tissue cells that are localized at external and internal surfaces of the body where they act in the first line of defense. MCs are found in all studied vertebrates and have also been identified in tunicates, an early chordate. To obtain a detailed insight into the biology of MCs, here we analyzed the transcriptome of MCs from different mouse organs by RNA-seq and PCR-based transcriptomics. We show that MCs at different tissue locations differ substantially in their levels of transcripts coding for the most abundant MC granule proteins, even within the connective tissue type, or mucosal MC niches. We also demonstrate that transcript levels for the major granule proteins, including the various MC-restricted proteases and the heparin core protein, can be several orders of magnitude higher than those coding for various surface receptors and enzymes involved in protease activation, as well as enzymes involved in the synthesis of heparin, histamine, leukotrienes, and prostaglandins. Interestingly, our analyses revealed an almost complete absence in MCs of transcripts coding for cytokines at baseline conditions, indicating that cytokines are primarily produced by activated MCs. Bone marrow-derived MCs (BMMCs) are often used as equivalents of tissue MCs. Here, we show that these cells differ substantially from tissue MCs with regard to their transcriptome. Notably, they showed a transcriptome indicative of relatively immature cells, both with respect to the expression of granule proteases and of various enzymes involved in the processing/synthesis of granule compounds, indicating that care should be taken when extrapolating findings from BMMCs to the in vivo function of tissue-resident MCs. Furthermore, the latter finding indicates that the development of fully mature tissue-resident MCs requires a cytokine milieu beyond what is needed for in vitro differentiation of BMMCs. Altogether, this study provides a comprehensive quantitative view of the transcriptome profile of MCs resident at different tissue locations that builds nicely on previous studies of both the mouse and human transcriptome, and form a solid base for future evolutionary studies of the role of MCs in vertebrate immunity.

Combined single-cell and spatial transcriptomics reveal the molecular, cellular and spatial bone marrow niche organization

The bone marrow constitutes the primary site for life-long blood production and skeletal regeneration. However, its cellular and spatial organization remains controversial. Here, we combine single-cell and spatially resolved transcriptomics to systematically map the molecular, cellular and spatial composition of distinct bone marrow niches. This allowed us to transcriptionally profile all major bone-marrow-resident cell types, determine their localization and clarify sources of pro-haematopoietic factors. Our data demonstrate that Cxcl12-abundant-reticular (CAR) cell subsets (Adipo-CAR and Osteo-CAR) differentially localize to sinusoidal and arteriolar surfaces, act locally as 'professional cytokine-secreting cells' and thereby establish peri-vascular micro-niches. Importantly, the three-dimensional bone-marrow organization can be accurately inferred from single-cell transcriptome data using the RNA-Magnet algorithm described here. Together, our study reveals the cellular and spatial organization of bone marrow niches and offers a systematic approach to dissect the complex organization of whole organs.

Inflammasome-Independent Role for NLRP3 in Controlling Innate Antihelminth Immunity and Tissue Repair in the Lung

Alternatively activated macrophages are essential effector cells during type 2 immunity and tissue repair following helminth infections. We previously showed that Ym1, an alternative activation marker, can drive innate IL-1R-dependent neutrophil recruitment during infection with the lung-migrating nematode, Nippostrongylus brasiliensis, suggesting a potential role for the inflammasome in the IL-1-mediated innate response to infection. Although inflammasome proteins such as NLRP3 have important proinflammatory functions in macrophages, their role during type 2 responses and repair are less defined. We therefore infected Nlrp3 -/- mice with N. brasiliensis Unexpectedly, compared with wild-type (WT) mice, infected Nlrp3 -/- mice had increased neutrophilia and eosinophilia, correlating with enhanced worm killing but at the expense of increased tissue damage and delayed lung repair. Transcriptional profiling showed that infected Nlrp3 -/- mice exhibited elevated type 2 gene expression compared with WT mice. Notably, inflammasome activation was not evident early postinfection with N. brasiliensis, and in contrast to Nlrp3 -/- mice, antihelminth responses were unaffected in caspase-1/11-deficient or WT mice treated with the NLRP3-specific inhibitor MCC950. Together these data suggest that NLRP3 has a role in constraining lung neutrophilia, helminth killing, and type 2 immune responses in an inflammasome-independent manner.

Persistent Toll-like receptor 7 stimulation induces behavioral and molecular innate immune tolerance

Toll-like receptors 7 and 8 (TLR7 and TLR8) are endosomal pattern recognition receptors that detect a variety of single-stranded RNA species. While TLR7/8 agonists have robust therapeutic potential, clinical utility of these agents is limited by sickness responses associated with treatment induction. To understand the kinetics and mechanism of these responses, we characterized the acute and chronic effects of TLR7 stimulation. Single-cell RNA-sequencing studies, RNAscope, and radiolabeled in situ hybridization demonstrate that central nervous system gene expression of TLR7 is exclusive to microglia. In vitro studies demonstrate that microglia are highly sensitive to TLR7 stimulation, and respond in a dose-dependent manner to the imidazoquinoline R848. In vivo, both intraperitoneal (IP) and intracerebroventricular (ICV) R848 induce acute sickness responses including hypophagia, weight loss, and decreased voluntary locomotor activity, associated with increased CNS pro-inflammatory gene expression and changes to glial morphology. However, chronic daily IP R848 resulted in rapid tachyphylaxis of behavioral and molecular manifestations of illness. In microglial in vitro assays, pro-inflammatory transcriptional responses rapidly diminished in the context of repeated R848. In addition to TLR7 desensitization, we found that microglia become partially refractory to lipopolysaccharide (LPS) following R848 pretreatment, associated with induction of negative regulators A20 and Irak3. Similarly, mice pre-treated with R848 demonstrate reduced sickness responses, hypothalamic inflammation, and hepatic inflammation in response to LPS. These data combined demonstrate that TLR7 stimulation induces acute behavioral and molecular evidence of sickness responses. Following prolonged dosing, R848 induces a refractory state to both TLR7 and TLR4 activation, consistent with induced immune tolerance.