Immune cells must accurately interpret environmental signals to make robust cell fate decisions and control inflammatory signalling. Many signals (e.g. Tumor Necrosis Factor alpha (TNFα) or interferon gamma (IFNγ)) converge on just a few key signalling systems such as Nuclear Factor kappa B (NF-κB) or Signal Transducers and Activators of Transcription (STAT), which exhibit complex activation dynamics that control patterns of downstream gene expression. Often, seemingly identical cells show heterogeneous or random behaviour to a common stimulus. Therefore, a key question is how can immune cells coordinate inflammatory signalling in the presence of this noise. The NF-kB system dynamics were studied in response to rapidly changing inflammatory signals. It was shown that pulsed TNFa cytokine stimulations induced digital single-cell NF-kB responses, with only a fraction of cells able to respond to repeated pulses. These responses appeared to be reproducible in individual cells, but heterogeneous in the population. Mathematical models of the NF-kB signalling network suggested that single-cell responses were governed through a refractory state potentially encoded via 'extrinsic' noise in the levels of signalling molecules related to the TNFa signal transduction pathway. Such signal processing enabled robust and reproducible single cell responses and maintained acute tissue-level signalling, with fewer cells responding to shorter pulsing intervals. The NF-kB system is involved in effector cytokine propagation in response to pathogen infection. It was shown that in macrophages, the dose of TLR4 stimulation (mimicking the pathogen infection) was encoded in graded (yet heterogeneous) NF-kB dynamics in single cells. This resulted in analogue inflammatory gene expression patterns in the population. However, individual cells substantially differed in their ability to encode TLR4 signal and to regulate TNFa expression, which was explained by extrinsic noise in the NF-kB system. Quantitative mathematical modelling showed that tissue-level environment modulates heterogeneous single-cell TNFa outputs; by effectively removing it from circulation. This may determine the interaction distance between tissue-resident immune cells to enable propagation of cellular inflammation. Heterogeneity of single cell macrophage signalling was also observed in NF-kB and STAT1 system responses to a range of IFN stimulation doses. Although each system showed substantial variability between cells, their responses were surprisingly well correlated in individual cells. It was however apparent (based on gene expression studies) that individual cells may not be able to precisely discriminate different IFNg doses. Overall, this work suggests that heterogeneity in the NF-kB (and other) regulatory networks might be a part of an inherent design motif in the inflammatory response, which enables robust control of the tissue-level inflammatory response by preventing homogeneous and thus potentially harmful activation.

A previous thoracic aortic aneurysm time course study from our laboratory determined that ascending aortic dilation was significantly increased by day 5, and reached a plateau by day 28 of angiotensin II (AngII) infusion. We also found that mice had hemorrhage localized to the ascending aortas by day 5 of AngII infusion. The purpose of these studies was to provide mechanistic insight into the development of AngII-induced ascending aortic hemorrhage. Male C57BL/6 mice fed normal diet were subcutaneously infused with either AngII (1000 ng/kg/min) or saline for 5 days. To examine cellular heterogeneity, hemorrhaged ascending aortas were collected and sectioned serially for histological staining and immunostaining. I was unable to identify an entry point for blood into the media of the aortic root and ascending aorta. However, I found incomplete intimo-medial dissection near the hemorrhaged regions that may potentially be contiguous with the blood. To investigate infiltration of immune cells during AngII infusion, immunohistochemistry of hemorrhaged ascending aortas was performed. The numbers of macrophages and neutrophils in AngII-infused aortas were increased in both medial and adventitial areas when compared with saline-infused aortas. Therefore, infiltration of immune cells at the point of dissection is associated with aortic hemorrhage during AngII infusion.

