Computational and Mathematical Methods in Medicine

Volume 2012, Article ID 135256, 7 pages

http://dx.doi.org/10.1155/2012/135256

## Relationship between Concentration Difference of Different Density Lipoproteins and Shear Stress in Atherosclerosis

^{1}Department of Thoracic and Cardiovascular Surgery, West China Hospital, Sichuan University, Chengdu 610041, China^{2}Institute of Biomedical Engineering, West China Medical Center, Sichuan University, Chengdu 610041, China^{3}State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China

Received 15 August 2011; Revised 16 October 2011; Accepted 8 November 2011

Academic Editor: Kyehan Rhee

Copyright © 2012 Wei Meng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Previous research has observed concentration polarization in LDL and HDL in the arterial system. However, there is no report that links this concentration polarization to the development of vascular atherosclerosis (AS). Therefore, the purpose of this study is to establish the relationship between concentration difference of LDL and HDL and shear stress using a carotid bifurcation vascular model. PTFE was employed to create the carotid bifurcation model. Endothelial cells were coated on the inner wall of the graft. In a recirculation system, HDL and LDL concentration were measured under two different ICA flow velocities at 5 different locations within our model. We report the following: (1) LDL and HDL concentration difference was observed in both high flow and low flow environments; (2) the degree of LDL and HDL concentration polarization varied depending of high flow and low flow environment; (3) absolute values of concentration difference between LDL and HDL at the inner wall surface decreased with the increase in shear stress when shear stress was more than 1.5 Pa. This variation trend would be more pronounced if shear stress were less than 0.5 Pa. Our study suggests that under the action of shear stress, concentration differences of LDL or HDL create a disturbance in the balance of atherogenic factors and anti-As factors, resulting in the occurrence of AS.

#### 1. Introduction

Epidemiology has demonstrated that concentration increment of low-density lipoprotein (LDL) and concentration decrement of high-density lipoprotein (HDL) in blood are independent risk factors for the development of atherosclerosis (AS) [1, 2]. They assume opposite roles: the former is a harmful lipoprotein while the latter is protective [1–4]. In the majority of cases, AS is found to be located at the bifurcations, bends, or stenosis sites of the artery. In these areas, blood flow is disturbed resulting in flow separation and vortex zone, which is termed localization of AS [5, 6]. Studies have revealed that reasons resulting in localization of AS include the following: (1) shear stress appearing at the surface of vascular wall and (2) areas of disturbed blood flow leads to its slowing and thus, gives LDL considerable time to interact with the inner wall. As a result, an increasing number of LDL can penetrate the vascular wall.

At present, research has observed concentration polarization in LDL and HDL in the arterial system; these studies also revealed that polarization degree is related to shear stress. Lipid concentration polarization plays an important role in AS development [7–10]. In our research on HDL concentration polarization, we found that under high shear stress circumstance, HDL concentration on the inner surface of vascular wall is inversely proportional to shear stress [10]. This result appeared to contradict with the accepted protective effect of HDL. Based on this observation, the following have been assumed: AS is a result of concerted action of multiple factors and a conclusion cannot be drawn from HDL concentration of the inner surface of vascular wall alone. Under high shear stress circumstance, HDL concentration absolute value will decrease; the protective effect for preventing AS development appears to decrease. However, when integrated consideration with LDL concentration is made, can it show the interactive function of concentration polarization of different density lipoproteins in AS development?

On the above-mentioned aspects, we studied the relationship between concentration difference of LDL and HDL on the inner surface of vascular wall and shear stress based on a carotid bifurcation vascular model which designed to mimic the configuration of the human carotid artery [10].

#### 2. Material and Methods

##### 2.1. Experimental Model

Semipermeable poly tetra fluoro ethylene (PTFE) was employed to create the carotid bifurcation model. The inner diameter of each region of a human carotid artery was measured by Doppler ultrasound and computerized tomography angiography (CTA) (Figures 2(c) and 2(d)). The measurements were amplified at a rate of 1 : 1.5 (Table 1) to create the bifurcation model of carotid artery that was subsequently coated with endothelial cells on the inner wall of the model to induce the endothelialization of the inner wall of the model as shown in Figure 1 [10, 11].

