Applying Mathematical Models to Study the Role of the Immune System in Chronic Myelogenous Leukemia

Although tyrosine kinase inhibitors (TKIs) such as imatinib have transformed chronic myelogenous leukemia (CML) into a chronic condition, these therapies are not curative in the majority of cases. Most patients must continue TKI therapy indefinitely, a requirement that is both expensive and that compromises a patient's quality of life. While TKIs are known to reduce leukemic cells' proliferative capacity and to induce apoptosis, their effects on leukemic stem cells, the immune system, and the microenvironment are not fully understood. A more complete understanding of their global therapeutic effects would help us to identify any limitations of TKI monotherapy and to address these issues through novel combination therapies. Mathematical models are a complementary tool to experimental and clinical data that can provide valuable insights into the underlying mechanisms of TKI therapy. Previous modeling efforts have focused on CML patients who show biphasic and triphasic exponential declines in BCR-ABL ratio during therapy. However, our patient data indicates that many patients treated with TKIs show fluctuations in BCR-ABL ratio yet are able to achieve durable remissions. To investigate these fluctuations, we construct a mathematical model that integrates CML with a patient's autologous immune response to the disease. In our model, we define an immune window, which is an intermediate range of leukemic concentrations that lead to an effective immune response against CML. While small leukemic concentrations provide insufficient stimulus, large leukemic concentrations actively suppress a patient's immune system, thus limiting it's ability to respond. Our patient data and modeling results suggest that at diagnosis, a patient's high leukemic concentration is able to suppress their immune system. TKI therapy drives the leukemic population into the immune window, allowing the patient's immune cells to expand and eventually mount an efficient response against the residual CML. This response drives the leukemic population below the immune window, causing the immune population to contract and allowing the leukemia to partially recover. The leukemia eventually reenters the immune window, thus stimulating a sequence of weaker immune responses as the two populations approach equilibrium. We hypothesize that a patient's autologous immune response to CML may explain the fluctuations in BCR-ABL ratio that are regularly seen during TKI therapy. These fluctuations may serve as a signature of a patient's individual immune response to CML. By applying our modeling framework to patient data, we are able to construct an immune profile that can then be used to propose patient-specific combination therapies aimed at further reducing a patient's leukemic burden. Our characterization of a patient's anti-leukemia immune response may be especially valuable in the study of drug resistance, treatment cessation, and combination therapy.

Analysis and modelling of immune cell behaviour in lymph nodes based on multi-photon imaging

This thesis presents quantitative studies of T cell and dendritic cell (DC) behaviour in mouse lymph nodes (LNs) in the naive state and following immunisation. These processes are of importance and interest in basic immunology, and better understanding could improve both diagnostic capacity and therapeutic manipulations, potentially helping in producing more effective vaccines or developing treatments for autoimmune diseases. The problem is also interesting conceptually as it is relevant to other fields where 3D movement of objects is tracked with a discrete scanning interval. A general immunology introduction is presented in chapter 1. In chapter 2, I apply quantitative methods to multi-photon imaging data to measure how T cells and DCs are spatially arranged in LNs. This has been previously studied to describe differences between the naive and immunised state and as an indicator of the magnitude of the immune response in LNs, but previous analyses have been generally descriptive. The quantitative analysis shows that some of the previous conclusions may have been premature. In chapter 3, I use Bayesian state-space models to test some hypotheses about the mode of T cell search for DCs. A two-state mode of movement where T cells can be classified as either interacting to a DC or freely migrating is supported over a model where T cells would home in on DCs at distance through for example the action of chemokines. In chapter 4, I study whether T cell migration is linked to the geometric structure of the fibroblast reticular network (FRC). I find support for the hypothesis that the movement is constrained to the fibroblast reticular cell (FRC) network over an alternative 'random walk with persistence time' model where cells would move randomly, with a short-term persistence driven by a hypothetical T cell intrinsic 'clock'. I also present unexpected results on the FRC network geometry. Finally, a quantitative method is presented for addressing some measurement biases inherent to multi-photon imaging. In all three chapters, novel findings are made, and the methods developed have the potential for further use to address important problems in the field. In chapter 5, I present a summary and synthesis of results from chapters 3-4 and a more speculative discussion of these results and potential future directions.

