Integrative Computational Modeling of the Lymph Node Stromal Cell Landscape

Adaptive immune responses develop in secondary lymphoid organs such as lymph nodes (LNs) in a well-coordinated series of interactions between migrating immune cells and resident stromal cells. Although many processes that occur in LNs are well understood from an immunological point of view, our understanding of the fundamental organization and mechanisms that drive these processes is still incomplete. The aim of systems biology approaches is to unravel the complexity of biological systems and describe emergent properties that arise from interactions between individual constituents of the system. The immune system is greater than the sum of its parts, as is the case with any sufficiently complex system. Here, we review recent work and developments of computational LN models with focus on the structure and organization of the stromal cells. We explore various mathematical studies of intranodal T cell motility and migration, their interactions with the LN-resident stromal cells, and computational models of functional chemokine gradient fields and lymph flow dynamics. Lastly, we discuss briefly the importance of hybrid and multi-scale modeling approaches in immunology and the technical challenges involved

B-cell Zone Reticular Cell Microenvironments Shape CXCL13 Gradient Formation

Through the formation of concentration gradients, morphogens drive graded responses to extracellular signals, thereby fine-tuning cell behaviors in complex tissues. Here we show that the chemokine CXCL13 forms both soluble and immobilized gradients. Specifically, CXCL13+ follicular reticular cells form a small-world network of guidance structures, with computer simulations and optimization analysis predicting that immobilized gradients created by this network promote B-cell trafficking. Consistent with this prediction, imaging analysis show that CXCL13 binds to extracellular matrix components in situ, constraining its diffusion. CXCL13 solubilization requires the protease cathepsin B that cleaves CXCL13 into a stable product. Mice lacking cathepsin B display aberrant follicular architecture, a phenotype associated with effective B cell homing to but not within lymph nodes. Our data thus suggest that reticular cells of the B cell zone generate microenvironments that shape both immobilized and soluble CXCL13 gradient

Computational Approach to 3D Modeling of the Lymph Node Geometry

In this study we present a computational approach to the generation of the major geometric structures of an idealized murine lymph node (LN). In this generation, we consider the major compartments such as the subcapsular sinus, B cell follicles, trabecular and medullar sinuses, blood vessels and the T cell zone with a primary focus on the fibroblastic reticular cell (FRC) network. Confocal microscopy data of LN macroscopic structures and structural properties of the FRC network have been generated and utilized in the present model. The methodology sets a library of modules that can be used to assemble a solid geometric LN model and subsequently generate an adaptive mesh model capable of implementing transport phenomena. Overall, based on the use of high-resolution confocal microscopy and morphological analysis of cell 3D reconstructions, we have developed a computational model of the LN geometry, suitable for further investigation in studies of fluid transport and cell migration in this immunologically essential organ

Origin and differentiation trajectories of fibroblastic reticular cells in the splenic white pulp

The splenic white pulp is underpinned by poorly characterized stromal cells that demarcate distinct immune cell microenvironments. Here we establish fibroblastic reticular cell (FRC)-specific fate-mapping in mice to define their embryonic origin and differentiation trajectories. Our data show that all reticular cell subsets descend from multipotent progenitors emerging at embryonic day 19.5 from periarterial progenitors. Commitment of FRC progenitors is concluded during the first week of postnatal life through occupation of niches along developing central arterioles. Single cell transcriptomic analysis facilitated deconvolution of FRC differentiation trajectories and indicated that perivascular reticular cells function both as adult lymphoid organizer cells and mural cell progenitors. The lymphotoxin-β receptor-independent sustenance of postnatal progenitor stemness unveils that systemic immune surveillance in the splenic white pulp is governed through subset specification of reticular cells from a multipotent periarterial progenitor cell. In sum, the finding that discrete signaling events in perivascular niches determine the differentiation trajectories of reticular cell networks explains the development of distinct microenvironmental niches in secondary and tertiary lymphoid tissues that are crucial for the induction and regulation of innate and adaptive immune processes

Critical Issues in Modelling Lymph Node Physiology

In this study, we discuss critical issues in modelling the structure and function of lymph nodes (LNs), with emphasis on how LN physiology is related to its multi-scale structural organization. In addition to macroscopic domains such as B-cell follicles and the T cell zone, there are vascular networks which play a key role in the delivery of information to the inner parts of the LN, i.e., the conduit and blood microvascular networks. We propose object-oriented computational algorithms to model the 3D geometry of the fibroblastic reticular cell (FRC) network and the microvasculature. Assuming that a conduit cylinder is densely packed with collagen fibers, the computational flow study predicted that the diffusion should be a dominating process in mass transport than convective flow. The geometry models are used to analyze the lymph flow properties through the conduit network in unperturbed- and damaged states of the LN. The analysis predicts that elimination of up to 60%–90% of edges is required to stop the lymph flux. This result suggests a high degree of functional robustness of the network