As shown in Figure 2(a), the model included the common carotid artery (CCA), external carotid artery (ECA), internal carotid artery (ICA), and internal carotid artery sinus (ICAS). The low shear stress core region and its margins were marked on the model. The low shear stress core region was marked as point 5, and the proximal and distal margins were marked as points 3 and 4 (Figure 2(b)). The measuring points of the inner diameter of the CCA and the inner diameter of the ICA were used as control points and marked as points 1 and 2 (Figure 2(b)). In addition, the final model is provided in the schematic diagram in Figure 3(a), and the schematic diagram of the experimental system was shown in Figure 3(b). The whole recirculate system was composed of the upstream reservoir, the model, the flow meter, a centrifugal pump, and the downstream reservoir, connected with tubes. The system is connected with tubing form a complete circuit. A centrifugal pump would pump the fluid from the downstream back to the upper reservoir to establish the circuit. (The fluid was conserved by a reservoir with adjustable height. Before entering the model, the fluid would pass through a segment of straight and horizontal tube. After leaving the model, the fluid flowed forwards to a downstream reservoir through tubes with flow meter. A pump would suck the fluid of the downstream reservoir back to the upper one in order to establish a circulation.) A threshold bar in the upper reservoir could automatically introduce extra fluid back into the downstream container to maintain a stable experimental water stage. The flow direction was shown by arrows (Figure 3(b)).

##### 2.2. Numerical Simulation Methodology

To simplify the analysis, the following assumptions were made: (1) the circulating liquid is incompressible Newtonian fluid; (2) the flow is a steady flow; (3) the wall of bifurcation model is permeable to plasma and has a filtration rate of the order of 10^{6} cm/s. If the volume power, heat exchange, and other physical and chemical factors are not considered then the equations provided in the footnote can be employed. The first formula is a continuity equation and the second formula is an equation of motion where is the velocity of flow field, is the fluid pressure, *ρ* is the fluid density and * μ* is the fluid viscosity. Because the basic equations mentioned above are intensive non-linear equations, the Finite Volume Method (FVM) was employed as it is the most commonly used numerical method for resolving this type of mathematical problem at present [12]. The FLUENT software is hydrodynamic calculating software based on the FVM and it is the CFD software applied most widely utilized.

*The boundary conditions are*

The bifurcation model was considered as a semipermeable membrane with a filtration rate of Vw.

##### 2.3. Hydrodynamic Parameters of Blood Flow

The parameters of the blood flow in the model were controlled by the particle image velocimetry (PIV) using a type PIV-400-10 (TSI Company, Shoreview, MN, USA) and numerical simulation (NS) [13]. The average flow velocity of ICA was set at 0.559 m/s, which is the average flow velocity within the ICA measured in the human body at 150 mmHg blood pressure.

A solution prepared with 7.5% glycerol with a viscosity of 0.782 mPa.s and a density of 1.005 × 10^{3} kg/m^{3} (for pre run and stabilize the equipment) was measured by the Low Shear 30 (CONTRAVES LOW SHEAR 30 ISCOMETER, Swiss). This viscosity was chosen because it was the same viscosity as the M199 culture medium used for the endothelial cells. Blood flow parameters determined were flow rate (mL/s) and velocity of the circulation liquid through the model (flow velocity, m/s).

##### 2.4. Separation of HDL and LDL

Human plasma lipoproteins were collected and separated using the one-time density gradient ultracentrifugation method described by Zhang and Liu [14]. HDL and LDL bands were collected from centrifuged samples using a long syringe needle and were dialyzed in a buffer containing 0.02 mol/L Tris-HCl, 0.85% NaCl, 0.01% EDTA, and 0.01% NaN_{3} at pH 7.6. Dialysis was performed at 4°C in the dark for 6 hours each time and completely repeated 4 times in order to remove sodium bromide. Collected lipoproteins were stored at 4°C following filtration (storage and dialysis were performed under a nitrogen atmosphere to avoid oxidation). Figure 4 shows the high purity of the isolated LDLs and HDLs. In total, 100 mL of circulation liquid was prepared for the experiments, including 80 ml M199 medium and 20 mL of separated human plasma lipoproteins (i.e., 10 mL LDL and 10 mL HDL). Lipoprotein concentrations were determined with an OLYMPUS automatic biochemical analyzer (OLYMPUS automatic biochemical analyzer AU2700, Japan). The concentrations of LDL and HDL in the bulk flow () were 0.575 mmol/L and 0.242 mmol/L, respectively.

