A method for the inference of cytokine interaction networks.

Cell-cell communication is mediated by many soluble mediators, including over 40 cytokines. Cytokines, e.g. TNF, IL1ß, IL5, IL6, IL12 and IL23, represent important therapeutic targets in immune-mediated inflammatory diseases (IMIDs), such as inflammatory bowel disease (IBD), psoriasis, asthma, rheumatoid and juvenile arthritis. The identification of cytokines that are causative drivers of, and not just associated with, inflammation is fundamental for selecting therapeutic targets that should be studied in clinical trials. As in vitro models of cytokine interactions provide a simplified framework to study complex in vivo interactions, and can easily be perturbed experimentally, they are key for identifying such targets. We present a method to extract a minimal, weighted cytokine interaction network, given in vitro data on the effects of the blockage of single cytokine receptors on the secretion rate of other cytokines. Existing biological network inference methods typically consider the correlation structure of the underlying dataset, but this can make them poorly suited for highly connected, non-linear cytokine interaction data. Our method uses ordinary differential equation systems to represent cytokine interactions, and efficiently computes the configuration with the lowest Akaike information criterion value for all possible network configurations. It enables us to study indirect cytokine interactions and quantify inhibition effects. The extracted network can also be used to predict the combined effects of inhibiting various cytokines simultaneously. The model equations can easily be adjusted to incorporate more complicated dynamics and accommodate temporal data. We validate our method using synthetic datasets and apply our method to an experimental dataset on the regulation of IL23, a cytokine with therapeutic relevance in psoriasis and IBD. We validate several model predictions against experimental data that were not used for model fitting. In summary, we present a novel method specifically designed to efficiently infer cytokine interaction networks from cytokine perturbation data in the context of IMIDs.

Analysis of cellular kinetic models suggest that physiologically based model parameters may be inherently, practically unidentifiable.

Physiologically-based pharmacokinetic and cellular kinetic models are used extensively to predict concentration profiles of drugs or adoptively transferred cells in patients and laboratory animals. Models are fit to data by the numerical optimisation of appropriate parameter values. When quantities such as the area under the curve are all that is desired, only a close qualitative fit to data is required. When the biological interpretation of the model that produced the fit is important, an assessment of uncertainties is often also warranted. Often, a goal of fitting PBPK models to data is to estimate parameter values, which can then be used to assess characteristics of the fit system or applied to inform new modelling efforts and extrapolation, to inform a prediction under new conditions. However, the parameters that yield a particular model output may not necessarily be unique, in which case the parameters are said to be unidentifiable. We show that the parameters in three published physiologically-based pharmacokinetic models are practically (deterministically) unidentifiable and that it is challenging to assess the associated parameter uncertainty with simple curve fitting techniques. This result could affect many physiologically-based pharmacokinetic models, and we advocate more widespread use of thorough techniques and analyses to address these issues, such as established Markov Chain Monte Carlo and Bayesian methodologies. Greater handling and reporting of uncertainty and identifiability of fit parameters would directly and positively impact interpretation and translation for physiologically-based model applications, enhancing their capacity to inform new model development efforts and extrapolation in support of future clinical decision-making.

Immunofibroblasts are pivotal drivers of tertiary lymphoid structure formation and local pathology.

Resident fibroblasts at sites of infection, chronic inflammation, or cancer undergo phenotypic and functional changes to support leukocyte migration and, in some cases, aggregation into tertiary lymphoid structures (TLS). The molecular programming that shapes these changes and the functional requirements of this population in TLS development are unclear. Here, we demonstrate that external triggers at mucosal sites are able to induce the progressive differentiation of a population of podoplanin (pdpn)-positive stromal cells into a network of immunofibroblasts that are able to support the earliest phases of TLS establishment. This program of events, that precedes lymphocyte infiltration in the tissue, is mediated by paracrine and autocrine signals mainly regulated by IL13. This initial fibroblast network is expanded and stabilized, once lymphocytes are recruited, by the local production of the cytokines IL22 and lymphotoxin. Interfering with this regulated program of events or depleting the immunofibroblasts in vivo results in abrogation of local pathology, demonstrating the functional role of immunofibroblasts in supporting TLS maintenance in the tissue and suggesting novel therapeutic targets in TLS-associated diseases.

Inducible ablation of CD11c+ cells to determine their role in skin wound repair.