La malària és encara avui en dia una malaltia que causa aproximadament un milió de morts a l'any a tot el món. La seva eradicació suposa un gran repte per a la humanitat i per a la comunitat científica, en particular. El cultiu in vitro del paràsit és essencial per al desenvolupament de nous medicaments. Els mètodes de cultiu actuals es basen en l'heurística i requereixen millores.En aquesta tesi es presenta una aproximació teòrica al procés d'infecció a eritròcits en cultius in vitro amb Plasmodium falciparum, un dels protozous paràsits causants de la malària. El treball està centrat en la construcció i avaluació de models d'una complexitat adequada per tractar els problemes específics detectats pels experts en l'àmbit, i inclou també la formulació d'algorismes de simulació i el disseny de protocols experimentals.Aquest tipus de treball requereix de la col·laboració multidisciplinària. La visió dels experts en malària es complementa amb la modelització i simulació, que permet la comprovació dels supòsits preestablerts, la comprensió de fenòmens observats i la millora dels mètodes de cultiu actuals. Així doncs, cal establir i desenvolupar eines que permetin crear, analitzar i compartir models amb grups que estudien la malària des d'altres perspectives. En aquesta tesi, s'ha optat per la modelització basada en l'individu (IbM) i orientada a la reproducció de múltiples patrons (PoM). El model s'ha formulat seguint l'ODD, un protocol estàndard en el camp de l'ecologia teòrica, que s'ha adaptat a la representació de comunitats microbianes.Els models basats en l'individu (IbMs) defineixen un conjunt de normes que regeixen el comportament de cada cèl·lula i les seves interaccions amb les altres cèl·lules i amb el seu entorn immediat. A partir d'aquestes regles, i tenint en compte una certa diversitat dins de la població i un cert grau d'aleatorietat en els processos individuals, els IbMs mostren explícitament el comportament emergent del sistema en conjunt. Complementàriament, s'han aplicat conceptes propis de la termodinàmica per tal d'entendrel'aparició de patrons macroscòpics a partir de l'estructura de la població (per exemple de la distribució de les fases d'infecció entre els glòbuls vermells infectats).Aquesta recerca ha comportat la la creació i aplicació del model i simulador INDISIM-RBC, que ha demostrat ser una bona eina per millorar la comprensió dels cultius estudiats. Es tracta d'un model mecanicista, basat en l'individu, que reprodueix quantitativament els patrons observats en cultius reals a diferents nivells de descripció, i que en prediu el comportament sota determinades condicions.Hem demostrat que INDISIM-RBC pot ser emprat per a estudiar en detall alguns aspectes del cultiu del paràsit causant de la malària que calia aclarir. Permet realitzar experiments virtuals i així impulsar noves línies de recerca i explorar noves tècniques de cultiu. En particular, INDISIM-RBC s'ha utilitzat per millorar els protocols experimentals actuals del cultius estàtics, definint la geometria òptima de l'hematòcrit i els protocols de subcultiu més adequats per als cultius continus.El treball realitzat en malària s'ha comparat amb la investigació duta a terme pel grup de recerca em relació amb d'altres comunitats microbianes. D'aquesta manera, podem estudiar les propietats emergents dels sistemes microbians en general en relació als efectes de la individualitat de la cèl·lula, la diversitat de les poblacions, l'heterogeneïtat en el medi, o el caràcter local de les interaccions, entre d'altres. Aquesta visió general proporciona eines conceptuals que poden ser emprades per refinar l'anàlisi dels processos d'infecció sota estudi. / Malaria is still a major burden that causes approximately one million deaths annually worldwide. Its eradication supposes a great challenge to the humanity and to the scientific community, in particular. In vitro cultivation of the parasite is essential for the development of new drugs. Current culture methods are based on heuristics and demand for specific improvements.The present thesis is a theoretical approach to in vitro cultivation of the protozoan parasite Plasmodium falciparum infecting human red blood cells. It mainly focuses on the process of building a model of appropriate complexity to deal with the specific demands above mentioned, but it also includes the formulation and implementation of algorithms, and the design and execution of experimental trials.This kind of work requires multidisciplinary collaboration: the insight of the experts in malaria research is complemented with modeling and simulation, which allows for checking settled assumptions, increasing the understanding on the system and improving the current culturing methods.The use of tools for building, analyzing and sharing models is an imperative to this end. In this thesis, Pattern-oriented Modeling (PoM) has been adopted as the most appropriate way for raising of models and the ODD protocol (Objectives, Design Concepts and Details) has been proposed as the standard tool for communicating them.Individual-based Modeling (IbM) has been used to tackle malaria culture systems. IbMs define a set of rules governing each cell, its interactions with others and with its immediate surroundings. From this set of rules, and taking into account diversity within the population and a certain degree of randomness in the individual processes, IbMs explicitly show the emerging behavior of the system as a whole. Methods from statistical thermodynamics have been applied to understand the emergence of macroscopic patterns from the population structure (e.g. distribution of infection stages among infected red blood cells).The research resulted in the development of the model and simulator INDISIM-RBC, which has proved to be a good tool to improve understanding of the cultures under study. It is a mechanistically rich individual-based model and it quantitatively reproduces and predicts several patterns observed in real cultures at different levels of description.We demonstrated that INDISIM-RBC can be used to study in detail several aspects of malaria cultivation that remained unclear, as well as to perform virtual experiments. Consequently, it can be used to open novel lines of research and to examine potential experimental techniques. INDISIM-RBC has also been used to improve the current experimental culturing protocols in static cultivation by obtaining the optimal geometry of the hematocrit layer and subcultivation periods in the continuous cultures.This study on malaria has been compared to the research carried out by the group regarding other microbial communities. Thereby studying general emerging properties of microbial systems in general, with regard to the effect of cell individuality, heterogeneity and diversity, the local nature of interactions; and biological and spatial complexity. In doing so, the acquired holistic view has been used to develop tools that allow for a better characterization and study of the infection process, in particular.

In mammalian hearts, initiation of the heartbeat occurs in a region of specialised pacemaker cells known as the sinoatrial node (SAN). The SAN is a highly complex spatially distributed structure which displays considerable cellular heterogeneity and is subject to complex electrotonic interactions with the surrounding atrial tissue. In this study, biophysically detailed ionic models of central and peripheral SAN pacemaker cells are described. These models are able to accurately reproduce experimental recordings of the membrane potential from central and peripheral SAN tissue. These models are used to investigate frequency synchronisation of electrically coupled cardiac pacemaker cells. Based on simulation results presented, it is proposed that cellular heterogeneity in the SAN plays an important role in achieving rhythm coordination and possibly contributes to the efficient activation of the surrounding atrial myocardium. This represents an important, previously unexplored, mechanism underlying pacemaker synchronisation and cardiac activation in vivo. A spatial-gradient model of action potential heterogeneity within the SAN is then formulated using a large-scale least squares optimisation technique. This model accurately reproduces the smooth spatial variation in action potential characteristics observed in the SAN. One and two dimensional models of the intact SAN are then formulated and three proposed models of SAN heterogeneity are investigated: 1) the discrete-region model, in which the SAN consists of a compact central region surrounded by a region of transitional pacemaker cells, 2) the gradient model, in which cells of the SAN exhibit a smooth variation in properties from the centre to the periphery of the SAN, and 3) the mosaic model, in which SAN and atrial cells are scattered throughout the SAN region with the proportion of atrial cells increasing towards the periphery. Simulation results suggest that the gradient model achieves frequency entrainment more easily than the other models of SAN heterogeneity. The gradient model also reproduces action potential waveshapes and a site of earliest activation consistent with experimental observations in the intact SAN. It is therefore proposed that the gradient model of SAN heterogeneity represents the most plausible model of SAN organisation.

