Growth rate and rupture rate of unruptured intracranial aneurysms: a population approach
© Jou and Mawad; licensee BioMed Central Ltd. 2009
Received: 30 January 2009
Accepted: 18 June 2009
Published: 18 June 2009
Understanding aneurysm growth rate allows us to predict not only the current rupture risk, but also accumulated rupture risk in the future. However, determining growth rate of unruptured intracranial aneurysms often requires follow-up of patients for a long period of time so that significant growth can be observed and measured. We investigate a relationship between growth rate and rupture rate and develop a theoretical model that can predict average behavior of unruptured intracranial aneurysms based on existing clinical data.
A mathematical model is developed that links growth rate and rupture rate. This model assumes a stable aneurysm size distribution so the number of aneurysm ruptures is balanced by the growth of aneurysms. Annual growth rates and growth profiles are calculated from a hypothetical size distribution and data from a previous clinical study.
Our model predicts a growth rate of 0.34–1.63 mm/yr for three different growth models when the rupture rate at 10 mm is 1%. The growth rate is 0.56–0.65 mm/yr if annual rupture rate averaged over all aneurysm sizes is assumed to be 2%. The peak of aneurysm size distribution coincides with a period of slow growth between 5 mm and 8 mm.
This mathematical model can be used to predict aneurysm growth rate, and the results are consistent with previous clinical studies. Predictions from both hypothetical and clinical cases agree very well. This model explains why some aneurysms may grow into a stable size and remain so without rupture.
Rupture of an intracranial aneurysm is devastating; however, prevention remains a challenge since most aneurysms are asymptomatic before rupture. Rupture risk is difficult to assess even if an aneurysm is detected. Each aneurysm behaves differently, so the only way to study this cerebrovascular disease is through its natural history. Natural history of intracranial aneurysms can help us to understand the course of aneurysm development and allow proper assessment of aneurysm rupture risk. Ultimately, patients can benefit from rational decisions regarding their treatments based on their individual risk.
A large clinical study on unruptured intracranial aneurysms, the International Study of Unruptured Intracranial Aneurysms (ISUIA), has studied more than 4000 aneurysm patients and raised serious interests on risk factors and rupture rate of unruptured aneurysms [1, 2]. They found that rupture rate increased with aneurysm size in both retrospective and prospective studies, and rupture rate was higher for patients with a history of subarachnoid hemorrhage. Rupture rate assesses immediate rupture risk, and growth rate allows prediction of accumulated rupture risk and evaluation of treatment timing. While large studies of this kind have provided critical information on rupture rate and risk factors, few studies have analyzed aneurysm growth rate in the same rigorous manner. Most studies analyzed angiographical images of aneurysms retrospectively and had very short follow-up periods [3, 4]. Ethical concerns and small growth rates are among the major reasons that prohibit a large prospective clinical study on growth rate. It often requires a few years for aneurysms to grow at least 1mm larger so that growth can be appreciated on angiograms. Intra- and inter-observers variation may also seriously affect the accuracy of growth rate [6, 7].
Koffijberg et al. found that intracranial aneurysms did not grow constantly and prediction of their behavior might be extremely difficult . Markov chain models have been used to study the natural history of unruptured intracranial aneurysms and benefit of prophylactic treatment . These reports all relied on estimated aneurysm initiation rate, growth rate, and rupture rate [9, 10].
Rupture rate and growth rate are two seemingly unrelated concepts. Will it be possible to determine the growth rate of unruptured intracranial aneurysms from rupture rates in previous clinical studies? In this report, we explore and investigate the possible mathematical relationship between rupture rate and growth rate. We employ a different approach and intend to derive aneurysm growth rate from clinical data. A simple mathematical model with only a few parameters can give us a reasonable description of average aneurysm behavior without organizing large clinical studies to understand the natural history of unruptured intracranial aneurysms.
where t is the time. The first term on the left hand side describes the increase of the number of aneurysms with time, and the second term on the left hand side is the net change of the number of aneurysms due to growth. A number of aneurysms will grow into a bigger size (>r) and a number of smaller aneurysms (<r) will grow to the size r, and the net result will be their difference, . Net change of these aneurysms is balanced by the number of aneurysm ruptures on the right hand side. Equation (1) essentially describes the conservation of unruptured aneurysms at the size r. Once an aneurysm bleeds, it will become a ruptured aneurysm, and its behavior is different from unruptured aneurysms. Ruptured aneurysms need to be dealt with differently since it involves not only the aneurysms that rupture now, but also the aneurysms that ruptured in the past; patients suffered from ruptured aneurysms may or may not survive. If we consider only the cases that rupture within a short time period Δt, then the distribution of ruptured aneurysms is , and the average annual rupture rate is . h will be used for this average annual rupture rate over al aneurysm sizes, and h is actually a weighted average of rupture rate f(r).
One can easily calculate the growth rate g(r) if the aneurysm distribution, N(r), aneurysm rupture rate, f(r), and increase of total number of aneurysms, , are known. Because N(r) is an aneurysm size distribution, integration of N(r) over all aneurysm sizes will be the unity, that is,
Equation (3) shows that initiation of new aneurysms CN(r = 0) g(r = 0) is tied to the average rupture rate (h) and temporal change of patient population ( ). Therefore, one cannot assume both the growth rate and rupture rate at the same time, or Equation (3) will be invalid. Initiation and growth of aneurysms are important in maintaining a stable aneurysm size distribution. The change of C also reflects the deaths of these patients from other diseases or life events.
Although g(r) is the percentage of aneurysms that grow, it is related to aneurysm size growth. Considering that the aneurysms at the size r is N(r), the average size of these aneurysms after a time period Δt is . CN(1-gΔt) is the number of aneurysms remaining at the size r, and CNgΔt is the number of aneurysms growing into the size r + Δr. So the change of average size for these aneurysms is g(r)ΔtΔr, which describes the increase of mean aneurysm size within a given time period. For example, there are 100 aneurysms in the 5 mm range, and the mean aneurysm size for these aneurysms is 5 mm. 20% of these aneurysms grow into the 6 mm range, then the mean aneurysm size for these 100 aneurysms is 5.2 mm (=(80*5+20*6)/100). Thus, 20% of aneurysms that grow lead to a growth rate of 0.2 mm/year. This argument works only when the time period (Δt) is small enough that no aneurysms grow into the 7 mm range. The average time needed for the aneurysms size r to grow Δr is and the time that it takes an aneurysm to grow from r = 0 to r = R is
These equations establish a mathematical model for our analysis. In the following, we will analyze average aneurysm growth behavior based on these relationships.
Our model basically provides a link between growth rate and rupture rate; however, the aneurysm size distribution is a pre-requisite to the calculation of growth rate. Thus, we assume a profile for the aneurysm size distribution: N(r) = ar/(1+kr2)2. This is a fitted curve of the data in . In the profile, k controls the peak of the profile, r m , and a is determined by An example for r m = 5.77 is shown in Figure 1. There is only one parameter in the size distribution so determining the size distribution is easy if the mode (r m ) can be found from clinical data. Nonetheless, the exact mathematical formula is not used in the simulations, and only the numerical values at each size group are used.
In the following, we discuss three different cases for which we derive the growth rates from rupture rates. In these cases, the rupture rate f is proportional to various power of aneurysm size.
A) Rupture rate proportional to aneurysm size
Because growth rate depends on the change of total number of aneurysms (dC/dt), a sensitivity analysis is conducted for different patient population growth rate (0.5%, 1%, and 2%). The patient population growth, dC/Cdt, is different from aneurysm growth rate. Aneurysm growth rate describes how fast an aneurysm grows into a larger size, and the patient population growth reflects the change of total number of aneurysm patients. Suppose that the aneurysm prevalence is 2% in the U.S., then aneurysm patient population growth rate will be 1%, which is the same as the general population growth rate in the U.S. Aneurysm patient population growth needs to be the same as the general population growth so the prevalence of aneurysms is a constant.
The sensitivity of aneurysm growth rate to patient population growth rate is shown in Figure 2B. The growth rate is 0.32 mm/yr, 0.37 mm/yr, 0.47 mm/yr at 20 mm for the patient population growth rates of 0.5%, 1%, and 2%, respectively. Patient population growth rate of 1% is used in all the following analyses unless stated otherwise.
B) Rupture rate proportional to the square of aneurysm size
In this case, rupture rate is proportional to the square of aneurysm size, so f(r) = αr2. α is adjusted so that annual rupture rate is 1% for 10 mm aneurysms. For spherical aneurysms, this rupture rate is proportional to aneurysm surface area.
(C) Rupture rate proportional to the cube of aneurysm size
A previous study hypothesized that the rupture rate of an aneurysm is proportional to the aneurysm volume . They argued that vessel thickness is proportional to vessel size and weakest sites of an aneurysm are proportional to the total cell mass (thickness × surface area), so rupture rate is proportional to the aneurysm volume. In this case, rupture rate is assumed to be proportional to the cube of aneurysm size: f(r) = αr3. α is adjusted so that rupture rate is 1% for 10 mm aneurysms. For spherical aneurysms, this rupture rate is proportional to aneurysm volume. Aneurysm growth in this case is obviously faster than previous two cases.
Growth curves in Figure 3B present the time it takes for an aneurysm to grow to a certain size. The analysis shows that it needs 20 and 9 years for an aneurysm to grow to 10 mm for Case B and C, respectively. There are local minimum and maximum growth rates for aneurysms from Figure 3C. For Case C, the growth rate varies from 0.5 mm/yr at 5 mm to 4 mm/yr at 30 mm.
Differences in the mean growth rate (mm/yr) for three different cases.
f 10 mm = 1%/yr
h = 2%/yr
Differences in the time (yrs) for an aneurysm to grow to 10 mm after initiation.
f 10 mm = 1%/yr
h = 2%/yr
Obtaining size distribution of unruptured aneurysms needs at least one angiographical examination for a large number of aneurysm patients . Determining rupture rate requires additional follow-up on these patients either by interviews or questionnaires for their current conditions [1, 2]. Even then, the results represent only the average rupture rate of a specific range of initial aneurysm size, assuming that there is no significant growth during the study period. At least two angiographical examinations are required to determine the growth rate for each patient [3–5]. Since only a small number of aneurysms actually grow, the effort leading to determination of aneurysm growth rate is enormous. We have derived a model that describes the relationship between growth rate and rupture rate of unruptured intracranial aneurysms. Based on the assumption that there is a stable aneurysm size distribution, we learn that one cannot assume both growth rate and rupture rate in an analysis, or the model is invalid.
Our model does not deal with any specific aneurysm or aneurysms at a particular location, but it predicts average behavior of entire aneurysm population. Based on our analysis, aneurysm size is no longer the only consideration for rupture risk. For Case A, an 8 mm aneurysm growing at 0.22 mm/yr behaves just like any other 8 mm aneurysms. However, a 5 mm aneurysm with a growth rate of 0.3 mm/yr is growing 50% faster than average 5 mm aneurysms, and this increase in growth rate may out-weight aneurysm size in assessment of rupture risk. Therefore, the 8 mm aneurysm may have a smaller rupture risk than a 5 mm aneurysm in this case; it becomes important when both aneurysms are found in the same individual and one of them needs to be treated immediately.
Clinical experiences show that intracranial aneurysms either grow to rupture after development or grow into a certain size without rupture. Our model explains why many aneurysms grow into a certain size and remain stable without rupture for a long period of time. Initial aneurysm development is often followed by rapid growth at early stage, and then the growth is slowed down between 5 and 8 mm, as depicted in Figure 6. This slow growth period coincides with the peak of unruptured aneurysm size distribution in Figure 1, implying that only a very small percentage of aneurysms sized between 5–8 mm are growing. Thus, more aneurysms are found at this size range and their growth is slow. Aneurysm growth regains its pace when aneurysm size is greater than 10 mm, and growth accelerates until it is larger than 30 mm. Fast growth and increasing rupture rate render fewer aneurysms at larger sizes. Note that the end of this slow growth period between aneurysm size of 5–8 mm happens to be near the critical aneurysm size , and this might also indicate a change of aneurysm growth behavior and beginning of an increasing rupture risk. Nevertheless, our model shows the same behavior both for a hypothetical model and from clinical data .
Growth rate of unruptured aneurysms can be related to rupture rate. Therefore, we can calculate growth rate given a specific rupture rate. Our model gives a growth rate of 0.34–1.63 mm/yr when f 10 mm = 1%/yr. However, our model shows a narrow range of growth rate (0.56–0.65 mm/yr) for h = 2%/yr. By assuming a 2% annual rupture rate averaged over entire aneurysm size range, the difference among three cases is smaller for small aneurysms, and their difference would be difficult to quantify based on clinical data. Juvela et al. reported a 0.31 mm/yr growth rate in 45% of patients with a mean follow-up of 19 years ; however, nearly 90% of patients in their study had subarachnoid hemorrhage at other locations before the study began. Phan et al. reported that 7% of unruptured intracranial aneurysms grew during an average follow-up period of 50-month, and all growing aneurysms were greater than 9 mm . We learn from Figure 6 that unruptured aneurysms smaller than 9 mm grow 0.18 mm/yr (or 0.75 mm in 50 months), and few aneurysms can grow 1 mm larger within 50 months, consistent with Phan's finding. Therefore, observation of growth of small aneurysms is infrequent , and only a long-term study can reveal significant growth .
Because of the connection between rupture rate and growth rate, high rupture rate can only be explained by faster aneurysm growth. Among the cases that we have studied, Case C, in which rupture rate is proportional to the cube of aneurysm size, gives the fastest growth. However, it generates an unusual ruptured aneurysm size distribution in Figure 3A. Mean ruptured aneurysm size is larger than 20 mm, and this clearly differs from previous data. Aneurysm growth of 10mm3/yr suggested in another study is, however, too restrictive since it yields little growth for large and giant aneurysms . Therefore, rupture rate probably will not increase as fast as in Case C.
Rupture risk involves a short-term rupture rate and an increase of rupture rate associated with aneurysm growth. Case B gives a 75% survival rate at 23 years, similar to 70% survival rate at 25 years reported by Juvela et al. . However, 84% of aneurysms survive without rupture after 25 years when the growth is not considered, and rupture risk is actually underestimated by 35%. Therefore, the rupture rate derived from clinical survival rate may need to be lowered to reflect the growth of aneurysms. Since the surface area represents the amount of tissue in an aneurysm, the assumption that rupture rate is proportional to the square of aneurysm size (Case B) seems to be logical as long as wall thinning is not significant.
Our theoretical model has several advantages. First, the model parameters can be determined from clinical data. Unruptured aneurysm distribution is available from many studies [1, 2, 11, 15, 18]. The parameter α can be computed from average rupture rate. The resulting model then predicts average behavior of unruptured aneurysms at any given size. Second, we can compute the growth rate from rupture rate easily. Growth of unruptured aneurysms is very difficult to study. Very few aneurysms grow, and they grow very slowly when they do grow. Therefore, our understanding of aneurysm growth often relies on a small number of cases [3–5]. Our model does help us in understanding general growth behavior of unruptured aneurysms.
Limitations exist in our model. We focus on gross behavior of unruptured aneurysms and ignore individual risk factors that might ultimately dictate rupture. The growth rate g and rupture rate f are averaged over all risk factors, so our model can be used to predict only the average behavior of unruptured aneurysms. The biological ages of these aneurysms or patients are not explicitly represented in the model. However, inclusion of biological age of an aneurysm is possible , but the final model may be so complex that model parameters cannot be obtained or verified from existing data. Since patient age is considered in every aspect of a clinical decision, the benefit from inclusion of patient age explicitly may be minimal.
We have formulated a theory by which growth rate of unruptured intracranial aneurysms can be determined from aneurysm size distribution and rupture rate. This mathematical model can be used to predict aneurysm growth rate, and the results are consistent with previous clinical studies. Aneurysm behaves similarly under hypothetical conditions and from clinical date. This model also explains why some aneurysms may grow into a certain size and remain so without rupture. Rupture risk may also be underestimated if aneurysm growth is not accounted for.
LDJ was partially supported by a grant from the Brain Aneurysm Foundation.
- ISUIA-Investigators: Unruptured Intracranial Aneurysms: Risks of Rupture and Risks of Surgical Intervention. N Engl J Med 1998, 339: 1725–1733. 10.1056/NEJM199812103392401View ArticleGoogle Scholar
- ISUIA-Investigators: Unruptured Intracranial Aneurysms: Natual History, Clinical Outcome, and Risks of Surgical and Endovascular Treatment. Lancet 2003, 362: 103–110. 10.1016/S0140-6736(03)13860-3View ArticleGoogle Scholar
- Phan TG, Houston J, Brown RD, Wiebers DO, Piepgras DG: Intracranial saccular aneurysm enlargement determined using serial magnetic resonance angiography. J Neurosurg 2002, 97: 1023–1028. 10.3171/jns.2002.97.5.1023View ArticleGoogle Scholar
- Miyazawa N, Akiyama I, Yamagata Z: Risk factors for growth of unruptured intracranial aneurysms: follow-up study by serial 0.5-T magnetic resonance angiography. Neurosurgery 2006, 58: 1047–1053. 10.1227/01.NEU.0000217366.02567.D2View ArticleGoogle Scholar
- Juvela S, Poussa K, Porras M: Factors Affectng Formation and Growth of Intracranial Aneurysms: A Long-term Follow-up Study. Stroke 2001, 32: 485–491.View ArticleGoogle Scholar
- White PM, Teasdale EM, Wardlaw JM, Easton V: Intracranial aneurysms: CT angiography and MR angiography for detection prospective blinded comparison in a large patient cohort. Radiology 2001, 219: 739–749.View ArticleGoogle Scholar
- Dispensa BP, Saloner DA, Acevedo-Bolton G, Archol AS, Jou LD, McCulloch RT, Higashida RT, Hallach VV, Ko NU, Lawton MT, Quinnine NJ, Young WL: Estimation of fusiform intracranial aneurysm growth by serial magnetic resonance imaging. J Magn Reson Imaging 2007, 26: 177–183. 10.1002/jmri.20944View ArticleGoogle Scholar
- Koffijberg H, Buskens E, Algra A, Wermer MJ, Rinkel GJ: Growth Rates of Intracranial Aneurysms: Exploring Constancy. J Neurosurg 2008, 109: 176–185. 10.3171/JNS/2008/109/8/0176View ArticleGoogle Scholar
- Chang HS: Simulation of the Natural History of Cerebral Aneurysms Based on Data From the International Study of Unruptured Intracranial Aneurysms. J Neurosurg 2006, 104: 188–194. 10.3171/jns.2006.104.2.188View ArticleGoogle Scholar
- Yoshimoto Y: A Mathematical Model of the Natural History of Intracranial Aneurysms: Quantification of the Benefit of Prophylactic Treatment. J Neurosurg 2006, 104: 195–200. 10.3171/jns.2006.104.2.195View ArticleGoogle Scholar
- Carter BS, Sheth S, Chang E, Sethl M, Ogilvy CS: Epidemiology of the Size Distribution of Intracranial Bifurcation Aneurysms: Smaller Size of Distal Aneurysms and Increasing Size of Unruptured Aneurysms with Age. Neurosurg 2006, 58: 217–223. 10.1227/01.NEU.0000194639.37803.F8View ArticleGoogle Scholar
- Juvela S, Porras M, Poussa K: Natural History of Unruptured Intracranial Aneurysms: Probability of and Risk Factors for Aneurysm Rupture. J Neurosurg 2000, 93: 379–387. 10.3171/jns.2000.93.3.0379View ArticleGoogle Scholar
- Dickey P, Kailasnath P: The Diameter-Cube Hypothesis: a New Biophyiscal Model of Aneurysm Rupture. Surg Neurol 2002, 58: 166–180. 10.1016/S0090-3019(02)00848-0View ArticleGoogle Scholar
- Rinkel GJ, Djibuti M, Algra A, van Gijn J: Prevalence and risk of rupture of intracranial aneurysms: a systematic review. Stroke 1998, 29: 251–256.View ArticleGoogle Scholar
- Weir B, Disney L, Karrison T: Sizes of Ruptured and Unruptured Aneurysms in Relation to Their Sites and the Ages of Patients. JNeurosurg 2002, 96: 64–70. 10.3171/jns.2002.96.1.0064View ArticleGoogle Scholar
- Ferguson GG, Peerless SJ, Drake CG: Natural history of intra cranial aneurysms. N England J Med 1981, 305: 99.Google Scholar
- Wermer MJ, Schaaf IC, Velthuis BK, Majoie CB, Albrecht KW, Rinkel GJ: Yield of short-term follow-up CT/MR angiography for small aneurysms detected at screening. Stroke 2006, 37: 414–418. 10.1161/01.STR.0000199077.06390.35View ArticleGoogle Scholar
- McCormick WF, Acosta-Rua GJ: The Size of Intrancranial Aneurysms: An Autopsy Study. J Neurosurg 1970, 33: 422–427. 10.3171/jns.1970.33.4.0422View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.