Spatial Autocorrelation in Health Loss of Kashin-Beck Disease and Relationship to Environmental Factors: A Cross-Sectional Study in Bin County, Shaanxi Province, China

doi:10.21203/rs.3.rs-86107/v1

Download PDF

Research

Spatial Autocorrelation in Health Loss of Kashin-Beck Disease and Relationship to Environmental Factors: A Cross-Sectional Study in Bin County, Shaanxi Province, China

https://doi.org/10.21203/rs.3.rs-86107/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Background: Kashin-Beck disease (KBD) is one of the major endemic diseases in China, which severely impacts the physical health and life quality of people. A better understanding of the spatial distribution of the health loss from KBD and its influencing factors will help to identify areas and populations at high risk so as to plan for targeted interventions.

Methods: The data of patients with KBD at village-level were collected to estimate and analyze the spatial pattern of health loss from KBD in Bin County, Shaanxi Province. The years lived with disability (YLDs) index was applied as a measure of health loss from KBD. Spatial autocorrelation methodologies, including Global Moran’s I and Local Moran’s I, were used to describe and map spatial clusters of the health loss from KBD. In addition, selenium concentrations in soil and wheat samples in Bin County were determined to detect their relationships with the distribution of health loss of KBD.

Results: The estimation of YLDs for KBD showed that patients with KBD of grade II and patients over 50 years old contributed the most to the health loss in Bin County. There was no significant difference between the two genders. The spatial patterns of YLDs and YLD rate of KBD were clustered significantly at both global and local scales. Villages in the southwestern and eastern regions revealed higher health loss, while those in the northern regions exhibited lower health loss. This clustering was found to be closely related to organically bound Se in soil and poverty rate of KBD patients.

Conclusions: Our results suggest that future treatment and prevention of KBD should focus on endemic areas with high organically bound Se in soil and poor economic conditions. The method of estimating the health loss of KBD with YLDs can be useful for KBD surveillance for public health officials.

Geographic Information Systems

Kashin-Beck disease (KBD)

health loss

spatial autocorrelation

selenium

Kashin-Beck disease (KBD) is a chronic, endemic, deformative osteoarthropathy, which is known for the formation of multi-joint hyperplasia bone changes [1]. The disease usually starts in childhood and attacks the growth of joint cartilage [2]. Patients with mild KBD have symptoms of joint thickening and deformation, muscle atrophy, and often accompanied by pain; while those with serious KBD manifest developmental disorders, short limbs and malformation, loss of labor capacity, and confined self-care ability [3]. KBD has been discovered since the sixteenth century, and it is distributed diagonally from northeastern China to Tibet in the southwest, with additional endemic regions in Siberia and North Korea [4]. In China, KBD has been effectively controlled and even eliminated in most affected areas, but the number of existing patients is still huge. According to the 2018 health statistics issued by the Chinese Ministry of Health [5], there were currently 177,018 individuals affected by KBD in 379 counties of 13 provinces or autonomous regions, and the largest number was found in Shaanxi Province (60,157 individuals), accounting for 34.0% of the total existing patients with KBD.

Although the prevalence of KBD has reached the level of control in Shaanxi Province [6], the harm of KBD to human health still exists. For a long time, the severity of KBD has been evaluated by the prevalence rate or X-ray positive rate of children, which can only count the number of cases and cannot value the harm of the non-fatal disability caused by different degrees of KBD to population. The disability adjusted life year (DALY) is a time-based metric, proposed by the World Bank and World Health Organization (WHO) to estimate the burden of disease [7]. DALY is the sum of years of life lost due to premature mortality (YLLs) and years lived with disability (YLDs), which quantify the health loss of both fatal and non-fatal consequences [8]. As far as the non-fatal KBD is concerned, its loss of healthy life years is mainly caused by YLDs, and there is no YLLs. DALYs have been widely used to reflect health gaps of diseases among age-sex groups and regions [9–12]. However, to the best of my knowledge, there have been no relevant research on KBD. Given its serious impacts on the health and life of local residents, an estimation of the health loss is quite necessary. Furthermore, a better understanding of the spatial distribution of the health loss of KBD will be helpful to identify areas and populations at high risk so as to plan for targeted interventions.

