The objective of the current study was to assess the secular trends in the prevalence of diabetes, prediabetes, and risk factors from two epidemiological surveys done 10 years apart in three adult populations of different geographic and socioeconomic backgrounds in Tamil Nadu, India.
This survey was conducted in 2016 using methodology similar to that used in 2006. Persons aged ≥20 years (n = 9,848) were screened for diabetes, prediabetes, and the risk variables. Fasting and 2-h plasma glucose, lipid profile, blood pressure, anthropometry, and socioeconomic and behavioral details were recorded. Comparative analyses of age-standardized prevalence were done. Prevalence ratios (PRs) between 2016 and 2006 of diabetes and also prediabetes were assessed using Poisson regression analyses.
Prevalence of diabetes increased from 18.6% (95% CI 16.6–20.5) to 21.9 (20.5–23.3) in the city, 16.4 (14.1–18.6) to 20.3 (18.9–21.6) in the town, and 9.2 (8.0–10.5) to 13.4 (11.9–14.8) in the periurban villages (PUVs) (P < 0.0001 in all). The PR showed a nonsignificant 8% rise in diabetes in the city, while significant increases had occurred in the town (39%) and PUVs (34%). Prevalence of prediabetes also increased. Age, family history of diabetes, and waist circumference were common risk determinants among the populations. Though general obesity and abdominal obesity increased, the latter was associated with the increased prevalence.
Prevalence of diabetes and prediabetes increased in all locations; the rise was significant only in the town and PUVs. Abdominal obesity is significantly associated with increased trend even among the villagers. Rural populations may be targeted for future public health measures to combat diabetes.
Introduction
The latest estimate by the International Diabetes Federation in 2017 shows that the prevalence of diabetes among the adults (20–79 years of age) in India is 72.9 million (1). A phenomenal increase in diabetes in the last 10 years, from 40.9 million in 2007 to 72.9 million in 2017 (1.78-fold increase), occurred. The prevalence is likely to increase further by 1.84-fold, i.e., to 134.3 million by 2045. Once again, India is likely to have the largest number of people with diabetes in the world (1). The prevalence of impaired glucose tolerance (IGT) is 2.9% (24 million) (1).
Asian countries, specifically, India and China, are undergoing rapid socioeconomic progress and are susceptible to many adverse consequences such as unhealthy lifestyle leading to high prevalence of diabetes and other noncommunicable diseases (2).
Epidemiological studies in India since 1980 have highlighted the changing trends in prevalence of diabetes, prediabetes, and associated risk factors (3–6). The prevalence of diabetes is higher in urban areas. The determinants of the rising prevalence vary among populations and with time. Unhealthy changes in lifestyle associated with urbanization, familial aggregation, and aging of the population are common risk associations. Large numbers of Indians have migrated to many Western and Asian countries, and the prevalence of diabetes among them is much higher than that of the host populations (2,7,8). In India (9–12) and other developing countries (13–15), lifestyle transition is a major determinant of the rising trends in diabetes, even in rural areas. The phenomenal increase in diabetes can be attributed to changes in dietary pattern, sedentary behavior, and obesity superimposed on the background of genetic and epigenetic susceptibility (16).
Regional disparities in dietary pattern, occupation, and demography are likely to influence the prevalence rates. The extent to which the risk factors are changing in different geographical areas in India and its influence on diabetes is not known. There are sparse data on the time-trend comparisons using standardized methods from low- and middle-income countries, which contribute to nearly 79% of the global prevalence of diabetes (1). Although there are studies showing the epidemiological transition in different Asian countries, only a limited number of studies (5,17) have compared the secular trends in diabetes within the same population at different time points using similar methodology. Such studies from the same region can help with understanding the temporal changes occurring in a specific population. Hence, the current study was done to assess the secular trends in the prevalence of diabetes and prediabetes and the changes in the associated risk factors in urban and rural areas, i.e., in a city, town, and periurban villages (PUVs) in the state of Tamil Nadu in South India, in a 10-year period (2006–2016). The census of 2011 for India indicated that 48.5% of Tamil Nadu is urbanized and has the highest rate of urbanization (18). We did cross-sectional epidemiological studies in the same city, town, and PUVs in 2006 (10) and in 2016 using similar methodology, and the changing trends are reported in this article.
Research Design and Methods
Study Design and Participants
The present epidemiological surveys were conducted from July 2016 to September 2017 in three geographical locations in Tamil Nadu, South India, where similar surveys had been done in 2006. Three locations, classified as city, town, and PUVs by the Indian census commission, were chosen in 2006. In India as per the census commission, a city has a population of >4 million and a town is smaller than a city, with a minimum population of 100,000 and at least 75% of the male working population engaged in nonagricultural pursuits. PUVs are defined as landscapes between a town and rural area. People in PUVs have low economic status and are mostly casual workers engaged in agricultural work. In 2016, the same locations were chosen: city, Chennai; town, Kanchipuram (80 km from Chennai); and PUVs, Panruti (186 km from Chennai). The methodologies used in selection of the study samples, screening, and biochemical procedures were similar in both surveys. For Chennai, the required sample size was 3,824, at 80% power with an α error of 5% and a design effect of 1.5 and an assumed increase in prevalence from 2006 to 2016 from 18.6 to 21.0%. In the town, the number was 3,547 assuming an increase from 16.4 to 19.0%, and the number for PUVs was 2,434 with an expected increase from 9.2 to 12.2%. Samples larger than the calculated numbers were studied. A total of 9,848 were tested (city, n = 3,850; town, n = 3,530; and PUVs, n = 2,468). Response rates were 86.5% in the city, 88.1% in the town, and 87.6% in the PUVs. A multistage random selection was done from streets with population characteristics similar to those in 2006, which gave a fair representation of all the socioeconomic strata. In the city, 5 corporation zones out of 10 marked in 2006 were randomly selected, in the town the selection was done from municipal wards, and in the PUVs panchayat census wards were used. We selected 18 villages in Panruti, which had a population ranging from 2,054 to 8,313 persons. All subjects aged ≥20 years were eligible to participate and were enumerated and invited by house visits. Families in which all eligible members gave written informed consent for participation were selected. There was no significant difference in the demographic characteristics of responders and nonresponders.
The study protocol was approved by the Ethics Committee of India Diabetes Research Foundation (IDRF) and Dr. A. Ramachandran’s Diabetes Hospitals. An independent data-monitoring committee reviewed the progress of the survey. Written permission from the Department of Public Health and Preventive Medicine, Government of Tamil Nadu, was obtained to conduct the survey in the selected locations. Periodic inspections were done by these authorities while the survey was being conducted. The study was registered on www.ClinicalTrials.gov (identifier: NCT03490136).
Demographic, Anthropometric, and Behavioral Assessment
Details of demography, anthropometry, education, monthly family income, occupation, physical activity, diet habits, smoking, and alcohol consumption were filled in a questionnaire by personal interview by trained team members. Height, weight, and waist circumference (WC) were measured by standard procedures (10). BMI was calculated as weight in kilograms divided by the square of height in meters. Generalized obesity was defined as BMI ≥23.0 kg/m2. BMI between ≥23.0 and 24.9 kg/m2 was considered overweight, and BMI ≥25.0 kg/m2 was defined as obese. Abdominal obesity was indicated by WC of ≥90 cm for men and ≥80 cm for women (19).
Nutrient intake was recorded by trained dietitians using the 24-h dietary recall method. The total energy intake (kcal) and components of individual food constituents (carbohydrates, proteins, and fat [grams] were calculated with an in-house dietary analysis program (visual basic programming tool) using the National Institute of Nutrition guidelines for India. Details of physical activity were recorded and quantified using a questionnaire originally used in U.K. Asian Indians (20) that was modified to suit the Indian environment (10,21,22). The quantification was based on the occupational activity and hours of moderate/vigorous activity and leisure-time activity. A score of 7–70 was used.
The economic status of the family was assessed by monthly family income in Indian rupees. The family income was divided in tertiles representing the lower-, middle-, and upper-income groups. The income data of the population in 2006 were also similarly categorized. Education level was categorized as no formal education, school education, and college education. Categorizes of occupation were laborers, homemakers/business, and executives/professional/clerical/students/retired.
Clinical Assessment
Each day 40–60 persons were tested. A clinician took down the medical history and measured the blood pressure in the sitting position using the electronic measuring device (Omron HEM 7132; Omron, Tokyo, Japan). An average of two readings taken at 5-min intervals was recorded. Self-reported history of diabetes and hypertension with details of diagnosis and medication were noted. Subjects with a history of hypertension and those with newly diagnosed diabetes with blood pressure readings ≥140/90 mmHg were categorized as hypertensive (23). Newly diagnosed patients with diabetes during the survey were referred to their physicians or the local public health center or hospital for management.
Biochemical Assessment
Blood samples were collected by door-to-door visits. Fasting samples were taken between 6:00 and 8:00 a.m. For known cases of diabetes, fasting and 2-h postprandial blood glucose was measured. All other subjects underwent an oral glucose tolerance test after 8–12 h of fasting. The fasting state was ensured by questioning the participants. Blood glucose was measured at fasting and at 2 h after consumption of 75 g oral glucose, using venous blood and Accu-Chek Performa (Roche Diagnostics, GmbH, Germany) by the hexokinase method and calibrated to read plasma glucose (24). In this survey, plasma glucose levels were measured using a glucometer to avoid errors that could occur while transporting the samples from distant places to the laboratory. In a subsample of 111 subjects, in addition to the glucose estimation by glucometer, venous blood collected in sodium fluoride was separated at the work site and plasma brought to the central laboratory. Glucose estimations were done within 6 h of collection. The correlation was as follows: reported plasma glucose = 17.1 + 0.887 × venous plasma glucose (glucometer), P < 0.0001. The glucometer values were corrected using the equation.
Diagnosis of diabetes and IGT was made using World Health Organization criteria (25). Diabetes was diagnosed if fasting plasma glucose (FPG) value was ≥7 mmol/L (≥126 mg/dL) or 2-h post glucose ≥11.1 mmol/L (≥200 mg/dL), impaired fasting glucose (IFG) if FPG was between 6.1 and <7 mmol/L (110–125 mg/dL), and IGT if the 2-h values were ≥7.8–11.0 mmol/L (≥140–199 mg/dL) with fasting values <6.1 mmol/L (<110 mg/dL).
All blood samples were centrifuged within an hour of collection, and aliquots were stored in freezers until they were sent to the IDRF laboratory. Fasting serum lipid profile was estimated by standard enzymatic procedures (Roche Diagnostics). HDL cholesterol was estimated by the direct assay method.
Statistical Analysis
Estimates of prevalence of diabetes and prediabetes were age standardized by a direct standardization method using the 2011 census data for the respective populations in Tamil Nadu. Means ± SD are reported for normally distributed variables. Student t test and one-way ANOVA were used for group comparisons. For skewed variables, median and interquartile ranges are given. Mann-Whitney U test or Kruskal-Wallis H test was used as relevant for comparisons. Trend χ2 test or z test was used for comparing categorical variables.
Poisson regression analyses with robust variance were done for 2006 and 2016 to identify the risk variables associated with diabetes versus normoglycemia and prediabetes versus normoglycemia. Independent risk variables included in the models were age (units of 10), BMI (units of 5), WC (units of 10), sex (reference: female), family history of diabetes (reference: negative), monthly family income (three categories [reference: low income]), total calorie consumption (three categories [reference: low calorie]), physical activity (four categories [reference: heavy activity]), education (three categories [reference: no formal education]), occupation (three categories [reference: laborers]), smoking (reference: nil), and alcohol (reference: nil). The prevalence ratios (PRs) (2016 vs. 2006) of diabetes and prediabetes were calculated using Poisson regression analyses with corrections for the significant confounders identified in the above analyses.
Results
Table 1 shows the characteristics of the current study cohorts in comparison with the corresponding data in 2006. Ratio of men and women was representative of the population in each location. Significant increases in age, BMI, and WC were seen in all of the three populations. Systolic blood pressure showed an increase in the urban areas, and diastolic blood pressure increased in all areas. Lipid levels increased in all locations (P < 0.0001). However, the total cholesterol–to–HDL cholesterol ratio improved in PUVs (P < 0.0001). In 2016, the majority of the populations (63–74%) had school education, and percentages with college education improved in all locations including the PUVs (P < 0.0001). The proportion of laborers increased and percentage of homemakers and numbers doing business decreased both in the urban and rural areas.
In 2006 and 2016, the mean age of people with diabetes was similar in the city and in PUVs (49.4 ± 11.5 vs. 49.8 ± 10.8 years, P = 0.560, and 47.5 ± 10.8 vs. 47.7 ± 11.9 years, P = 0.531, respectively). In the town, it was higher in 2016 (51.9 ± 10.7 vs. 47.7 ± 9.7 years, P < 0.0001). The mean age of those with newly diagnosed diabetes in Chennai (45.2 ± 11.3 vs. 46.6 ± 10.9 years, P = 0.190) and PUV (43.9 ± 11.4 vs. 45.6 ± 11.9 years, P = 0.262) were similar in the 2006 and 2016 surveys, respectively. In the town, there was an increase in the mean age (44.6 ± 9.9 vs. 48.1 ± 10.9 years, P = 0.003).
Total calorie consumption increased significantly (P < 0.0001) in urban areas. In the city, fat consumption increased. No changes were seen in the dietary pattern in the PUVs. Moderate and heavy levels of physical activity were higher among persons in the PUVs (P < 0.0001).
Table 2 shows the comparison of glucose intolerance and prevalence of generalized obesity, abdominal obesity, and hypertension in the study periods.
Diabetes
The prevalence of diabetes was 21.9% (95% CI 20.5–23.3), 20.3 (18.9–21.6), and 13.4 (11.9–14.8) in the city, town, and PUVs, respectively. The prevalence had increased significantly compared with the previous survey (P < 0.003 in the city and P < 0.0001 in other locations). The increase was statistically significant among men only in the PUVs (P < 0.0001), but significant increase occurred in women in all locations (P < 0.0001). Self-reported diabetes increased in all locations. Percentage of new diabetes increased in both the urban and rural areas (P < 0.0001).
Prediabetes
Increasing trends in total prediabetes (isolated IFG [iIFG], isolated IGT [iIGT], and IFG + IGT) and its individual categories were noted in all locations (P < 0.0001). Prevalence of iIGT was more common than other forms of prediabetes. Maximum increase in prediabetes occurred in the town.
Compared with the data of 2006, significant increases in generalized obesity and abdominal obesity occurred in all three populations; the highest percentage of increase was observed in the PUVs (Table 2). In the urban areas, prevalence of hypertension increased significantly. Generalized obesity increased by 16.3, 26.4, and 45.4% in the city, town, and PUVs, respectively. Higher percentages of increase in abdominal obesity was seen in all populations (34, 35.2, and 57.2%). The highest increase in both variables occurred in the PUVs (P < 0.0001).
Table 3 shows the PRs of diabetes (2016 vs. 2006) and the risk associations for diabetes versus normoglycemia, analyzed by Poisson regression in the three locations. The PR for diabetes for the city (2016 vs. 2006) was not significant (PR 1.08 [95% CI −0.01 to 0.17]; P = 0.08), whereas the ratios were significant for the town (PR 1.39 [95% CI 1.25–1.55], P < 0.0001) and PUV (PR 1.34 [95% CI 1.15–1.54], P < 0.0001).
In all locations, age, positive family history of diabetes, and WC were associated with the increasing trend in diabetes. Male sex showed an association with the increase in the city and PUVs. BMI showed an association only in the city. In the city, people with higher education level had lower prevalence and persons having high-income occupations had higher prevalence. In PUVs, sedentary activity was associated with diabetes. Total dietary calorie consumption was not independently associated with diabetes. Smoking and alcohol consumption were not associated with diabetes.
Table 3 also shows the PR for prediabetes in 2016 versus 2006. The ratio increased significantly (P < 0.0001) in all populations; the highest increase occurred in the town (PR 4.01 [95% CI 3.31–4.86], P < 0.0001) followed by the PUVs (PR 1.87 [95% CI 1.59–2.20], P < 0.0001). Age was the only common risk association, and BMI was associated with the rising trend in the town. Family history and BMI were associated with the trend in the city. In the PUVs, age and WC showed a positive association and male sex had a higher prevalence.
Figure 1 shows that compared with prevalence in 2006, prevalence of diabetes increased in all age-groups in all three populations. The takeoff point for diabetes was between 25 and 34 years. In the city, the age-specific increase was not significant. In the town, an increasing trend was observed at all ages. At ≥64 years of age, the increase was highly significant (P < 0.0001), and the maximum increase was at ≥64 years (P < 0.022). A significant fall in the prevalence after the age of 64 years was present in the town in 2006. In the PUVs, the increase in 2016 was significant until the age of 64 years.
Age-specific prevalence of prediabetes increased in 2016 in all populations (Fig. 1). The maximum change was in the town, followed by the PUV. The prevalence started increasing at 25–34 years in all populations. In the city, the increasing trend was present until 64 years; in the town, the maximum increase was at 45–54 years; and in the PUV, the rising trend continued until ≥64 years.
The percentage of smokers was 22.5, 15.0, and 36.8% and alcohol consumption was reported by 30.7, 25.3, and 39% in the city, town, and PUVs, respectively.
Conclusions
In this study from Tamil Nadu in India, we observed that in the last 10 years the prevalence of diabetes and prediabetes increased in all populations studied. However, varied secular changes were noted among the populations. Significant variations in the PR were observed. In the city, the prevalence (21.9%) had increased only by 8% in 10 years (PR 1.08 [95% CI 0.99–1.19], P = 0.08), indicating the probability of the prevalence reaching a plateau. The greatest increase, by 39%, was seen in the town (PR 1.39 [95% CI 1.25–1.55], P < 0.0001), and the prevalence (20.3%) was almost similar to that in the city. In the PUVs, a 34% increase in the prevalence in 10 years, from 9.2 to 13.4%, was seen in 2016 (PR 1.34 [95% CI 1.15–1.54], P < 0.0001). The prevalence was lower in the PUVs compared with the urban areas. It is important to note that a significant increase occurred among the populations in 2016, with much lower BMI and WC compared with the urban populations.
Prevalence of diabetes in urban India seems to be higher among the states that are economically stronger (1). At present, the total prevalence of diabetes is chiefly contributed by the urban population. Our study shows that the prevalence of diabetes in 2016 is 21.9% in the city of Chennai. The estimate was similar to that reported by the Center for cArdio-metabolic Risk Reduction in South Asia (CARRS) Study (26) in 2015, which was 22.8%. In that study, diagnosis of diabetes was made using the FPG and/or HbA1c values. The survey also reported that the prevalence in Delhi was 16.8%.
A series of studies in Tamil Nadu, in South India, had shown the changing scenario in lifestyle with associated increase in the prevalence of diabetes and prediabetes in the rural areas (6,10,27). Reports from other South Indian states, Pondicherry (28), Andhra Pradesh (29), Karnataka (30), and Kerala (31), had also highlighted the rising trend in diabetes in rural parts. Studies from other parts of India also showed a similar scenario (11). The recent INdia DIABetes (INDIAB) study (9) reported an increase in rural prevalence in many parts of India.
A systematic review and meta-analysis of the global rural diabetes prevalence covering 1990–2012 showed that in low- and middle-income countries, the rural prevalence had increased from 1.8 to 7.5% (32). High prevalence of diabetes in rural areas is reported from other South Asian countries, such as Bangladesh (11), Pakistan (33), and Sri Lanka (13).
In the earlier studies, diabetes prevalence in Chennai was higher in high socioeconomic stratum (4,10). Presently this association is not observed in the urban areas. In 2016, income did not show an association with diabetes prevalence in any of the three locations, whereas in the 2006 survey in the city, high income was associated with a higher prevalence of diabetes (10). The current study suggests a rising trend of diabetes in the low- and middle-income groups. The INDIAB study (9) also observed a transition occurring in India, with diabetes prevalence increasing in rural India as well as among the people of lower socioeconomic status living in the urban areas.
The strong associations of generalized obesity and abdominal obesity with diabetes are well-known. Interestingly, this study showed that abdominal obesity had a higher increase than generalized obesity in all the populations and its association with diabetes was more significant. BMI was found to be a determinant only in the city. The significance of the Asian phenotype of obesity, namely, high abdominal obesity with normal BMI, was seen in all the populations studied, particularly in the PUVs. The rise in WC was more predominant than that of BMI in the rural population (57.2%) than in the city (34%) and in the town (35.2%). Since nearly 60% of India’s population lives in rural areas, the major contribution to the future increase in the prevalence of diabetes will be by the rural population.
The present increase in prevalence of diabetes among Indian populations, especially in the rural areas, could have been partially influenced by the complex interaction between genes and environmental factors. Epigenetic changes occur during life and commonly in response to environmental stimuli, providing gene-environmental interaction (34), and are known to induce persistent changes in gene expression patterns and alter physiology (35). Intrauterine exposure to diabetes itself increases the risk of diabetes and obesity in the offspring (36). Several animal and human studies suggest that DNA methylation, a class of epigenetic modification, is sensitive to early environmental factors including maternal diet (37). It has been shown that Indian newborn babies, especially among the rural population, have low birth weight and lean body mass but higher subcutaneous adiposity (38), which is related to an intrauterine imbalance between two related micronutrients (vitamin B12 and folate) (39). Further, B12 supplementation has been shown to influence methylation of genes associated with diabetes and its intermediate traits. Overall, there is strong evidence that in addition to the “thrifty phenotype” genes, a combination of antenatal, epigenetic, and postnatal influences renders a high risk of diabetes in the South Asian populations (40).
The mean age of newly diagnosed diabetes was similar in the city and the PUVs, but it showed an increase in the town in 2016. No reduction in the age of people with new diabetes was seen. Self-reported diabetes had increased in all populations, indicating better awareness and improved health seeking behavior, even among the rural population.
In 2006 a significant reduction in IGT was noted in Chennai and the PUVs compared with corresponding prevalence in previous years (10). During the same period, Mohan et al. (5) observed a rapid conversion of IGT to diabetes in Chennai with a reduction in the prevalence of IGT. In the current study, the prevalence of prediabetes increased in urban and rural populations.
Presently, the rise in prevalence of diabetes in India (2) and other South Asian countries (14,15) is contributed chiefly by the urban population. Rural to urban migration was shown to be associated with increase in obesity and diabetes (12). The recent studies in India (9,10,14,15) and other South Asian countries (7) indicate that adoption of urban lifestyle within the rural habitats produces similar detrimental changes resulting in health hazards. The trend is likely to change, since many of the Asian countries show rising trends in the prevalence of diabetes in the rural populations as well. Our present study has demonstrated definite evidence of this phenomenon. The population living in rural areas is higher in India. In the future the major contribution to population with diabetes is likely to be due to the rising rural prevalence of diabetes. Although we have demonstrated this phenomenon in the southern state of Tamil Nadu, it is possible that a similar phenomenon will be occurring in rest of the country.
The current study has the merit that the survey was done in the same three populations of varied urbanization, by virtue of which the transitions occurring in a decade in the pattern of diabetes and its risk associations could be assessed. Moreover, the good response rates of participants in all the study locations give strength to our findings.
The limitations of our study included subjective assessments of dietary habits and physical activity. We have used similar quantification of dietary habits and statistical methodology in both the surveys so that the comparative assessments are more valid. Objective measurements of physical activity are time-consuming and may not be acceptable to the large number of people surveyed. Plasma glucose levels were measured using a glucometer to avoid the errors that could occur if the blood samples had to be transported from distant places to the laboratory. These readings correlated well with those obtained in the laboratory. A correlation factor has also been applied.
In summary, the findings of the 2016 survey highlight a transition in the secular trends of diabetes and prediabetes in different populations studied in Tamil Nadu. In a decade of time, the prevalence of diabetes in the city showed only a nonsignificant increase, whereas a higher percentage increase occurred in the town and PUVs. The prevalence in the town is closer to that in the city.
Recent regional and national studies (9,10,27–31) show that the prevalence of diabetes is increasing in all regions of India including the rural areas. There have been attempts by governmental and nongovernmental agencies to promote preventive endeavors at the national level. However, more urgent focus on prevention is needed for India and also other developing nations to curb the diabetes epidemic. National policies have to be enforced in all countries.
Clinical trial reg. no. NCT03490136, clinicaltrials.gov
Article Information
Acknowledgments. The authors acknowledge the help rendered by the epidemiology and laboratory team of the IDRF (C.K. Sathish Kumar, M. Karthikeyan, A. Catherin Seeli, R. Sekar, and G. Sridhar). The authors thank D.T. Arasu and his team, of I5 Clinical Research Private Limited, Chennai, India, for helping with the epidemiological screening and data entry. The authors thank Stephen Sharp, MRC Epidemiology Unit, Cambridge, U.K., for advice on statistical analyses. The authors thank the participants of the epidemiology survey.
Funding. The study was funded by the IDRF.
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. A.N., C.S., and A. Ram. developed the study plan, designed and supervised the conduct of the study, and prepared the manuscript. A.N., C.S., and A. Ram. analyzed data and interpreted the results. K.S., P.S., and M.S. participated in fieldwork, study coordination, and data collection. L.V. helped with data management and statistical analyses. A. Rag. and R.V. contributed to discussion and reviewed the manuscript. All authors read and approved the final manuscript. A. Ram. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.