Whether resident and recruited myeloid cells may impair or aid healing of acute skin wounds remains a debated question. To begin to address this, we examined the importance of CD11c+ myeloid cells in the early activation of skin wound repair. We find that an absence of CD11c+ cells delays wound closure and epidermal proliferation, likely due to defects in the activation of the IL-23-IL-22 axis that is required for wound healing.

MicroRNA-155 induction via TNF-α and IFN-γ suppresses expression of programmed death ligand-1 (PD-L1) in human primary cells.

Programmed death ligand-1 (PD-L1) is a critical regulator of T cell function contributing to peripheral immune tolerance. Although it has been shown that posttranscriptional regulatory mechanisms control PD-L1 expression in cancer, it remains unknown whether such regulatory loops operate also in non-transformed cells. Here we studied PD-L1 expression in human dermal lymphatic endothelial cells (HDLECs), which play key roles in immunity and cancer. Treatment of HDLECs with the pro-inflammatory cytokines IFN-γ and TNF-α synergistically up-regulated PD-L1 expression. IFN-γ and TNF-α also affected expression of several microRNAs (miRNAs) that have the potential to suppress PD-L1 expression. The most highly up-regulated miRNA following IFN-γ and TNF-α treatment in HDLECs was miR-155, which has a central role in the immune system and cancer. Induction of miR-155 was driven by TNF-α, the effect of which was significantly enhanced by IFN-γ. The PD-L1 3'-UTR contains two functional miR-155-binding sites. Endogenous miR-155 controlled the kinetics and maximal levels of PD-L1 induction upon IFN-γ and TNF-α treatments. We obtained similar findings in dermal fibroblasts, demonstrating that the IFN-γ/TNF-α/miR-155/PD-L1 pathway is not restricted to HDLECs. These results reveal miR-155 as a critical component of an inflammation-induced regulatory loop controlling PD-L1 expression in primary cells.

ASPASIA: A toolkit for evaluating the effects of biological interventions on SBML model behaviour.

A calibrated computational model reflects behaviours that are expected or observed in a complex system, providing a baseline upon which sensitivity analysis techniques can be used to analyse pathways that may impact model responses. However, calibration of a model where a behaviour depends on an intervention introduced after a defined time point is difficult, as model responses may be dependent on the conditions at the time the intervention is applied. We present ASPASIA (Automated Simulation Parameter Alteration and SensItivity Analysis), a cross-platform, open-source Java toolkit that addresses a key deficiency in software tools for understanding the impact an intervention has on system behaviour for models specified in Systems Biology Markup Language (SBML). ASPASIA can generate and modify models using SBML solver output as an initial parameter set, allowing interventions to be applied once a steady state has been reached. Additionally, multiple SBML models can be generated where a subset of parameter values are perturbed using local and global sensitivity analysis techniques, revealing the model's sensitivity to the intervention. To illustrate the capabilities of ASPASIA, we demonstrate how this tool has generated novel hypotheses regarding the mechanisms by which Th17-cell plasticity may be controlled in vivo. By using ASPASIA in conjunction with an SBML model of Th17-cell polarisation, we predict that promotion of the Th1-associated transcription factor T-bet, rather than inhibition of the Th17-associated transcription factor RORγt, is sufficient to drive switching of Th17 cells towards an IFN-γ-producing phenotype. Our approach can be applied to all SBML-encoded models to predict the effect that intervention strategies have on system behaviour. ASPASIA, released under the Artistic License (2.0), can be downloaded from http://www.york.ac.uk/ycil/software.

A Stromal Cell Niche for Human and Mouse Type 3 Innate Lymphoid Cells.

Adaptive immunity critically depends on the functional compartmentalization of secondary lymphoid organs. Mesenchymal stromal cells create and maintain specialized niches that support survival, activation, and expansion of T and B cells, and integrated analysis of lymphocytes and their niche has been instrumental in understanding adaptive immunity. Lymphoid organs are also home to type 3 innate lymphoid cells (ILC3), innate effector cells essential for barrier immunity. However, a specialized stromal niche for ILC3 has not been identified. A novel lineage-tracing approach now identifies a subset of murine fetal lymphoid tissue organizer cells that gives rise exclusively to adult marginal reticular cells. Moreover, both cell types are conserved from mice to humans and colocalize with ILC3 in secondary lymphoid tissues throughout life. In sum, we provide evidence that fetal stromal organizers give rise to adult marginal reticular cells and form a dedicated stromal niche for innate ILC3 in adaptive lymphoid organs.