Notre équipe a développé un protocole d’ingénierie de l’organe dentaire basé sur la biomimétique et l’utilisation de cellules dentaires embryonnaires dissociées. La recherche de ressources cellulaires permettant d’éviter le recours aux cellules embryonnaire reste un défi majeur, et nécessite une meilleure connaissance des paramètres limitants. Nous avons testé les potentialités odontogènes de lignées cellulaires, dentaires ou non, embryonnaires ou adultes. Le développement dentaire étant contrôlé par des interactions réciproques entre ectomésenchyme dérivé des crêtes neurales et épithélium, ces cellules ont été réassociés à un épithélium dentaire compétent. Nous avons recherché la formation de dents in vitro et/ou après implantation chez la souris adulte et étudié un certain nombre de paramètres biologiques et techniques. Ainsi, nous avons étudié l’impact de l’âge, de la mise en culture, et de l’hétérogénéité cellulaire sur les potentialités odontogènes des cellules mésenchymateuses. Nos résultats montrent que le potentiel odontogène des différentes lignées mésenchymateuses testées pouvait être lié à l’âge des cellules et qu’il est perdu lorsque les cellules mésenchymateuses sont cultivées avant d’être ré-associées. Ceci pouvant s’expliquer par un changement phénotypique, nous avons testé un certain nombre de gènes essentiels au développement dentaire, et suivi l’expression de marqueurs de surface. Les changements observés peuvent être liés à une sélection cellulaire in vitro pouvant conduire à des modifications de l’hétérogénéité des cellules en monocouche. / Our laboratory has developed a protocol for tooth organ engineering, based on biomimetic and the use of dissociated embryonic dental cells. Searching for mesenchymal cell sources avoiding the use of embryonic cells still remains a major challenge. We tested the odontogenic potential of several cell lines, dental or not, embryonic or adult. These cells were re-associated with a competent intact dental epithelium. We searched for tooth formation in vitro and/or after implantation in adult mice and studied different biological and technical parameters. For this purpose we analyzed the effects of the age, of a pre-culture step, and of the cellular heterogeneity on the odontogenetic potential of mesenchymal cells. To test the heterogeneity, we compared the patterns of expression of cell surface markers in cultured and implanted re-association with those observed during tooth development..Our results show that the absence of odontogenetic potential in the different cell lines that were tested, in part depends on the age of cells and that it is lost when mesenchymal cells are cultured in monolayer before their re-association. This could be explained by phenotypical changes as shown by testing several genes involved in tooth development, and tracing cell surface markers expression. Changes were observed, wich could be related to a cell selection in vitro, leading to variations in cellular heterogeneity. Indeed, pulpal cellular heterogeneity shows specific patterns of expression, that are space-time defined during tooth development, and is controlled by epithelial-mesenchymal interactions.