The health loss of KBD is directly related to the severity of KBD. A great number of studies have reported that the prevalence of KBD is closely associated with the low-selenium (Se) environment [13–15]. Se in soil and crops in KBD endemic areas are generally lower than those in non-KBD endemic areas. Nevertheless, little is known about the relationship between the health loss of KBD and Se in the environment. The exploration of their relationships may provide a new perspective for the etiologic research of KBD. In addition, socio-economic factors such as educational attainment [16, 17] and family income [18, 19] can also affect the health loss from diseases in varying degrees. Although a few studies have investigated the key factors (i.e., age, educational attainment, severity of KBD, economic level, etc.) influencing the health-related quality of life of adult patients with KBD in Shaanxi Province [20, 21], the knowledge of factors contributing to the health loss of KBD is still limited.

Bin County in Shaanxi Province, as one of the most serious KBD endemic areas, was chosen as the study area. In this study, basic conditions of KBD areas and patients in Bin County were investigated at the village-scale. Corresponding environmental samples including cultivated topsoil and wheat were collected to analyze concentrations of total Se and soil Se speciation. The main objectives were to 1) quantify the health loss of KBD based on the YLDs metric; 2) analyze the spatial distribution of YLDs and YLD rate; and 3) explore factors influencing the health loss of KBD from natural and social environment. The results will provide theoretical basis for assisting public health officers to optimize the allocation of health resources and to prevent and control endemic diseases.

Study area

Bin County is located in the central and western part of Shaanxi Province and the middle of the natural Se deficiency belt in China (Fig. 1). It is a national surveillance site for KBD. KBD in Bin County was ever very serious, whose prevalence rate had been the highest in China from 1992 to 1995 [22]. With the measures of selenium supplement, water source changes and economic development, the prevalence of KBD was under control basically. In 2007, the X-ray detection rate in Bin County was only 0.43%, which had reached the national standard [6]. There have been no new cases discovered in children in recent years, but KBD in adults is still very serious.

Data sources

The epidemiological data used for the calculation of YLDs were obtained from the general survey of KBD in Bin County, which was carried out by the Binxian Center for Disease Prevention and Control in 2018. This survey focused on the prevalence status of KBD at each village in 16 townships of Bin County, including the number of existing patients, the degree of KBD. The detailed individual information (such as age, gender, education level, occupation, place of residence, medical history, etc.) for each patient with KBD was recorded. As part of the survey, basic demographic and socioeconomic information were collected at each village as well.

Sampling and analyses

Se contents in the environment in Bin County were analyzed to explore factors influencing the health loss of KBD. According to the distribution of KBD endemic areas in Bin County, 84 villages were selected to collect cultivated topsoil samples (0–20 cm). The location of sampling sites is shown in Fig. 1. Three sub-samples were taken from the farmland with relatively uniform distance in each village, which were then mixed together to create a composite sample. At each sampling site, wheat samples were randomly collected from a household. The determination of total Se referred to national standards on the determination of Se in foods (GB/T 5009.93–2003) and soil (NY/T 1104–2006). Furthermore, three different fractions of soil Se (water-soluble, exchangeable, and alkali-soluble organically bound) were extracted successively. The details of the method were represented in our previously published paper [23]. Se in soil and wheat samples were determined by hydrogen generation-atomic fluorescence spectrometry (AFS-9780, HaiGuang Instruments, Beijing, China), whose detection limit is 0.02 ng/ml and RSD is less than 1.0%. Reagent blanks, duplicated samples, and national standard reference materials (GBW10011 for wheat and GBW07410 for Tibetan soil) were used for analytical quality control.

Calculation of YLDs and YLD rate

The YLDs can be calculated from either an incidence perspective or a prevalence perspective [9]. The former is the product of incidence, disability weights and average duration of disease; the latter is the product of prevalence of disease and disability weights, which is convenient for comparison with the recent GBD studies [24]. In this study, the prevalence-based YLDs were calculated for the analysis. Discounting and age weighting were not applied. The basic formula is as follows:

