Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2759354

Formats

Article sections

Authors

Related links

Int J Cardiol. Author manuscript; available in PMC 2010 June 12.

Published in final edited form as:

Published online 2008 July 1. doi: 10.1016/j.ijcard.2008.03.087

PMCID: PMC2759354

NIHMSID: NIHMS125602

Hui Zhu, PhD,^{a} Zhaohua Ding, PhD,^{b} Robert N. Piana, MD,^{c} Thomas R. Gehrig, MD,^{d} and Morton H. Friedman, PhD^{a}

Correspondence to: Professor Morton H. Friedman, Department of Biomedical Engineering, 136 Hudson Hall, Box 90281, Duke University, Durham, NC 27708-0281 USA, Fax: 919-6844488; Phone: 919-6605154; Email: ude.ekud@namdeirf.trom

The publisher's final edited version of this article is available at Int J Cardiol

See other articles in PMC that cite the published article.

The nonuniform distribution of atherosclerosis in the human vasculature suggests that local fluid dynamics or wall mechanics may be involved in atherogenesis. Thus certain aspects of vascular geometry, which mediates both fluid dynamics and wall mechanics, might be risk factors for coronary atherosclerosis. Cataloguing the geometry of normal human coronary arteries and its variability is a first step toward identifying specific geometric features that increase vascular susceptibility to the disease.

Images of angiographically normal coronary arteries, including 32 left anterior descending (LAD) and 35 right coronary arteries (RCA), were acquired by clinical biplane cineangiography from 52 patients. The vessel axes in end diastole were reconstructed and geometric parameters that included measures of curvature, torsion and tortuosity were quantified for the proximal, middle and distal segments of the arteries.

Statistical analysis shows that (1) in the LAD, curvature, torsion and tortuosity are generally highest in the distal portion, (2) in the RCA, these parameters are smallest in the middle segment, (3) the LAD exhibits significant higher torsion than the RCA (*P*<0.005), and (4) >80% of the variability of coronary arterial geometry can be expressed in terms of two factors, one dominated by the curvature measures and tortuosity, and the other emphasizing the torsion parameters.

This study has comprehensively documented the normal arterial geometry of the LAD and RCA in end diastole. This information may be used to guide the identification of geometric features that might be atherogenic risk factors.

It has been well documented through extensive epidemiological studies that human coronary atherosclerosis shows a non-uniform distribution in the coronary vasculature [1–4]. The discrepancy in disease occurrence exists not only among different coronary arteries [5], but along different sites of the same vessel as well [6]. This phenomenon cannot be explained by conventionally established risk factors, such as hypercholesterolemia, hypertension, diabetes, smoking, and male gender, which are all systemic and presumably influence the entire coronary vasculature and hence provide no mechanisms for disease localization.

The preferential disease distribution has led to a widespread speculation that local mechanical forces may play an important role in the initiation and localization of atherosclerosis. In fact, this speculation has been supported by ample evidence gathered from laboratory investigations that include in vivo or vitro experiments, cell culture studies and gene expression profiling. For instance, experiments on vessel wall cells in culture have demonstrated that they respond morphologically, functionally, metabolically, and at the molecular level to a variety of hemodynamic stresses (eg. [7–10]); these stresses vary regionally, and some of the responses they elicit are likely to be involved in the atherosclerotic process. Similar results have also been reported from in vivo studies (eg. [11–13]). It has also been argued that mechanical stresses in the vessel wall might act as mediators of the atherosclerotic process. Such stresses, which also vary from site to site, have been claimed to correlate with the location of atherosclerotic disease [14]. Furthermore, vascular wall cells in culture and in vivo have been shown to respond to imposed strains by activating signaling pathways and remodeling (eg. [15–18]). Therefore, it seems plausible that the non-uniformity in disease distribution is due, at least in part, to the presence of certain mechanical forces that vary spatially; these forces might mediate disease development in concert with some of the established systemic risk factors [19].

While the notion that biomechanics may be involved in the atherosclerotic process is by now widely accepted, to date the precise mechanism of its involvement remains unknown. The biomechanical factors that have been proposed to play a role in atherogenesis and atherosclerotic development include low fluid dynamic shear stress and impaired mass transfer, concentrations of mechanical stress in the wall, long near-wall residence times, oscillatory or reversing shear stress, high spatial or temporal gradients in shear, and an imprecisely defined "disturbed flow" characterized by oscillatory and low mean stress, and possibly including true turbulence.

