A mathematical model of aortic aneurysm formation

Abstract

Abdominal aortic aneurysm (AAA) is a localized enlargement of the abdominal aorta, such that the diameter exceeds 3 cm. The natural history of AAA is progressive growth leading to rupture, an event that carries up to 90% risk of mortality. Hence there is a need to predict the growth of the diameter of the aorta based on the diameter of a patient’s aneurysm at initial screening and aided by non-invasive biomarkers. IL-6 is overexpressed in AAA and was suggested as a prognostic marker for the risk in AAA. The present paper develops a mathematical model which relates the growth of the abdominal aorta to the serum concentration of IL-6. Given the initial diameter of the aorta and the serum concentration of IL-6, the model predicts the growth of the diameter at subsequent times. Such a prediction can provide guidance to how closely the patient’s abdominal aorta should be monitored. The mathematical model is represented by a system of partial differential equations taking place in the aortic wall, where the media is assumed to have the constituency of an hyperelastic material.

About the Authors:

Wenrui Hao

* E-mail: [email protected]

Affiliation: Department of Mathematics, Pennsylvania State University, State College, PA, United States of America

Shihua Gong

Affiliation: Beijing International Center for Mathematical Research, Peking University, Beijing, China

Shuonan Wu

Affiliation: Department of Mathematics, Pennsylvania State University, State College, PA, United States of America

Jinchao Xu

Affiliation: Department of Mathematics, Pennsylvania State University, State College, PA, United States of America

Michael R. Go

Affiliation: Division of Vascular Diseases & Surgery, The Ohio State University Wexner Medical Center, Columbus, OH, United States of America

Avner Friedman

Affiliation: Department of Mathematics, The Ohio State University, Columbus, OH, United States of America

Dai Zhu

Affiliation: Yunnan Cocreative Scientific Computing and Datamining Center, Kunming, China

Introduction

AAA is an abnormal dilatation most commonly of the infrarenal aorta. One definition of AAA is a diameter greater than 3 cm [1]. The clinical significance of AAA stems from the high mortality associated with rupture. Approximately 60% of patients with ruptured AAA die before reaching the hospital [2] and mortality rates of emergency surgical repair are as high as 35–70% [3]. The risk of rupture increases with AAA diameter. The pathogenesis of AAA is largely unknown and likely multifactorial. Diameter remains the only clinically useful and available marker for risk of rupture. Surgical repair is recommended for aneurysms measuring greater than 5.5 cm [4, 5] and there is insufficient evidence to recommend surgery for all patients with smaller AAA [2]. However, even elective surgical repair of AAA remains a major operation, incurring significant morbidity and mortality particularly in older, sicker patients where there is a higher prevalence of AAA. Furthermore, patients die every year from rupture of aneurysms smaller than 5.5 cm, and up to 60% of AAA larger than 5 cm remain stable [6]. Thus, some patients with smaller AAA are denied lifesaving surgery and others with larger AAA undergo unnecessary major surgery. The reasons why some smaller aneurysms go on to rupture while some larger ones remain stable are not understood. Techniques that provide early identification of small AAA with increased risk for rupture and large AAA with low risk for rupture will improve overall mortality by prompting personalized treatment plans for AAA.

There are a number of mathematical papers that describe the dynamics of the weakening and dilation of the arterial wall and the risk of rupture. The more comprehensive models include the nonlocal, nonlinear elastic nature of the arterial wall in response to wall shear stress [7–10]. However these models do not explain the clinical fact that some abdominal aortae rupture when the diameter of the cross section of the dilated aorta is less than 5.5 cm while others do not rupture even when much larger than 5.5 cm. There is a need to discover non-invasive biomarkers that will provide the following prognosis for a patient undergoing initial screening for AAA:

* Given initial diameter of the cross section of the aorta, say R₀, how will the diameter, say R(t), evolve with time?

* How will the risk of rupture increase with time?

The present paper develops a mathematical model that addresses the first question, with IL-6 as the serum biomarker. IL-6 has already been suggested as prognostic biomarker for AAA in [11, 12]; additional potential biomarkers are reviewed in [13, 14].

The mathematical model consists of two parts. The first part focuses on the biology. It introduces a network of cells and cytokines which, in disease state, lead to the destruction of smooth muscle cells (SMCs) and elastin depletion in the media, and disruption of the extracellular matrix (ECM) in the adventitia. The biological predictions of this part of the model, in terms of the expression levels of the proteins involved, are in agreement with patients data [15]. The second part of the model develops the mechanics of diameter growth of the artery as a result of deficiency in elastin and disruptions in the ECM. Here, instead of the approach taken in [7–10], we use the fact that the arterial wall behaves like a hyperelastic material [16]. Hence its elastic strength and displacement under stress depend on the smooth muscle cells in the aortic wall. We can therefore relate the wall dilation to several cytokines whose expression is upregulated in disease state. Such cytokines could be used as potential biomarkers for predicting the arterial wall dilation. We shall focus on IL-6, which is upregulated in AAA patients [11, 12].

