Author Correction: Commensal-driven immune zonation of the liver promotes host defence.

A Correction to this paper has been published: https://doi.org/10.1038/s41586-021-03346-0.

Commensal-driven immune zonation of the liver promotes host defence.

The liver connects the intestinal portal vasculature with the general circulation, using a diverse array of immune cells to protect from pathogens that translocate from the gut1. In liver lobules, blood flows from portal triads that are situated in periportal lobular regions to the central vein via a polarized sinusoidal network. Despite this asymmetry, resident immune cells in the liver are considered to be broadly dispersed across the lobule. This differs from lymphoid organs, in which immune cells adopt spatially biased positions to promote effective host defence2,3. Here we used quantitative multiplex imaging, genetic perturbations, transcriptomics, infection-based assays and mathematical modelling to reassess the relationship between the localization of immune cells in the liver and host protection. We found that myeloid and lymphoid resident immune cells concentrate around periportal regions. This asymmetric localization was not developmentally controlled, but resulted from sustained MYD88-dependent signalling induced by commensal bacteria in liver sinusoidal endothelial cells, which in turn regulated the composition of the pericellular matrix involved in the formation of chemokine gradients. In vivo experiments and modelling showed that this immune spatial polarization was more efficient than a uniform distribution in protecting against systemic bacterial dissemination. Together, these data reveal that liver sinusoidal endothelial cells sense the microbiome, actively orchestrating the localization of immune cells, to optimize host defence.

IBEX: A versatile multiplex optical imaging approach for deep phenotyping and spatial analysis of cells in complex tissues.

The diverse composition of mammalian tissues poses challenges for understanding the cell-cell interactions required for organ homeostasis and how spatial relationships are perturbed during disease. Existing methods such as single-cell genomics, lacking a spatial context, and traditional immunofluorescence, capturing only two to six molecular features, cannot resolve these issues. Imaging technologies have been developed to address these problems, but each possesses limitations that constrain widespread use. Here we report a method that overcomes major impediments to highly multiplex tissue imaging. "Iterative bleaching extends multiplexity" (IBEX) uses an iterative staining and chemical bleaching method to enable high-resolution imaging of >65 parameters in the same tissue section without physical degradation. IBEX can be employed with various types of conventional microscopes and permits use of both commercially available and user-generated antibodies in an "open" system to allow easy adjustment of staining panels based on ongoing marker discovery efforts. We show how IBEX can also be used with amplified staining methods for imaging strongly fixed tissues with limited epitope retention and with oligonucleotide-based staining, allowing potential cross-referencing between flow cytometry, cellular indexing of transcriptomes and epitopes by sequencing, and IBEX analysis of the same tissue. To facilitate data processing, we provide an open-source platform for automated registration of iterative images. IBEX thus represents a technology that can be rapidly integrated into most current laboratory workflows to achieve high-content imaging to reveal the complex cellular landscape of diverse organs and tissues.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes.

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

RNA-Seq identifies SPGs as a ventral skeletal patterning cue in sea urchins.

The sea urchin larval skeleton offers a simple model for formation of developmental patterns. The calcium carbonate skeleton is secreted by primary mesenchyme cells (PMCs) in response to largely unknown patterning cues expressed by the ectoderm. To discover novel ectodermal cues, we performed an unbiased RNA-Seq-based screen and functionally tested candidates; we thereby identified several novel skeletal patterning cues. Among these, we show that SLC26a2/7 is a ventrally expressed sulfate transporter that promotes a ventral accumulation of sulfated proteoglycans, which is required for ventral PMC positioning and skeletal patterning. We show that the effects of SLC perturbation are mimicked by manipulation of either external sulfate levels or proteoglycan sulfation. These results identify novel skeletal patterning genes and demonstrate that ventral proteoglycan sulfation serves as a positional cue for sea urchin skeletal patterning.

Neuronal maturation reduces the type I IFN response to orthobunyavirus infection and leads to increased apoptosis of human neurons.