Collagen-Poly(N-isopropylacrylamide) Hydrogels with Tunable Properties.

There is a lack of hydrogel materials whose properties can be tuned at the point of use. Biological hydrogels, such as collagen, gelate at physiological temperatures; however, they are not always ideal as scaffolds because of their low mechanical strength. Their mechanics can be improved through cross-linking and chemical modification, but these methods still require further synthesis. We have demonstrated that by combining collagen with a thermoresponsive polymer, poly(N-isopropylacrylamide) (PNIPAM), the mechanical properties can be improved while maintaining cytocompatibility. Furthermore, different concentrations of this polymer led to a range of hydrogels with shear moduli ranging from 10(5) Pa down to less than 10(2) Pa, similar to the soft tissues in the body. In addition to variable mechanical properties, the hydrogel blends have a range of micron-scale structures and porosities, which caused adipose-derived stromal cells (ADSCs) to adopt different morphologies when encapsulated within and may therefore be able to direct cell fate.

Blood flukes exploit Peyer's Patch lymphoid tissue to facilitate transmission from the mammalian host.

Schistosomes are blood-dwelling parasitic helminths which produce eggs in order to facilitate transmission. Intestinal schistosomes lay eggs in the mesenteries, however, it is unclear how their eggs escape the vasculature to exit the host. Using a murine model of infection, we reveal that Schistosoma mansoni exploits Peyer's Patches (PP) gut lymphoid tissue as a preferential route of egress for their eggs. Egg deposition is favoured within PP as a result of their more abundant vasculature. Moreover, the presence of eggs causes significant vascular remodeling leading to an expanded venule network. Egg deposition results in a decrease in stromal integrity and lymphoid cellularity, including secretory IgA producing lymphocytes, and the focal recruitment of macrophages. In mice lacking PP, egg excretion is significantly impaired, leading to greater numbers of ova being entrapped in tissues and consequently, exacerbated morbidity. Thus, we demonstrate how schistosomes directly facilitate transmission from the host by targeting lymphoid tissue. For the host, PP-dependency of egg egress represents a trade-off, as limiting potentially life-threatening morbidity is balanced by loss of PP structure and perturbed PP IgA production.

Tyrosine kinase receptor RET is a key regulator of Peyer's patch organogenesis.

Normal organogenesis requires co-ordinate development and interaction of multiple cell types, and is seemingly governed by tissue specific factors. Lymphoid organogenesis during embryonic life is dependent on molecules the temporal expression of which is tightly regulated. During this process, haematopoietic 'inducer' cells interact with stromal 'organizer' cells, giving rise to the lymphoid organ primordia. Here we show that the haematopoietic cells in the gut exhibit a random pattern of motility before aggregation into the primordia of Peyer's patches, a major component of the gut-associated lymphoid tissue. We further show that a CD45+CD4-CD3-Il7Ralpha-c-Kit+CD11c+ haematopoietic population expressing lymphotoxin has an important role in the formation of Peyer's patches. A subset of these cells expresses the receptor tyrosine kinase RET, which is essential for mammalian enteric nervous system formation. We demonstrate that RET signalling is also crucial for Peyer's patch formation. Functional genetic analysis revealed that Gfra3-deficiency results in impairment of Peyer's patch development, suggesting that the signalling axis RET/GFRalpha3/ARTN is involved in this process. To support this hypothesis, we show that the RET ligand ARTN is a strong attractant of gut haematopoietic cells, inducing the formation of ectopic Peyer's patch-like structures. Our work strongly suggests that the RET signalling pathway, by regulating the development of both the nervous and lymphoid system in the gut, has a key role in the molecular mechanisms that orchestrate intestine organogenesis.

EphB-ephrin-B2 interactions are required for thymus migration during organogenesis.

Thymus organogenesis requires coordinated interactions of multiple cell types, including neural crest (NC) cells, to orchestrate the formation, separation, and subsequent migration of the developing thymus from the third pharyngeal pouch to the thoracic cavity. The molecular mechanisms driving these processes are unclear; however, NC-derived mesenchyme has been shown to play an important role. Here, we show that, in the absence of ephrin-B2 expression on thymic NC-derived mesenchyme, the thymus remains in the cervical area instead of migrating into the thoracic cavity. Analysis of individual NC-derived thymic mesenchymal cells shows that, in the absence of ephrin-B2, their motility is impaired as a result of defective EphB receptor signaling. This implies a NC-derived cell-specific role of EphB-ephrin-B2 interactions in the collective migration of the thymic rudiment during organogenesis.