The biological network introduced in the present paper includes macrophages, T cells, SMCs and fibroblasts, as well as cytokines such as IL-6, IL-10, IL-12, TNF-α and IFN-γ; MMP and TIMP are also included. The densities of the cells and the concentrations of the proteins satisfy a coupled system of partial differential equations (PDEs) in the arterial wall.

In the early stage of aneurysm, endothelial cells, under hemodynamic shear stress, secrete major attractant protein, MCP-1, and IL-6 [17]. IL-6 is also produced by SMCs (in conjunction with IL-1) [18]. This triggers recruitment of monocytes from the blood into the media layer of the arterial wall [19]. The monocytes mature into macrophages, while, at the same time, additional macrophages from the adventitia are also chemoattracted by MCP-1, IL-6 and IL-8 [20–22]. More precisely, IL6/MCP-1 combined, and MCP-1 and IL-8 separately chemoattract monocytes; for simplicity, we shall represent the chemoattractants by MCP-1 and IL-6. T cells are activated by contact with macrophages in the presence of IL-12, and macrophages are activated by IFN-γ produced by the T cells [23]. Macrophages produce TNF-α [24], MMP and TIMP [24, 25], and IL-10, IL-12 and IL-6 [26], but MMP is also produced by SMCs [27]. Fibroblasts produce collagen [28], and the collection of MMP, TIMP and collagen weaken the ability of the adventitia layer to withstand stress. Macrophages are known to cause apoptosis in SMCs [19], and this leads to reduction in elastin [29], thus weakening the elastic strength of the media; The apoptotic SMCs are known to produce MCP-1 [19].

Fig 1 shows SMCs residing in the media, fibroblasts residing in the adventitia, while macrophages and T cells are present in both layers of the aortic wall.

[Figure omitted. See PDF.]

Fig 1. A schematic of the aortic wall including the aortic lumen, the intima, media, and adventitia.

https://doi.org/10.1371/journal.pone.0170807.g001

Fig 2 is a schematics of the network within the arterial wall during aneurysm; all the cytokines and the ECM (elastin and collagen) are present in both the media and the adventitia. A schematics of the balloon-like bulge geometry in the aorta is shown in Fig 3.

[Figure omitted. See PDF.]

Fig 2. Interaction network among cells and cytokines in their respective layers.

SMCs remain in the media, and fibroblasts remain in the adventitia.

https://doi.org/10.1371/journal.pone.0170807.g002

[Figure omitted. See PDF.]

Fig 3. Schematic representation of a 2D section of the computational domain.

Ω_B is blood vessel; Ω_M represents the media and Ω_A the adventitia; Γ_A, Γ_B and Γ_M are three free boundaries. Γ_B is the inner surface of the artery, Γ_A is the outer surface, and Γ_M is the surface between the media and adventitia.

https://doi.org/10.1371/journal.pone.0170807.g003

1 Mathematical model of the biology

In this section we develop a mathematical model of aneurysm based on the diagram shown in Fig 2. The model, represented by a system of PDEs, includes the variables listed in Table 1. We assume that all cells are moving with a common velocity v; the velocity is the result of movement of macrophages, T cells and SMCs into the media and adventitia. We also assume that all species are diffusing with appropriate diffusion coefficients. The equation for each species of cells, X, has a formwhere the expression on the left-hand side includes advection and diffusion, and F_X accounts for various biochemical reactions, and chemotaxis. The equations for the chemical species are the same as for cells, but without the advection term, which is relatively very small compared to their large diffusion coefficients. Fig 3 shows a 2D projection of a blood vessel Ω_B with an aortic aneurysm Ω_M ∪ Ω_A.

[Figure omitted. See PDF.]

Table 1. The variables of the model; concentration and densities are in units of g/cm³.

https://doi.org/10.1371/journal.pone.0170807.t001

Equation for MCP-1 (P)

The equation for MCP-1 is given in the whole domain (media Ω_M, and adventitia Ω_A) by(1)The first term on the right-hand side is the production of P by apoptotic SMCs, which is induced by macrophages [19]. The second term represents lost of MCP-1 as it binds to, and internalized by, macrophages which are chemoattracted to MCP-1; the internalization of MCP-1 may be limited due to the limited rate of receptor recycling. The last term in Eq (1) accounts for degradation of MCP-1.

Equation for macrophages (M)