Recognizing that vessel geometry has a major influence on the mechanical environment of the vessel wall, Friedman *et al* proposed the concept of “geometric risk factors” for the development of atherosclerosis [20]. It is hypothesized that certain geometric features of the human vasculature may be risk factors that increase vessel’s susceptibility to atherosclerotic disease by mediating an adverse mechanical environment in its vicinity. Indeed, this hypothesis has been supported by many studies based on *in vivo* observations or post-mortem examination [21–25]. One of the appealing aspects of the hypothesis is that it does not rely on any particular biomechanical mechanism of atherogenesis, since all of them depend quantitatively, and in some cases qualitatively, on the geometry of the conduit. Because the candidate mechanisms depend differently on geometry, the identification of geometric features that predispose to disease can help elucidate the mechanisms by which biomechanics participates in the pathological development of atherosclerosis.

To further evaluate the geometric risk factor hypothesis, it is necessary to understand the normal geometry of vascular regions of interest, and the variability, both along the vessel length and among different individuals, of geometric features that might affect the mechanical stress experienced by the wall. We begin here with a description of the static parameters of the human left anterior descending (LAD) and right coronary arteries (RCA), selected because of their clinical significance. Our long-term goal is to use these data and dynamic data to be derived subsequently, in relation to epidemiological data regarding the distribution of disease, to identify atherogenic geometric features that can predict risk and suggest mechanisms of atherogenesis. This work extends our earlier research on the dynamics of a limited number of coronary artery segments to include a larger number of cases, well defined arterial segments, and more comprehensive statistical analysis for the end diastolic geometry of the two coronary vessels.

Patients > 18 years of age undergoing elective coronary angiography at the cardiac catheterization laboratories of Duke University Medical Center and Vanderbilt University Medical Center were approached for enrollment in this study if they did not have a history or presence of acute or complicated myocardial infarction, renal insufficiency, uncontrolled hypertension, uncontrolled diabetes mellitus, or known severe coronary artery disease. Informed consent was obtained from each participating patient prior to angiography. The images were reviewed after acquisition, and those arteries with evident stenosis or unsatisfactory image quality, were excluded from further evaluation. Sixty-seven angiographically normal adult human coronary arteries (32 LADs and 35 RCAs) from 52 patients were finally selected for characterization of coronary geometry. The demographics of these patients are summarized in Table 1. The recruitment and study protocols were approved by the Institutional Review Boards of both medical centers.

The imaging and analysis procedures are similar to those described previously [26,27] and are summarized as follows.

Left or right coronary cineangiograms were obtained from recruited subjects using clinical biplane systems. Both anterior-posterior and lateral projections were used to image the vessel of interest. The projection angles were adjusted so that the entire vessel course was well visualized and imaged with minimal overlap with other vessels. The pixel dimensions were approximately 0.3×0.3 mm^{2}, and the rate of acquisition was 15 or 30 frames per second.

To allow reconstruction of the three dimensional (3-D) vessel course, the imaging geometry was derived from images of a Plexiglas calibration cube that were obtained subsequent to patient imaging. Embedded within the cube were twelve radio-opaque beads whose 3-D coordinates were known in relationship to one another. The calibration cube was placed approximately in the former position of the heart and imaged with the same imaging geometry as used with the patient. Based on the known spatial relationship of the 3-D coordinates of the beads and their corresponding 2-D coordinates in the image planes, two projection matrices, one for each projection direction, were derived. These matrices allow the coordinates of a 3-D point to be computed from those in the image planes.

Three-dimensional reconstruction was performed using end diastolic images, with an epipolar line based technique (See [28] for technical details). Briefly, two paired cineangiographic frames, one from each projection and obtained at end diastole, were identified. Then, using the projection matrices obtained as described above and the epipolar line technique, the 3-D coordinates of several salient landmarks that could be well visualized in both projections were reconstructed. The initial reconstruction, including the landmarks and connecting "strings", was iteratively refined until a 3-D vessel axis that minimized a cost function was found. The refinement was based on an extension of the well established 2-D “snake” model [29] into 3-D space, with a smoothness constraint on the 3-D vessel axis and a requirement that the image intensity in both projections be simultaneously minimized [30]. Fig. 1 shows the reconstruction results for an LAD and an RCA.

The geometry of the reconstructed 3-D vessel axis was characterized on the basis of curvature (κ), torsion (τ) and tortuosity. These measures were selected because they may be related to atherosusceptibility [22,26], and they are invariant under translation and rotation.

Curvature and torsion are defined according to the following equations [31]:

$$\kappa ={\left(\frac{{\mathit{(}y\prime z\u2033-y\u2033z\prime \mathit{)}}^{\mathit{2}}+{\mathit{(}z\prime x\u2033-z\u2033x\prime \mathit{)}}^{\mathit{2}}+{\mathit{(}x\prime y\u2033-x\u2033y\prime \mathit{)}}^{\mathit{2}}}{{\mathit{(}{x\prime}^{\mathit{2}}+{y\prime}^{\mathit{2}}+{z\prime}^{\mathit{2}}\mathit{)}}^{\mathit{3}}}\right)}^{\mathit{1}/\mathit{2}},$$

