Endothelialized Bronchioalveolar Lung Organoids Model Endothelial Cell Responses to Injury.

Organoid 3D systems are powerful platforms to study development and disease. Recently, the complexity of lung organoid models derived from adult mouse and human stem cells has increased substantially in terms of cellular composition and structural complexity. However, a murine lung organoid system with a clear integrated endothelial compartment is still missing. Here, we describe a novel method that adds another level of intricacy to our published bronchioalveolar lung organoid (BALO) model by microinjection of FACS-sorted lung endothelial cells (ECs) into differentiated organoid cultures. Before microinjection, ECs obtained from the lung homogenate (LH) of young mice expressed typical ECs markers such as CD31 and vascular endothelial (VE)-Cadherin and showed tube formation capacity. Following microinjection, ECs surrounded BALO´s alveolar-like compartment aligning with both alveolar epithelial cells type I (AECI) and type II (AECII), as demonstrated by confocal and electron microscopy. Notably, expression of Car4 and Aplnr was as well detected, suggesting presence of EC microvascular phenotypes in the cultured ECs. Moreover, upon epithelial cell injury by lipopolysaccharides (LPS) and influenza A virus (IV), endothelialized BALO (eBALO) released proinflammatory cytokines leading to the upregulation of the intercellular adhesion molecule 1 (ICAM-1) in ECs. In summary, we characterized for the first time a organoid model that incorporates ECs into the alveolar structures of lung organoids, not only increasing our previous model ́s cellular and structural complexity but also providing a suitable niche to model lung endothelium responses to injury ex vivo.

Design patterns for the construction of computational biological models.

Computational biological models have proven to be an invaluable tool for understanding and predicting the behaviour of many biological systems. While it may not be too challenging for experienced researchers to construct such models from scratch, it is not a straightforward task for early stage researchers. Design patterns are well-known techniques widely applied in software engineering as they provide a set of typical solutions to common problems in software design. In this paper, we collect and discuss common patterns that are usually used during the construction and execution of computational biological models. We adopt Petri nets as a modelling language to provide a visual illustration of each pattern; however, the ideas presented in this paper can also be implemented using other modelling formalisms. We provide two case studies for illustration purposes and show how these models can be built up from the presented smaller modules. We hope that the ideas discussed in this paper will help many researchers in building their own future models.

GLI1+ Cells Contribute to Vascular Remodeling in Pulmonary Hypertension.

BACKGROUND: The precise origin of newly formed ACTA2+ (alpha smooth muscle actin-positive) cells appearing in nonmuscularized vessels in the context of pulmonary hypertension is still debatable although it is believed that they predominantly derive from preexisting vascular smooth muscle cells (VSMCs). METHODS: Gli1Cre-ERT2; tdTomatoflox mice were used to lineage trace GLI1+ (glioma-associated oncogene homolog 1-positive) cells in the context of pulmonary hypertension using 2 independent models of vascular remodeling and reverse remodeling: hypoxia and cigarette smoke exposure. Hemodynamic measurements, right ventricular hypertrophy assessment, flow cytometry, and histological analysis of thick lung sections followed by state-of-the-art 3-dimensional reconstruction and quantification using Imaris software were used to investigate the contribution of GLI1+ cells to neomuscularization of the pulmonary vasculature. RESULTS: The data show that GLI1+ cells are abundant around distal, nonmuscularized vessels during steady state, and this lineage contributes to around 50% of newly formed ACTA2+ cells around these normally nonmuscularized vessels. During reverse remodeling, cells derived from the GLI1+ lineage are largely cleared in parallel to the reversal of muscularization. Partial ablation of GLI1+ cells greatly prevented vascular remodeling in response to hypoxia and attenuated the increase in right ventricular systolic pressure and right heart hypertrophy. Single-cell RNA sequencing on sorted lineage-labeled GLI1+ cells revealed an Acta2high fraction of cells with pathways in cancer and MAPK (mitogen-activated protein kinase) signaling as potential players in reprogramming these cells during vascular remodeling. Analysis of human lung-derived material suggests that GLI1 signaling is overactivated in both group 1 and group 3 pulmonary hypertension and can promote proliferation and myogenic differentiation. CONCLUSIONS: Our data highlight GLI1+ cells as an alternative cellular source of VSMCs in pulmonary hypertension and suggest that these cells and the associated signaling pathways represent an important therapeutic target for further studies.

