Total cardiovascular risk for next 10 years among rural population of Nepal using WHO/ISH risk prediction chart
© The Author(s) 2017
Received: 16 November 2016
Accepted: 24 February 2017
Published: 7 March 2017
Cardiovascular diseases (CVD) are the leading cause of morbidity and mortality globally. Primary prevention of CVD based on total CVD risk approach using WHO/ISH risk prediction chart would be more effective to stratify population under different risk levels, prioritize and utilize the scarce resources of low and middle-income countries. This study estimated total 10-year CVD risk and determined the proportion of population who need immediate drug therapy among the rural population of Nepal.
A community based cross-sectional study conducted among 345 participants aged 40–80 years in rural villages of Lamjung District of Nepal. They were selected randomly from total eighteen wards. Data were collected using WHO STEPS questionnaires. WHO/ISH risk prediction chart for SEAR D was used to estimate total cardiovascular risk. Chi-square and independent t-test were used to test significance at the level of p < 0.05 in SPSS version 16.0.
Of the total participants, 55.4% were female. The mean age (standard deviation) of the participants was 53.5 ± 10.1 years. According to WHO/ISH chart proportions of low, moderate and high CVD risk were 86.4%, 9.3%, and 4.3%, respectively. Eleven percent of participants were in need of immediate pharmacotherapy. Age (p = 0.001), level of education (p = 0.01) and occupation (p = 0.001) were significantly associated with elevated CVD risk.
A large proportion of Nepalese rural population is at moderate and high CVD risk. Immediate pharmacological interventions are warranted for at least one in every ten individuals along with lifestyle interventions. Both population-wise and high-risk approaches are required to minimize CVD burden in the future.
KeywordsCardiovascular diseases WHO/ISH chart Risk prediction Total 10-year CVD risk Rural population
Cardiovascular diseases (CVD) are the major public health problem globally with disproportionately higher occurrence and consequences in developing countries. Data of 2012 showed estimated deaths due to CVD was 17.5 million (31% of total death) with 37% of premature deaths below the age of 70, globally. Almost 80% of such death occurred in low and middle-income countries (LMIC) with the high disparity of premature deaths among such LMIC . In another way, it has been projected that 23.3 million people would die by 2030 only due to CVD . Even in South–East Asia Region (SEAR), CVD accounted for a quarter (25%) of all deaths due to Noncommunicable Diseases (NCD) . Likewise, in Nepal, the prevalence of CVD was 40% among all NCD in-patients in non-specialty hospitals .
Population-wise strategies and high-risk strategies are two important methods for prevention of CVD complementing each other . The high-risk strategy identifies individuals who need immediate and sustained care to prevent the development of CVD . The high-risk strategy can be applied through two distinct approaches, first one is targeting and managing each of the single risk factors like hypertension, hyper-cholesterol separately. Another is accounting the total risk of CVD by measuring the cumulative effect of relative risk of all risk factors . However, later one seems to be accurate for identifying a high-risk individual and can enable to guide limited healthcare resources of LMIC cost-effectively [7–9].
Established tools to estimate the total CVD risk are Framingham, SCORE, QRISK, ASSIGN. These tools were developed considering the specific population and might not be applicable to others [10–15]. However, World Health Organization/International Society of Hypertension (WHO/ISH) risk prediction charts have been developed using information on the population distribution of risk factors from WHO Comparative Risk Assessment study of different sub-regions focusing non western developing countries [16, 17]. The chart estimates the total risk of developing fatal and non-fatal cardiovascular diseases (coronary artery disease and stroke) in following 10 years among those people whose CVD were not manifested clinically, categorizes an individual in different risk groups so that the person can be managed only by lifestyle modification or need drug therapy also . The chart is easy to utilize by all physicians and other health workers at all levels of health care. This tool not only helps to make a clinical decision but also estimates and monitor population distribution of total CVD risk [16, 18]. However, there is paucity of evidence on the use of WHO/ISH risk prediction chart in Nepal. Therefore, this study aimed to estimate total 10-year CVD risk and to determine the proportion of population who need immediate drug therapy among the rural community of Nepal using WHO/ISH chart.
Study design and settings
This was a community-based cross-sectional study conducted at Bhotewodar and Sunderbazar Village Development Committee (VDC) of Lamjung between October and November 2014. These villages are about 200 km west from capital city Kathmandu of Nepal and consist of 14,275 population living in 3,869 households. These are the hilly rural areas where agriculture is the main source of livelihood. The sites are the inhabitants of culturally rich people like Gurung, Tamang, Magar, Newar and Dalit, most of them are socially deprived .
All eligible men and women were included as the study population. Eligible criteria were age 40–80 years, without having self-reported angina, myocardial infarction, stroke, and intermittent claudication. Pregnant women, mentally retarded and bed-ridden patients were excluded from the study.
Sample size and sampling technique
The sample size was calculated by using the formula for estimation of proportion for one sample situation. To detect prevalence of 9.5% of moderate and high CVD risk population, as determined by a study conducted in sub-urban municipalities of Nepal , the minimum sample size required was 292 with allowable error of 4% and a design effect of 1.5. Taking into consideration a non-response rate of 15%, we decided to select 343 subjects from all male and female population aged 40–80 years (N = 3712).
A cluster sampling technique was used to recruit the required number of participants for the study. Firstly, four wards were selected out of 18 wards using probability proportion to size (PPS). In each ward, we selected the first household from a randomly selected point at the entry of the main road. Then, every 6th household was taken as the study household. All eligible members of each household were listed and one subject was selected randomly for interview by using Kish Table .
Five health professionals who had at least intermediate level of health-related education were trained for one week on study protocol, tools and data collection techniques of the study. Data were collected using standard STEPS instrument version 2.2 . The survey tool encompassed questions related to demography and socio-economy, tobacco and alcohol use, dietary intake, physical activity level, and history of high blood pressure and diabetes. The questionnaires were translated into Nepali language and revised after conducting a pilot survey among 20 participants. Blood pressure was measured using doctor’s aneroid sphygmomanometer and stethoscope with an appropriately sized cuff on the unclothed left arm. First reading was taken after 15 min rest-crossing leg and second and third readings were recorded at 5 min interval. Trained technical person collected venous blood, after conforming overnight 12 h fasting to measure blood sugar and lipid profile. All biochemical measurements were done using semi-automated machine of Erba Mannheim (Germany), chem5 v3 model.
WHO/ISH risk prediction chart
WHO/ISH chart for SEAR D of WHO epidemiological sub-region  was used to estimate the total 10-year risk of CVD of all participants. Age, gender, smoking status, systolic blood pressure, presence or absence of diabetes and total serum cholesterol in millimole/liter (mmol/l) were used to compute the total CVD risk. The chart can stratify an individual into low (<10%), moderate (10% to <20%), high (20% to <30%), and very high (>30%) risk groups .
Definitions of variables
A smoker was the one who has smoked currently and who quit smoking less than one year before the assessment . Systolic blood pressure was mean of last two readings. A respondent with diabetes mellitus was defined as the one who had fasting venous blood sugar level 7 mmol/l or who was taking oral hypoglycemic drug or insulin [23, 24]. Caste was categorized into upper caste (Bramhan and Kshetri), Janajati (Gurung, Tamang, Newar, Magar and Thakali) and Dalit (low caste). For occupation, respondent working in paid governmental and nongovernmental organization was labeled as employed, running own business as self-employed, and student or housemaker or non-paid worker as unemployed. Economic status was categorized as poor (<8498) and medium to high (≥8498) economic status according to family income per person per month in Nepalese Rupees (NR) .
Data were compiled and entered in Epidata version 3.1. Duplication and omissions of data were checked and edited before coding and exporting into SPSS V.16.0 for further analysis. Descriptive statistics were used to compute different risk categories with 95% confidence interval (CI) for proportion. Mean was calculated with standard deviation (SD) for continuous variables. Chi-square and independent t-tests were conducted to compare categorical and continuous variables, respectively. All tests were two-tail test and p < 0.05 was considered statistically significant.
Ethical clearance was obtained from ethical review committee of Nepal Health Research Council. Written informed consent or thumb impression (if unable to write) were taken from all participants. Study objectives, data collection procedures, benefits and risks of the study, confidentiality, and anticipated use of the results were explained to all participants in detail before taking consent.
Socio-demographic characteristics of respondents
Male, n (%)
Female, n (%)
Total, n (%)
Age in years
Mean ± SDb
53.3 ± 9.9
53.7 ± 10.3
53.5 ± 10.1
Widow or widower
Level of educationa
No formal education
Medium and high
Total 10-year CVD risk using WHO/ISH chart
Distribution of study population into low, moderate and high cardiovascular risk categories, n (%)
Male and Female
Socio-demographic factors associated with elevated 10-year CVD risk, n (%)
<10% risk (low)
>10% risk (elevated)
Level of education
No formal education
Widow or widower
Poor (≤NRs 8498)
Higher (>NRs 8498)
Socio-demographic factors associated with elevated 10-year CVD risk
CVD risk was significantly (p = 0.01) associated with the level of education. Of participants without formal education, 18.4% had elevated CVD risk. However, it was 10% and 6.7% of participants who had a primary and secondary level of education. Elevated CVD risk was higher among Janajati and Dalit (14.7% and 14.3%) in comparison to upper caste participants (12%). Participants whose partner was not alive (30%) had higher CVD risk compared to married (12.1%). Occupation was significantly (p = 0.001) associated with CVD risk. Of students, homemaker or non-paid workers, 16% had elevated CVD risk. Similarly, among retired participants, 37% had more than 10% CVD risk. However, the elevated risk was only 1.3% and 10% among self-employed and employed participants, respectively (Table 3).
Proportion of population who need immediate drug therapy
Study population requiring pharmacotherapy using WHO/ISH chart at different threshold levels and single risk factors approach, n (%)
Persons requiring drugs
Threshold CV risk >30%
Threshold CV risk >20%
Chart + BPa
Single risk factor
BP ≥ 140/90
TC ≥ 6
Our study had an objective to estimate total 10-year CVD risk among study population. This study found that large proportion of study participants had moderate to high risk (13.6%) which was consistent with the findings reported in national and international studies conducted in similar settings. This proportion of population at elevated risk of CVD was almost same in one sample of Nepal (13.4%) reported by Mendis et al. , while higher than the finding of another study conducted in a sub-urban population of Nepal (9.5%) . However, the proportion of raised risk in current study was slightly lower than that of Kathmandu (14.6%)  and rural India (16%) . Regarding the proportion of population under high risk of CVD, it was slightly higher in our study (4.3%) than that of study conducted in Malaysia (2.3%) and Cambodia (1.3%) and lower than that of Mangolian study (6%) . Similar heterogeneous results were observed in different countries (China 1.1%, Iran 1.7%, Sri Lanka 2.2%, Cuba 2.8%, Nigeria 5.0%, Georgia 9.6%, Pakistan 10.0%) . Very high-risk (>30%) was 1.7% in our study, which was slightly lower than the estimated for age 40–65 years throughout the country of Nepal (3.2%) . These different data revealed that the population under risk for development of stroke and coronary heart disease were high even in a rural community of Nepal. Application of WHO/ISH chart helps to categorize population into different risk levels. Population-based lifestyle modification strategy can be applied to low-risk population while individual counseling and frequent follow-up assessment need to be added for moderate risk population. However, more rigorous treatment strategies would be needed for the population at high and very high-risk category .
CVD risk was significantly associated with the level of education. A study conducted in Vienna echoed the same finding that education is associated with the level of understanding on risk factors . Evidence also suggests that increasing knowledge on risk factors can substantially motivate people for changing risk behaviors . This fact suggests that we need to improve awareness and knowledge about prevention of cardiovascular diseases. Occupation was significantly associated with elevated CVD risk. A higher proportion of unemployed and retired participants had raised risk. One study of Mauritius also proved the association of occupation and cardiovascular risk, mostly related to physical activities and economic status . These findings provide the evidence that we should target unemployed population for early prevention of cardiovascular disease. As expected, age-related increase in CVD risk was associated with an increase in risk factor level .
As explained in the WHO guideline for assessment and management of cardiovascular disease , individuals who have blood pressure of ≥160/100 mm Hg need appropriate drug therapy even though they are in low and moderate CVD risk category. However, there is still no universal consensus for appropriate threshold (30% vs 20%) of total CVD risk at which drug treatment are initiated [7, 16]. It depends on the availability of resource of a country and impact of the specific intervention. Current data reveals that we can prevent large numbers of ischemic heart disease and stroke incidence if threshold is reduced from >30 to >20% risk. Our result demonstrated that about 11% would be treated by drug if >20% risk was taken as threshold while about 8.6% would need drug at threshold of >30% risk. However, about threefold more individual would require pharmacotherapy if we target single risk factors like hypertension and raised cholesterol. Our estimation of proportion of the population in need of drug therapy was slightly lower than that of Bangladesh at same threshold level (8.3% and 4%). However, as in our study, population who require drug were almost threefold higher by targeting single risk factors (28.2%) . Our estimation was comparable to another study conducted at Cuba , where sample population requiring drug at aforementioned threshold level were 12.5% and 11.5% with total risk approach restricting about threefold of the population from over prescription of drug than single risk factor approach (43.8%). Similar result was found in a study in Africa [7, 9]. Therefore, high proportion of the population were in need of pharmacotherapy and it would prevent over prescription of the drug if total CVD risk approach applied to prevent CVD.
Due to lack of data related to cardiovascular morbidity and mortality, WHO/ISH charts have not validated yet in Nepal. The chart might have underestimated the total CVD risk, as it did not include the individual with antihypertensive therapy, raised triglyceride, low HDL, family history of CVD and obesity. Higher number of janajati participants were involved in the study. Therefore, findings of the study should cautiously be generalized to other rural population of Nepal.
A large proportion of Nepalese rural population is at moderate and high CVD risk. Immediate pharmacological interventions are warranted for at least one in every ten individuals along with lifestyle interventions. This situation underlines the importance of immediate designing and implementation of both population-wise and high-risk approaches to minimize CVD burden in the future. Total risk approach would be more effective for initiation of preventive pharmacotherapy compared to the single risk factor approach to allocate the limited resources more efficiently.
Low and Middle Income Countries
World Health Organization/International Society for Hypertension
Stepwise Approach to Surveillance
South East Asia Region
MKK designed the main thesis, trained field enumerators, analysed data, interpreted data and drafted the manuscript; MSAMA oversaw the quality of the main thesis and reviewed the manuscript; MZ interpreted results critically, conceptualized the manuscript and critically reviewed it; PCB and PB collected, compiled, analyzed, interpreted the data and involved in drafting the manuscript; RD interpreted the data, reviewed the draft and provided substantial comments; SD and AS critically revised the main thesis and manuscript. All authors read and approved the final manuscript.
This article is based on the thesis which was submitted to Bangladesh Institute of Health Sciences (BIHS) as a partial fulfillment of MPhil degree in Noncommunicable Diseases under Dhaka University. We would like to thank the study teams and study participants. We acknowledge Mr. Ramchandra Adhikari, managing director of Karmada Hospital Pvt. Ltd. for providing laboratory, logistics and transportation facilities. We would like to thank lab technician Mr. Bikash Siluwal for blood sample collection and analysis. Mr. Shailendra Pandit, Miss Pratiksha Adhikari, Miss Yem Kumari Gurung, Dr. Anil Koju, Dr. Krishna Malakar, Mr. Dwarika Mishra, Janaklal Shreshtha, Dr. Sabrina Ahmed, Dr. Debdeep Mukherjee, Meena Iqbal and secretaries of Bhowtewodar and Sundarbazar VDC who played key role in data collection and data management.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Ethical clearance was obtained from the ethical review committee of Nepal Health Research Council. Written informed consent (or thumb impression if unable to write) were taken using a consent form in Nepali language from all participants.
This research work (M.Phil Thesis) was funded by the ‘Norad’s Programme for Master Studies (NOMA)’ grant of the University of Oslo, Norway.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- World Health Organization. Global status report on noncommunicable diseases 2014. Geneva: World Health Organization; 2014.Google Scholar
- Mathers CD, Loncar D. Projections of global mortality and burden of disease from 2002 to 2030. PLoS Med. 2006;3(11):e442.View ArticlePubMedPubMed CentralGoogle Scholar
- World Health Organization. Non communicable diseases in the South East Asia region: situation and response 2011. Geneva: World Health Organization; 2011.Google Scholar
- Bhandari G, Angdembe M, Dhimal M, Neupane S, Bhusal C. State of non-communicable diseases in Nepal. BMC Public Health. 2014;14(1):23.View ArticlePubMedPubMed CentralGoogle Scholar
- Emberson J, Whincup P, Morris R, Walker M, Ebrahim S. Evaluating the impact of population and high-risk strategies for the primary prevention of cardiovascular disease. Eur Heart J. 2004;25(6):484–91.View ArticlePubMedGoogle Scholar
- Hobbs FDR. Cardiovascular disease: different strategies for primary and secondary prevention? Heart. 2004;90(10):1217–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Mendis S, Lindholm LH, Anderson SG, Alwan A, Koju R, Onwubere BJ, Kayani AM, Abeysinghe N, Duneas A, Tabagari S, et al. Total cardiovascular risk approach to improve efficiency of cardiovascular prevention in resource constrain settings. J Clin Epidemiol. 2011;64(12):1451–62.View ArticlePubMedGoogle Scholar
- Batsis J, Lopez-Jimenez F. Cardiovascular risk assessment—from individual risk prediction to estimation of global risk and change in risk in the population. BMC Med. 2010;8(1):29.View ArticlePubMedPubMed CentralGoogle Scholar
- Ndindjock R, Gedeon J, Mendis S, Paccaud F, Bovet P. Potential impact of single-risk-factor versus total risk management for the prevention of cardiovascular events in Seychelles. Bull World Health Organ. 2011;89(4):286–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson PWF, D’Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB. Prediction of coronary heart disease using risk factor categories. Circulation. 1998;97(18):1837–47.View ArticlePubMedGoogle Scholar
- Kannel WB, McGee D, Gordon T. A general cardiovascular risk profile: the Framingham Study. Am J Cardiol. 1976;38(1):46–51.View ArticlePubMedGoogle Scholar
- D’Agostino RB Sr, Vasan RS, Pencina MJ, Wolf PA, Cobain M, Massaro JM, Kannel WB. General cardiovascular risk profile for use in primary care: the Framingham Heart Study. Circulation. 2008;117(6):743–53.View ArticlePubMedGoogle Scholar
- Conroy RM, Pyorala K, Fitzgerald AP, Sans S, Menotti A, De Backer G, De Bacquer D, Ducimetiere P, Jousilahti P, Keil U, et al. Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE project. Eur Heart J. 2003;24(11):987–1003.View ArticlePubMedGoogle Scholar
- Julia H-C, Carol C, John R, Peter B. Derivation, validation, and evaluation of a new QRISK model to estimate lifetime risk of cardiovascular disease: cohort study using QResearch database. BMJ. 2010;341:c6624.View ArticleGoogle Scholar
- Woodward M, Brindle P, Tunstall-Pedoe H. Adding social deprivation and family history to cardiovascular risk assessment: the ASSIGN score from the Scottish Heart Health Extended Cohort (SHHEC). Heart. 2007;93(2):172–6.View ArticlePubMedGoogle Scholar
- World Health Organization. Prevention of cardiovascular disease: Guideline for assessment and management of cardiovascular risk. Geneva: World Health Organization; 2007.Google Scholar
- Mendis S, Lindholm LH, Mancia G, Whitworth J, Alderman M, Lim S, Heagerty T. World Health Organization (WHO) and International Society of Hypertension (ISH) risk prediction charts: assessment of cardiovascular risk for prevention and control of cardiovascular disease in low and middle-income countries. J Hypertens. 2007;25(8):1578–82.View ArticlePubMedGoogle Scholar
- Otgontuya D, Oum S, Buckley B, Bonita R. Assessment of total cardiovascular risk using WHO/ISH risk prediction charts in three low and middle income countries in Asia. BMC Public Health. 2013;13(1):539.View ArticlePubMedPubMed CentralGoogle Scholar
- Central Bureau of Statistics. National Population and Housing Census 2011, vol. 05, 1st edn. Kathmandu: Central Bureau of Statistics; 2014.Google Scholar
- Koju R, Gurung R, Pant P, Humagain S, Yogol CM, Koju A, Manandhar K, Karmacharya B, Bedi TRS. Prediction of cardiovascular disease in suburban population of 3 municipalities in Nepal. Nepal Heart J. 2011;8(1):3–7.Google Scholar
- World Health Organization. WHO STEPS surveillance manual: the WHO STEPwise approach to chronic disease risk factor surveillance/noncommunicable diseases and mental health. Geneva: World Health Organization; 2005.Google Scholar
- WHO/ISH risk prediction charts for 14 WHO epidemiological sub-regions. http://ish-world.com/downloads/activities/colour_charts_24_Aug_07.pdf.
- World Health Organization. Prevention of cardiovascular diseases: pocket guidelines for assessment and management of cardiovascular risk. Geneva: World Health Organization; 2007.Google Scholar
- World Health Organization. Definition and diagnosis of diabetes mellitus and intermediate hyperglycemia: report of a WHO/IDF consultation. Geneva: World Health Organization; 2006.Google Scholar
- Central Bureau of Statistics. Nepal living standards survey 2010/11. Kathmandu: Central Bureau of Statistics; 2011.Google Scholar
- Dhungana RR, Khanal MK, Pandey AR, Thapa P, Devkota S, Mumu SJ, Shayami A, Ali L. Assessment of short term cardiovascular risk among 40 years and above population in a selected community of Kathmandu, Nepal. J Nepal Health Res Counc. 2015;13(29):66–72.PubMedGoogle Scholar
- Ghorpade AG, Shrivastava SR, Kar SSSS, Majgi SM, Roy G. Estimation of the cardiovascular risk using World Health Organization/International Society of Hypertension (WHO/ISH) risk prediction charts in a rural population of South India. Int J Health Policy Manag. 2015;4(8):1–6.View ArticleGoogle Scholar
- Aryal KN, Mehata S, Vaidya A, Singh S, Paulin F, Madanlal RG, Riley LM, Cowan M, Guthold R, Singh SP, Bhusal CL, Lohani GR. Non communicable diseases risk factors: STEPS Survey Nepal 2013. Kathmandu; 2014.
- Samal D, Greisenegger S, Auff E, Lang W, Lalouschek W. The relation between knowledge about hypertension and education in hospitalized patients with stroke in Vienna. Stroke. 2007;38(4):1304–8.View ArticlePubMedGoogle Scholar
- McDermott MM, Mandapat AL, Moates A, Albay M, Chiou E, Celic L, Greenland P. Knowledge and attitudes regarding cardiovascular disease risk and prevention in patients with coronary or peripheral arterial disease. Arch Intern Med. 2003;163(18):2157–62.View ArticlePubMedGoogle Scholar
- Pereira MA, Kriska AM, Collins VR, Dowse GK, Tuomilehto J, Alberti KG, Gareeboo H, Hemraj F, Purran A, Fareed D, et al. Occupational status and cardiovascular disease risk factors in the rapidly developing, high-risk population of Mauritius. Am J Epidemiol. 1998;148(2):148–59.View ArticlePubMedGoogle Scholar
- Jousilahti P, Vartiainen E, Tuomilehto J, Puska P. Sex, age, cardiovascular risk factors, and coronary heart disease: a prospective follow-up study of 14,786 middle-aged men and women in Finland. Circulation. 1999;99(9):1165–72.View ArticlePubMedGoogle Scholar
- Fatema K, Zwar NA, Milton AH, Rahman B, Ali L. Application of two versions of the WHO/international society of hypertension absolute cardiovascular risk assessment tools in a rural Bangladeshi population. BMJ Open. 2015;5(10):e008140.View ArticlePubMedPubMed CentralGoogle Scholar
- Nordet P, Mendis S, Duenas A, de la Noval R, Armas N, de la Noval IL, Pupo H. Total cardiovascular risk assessment and management using two prediction tools, with and without blood cholesterol. MEDICC Rev. 2013;15(4):36–40.PubMedGoogle Scholar