(1)

$$\tau =\frac{abs(|\begin{array}{ccc}x\prime \hfill & y\prime \hfill & z\prime \hfill \\ x\u2033\hfill & y\u2033\hfill & z\u2033\hfill \\ {x}^{}{y}^{}{z}^{}\hfill \hfill \hfill & |\end{array}){\mathit{(}y\prime z\u2033-y\u2033z\prime \mathit{)}}^{\mathit{2}}+{\mathit{(}z\prime x\u2033-z\u2033x\prime \mathit{)}}^{2}+{\mathit{(}x\prime y\u2033-x\u2033y\prime \mathit{)}}^{2},}{}$$

(2)

where primes denote derivatives of the coordinates of the curve with respect to distance along the curve. The shape of a 3-D curve is uniquely determined by these two parameters. Curvature and torsion were calculated for each point on the vessel axis, using a cubic polynomial fitting method and derivatives that were computed algebraically [31]. Points close to the ends of the axis were excluded because there were too few adjacent points to permit reliable curve fitting.

The tortuosity between two points was defined as:

$$Tortuosity=1-\text{chord length}\phantom{\rule{thinmathspace}{0ex}}/\text{arc length}.$$

(3)

The tortuosity defined as above has a range between zero and unity. It is zero when the vessel is straight, and increases when the vessel is more tortuous. It is unity when the two endpoints of the curve are coincident.

As disease susceptibility varies from region to region, the geometry of each vessel axis was characterized on the basis of the three major segments commonly referred to at clinic and recommended by the American Heart Association [32]. The region of the LAD that we study here extends from its origin at the bifurcation of the left main coronary artery to the apex of the heart. This region is divided into the proximal, middle, and distal segments (LADp, LADm, and LADd) by two landmarks: the origins of the first septal perforator branch and the second diagonal branch. For the RCA, the study region extends from the RCA ostium to the posterior descending branch. It is also separated into three segments (RCAp, RCAm, and RCAd) by the ostia of the right ventricular and acute marginal branches. If there was more than one ventricular branch, the one closest to a point midway between the RCA ostium and the acute marginal branch was used. If there were no ventricular branches, the midpoint was used. For each major vessel segment, seven geometric parameters were calculated: the mean, maximum, and standard deviation (SD) of curvature (*MeC*, *MaC*, and *SDC*) and torsion (*MeT*, *MaT*, and *SDT*), and the tortuosity (*TTS*).

The LAD and RCA data were analyzed separately. The distributions among the cases of all seven geometric parameters were skewed. Therefore, the median (M) and interquartile range [IQR, the range between the first (Q_{1}) and third quartile (Q_{3})] were used to summarize the population data instead of the mean and standard deviation (SD). The variation of each parameter was measured by the coefficient of quartile variation (CQV), which is defined as (Q_{3}-Q_{1})/(Q_{3}+Q_{1}). Paired *t*-tests were used to compare each parameter among the three segments of the LAD and the three segments of the RCA. The difference of these parameters between the LAD and RCA was examined by Wilcoxon rank-sum tests. A difference was considered statistically significant if *P*<0.05.

Prior to multivariate analysis, the skewed distributions of the parameters were transformed to approximately normal distributions. A logarithmic transformation [to base 10, lg(1+*x*)] was applied to all parameters except tortuosity, for which a square root transformation was used. Accordingly, the transformed parameters were named *LMeC*, *LMaC*, *LSDC*, *LMeT*, *LMaT*, *LSDT*, and *RTTS* respectively. The correlation matrix **R** of the seven transformed geometric parameters indicated that they were not independent of one another. Therefore, exploratory factor analysis (FA) based on principal components analysis (PCA) was used to summarize the geometric parameters by a smaller number of independent factors; outliers identified on the basis of their Mahalanobis distance were excluded from the analysis. There were three steps for the FA [33]: first, PCA was applied to **R**, and principal components (PCs) with eigenvalues greater than unity were retained; second, these PCs were then orthogonally rotated by the *varimax* method to maximize high correlations and minimize low ones between each PC and the seven geometric parameters. The rotated PCs were the solution of the factors; last, to find out the factor scores for each case, the score coefficients which transformed raw scores to factor scores were calculated. The coefficients can be used (see Appendix) to calculate the score of a given case in the factor space and determine whether it is statistically consistent with the population of previously characterized cases.

All statistical analyses were performed using JMP 6.0 (SAS Institute Inc., Cary, NC).

The distributions of the seven geometric parameters (the mean, maximum, and SD of curvature and torsion, and tortuosity) of each of the three segments of the 32 LADs and 35 RCAs are summarized in Table 2, including the median, interquartile range, and coefficient of quartile variation. Pairwise comparisons of the geometric parameters among the three major segments of the LAD and among those of the RCA, using paired t-tests, are given in Table 3.