According to [26] the two phenotype of macrophages, M1 and M2, are both present in aneurysm; M1 secretes inflammatory cytokines such as IL-12 and TNF-α, while M2 secretes immunoregulatory cytokines such as IL-10. For simplicity, we do not differentiate between the two phenotypes, but take into account that the productions of IL-12 and TNF-α are resisted by IL-10 [30].

The evolution of the macrophage density in Ω_M ∪ Ω_A is modeled by(2)The first two terms on right-hand side account for the recruitment of macrophages by MCP-1 [19] and IL-6 [20–22]. The third term accounts for the activation of macrophages by IFN-γ [23], which is enhanced by TNF-α [31].

Equation for T cells (T)

The density of T cells in Ω_M ∪ Ω_A satisfies the equation(3)We assume that T cells are activated by contact with macrophages in IL-12 environment, while IL-12 is inhibited by IL-10 [30].

Equation for SMCs (S)

The equation of the SMCs density in Ω_M is given by(4)The second term of the right-hand side accounts for apoptosis of SMCs, caused by macrophages through a FasL/Fas-Caspase8-RIP1 mediated mechanism [19].

Equations for cytokines (in Ω_M ∪ Ω_A)

IL-6 is produced by macrophages [26], and by SMCs (in conjunction with IL-1) [18]; hence(5)

IL-10 is produced by macrophages [26], so that(6)

IL-12 is produced by macrophages [26], a process inhibited by IL-10. The equation for I₁₂ is given by(7)

Similarly, TNF-α is produced by macrophages, a process inhibited by IL-10 [32], so that(8)

IFN-γ is produced by T cells [33], hence(9)

MMP is produced by macrophages [25] and by SMCs [34], a process enhanced by TNF-α [24], and MMP is lost by binding with TIMP [25]. Accordingly, Q satisfies the following equation:(10)TIMP is produced by macrophages [25], and is lost by binding to MMP, so that(11)

Equation for ECM (ρ)

SMCs produce elastin [29] and fibroblasts produce collagen [28]. We assume that the density of fibroblasts is approximately constant, and that the ECM concentration is proportional to the combined concentrations of collagen and elastin, which are produced by fibroblasts and SMCs, respectively. Hence,(12)where the term d_ρQ Qρ represents ECM degradation by MMP. Here we used the notations: X⁺ = X if X > 0, X⁺ = 0 if X ≤ 0, and χ_A = 1 in Ω_A, χ_A = 0 in Ω_M.

Boundary conditions

Fig 3 shows a 2-D projection of the media and adventitia layers of the arterial wall. We shall take as our computational domain the right region lying between y = y₁ and y = y₂, and impose periodic boundary condition at y = y₁ and y = y₂.

T cells, which are primarily CD4+ Th1 cells [35], migrate to the adventitia and intima [36, 37], and so do macrophages [22, 37].

The concentration ρ is continuous across the media/adventitia membrane, while S satisfies no-flux boundary condition, , at this membrane. All the other variables satisfy the following flux conditions at the interface Γ_M:where X_M and X_A are the concentrations of cytokines or cell densities in media and adventitia, respectively. We take γ = 50 for cells, and γ = 500 for the cytokines [38].

On Γ_B and Γ_A, we havewhere .

Both MCP-1 and IL-6 are produced by endothelial cells [17], which lie near the inner boundary of the media, Γ_B. Hence(13)where P₀ and I₆₀ are the concentrations of MCP-1 and I₆ in the blood. We assume no-flux boundary conditions on Γ_B for all remaining variables. We also assume no-flux conditions on Γ_A for all variables (except M and T, as stated above)

Initial conditions

We take initially , where I₆₀ = 6 × 10⁻⁹ gcm⁻³ in the source of influx of I₆ from the endothelial cells [39], and as in [40],All other cell densities and cytokine concentrations are taken to be zero initially. We prescribe the initial geometry of aneurysm by

Mechanical hyperelasticity model

We set Ω = Ω_M ∪ Γ_M ∪ Ω_A and refer to this region, briefly, as the aortic wall. Aortic wall is assumed to have nonlinear, hyperelastic material properties [27, 41, 42]. Specifically, it is an ‘almost’ incompressible, homogeneous, and isotropic material with energy density function W of the form(14)where β₁ and β₂ are forces that represent elastic coefficients; I_B is the first invariant of the Left Cauchy-Green tensor B (namely I_B = tr(B)), B = F F^T, while F is the deformation gradient tensor. In terms of the displacement vector u = (u_i), we can write , where I is the identity tensor.

The Cauchy stress tensor is given by(15)where p is the hydrostatic pressure. The momentum equation is then(16)where f is the body force and d is the tissue density. In the sequel we assume no body force and neglect tissue inertia, so that(17)