where P is the number of prevalent cases of KBD; DW is the disability weight of different degrees of KBD. DW is usually estimated based on evaluation scales [9, 25] or referenced from the results of Global Burden of Disease (GBD) studies [12, 26]. However, no paper has previously reported disability weights of KBD. Considering that KBD has many similarities with rheumatoid arthritis in clinical manifestations [27], disability weights of rheumatoid arthritis in the latest GBD 2017 study were directly adopted in this study. Thus, disability weights of KBD in different grades were assigned as 0.117 (grade I), 0.317 (grade II) and 0.581 (grade III) in the light of the sequela of rheumatoid arthritis [28].

The YLD rate (YLDs per 1000 population) is calculated from the YLDs in different cohorts divided by the total target population and then multiplied by 1000.

Spatial autocorrelation analysis

Clustering in the health loss from KBD was analyzed using both global and local spatial autocorrelation statistics. First, the global Moran’s I test statistic was computed to test the null hypothesis of no significant clustering of YLDs and YLD rate in the entire study region [29]. The values of Moran’s I range from − 1 (dispersed) to 1 (clustered). The threshold value of Moran’s I index is 0, indicating complete spatial randomness. The statistical significance for the spatial autocorrelation relationship is determined by standardizing the statistic Z value [30]. At a confidence level of 0.05, |Z|=1.96; at a confidence level of 0.01, |Z|=2.58. The significant or highly significant level was set when |Z|>1.96 or |Z|>2.58.

Second, Anselin Local Moran’s I statistic was applied to examine disease spatial autocorrelation at the local level. Unlike the global Moran’s I, the expected value of local Moran’s I varies for each sampling village because it is calculated in relation to its particular set of neighbours [31]. The local Moran’s I index identified locations of clusters or hotspots where the value of the index was extremely pronounced across localities, as well as those of spatial outliers [32]. The significance of the local Moran’s I was calculated using a randomization test on the Z-score value [33]. A positive Z-score value indicates that the health loss of KBD in one village is surrounded by similar health loss in neighboring locations (high-high or low-low), thus forming a spatial cluster. A negative Z-score value indicates that the high health loss of KBD in one village is surrounded by low neighbors (high-low) and vice versa (low-high). Similarly, the significance level was set when |Z|>1.96. |Z|≤1.96 indicates presence of a random distribution. The results were mapped to display the specific locations of clusters (high-high and low-low) and potential outliers (high-low and low-high).

Data processing and chart production were mainly done using SPSS 23.0, ArcGIS 10.5, and Origin 8.0. Correlation analysis and T test were performed using SPSS 23.0. Spatial autocorrelation analysis was conducted in ArcGIS 10.5 using spatial statistics tools.

Study population

296,770 of people (17.3% were children aged under 12 years) in 244 villages were examined in Bin County in 2018. The gender ratio of the investigated population was 1.12:1 (male to female). Overall, 1.34% were reported suffering from KBD, with no new cases discovered in children. The age of the patients ranged from 19 to 97 years old. Amongst patients with KBD, 52.9% were males and 47.1% were females. Patients with KBD of grade I, II and III accounted for 57.4%, 37.4% and 5.3% of the total amount, respectively. 44.6% of them had been taking medical treatment for a long-term. Farmers were the main occupation for KBD patients. 39.9% never received education, and 19.2% were from poor poverty-stricken households whose annual net income per capita was lower than 2950 RMB.

Age-sex distribution of YLDs for KBD

Table 1 showed the calculated YLDs in different KBD grades by gender. The total health loss from KBD in Bin County in 2018 was estimated at 858.78 YLDs (2.89 YLD per 1000 population, 53.8% for males and 46.2% for females). YLDs for males were higher in all grades of KBD than those for females, but with no statistical significance (p > 0.05). Among different KBD grades, KBD of grade II contributed most to the YLDs, followed by KBD of grade I, accounting for 54.4% and 31.0% of the total YLDs respectively. The same trends were observed for males and females. When compared with the prevalence rate of KBD, it was found that there was no consistent corresponding relationship. Although the largest contribution of YLDs was KBD of grade II (1.58 YLD per 1000 population), the highest prevalence rate was observed in KBD of grade I (0.77%).

Table 1

Years lived with disability (YLDs) for KBD in Bin County (2018)
KBD grades	Male	Female	Total		Prevalence (%)
KBD grades	YLDs	YLDs	YLDs	YLD/1000	Prevalence (%)
I	136.42	129.87	266.29	0.90	0.77
II	258.99	208.59	467.58	1.58	0.50
III	66.82	58.10	124.92	0.42	0.07
Total	462.23	396.56	858.78	2.89	1.34

Figure 2 displayed the age-specific distribution of YLDs between genders and KBD grades in Bin County. According to the distribution of the number of patients, age was divided into five groups: 19~, 40~, 50~, 60~, 70 years and over. The age group of 0 ~ 18 was not included in Fig. 2 as there were no cases under the age of 18 years. The results showed that YLDs increased with age and decreased at the highest 50 ~ age group (330.49 YLDs, 38.5% of total YLDs). Approximately 85.6% of the total YLDs were from the age group of 50 years and above. YLDs in KBD of grade II were the highest in all age groups, followed by KBD of grade I. There were no significant differences in YLDs among KBD grades (p > 0.05). These trends were generally consistent in the two genders, only with relatively high YLDs for males in all age groups.

Spatial distribution of health loss from KBD

Figure 3 presented the spatial distribution of healthy life years loss from KBD in Bin County at the village-scale. Obvious regional variations were observed in both YLDs and YLD rate, with higher values in the southern regions and lower values in the northern regions of Bin County. Areas with the most serious health loss were mainly clustered in the southwest. The loss of healthy life years from KBD in villages of Bin County ranged from 0 to 25.50 YLDs, while their YLD rates had a larger variation, ranging from 0 to 40.47 YLD/1000. The top five townships with the highest YLD rates were Hanjia (19.83 YLD/1000, 97.16 YLDs), Xiangmiao (6.04 YLD/1000, 78.49 YLDs), Tandian (5.74 YLD/1000, 88.97 YLDs), Xinbaozi (5.54 YLD/1000, 73.43 YLDs), and Chejiazhuang (4.95 YLD/1000, 69.39 YLDs).

Spatial autocorrelation and hotspot detection

To detect the spatial autocorrelations of YLDs and YLD rate, global Moran’s I statistics were calculated at village-level in Bin County. It was found that the Moran’s I index values for YLDs and YLD rate were 0.27 and 0.37, indicating that the distribution of health loss from KBD have positive correlations. The Z score, which is a standardized statistic, for YLDs and YLD rate were 7.607 (p < 0.0001) and 10.511 (p < 0.0001), respectively. Apparently, the global spatial autocorrelations of YLDs and YLD rate presented significant spatial clustering. Villages with KBD cases were not randomly distributed in space among all villages in Bin County. Similarly, the YLDs and YLD rate of different gender groups also showed significant spatial clustering at village-level. The Moran’s I index values of YLDs and YLD/1000 for males were 0.15 and 0.21 (p < 0.01), those for females were 0.12 and 0.17 (p < 0.01), respectively.

To further identify the locations of significant clusters (hot and cold spots), Anselin local Moran’s I index was applied at village-level in Bin County. Maps in Fig. 4 showed the locations with significant local Moran’s I statistics and classified those locations by type of association (LISA cluster map). The high-high and low-low clusters are indications of high values surrounded by high values and low values surrounded by low values. In contrast, the high-low and low-high locations are indications of spatial outliers. There were some outstanding spatial clusters of both YLDs and YLD/1000 observed in Bin County. The clustered villages with high YLDs and YLD/1000 (hotspots) were found to cover most areas of Hanjia and Chejiazhuang Townships in the southwest and parts of Tandian and Xiangmiao Townships in the east of Bin County. The clustered villages with low YLDs and YLD/1000 (cold spots) were mainly found in the north and central of Bin County, covering Xipo, Yongle, Xinmin, Xiaozhang and Chengguan Townships. The hotspots of YLDs were slightly different from that of YLD/1000. Several clustered villages with high YLDs were also observed in Xinbaozi and Longgao Townships, discontinuously distributed in the southeast of Bin County.

Effect of environmental factors on health loss