Summary of the Geometric Parameters: Median (M), Interquartile Range (Q_{1} to Q_{3}), and Coefficient of Quartile Variation (CQV)

Pairwise Comparisons of Geometric Parameters among the Three Major Segments of the LADs and RCAs^{*}

In the LADs, the variation of *MeC* among the three segments is small, with the median values ranging from 0.50 cm^{−1} for the middle segment to 0.57 cm^{−1} for the distal segment; only the difference between the middle and distal segments is statistically significant (*P*= 0.01). *MaC*, however, exhibits large variations among the three segments, with the median value increasing from 0.87 cm^{−1} for the proximal segment to 1.59 cm^{−1} for the distal segment; pairwise differences among the three segments are all statistically significant. The variation in *MaC* is accompanied by a similar trend in *SDC*, whose median value increases from 0.23 cm^{−1} for the proximal segment to 0.37 cm^{−1} for the distal segment; however, the increment in this parameter between the proximal and middle segments does not reflect a statistically significant differnce between the two populations. The finding that *MaC* and *SDC* tend to increase distally probably reflects the increasing sinuosity of the more distal portions of the artery.

In the RCAs, the median values of *MeC* of the proximal (0.66 cm^{−1}) and distal (0.62 cm^{−1}) segments are higher than that of the middle segment (0.58 cm^{−1}), and there is a significant difference between the proximal and middle segments (*P* = 0.04). The median values of *MaC* of the proximal (1.27 cm^{−1}) and distal segments (1.51 cm^{−1}) are also higher than that of the middle segment (1.09 cm^{−1}), although none of the pairwise differences is significant. Similarly, the median values of *SDC* of both the proximal (0.36 cm^{−1}) and distal (0.38 cm^{−1}) segments is higher than that of the middle segment (0.27 cm^{−1}), and the *P*-value for the difference between the middle and distal segments is significant (*P* = 0.04). Overall, the middle segment has the lowest curvature measures among the three segments.

The torsion measures share some similarities with the curvature measures in the LADs. The median values of *MeT* are comparable among the three segments, with values of 3.31, 3.36 and 3.82 cm^{−1} for the proximal, middle and distal segments, respectively. Similarly to *MaC*, the median of *MaT* increases distally, from 10.4 cm^{−1} for the proximal segment to 26.4 cm^{−1} for the distal segment, but the difference between the middle and distal segments is not significant. The median *SDT* of the distal segment is higher than that of the proximal and middle segments, although the differences between the corresponding populations are not statistically significant at the *P*=0.05 level.

In the RCAs, the torsion measures follow the pattern of the curvature measures, but with more significant differences; the distal segment exhibits the highest values. The median value of *MeT* of the middle segment (2.03 cm^{−1}) is lower than that of the proximal segment (2.46 cm^{−1}) and distal segment (2.91 cm^{−1}), and there is a marginally significant difference (*P* = 0.05) between the middle and distal segments. The median value of *MaT* of the distal segment (15.25 cm^{−1}) is considerably higher than that of the proximal (10.3 cm^{−1}) and middle segments (7.27 cm^{−1}), and this is reflected in the pairwise tests. The same is true for *SDT*, though the differences among the medians are not so large.

In the LADs, the median value of *TTS* varies from 0.05 for the middle segment to 0.15 for the distal segment; the pairwise *t*-tests show that *TTS* is significantly higher in the distal segment than in the proximal or middle segments. The variation of *TTS* is lowest in the distal segment, suggesting that the high tortuosity of this segment has a low variability among individuals.

In the RCAs, the median value of *TTS* of the proximal segment (0.16) is similar to that of the distal segment (0.17), and both are considerably higher than that of the middle segment (0.06). This is paralleled by the results of the pairwise tests, which show significant differences in both cases. This result, and the curvature and torsion measures summarized above, consistently indicate that the course of the RCA is least tortuous in the middle segment. The CQV of *TTS*, however, is higher in the middle segment (59.7%) than in the proximal (34.0%) or distal (50.1%) segments.