The coefficients β₁ and β₂ are functions of the SMCs which we take to represent the elasticity coefficients. We shall use the following specific relations in Eq (14):(18)The pressure, p, is a function of ECM, which we take to have the form(19)for some parameters β_p, p*.

For boundary conditions, we take(20) (21)where σ_B is the stress tensor from the blood and κ is the mean curvature.

Results

Fig 4 shows the profile of average of cell densities and cytokines concentrations in the first 500 days, and Fig 5 shows how the aneurysm deforms at day 500.

[Figure omitted. See PDF.]

Fig 4. Simulation results over 500 days.

The parameter used are as in Tables 2 and 3 with I₆₀ = 6 × 10⁻⁹ g/ml.

https://doi.org/10.1371/journal.pone.0170807.g004

[Figure omitted. See PDF.]

Fig 5.

Aneurysm deformation at day 500 with initial deformation shown in (A). The parameters used are as in Tables 2 and 3 with (B) I₆₀ = 6 × 10⁻⁹ g/ml, and (C) I₆₀ = 6 × 10⁻⁸ g/ml. The diameter grows from 2 cm to 2.5 cm in case (B), and to 4.5 cm in case (C).

https://doi.org/10.1371/journal.pone.0170807.g005

From Fig 4, we see that, as AAA progresses, the average densities of SMCs and ECM decrease. The concentration of the inflammatory cytokines TNF-α, IFN-γ and MCP-1 are increasing. On the other hand, the average concentration of MMP is increasing only for the first 200 days, and it then decreases as the concentration of TIMP keeps growing. The average concentration of IL-6 in the tissue decreases and nearly stabilizes after 500 days. Fig 5 shows growth of the aorta bulge after 500 days: from diameter of 2 cm to diameter of nearly 2.5 cm when I₆₀ = 6 × 10⁻⁹ g/ml and nearly 4.5 cm when I₆₀ = 6 × 10⁻⁸ g/ml.

In Figs 4 and 5, the initial values of the pro-inflammatory cytokines were taken to be below their steady state. On the other hand the SMCs and ECM densities were taken to be above their steady state. Other choices of the initial conditions yield similar profiles (not shown here).

If at the initial screening for AAA we measure both the aorta diameter, R₀, and the concentration of IL-6 in the blood, I₆₀, then by simulating the model with these values of R(0) = R₀ and I₆₀ (in Eq (13)), we can predict the AAA diameter at any future time T. Figs 6 and 7 simulate the cases T = 300 and T = 500 days. Thus, for every pair of I₆₀ (on the horizontal axis) and R₀ (on the vertical axis) the color on the column shows the diameter corresponding to the point (I₆₀, R₀) at day T = 300 (upper) and T = 500 (lower).

[Figure omitted. See PDF.]

Fig 6. AAA growth.

The x-axis scales I₆₀ from 6 × 10⁻⁹ g/ml to 6 × 10⁻⁸ g/ml; the y-axis scales R₀ from 2 cm to 3 cm; The color represents the diameter of aortic bulge at 300 day.

https://doi.org/10.1371/journal.pone.0170807.g006

[Figure omitted. See PDF.]

Fig 7. AAA growth.

The x-axis scales I₆₀ from 6 × 10⁻⁹ g/ml to 6 × 10⁻⁸ g/ml; the y-axis scales R₀ from 2 cm to 3 cm; The color represents the diameter of aortic bulge at 500 day.

https://doi.org/10.1371/journal.pone.0170807.g007

Conclusions

Abdominal aortic aneurysm is a localized enlargement of the abdominal aorta, which may lead to rupture of the aorta. The disease is asymptomatic until rupture, which is nearly always fatal. The risk of rupture associated with the increased diameter of the aorta varies among people. Hence it is difficult to determine, on the basis of just the current diameter (R₀) of a patient’s aortic wall, how closely this patient’s abdominal aorta should be monitored. In order to make such a determination, we need address the following questions:

* how fast the diameter will grow over time?

* how does the risk of rupture depend on the growing diameter?

Several non-invasive prognostic biomarkers for AAA have been suggested [13, 14]; specifically, IL-6 [11, 12].

The present paper addresses the first question with a mathematical model: Given R₀ and the serum concentration of IL-6, the model predicts the dynamical growth of the diameter, R(t), for any future time t. The mathematical model includes the basic biology underlying AAA formation as a disease that degrades the elastic strength of the aortic wall. Simulations of the model can show how the initial diameter R₀ and the serum concentration of IL-6 determine the diameter R after a given period of time T; This is done in Fig 6 in the cases of T = 300 days and T = 500 days. Thus, the level of IL-6 concentration in the blood can be used to suggest how soon to schedule the next screening of the patient’s abdominal aorta.