BACKGROUND: La Crosse virus (LACV) is the leading cause of pediatric arboviral encephalitis in the USA. LACV encephalitis can result in learning and memory deficits, which may be due to infection and apoptosis of neurons in the brain. Despite neurons being the primary cell infected in the brain by LACV, little is known about neuronal responses to infection. METHODS: Human cerebral organoids (COs), which contain a spectrum of developing neurons, were used to examine neuronal responses to LACV. Plaque assay and quantitative reverse transcription (qRT) PCR were used to determine the susceptibility of COs to LACV infection. Immunohistochemistry, flow cytometry, and single-cell transcriptomics were used to determine specific neuronal subpopulation responses to the virus. RESULTS: Overall, LACV readily infected COs causing reduced cell viability and increased apoptosis. However, it was determined that neurons at different stages of development had distinct responses to LACV. Both neural progenitors and committed neurons were infected with LACV, however, committed neurons underwent apoptosis at a higher rate. Transcriptomic analysis showed that committed neurons expressed fewer interferon (IFN)-stimulated genes (ISGs) and genes involved IFN signaling in response to infection compared to neural progenitors. Furthermore, induction of interferon signaling in LACV-infected COs by application of recombinant IFN enhanced cell viability. CONCLUSIONS: These findings indicate that neuronal maturation increases the susceptibility of neurons to LACV-induced apoptosis. This susceptibility is likely due, at least in part, to mature neurons being less responsive to virus-induced IFN as evidenced by their poor ISG response to LACV. Furthermore, exogenous administration of recombinant IFN to LACV COs rescued cellular viability suggesting that increased IFN signaling is overall protective in this complex neural tissue. Together these findings indicate that induction of IFN signaling in developing neurons is an important deciding factor in virus-induced cell death.

HoTResDB: host transcriptional response database for viral hemorrhagic fevers.

SUMMARY: High-throughput screening of the host transcriptional response to various viral infections provides a wealth of data, but utilization of microarray and next generation sequencing (NGS) data for analysis can be difficult. The Host Transcriptional Response DataBase (HoTResDB), allows visitors to access already processed microarray and NGS data from non-human primate models of viral hemorrhagic fever to better understand the host transcriptional response. AVAILABILITY: HoTResDB is freely available at http://hotresdb.bu.edu.

T-Cell Receptor Diversity and the Control of T-Cell Homeostasis Mark Ebola Virus Disease Survival in Humans.

Differences in T-cell phenotype, particularly the expression of markers of T-cell homeostasis, have been observed in fatal and nonfatal Ebola virus disease (EVD). However, the relationship between these markers with T-cell function and virus clearance during EVD is poorly understood. To gain biological insight into the role of T cells during EVD, combined transcriptomics and T-cell receptor sequencing was used to profile blood samples from fatal and nonfatal EVD patients from the recent West African EVD epidemic. Fatal EVD was characterized by strong T-cell activation and increased abundance of T-cell inhibitory molecules. However, the early T-cell response was oligoclonal and did not result in viral clearance. In contrast, survivors mounted highly diverse T-cell responses, maintained low levels of T-cell inhibitors, and cleared Ebola virus. Our findings highlight the importance of T-cell immunity in surviving EVD and strengthen the foundation for further research on targeting of the dendritic cell-T cell interface for postexposure immunotherapy.

Comparative Transcriptomics in Ebola Makona-Infected Ferrets, Nonhuman Primates, and Humans.

The domestic ferret is a uniformly lethal model of infection for 3 species of Ebolavirus known to be pathogenic in humans. Reagents to systematically analyze the ferret host response to infection are lacking; however, the recent publication of a draft ferret genome has opened the potential for transcriptional analysis of ferret models of disease. In this work, we present comparative analysis of longitudinally sampled blood taken from ferrets and nonhuman primates infected with lethal doses of the Makona variant of Zaire ebolavirus. Strong induction of proinflammatory and prothrombotic signaling programs were present in both ferrets and nonhuman primates, and both transcriptomes were similar to previously published datasets of fatal cases of human Ebola virus infection.

Preservation of scRNA-Seq Libraries Using Existing Inactivation Protocols.

Single-cell RNA sequencing has soared in popularity in recent years. The ability to deeply profile the states of individual cells during the course of disease or infection has helped to expand our knowledge of coordinated responses. However, significant challenges arise when performing this analysis in high containment settings such as biosafety level 3 (BSL-3), BSL-3+ and BSL-4. Working in containment is necessary for many important pathogens, such as Ebola virus, Marburg virus, Lassa virus, Nipah and Hendra viruses. Since standard operating procedures (SOPs) for inactivation are extensive and may compromise sample integrity, we tested whether the removal of single-cell sequencing libraries from containment laboratories using existing inactivation protocols for nucleic acid extraction (Trizol, RLT buffer, or AVL buffer) was feasible. We have demonstrated that the inactivation does not affect sample quality and can work with existing methods for inactivation.

Rebalancing Viral and Immune Damage versus Tissue Repair Prevents Death from Lethal Influenza Infection.