We computed Spearman’s correlation coefficients and its 95% confidence intervals to assess the relationships between the health loss and variables of interest (Table 2). Natural factors related to environmental Se contents and social factors including the poverty rate and educational attainment were selected for the correlation analysis. The results showed that YLDs and YLD rate of KBD had no significant correlations with the total Se contents in soil and wheat, but were positively and significantly correlated with organically bound Se in soil (r_YLDs = 0.216, r_YLD/1000 = 0.217, p < 0.05, N = 83). By contrast, the prevalence of KBD had no significant correlations with any environmental Se factors. As for social environmental factors, only the poverty rate of patients with KBD showed significantly positive correlations with YLD rate (r = 0.267, p < 0.05, N = 83) and prevalence rate (r = 0.264, p < 0.05, N = 83).

Table 2

Spearman’s correlation coefficients for YLDs, YLD rate, prevalence and environmental factors
	Natural environmental factors				Socio-economic factors
	Soil Se	Wheat Se	Available Se^a	Organic Se	Poverty	Education
YLDs	0.153	-0.008	-0.052	0.216*	0.203	0.223
YLD/1000	0.150	0.064	-0.070	0.217*	0.267*	0.220
Prevalence	0.140	0.089	-0.060	0.213	0.264*	0.213
^a Available Se represents water-soluble and exchangeable fractions of Se in soil.
* indicates significant (p < 0.05).

One of the relevant contributions of using YLDs is that they give disability weights in different degrees according to the severity of KBD, while simultaneously taking its prevalence and severity into account. Consequently, the ranking of the health loss caused by KBD differs from the ranking based on prevalence rates (Table 1). Compared with traditional index (e.g. prevalence), YLDs is more conductive to the measurement of these non-fatal consequences and the exploration of relationships between KBD and low-Se environment.

The study showed that spatial distribution patterns of YLDs and YLD rate of KBD were significantly clustered, and identified their hotspots and cold spots in Bin County. Factors that might affect this spatial clustering were analyzed from natural and social environment. In the aspect of natural factors, correlations of YLDs and YLD rate with environmental Se were analyzed. The results found that organically bound Se in soil significantly and positively influenced the distribution of YLDs and YLD rate. Organically bound Se is known as the unavailable fraction in soil, which is mainly found in the fulvic acid and humic acid of soil humus [23]. Previous studies have revealed that the organic matter (mainly fulvic acid) in drinking water may be an etiological factor of KBD [34]. They found that the total amount of organic matter and humic acid in drinking water from KBD endemic areas were significantly higher than those in non-endemic areas, and changes in water sources effectively prevented the occurrence of KBD [35, 36]. It is well-known that soluble or suspended compounds (including organic matter) in soil can infiltrate into water through leaching, finally leading to accumulation in water. Consequently, the distribution of the severity of KBD may be directly associated with that of organic matters (mainly fuvic acid) in soil. The non-significant relationship between prevalence rate and organically bound Se is possibly due to the limitation and one-sidedness of this index, which only counts the number of cases and fails to reflect the harm of different KBD grades comprehensively. With regard to total Se contents in soil and wheat, their non-significant correlations with both YLD rate and prevalence rate indicate that the dependence of local residents on low-Se environment is weakened. With the improvement of living standards and measures of Se supplement, the food sources of Se intake increase. Residents do not have to entirely rely on local crops.

In aspect of social factors, positive and significant correlations were observed for YLD rate and prevalence with poverty rate. This is consistent with previous findings of the impact of economic levels on life quality of patients with KBD [21]. KBD tends to occur in remote areas with weak infrastructure and inconvenient transportation. These regions are usually poverty-stricken areas. Poor economic conditions restrict the improvement of living and nutritional levels of local residents, thus increasing the risk of occurrence and development of KBD. Moreover, positive but non-significant correlations were observed for the indexes with educational attainment, indicating that people with higher education levels suffered more health loss. This result does not accord with previous findings. According to Chen et al. [20], those with higher education levels tended to easily acquire and accept health knowledge of KBD, thus early diagnosing and treating the disease. By comparison, the positive correlation of this study may be attributed to the generally poor education attainment of patients in Bin County. Among the educated population, 72.0% only received primary education, which may not be substantively helpful with the prevention of KBD.