Cellular characterisation of advanced osteoarthritis knee synovium.

OBJECTIVES: Osteoarthritis (OA) is increasingly recognised as a whole joint disease, with an important role for synovium. However, the repertoire of immune cells and fibroblasts that constitute OA synovium remains understudied. This study aims to characterise the cellular composition of advanced OA synovium and to explore potential correlations between different cell types and patient demographics or clinical scores. METHODS: Synovium, collected from 10 patients with advanced OA during total knee replacement surgery, was collagenase-digested, and cells were stained for flow cytometry analysis. Formalin-fixed paraffin-embedded synovium was sectioned, stained with immunofluorescence, and imaged using the multiplex Cell DIVE platform. Patient demographics and clinical scores were also collected. RESULTS: The proportion of immune cells in OA synovium varied between patients (8-38% of all cells). Macrophages and T cells were the dominant immune cell populations, together representing 76% of immune cells. Age positively correlated with the proportion of macrophages, and negatively correlated with T cells. CCR6+ T cells were found in 6/10 patients; these patients had a higher mean Kellgren-Lawrence grade across the three knee compartments. Immunofluorescence staining showed that macrophages were present in the lining as well as distributed throughout the sublining, while T and B cells were mainly localised near vessels in the sublining. Fibroblast subsets (CD45-PDPN+) based on the expression of CD34/CD90 or FAP/CD90 were identified in all patient samples, and some populations correlate with the percentage of immune cells or clinical scores. Immunofluorescence staining showed that FAP expression was particularly strong in the lining layer, but also present throughout the sublining layer. CD90 expression was exclusively found around vessels in the sublining, while CD34 was mostly found in the sublining but also occasionally in the lining layer. CONCLUSIONS: There are significant differences in the relative proportions and subsets of immune cells in OA synovium; exploratory correlative analyses suggest that these differences might be correlated with age, clinical scores, or fibroblast subsets. Additional studies are required to understand how different cell types affect OA pathobiology, and if the presence or proportion of cell subsets relates to disease phenotypes.

De-risking clinical trial failure through mechanistic simulation.

Drug development typically comprises a combination of pre-clinical experimentation, clinical trials, and statistical data-driven analyses. Therapeutic failure in late-stage clinical development costs the pharmaceutical industry billions of USD per year. Clinical trial simulation represents a key derisking strategy and combining them with mechanistic models allows one to test hypotheses for mechanisms of failure and to improve trial designs. This is illustrated with a T-cell activation model, used to simulate the clinical trials of IMA901, a short-peptide cancer vaccine. Simulation results were consistent with observed outcomes and predicted that responses are limited by peptide off-rates, peptide competition for dendritic cell (DC) binding, and DC migration times. These insights were used to hypothesise alternate trial designs predicted to improve efficacy outcomes. This framework illustrates how mechanistic models can complement clinical, experimental, and data-driven studies to understand, test, and improve trial designs, and how results may differ between humans and mice.

Alveolar macrophages transport pathogens to lung draining lymph nodes.

The first step in inducing pulmonary adaptive immunity to allergens and airborne pathogens is Ag acquisition and transport to the lung draining lymph nodes (dLN). Dendritic cells (DC) sample the airways, and active transfer of Ag to the lung dLN is considered an exclusive property of migratory DC. However, alveolar macrophages (AM) are the first phagocytes to contact inhaled particulate matter. Although having well-defined immunoregulatory capabilities, AM are generally considered as restricted to the alveoli. We show that murine AM constitutively migrate from lung to dLN and that following exposure to Streptococcus pneumoniae, AM rapidly transport bacteria to this site. Thus AM, and not DC, appear responsible for the earliest delivery of these bacteria to secondary lymphoid tissue. The identification of this novel transport pathway has important consequences for our understanding of lung immunity and suggests more widespread roles for macrophages in the transport of Ags to lymphoid organs than previously appreciated.

Simulating CXCR5 Dynamics in Complex Tissue Microenvironments.