Maintaining tissue function while eliminating infected cells is fundamental to host defense. Innate inflammatory damage contributes to lethal influenza and COVID-19, yet other than steroids, immunomodulatory drugs have modest effects. Among more than 50 immunomodulatory regimes tested in mouse lethal influenza infection, only the previously reported early depletion of neutrophils showed efficacy, suggesting that the infected host passes an early tipping point in which limiting innate immune damage alone cannot rescue physiological function. To re-balance the system late in infection, we investigated whether partial limitation of viral spread using oseltamivir (Tamiflu) together with enhancement of epithelial repair by blockade of interferon signaling or the limitation of further epithelial cell loss mediated by cytotoxic CD8 + T cells would prevent death. These treatments salvaged a large fraction of infected animals, providing new insight into the importance of repair processes and the timing of adaptive immune responses in survival of pulmonary infections.

Spatial Patterning Analysis of Cellular Ensembles (SPACE) discovers complex spatial organization at the cell and tissue levels.

Spatial patterns of cells and other biological elements drive both physiologic and pathologic processes within tissues. While many imaging and transcriptomic methods document tissue organization, discerning these patterns is challenging, especially when they involve multiple elements in complex arrangements. To address this challenge, we present Spatial Patterning Analysis of Cellular Ensembles (SPACE), an R package for analysis of high-plex spatial data. SPACE is compatible with any data collection modality that records values (i.e., categorical cell/structure types or quantitative expression levels) at fixed spatial coordinates (i.e., 2d pixels or 3d voxels). SPACE detects not only broad patterns of co-occurrence but also context-dependent associations, quantitative gradients and orientations, and other organizational complexities. Via a robust information theoretic framework, SPACE explores all possible ensembles of tissue elements - single elements, pairs, triplets, and so on - and ranks the most strongly patterned ensembles. For single images, rankings reflect patterns that differ from random assortment. For sets of images, rankings reflect patterns that differ across sample groups (e.g., genotypes, treatments, timepoints, etc.). Further tools then thoroughly characterize the nature of each pattern for intuitive interpretation. We validate SPACE and demonstrate its advantages using murine lymph node images for which ground truth has been defined. We then use SPACE to detect new patterns across varied datasets, including tumors and tuberculosis granulomas.

Multi-omic profiling of follicular lymphoma reveals changes in tissue architecture and enhanced stromal remodeling in high-risk patients.

Follicular lymphoma (FL) is a generally incurable malignancy that evolves from developmentally blocked germinal center (GC) B cells. To promote survival and immune escape, tumor B cells undergo significant genetic changes and extensively remodel the lymphoid microenvironment. Dynamic interactions between tumor B cells and the tumor microenvironment (TME) are hypothesized to contribute to the broad spectrum of clinical behaviors observed among FL patients. Despite the urgent need, existing clinical tools do not reliably predict disease behavior. Using a multi-modal strategy, we examined cell-intrinsic and -extrinsic factors governing progression and therapeutic outcomes in FL patients enrolled onto a prospective clinical trial. By leveraging the strengths of each platform, we identify several tumor-specific features and microenvironmental patterns enriched in individuals who experience early relapse, the most high-risk FL patients. These features include stromal desmoplasia and changes to the follicular growth pattern present 20 months before first progression and first relapse.

The IBEX Knowledge-Base: Achieving more together with open science.

Iterative Bleaching Extends multipleXity (IBEX) is a versatile method for highly multiplexed imaging of diverse tissues. Based on open science principles, we created the IBEX Knowledge-Base, a resource for reagents, protocols and more, to empower innovation.

Understanding virus-host interactions in tissues.

Although virus-host interactions are usually studied in a single cell type using in vitro assays in immortalized cell lines or isolated cell populations, it is important to remember that what is happening inside one infected cell does not translate to understanding how an infected cell behaves in a tissue, organ or whole organism. Infections occur in complex tissue environments, which contain a host of factors that can alter the course of the infection, including immune cells, non-immune cells and extracellular-matrix components. These factors affect how the host responds to the virus and form the basis of the protective response. To understand virus infection, tools are needed that can profile the tissue environment. This Review highlights methods to study virus-host interactions in the infection microenvironment.

IBEX: an iterative immunolabeling and chemical bleaching method for high-content imaging of diverse tissues.