The occurrences of extreme (highest or lowest) values of the geometric parameters among the major segments of the LAD and RCA are summarized in Table 4, which gives for each major segment the number of cases in which that segment exhibited the highest maximum curvature (*MaC _{H}*), highest and lowest mean curvature (

Region-by-region analysis in the preceding sections demonstrates that the LAD generally has higher values of curvature, torsion and tortuosity at the distal portion, whereas the RCA consistently has lower values of these measures in the middle segment. To compare the overall geometries of the LAD and RCA, the seven geometric parameters were calculated for the entire length of the vessels. The comparison of these parameters is shown in Table 5. It can be seen that there are no significant differences in curvature measures between the LAD and the RCA; the LAD exhibits significantly higher torsion measures, but lower tortuosity, than the RCA. These findings are in keeping with our earlier report [34].

Seven geometric parameters were studied individually in the above analysis. The results show that there are quite large variations in LAD and RCA geometry among both individuals and vascular segments. Correlation analyses of these parameters revealed that they were not independent of one another. Therefore, FA was applied to derive a smaller number of independent factors that account for most of the variation of the seven parameters.

The transformed parameters *LMeC*, *LMaC*, *LSDC*, *LMeT*, *LMaT*, *LSDT*, and *RTTS*, which are summarized in Table 6, were used in the PCA, which strictly applies to normally distributed variables. Three segments in two LAD's and two segments in two RCA's were excluded from the analysis as outliers.

Table 7 gives the results of FA of the seven geometric parameters of the LAD and RCA segments. As shown in the table, for all six vessel segments, only two factors are needed to account for 82% to 89% of the total variance in each segment. It is interesting to notice that the two independent factors for each vessel segment of both the LAD and RCA are very similar, one (denoted as F_{C}) with dominant factor loadings from the three curvature measures (mean, maximum and SD) and tortuosity, and the other one (denoted as F_{T}) associated with the variation in the three torsion parameters. The score coefficients for calculating the factor scores of each case are listed in Table 8. In the Appendix, we will demonstrate how to use the results from the FA.

The geometry of the human vasculature and its implications for atherogenesis have been addressed in many studies [21–25,35–43]. Most previous characterizations of vascular geometry were focused on a number of geometric features in disease prone regions, and provided valuable insights into the role of vessel geometry in atherogenesis. Earlier studies of the entire coronary vasculature were focused on anatomical aspects such as the locations and dimensions of coronary artery segments [44–46]. The present study catalogues the axial geometry of 67 coronary arteries in a region-specific manner, defines a set of parameters that characterize the end diastolic axial geometry of the LAD and RCA, and describes their variations among individuals and vascular regions. This represents a first comprehensive documentation of this kind. It complements the former studies, and has the potential of adding new insights and providing useful tools for studying coronary disease by exploring the causal relationships among vascular geometry, biomechanics and atherogenesis. Further study will incorporate the coronary geometry at systolic phase so that dynamic nature of the heart can be reflected. The geometric features characterized in this study include curvature, torsion and tortuosity, with which the tortuous nature of the coronary artery can be described. This tortuous nature is, in fact, closely tied with unusual fluid dynamic and wall mechanical forces, particularly in the presence of the phasic cardiac motion that is imposed on these vessels. Presumably, the unusual biomechanics is among the primary factors that may put the vessel at a greater risk for atherosclerosis.

A key premise of the concept of geometric risk factors is that, for a particular geometric feature to mediate atherosclerotic risk, it must exhibit large variability among individuals or vascular regions, sufficient to cause a significant variability in mechanical stresses. In this study, characterization of 32 LADs and 35 RCAs indeed confirmed the presence of considerable geometric variability in the human coronary vasculature. For instance, in the LAD, there is nearly a two-fold difference between segments in maximum curvature and a 2.5-fold difference in maximum torsion. The CQV of mean curvature is 25% or more in all segments, and that of maximum torsion is nearly 60% in the middle segment, demonstrating that some geometric variables vary considerably among individuals. A geometric parameter that exhibits a large variability in the population may induce a corresponding variability in the mechanical environment of the vessel, and therefore would be more likely to be responsible for individual differences in disease susceptibility and localization.

Factor analysis reveals that, in every segment examined, much of the geometric variation of the coronary arteries is attributable to two factors, to which curvature-tortuosity and torsion measures are the dominant contributors. As can be seen in Table 7, the two factors explain more than 80% of the total variance in the transformed geometric parameters. The contributions of the two are similar; differing by no more than 13% in the LAD distal segment. Therefore, these two factors can form a basis for characterizing normal coronary artery geometry. The results in Table 6 and Table 8 provide a good starting point for this kind of study. The table entries can be updated as more cases are examined. Furthermore, when a certain number of cases has been analyzed, the effects of gender, age, and race *et al*. on coronary geometry can be examined. When norms have been established, it should be possible to identify individuals whose geometry is abnormal and whose risk of atherosclerotic disease may be more or less than expected on the basis of traditional risk factors.

Comparisons of geometric parameters between the LAD and RCA show that, overall, the curvature measures in the two vessels are comparable, although their distributions along the vessel are different. The torsion measures, however, are more than 30% higher in the LAD than in the RCA. Further investigation is needed to determine whether the difference in torsion measures is responsible for the greater disease susceptibility of the LAD [47], and if so, how torsion interacts with hemodynamic and wall mechanical parameters. For example, torsion may play a more important role than curvature in generating helical flow patterns in the vessel; helical flow has been found to relate to coronary disease [48].

The goal of this study is to catalogue the normal values and variabilities of coronary artery geometric parameters to guide future exploration and identification of geometric risk factors. The identification of vessel geometric risk factors will not only illuminate the mechanical mechanisms of atherogenesis, but potentially can be valuable to clinical practice as well. For instance, measurements of coronary geometry may be used prognostically in a patient-specific manner, as noted above. Knowledge of adverse geometric features can support surgical planning to achieve a satisfactory post-surgical geometry. Technically, measuring the geometry of coronary arteries is quite convenient at clinic, given the availability of biplane cineangiography and the increasing applicability of less invasive procedures such as x-ray computed tomography and magnetic resonance angiography.

It is worth noting that the axial geometry of the coronary arteries, the main focus of this study, is not the sole determinant of the mechanical environment of the vessel. Other geometric parameters, most notably those describing the vessel cross-sectional geometry, such as the vessel diameter and wall thickness, are also important determinants. A combination of axial and cross-sectional geometric parameters may be most useful for studying geometric risk factors. The inclusion of wall thickness would require another catheter based intervention, presumably intravascular ultrasound; however, lumen diameter data can be obtained from the angiograms. A mechanically appropriate way to introduce lumen diameter is to use it to multiplicatively nondimensionalize curvature and torsion. The nondimensional product of tube radius and axial curvature (or torsion) is the natural variable for parameterizing flow in a bend (or helix), and these variables will also be included in future work.

Finally, it should be pointed out that, in this study, the geometry of the coronary arteries was characterized at the end diastolic phase of the cardiac cycle, similar to that in [45,46]. Indeed, with the techniques used here, it is possible to characterize the coronary arterial geometry throughout the cardiac cycle. Measurements at specific times during the cycle, such as end diastole and end systole, will allow the dynamics of coronary geometry to be characterized. This effort is under way, with the hope that it will add a dynamic dimension to the definition of geometric risk factors.

Grant support: NIH HL058856

This study is supported by the NIH grant HL058856.

We can use the data from the FA to evaluate the coronary geometry of a new case. For example, after 3-D reconstruction of the axis, the proximal segment of an LAD had the original geometric measures *MeC*, *MaC*, *SDC*, *MeT*, *MaT*, *SDT*, and *TTS* as 0.41cm^{−1}, 0.76 cm^{−1}, 0.24 cm^{−1}, 4.02 cm^{−1}, 23.52 cm^{−1}, 7.21 cm^{−1}, and 0.046. After the logarithmic transformation for the first six parameters and square root transformation of the last one, the mean and SD values of LAD proximal segment in Table 6 were used to standardize the data, resulted in −0.521, −0.221, 0.201, 0.290, 1.186, 0.974, and −0.287. By the product of these values and the factor score coefficients in Table 8, the scores on F_{C} and F_{T} of this case were obtained as −0.16 and 0.83, which means they are within about 0.2 and 0.8 standard deviation of the two factors’ population distributions respectively. So the geometry of this case is in a very normal range of population.

**Publisher's Disclaimer: **This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

1. Monteneg MR, Eggen DA. Topography of atherosclerosis in coronary arteries. Lab Invest. 1968;18(5):586–593. [PubMed]

2. Debakey ME, Lawrie GM, Glaeser DH. Patterns of atherosclerosis and their surgical significance. Ann Surg. 1985;201(2):115–131. [PubMed]

3. Cornhill JF, Herderick EE, Stary HC. Topography of human aortic sudanophilic lesions. Monogr Atheroscler. 1990;15:13–19. [PubMed]

4. Schwartz CJ, Chandler AB. Coronary atherosclerotic disease: Pathological background. In: Fuster V, editor. Syndromes of atherosclerosis: Correlations of clinical imaging and pathology. Armonk, NY: Futura; 1996.

5. Schmermund A, Baumgart D, Mohlenkamp S, Kriener P, Pump H, Gronemeyer D, Seibel R, Erbel R. Natural history and topographic pattern of progression of coronary calcification in symptomatic patients - an electron-beam ct study. Arterioscler Thromb Vasc Biol. 2001;21(3):421–426. [PubMed]

6. Halon DA, Sapoznikov D, Lewis BS, Gotsman MS. Localization of lesions in the coronary circulation. Am J Cardiol. 1983;52(8):921–926. [PubMed]

7. Resnick N, Collins T, Atkinson W, Bonthron DT, Dewey CF, Jr, Gimbrone MA., Jr Platelet-derived growth factor b chain promoter contains a cis-acting fluid shear-stress-responsive element. Proc Natl Acad Sci USA. 1993;90(10):4591–4595. [PubMed]

8. Davies PF. The convergence of haemodynamics, genomics, and endothelial structure in studies of the focal origin of atherosclerosis. Biorheology. 2002;39(3):299–306. [PMC free article] [PubMed]

9. Chen BP, Li YS, Zhao Y, Chen KD, Li S, Lao J, Yuan S, Shyy JY, Chien S. DNA microarray analysis of gene expression in endothelial cells in response to 24-h shear stress. DNA microarray analysis of gene expression in endothelial cells in response to 24-h shear stress. 2001;7(1):55–63. [PubMed]

10. McCormick SM, Eskin SG, McIntire LV, Teng CL, Lu CM, Russell CG, Chittur KK. DNA microarray reveals changes in gene expression of shear stressed human umbilical vein endothelial cells. Proc Natl Acad Sci USA. 2001;98(16):8955. [PubMed]

11. Berceli SA, Warty VS, Sheppeck RA, Mandarino WA, Tanksale SK, Borovetz HS. Hemodynamics and low density lipoprotein metabolism: Rates of low density lipoprotein incorporation and degradation along medial and lateral walls of the rabbit aorto-iliac bifurcation. Arterioscler Thromb Vasc Biol. 1990;10(5):686–694. [PubMed]

12. Gibson CM, Diaz L, Kandarpa K, Sacks FM, Pasternak RC, Sandor T, Feldman C, Stone PH. Relation of vessel wall shear stress to atherosclerosis progression in human coronary arteries. Arterioscler Thromb Vasc Biol. 1993;13(2):310–315. [PubMed]

13. Himburg HA, Grzybowski DM, Hazel AL, LaMack JA, Li XM, Friedman MH. Spatial comparison between wall shear stress measures and porcine arterial endothelial permeability. Am J Physiol Heart Circ Physiol. 2004;286(5):1916–1922. [PubMed]

14. Thubrikar MJ, Robicsek F. Pressure-induced arterial wall stress and atherosclerosis. The Annals of Thoracic Surgery: Soc Thorac Surgeons. 1995:1594–1603. [PubMed]

15. Iba T, Mills I, Sumpio BE. Intracellular cyclic amp levels in endothelial cells subjected to cyclic strain in vitro. Intracellular cyclic AMP levels in endothelial cells subjected to cyclic strain in vitro. 1992;52(6):625–630. [PubMed]

16. Zhao S, Suciu A, Ziegler T, Moore JE, Burki E, Meister JJ, Brunner HR. Synergistic effects of fluid shear stress and cyclic circumferential stretch on vascular endothelial cell morphology and cytoskeleton. Arterioscler Thromb Vasc Biol. 1995;15(10):1781–1786. [PubMed]

17. Lee RT, Yamamoto C, Feng Y, Potter-Perigo S, Briggs WH, Landschulz KT, Turi TG, Thompson JF, Libby P, Wight TN. Mechanical strain induces specific changes in the synthesis and organization of proteoglycans by vascular smooth muscle cells. J Biol Chem. 2001;276(17):13847–13851. [PubMed]

18. Li W, Chen Q, Mills I, Sumpio BE. Involvement of s 6 kinase and p 38 mitogen activated protein kinase pathways in strain-induced alignment and proliferation of bovine aortic smooth muscle cells. J Cell Physiol. 2003;195(2):202–209. [PubMed]

19. VanderLaan PA, Reardon CA, Getz GS. Site specificity of atherosclerosis: Site-selective responses to atherosclerotic modulators. Arterioscler Thromb Vasc Biol. 2004;24(1):12–22. [PubMed]

20. Friedman MH, Deters OJ, Mark FF, Bargeron CB, Hutchins GM. Arterial geometry affects hemodynamics - a potential risk factor for atherosclerosis. Atherosclerosis. 1983;46(2):225–231. [PubMed]

21. Fisher M, Fieman S. Geometric factors of the bifurcation in carotid atherogenesis. Stroke. 1990;21(2):267–271. [PubMed]

22. Smedby O. Geometrical risk factors for atherosclerosis in the femoral artery: A longitudinal angiographic study. Ann Biomed Eng. 1998;26(3):391–397. [PubMed]

23. Friedman MH, Baker PB, Ding ZH, Kuban BD. Relationship between the geometry and quantitative morphology of the left anterior descending coronary artery. Atherosclerosis. 1996;125(2):183–192. [PubMed]

24. Ding ZH, Biggs T, Seed WA, Friedman MH. Influence of the geometry of the left main coronary artery bifurcation on the distribution of sudanophilia in the daughter vessels. Arterioscler Thromb Vasc Biol. 1997;17(7):1356–1360. [PubMed]

25. Wood NB, Zhao SZ, Zambanini A, Jackson M, Gedroyc W, Thom SA, Hughes AD, Xu XY. Curvature and tortuosity of the superficial femoral artery: A possible risk factor for peripheral arterial disease. J Appl Physiol. 2006;101(5):1412–1418. [PubMed]

26. Ding ZH, Zhu H, Friedman MH. Coronary artery dynamics in vivo. Ann Biomed Eng. 2002;30(4):419–429. [PubMed]

27. Ding ZH, Friedman MH. Quantification of 3-d coronary arterial motion using clinical biplane cineangiograms. Int J Card Imaging. 2000;16(5):331–346. [PubMed]

28. Guggenheim N, Doriot PA, Dorsaz PA, Descouts P, Rutishauser W. Spatial reconstruction of coronary-arteries from angiographic images. Phys Med Biol. 1991;36(1):99–110. [PubMed]

29. Kass M, Witkin A, Terzopoulos D. Snakes - active contour models. Snakes - Active Contour Models. 1987;1(4):321–331.

30. Zhu H, Friedman MH. Tracking 3-d coronary artery motion with biplane angiography. Proc 2002 IEEE Int Symp Biomed Imaging: From Nano to Macro. 2002:605–608.

31. Pao YC, Lu JT, Ritman EL. Bending and twisting of an invivo coronary-artery at a bifurcation. J Biomech. 1992;25(3):287–295. [PubMed]

32. Austen WG, Edwards JE, Frye RL, Gensini GG, Gott VL, Griffith LS, McGoon DC, Murphy ML, Roe BB. A reporting system on patients evaluated for coronary artery disease. Report of the ad hoc committee for grading of coronary artery disease, council on cardiovascular surgery, american heart association. Circulation. 1975;51(4 Suppl):5. [PubMed]

33. Tabachnick BG, Fidell LS. Boston: Allyn and Bacon; 2001. Using multivariate statistics.

34. Ding ZH, Friedman MH. Dynamics of human coronary arterial motion and its potential role in coronary atherogenesis. J Biomech Eng. 2000;122(5):488–492. [PubMed]

35. Demirbag R, Yilmaz R. Effects of the shape of coronary arteries on the presence, extent, and severity of their disease. Heart Vessels. 2005;20(5):224–229. [PubMed]

36. Ikeda U, Kuroki M, Ejiri T, Hosoda S, Yaginuma T. Stenotic lesions and the bifurcation angle of coronary-arteries in the young. Jpn Heart J. 1991;32(5):627–633. [PubMed]

37. Sakata N, Takebayashi S. Localization of atherosclerotic lesions in the curving sites of human internal carotid arteries. Biorheology. 1988;25(3):567–578. [PubMed]

38. Endoh R, Homma T, Furihata Y, Sasaki Y, Fukushima T. A morphometric study of the distribution of early coronary atherosclerosis using arteriography. Artery. 1988;15(4):192–202. [PubMed]

39. Smedby O, Hogman N, Nilsson S, Erikson U, Olsson AG, Walldius G. 2-dimensional tortuosity of the superficial femoral-artery in early atherosclerosis. J Vasc Res. 1993;30(4):181–191. [PubMed]

40. Friedman MH, Ding ZH. Relation between the structural asymmetry of coronary branch vessels and the angle at their origin. J Biomech. 1998;31(3):273–278. [PubMed]

41. Friedman MH, Ding ZH. Variability of the planarity of the human aortic bifurcation. Med Eng Phys. 1998;20(6):469–472. [PubMed]

42. Dvir D, Kornowski R, Ben-Gal T, Berman M, Vidne B, Aravot D. Relation of amounts of narrowing to the length of the right coronary artery. Am J Cardiol. 2002;90(1):46–48. [PubMed]

43. Dvir D, Kornowski R, Gurevich J, Orlov B, Aravot D. Degrees of severe stenoses in sigma-shaped versus c-shaped right coronary arteries. Am J Cardiol. 2003;92(3):294–298. [PubMed]

44. Zamir M, Chee H. Segment analysis of human coronary-arteries. Blood Vessels. 1987;24(1–2):76–84. [PubMed]

45. Dodge JT, Brown BG, Bolson EL, Dodge HT. Intrathoracic spatial location of specified coronary segments on the normal human-heart - applications in quantitative arteriography, assessment of regional risk and contraction, and anatomic display. Circulation. 1988;78(5):1167–1180. [PubMed]

46. Dodge JT, Brown BG, Bolson EL, Dodge HT. Lumen diameter of normal human coronary-arteries - influence of age, sex, anatomic variation, and left-ventricular hypertrophy or dilation. Circulation. 1992;86(1):232–246. [PubMed]

47. Chatzizisis YS, Giannoglou GD, Parcharidis GE, Louridas GE. Is left coronary system more susceptible to atherosclerosis than right? A pathophysiological insight. Int J Cardiol. 2007;116(1):7–13. [PubMed]

48. Van Langenhove G, Wentzel JJ, Krams R, Slager CJ, Hamburger JN, Serruys PW. Helical velocity patterns in a human coronary artery: A three-dimensional computational fluid dynamic reconstruction showing the relation with local wall thickness. Circulation. 2000;102(3):E22. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |