ABSTRACT
Colorectal cancer morphogenesis begins at the cellular level from cell mutations in the intestinal epithelium cavities called crypts. These mutations lead to a pressure difference in the epithelium crypt walls, which can cause deformation and generate visible abnormalities in the epithelium. The geometrical modeling of these crypts and the mathematical modeling of the biomechanical process that leads to deformations can be simulated by using a Finite Element Method. The method solves numerically the system of partial differential equations (PDEs) that governs this phenomenon and permits to estimate the deformations of the crypt walls. In this work we simulate the crypt deformation when the cell mutations appear in several regions of the crypt epithelium.
Keywords:
colorectal cancer; crypts; finite element method; computer simulation
1 INTRODUCTION
In the process of the development of colorectal cancer, called carcinogenesis, cell mutations can arise in colon cavities, called crypts, containing transit and stem cells 3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. Mutant cells usually show higher proliferation rate than normal cells, which causes a larger pressure on the crypts walls than that generated by normal cells. This leads to a deformation of the crypt walls and its top orifice. If mutant cells fill adjacent crypts, a set of abnormal and deformed crypts appear. Thereafter there is a growth of the epithelium tissue caused by these abnormal crypts inside the colon that evolves to an adenoma. Adenomas in the colon are characterized in early stages by the formation of larger and deformed orifices of the epithelium crypts, which have a displacement in the intestinal lumen 1414 K. Drasdo & M. Loeffler . Individual-based models to growth and folding in one-layered tissues: intestinal crypts and early development. Nonlinear Analysis: Theory, Methods & Applications, 47(1) (2001), 245-256.), (3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181..
Regarding the cell and tissue growth modeling, we present in Section 2 a spatial continuous model based on PDEs, similar to that used in 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221., that determines the density and adhesion cell-cell pressure of fully differentiated and proliferative transit cells in a single crypt.
The exact morphology resulting from an abnormal proliferation of cells at some locations along the crypt axis is not known. However this leads to a deformation that can proceed to a budding, a orifice enlargement or an abnormal crypt fission 11 A.A. Almet, H.M. Byrne, P.K. Maini & D.E. Moulton. The role of mechanics in the growth and homeostasis of the intestinal crypt. Journal of Theoretical Biology, 20 (2021), 585-608.), (1515 C.M. Edwards & S.J. Chapman . Biomechanical Modelling of Colorectal Crypt Budding and Fission. Bulletin of Mathematical Biology , 69(6) (2007), 1927-1942.. We follow a model similar of that introduced in 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221., but instead of using a biomechanical model considering only the crypt orifice deformation, we consider the entire crypt walls deformation which permits to simulate the colon epithelium tissue growth generated by abnormal crypt cells located in any position in relation to the crypt base. This allows to be more effective in understanding how orifice enlargement and deformation observed in usual colonoscopy images can be originated by abnormal cell proliferation along the usual colonoscopy camera hidden crypt axis. Understanding such origin of crypt morphology permits to associate different orifice deformation shapes to particular locations of the abnormal cells along the crypt axis. It is then feasible to apply a specific target therapy to these locations to block the adenoma formation since cells at different crypt axis quotes have different level of Wnt/-catenin/TCF4 signaling pathway and APC protein. Wnt pathway is for instance more active in the lower part of the crypt and decreases moving upwards up to be null at the crypt top, whereas APC protein has an opposite gradient along the crypt axis 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.), (2323 P. Hogeweg. Evolving Mechanisms of Morphogenesis: on the Interplay between Differential Adhesion and Cell Differentiation. Journal of Theoretical Biology , 203(4) (2000), 317-333.. Such a model based on longitudinal crypt section can allow a model calibration based on real data of a three dimensional crypt as that obtained by the scanning electron microscope SEM as that used in 44 K. Araki, T. Ogata, M. Kobayashi & R. Yatani. A Morphological Study on the Histogenesis of Human Colorectal Hyperplastic Polyps. Gastroenterology, 109 (1995), 1468-1474.), (2525 Y. Kuratani, S. Tamura, Y. Furuya & S. Onishi. Morphogenesis of a colorectal neoplasm with a type III S pit pattern inferred from isolated crypts. J Gastroenterol, 43 (2008), 597-602..
In literature, there are many spatial models describing the cell spatial location, also called cellbased models, including in-lattice or grid-models 2323 P. Hogeweg. Evolving Mechanisms of Morphogenesis: on the Interplay between Differential Adhesion and Cell Differentiation. Journal of Theoretical Biology , 203(4) (2000), 317-333.),(3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.),(3131 S.Y. Wong, K.H. Chiam, C.T. Lim & P. Matsudaira. Computational model of cell positioning: directed and collective migration in the intestinal crypt epithelium. Journal of The Royal Society Interface, 7(3) (2010), S351-S363. and off-lattice or lattice-free models 22 A.A. Almet , B.D. Hughes, K.A. Landman, I.S. Näthke & J.M. Osborne. A Multicellular Model of Intestinal Crypt Buckling and Fission. Bulletin of Mathematical Biology, 80(2) (2018), 335-359.), (99 P. Buske, J. Galle, N. Barker, G. Aust, H. Clevers & M. Loeffler. A Comprehensive Model of the Spatio-Temporal Stem Cell and Tissue Organisation in the Intestinal Crypt. PLoS Computational Biology, 7(1) (2011), e1001045.), (1919 J. Galle , M. Hoffmann & G. Aust . From single cells to tissue architecture -a bottom-up approach to modelling the spatio-temporal organisation of complex multi-cellular systems. Journal of Mathematical Biology , 58(1) (2009), 261-283.. The choice to preferring a spatial continuum model with respect to a cell-based model is in the possibility of extend its application to simulate cell dynamics in millions of crypts that can be done using a multiscale strategy, as done in 1616 I.N. Figueiredo, C. Leal, G. Romanazzi & B. Engquist. Homogenization model for aberrant crypt foci. SIAM Journal on Applied Mathematics, 76(3) (2016), 1152-1177..
For the tissue growth modeling, as highlighted in reviews 33 A.A. Almet , P.K. Maini , D.E. Moulton & H.M. Byrne . Modeling perspectives on the intestinal crypt, a canonical system for growth, mechanics, and remodeling. Current Opinion in Biomedical Engineering, 15 (2020), 32-39.),(1111 G. De Matteis, A. Graudenzi & M. Antoniotti. A review of spatial computational models for multicellular systems, with regard to intestinal crypts and colorectal cancer development. Journal of Mathematical Biology, 66(7) (2013), 1409-1462.),(2424 S.K. Kershaw, H.M. Byrne , D.J. Gavaghan & J.M. Osborne . Colorectal cancer through simulation and experiment. IET Systems Biology, 7(3) (2013), 57-73., many possible models based on cell dynamics in the intestinal crypts have been considered, among them: continuum mechanics models 1515 C.M. Edwards & S.J. Chapman . Biomechanical Modelling of Colorectal Crypt Budding and Fission. Bulletin of Mathematical Biology , 69(6) (2007), 1927-1942.), (2727 M.R. Nelson, J.R. King & O.E. Jensen. Buckling of a growing tissue and the emergence of twodimensional patterns. Mathematical Biosciences, 246(2) (2013), 229-241.; compartmental models describing the transition between cell types, regardless of their position 77 M. Bjerknes. Expansion of Mutant Stem Cell Populations in the Human Colon. Journal of Theoretical Biology , 178(4) (1996), 381-385.), (88 B.M. Boman, J.Z. Fields, O. Bonham-Carter & O.A. Runquist. Computer Modeling Implicates Stem Cell Overproduction in Colon Cancer Initiation. Cancer Research, 61(23) (2001), 8408-8411.; and nonspatial stochastic models based on the use of epidemiological data to predict cancer formation or some genetic instability 1212 A. Di Garbo, M.D. Johnston, S.J. Chapman & P.K. Maini . Variable renewal rate and growth properties of cell populations in colon crypts. Physical Review E, 81(6) (2010), 061909.), (1313 A. d’Onofrio & I.P.M. Tomlinson. A nonlinear mathematical model of cell turnover, differentiation and tumorigenesis in the intestinal crypt. Journal of Theoretical Biology , 244(3) (2007), 367-374.. Our model determines the displacement and deformation of the epithelium tissue along the crypt based on a continuum equilibrium relation, this allows us to use directly the pressure gradient between cells as a force acting on the tissue.
2 MODELING CRYPT CELL DYNAMICS AND DEFORMATION
In this Section we present the geometrical domain used to represent a crypt and two half crypts at its left and right in a cross section of the colon epithelium. Then we introduce the PDE system that describes the cell dynamics in this domain, see Subsection 2.2. Later in the Subsection 2.3 we define a viscoelastic model responsible for the tissue deformation along the crypt, which coupled with the cell dynamics model allows to deform the geometrical domain of the crypt and to observe, in some abnormal cases, tissue growth as it is shown in Section 4.
2.1 Geometrical Domain
The colon epithelium is formed by a periodic distribution of crypts that in three dimension can be locally represented as in Figure 1. We study and analyze the deformation occurring in a small region Ω of the colon epithelium during a time interval [0, T]. Its external crypt subdomain ΩE is formed by the cross section of the external epithelium tissue, which has one crypt in its center and two walls of two adjacent crypts at its left and right, see Figure 2.
Domain ΩE formed by a cross section of a single crypt in the center with two adjacent crypt walls at left and right.
The domain Ω has also a tissue region ΩC that is formed by the connective tissue of the colon epithelium between the crypts. Thus we model the deformations and cell dynamics in , see Figure 3. In ΩE , we use a diffusive-convective model for the dynamics of proliferative (also called transit) cells and fully differentiated cells. Both ΩE and ΩC have elastic structures that can move during the time interval [0, T]. Figure 4 illustrates the boundary domain subdivisions and their notations. Region ΩE has an external boundary formed by Γr that contains its upper circular borders, Γ3 is the vertical boundary of the ΩE bottom and Γtop is the horizontal boundary of the ΩE top. We note that normal mature (fully differentiated) cells are released in the colon lumen through this last boundary Γtop3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. The remaining external ΩE boundary is denoted by Γ1. In the interface between ΩE and ΩC we localize Γbot that represents the circular bottom boundary, and Γ2 that is the remaining interface boundary. Furthermore we denoted by Γ4 the lower ΩC boundary where the tissue is supposed to be fixed.
Full Problem domain Ω with the external crypt tissue ΩE and the internal connectivity tissue ΩC .
Since we deform such domains in time, we denote by the geometrical domain of the crypt at time that is deformed in relation to the original domain Ω represented in the figures above. Similarly we denote by and the deformed external and connective tissue domain at time t.
2.2 PDE model for cell dynamics
We suppose, as done in 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.), (1818 I.N. Figueiredo , C. Leal , G. Romanazzi , B. Engquist & P.N. Figueiredo. A convection-diffusionshape model for aberrant colonic crypt morphogenesis. Computing and Visualization in Science, 14(4) (2011), 157-166., that a crypt contains two families of cells: proliferative transit cells with density and fully differentiated living cells with density where and . It is also supposed that N 1 and N 2 satisfy the overall density condition and that cell families have the same diffusion coefficient D as it is asserted in 2828 T. Roose, S.J. Chapman & P.K. Maini . Mathematical Models of Avascular Tumor Growth. SIAM REVIEW, 49(2) (2007), 179-208.. According to the model used in 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.),(1818 I.N. Figueiredo , C. Leal , G. Romanazzi , B. Engquist & P.N. Figueiredo. A convection-diffusionshape model for aberrant colonic crypt morphogenesis. Computing and Visualization in Science, 14(4) (2011), 157-166., N 1 and N 2 satisfy in for the following PDE system
where α is the proliferation rate of N 1, β is the proliferation rate of N 2, p is the cell-cell adhesion pressure responsible for the upwards cell transition along the crypt axis. Since , if we sum the first two equations of (2.1) and take , the temporal and spatial derivatives cancel. Thus we obtain that and satisfy in for :
The boundary conditions for the pressure are
In the second equation of (2.3), η is the modulus of y-coordinate of outward normal vector in Γr . This is a Robin boundary condition that allows the transition between the Dirichlet and Neumann conditions in Γ1 and Γtop .
For the density, the boundary conditions are
The proliferation rate α is higher at the crypt bottom and becomes null at the top, normally in the upper third of the crypt height where most of the fully differentiated cells, that are not proliferative, reside 3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. On the other hand, the transformation rate β behaves in the opposite way. It means that it is higher at the crypt top and null at the bottom, thus we use as in 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.), (1818 I.N. Figueiredo , C. Leal , G. Romanazzi , B. Engquist & P.N. Figueiredo. A convection-diffusionshape model for aberrant colonic crypt morphogenesis. Computing and Visualization in Science, 14(4) (2011), 157-166. the rates
where h denotes the crypt height and is a constant that is determined by using some biological information of the analyzed crypt such as the cell size and the cell convective velocity along the crypt axis.
The full model (2.2) with the boundary conditions (2.3)-(2.4) guarantees that transit cells diffuse in the crypt, proliferate with rate α and differentiate with rate β . Moreover these cells have a convective upward velocity . This velocity increases along the vertical crypt axis up to reach the crypt top, where fully differentiated cells are released in the lumen, as it has been observed experimentally 3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. We note in fact that the boundary conditions for the pressure and density in Γ permit cells to move upwards, and to be released in the colon lumen across the Γtop boundary.
2.3 The Viscoelastic Model
As discussed in 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221., the crypt domain Ω can be considered as a viscoelastic structure. Then it can have a displacement in relation to its original position, where u 1 and u 2 are the displacement along the x−axis and y−axis of point x at time t, respectively. Since connective tissue domain base is normally considered as fixed we allow only Γ1, Γ2, Γbot , Γtop , Γr to be free boundaries whereas Γ3, Γ4 remain fixed. The gradient of pressure differences acts as a force on the external crypt boundary , which causes a deformation of Ω, with a displacement u. The stress tensor has elastic and viscous components, so that . Since ΩE contains only cells it relaxes more linearly than ΩC and so it can be represented as follows with an elastic structure, see 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.:
In (2.6) λ E and µ E are the Lamé coefficients, which describe mechanical properties of elastic materials 2929 J. Salençon. “Handbook of Continuum Mechanics: General Concepts, Thermoelasticy”. Springer Science & Business Media (2012).. On the other hand, in ΩC , we have a viscoelastic structure that has elastic and viscous components. We use a Standard Linear Solid (SLS) model for modeling the viscoelastic tissue domain ΩC . This model is largely used in the representation of biological tissues 2626 N.K. Kyslstad. “Simulating the viscoelastic response of the spinal cord”. Ph.D. thesis, Faculty of Mathematics and Natural Sciences, University of Oslo, Oslo (2014).. Then, using the Einstein summation convention, we have
where
with
where with , see 1010 E. Carniel, M. Mencattelli & G.e.a. Bonsignori. Analysis of the structural behaviour of colonic segments by inflation tests: Experimental activity and physio-mechanical model. Proc Inst Mech Eng H., 229(11) (2015), 794-803.), (1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221..
At the equilibrium, the crypt structure displacements satisfy for for all
where and is the outward normal vector in Γ1. Furthermore, we have , with being the pressure exerted between the cells in the normal case and p the exerted pressure in the abnormal case, that causes the crypt deformation, and γ is an adimensional term related with the current geometry of the tissue domain and cell dimensions.
3 METHODOLOGY
We simulate the crypt deformations and the associated tissue growth during the first stages of colorectal cancer in which abnormalities occur in a single crypt as shown in Section 4. We use a classic Finite Element Method with piecewise linear basis to solve the differential problem in different time instants and we advance in time with a fixed time interval ∆t by using a Backward Euler’s method. Once the crypt deforms as the process goes forward, the domain mesh changes in each time iteration. The computational code is implemented using the FreeFem++ software 2222 F. Hecht . New development in FreeFem++. J. Numer. Math., 20(3-4) (2012), 251-265. URL https://freefem.org/.
https://freefem.org/...
, that solves numerically the density N, pressure p, and displacement u by using a piecewise linear finite element basis of the domain at each .
4 RESULTS AND DISCUSSIONS
In order to simulate different crypt deformations and the associated tissue growth due to abnormal crypt cell proliferation, we need to assert the normal case and how it differs with respect the abnormal cases.
4.1 Stable Solution for N
In order to have a stable solution of system (2.2) with boundary conditions we find an expression for the transformation rate β, for a given α and a stable solution N. We suppose that
in which τ is the relaxation constant, and the stable solution is
Where , in which yt is the minimum (vertical) coordinate y of the top region and y b is the maximum (vertical) coordinate of the bottom region. Since the solution N is stable, we have by hypothesis that . Then, from the first equation of (2.2), we have
But, from the second equation, , with D constant, we obtain in
Moreover since and , we obtain
Since on and on Γtop , we assume that in order to get a simpler solution. Then, recalling that , we have
Integrating both sides of the equation we obtain
for all .
We know that, at the crypt top, position/hour = 0.85 h cell µm/h3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. Moreover, since a crypt has a height 433µm on average, and cells have height 5.9µm2020 D.V. Guebel & N.V. Torres. A computer model of oxygen dynamics in human colon mucosa: implications in normal physiology and early tumor development. Journal of Theoretical Biology , 250(3) (2008), 389-409., then for a crypt in the model with height h crypt , we have
Thus . Then
If , we have
On the other hand, if , we have
Thus, we get
if and
if . At last, if we define β only for with and recall that , we obtain
where is a solution of the problem
Figures 5 and 6 show the stable solution for density and the pressure solution of the equation .
4.2 Numerical Simulations
We tested many abnormal examples each leading to a different cell dynamics and crypt deformation. In each test the abnormality, characterized by a region with an unexpected high proliferation rate α, is located in a different region. We use the following parameter values in each simulation: crypt height (µm)2/hour, hour, hour−1(µm)−2 and µm, with H max being the maximum diameter of the mesh triangles. Moreover we use the boundary conditions for the pressure and density (2.3)-(2.4) and the initial condition given in (4.2) with and .
In all simulations, meshes are automatically generated by the software FreeFem++. The finite element mesh used at the beginning of the simulations, at time t=0, have 17174 triangles and 9307 vertices. After each time step, the mesh is updated based on the calculated displacements using the “movemesh” routine. If this update causes overlapping of the elements, FreeFem++ adapts the mesh using the “adaptmesh” routine based on the BAMG mesh adaptive method 2121 F. Hecht. BAMG: bidimensional anisotropic mesh generator. User Guide. Technical report, INRIA, Rocquencourt (1998)..
4.2.1 Test I
In the first test the abnormality is centered in the lower region of the crypt, as shown by the black circle in Figure 7.
Initial mesh and the abnormality location marked by a black circle located in the crypt bottom.
After 15 time steps, the solutions u 1 and u 2 calculated on the updated mesh are shown in 8 and 9. Moreover, solutions N and are plotted respectively in Figures 10 and 11.
We observe that the displacements along the x-axis are symmetric since they are the same in the right and left crypt walls, according to the expectations. It results then in a symmetric deformation problem. Moreover, we note a relevant opening of the crypt orifice along the time, that characterizes microadenomas such as the Aberrant Crypt Foci1 1 Aberrant crypt foci (ACF) 5), (6 are clusters of crypts in the colon epithelium, containing cells with a deviant behavior with respect to the normal ones such as those hyper proliferating that due to their action can deform the crypt structure giving rise to aberrant crypt orifices. It is believed that ACF are precursors of colorectal cancer morphogenesis. , see 1616 I.N. Figueiredo, C. Leal, G. Romanazzi & B. Engquist. Homogenization model for aberrant crypt foci. SIAM Journal on Applied Mathematics, 76(3) (2016), 1152-1177.), (1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.), (3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. Besides that, we have a large pressure difference in the region of the abnormality. This pressure difference decreases as we ascend the crypt. Regarding cell density N, we note that it is equal to 1 at the crypt base and decreases its value moving upwards in the direction of the crypt top where it has a null value.
4.2.2 Test II
In the second simulation the abnormality is located in the left branch (wall) of the crypt, in a centered position in relation to the y-axis, as it can be seen by the black circle in Figure 12.
After 20 time steps, the displacements u 1 and u 2 computed on the updated mesh are shown in Figures 13 and 14. N and are respectively shown in Figures 15 and 16.
In this simulation we observe that there is a larger displacement in the left crypt wall in relation to the portion on the right. This behavior agrees with what is expected because the abnormality is on the left side and, therefore, exerts larger influence in this region than on the right. Once again the pressure difference is large close to the abnormality and this difference decreases as we move away from this region. In this case, the density behavior is similar to that observed in test 1: it varies from one to zero, from the bottom to the top.
4.2.3 Test III
In the third situation, the abnormality is on the left branch of the crypt, in a higher position, as shown by Figure 17.
After 20 time steps, the solutions to u 1 and u 2 calculated on the updated mesh are shown in Figures 18 and 19. Besides that, solutions to N and are respectively shown in Figures 20 and 21.
In this simulation, the region on the left of the crypt moves, forming an opening in this direction. The right-hand portion of the crypt has a very small displacement. This occurs because the difference in pressure in the right portion is very close to zero, since the abnormality is found in the left portion. Again, cell density decreases along the crypt vertical axis, from one at the bottom to zero at the top.
4.2.4 Test IV
In the fourth simulation, the abnormality is centered in the upper right portion of the crypt, as shown by the black circle in Figure 22.
After 10 time steps, the solutions to u 1 and u 2 calculated on the updated mesh are shown in Figures 23 e 24. Moreover, the solutions N and are respectively shown in Figures 25 and 26.
In this test we observe an elongation of the left domain region of the epithelium due to the presence of upward-oriented forces in this region. Such forces come from the pressure difference in the abnormality region. The displacements in y-axis, shown by u 2, are larger than those in x-axis. Moreover, the pressure difference is larger in the upper left portion of the crypt, being practically null in the entire right branch. Cell density presents a behavior similar to that of the other tests, decreasing from the crypt bottom to the crypt top.
5 CONCLUSIONS
The process of colorectal cancer development arises from the presence of abnormal cells in the intestinal crypts. This leads to a pressure difference, which generates tissue displacements. Using a PDE model, it was possible to numerically simulate the tissue deformation over time by using the Finite Element Method. This is a problem in which the domain is updated in each time iteration.
We analyzed how different positions of the abnormality in the tissue affects the process evolution. The displacements occurs in different ways depending on the inital location of abnormal proliferative cells. This is in line with what is biologically expected since abnormal cell proliferation can lead to forces acting in the tissue epithelium 11 A.A. Almet, H.M. Byrne, P.K. Maini & D.E. Moulton. The role of mechanics in the growth and homeostasis of the intestinal crypt. Journal of Theoretical Biology, 20 (2021), 585-608.), (1515 C.M. Edwards & S.J. Chapman . Biomechanical Modelling of Colorectal Crypt Budding and Fission. Bulletin of Mathematical Biology , 69(6) (2007), 1927-1942..
The crypt orifice deformation is an early sign of the formation of adenomas in the colon 1717 I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.), (3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. Therefore, with this research, we were able to model how abnormal cells characterized by excessive proliferation can lead to crypt deformation that can occur before cancer appears 3030 I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.. This work suggests that we can understand where the abnormal cells are acting, observing the deformed orifice and crypt walls. This can be useful to force target drugs to act in specific colon epithelium regions in order to block the colorectal cancer formation or reducing its growth.
REFERENCES
-
1A.A. Almet, H.M. Byrne, P.K. Maini & D.E. Moulton. The role of mechanics in the growth and homeostasis of the intestinal crypt. Journal of Theoretical Biology, 20 (2021), 585-608.
-
2A.A. Almet , B.D. Hughes, K.A. Landman, I.S. Näthke & J.M. Osborne. A Multicellular Model of Intestinal Crypt Buckling and Fission. Bulletin of Mathematical Biology, 80(2) (2018), 335-359.
-
3A.A. Almet , P.K. Maini , D.E. Moulton & H.M. Byrne . Modeling perspectives on the intestinal crypt, a canonical system for growth, mechanics, and remodeling. Current Opinion in Biomedical Engineering, 15 (2020), 32-39.
-
4K. Araki, T. Ogata, M. Kobayashi & R. Yatani. A Morphological Study on the Histogenesis of Human Colorectal Hyperplastic Polyps. Gastroenterology, 109 (1995), 1468-1474.
-
5R.P. Bird. Role of aberrant crypt foci in understanding the pathogenesis of colon cancer. Cancer Letters, 93 (1995), 55-71.
-
6R.P. Bird & C.K. Good. The significance of aberrant crypt foci in understanding the pathogenesis of colon cancer. Toxicology Letters, 112-113 (2000), 395-402.
-
7M. Bjerknes. Expansion of Mutant Stem Cell Populations in the Human Colon. Journal of Theoretical Biology , 178(4) (1996), 381-385.
-
8B.M. Boman, J.Z. Fields, O. Bonham-Carter & O.A. Runquist. Computer Modeling Implicates Stem Cell Overproduction in Colon Cancer Initiation. Cancer Research, 61(23) (2001), 8408-8411.
-
9P. Buske, J. Galle, N. Barker, G. Aust, H. Clevers & M. Loeffler. A Comprehensive Model of the Spatio-Temporal Stem Cell and Tissue Organisation in the Intestinal Crypt. PLoS Computational Biology, 7(1) (2011), e1001045.
-
10E. Carniel, M. Mencattelli & G.e.a. Bonsignori. Analysis of the structural behaviour of colonic segments by inflation tests: Experimental activity and physio-mechanical model. Proc Inst Mech Eng H., 229(11) (2015), 794-803.
-
11G. De Matteis, A. Graudenzi & M. Antoniotti. A review of spatial computational models for multicellular systems, with regard to intestinal crypts and colorectal cancer development. Journal of Mathematical Biology, 66(7) (2013), 1409-1462.
-
12A. Di Garbo, M.D. Johnston, S.J. Chapman & P.K. Maini . Variable renewal rate and growth properties of cell populations in colon crypts. Physical Review E, 81(6) (2010), 061909.
-
13A. d’Onofrio & I.P.M. Tomlinson. A nonlinear mathematical model of cell turnover, differentiation and tumorigenesis in the intestinal crypt. Journal of Theoretical Biology , 244(3) (2007), 367-374.
-
14K. Drasdo & M. Loeffler . Individual-based models to growth and folding in one-layered tissues: intestinal crypts and early development. Nonlinear Analysis: Theory, Methods & Applications, 47(1) (2001), 245-256.
-
15C.M. Edwards & S.J. Chapman . Biomechanical Modelling of Colorectal Crypt Budding and Fission. Bulletin of Mathematical Biology , 69(6) (2007), 1927-1942.
-
16I.N. Figueiredo, C. Leal, G. Romanazzi & B. Engquist. Homogenization model for aberrant crypt foci. SIAM Journal on Applied Mathematics, 76(3) (2016), 1152-1177.
-
17I.N. Figueiredo , C. Leal , G. Romanazzi & B. Engquist . Biomathematical model for simulating abnormal orifice patterns in colonic crypts. Mathematical biosciences, 315 (2019), 108221.
-
18I.N. Figueiredo , C. Leal , G. Romanazzi , B. Engquist & P.N. Figueiredo. A convection-diffusionshape model for aberrant colonic crypt morphogenesis. Computing and Visualization in Science, 14(4) (2011), 157-166.
-
19J. Galle , M. Hoffmann & G. Aust . From single cells to tissue architecture -a bottom-up approach to modelling the spatio-temporal organisation of complex multi-cellular systems. Journal of Mathematical Biology , 58(1) (2009), 261-283.
-
20D.V. Guebel & N.V. Torres. A computer model of oxygen dynamics in human colon mucosa: implications in normal physiology and early tumor development. Journal of Theoretical Biology , 250(3) (2008), 389-409.
-
21F. Hecht. BAMG: bidimensional anisotropic mesh generator. User Guide. Technical report, INRIA, Rocquencourt (1998).
-
22F. Hecht . New development in FreeFem++. J. Numer. Math., 20(3-4) (2012), 251-265. URL https://freefem.org/
» https://freefem.org/ -
23P. Hogeweg. Evolving Mechanisms of Morphogenesis: on the Interplay between Differential Adhesion and Cell Differentiation. Journal of Theoretical Biology , 203(4) (2000), 317-333.
-
24S.K. Kershaw, H.M. Byrne , D.J. Gavaghan & J.M. Osborne . Colorectal cancer through simulation and experiment. IET Systems Biology, 7(3) (2013), 57-73.
-
25Y. Kuratani, S. Tamura, Y. Furuya & S. Onishi. Morphogenesis of a colorectal neoplasm with a type III S pit pattern inferred from isolated crypts. J Gastroenterol, 43 (2008), 597-602.
-
26N.K. Kyslstad. “Simulating the viscoelastic response of the spinal cord”. Ph.D. thesis, Faculty of Mathematics and Natural Sciences, University of Oslo, Oslo (2014).
-
27M.R. Nelson, J.R. King & O.E. Jensen. Buckling of a growing tissue and the emergence of twodimensional patterns. Mathematical Biosciences, 246(2) (2013), 229-241.
-
28T. Roose, S.J. Chapman & P.K. Maini . Mathematical Models of Avascular Tumor Growth. SIAM REVIEW, 49(2) (2007), 179-208.
-
29J. Salençon. “Handbook of Continuum Mechanics: General Concepts, Thermoelasticy”. Springer Science & Business Media (2012).
-
30I.M.M. Van Leeuwen, H.M. Byrne , O.E. Jensen & J.R. King . Crypt dynamics and colorectal cancer: advances in mathematical modelling. Cell proliferation, 39(3) (2006), 157-181.
-
31S.Y. Wong, K.H. Chiam, C.T. Lim & P. Matsudaira. Computational model of cell positioning: directed and collective migration in the intestinal crypt epithelium. Journal of The Royal Society Interface, 7(3) (2010), S351-S363.
-
1
Aberrant crypt foci (ACF) 55 R.P. Bird. Role of aberrant crypt foci in understanding the pathogenesis of colon cancer. Cancer Letters, 93 (1995), 55-71.), (66 R.P. Bird & C.K. Good. The significance of aberrant crypt foci in understanding the pathogenesis of colon cancer. Toxicology Letters, 112-113 (2000), 395-402. are clusters of crypts in the colon epithelium, containing cells with a deviant behavior with respect to the normal ones such as those hyper proliferating that due to their action can deform the crypt structure giving rise to aberrant crypt orifices. It is believed that ACF are precursors of colorectal cancer morphogenesis.
Publication Dates
-
Publication in this collection
04 Apr 2022 -
Date of issue
Jan-Mar 2022
History
-
Received
21 Mar 2021 -
Accepted
24 Aug 2021