Alveolar macrophage-expressed Plet1 is a driver of lung epithelial repair after viral pneumonia.

Influenza A virus (IAV) infection mobilizes bone marrow-derived macrophages (BMDM) that gradually undergo transition to tissue-resident alveolar macrophages (TR-AM) in the inflamed lung. Combining high-dimensional single-cell transcriptomics with complex lung organoid modeling, in vivo adoptive cell transfer, and BMDM-specific gene targeting, we found that transitioning ("regenerative") BMDM and TR-AM highly express Placenta-expressed transcript 1 (Plet1). We reveal that Plet1 is released from alveolar macrophages, and acts as important mediator of macrophage-epithelial cross-talk during lung repair by inducing proliferation of alveolar epithelial cells and re-sealing of the epithelial barrier. Intratracheal administration of recombinant Plet1 early in the disease course attenuated viral lung injury and rescued mice from otherwise fatal disease, highlighting its therapeutic potential.

Protocol for biomodel engineering of unilevel to multilevel biological models using colored Petri nets.

Biological systems inherently span multiple levels, which can pose challenges in spatial representation for modelers. We present a protocol that utilizes colored Petri nets to construct and analyze biological models of systems, encompassing both unilevel and multilevel scenarios. We detail a modeling workflow exploiting the PetriNuts platform comprising a set of tools linked together via common file formats. We describe steps for modeling preparation, component-level modeling and analysis, followed by system-level modeling and analysis, and model use.

FGF10 Triggers De Novo Alveologenesis in a Bronchopulmonary Dysplasia Model: Impact on Resident Mesenchymal Niche Cells.

Bronchopulmonary dysplasia (BPD) is a neonatal lung disease developing in premature babies characterized by arrested alveologenesis and associated with decreased Fibroblast growth factor 10 (FGF10) expression. One-week hyperoxia (HYX) exposure of newborn mice leads to a permanent arrest in alveologenesis. To test the role of Fgf10 signaling to promote de novo alveologenesis following hyperoxia, we used transgenic mice allowing inducible expression of Fgf10 and recombinant FGF10 (rFGF10) protein delivered intraperitoneally. We carried out morphometry analysis, and IF on day 45. Alveolospheres assays were performed co-culturing AT2s from normoxia (NOX) with FACS-isolated Sca1Pos resident mesenchymal cells (rMC) from animals exposed to NOX, HYX-PBS, or HYX-FGF10. scRNAseq between rMC-Sca1Pos isolated from NOX and HYX-PBS was also carried out. Transgenic overexpression of Fgf10 and rFGF10 administration rescued the alveologenesis defects following HYX. Alveolosphere assays indicate that the activity of rMC-Sca1Pos is negatively impacted by HYX and partially rescued by rFGF10 treatment. Analysis by IF demonstrates a significant impact of rFGF10 on the activity of resident mesenchymal cells. scRNAseq results identified clusters expressing Fgf10, Fgf7, Pdgfra, and Axin2, which could represent the rMC niche cells for the AT2 stem cells. In conclusion, we demonstrate that rFGF10 administration is able to induce de novo alveologenesis in a BPD mouse model and identified subpopulations of rMC-Sca1Pos niche cells potentially representing its cellular target.

Hybrid modelling of biological systems: current progress and future prospects.

Integrated modelling of biological systems is becoming a necessity for constructing models containing the major biochemical processes of such systems in order to obtain a holistic understanding of their dynamics and to elucidate emergent behaviours. Hybrid modelling methods are crucial to achieve integrated modelling of biological systems. This paper reviews currently popular hybrid modelling methods, developed for systems biology, mainly revealing why they are proposed, how they are formed from single modelling formalisms and how to simulate them. By doing this, we identify future research requirements regarding hybrid approaches for further promoting integrated modelling of biological systems.

Coloured fuzzy Petri nets for modelling and analysing membrane systems.

Membrane systems are a very powerful computational modelling formalism inspired by the internal organisation of living cells. Modelling of membrane systems is challenged by composing many structurally similar components, which may result in very large models. Furthermore, some components may suffer from a lack of precise kinetic parameters. FPNC (coloured fuzzy Petri nets) combine coloured Petri nets with fuzzy kinetic parameters, and thus offer an approach to address these challenges. In this paper, we use FPNC to model and simulate membrane systems which are enriched by fuzzy kinetic parameters. We also introduce a methodology and workflow utilising FPNC for modelling and simulating general biological systems which have to cope with incomplete knowledge of their kinetic data.

A Petri nets-based framework for whole-cell modeling.

Whole-cell modeling aims to incorporate all main genes and processes, and their interactions of a cell in one model. Whole-cell modeling has been regarded as the central aim of systems biology but also as a grand challenge, which plays essential roles in current and future systems biology. In this paper, we analyze whole-cell modeling challenges and requirements and classify them into three aspects (or dimensions): heterogeneous biochemical networks, uncertainties in components, and representation of cell structure. We then explore how to use different Petri net classes to address different aspects of whole-cell modeling requirements. Based on these analyses, we present a Petri nets-based framework for whole-cell modeling, which not only addresses many whole-cell modeling requirements, but also offers a graphical, modular, and hierarchical modeling tool. We think this framework can offer a feasible modeling approach for whole-cell model construction.

Characterization in mice of the resident mesenchymal niche maintaining AT2 stem cell proliferation in homeostasis and disease.

Resident mesenchymal cells (rMCs defined as Cd31Neg Cd45Neg EpcamNeg ) control the proliferation and differentiation of alveolar epithelial type 2 (AT2) stem cells in vitro. The identity of these rMCs is still elusive. Among them, Axin2Pos mesenchymal alveolar niche cells (MANCs), which are expressing Fgf7, have been previously described. We propose that an additional population of rMCs, expressing Fgf10 (called rMC-Sca1Pos Fgf10Pos ) are equally important to maintain AT2 stem cell proliferation. The alveolosphere model, based on the AT2-rMC co-culture in growth factor-reduced Matrigel, was used to test the efficiency of different rMC subpopulations isolated by FACS from adult murine lung to sustain the proliferation and differentiation of AT2 stem cells. We demonstrate that rMC-Sca1Pos Fgf10Pos cells are efficient to promote the proliferation and differentiation of AT2 stem cells. Co-staining of adult lung for Fgf10 mRNA and Sftpc protein respectively, indicate that 28% of Fgf10Pos cells are located close to AT2 cells. Co-ISH for Fgf7 and Fgf10 indicate that these two populations do not significantly overlap. Gene arrays comparing rMC-Sca1Pos Axin2Pos and rMC-Sca1Pos Fgf10Pos support that these two cell subsets express differential markers. In addition, rMC function is decreased in obese ob/ob mutant compared to WT mice with a much stronger loss of function in males compared to females. In conclusion, rMC-Sca1Pos Fgf10Pos cells play important role in supporting AT2 stem cells proliferation and differentiation. This result sheds a new light on the subpopulations of rMCs contributing to the AT2 stem cell niche in homeostasis and in the context of pre-existing metabolic diseases.

Identification of a novel subset of alveolar type 2 cells enriched in PD-L1 and expanded following pneumonectomy.

Alveolar type 2 (AT2) cells are heterogeneous cells, with specialised AT2 subpopulations within this lineage exhibiting stem cell properties. However, the existence of quiescent, immature cells within the AT2 lineage that are activated during lung regeneration is unknown.SftpcCreERT2/+;tdTomatoflox/flox mice were used for the labelling of AT2 cells and labelled subpopulations were analysed by flow cytometry, quantitative PCR, assay for transposase-accessible chromatin using sequencing (ATAC-seq), gene arrays, pneumonectomy and culture of precision-cut lung slices. Single-cell RNA-sequencing (scRNA-seq) data from human lungs were analysed.In mice, we detected two distinct AT2 subpopulations, with low tdTomato level (TomLow) and high tdTomato level (TomHigh). TomLow cells express lower levels of the AT2 differentiation markers Fgfr2b and Etv5, while TomHigh, as bona fide mature AT2 cells, show higher levels of Sftpc, Sftpb, Sftpa1, Fgfr2b and Etv5 expression. ATAC-seq analysis indicates that TomLow and TomHigh cells constitute two distinct cell populations, with specific silencing of Sftpc, Rosa26 and cell cycle gene loci in the TomLow population. Upon pneumonectomy, the number of TomLow but not TomHigh cells increases and TomLow cells show upregulated expression of Fgfr2b, Etv5, Sftpc, Ccnd1 and Ccnd2 compared to Sham. TomLow cells overexpress programmed cell death 1 ligand 1 (PD-L1), an immune inhibitory membrane receptor ligand, which is used by flow cytometry to differentially isolate these two subpopulations. In the human lung, data mining of a recent scRNA-seq AT2 data set demonstrates the existence of a PD-L1 Pos population. Therefore, we have identified a novel population of AT2 quiescent, immature progenitor cells in mouse that expand upon pneumonectomy and we have provided evidence for the existence of such cells in human.

Hybrid modelling of biological systems using fuzzy continuous Petri nets.

Integrated modelling of biological systems is challenged by composing components with sufficient kinetic data and components with insufficient kinetic data or components built only using experts' experience and knowledge. Fuzzy continuous Petri nets (FCPNs) combine continuous Petri nets with fuzzy inference systems, and thus offer an hybrid uncertain/certain approach to integrated modelling of such biological systems with uncertainties. In this paper, we give a formal definition and a corresponding simulation algorithm of FCPNs, and briefly introduce the FCPN tool that we have developed for implementing FCPNs. We then present a methodology and workflow utilizing FCPNs to achieve hybrid (uncertain/certain) modelling of biological systems illustrated with a case study of the Mercaptopurine metabolic pathway. We hope this research will promote the wider application of FCPNs and address the uncertain/certain integrated modelling challenge in the systems biology area.

Multilineage murine stem cells generate complex organoids to model distal lung development and disease.

Organoids derived from mouse and human stem cells have recently emerged as a powerful tool to study organ development and disease. We here established a three-dimensional (3D) murine bronchioalveolar lung organoid (BALO) model that allows clonal expansion and self-organization of FACS-sorted bronchioalveolar stem cells (BASCs) upon co-culture with lung-resident mesenchymal cells. BALOs yield a highly branched 3D structure within 21 days of culture, mimicking the cellular composition of the bronchioalveolar compartment as defined by single-cell RNA sequencing and fluorescence as well as electron microscopic phenotyping. Additionally, BALOs support engraftment and maintenance of the cellular phenotype of injected tissue-resident macrophages. We also demonstrate that BALOs recapitulate lung developmental defects after knockdown of a critical regulatory gene, and permit modeling of viral infection. We conclude that the BALO model enables reconstruction of the epithelial-mesenchymal-myeloid unit of the distal lung, thereby opening numerous new avenues to study lung development, infection, and regenerative processes in vitro.

Fuzzy Petri nets for modelling of uncertain biological systems.

The modelling of biological systems is accompanied with epistemic uncertainties that range from structural uncertainty to parametric uncertainty due to such limitations as insufficient understanding of the underlying mechanism and incomplete measurement data of a system. Fuzzy logic approaches such as fuzzy Petri nets (FPNs) are effective in addressing these issues. In this paper, we review FPNs that have been used for modelling uncertain biological systems, which we classify in three categories: basic fuzzy Petri nets, fuzzy quantitative Petri nets and Petri nets with fuzzy kinetic parameters. For each category of these FPNs, we summarize its modelling capabilities and current applications, discuss its merits and drawbacks and give suggestions for further research. This understanding on how to use FPNs for modelling uncertain biological systems will assist readers in selecting appropriate FPN classes for specific modelling circumstances. This review may also promote the extensive research and application of FPNs in the systems biology area.

Regulatory Dynamics of Cell Differentiation Revealed by True Time Series From Multinucleate Single Cells.

Dynamics of cell fate decisions are commonly investigated by inferring temporal sequences of gene expression states by assembling snapshots of individual cells where each cell is measured once. Ordering cells according to minimal differences in expression patterns and assuming that differentiation occurs by a sequence of irreversible steps, yields unidirectional, eventually branching Markov chains with a single source node. In an alternative approach, we used multi-nucleate cells to follow gene expression taking true time series. Assembling state machines, each made from single-cell trajectories, gives a network of highly structured Markov chains of states with different source and sink nodes including cycles, revealing essential information on the dynamics of regulatory events. We argue that the obtained networks depict aspects of the Waddington landscape of cell differentiation and characterize them as reachability graphs that provide the basis for the reconstruction of the underlying gene regulatory network.

Agent-based modeling and bifurcation analysis reveal mechanisms of macrophage polarization and phenotype pattern distribution.

Macrophages play a key role in tissue regeneration by polarizing to different destinies and generating various phenotypes. Recognizing the underlying mechanisms is critical in designing therapeutic procedures targeting macrophage fate determination. Here, to investigate the macrophage polarization, a nonlinear mathematical model is proposed in which the effect of IL4, IFNγ and LPS, as external stimuli, on STAT1, STAT6, and NFκB is studied using bifurcation analysis. The existence of saddle-node bifurcations in these internal key regulators allows different combinations of steady state levels which are attributable to different fates. Therefore, we propose dynamic bifurcation as a crucial built-in mechanism of macrophage polarization. Next, in order to investigate the polarization of a population of macrophages, bifurcation analysis is employed aligned with agent-based approach and a two-layer model is proposed in which the information from single cells is exploited to model the behavior in tissue level. Also, in this model, a partial differential equation describes the diffusion of secreted cytokines in the medium. Finally, the model was validated against a set of experimental data. Taken together, we have here developed a cell and tissue level model of macrophage polarization behavior which can be used for designing therapeutic interventions.

Spatial quorum sensing modelling using coloured hybrid Petri nets and simulative model checking.

BACKGROUND: Quorum sensing drives biofilm formation in bacteria in order to ensure that biofilm formation only occurs when colonies are of a sufficient size and density. This spatial behaviour is achieved by the broadcast communication of an autoinducer in a diffusion scenario. This is of interest, for example, when considering the role of gut microbiota in gut health. This behaviour occurs within the context of the four phases of bacterial growth, specifically in the exponential stage (phase 2) for autoinducer production and the stationary stage (phase 3) for biofilm formation. RESULTS: We have used coloured hybrid Petri nets to step-wise develop a flexible computational model for E.coli biofilm formation driven by Autoinducer 2 (AI-2) which is easy to configure for different notions of space. The model describes the essential components of gene transcription, signal transduction, extra and intra cellular transport, as well as the two-phase nature of the system. We build on a previously published non-spatial stochastic Petri net model of AI-2 production, keeping the assumptions of a limited nutritional environment, and our spatial hybrid Petri net model of biofilm formation, first presented at the NETTAB 2017 workshop. First we consider the two models separately without space, and then combined, and finally we add space. We describe in detail our step-wise model development and validation. Our simulation results support the expected behaviour that biofilm formation is increased in areas of higher bacterial colony size and density. Our analysis techniques include behaviour checking based on linear time temporal logic. CONCLUSIONS: The advantages of our modelling and analysis approach are the description of quorum sensing and associated biofilm formation over two phases of bacterial growth, taking into account bacterial spatial distribution using a flexible and easy to maintain computational model. All computational results are reproducible.

Coloured Petri nets for multilevel, multiscale and multidimensional modelling of biological systems.

Owing to the availability of data of one biological phenomenon at different levels/scales, modelling of biological systems is moving from single level/scale to multiple levels/scales, which introduces a number of challenges. Coloured Petri nets (ColPNs) have been successfully applied to multilevel, multiscale and multidimensional modelling of some biological systems, addressing many of these challenges. In this article, we first review the basics of ColPNs and some popular extensions, and then their applications for multilevel, multiscale and multidimensional modelling of biological systems. This understanding of how to use ColPNs for modelling biological systems will assist readers in selecting appropriate ColPN classes for specific modelling circumstances.

Towards dynamic genome-scale models.

The analysis of the dynamic behaviour of genome-scale models of metabolism (GEMs) currently presents considerable challenges because of the difficulties of simulating such large and complex networks. Bacterial GEMs can comprise about 5000 reactions and metabolites, and encode a huge variety of growth conditions; such models cannot be used without sophisticated tool support. This article is intended to aid modellers, both specialist and non-specialist in computerized methods, to identify and apply a suitable combination of tools for the dynamic behaviour analysis of large-scale metabolic designs. We describe a methodology and related workflow based on publicly available tools to profile and analyse whole-genome-scale biochemical models. We use an efficient approximative stochastic simulation method to overcome problems associated with the dynamic simulation of GEMs. In addition, we apply simulative model checking using temporal logic property libraries, clustering and data analysis, over time series of reaction rates and metabolite concentrations. We extend this to consider the evolution of reaction-oriented properties of subnets over time, including dead subnets and functional subsystems. This enables the generation of abstract views of the behaviour of these models, which can be large-up to whole genome in size-and therefore impractical to analyse informally by eye. We demonstrate our methodology by applying it to a reduced model of the whole-genome metabolism of Escherichia coli K-12 under different growth conditions. The overall context of our work is in the area of model-based design methods for metabolic engineering and synthetic biology.

Emerging ensembles of kinetic parameters to characterize observed metabolic phenotypes.

BACKGROUND: Determining the value of kinetic constants for a metabolic system in the exact physiological conditions is an extremely hard task. However, this kind of information is of pivotal relevance to effectively simulate a biological phenomenon as complex as metabolism. RESULTS: To overcome this issue, we propose to investigate emerging properties of ensembles of sets of kinetic constants leading to the biological readout observed in different experimental conditions. To this aim, we exploit information retrievable from constraint-based analyses (i.e. metabolic flux distributions at steady state) with the goal to generate feasible values for kinetic constants exploiting the mass action law. The sets retrieved from the previous step will be used to parametrize a mechanistic model whose simulation will be performed to reconstruct the dynamics of the system (until reaching the metabolic steady state) for each experimental condition. Every parametrization that is in accordance with the expected metabolic phenotype is collected in an ensemble whose features are analyzed to determine the emergence of properties of a phenotype. In this work we apply the proposed approach to identify ensembles of kinetic parameters for five metabolic phenotypes of E. Coli, by analyzing five different experimental conditions associated with the ECC2comp model recently published by Hädicke and collaborators. CONCLUSIONS: Our results suggest that the parameter values of just few reactions are responsible for the emergence of a metabolic phenotype. Notably, in contrast with constraint-based approaches such as Flux Balance Analysis, the methodology used in this paper does not require to assume that metabolism is optimizing towards a specific goal.