##### 2.5. Experimental Procedure

HDL and LDL concentration were measured under two different ICA flow velocities at 5 different locations in the model. The low-speed group had an average ICA flow velocity of 0.559 m/s (the average flow velocity of ICA in human body under 150 mmHg blood pressure), while the high-speed group had an average ICA flow velocity of 1.451 m/s (the peak flow velocity of ICA in human body under 90 mmHg blood pressure). Hydrodynamic parameters (i.e., blood flow in the low- and high-speed groups) were measured and have been summarized in Table 2.

After allowing the model to stabilize for 30 minutes, 50 *μ*L samples were sequentially collected from each of the 5 locations and 5 samples were collected consecutively from each location. The samples were collected 15 minutes apart to ensure that samples were collected at a constant flow. The collected samples were individually placed in polyethylene tubes and stored protected from light in brown bottles at 4°C. Lipoprotein concentrations were measured within 4 hours of collection. The ratio of the concentration of LDL at the surface () to the concentration in bulk was used as an index for concentration polarization of LDL. The ratio of the concentration of HDL at the surface () to the concentration in bulk was used as an index for concentration polarization of HDL. Polarization of LDL and HDL was considered to have occurred if the ratio was greater than 1.000.

##### 2.6. Statistical Analysis

The experimental result is shown in mean number ± standard deviation . Comparison between multiple groups was performed using one-factor analysis of variance, whereas comparison between two groups involved the *t*-test. Data correlation contrast analysis was performed using linear regression. The result was deemed as a significant difference when *P* was below 0.05. All statistical analyses were computed using the SPSS11.5 statistics package.

#### 3. Results

##### 3.1. The Surface Concentration of LDL and HDL at Low-Speed Flow Group

See Table 3.

##### 3.2. The Surface Concentration of LDL and HDL at High-Speed Flow Group

From Tables 3 and 4, LDL and HDL concentration polarization on the inner surface of carotid artery bifurcation model was observed at every sampling point in both two groups.

##### 3.3. Relationship between Concentration Difference of HDL and LDL and Shear Stress

LDL and HDL concentration difference on the inner surface of carotid artery bifurcation model was observed.

Absolute values of concentration difference between LDL and HDL on each sampling point of low- and high-speed flow groups (Table 5) were compared. No statistic difference was observed when absolute values of concentration difference between LDL and HDL on Sampling Points 1 and 2 were compared. However, statistical difference was noted by comparing absolute values of concentration difference between LDL and HDL on Sampling Points 3, 4, and 5; value was less than 0.05.

Comparing the ratio of LDL concentration at the inner wall surface on each sampling point of low- and high-speed flow group to concentration in bulk flow and the ratio of HDL concentration at the inner wall surface to concentration in bulk flow, the difference value showed the following: ratio of LDL concentration at the inner wall surface in low-speed flow to concentration in bulk flow was significantly higher than HDL; the difference between these two values was a positive value. However, the ratio of LDL concentration at the inner wall surface in high-speed flow to concentration in bulk flow was significantly lower than HDL; the difference between these two values was a negative value. Statistical difference was observed by comparing difference between the ratios of LDL and HDL concentration at the inner wall surface from these two groups to concentration in bulk flow.

As demonstrated in Figure 5, absolute values of concentration difference between LDL and HDL at the inner wall surface decreased with the increase in shear stress when shear stress was less than 1.5 Pa. This variation trend would be more pronounced if shear stress were less than 0.5 Pa. After shear stress increased to 1.5 Pa, absolute values of concentration difference between LDL and HDL at the inner surface of wall remained unchanged according to shear stress. Statistics revealed that on the inner wall surface of ICAS, absolute values of concentration difference between LDL and HDL at the inner wall surface were negatively correlated with shear stress (, ). After shear stress increased to 1.5 Pa, concentration difference between LDL and HDL at the inner wall surface was maintained at a relatively constant lower level. No negative correlation was observed.

Results revealed a concentration difference between LDL and HDL at the inner wall surface, particularly on the inner wall surface of ICAS.

#### 4. Discussion