Some limitations should be noted. First, reliable sources of disability weight are required to calculate YLDs. Disability weights of different sequela of KBD in the present study were based on the results of rheumatoid arthritis in GBD 2017. It was assumed that these data would be acceptable as the two diseases have similar symptoms according to the health state lay descriptions in GBD 2017 disability weights dataset [28]. Second, discounting and age weighting were not considered when calculating YLDs. Given the long duration, non-fatal outcomes and late-onset serious symptoms of KBD, which is quite similar to other endemic diseases such as schistosomiasis and goiter, we adopted the same strategy carried out in the GBD study of these two diseases [7, 37]. Namely, the duration of the disease was hypothesized as one year. Consequently, there is no discounting problem. With regard to age weights, the original GBD study weighted a healthy life lived at very young and old ages lower than other ages [38]. However, this is still controversial, since not all such studies agree that the youngest and oldest ages should be given less weight, nor do they agree on the relative magnitude of the differences [10, 39–40]. Thus, the social values were not considered in this study. Third, given that the income variable of each patient or family is not available in the dataset, we instead used the poverty rate of each village, which failed to distinguish the effect of specific economic levels. Finally, the effect of population age structure on the distribution of health loss of KBD was not considered in the analysis mainly due to the insufficiency of these data in the study area. It should be noted that 85.6% of the total YLDs in Bin County were from the age group of 50 years and above (Fig. 2), indicating that a higher proportion of aging population may lead to greater loss of healthy life years.

The present study found that the health loss from KBD in Bin County was significantly clustered in spatial distribution. This clustering was collectively affected by environmental Se and socio-economic factors. In particular, organically bound Se in soil and poverty rate might be the most influential natural and social factors. Our results suggest that future treatment and prevention of KBD should focus on endemic areas with high organically bound Se in soil and poor economic conditions. Furthermore, this study exhibited the method of estimating the health loss of KBD with YLDs, which can be useful for KBD surveillance for public health officials.

KBD: Kashin-Beck disease; DALY: disability adjusted life year; YLDs: years lived with disability.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

The datasets collected and analyzed during the current study are not publicly available due to confidentiality requirements, but are available from the corresponding author on reasonable request and with permission of the Binxian Center for Disease Prevention and Control.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was funded by the National Natural Science Foundation of China [Project Number: 41907390 and 41671500].

Author’s contributions

JW performed data analysis and developed the initial draft of the manuscript. XYW carried out the field sampling and chemical analysis. HRL contributed to revising the paper. LSY conceived the presented idea. YCL and CK helped with the field investigation and sampling. All authors read and approved the final manuscript.

Acknowledgements

Special thanks to the Binxian Center for Disease Prevention and Control for providing the data of patients with KBD.