High-content imaging is needed to catalog the variety of cellular phenotypes and multicellular ecosystems present in metazoan tissues. We recently developed iterative bleaching extends multiplexity (IBEX), an iterative immunolabeling and chemical bleaching method that enables multiplexed imaging (>65 parameters) in diverse tissues, including human organs relevant for international consortia efforts. IBEX is compatible with >250 commercially available antibodies and 16 unique fluorophores, and can be easily adopted to different imaging platforms using slides and nonproprietary imaging chambers. The overall protocol consists of iterative cycles of antibody labeling, imaging and chemical bleaching that can be completed at relatively low cost in 2-5 d by biologists with basic laboratory skills. To support widespread adoption, we provide extensive details on tissue processing, curated lists of validated antibodies and tissue-specific panels for multiplex imaging. Furthermore, instructions are included on how to automate the method using competitively priced instruments and reagents. Finally, we present a software solution for image alignment that can be executed by individuals without programming experience using open-source software and freeware. In summary, IBEX is a noncommercial method that can be readily implemented by academic laboratories and scaled to achieve high-content mapping of diverse tissues in support of a Human Reference Atlas or other such applications.

Age-related differences in immune dynamics during SARS-CoV-2 infection in rhesus macaques.

Advanced age is a key predictor of severe COVID-19. To gain insight into this relationship, we used the rhesus macaque model of SARS-CoV-2 infection. Eight older and eight younger macaques were inoculated with SARS-CoV-2. Animals were evaluated using viral RNA quantification, clinical observations, thoracic radiographs, single-cell transcriptomics, multiparameter flow cytometry, multiplex immunohistochemistry, cytokine detection, and lipidomics analysis at predefined time points in various tissues. Differences in clinical signs, pulmonary infiltrates, and virus replication were limited. Transcriptional signatures of inflammation-associated genes in bronchoalveolar lavage fluid at 3 dpi revealed efficient mounting of innate immune defenses in both cohorts. However, age-specific divergence of immune responses emerged during the post-acute phase. Older animals exhibited sustained local inflammatory innate responses, whereas local effector T-cell responses were induced earlier in the younger animals. Circulating lipid mediator and cytokine levels highlighted increased repair-associated signals in the younger animals, and persistent pro-inflammatory responses in the older animals. In summary, despite similar disease outcomes, multi-omics profiling suggests that age may delay or impair antiviral cellular immune responses and delay efficient return to immune homeostasis.

Pulmonary infection induces persistent, pathogen-specific lipidomic changes influencing trained immunity.

Resolution of infection results in development of trained innate immunity which is typically beneficial for defense against unrelated secondary infection. Epigenetic changes including modification of histones via binding of various polar metabolites underlie the establishment of trained innate immunity. Therefore, host metabolism and this response are intimately linked. However, little is known regarding the influence of lipids on the development and function of trained immunity. Utilizing two models of pulmonary bacterial infection combined with multi-omic approaches, we identified persistent, pathogen-specific changes to the lung lipidome that correlated with differences in the trained immune response against a third unrelated pathogen. Further, we establish the specific cellular populations in the lung that contribute to this altered lipidome. Together these results expand our understanding of the pulmonary trained innate immune response and the contributions of host lipids in informing that response.

Single-cell RNA sequencing reveals SARS-CoV-2 infection dynamics in lungs of African green monkeys.

Detailed knowledge about the dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is important for uncovering the viral and host factors that contribute to coronavirus disease 2019 (COVID-19) pathogenesis. Old-World nonhuman primates recapitulate mild to moderate cases of COVID-19, thereby serving as important pathogenesis models. We compared African green monkeys inoculated with infectious SARS-CoV-2 or irradiated, inactivated virus to study the dynamics of virus replication throughout the respiratory tract. Genomic RNA from the animals inoculated with the irradiated virus was found to be highly stable, whereas subgenomic RNA, an indicator of viral replication, was found to degrade quickly. We combined this information with single-cell RNA sequencing of cells isolated from the lung and lung-draining mediastinal lymph nodes and developed new analysis methods for unbiased targeting of important cells in the host response to SARS-CoV-2 infection. Through detection of reads to the viral genome, we were able to determine that replication of the virus in the lungs appeared to occur mainly in pneumocytes, whereas macrophages drove the inflammatory response. Monocyte-derived macrophages recruited to the lungs, rather than tissue-resident alveolar macrophages, were most likely to be responsible for phagocytosis of infected cells and cellular debris early in infection, with their roles switching during clearance of infection. Together, our dataset provides a detailed view of the dynamics of virus replication and host responses over the course of mild COVID-19 and serves as a valuable resource to identify therapeutic targets.