Vascular smooth muscle cell (VSMC) accumulation is a hallmark of atherosclerosis and vascular injury. However, fundamental aspects of proliferation and the phenotypic changes within individual VSMCs, which underlie vascular disease remain unresolved. In particular, it is not known if all VSMCs proliferate and display plasticity, or whether individual cells can switch to multiple phenotypes. To assess whether proliferation and plasticity in disease is a general characteristic of VSMCs or a feature of a subset of cells, multi-colour lineage labelling is used to demonstrate that VSMCs in injury-induced neointimal lesions and in atherosclerotic plaques are oligo-clonal, derived from few expanding cells, within mice. Lineage tracing also revealed that the progeny of individual VSMCs contribute to both alpha Smooth muscle actin (aSma)-positive fibrous cap and Mac-3-expressing macrophage-like plaque core cells. Co-staining for phenotypic markers further identified a double-positive aSma+ Mac3+ cell population, which is specific to VSMC-derived plaque cells. In contrast, VSMC-derived cells generating the neointima after vascular injury generally retained expression of VSMC markers and upregulation of Mac3 was less pronounced. Monochromatic regions in atherosclerotic plaques and injury-induced neointima did not contain VSMC-derived cells expressing a different fluorescent reporter protein, suggesting that proliferation-independent VSMC migration does not make a major contribution to VSMC accumulation in vascular disease. Similarly, VSMC proliferation was examined in an Angiotensin II perfusion model of aortic aneurysm in mice, oligo-clonal proliferation was observed in remodelling regions of the vasculature, however phenotypic changes were observed in a large proportion of VSMCs, suggesting that the majority of VSMCs have some potential to modulate their phenotype. To understand the mechanisms behind the inherent VSMC heterogeneity and observed functionality, the single cell transcriptomic techniques Smart-seq2 and the Chromium 10X system were optimized for use on VSMCs. The work within this thesis suggests that extensive proliferation of a low proportion of highly plastic VSMCs results in the observed VSMC accumulation after injury, and the atherosclerotic and aortic aneurysm models of cardiovascular disease.