Chemokine Receptor Expression on Normal Blood CD56(+) NK-Cells Elucidates Cell Partners That Comigrate during the Innate and Adaptive Immune Responses and Identifies a Transitional NK-Cell Population

Studies of chemokine receptors (CKR) in natural killer- (NK-) cells have already been published, but only a few gave detailed information on its differential expression on blood NK-cell subsets. We report on the expression of the inflammatory and homeostatic CKR on normal blood CD56(+low) CD16(+) and CD56(+high)  CD16(-/+low) NK-cells. Conventional CD56(+low) and CD56(+high) NK-cells present in the normal PB do express CKR for inflammatory cytokines, although with different patterns CD56(+low) NK-cells are mainly CXCR1/CXCR2(+) and CXCR3/CCR5(-/+), whereas mostly CD56(+high) NK-cells are CXCR1/CXCR2(-) and CXCR3/CCR5(+). Both NK-cell subsets have variable CXCR4 expression and are CCR4(-) and CCR6(-). The CKR repertoire of the CD56(+low) NK-cells approaches to that of neutrophils, whereas the CKR repertoire of the CD56(+high) NK-cells mimics that of Th1(+) T cells, suggesting that these cells are prepared to migrate into inflamed tissues at different phases of the immune response. In addition, we describe a subpopulation of NK-cells with intermediate levels of CD56 expression, which we named CD56(+int) NK-cells. These NK-cells are CXCR3/CCR5(+), they have intermediate levels of expression of CD16, CD62L, CD94, and CD122, and they are CD57(-) and CD158a(-). In view of their phenotypic features, we hypothesize that they correspond to a transitional stage, between the well-known CD56(+high) and CD56(+low) NK-cells populations.

Idiotypic Immune Networks in Mobile Robot Control

Jerne's idiotypic network theory postulates that the immune response involves inter-antibody stimulation and suppression as well as matching to antigens. The theory has proved the most popular Artificial Immune System (AIS) model for incorporation into behavior-based robotics but guidelines for implementing idiotypic selection are scarce. Furthermore, the direct effects of employing the technique have not been demonstrated in the form of a comparison with non-idiotypic systems. This paper aims to address these issues. A method for integrating an idiotypic AIS network with a Reinforcement Learning based control system (RL) is described and the mechanisms underlying antibody stimulation and suppression are explained in detail. Some hypotheses that account for the network advantage are put forward and tested using three systems with increasing idiotypic complexity. The basic RL, a simplified hybrid AIS-RL that implements idiotypic selection independently of derived concentration levels and a full hybrid AIS-RL scheme are examined. The test bed takes the form of a simulated Pioneer robot that is required to navigate through maze worlds detecting and tracking door markers.

Investigating the CD4+ T cell response to infection with Trypanosoma brucei

African trypanosomes are a protozoan parasite with the ability to cause disease states in both humans and animals; rendering them an important and relevant subject of study. These diseases have a vast socioeconomic impact upon the African continent and are perpetuated in part by the parasite's ability to evade the adaptive immune response. This thesis describes CD4+ T cell response to trypanosome infection, through the use of murine infection models and in vitro assays.

Idiotypic Immune Networks in Mobile Robot Control

Jerne's idiotypic network theory postulates that the immune response involves inter-antibody stimulation and suppression as well as matching to antigens. The theory has proved the most popular Artificial Immune System (AIS) model for incorporation into behavior-based robotics but guidelines for implementing idiotypic selection are scarce. Furthermore, the direct effects of employing the technique have not been demonstrated in the form of a comparison with non-idiotypic systems. This paper aims to address these issues. A method for integrating an idiotypic AIS network with a Reinforcement Learning based control system (RL) is described and the mechanisms underlying antibody stimulation and suppression are explained in detail. Some hypotheses that account for the network advantage are put forward and tested using three systems with increasing idiotypic complexity. The basic RL, a simplified hybrid AIS-RL that implements idiotypic selection independently of derived concentration levels and a full hybrid AIS-RL scheme are examined. The test bed takes the form of a simulated Pioneer robot that is required to navigate through maze worlds detecting and tracking door markers.

Novel strategies for targeting innate immune responses to influenza.

NlmCategory="UNASSIGNED">We previously reported that TLR4(-/-) mice are refractory to mouse-adapted A/PR/8/34 (PR8) influenza-induced lethality and that therapeutic administration of the TLR4 antagonist Eritoran blocked PR8-induced lethality and acute lung injury (ALI) when given starting 2 days post infection. Herein we extend these findings: anti-TLR4- or -TLR2-specific IgG therapy also conferred significant protection of wild-type (WT) mice from lethal PR8 infection. If treatment is initiated 3 h before PR8 infection and continued daily for 4 days, Eritoran failed to protect WT and TLR4(-/-) mice, implying that Eritoran must block a virus-induced, non-TLR4 signal that is required for protection. Mechanistically, we determined that (i) Eritoran blocks high-mobility group B1 (HMGB1)-mediated, TLR4-dependent signaling in vitro and circulating HMGB1 in vivo, and an HMGB1 inhibitor protects against PR8; (ii) Eritoran inhibits pulmonary lung edema associated with ALI; (iii) interleukin (IL)-1β contributes significantly to PR8-induced lethality, as evidenced by partial protection by IL-1 receptor antagonist (IL-1Ra) therapy. Synergistic protection against PR8-induced lethality was achieved when Eritoran and the antiviral drug oseltamivir were administered starting 4 days post infection. Eritoran treatment does not prevent development of an adaptive immune response to subsequent PR8 challenge. Overall, our data support the potential of a host-targeted therapeutic approach to influenza infection.Mucosal Immunology advance online publication 27 January 2016; doi:10.1038/mi.2015.141.

Distinct immune responses of juvenile and adult oysters (Crassostrea gigas) to viral and bacterial infections

Since 2008, massive mortality events of Pacific oysters (Crassostrea gigas) have been reported worldwide and these disease events are often associated with Ostreid herpesvirus type 1 (OsHV-1). Epidemiological field studies have also reported oyster age and other pathogens of the Vibrio genus are contributing factors to this syndrome. We undertook a controlled laboratory experiment to simultaneously investigate survival and immunological response of juvenile and adult C. gigas at different time-points post-infection with OsHV-1, Vibrio tasmaniensis LGP32 and V. aestuarianus. Our data corroborates epidemiological studies that juveniles are more susceptible to OsHV-1, whereas adults are more susceptible to Vibrio. We measured the expression of 102 immune-genes by high-throughput RT-qPCR, which revealed oysters have different transcriptional responses to OsHV-1 and Vibrio. The transcriptional response in the early stages of OsHV-1 infection involved genes related to apoptosis and the interferon-pathway. Transcriptional response to Vibrio infection involved antimicrobial peptides, heat shock proteins and galectins. Interestingly, oysters in the later stages of OsHV-1 infection had a transcriptional response that resembled an antibacterial response, which is suggestive of the oyster's microbiome causing secondary infections (dysbiosis-driven pathology). This study provides molecular evidence that oysters can mount distinct immune response to viral and bacterial pathogens and these responses differ depending on the age of the host.

Early growth response genes 2 and 3 regulate the expression of Bcl6 and differentiation of T follicular helper cells

T follicular helper (Tfh) cells support differentiation of B cells to plasma cells and high affinity antibody production in germinal centers (GC) and Tfh differentiation requires the function of B cell lymphoma 6 (Bcl6). We have now discovered that early growth response gene (Egr) 2 and 3 directly regulate the expression of Bcl6 in Tfh cells which is required for their function in regulation of GC formation. In the absence of Egr2 and 3, the expression of Bcl6 in Tfh cells is defective leading to impaired differentiation of Tfh cells resulting in a failure to form GCs following virus infection and defects in production of anti-viral antibodies. Enforced expression of Bcl6 in Egr2/3 deficient CD4 T cells partially restored Tfh differentiation and GC formation in response to virus infection. Our findings demonstrate a novel function of Egr2/3 which is important for Tfh cell development and Tfh cell mediated B cell immune responses.