Concentration increment of LDL and concentration decrement of HDL are independent risk factors for the development of AS. They assume opposite roles. While the former is a harmful lipoprotein, the latter is a protective one. The major reasons for the anti-AS mechanism of HDL are as follows: participating reverse cholesterol transcription, anti-LDL oxidation, and the protection of endothelial cells [15–22]. In the atherogenic mechanism of LDL, there are more ROS (, H_{2}O_{2} and ect) in the area with low shear stress because of increased expression of NADPH oxidase and defect in oxygen transfer. This accelerates LDL oxidation and generates Ox-LDL. Ox-LDL can induce expression of vascular endothelial cell MCP-1 to result in adhesion and migration of monocyte. Similarly, it can reduce expressions of antiapoptotic proteins, Bcl-2, and c-IAP-1, by combining specific receptors LOX-1 to induce cell apoptosis. Thus, concentration polarization of LDL provides pathogenic Ox-LDL with sufficient substrate, and in the setting of low shear stress environment, makes the endothelial cells in these areas more vulnerable to the damage, and thus facilitates occurrence of atherosclerosis [23–26].

Concentration polarization degree of HDL and LDL is correlated with shear stress. However, the degree of HDL polarization and LDL polarization is different under different shear stress. As evidenced in Table 5, under high-speed flow circumstance, shear stress at the inner wall surface of ICAS (sampling point 5) increased when low-speed flow was compared, and the LDL concentration polarization degree was lower than HDL. This indicated that concentration polarization degree of lipoprotein with different molecular weight and size is different under different shear stress. We have furthered the study on this phenomenon and proposed the concept of “concentration difference” between HDL and LDL.

According to the phenomenon observed in experiments, we compared the absolute values of LDL-HDL concentration difference at the wall surfaces of every point in low- and high-speed flow with difference value of ratio of concentration in bulk flow. Statistical results revealed that in CCA and ICA (sampling points 1 and 2), no obvious difference in absolute values of concentration difference was observed. However, in ICAS, absolute value of concentration difference between LDL and HDL at the inner wall surface in low-speed flow was significantly higher than that in high-speed flow; in the ICAS area with wall shear stress less than 0.5 Pa, the absolute value of concentration difference between LDL and HDL at the inner wall surface was negatively correlated with shear stress. In ICAS, a significant difference was observed between the difference value of ratios of LDL and HDL concentration at the inner wall surface in high- and low-speed flow to concentration in bulk flow. In low-speed flow, wall shear stress in ICAS part was lower, the difference value between the ratio of concentration of LDL at the inner wall surface to concentration in bulk flow and ratio of concentration of HDL to concentration in bulk flow was a positive value. In high-speed flow, wall shear stress in ICAS was obviously increased, the difference value between the ratio of concentration of LDL at the inner wall surface to concentration in bulk flow and ratio of concentration of HDL to concentration in bulk flow was a negative value. This indicated that concentration polarization degree of LDL was higher than HDL in low-speed flow and low shear stress, and concentration polarization degree of HDL was higher than LDL in high-speed flow and high shear stress.

We first proposed that in the situation with certain low shear stress in the flow field, concentration polarization at the inner wall surface of lipoproteins with different molecular weights would change along with change in shear stress, resulting in inconsistent concentration profile; this is referred to as concentration difference. According to our study on concentration difference, in low stress situations, the increase in harmful LDL outweighs the increase in protective HDL. Other studies have suggested that the occurrence probability of AS angiocardiopathy would increase by 1%-2% if LDL level increased by 1% or if HDL level reduced by 1% [27, 28].

#### 5. Conclusion

In conclusion, we propose that the formation of AS may be secondary to an imbalance between the proatherosclerotic LDL and the antiatherosclerotic HDL. We demonstrated that in low stress environment, the increase in LDL outweighs the increase in HDL leading to AS. On the other hand, in high stress environment, the increases in HDL outweighs the increase in LDL which leads to a protective effect on atherosclerosis.

#### Authors’ Contribution

Y. Guo and Y. Shi contributed equally to this paper.

#### Acknowledgment

This Work Was Supported by the National Natural Science Foundation of China (81170288).

#### References