Der TGFβ-Signalweg ist ein multifunktionales System, das zelluläre Prozesse reguliert, die von Proliferation und Migration bis zu Differenzierung und Zelltod reichen. Nach Ligandenbindung und Rezeptoraktivierung translozieren SMAD-Proteine zum Zellkern und induzieren die Expression zahlreicher Zielgene. Während viele Komponenten des TGFβ-Signalweges identifiziert wurden, verstehen wir noch nicht genau, wie die Aktivierung des Signalwegs in verschiedene zelluläre Antworten übersetzt wird. Da die zelluläre Antwort auf einen gegebenen Stimulus oft sogar in genetisch identischen Zellen variiert, konzentrierte ich mich auf die Messung der Signalwegaktivität auf der Einzelzellebene. Durch die Kombination fluoreszierender Reporterzelllinien mit Zeitraffer-Lebendzellmikroskopie und automatisierter Bildanalyse beobachtete ich die zytoplasmatische und nukleäre Translokation von SMADs mit hoher zeitlicher und räumlicher Auflösung in Hunderten einzelner Zellen. Unsere Experimente zeigten, dass die Signalwegaktivität in eine erste synchrone Phase der SMAD-Translokation, gefolgt von einer Adaption und einer zweiten Signalphase mit hoher Variabilität in Stärke und Dauer der nuklearen Akkumulation unterteilt werden kann. Darüber hinaus beobachtete ich, dass Zellen, die aufgrund ihrer dynamischen Eigenschaften in Subpopulationen gruppiert sind, unterschiedliche phänotypische Reaktionen zeigen. Ich war nun daran interessiert, die Netzwerkinteraktionen zu identifizieren, die diese Dynamiken formen und fokussierte mich auf den Crosstalk mit nicht-kanonischen Komponenten des TGFβ-Signalweges. Ich konnte zeigen, dass die Hemmung der MAP Kinasen p38 und ERK die zweite Signalphase spezifisch aufhebt. Diese dynamische Remodellierung führt zu Veränderungen in der Zielgenexpression und den Zellschicksalen. Dies wird zu einem tieferen Verständnis der molekularen Netzwerke führen, die die TGFβ-Signaltransduktion regulieren und Möglichkeiten eröffnen, es in erkrankten Zellen zu modulieren. / The TGFβ pathway is a multi-functional signaling system regulating cellular processes ranging from proliferation and migration to differentiation and cell death. Upon ligand binding and receptor activation, SMAD proteins translocate to the nucleus and induce expression of numerous target genes. While many components of the TGFβ pathway have been identified, we are still challenged to understand how pathway activation is translated into distinct cellular responses. As the cellular response to a given stimulus often varies even in genetically identical cells, I focused on measuring pathway activity on the single cell level. By combining fluorescent reporter cell lines with time-lapse live-cell microscopy and automated image analysis, I monitored the cytoplasmic to nuclear translocation of SMADs with high temporal and spatial resolution in hundreds of individual cells. Our experiments demonstrated that pathway activity can be divided into a first synchronous phase of SMAD translocation, followed by adaptation and a second signaling phase with high variability in the extent and duration of nuclear accumulation. Furthermore, I observed that cells clustered into subpopulations according to their dynamic features show different phenotypic responses. I was interested in identifying the network interactions that shape these dynamics and focus on crosstalk with non-canonical components of the TGFβ pathway. I could show that inhibition of the MAP kinases p38 and ERK specifically abrogates the second signaling phase. This dynamic remodeling led to changes in target gene expression and cell fate decisions. I explored the molecular mechanisms underlying this interaction of the canonical and non-canonical pathways. This will provide a deeper understanding of the molecular networks regulating TGFβ signaling and open opportunities to modulate it in diseased cells.