Correlation between TH1 response standard cytokines as biomarkers in patients with the delta virus in the western Brazilian Amazon

Hepatitis D virus (HDV) is endemic in the Amazon Region and its pathophysiology is the most severe among viral hepatitis. Treatment is performed with pegylated interferon and the immune response appears to be important for infection control. HDV patients were studied: untreated and polymerase chain reaction (PCR) positive (n = 9), anti- HDV positive and PCR negative (n = 8), and responders to treatment (n = 12). The cytokines, interleukin (IL)-2 (p = 0.0008) and IL-12 (p = 0.02) were differentially expressed among the groups and were also correlated (p = 0.0143). Future studies will be conducted with patients at different stages of treatment, associating the viral load with serum cytokines produced, thereby attempting to establish a prognostic indicator of the infection.

Estudio de la capacidad inmunoprofiláctica de proteínas recombinantes de Leishmania infantum en el modelo murino

En esta Tesis Doctoral, cuatro proteínas que están sobre-expresadas en las fases infectivas del ciclo celular de L. infantum (STPKA, STPK, PP2B and PABP3) fueron purificadas para inmunizar ratones BALB/c, con el interés de estudiar su posible capacidad protectora contra la VL. Tres de ellas indujeron una fuerte respuesta humoral en los ratones, rLiSTPKA, rLiSTPK y rLiPABP3, a diferencia de la proteína rLiPP2B. Además, la rLiPABP3 fue la única que indujo una disminución significativa de la carga parasitara tanto en bazo como en hígado 2 meses después de la infección experimental con el parásito. Con el objetivo de determinar los efectos que la inmunización con rLiPABP3 tenía sobre el sistema inmune murino, se determinó la expresión génica de 106 genes relacionados con el sistema inmune en ratones control, inmunizados e infectados. Los niveles de expresión génica fueron comparados con los obtenidos para los grupos control. Los resultados mostraron que durante todo el experimento la proteína rLiPABP3 promueve la inhibición de la respuesta inmune inflamatoria en el bazo de los ratones infectados. Además, no se observa una respuesta adaptiva claramente polarizada hacia un tipo concreto. Un mes después de la infeccion, en los animales previamente inmunizados se observa la sobre-expresion de Il2rb lo que nos hace pensar que rLiPABP3 podría estimular la presencia de linfocitos T memoria. En fases más avanzadas de la infección, a pesar de que no observamos una clara diferenciación de una población concreta de linfocitos T efectores, sí observamos la infra-expresión del gen codificante del TNF-alfa, lo que unido a la sobre-expresión del gen Cxcr4 y la ausencia de cambios de la expresión del gen Ccr7 nos hace pensar que la inmunización no sólo mantiene la micro-arquitectura del bazo, sino que también promueve la correcta migración de las DCs desde la MZ hasta el PALS.

Transcriptional analysis of PRRSV-infected porcine dendritic cell response to Streptococcus suis infection reveals up-regulation of inflammatory-related genes expression

The porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most important swine pathogens and often serves as an entry door for other viral or bacterial pathogens, of which Streptococcus suis is one of the most common. Pre-infection with PRRSV leads to exacerbated disease caused by S. suis infection. Very few studies have assessed the immunological mechanisms underlying this higher susceptibility. Since antigen presenting cells play a major role in the initiation of the immune response, the in vitro transcriptional response of bone marrow-derived dendritic cells (BMDCs) and monocytes in the context of PRRSV and S. suis co-infection was investigated. BMDCs were found to be more permissive than monocytes to PRRSV infection; S. suis phagocytosis by PRRSV-infected BMDCs was found to be impaired, whereas no effect was found on bacterial intracellular survival. Transcription profile analysis, with a major focus on inflammatory genes, following S. suis infection, with and without pre-infection with PRRSV, was then performed. While PRRSV pre-infection had little effect on monocytes response to S. suis infection, a significant expression of several pro-inflammatory molecules was observed in BMDCs pre-infected with PRRSV after a subsequent infection with S. suis. While an additive effect could be observed for CCL4, CCL14, CCL20, and IL-15, a distinct synergistic up-regulatory effect was observed for IL-6, CCL5 and TNF-α after co-infection. This increased pro-inflammatory response by DCs could participate in the exacerbation of the disease observed during PRRSV and S. suis co-infection.