The model is represented by a system of partial differential equations, based on some simplified assumptions. In particular, the model assumes that the disease is associated with initial inflammation, and that the aortic wall is a ‘pure’ hyperelastic material. Furthermore, some of the parameters were estimated crudely, for the lack of available data. The authors plan to follow up this study with patients data, which are unavailable at this time. The patients will be followed every 6 months with CT scan and measurement of their aneurysm, and of course any patient who comes for surgery or rupture will be recorded. When patients data become available, the model could then be refined, in particular by adjusting some parameters for specific groups of patients. The model could then be used as a predictive tool for monitoring the enlargement of the abdominal aorta.

Parameter estimates

All the parameter values are listed in Tables 2 and 3. Most parameters are taken from the literature while some are estimated below. A few remaining parameters are fitted.

[Figure omitted. See PDF.]

Table 2. Parameters’ description and value.

https://doi.org/10.1371/journal.pone.0170807.t002

[Figure omitted. See PDF.]

Table 3. Parameters’ description and value.

https://doi.org/10.1371/journal.pone.0170807.t003

λ_I6 M and λ_I6 S

According to [16], 10⁶ macrophages produced 10 ng/ml of IL-6. Assuming steady state, we haveSince d_I6 = 0.173 day⁻¹ [43], we get λ_I6 M = 1.73 × 10⁻⁶ day⁻¹.

According to [18], 4 × 10⁴ cells/ml SMCs produced 4 ng/ml of IL-6. Assuming steady state, we haveSince d_I6 = 0.173 day⁻¹, we get λ_I6 S = 1.73 × 10⁻⁵ day⁻¹.

β₁₀, β₂₀, k₁ and k₂

By [42], β₁ = 0.022 N/mm² and β₁₀ = 0.174 N/mm². We assume that in aneurysm S is approximately half its normal density, i.e., . Then by Eq (17),which gives k₁ = 30.4 N/cm².

Similarly, by [42], β₂ = 1.461 N/mm² and β₂₀ = 1.881 N/mm² so thathence k₂ = 84 N/cm².

β_p

By [44], the peak wall stress measured for normal people and AAA patients were 30-44 N/cm² and 28-37 N/cm², respectively. We assume that ECM in AAA patients is reduced to 75% of its value in healthy individuals. Hencewhich implies, by Eq (18), that β_p = 18 N/cm².

1.1 σ_B

We assume σ_B = −p_B I, where p_B is average central aortic systolic pressure (100-120 mm Hg) [45]; we take p_B = 110 mm Hg.

Numerical methods

Weak formulation

In this subsection, we rewrite in compact form the coupled biological and mechanical systems. We denote by Z the vector of the densities or concentrations of cells and cytokines, i.e. Z = (P, M, T, S, I₆, I₁₀, I₁₂, T_α, I_γ, Q, Q_r, ρ)^T. The entries of Z satisfy advection-diffusion equations with nonlinear source term, as seen from Eqs (1)–(12):(22)where the velocity v is defined by v = D_t u, where D_t denotes the material derivative,In the equation for Z, F(Z) is a source term and D_Z is a diagonal matrix with the diagonal terms for the diffusion coefficients of the biological species.

Some of the components of Z live only in either Ω_M(t) or Ω_A(t). Since there is a jump for Z at the media/adventitia interface, Γ_M(t), we split Z into two partsThe boundary conditions can then be written in the formfor an appropriate coefficient , some of its component being zero.

The governing equation for the displacement u can be written as(23)where the stress σ is is defined in terms of the displacement bywith the boundary conditions

Multiplying the coupled System (22) and (23) by arbitrary test functions (W_M, W_A,w), and performing integration by parts using the boundary conditions, we get the weak formulation:(24)Here we used (⋅, ⋅)_D and ⟨⋅, ⋅⟩_∂D to denote the L² inner product on the domain D and on the boundary ∂D, respectively.

Lagrangian description for hyperelastic equation

Ω_A(t), Ω_M(t) and their boundaries are moving in time. We set Ω(t) = Ω_M(t) ∪ Ω_A(t) ∪ Γ_M(t) and , , . To characterize the motion of the media and adventitia, we introduce a flow map for the position of the particle at time t.

Then each function can be described in Lagrangian coordinates, namely for each function f(x, t), its Lagrangian description isWe also introduce the following variables in Lagrangian coordinates: The displacement , the deformation tensor , and its determinant .

By change of variables, we can write the weak form for the equilibrium equation for w in Lagrangian coordinates,(25)where

As for the surface tension term, we can also write it in the Lagrangian coordinates as follows:(26)

By substituting this expression into the hyperelasitcity equation in Lagrangian coordinates Eq (25), we can get a nonlinear equation:(27)