Das normale Kolonepithel weist eine strenge Zellhierarchie auf, die aus bekannten Zelltypen besteht. Bei Darmkrebs (CRC) ist die Struktur weniger konserviert und nicht gut verstanden. Krebsauslösende Mutationen können die Prävalenz von Zelltypen verändern, und Zellen können sich auch dedifferenzieren, um einer gezielten Krebstherapie zu entgehen. Mein Ziel ist es, die Existenz heterogener Zelltypen in Organoiden zu bestätigen und Signalnetzwerke in CRC zu untersuchen, indem ich mit pharmakologischen Eingriffen spezifische Signalwege inhibiere, die Zellhierarchien im normalen Darm kontrollieren. Strategisch ausgewählte Medikamente wurden eingesetzt, um Knotenpunkte in verschiedenen Signalwegen zu hemmen, die für das Fortschreiten von Darmkrebs relevant sind. Ich untersuchte, ob die Inhibition von Signalwegen die Zusammensetzung der Zelltypen und den Differenzierungszustand verändert oder welche Kombinationen von Inhibitoren Plastizität oder Apoptose auslösen könnten. Von Patienten stammende Organoide mit verschiedenen onkogenen Treibermutationen wurden kultiviert und 48 Stunden lang mit einer Reihe von Inhibitoren und Inhibitorkombinationen behandelt. Diese Organoide wurden hauptsächlich auf zwei Ebenen untersucht: durch scRNA seq zur Ermittlung ihres Transkriptoms und durch CyTOF, das die Proteinhäufigkeit pro Zelle misst, um die Aktivität von Signaltransduktionskaskaden zu beurteilen. Beide Methoden wurden eingesetzt, um die Heterogenität des CRC zu quantifizieren. Ich konnte feststellen, dass sich Organoide mit denselben Treibermutationen ähnlicher verhalten und dass die molekularen Grundlagen der verschiedenen Linien Unterschiede im Therapieerfolg bedingen. Heterogene Transkriptome und Proteinexpression wurden durch einen Differenzierungsgradienten beeinflusst und konnten durch die Zugabe von Inhibitoren verändert werden. Die MAPK-Aktivität folgt diesem Differenzierungsgradienten und eine MAPK-Inhibition verringerte die Zellheterogenität und führte zu Plastizität. Darüber hinaus stellte ich fest, dass ein Teil der Zellen in Apoptose geht und die verbleibenden Zellen einen nicht-proliferativen Stammzellzustand annehmen, der es den Zellen ermöglicht, sich nach Aussetzung der Behandlung zu erholen. Es wurden in silico und in vitro Analysen durchgeführt, um neuartige Inhibitorkombinationen zur Maximierung der Apoptose in CRC-Organoiden zu finden, um die Entstehung therapieresistenter Subpopulationen weiter zu reduzieren. Wirksame Behandlungskombinationen bleiben jedoch zelllinienabhängig. Durch die getrennte Analyse des Zelldifferenzierungszustands und des Zellsignalisierungszustands habe ich dazu beigetragen zu verstehen, wie Tumorzellen einer gezielten Therapie durch nicht-genetische Resistenzmechanismen entgehen können. Die MAPK-Inhibition zur Verringerung der Zellheterogenität in Kombination mit anderen Inhibitoren könnte in Zukunft zur Optimierung des Therapieerfolgs eingesetzt werden. / Normal colon epithelium has a strict cell hierarchy consisting of well-known cell types. In colorectal cancer (CRC) the structure is less conserved and poorly understood. Cancer driver mutations may modulate the prevalence of cell types, and cells may also dedifferentiate to overcome targeted cancer therapy. My aim is to confirm the existence of heterogeneous cell types in organoids and investigate signaling networks in CRC by targeting specific signaling pathways with pharmacological intervention, which control cell hierarchies in the normal intestine. Strategically selected drugs were used to inhibit nodes in different signaling pathways relevant to the progression of CRC. I explored whether signaling inhibition changes cell type composition and differentiation state, or which inhibitor combinations might induce plasticity or apoptosis. Patient-derived organoids with different oncogenic diver mutations were cultured and treated with a panel of inhibitors and inhibitor combinations for 48 hours. These organoids were mainly examined on two levels: by scRNA seq to assess their transcriptome and by CyTOF, which measures protein abundance to assess the activity of pathways. Both methods were used to quantify CRC heterogeneity. I was able to see that organoids with the same driver mutations behave more similarly and that the molecular underpinnings of the different lines drive differences in therapy response. Heterogeneous transcriptomes and protein expression were affected by a differentiation gradient and could be altered by inhibitor addition. MAPK activity was graded along this differentiation gradient, and MAPK inhibition reduced cell heterogeneity and induced plasticity. Additionally, I found that a fraction of cells undergo apoptosis, and the remaining cells adopt a non-proliferative stem cell state, which allows cells to recover after suspension of treatment. \textit{In silico} and \textit{in vitro} analyses were performed to find novel inhibitor combinations to maximize apoptosis in CRC organoids to further reduce the emergence of therapy-resistant subpopulations. However, effective treatment combinations remain cell-line dependent. By separately analyzing cell differentiation state and cell signaling state I contributed to our understanding of how tumor cells can evade targeted therapy by non-genetic resistance mechanisms. Using MAPK inhibition to reduce cell heterogeneity in combination with other inhibitors may be used in the future to optimize therapy success.