- A. M. Gotto Jr. and E. A. Brinton, “Assessing low levels of high-density lipoprotein cholesterol as a risk factor in coronary heart disease: a working group report and update,”
*Journal of the American College of Cardiology*, vol. 43, no. 5, pp. 717–724, 2004. View at Publisher · View at Google Scholar · View at Scopus - A. von Eckardstein and G. Assmann, “Prevention of coronary heart disease by raising high-density lipoprotein cholesterol?”
*Current Opinion in Lipidology*, vol. 11, no. 6, pp. 627–637, 2000. View at Publisher · View at Google Scholar · View at Scopus - J. Chen, J. L. Mehta, N. Haider, X. Zhang, J. Narula, and D. Li, “Role of caspases in Ox-LDL-induced apoptotic cascade in human coronary artery endothelial cells,”
*Circulation Research*, vol. 94, no. 3, pp. 370–376, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. R. Nofer, M. Walter, and G. Assmann, “Current understanding of the role of high-density lipoproteins in atherosclerosis and senescence,”
*Expert Review of Cardiovascular Therapy*, vol. 3, no. 6, pp. 1071–1086, 2005. View at Publisher · View at Google Scholar · View at Scopus - C. G. Caro, J. M. Fitz-Gerald, and R. C. Schroter, “Arterial wall shear and distribution of early atheroma in man,”
*Nature*, vol. 223, no. 5211, pp. 1159–1161, 1969. View at Publisher · View at Google Scholar · View at Scopus - C. G. Caro, J. M. Fitz-Gerald, and R. C. Schroter, “Atheroma and arterial wall shear. Observation, correlation and proposal of a shear dependent mass transfer mechanism for atherogenesis,”
*Proceedings of the Royal Society of London B*, vol. 177, no. 46, pp. 109–159, 1971. View at Google Scholar · View at Scopus - S. Wada and T. Karino, “Theoretical study on flow-dependent concentration polarization of low density lipoproteins at the luminal surface of a straight artery,”
*Biorheology*, vol. 36, no. 3, pp. 207–223, 1999. View at Google Scholar · View at Scopus - S. Wada and T. Karino, “Theoretical prediction of low-density lipoproteins concentration at the luminal surface of an artery with a multiple bend,”
*Annals of Biomedical Engineering*, vol. 30, no. 6, pp. 778–791, 2002. View at Publisher · View at Google Scholar · View at Scopus - G. Wang, X. Deng, and R. Guidoin, “Concentration polarization of macromolecules in canine carotid arteries and its implication for the localization of atherogenesis,”
*Journal of Biomechanics*, vol. 36, no. 1, pp. 45–51, 2003. View at Publisher · View at Google Scholar · View at Scopus - W. Meng, F. Yu, H. Chen et al., “Concentration polarization of high-density lipoprotein and its relation with shear stress in an in vitro model,”
*Journal of Biomedicine and Biotechnology*, vol. 2009, Article ID 695838, 8 pages, 2009. View at Publisher · View at Google Scholar · View at Scopus - P. Zilla, M. Deutsch, J. Meinhart et al., “Clinical in vitro endothelialization of femoropopliteal bypass grafts: an actuarial follow-up over three years,”
*Journal of Vascular Surgery*, vol. 19, no. 3, pp. 540–548, 1994. View at Google Scholar · View at Scopus - Y. Mohammed and J. F. Verhey, “A finite element method model to simulate laser interstitial thermotherapy in anatomical inhomogeneous regions,”
*BioMedical Engineering Online*, vol. 4, article 2, 2005. View at Publisher · View at Google Scholar · View at Scopus - F. Yu, Y. Shi, W. Deng, H. Chen, Q. An, and Y. Guo, “Particle image velocimetry in measuring the flow fields distribution in carotid artery bifurcation model,”
*Sheng Wu Yi Xue Gong Cheng Xue Za Zhi*, vol. 24, no. 1, pp. 104–109, 2007 (Chinese). View at Google Scholar · View at Scopus - L. H. Zhang and B.W. Liu, “Isolation of human plasma lipoprotiens by one-time density gradient ultracentrifugation method,”
*Acta Biochimica et Biophysica Sinica*, vol. 21, pp. 257–260, 1989. View at Google Scholar - P. J. Barter and K. A. Rye, “Molecular mechanisms of reverse cholesterol transport,”
*Current Opinion in Lipidology*, vol. 7, no. 2, pp. 82–87, 1996. View at Google Scholar · View at Scopus - E. Boisfer, D. Stengel, D. Pastier et al., “Antioxidant properties of HDL in transgenic mice overexpressing human apolipoprotein A-II,”
*The Journal of Lipid Research*, vol. 43, no. 5, pp. 732–741, 2002. View at Google Scholar · View at Scopus - P. Dimayuga, J. Zhu, S. Oguchi et al., “Reconstituted HDL containing human apolipoprotein A-1 reduces VCAM-1 expression and neointima formation following periadventitial cuff-induced carotid injury in apoE null mice,”
*Biochemical and Biophysical Research Communications*, vol. 264, no. 2, pp. 465–468, 1999. View at Publisher · View at Google Scholar · View at Scopus - H. Hattori, T. Kujiraoka, T. Egashira et al., “Association of coronary heart disease with pre-
*β*-HDL concentrations in Japanese men,”*Clinical Chemistry*, vol. 50, no. 3, pp. 589–595, 2004. View at Publisher · View at Google Scholar · View at Scopus - X. P. Li, S. P. Zhao, X. Y. Zhang, L. Liu, M. Gao, and Q. C. Zhou, “Protective effect of high density lipoprotein on endothelium-dependent vasodilatation,”
*International Journal of Cardiology*, vol. 73, no. 3, pp. 231–236, 2000. View at Publisher · View at Google Scholar · View at Scopus - M. Navab, S. Y. Hama, G. M. Anantharamaiah et al., “Normal high density lipoprotein inhibits three steps in the formation of mildly oxidized low density lipoprotein: steps 2 and 3,”
*The Journal of Lipid Research*, vol. 41, no. 9, pp. 1495–1508, 2000. View at Google Scholar · View at Scopus - J. R. Nofer, B. Levkau, I. Wolinska et al., “Suppression of endothelial cell apoptosis by high density lipoproteins (HDL) and HDL-associated lysosphingolipids,”
*The Journal of Biological Chemistry*, vol. 276, no. 37, pp. 34480–34485, 2001. View at Publisher · View at Google Scholar · View at Scopus - I. S. Yuhanna, Y. Zhu, B. E. Cox et al., “High-density lipoprotein binding to scavenger receptor-BI activates endothelial nitric oxide synthase,”
*Nature Medicine*, vol. 7, no. 7, pp. 853–857, 2001. View at Publisher · View at Google Scholar · View at Scopus - J. Hwang, M. H. Ing, A. Salazar et al., “Pulsatile versus oscillatory shear stress regulates NADPH oxidase subunit expression: implication for native LDL oxidation,”
*Circulation Research*, vol. 93, no. 12, pp. 1225–1232, 2003. View at Publisher · View at Google Scholar · View at Scopus - D. W. Crawford, L. H. Back, and M. A. Cole, “In vivo oxygen transport in the normal rabbit femoral arterial wall,”
*The Journal of Clinical Investigation*, vol. 65, no. 6, pp. 1498–1508, 1980. View at Google Scholar · View at Scopus - J. A. Berliner, M. C. Territo, A. Sevanian et al., “Minimally modified low density lipoprotein stimulates monocyte endothelial interactions,”
*The Journal of Clinical Investigation*, vol. 85, no. 4, pp. 1260–1266, 1990. View at Google Scholar · View at Scopus - J. Chen, J. L. Mehta, N. Haider, X. Zhang, J. Narula, and D. Li, “Role of caspases in Ox-LDL-induced apoptotic cascade in human coronary artery endothelial cells,”
*Circulation Research*, vol. 94, no. 3, pp. 370–376, 2004. View at Publisher · View at Google Scholar · View at Scopus - D. J. Gordon, J. L. Probstfield, R. J. Garrison et al., “High-density lipoprotein cholesterol and cardiovascular disease. Four prospective American studies,”
*Circulation*, vol. 79, no. 1, pp. 8–15, 1989. View at Google Scholar · View at Scopus - P. W. Wilson, K. M. Anderson, W. P. Castelli, and W. B. Kannel, “Twelve-year incidence of coronary heart disease in middle-aged adults during the era of hypertensive therapy: the Framingham Offspring Study,”
*The American Journal of Medicine*, vol. 90, no. 1, pp. 11–16, 1991. View at Publisher · View at Google Scholar · View at Scopus