Finite element discretization on the moving mesh

In order to deal with the free boundary and the jump at the interface, we use the moving mesh method.

The initial domains and are discretized by triangulations, denoted by and , respectively. We consider the case that and are matching on the interface . We use piecewise linear and continuous function to approximate the displacement. Then the triangulation for the current domain can be obtained by moving the grid nodes of according to the flow map .

We next introduce the following time discretization for the material derivativeand forwhere we use the superscript n to indicate the function at time nk, where k is the time step. Note that no interpolation is needed for evaluating the material derivative D_t,k at grid points, since the grid points are moving according to the flow map.

For the space discretization, we use piecewise linear and continuous function spaces and to approximate the solutions. Let , and be the bases of and , respectively, where N, N_M and N_A are the numbers of basis functions of the corresponding spaces. The approximation solutions at time nk can be written as summations with coefficients to be determined by the discrete system, corresponding to the System (24),(28)for any and . Here we use the explicit discretization of the coefficients for the elasticity, and the source terms for the biological species. We denote the first two equations in Eulerian coordinates of the System (28) by

By this discretization, the elastic equation and the biological equation are decoupled. Given the solutions in last time step, we use the Newton’s nonlinear iteration method to solve the system a_L = 0 in the initial domain. Then we update the triangulations for the current domain by the approximate flow map . Following the standard technique of finite element method, we then solve the variational problem a_E = 0 for the biological species. We thus get the solution of current time step.

Author Contributions

1. Conceptualization: WH AF MG JX.

2. Funding acquisition: DZ.

3. Investigation: WH AF MG.

4. Methodology: SW SG JX.

5. Project administration: WH JX.

6. Resources: MG.

7. Software: DZ.

8. Validation: WH AF MG.

9. Visualization: WH SW SG.

10. Writing – original draft: WH AF SW SG.

11. Writing – review & editing: WH AF SW SG AF MG JX DZ.

Citation: Hao W, Gong S, Wu S, Xu J, Go MR, Friedman A, et al. (2017) A mathematical model of aortic aneurysm formation. PLoS ONE 12(2): e0170807. https://doi.org/10.1371/journal.pone.0170807

References

1. McGregor JC, Pollock JG, Anton HC. The diagnosis and assessment of abdominal aortic aneurysms by ultrasonography. Ann R Coll Surg Engl. 1976;58(5):388–392. pmid:970885

2. Filardo G, Powell JT, Martinez MA, Ballard DJ. Surgery for small asymptomatic abdominal aortic aneurysms. Cochrane Database Syst Rev. 2012;3:CD001835. pmid:22419281

3. Veith FJ, Ohki T, Lipsitz EC, Suggs WD, Cynamon J. Treatment of ruptured abdominal aneurysms with stent grafts: a new gold standard? Semin Vasc Surg. 2003;16(2):171–175. pmid:12920689

4. No authors listed. Mortality results for randomised controlled trial of early elective surgery or ultrasonographic surveillance for small abdominal aortic aneurysms.. The UK Small Aneurysm Trial Participants Lancet. 1998;352(9141):1649–1655. pmid:9853436

5. Bush RL, Johnson ML, Collins TC, Henderson WG, Khuri SF, Yu HJ, et al. Open versus endovascular abdominal aortic aneurysm repair in VA hospitals. J Am Coll Surg. 2006;202(4):577–587. pmid:16571424

6. Nicholls SC, Gardner JB, Meissner MH, Johansen HK. Rupture in small abdominal aortic aneurysms. J Vasc Surg. 1998;28(5):884–888. pmid:9808857

7. Grytsan A, Watton PN, Holzapfel GA. A thick-walled fluid-solid-growth model of abdominal aortic aneurysm evolution: application to a patient-specific geometry. J Biomech Eng. 2015;137(3). pmid:25473877

8. Schmid H, Grytsan A, Poshtan E, Watton PN, Itskov M. Influence of differing material properties in media and adventitia on arterial adaptation–application to aneurysm formation and rupture. Comput Methods Biomech Biomed Engin. 2013;16(1):33–53. pmid:22149119

9. Volokh KY. Modeling failure of soft anisotropic materials with application to arteries. J Mech Behav Biomed Mater. 2011;4(8):1582–1594. pmid:22098860

10. Watton PN, Hill NA, Heil M. A mathematical model for the growth of the abdominal aortic aneurysm. Biomech Model Mechanobiol. 2004;3(2):98–113. pmid:15452732

11. Jones KG, Brull DJ, Brown LC, Sian M, Greenhalgh RM, Humphries SE, et al. Interleukin-6 (IL-6) and the prognosis of abdominal aortic aneurysms. Circulation. 2001;103(18):2260–2265. pmid:11342474