Intestinal fibroblastic reticular cell niches control innate lymphoid cell homeostasis and function

Innate lymphoid cells (ILCs) govern immune cell homeostasis in the intestine and protect the host against microbial pathogens. Various cell-intrinsic pathways have been identified that determine ILC development and differentiation. However, the cellular components that regulate ILC sustenance and function in the intestinal lamina propria are less known. Using single-cell transcriptomic analysis of lamina propria fibroblasts, we identify fibroblastic reticular cells (FRCs) that underpin cryptopatches (CPs) and isolated lymphoid follicles (ILFs). Genetic ablation of lymphotoxin-beta receptor expression in Ccl19-expressing FRCs blocks the maturation of CPs into mature ILFs. Interactome analysis shows the major niche factors and processes underlying FRC-ILC crosstalk. In vivo validation confirms that a sustained lymphotoxin-driven feedforward loop of FRC activation including IL-7 generation is critical for the maintenance of functional ILC populations. In sum, our study indicates critical fibroblastic niches within the intestinal lamina propria that control ILC homeostasis and functionality and thereby secure protective gut immunity. Fibroblastic reticular cells (FRCs) support localisation of immune cells in secondary lymphoid tissues but less is known about the lamina propria. Here the authors use scRNA-seq and intestinal infection to characterise FRCs in the intestinal lamina propria and show specialised niches that foster innate lymphoid cells during homeostasis and infection

Topological Small-World Organization of the Fibroblastic Reticular Cell Network Determines Lymph Node Functionality.

Fibroblastic reticular cells (FRCs) form the cellular scaffold of lymph nodes (LNs) and establish distinct microenvironmental niches to provide key molecules that drive innate and adaptive immune responses and control immune regulatory processes. Here, we have used a graph theory-based systems biology approach to determine topological properties and robustness of the LN FRC network in mice. We found that the FRC network exhibits an imprinted small-world topology that is fully regenerated within 4 wk after complete FRC ablation. Moreover, in silico perturbation analysis and in vivo validation revealed that LNs can tolerate a loss of approximately 50% of their FRCs without substantial impairment of immune cell recruitment, intranodal T cell migration, and dendritic cell-mediated activation of antiviral CD8+ T cells. Overall, our study reveals the high topological robustness of the FRC network and the critical role of the network integrity for the activation of adaptive immune responses

Alternative NF-κB signaling regulates mTEC differentiation from podoplanin-expressing presursors in the cortico-medullary junction

The thymic epithelium forms specialized niches to enable thymocyte differentiation. While the common epithelial progenitor of medullary and cortical thymic epithelial cells (mTECs and cTECs) is well defined, early stages of mTEC lineage specification have remained elusive. Here, we utilized in vivo targeting of mTECs to resolve their differentiation pathways and to determine whether mTEC progenitors participate in thymocyte education. We found that mTECs descend from a lineage committed, podoplanin (PDPN)-expressing progenitor located at the cortico-medullary junction. PDPN(+) junctional TECs (jTECs) represent a distinct TEC population that builds the thymic medulla, but only partially supports negative selection and thymocyte differentiation. Moreover, conditional gene targeting revealed that abrogation of alternative NF-κB pathway signaling in the jTEC stage completely blocked mTEC development. Taken together, this study identifies jTECs as lineage-committed mTEC progenitors and shows that NF-κB-dependent progression of jTECs to mTECs is critical to secure central tolerance

Remodeling of light and dark zone follicular dendritic cells governs germinal center responses

Efficient generation of germinal center (GC) responses requires directed movement of B cells between distinct microenvironments underpinned by specialized B cell-interacting reticular cells (BRCs). How BRCs are reprogrammed to cater to the developing GC remains unclear, and studying this process is largely hindered by incomplete resolution of the cellular composition of the B cell follicle. Here we used genetic targeting of Cxcl13-expressing cells to define the molecular identity of the BRC landscape. Single-cell transcriptomic analysis revealed that BRC subset specification was predetermined in the primary B cell follicle. Further topological remodeling of light and dark zone follicular dendritic cells required CXCL12-dependent crosstalk with B cells and dictated GC output by retaining B cells in the follicle and steering their interaction with follicular helper T cells. Together, our results reveal that poised BRC-defined microenvironments establish a feed-forward system that determines the efficacy of the GC reaction