To effectively navigate complex tissue microenvironments, immune cells sense molecular concentration gradients using G-protein coupled receptors. However, due to the complexity of receptor activity, and the multimodal nature of chemokine gradients in vivo, chemokine receptor activity in situ is poorly understood. To address this issue, we apply a modelling and simulation approach that permits analysis of the spatiotemporal dynamics of CXCR5 expression within an in silico B-follicle with single-cell resolution. Using this approach, we show that that in silico B-cell scanning is robust to changes in receptor numbers and changes in individual kinetic rates of receptor activity, but sensitive to global perturbations where multiple parameters are altered simultaneously. Through multi-objective optimization analysis we find that the rapid modulation of CXCR5 activity through receptor binding, desensitization and recycling is required for optimal antigen scanning rates. From these analyses we predict that chemokine receptor signaling dynamics regulate migration in complex tissue microenvironments to a greater extent than the total numbers of receptors on the cell surface.

Quantifying the limits of CAR T-cell delivery in mice and men.

CAR (Chimeric Antigen Receptor) T cells have demonstrated clinical success for the treatment of multiple lymphomas and leukaemias, but not for various solid tumours, despite promising data from murine models. Lower effective CAR T-cell delivery rates to human solid tumours compared to haematological malignancies in humans and solid tumours in mice might partially explain these divergent outcomes. We used anatomical and physiological data for human and rodent circulatory systems to calculate the typical perfusion of healthy and tumour tissues, and estimated the upper limits of immune cell delivery rates across different organs, tumour types and species. Estimated maximum delivery rates were up to 10 000-fold greater in mice than humans yet reported CAR T-cell doses are typically only 10-100-fold lower in mice, suggesting that the effective delivery rates of CAR T cells into tumours in clinical trials are far lower than in corresponding mouse models. Estimated delivery rates were found to be consistent with published positron emission tomography data. Results suggest that higher effective human doses may be needed to drive efficacy comparable to mouse solid tumour models, and that lower doses should be tested in mice. We posit that quantitation of species and organ-specific delivery and homing of engineered T cells will be key to unlocking their potential for solid tumours.

The Cellular Composition of the Uveal Immune Environment.

The uveal tract consists of the iris, the ciliary body and the choroid; these three distinct tissues form a continuous layer within the eye. Uveitis refers to inflammation of any region of the uveal tract. Despite being grouped together anatomically, the iris, ciliary body and choroid are distinct functionally, and inflammatory diseases may affect only one part and not the others. Cellular structure of tissues direct their function, and understanding the cellular basis of the immune environment of a tissue in health, the "steady state" on which the perturbations of disease are superimposed, is vital to understanding the pathogenesis of those diseases. A contemporary understanding of the immune system accepts that haematopoietic and yolk sac derived leukocytes, though vital, are not the only players of importance. An array of stromal cells, connective tissue cells such as fibroblasts and endothelial cells, may also have a role in the inflammatory reaction seen in several immune-mediated diseases. In this review we summarise what is known about the cellular composition of the uveal tract and the roles these disparate cell types have to play in immune homeostasis. We also discuss some unanswered questions surrounding the constituents of the resident leukocyte population of the different uveal tissues, and we look ahead to the new understanding that modern investigative techniques such as single cell transcriptomics, multi-omic data integration and highly-multiplexed imaging techniques may bring to the study of the uvea and uveitis, as they already have to other immune mediated inflammatory diseases.

Predicted limited redistribution of T cells to secondary lymphoid tissue correlates with increased risk of haematological malignancies in asplenic patients.

The spleen, a secondary lymphoid tissue (SLT), has an important role in generation of adaptive immune responses. Although splenectomy remains a common procedure, recent studies reported poor prognosis and increased risk of haematological malignancies in asplenic patients. The high baseline trafficking of T lymphocytes to splenic tissue suggests splenectomy may lead to loss of blood-borne malignant immunosurveillance that is not compensated for by the remaining SLT. To date, no quantitative analysis of the impact of splenectomy on the human T cell trafficking dynamics and tissue localisation has been reported. We developed a quantitative computational model that describes organ distribution and trafficking of human lymphocytes to explore the likely impact of splenectomy on immune cell distributions. In silico splenectomy resulted in an average reduction of T cell numbers in SLT by 35% (95%CI 0.12-0.97) and a comparatively lower, 9% (95%CI 0.17-1.43), mean decrease of T cell concentration in SLT. These results suggest that the surveillance capacity of the remaining SLT insufficiently compensates for the absence of the spleen. This may, in part, explain haematological malignancy risk in asplenic patients and raises the question of whether splenectomy has a clinically meaningful impact on patient responses to immunotherapy.