12. Takagi H, Watanabe T, Mizuno Y, Kawai N, Umemoto T. Circulating interleukin-6 levels are associated with abdominal aortic aneurysm presence: a meta-analysis and meta-regression of case-control studies. Ann Vasc Surg. 2014;28(8):1913–1922. pmid:25011090

13. Moris D, Mantonakis E, Avgerinos E, Makris M, Bakoyiannis C, Pikoulis E, et al. Novel biomarkers of abdominal aortic aneurysm disease: identifying gaps and dispelling misperceptions. Biomed Res Int. 2014;2014:925840. pmid:24967416

14. Golledge J, Tsao PS, Dalman RL, Norman PE. Circulating markers of abdominal aortic aneurysm presence and progression. Circulation. 2008;118(23):2382–2392. pmid:19047592

15. Abdul-Hussien H, Hanemaaijer R, Kleemann R, Verhaaren BF, van Bockel JH, Lindeman JH. The pathophysiology of abdominal aortic aneurysm growth: corresponding and discordant inflammatory and proteolytic processes in abdominal aortic and popliteal artery aneurysms. J Vasc Surg. 2010;51(6):1479–1487. pmid:20488324

16. Wang Q, Ni H, Lan L, Wei X, Xiang R, Wang Y. Fra-1 protooncogene regulates IL-6 expression in macrophages and promotes the generation of M2d macrophages. Cell Res. 2010;20(6):701–712. pmid:20386569

17. Volk T, Hensel M, Schuster H, Kox WJ. Secretion of MCP-1 and IL-6 by cytokine stimulated production of reactive oxygen species in endothelial cells. Mol Cell Biochem. 2000;206(1-2):105–112. pmid:10839200

18. Loppnow H, Bil R, Hirt S, Schonbeck U, Herzberg M, Werdan K, et al. Platelet-derived interleukin-1 induces cytokine production, but not proliferation of human vascular smooth muscle cells. Blood. 1998;91(1):134–141. pmid:9414277

19. Wang Q, Ren J, Morgan S, Liu Z, Dou C, Liu B. Monocyte chemoattractant protein-1 (MCP-1) regulates macrophage cytotoxicity in abdominal aortic aneurysm. PLoS ONE. 2014;9(3):e92053. pmid:24632850

20. Koch AE, Kunkel SL, Pearce WH, Shah MR, Parikh D, Evanoff HL, et al. Enhanced production of the chemotactic cytokines interleukin-8 and monocyte chemoattractant protein-1 in human abdominal aortic aneurysms. Am J Pathol. 1993;142(5):1423–1431. pmid:8494046

21. Middleton RK, Bown MJ, Lloyd GM, Jones JL, London NJ, Sayers RD. Characterisation of Interleukin-8 and monocyte chemoattractant protein-1 expression within the abdominal aortic aneurysm and their association with mural inflammation. Eur J Vasc Endovasc Surg. 2009;37(1):46–55. pmid:19008128

22. Tieu BC, Lee C, Sun H, Lejeune W, Recinos A, Ju X, et al. An adventitial IL-6/MCP1 amplification loop accelerates macrophage-mediated vascular inflammation leading to aortic dissection in mice. J Clin Invest. 2009;119(12):3637–3651. pmid:19920349

23. Xiong W, Zhao Y, Prall A, Greiner TC, Baxter BT. Key roles of CD4+ T cells and IFN-gamma in the development of abdominal aortic aneurysms in a murine model. J Immunol. 2004;172(4):2607–2612. pmid:14764734

24. Xiong W, MacTaggart J, Knispel R, Worth J, Persidsky Y, Baxter BT. Blocking TNF-alpha attenuates aneurysm formation in a murine model. J Immunol. 2009;183(4):2741–2746. pmid:19620291

25. Elmore JR, Keister BF, Franklin DP, Youkey JR, Carey DJ. Expression of matrix metalloproteinases and TIMPs in human abdominal aortic aneurysms. Ann Vasc Surg. 1998;12(3):221–228. pmid:9588507

26. Hasan D, Chalouhi N, Jabbour P, Hashimoto T. Macrophage imbalance (M1 vs. M2) and upregulation of mast cells in wall of ruptured human cerebral aneurysms: preliminary results. J Neuroinflammation. 2012;9:222. pmid:22999528

27. Vorp DA, Vande Geest JP. Biomechanical determinants of abdominal aortic aneurysm rupture. Arterioscler Thromb Vasc Biol. 2005;25(8):1558–1566. pmid:16055757

28. Annabi B, Shedid D, Ghosn P, Kenigsberg RL, Desrosiers RR, Bojanowski MW, et al. Differential regulation of matrix metalloproteinase activities in abdominal aortic aneurysms. J Vasc Surg. 2002;35(3):539–546. pmid:11877705