Kraus VB. 185 - Rare osteoarthritis: Ochronosis and Kashin-Beck disease. Rheumatology. 2015;2:1536-47.
Farooq U, Guo X, lrshad R, Yaqoob S. Pattern of joints involvement in Kashin-Beck disease: a local osteochondropathy in China. J Ayub Med Coll Abbottabad. 2010;22(1):97-100.
Guo X, Ma WJ, Zhang F, Ren FL, Qu CJ, Lammi MJ. Recent advances in the research of an endemic osteochondropathy in China: Kashin-Beck disease. Osteoarthr and Cartilage. 2014;22:1774-83.
Lian W, Liu H, Song Q, Liu YQ, Sun LY, Deng Q, et al. Prevalence of hand osteoarthritis and knee osteoarthritis in Kashin-Beck disease endemic areas and non Kashin-Beck disease endemic areas: A status survey. PLoS One. 2018;13(1): e0190505.
National Health Commission. China Health Statistics Yearbook 2019. Beijing: Beijing Union Medical University Press; 2019.
Lv XY, Xu GY, Chen P. Report on monitoring results of Kaschin-Beck disease in Shaanxi Province in 2007. Chin J Ctrl Endem Dis. 2009;24(4):309-10.
Mathers CD, Vos T, Lopez AD, Salomon J, Ezzati M. National Burden of Disease Studies: A Practical Guide. 2nd ed. Global Program on Evidence for Health Policy. Geneva: World Health Organization; 2001.
Murray CJL, Vos T, Lozano R, Naghavi M, Flaxman AD, Michaud C, et al. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2197-223.
Maredza M, Bertram MY, Tollman SM. Disease burden of stroke in rural South Africa: an estimate of incidence, mortality and disability adjusted life years. BMC Neurol. 2015;15:54.
Dodhia H, Phillips K. Measuring burden of disease in two inner London boroughs using disability adjusted life years. J Public Health. 2008;30(3):313-21.
Gènova-Maleras R, Álvarez-Martín E, Morant-Ginestar C, Larrea-Baz N, Catalá-López F. Measuring the burden of disease and injury in Spain using disability-adjusted life years: an updated and policy-oriented overview. Public Health. 2012;126:1024-31.
Yang LP, Liang SY, Wang XJ, Li XJ, Wu YL, Ma W. Burden of disease measured by disability-adjusted life years and a disease forecasting time series model of scrub typhus in Laiwu, China. PLoS Negl Trop Dis. 2015;9(1):e3420.
Li SJ, Li W, Hu X, Yang LS, Xirao RD. Soil selenium concentration and Kashin-Beck disease prevalence in Tibet, China. Front Environ Sci En. 2009;3:62-8.
Zhang BJ, Yang LS, Wang WY, Li YH. Environmental selenium in the Kaschin–Beck disease area, Tibetan Plateau, China. Environ Geochem Health. 2011;33(5):495-501.
Wang J, Li HR, Yang LS, Li YH, Wei BG, Yu JP, Feng FJ. Distribution and translocation of selenium from soil to highland barley in the Tibetan Plateau Kashin-Beck disease area. Environ Geochem Health. 2017;39:221-9.
Agardh EE, Sidorchuk A, Hallqvist J, Ljung R, Peterson S, Moradi T, Allebeck P. Burden of type 2 diabetes attributed to lower educational levels in Sweden. Popul Health Metr. 2011;9:60.
Zadro JR, Shirley D, Pinheiro MB, Sánchez-Romera JF, Pérez-Riquelme F, Ordoñana JR, Ferreira PH. Does educational attainment increase the risk of low back pain when genetics are considered? A population-based study of Spanish twins. Spine J. 2017;17(4):518-30.
Jia TW, Sun LP, Hong QB, Deng Y, Zhang GH, Yi P, Guo JG, Zhou XN. Burden of disease in schistosomiasis japonica IV Community socioeconomic determinants of chronic schistosomiasis. Chin J Schisto Control. 2012;24(1):10-3.
Williams JS, Ng N, Peltzer K, Yawson A, Biritwum R, Maximova T, Wu F, Arokiasamy P, Kowal P, Chatterji S. Risk factors and disability associated with low back pain in older adults in low- and middle-income countries. Results from the WHO study on Global Ageing and Adult Health (SAGE). PLoS One. 2015;10(6):e0127880.
Chen P, Dai HX, Shi Z, Chen G, Feng QH, Yang PR. An investigation and analysis of health-related quality of life in patients with Kashin-Beck disease in Shaanxi Province. Chin J Endemiol. 2018;37(9):754-9.
Li XX, Fang H, Yuan HJ, Yan H, Xue SW. Relationship between income level and quality of life in patients with Kashin-Beck disease. Chin J Endemiol. 2020;39(4):285-8.
Wang ZL, Chen JH, Xiong YM, Liu WJ, Tan XW, Yang ZT, Xue L, Geng D. The prevalent growth and decline studies on results of national Kashin-Beck disease surveillance in Shaanxi province in 1992 to 2002. Chin J Ctrl Endem Dis. 2003;18(3):137-9.
Wang J, Li HR, Li YH, Yu JP, Yang LS, Feng FJ, Chen Z. Speciation, distribution and bioavailability of soil selenium in the Tibetan Plateau Kashin–Beck disease area—A case study in Songpan County, Sichuan Province, China. Biol Trace Elem Res. 2013;156:367–75.
Vos T, Flaxman AD, Naghavi M, Lozano R, Michaud C, Ezzati M, et al. Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2163-96.
Jia TW, Zhou XN, Wang XH, Utzinger J, Steinmann P, Wu XH. Assessment of the age-specific disability weight of chronic schistosomiasis japonica. Bull World Health Organ. 2007;85:458-65.
Cai L, Bi WH, Wan CH, Huang WX, Li JM, Yang R. Study on burden of stroke in Shilin county of Kunming. Chin J Public Health. 2007;23(1):45-6.
Liu D, Wang ZL. Differential diagnosis of Kashin-Beck disease and rheumatoid arthritis. Chin J Ctrl Endem Dis. 2009;24:195-7.
Global Burden of Disease Collaborative Network. Global Burden of Disease Study 2017 (GBD 2017) Disability Weights. Seattle, United States: Institute for Health Metrics and Evaluation (IHME); 2018. http://ghdx.healthdata.org/record/ihme-data/gbd-2017-disability-weights. Accessed 16 Sep 2020.
Tsai PJ, Lin ML, Chu CM, Perng CH. Spatial autocorrelation analysis of health care hotspots in Taiwan in 2006. BMC Public Health. 2009;9:464.
Gruebner O, Khan MM, Lautenbach S, Müller D, Kraemer A, Lakes T, Hostert P. A spatial epidemiological analysis of self-rated mental health in the slums of Dhaka. Int J Health Geogr. 2011;10:36.
Chaikaew N, Tripathi NK, Souris M. Exploring spatial patterns and hotspots of diarrhea in Chiang Mai, Thailand. Int J Health Geogr. 2009;8:36.
Anselin L. Local indicators of spatial associations-LISA. Geogr Anal. 1995;27:93-115.
Yeshiwondim AK, Gopal S, Hailemariam AT, Dengela DO, Patel HP. Spatial analysis of malaria incidence at the village level in areas with unstable transmission in Ethiopia. Int J Health Geogr. 2009;8:5.
Peng A, Wang WH, Wang CX, Wang ZJ, Rui HF, Wang WZ, Yang ZW. The role of humic substances in drinking water in Kashin-Beck disease in China. Environ Health Persp. 1999;107(4):293-6.
Guo X. Progression and prospect of etiology and pathogenesis of Kashin-Beck disease. Journal of Xi’an Jiaotong University (Medical Sciences). 2008;29(5):481-8.
Hou X, Zhai SS, Han JY, Xu W, Kong XZ. A retrospective survey on the effect of Kashin-Beck disease prophylaxis by a change of water sources. J Environ Sci (China). 1991;3(4):67-71.
Jia TW, Sun LP, Hong QB, Deng Y, Zhang GH, Wang H, Yi P, Guo JG, Zhou XN. Burden of disease in schistosomiasis japonica I Calculation and evaluation of years lived with disability of chronic schistosomiasis. Chin J Schisto Control. 2011;23(3): 243-8.
Murray CJL, Lopez AD, World Health Organization, World Band, Harvard School of Public Health. The Global Burden of Disease: a comprehensive assessment of mortality and disability from diseases, injuries and risk factors in 1990 and projected to 2020. Boston: Harvard University Press. 1996. https://extranet.who.int/iris/restricted/handle/10665/41864. Accessed 16 Sep 2020.
Kominski GF, Simon PA, Ho A, Luck J, Lim YW, Fielding JE. Assessing the burden of disease and injury in Los Angeles County using disability-adjusted life years. Public Health Rep. 2002;117:185-91.
Anand S, Hanson K. Disability-adjusted life years: a critical review. J Health Econ. 1997;16:685-702.

Download PDF

Version 1

posted

You are reading this latest preprint version

Spatial Autocorrelation in Health Loss of Kashin-Beck Disease and Relationship to Environmental Factors: A Cross-Sectional Study in Bin County, Shaanxi Province, China

Status:

Version 1

Abstract

Figures

Background

Methods

Study area

Data sources

Sampling and analyses

Calculation of YLDs and YLD rate

Spatial autocorrelation analysis

Results

Study population

Age-sex distribution of YLDs for KBD

Spatial distribution of health loss from KBD

Spatial autocorrelation and hotspot detection

Effect of environmental factors on health loss

Discussion

Conclusions

Abbreviations

Declarations

References

Status:

Version 1