Transcriptional analysis of PRRSV-infected porcine dendritic cell response to Streptococcus suis infection reveals up-regulation of inflammatory-related genes expression

The porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most important swine pathogens and often serves as an entry door for other viral or bacterial pathogens, of which Streptococcus suis is one of the most common. Pre-infection with PRRSV leads to exacerbated disease caused by S. suis infection. Very few studies have assessed the immunological mechanisms underlying this higher susceptibility. Since antigen presenting cells play a major role in the initiation of the immune response, the in vitro transcriptional response of bone marrow-derived dendritic cells (BMDCs) and monocytes in the context of PRRSV and S. suis co-infection was investigated. BMDCs were found to be more permissive than monocytes to PRRSV infection; S. suis phagocytosis by PRRSV-infected BMDCs was found to be impaired, whereas no effect was found on bacterial intracellular survival. Transcription profile analysis, with a major focus on inflammatory genes, following S. suis infection, with and without pre-infection with PRRSV, was then performed. While PRRSV pre-infection had little effect on monocytes response to S. suis infection, a significant expression of several pro-inflammatory molecules was observed in BMDCs pre-infected with PRRSV after a subsequent infection with S. suis. While an additive effect could be observed for CCL4, CCL14, CCL20, and IL-15, a distinct synergistic up-regulatory effect was observed for IL-6, CCL5 and TNF-α after co-infection. This increased pro-inflammatory response by DCs could participate in the exacerbation of the disease observed during PRRSV and S. suis co-infection.

Chlamydial infection and spatial ascension of the female genital tract : a novel hybrid cellular automata and continuum mathematical model

Sexually transmitted chlamydial infection initially establishes in the endocervix in females, but if the infection ascends the genital tract, significant disease, including infertility, can result. Many of the mechanisms associated with chlamydial infection kinetics and disease ascension are unknown. We attempt to elucidate some of these processes by developing a novel mathematical model, using a cellular automata–partial differential equation model. We matched our model outputs to experimental data of chlamydial infection of the guinea-pig cervix and carried out sensitivity analyses to determine the relative influence of model parameters. We found that the rate of recruitment and action of innate immune cells to clear extracellular chlamydial particles and the rate of passive movement of chlamydial particles are the dominant factors in determining the early course of infection, magnitude of the peak chlamydial time course and the time of the peak. The rate of passive movement was found to be the most important factor in determining whether infection would ascend to the upper genital tract. This study highlights the importance of early innate immunity in the control of chlamydial infection and the significance of motility-diffusive properties and the adaptive immune response in the magnitude of infection and in its ascension.

Characterisation of mesenchymal cells from osteophytes in osteoarthritis