29. Tian X, Fan J, Yu M, Zhao Y, Fang Y, Bai S, et al. Adipose stem cells promote smooth muscle cells to secrete elastin in rat abdominal aortic aneurysm. PLoS ONE. 2014;9(9):e108105. pmid:25243605

30. Rahim SS, Khan N, Boddupalli CS, Hasnain SE, Mukhopadhyay S. Interleukin-10 (IL-10) mediated suppression of IL-12 production in RAW 264.7 cells also involves c-rel transcription factor. Immunology. 2005;114(3):313–321. pmid:15720433

31. Parameswaran N, Patial S. Tumor necrosis factor-alpha signaling in macrophages. Crit Rev Eukaryot Gene Expr. 2010;20(2):87–103. pmid:21133840

32. Baer M, Dillner A, Schwartz RC, Sedon C, Nedospasov S, Johnson PF. Tumor necrosis factor alpha transcription in macrophages is attenuated by an autocrine factor that preferentially induces NF-kappaB p50. Mol Cell Biol. 1998;18(10):5678–5689. pmid:9742085

33. Teixeira LK, Fonseca BP, Vieira-de Abreu A, Barboza BA, Robbs BK, Bozza PT, et al. IFN-gamma production by CD8+ T cells depends on NFAT1 transcription factor and regulates Th differentiation. J Immunol. 2005;175(9):5931–5939. pmid:16237086

34. Patel MI, Melrose J, Ghosh P, Appleberg M. Increased synthesis of matrix metalloproteinases by aortic smooth muscle cells is implicated in the etiopathogenesis of abdominal aortic aneurysms. J Vasc Surg. 1996;24(1):82–92. pmid:8691532

35. Galle C, Schandene L, Stordeur P, Peignois Y, Ferreira J, Wautrecht JC, et al. Predominance of type 1 CD4+ T cells in human abdominal aortic aneurysm. Clin Exp Immunol. 2005;142(3):519–527. pmid:16297165

36. Maiellaro K, Taylor WR. The role of the adventitia in vascular inflammation. Cardiovasc Res. 2007;75(4):640–648. pmid:17662969

37. Majesky MW, Dong XR, Hoglund V, Mahoney WM, Daum G. The adventitia: a dynamic interface containing resident progenitor cells. Arterioscler Thromb Vasc Biol. 2011;31(7):1530–1539. pmid:21677296

38. Kim Y, Friedman A. Interaction of tumor with its micro-environment: A mathematical model. Bull Math Biol. 2010;72(5):1029–1068. pmid:19908100

39. Yokoyama U, Ishiwata R, Jin MH, Kato Y, Suzuki O, et al. Inhibition of EP4 signaling attenuates aortic aneurysm formation. PLoS ONE. 2012;7(5):e36724. pmid:22570740

40. Hao W, Friedman A. The LDL-HDL profile determines the risk of atherosclerosis: a mathematical model. PLoS ONE. 2014;9(3):e90497. pmid:24621857

41. Raghavan ML, Vorp DA. Toward a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability. J Biomech. 2000;33(4):475–482. pmid:10768396

42. Zelaya JE, Goenezen S, Dargon PT, Azarbal AF, Rugonyi S. Improving the efficiency of abdominal aortic aneurysm wall stress computations. PLoS ONE. 2014;9(7):e101353. pmid:25007052

43. Lu ZY, Brailly H, Wijdenes J, Bataille R, Rossi JF, Klein B. Measurement of whole body interleukin-6 (IL-6) production: prediction of the efficacy of anti-IL-6 treatments. Blood. 1995;86(8):3123–3131. pmid:7579407

44. Wang DH, Makaroun MS, Webster MW, Vorp DA. Effect of intraluminal thrombus on wall stress in patient-specific models of abdominal aortic aneurysm. J Vasc Surg. 2002;36(3):598–604. pmid:12218961

45. Sule AA, Hwang TH, Chin TJ. Very high central aortic systolic pressures in a young hypertensive patient on telmisartan: Is central aortic systolic pressure associated with white coat hypertension? Int J Angiol. 2010;19(4):e132–134. pmid:22479144

46. Friedman A, Hao W. A Mathematical Model of Atherosclerosis with Reverse Cholesterol Transport and Associated Risk Factors. Bull Math Biol. 2014;. pmid:25205457

47. Hao W, Crouser ED, Friedman A. Mathematical model of sarcoidosis. Proc Natl Acad Sci USA. 2014;. pmid:25349384

48. Hao W, Rovin BH, Friedman A. Mathematical model of renal interstitial fibrosis. Proc Natl Acad Sci USA. 2014;111(39):14193–14198. pmid:25225370