Osteophytes form through the process of chondroid metamorphosis of fibrous tissue followed by endochondral ossification. Osteophytes have been found to consist of three different mesenchymal tissue regions including endochondral bone formation within cartilage residues, intra-membranous bone formation within fibrous tissue and bone formation within bone marrow spaces. All these features provide evidence of mesenchymal stem cells (MSC) involvement in osteophyte formation; nevertheless, it remains to be characterised. MSC from numerous mesenchymal tissues have been isolated but bone marrow remains the “ideal” due to the ease of ex vivo expansion and multilineage potential. However, the bone marrow stroma has a relatively low number of MSC, something that necessitates the need for long-term culture and extensive population doublings in order to obtain a sufficient number of cells for therapeutic applications. MSC in vitro have limited proliferative capacity and extensive passaging compromises differentiation potential. To overcome this barrier, tissue derived MSC are of strong interest for extensive study and characterisation, with a focus on their potential application in therapeutic tissue regeneration. To date, no MSC type cell has been isolated from osteophyte tissue, despite this tissue exhibiting all the hallmark features of a regenerative tissue. Therefore, this study aimed to isolate and characterise cells from osteophyte tissues in relation to their phenotype, differentiation potential, immuno-modulatory properties, proliferation, cellular ageing, longevity and chondrogenesis in in vitro defect model in comparison to patient matched bone marrow stromal cells (bMSC). Osteophyte derived cells were isolated from osteophyte tissue samples collected during knee replacement surgery. These cells were characterised by the expression of cell surface antigens, differentiation potential into mesenchymal lineages, growth kinetics and modulation of allo-immune responses. Multipotential stem cells were identified from all osteophyte samples namely osteophyte derived mesenchymal stem cells (oMSC). Extensively expanded cell cultures (passage 4 and 9 respectively) were used to confirm cytogenetic stability and study signs of cellular aging, telomere length and telomerase activity. Cultured cells at passage 4 were used to determine 84 pathway focused stem cell related gene expression profile. Micro mass pellets were cultured in chondrogenic differentiation media for 21 days for phenotypic and chondrogenic related gene expression. Secondly, cell pellets differentiated overnight were placed into articular cartilage defects and cultured for further 21 days in control medium and chondrogenic medium to study chondrogenesis and cell behaviour. The surface antigen expression of oMSC was consistent with that of mesenchymal stem cells, such as lacking the haematopoietic and common leukocyte markers (CD34, CD45) while expressing those related to adhesion (CD29, CD166, CD44) and stem cells (CD90, CD105, CD73). The proliferation capacity of oMSC in culture was superior to that of bMSC, and they readily differentiated into tissues of the mesenchymal lineages. oMSC also demonstrated the ability to suppress allogeneic T-cell proliferation, which was associated with the expression of tryptophan degrading enzyme indoleamine 2,3 dioxygenase (IDO). Cellular aging was more prominent in late passage bMSC than in oMSC. oMSC had longer telomere length in late passages compared with bMSC, although there was no significant difference in telomere lengths in the early passages in either cell type. Telomerase activity was detectable only in early passage oMSC and not in bMSC. In osteophyte tissues telomerase positive cells were found to be located peri vascularly and were Stro-1 positive. Eighty-four pathway-focused genes were investigated and only five genes (APC, CCND2, GJB2, NCAM and BMP2) were differentially expressed between bMSC and oMSC. Chondrogenically induced micro mass pellets of oMSC showed higher staining intensity for proteoglycans, aggrecan and collagen II. Differential expression of chondrogenic related genes showed up regulation of Aggrecan and Sox 9 in oMSC and collagen II in bMSC. The in vitro defect models of oMSC in control medium showed rounded and aggregated cells staining positively for proteoglycan and presence of some extracellular matrix. In contrast, defects with bMSC showed fragmentation and loss of cells, fibroblast-like cell morphology staining positively for proteoglycans. For defects maintained in chondrogenic medium, rounded, aggregated and proteoglycan positive cells were found in both oMSC and bMSC cultures. Extracellular matrix and cellular integration into newly formed matrix was evident only in oMSC defects. For analysis of chondrocyte hypertrophy, strong expression of type X collagen could be noticed in the pellet cultures and transplanted bMSC. In summary, this study demonstrated that osteophyte derived cells had similar properties to mesenchymal stem cells in the expression of antigen phenotype, differential potential and suppression of allo-immune response. Furthermore, when compared to bMSC, oMSC maintained a higher proliferative capacity due to a retained level of telomerase activity in vitro, which may account for the relatively longer telomeres delaying growth arrest by replicative senescence compared with bMSC. oMSC behaviour in defects supported chondrogenesis which implies that cells derived from regenerative tissue can be an alternative source of stem cells and have a potential clinical application for therapeutic stem cell based tissue regeneration.