key: cord-0804107-bgkusmi7 authors: Huang, Li-Juan; Zha, Jun-Jing; Cao, Nv-Wei; Zhou, Hao-Yue; Chu, Xiu-Jie; Wang, Hua; Li, Xian-Bao; Li, Bao-Zhu title: Temperature might increase the hospital admission risk for rheumatoid arthritis patients in Anqing, China: a time-series study date: 2021-10-31 journal: Int J Biometeorol DOI: 10.1007/s00484-021-02207-9 sha: abdb991d46f3fa1af696fd470eb930030aa575de doc_id: 804107 cord_uid: bgkusmi7 Temperature has been studied in relation to many health outcomes. However, few studies have explored its effect on the risk of hospital admission for rheumatoid arthritis (RA). A distributed lag non-linear model (DLNM) was used to analyze associations between mean temperature, diurnal temperature range (DTR), temperature change between neighboring days (TCN), and daily admissions for RA from 2015 to 2019 in Anqing, China. Subgroup analyses based on age, gender, rheumatoid factors, and admission route were performed. In total, 1456 patients with RA were hospitalized. Regarding the cumulative-lag effects of extreme cold temperature (5th percentile = 3℃), the risks of admissions for RA were increased and highest at lag 0–11 (RR = 2.68, 95% CI: 1.23–5.86). Exposing to low (5th percentile = 1.9℃) and high (95th percentile = 14.2℃) DTRs both had increased risks of RA admission, with highest RRs of 1.40 (95% CI: 1.03–1.91) and 1.24 (95% CI: 1.0–1.53) at lag 0 day, respectively. As for TCN, the marginal risk of admission in RA patients was found when exposed to high TCN (95th percentile = 2.9℃) with the largest single-day effect at lag 10 (RR = 1.11, 95% CI: 1.01–1.23). In subgroup analyses, females were more susceptible to extreme cold temperature, low and high DTRs, and high TCN. In regard to extreme cold temperature, significant risk of hospital admission in females only appeared at lag 2 (RR = 1.48, 95% CI: 1.02–2.15) and lag 0–2 (RR = 2.35, 95% CI: 1.11–4.95). It is clear that RA patients exposed to changing temperature may increase risks of admission. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00484-021-02207-9. Rheumatoid arthritis (RA) is an autoimmune disease, which troubles patients with pain and swelling in multiple joints (Cross et al. 2014) . It affects 0.5-1% of the adult population (Scott et al. 2010) , and commonly occurred in people aged 40-70 (Lee and Weinblatt 2001) . Previous studies have demonstrated that prevalence and annual incidence of RA increased from 1990 to 2017 (Safiri et al. 2019) . Besides, the number of crude disability adjusted of life years of RA patients increased by 25% from 1990 (Moradi-Lakeh et al. 2017 ). However, the causes for RA progression remain unclear. It is generally believed that genetic factors (Karami et al. 2019; Constantin et al. 2001 ) and environmental factors, such as smoking (Liu et al. 2019b) , infection (Arleevskaya et al. 2018) , and microbiota (Zhang et al. 2015) , might affect the progression. Previous studies have shown that air pressure and seasonal changes have harmful effects on disease activity and joints of RA (Terao et al. 2014; Mori et al. 2019) . Thus, the fluctuant meteorological factors have attracted much attention. Temperature, as a fluctuant environmental factor, has been studied in relation to many health outcomes (Malig et al. 2019; Schreier et al. 2013; Wang et al. 2021) . Evidence have shown that RA was particularly sensitive to weather fluctuations, including altering temperature. Rheumatic disorders were reported to be correlated with the mean temperature (Patberg and Rasker 2004; Vergés et al. 2004) . Joint pain is a common clinical manifestation in patients with RA, and the relationship between joint pain and temperature has been reported in previous studies. As early as the 1990s, Daniela and Abraham found that a change in the weather could exacerbate the symptoms of arthritis in most patients, and temperature was found to have a positive effect on pain (Guedj and Weinberger 1990) . At the beginning of the twentieth century, Strusberg et al. (2002) followed up 82 RA patients, and the results showed that these patients were sensitive to the changes in temperature, and low temperature had a negative impact on their pain. Hamida reported that with a decrease in temperature in the winter, RA patients were more likely to experience pain in their joints (Azzouzi and Ichchou 2020) . However, inconsistent results for the association between RA and temperature were demonstrated. Neogi et al. (2014) suggested that high temperature was related to the higher risk of pain. But Dorleijn et al. (2014) failed to found the association between temperature and pain. In a comprehensive systematic review conducted by Beukenhorst et al. (2020) , no link or opposite associations between temperatures and RA pain were found in different studies, and potential effects in different directions of temperature appeared in different populations. Furthermore, Elcik et al. (2020) stated distinguishing results noted for RA, that warmer air masses contributed to the greatest number of emergency room visits, while colder air masses led to the fewest. Taken together, these findings show that changes in temperature might be associated with the risk of development or deterioration of RA. Therefore, it is worthwhile to explore the relationship between changes in temperature and hospital admissions for RA. By understanding the impact of temperature changes on the risk of admission of patients with RA, timely measures can be taken to manage and alleviate the deterioration of the disease caused by temperature changes. Therefore, we conducted a time-series study to investigate the association between changing temperature metrics and hospital admissions for RA in Anqing, China, and to assess the risk of hospital admissions for RA. Anqing (30.52° N, 117 .05° E) is a prefecture level city in in the southwest of Anhui Province, China, with a total area of 13,589.99 km 2 and a permanent resident population of 4.723 million. Anqing is located in the northern subtropical humid climate zone, with obvious monsoons, four distinct seasons, a humid climate, and sufficient sunlight. The total amount of solar radiation is 112-117 kcal/cm 2 , and the annual average temperature is 14.5-16.6 °C, which is slightly lower in the autumn than in the spring (http:// www. anqing. gov. cn/ zjaq/ zrhj/ index. html). Daily admission data of patients with RA, including data on basic information, rheumatoid factors (RF) and home address of RA patients, admission date, and admission reasons, from January 1, 2015 to December 31, 2019, were obtained from Anqing Municipal Hospital. All patients were diagnosed as RA by rheumatologists according to the RA classification criteria revised by the American College of Rheumatology (Arnett et al. 1988) . Cases from outside Anqing were excluded from this study based on home address of the patients. Daily meteorological information, including data on relative humidity (RH) and maximum and minimum temperatures (℃), was retrieved from the China Meteorological Data Network (http:// data. cma. cn/). Daily mean temperature (℃) was calculated based on maximum and minimum temperatures. Besides, diurnal temperature range (DTR) and temperature change between neighboring days (TCN), as important metrics evaluating the associations between temperature changes and human health, were also calculated (Vicedo-Cabrera et al. 2016; Zhang et al. 2018) . DTR was the range subtracting the daily minimum temperature from the daily maximum temperature. TCN was the change by subtracting the mean temperature of the day from that of the previous day. Daily air quality data, including data on the mean daily concentrations of particulate matter with an aerodynamic diameter of < 10 μm (PM 10 ), particulate matter with an aerodynamic diameter of < 2.5 μm (PM 2.5 ), sulfur dioxide (SO 2 ), nitrogen dioxide (NO 2 ), carbon monoxide (CO), and ozone (O 3 ), were obtained from the National Environmental Monitoring Center of China (https:// www. aqist udy. cn/ histo rydata/). This study was approved by the Ethics Committee of the Anhui Medical University, China (approval no: 2020H007). The correlations between meteorological parameters, air quality, and daily admissions for RA were assessed using Spearman rank correlation analysis. Daily admission for RA was considered a small-probability event and showed an approximate Poisson distribution. Previous studies have shown that different health events have a non-linear relationship with temperature and air quality (Wang et al. 2021; Marinaccio et al. 2019; Liu et al. 2019a) . Environmental factors usually have a delay effect on a disease; hence, it is necessary to use a flexible statistical model to describe the additional time dimension of the exposure-outcome relationship (Gasparrini et al. 2010) . A distributed lag non-linear model (DLNM) considers the potential non-linearity of the dose-response curve and the delayed effect of exposure on the results (Gasparrini 2014; Gasparrini and Leone 2014) . In this study, combined with a quasi-Poisson generalized linear model, a DLNM was used to analyze the associations between meteorological factors, air pollution, and daily admissions for RA. Based on the results of Spearman's correlation analysis, RH, PM 10 , PM 2.5 , SO 2 , NO 2 , CO, and O 3 , all with P value < 0.05, were included in the final regression model as confounding factors ( Supplementary Fig. S1 ). After model debugging and fitting, the final models for three temperature variables evaluating the association between temperature and RA admissions in Anqing were as follow: In this model, t represents the number of observation days; μ t is the number of RA cases admitted on day t; is the intercept of the model, which is a constant; mt represents the daily mean temperature; DTR and TCN are two kinds of temperature changes; β is the matrix coefficient of tm; l is the lag days; and ns denotes the natural cubic spline of this model. Time is the observation days (1-1826 days), which was used to control for the time-varying factor (Bhaskaran et al. 2013) . DOW is a dummy variable in this model to control confounding effect of the day of week. γ is the regression coefficient of DOW. Three degrees of freedom (dfs) were used to adjust for the confounding effect of relative humidity and air pollutions in this model. The quasi-Poisson Akaike information criterion (Q-AIC) was used to assess the model fit and ensure the optimal dfs for different variables (Gasparrini et al. 2010) . Therefore, 5*year is selected to control for the peak or valley of the temperature effect curve (Ma et al. 2013 ). Lags of 12, 12, and 10 days were set up to quantify the association between mean temperature, DTR, TCN, and admissions for RA with the lowest Q-AIC values (Supplementary Fig. S2 ). Subgroup analyses according to gender and age were performed to identify susceptible population. Patients were divided into groups based on gender (male or female), age (< 65 years old or ≥ 65 years old), RF index (positive or negative), and admission route (consulting room visits or emergency room visits). The relative risk (RR) and 95% confidence interval (95% CI) were calculated to assess the lag risk and cumulative risk of admissions for RA. In this model, as for mean temperature and DTR, the median annual temperature was selected as the first reference temperature. Then, the minimum RR is obtained by model prediction, and this temperature was deemed as the most suitable temperature and considered the reference temperature. As for TCN, 0℃ was considered the reference temperature (Cheng et al. 2016) . The 5th percentile of annual temperature was considered the extreme cold temperature, and the 95th percentile of annual temperature was considered the extreme hot temperature (Li et al. 2017; Tobias et al. 2012; Zhao et al. 2019 ). The effect of very high temperature and very low temperature and the cumulative effect of the later lag on hospital admissions for RA were assessed. The packages "dlnm," "splines," and "ggplot2" in the R software version 4.0.3 were used for model analysis and result visualization. In this study, a P value < 0.05 was considered statistically significant. Sensitivity analysis was conducted to evaluate the robustness of the model by changing dfs for ns, including time variable (5-7 dfs) and air pollutant variables (3-5 dfs). If no significant changes were found in the adjusted risk estimate of RA admissions, the model was considered stable. Basic information, rheumatoid factors (RF), and admission reasons of RA patients are shown in Table 1 . In total, 1456 patients with RA, comprising 1053 (72.32%) females and 403 (27.68%) males, were hospitalized in Anqing from 2015 to 2019. Among all patients, 545 (37.43%) were aged below 65 years, and the remaining 911 (62.57%) patients were aged 65 years and above; the mean age of the study population was 58.71 years. Of all patients, 1165 (80.01%) were detected with positive RF, while 291 (19.99%) with negative RF. The main reason for hospitalization of RA patients was pejorative joint pain accounting for 84.62%, followed with pneumonia (9.20%) and regular drug injections (5.98%). In addition, the data on daily meteorological parameters and air quality are shown in Table 2 . The daily counts of admissions for RA ranged from 0 to 8, with the maximum count on September 16, 2019. In Anqing City, the mean temperature was 17.29℃ (from − 4.8 to 33.6℃), the mean DTR was 7.66℃ (from 0.6 to 18.4℃), and the mean TCN was 0℃ (from − 11.6 to 6.4℃). The time-series distributions of total admissions for RA and the minimum, maximum, and mean temperatures from January 1, 2015 to December 31, 2019 are presented in Fig. 1 . The average PM 2.5 , PM 10 , SO 2 , NO 2 , CO, and O 3 levels were 50.43 µg/m 3 , 70.52 µg/m 3 , 14.33 µg/ m 3 , 32.34 µg/m 3 , 0.79 µg/m 3 , and 86.01 µg/m 3 , respectively. Supplementary Fig. S3A and Fig. A2 and B present the exposure-response relationship between mean temperature and the risk of hospital admissions for RA at various lag days. Patients with RA exposing to low environmental temperature were found to have an increased risk of admission. The model showed that the risk of admission in patients with RA was the lowest when the mean temperature is − 4.5℃. Therefore, this temperature was regarded as the reference temperature. The single-day lag and cumulativelag associations of extreme heat temperature (95th percentile = 30℃) and extreme cold temperature (5th percentile = 3℃) with hospital admissions for RA are presented in Table 3 . On day 1, the single-day lags of extreme cold temperature and extreme heat temperature showed significant differences according to hospital admissions for RA with RRs of = 2.20 (95% CI: 1.14-4.21) and 3.46 (95% CI: 1.31-9.15), respectively. Regarding the cumulative-lag association, the risk of admission for RA related to extreme cold temperature increased from lag 0-1 (RR = 1.99, 95% CI: 1.05-3.77) to lag 0-3 (RR = 2.23, 95% CI: 1.19-4.18) and from lag 0-8 (RR = 2.16, 95% CI: 1.03-4.54) to lag 0-12 (RR = 2.59, 95% CI: 1.13-5.90), with the highest risk at lag 0-11 (RR = 2.68, 95% CI: 1.23-5.86). As for extreme heat temperature, only at lag 0-2 (RR = 3.29, 95% CI: 1.17-9.30), the risk of the hospital admission for RA was significant. Figure 3A shows the visual cumulative-lag associations between RA admissions, extreme cold temperature, and extreme heat temperature. The associations between extreme cold temperature and RA hospital admission in different subgroups are presented in Supplementary Table S1. No significant association was found between exposure to extreme cold temperature and the risk of admission in male patients with RA for single-day lag. However, as for multi-day lag, men had increased risks of hospital admission from lag 0-5 (RR = 4.20, 95% CI: 1.01-17.50) to lag 0-12 (RR = 8.33, 95% CI: 1.56-44.58), with the most robust risk at lag 0-10 (RR = 9.83, 95% CI: 1.88-51.40). In female patients exposed to extreme cold temperature, a significant risk of hospital admission was seen only at lag 2 (RR = 1.48, 95% CI: 1.02-2.15) with respect to single-day lag and lag 0-2 (RR = 2.35, 95% CI: 1.11-4.95) with respect to cumulative lag. In patients with RA age < 65 years, the exposure-lag-response association significantly increased at lag 0-9 (RR = 3.75, 95% CI: 1.08-12.98) to lag 0-12 (RR = 4.11, 95% CI: 1.10-15.28); however, the risk related to single-day lag was not significant. In the subgroup of hospitalized patients over 64 years old, no association was found according to both single-day lag and multi-day lag. At the subgroups of RF index, the single-lag effect was only significantly increased in RF-positive population at lag 9 (RR = 1.16, 95% CI: 1.00-1.35), while not significant in RF-negative population. The cumulative RRs increased and got the peak at lag 0-11 (RR = 2.61, 95% CI: 1.11-6.14) for RF-positive group, while RR peaked at lag 0-3 for RF-negative group (RR = 6.05, 95% CI: 1.28-28.66). Furthermore, regarding different admission routes, patients from consulting room had increased risk for hospitalization at extreme cold temperature (RR = 1.19, 95% CI: 1.01-1.40, lag 9). Regarding cumulative risks, hospitalized patients from consulting room and emergency room were both significantly affected with RRs > 1. Supplementary Fig. S3B and Fig. 2C and D present the exposure-response relationship between DTR and the risk of hospital admissions for RA at various lag days. Exposing to low (5th percentile = 1.9℃) and high (95th percentile = 14.2℃) DTRs both had increased risks of RA admission. Fig. 1 The time-series distributions of mean temperature, diurnal temperature range, temperature change between neighboring days, and total RA admission in Anqing from 2015 to 2019. MT: mean temperature; DTR: diurnal temperature range; TCN: temperature change between neighboring days; TA: total admission The DTR model showed that the risk of admission in RA patients was the lowest at DTR = 18℃. Thus, it was considered the reference temperature. The single-day lag and cumulative-lag associations of low DTR (5th percentile = 1.9℃) and high DTR (95th percentile = 14.2℃) with hospital admissions for RA are presented in Table 4 . On day 1 to day 5, the single-day lags of low and high DTRs both showed increased risks for RA hospital admissions, with highest RRs of 1.40 (95% CI: 1.03-1.91) and 1.24 (95% CI: 1.01-1.53) at lag 0 day, respectively. The cumulative-lag effects of low and high DTRs are visualized in Fig. 3B . The risks of admission for RA related to low and high DTRs increased from lag 0-1 to lag 0-12 with the highest RRs of 8.58 (95% CI: 1.82-40.48) and 5.03 (95% CI: 1.60-15.86), respectively. The associations between low DTR, high DTR, and RA hospital admission in different subgroups are presented in Supplementary Table S2 and S3. The subgroup analyses of different gender show that exposure to both low DTR and high DTR increased the hospitalized risk for female with highest value at lag 12 (RR = 1.97, 95% CI: 1.03-3.75; RR = 2.02, 95% CI: 1.03-3.96). For the population from emergency room, the significant and lowest risk of hospitalization appeared at lag 6 (RR = 1.45, 95% CI: 1.06-1.98). Besides, the cumulative risk for these patients was lowest at lag 0-4 (RR = 4.69, 95% CI: 1.15-19.11) and remained Fig. 2 The association between temperature changes and the risk of hospital admissions for RA. A The contour plot of exposure-lag-response risk of daily mean temperature for RA admission. B The overall cumulative exposure-response associations between daily mean temperature and admission risk for RA. C The contour plot of exposure-lag-response risk of DTR for RA admission. D The overall cumulative exposureresponse associations between DTR and admission risk for RA. E The contour plot of exposurelag-response risk of TCN for RA admission. F The overall cumulative exposure-response associations between TCN and admission risk for RA. DTR: diurnal temperature range; TCN: temperature change between neighboring days significantly increased. Only the emergency room population had significantly increased risk of hospitalization at high DTR in both single and cumulative-lag days. Supplementary Fig. S3C and Fig. 2E and F present the exposure-response relationship between TCN and the risk of hospital admissions for RA at various lag days. Marginal association between TCN and RA admissions was found. The marginal risk of admission in RA patients was found when exposed to high TCN (95th percentile = 2.9℃). The largest single-day effect of high TCN was observed at lag 10 (RR = 1.11, 95% CI: 1.01-1.23) (Table 5 ). However, no significant cumulative effect of TCN was found (Fig. 3C) . After stratifying by gender, adverse effects of high TCN were found on patients who aged < 65 (RR = 1.22, 95% CI: 1.03-1.45, lag 9; RR = 3.88, 95% CI: 1.05-14.29, lag 0-10). Surprisingly, opposite results were found in the subgroup analyses of different admission routes. Increased risk for RA hospitalization was found in population from consulting Table S4 ). The sensitivity analysis was conducted by changing the dfs for the time variable (5-7 dfs) and other environmental factors (3-5 dfs). When controlling the dfs of variables in the mean temperature model, no significant changes were appeared (Supplementary Fig. S4 and S5) . Besides, after adjusting the dfs in the DTR model, there were slightly changes in the results without statistical significance (Supplementary Fig. S6 and S7 ). With changing dfs of variables, the effects of TCN on RA admissions did not sharply change and appeared stable tendency with the original one (Supplementary Fig. S8 and S9 ). The results of sensitivity analyses indicate robust association between temperature and hospital admissions for RA. To our knowledge, few studies have evaluated the association between temperature and the risk of admission for RA. In this study, temperature changes were considered factors that might increase the risk of admission for RA with significantly delayed effects. According to the results of the subgroup analysis, males were more likely to be vulnerable in a low-temperature environment and admitted to the hospital for further treatment. Moreover, patients with RA aged < 65 years had a higher risk of admission in the cold environment than those aged ≥ 65 years. This study found a significant association between low temperature and admission for RA. Based on their experiences in daily life, people have maintained an old belief that weather can affect joint pain. Previous studies have indicated that exposure to cold environments might be associated with a risk of musculoskeletal disorders, such as back and neck (Burström et al. 2013) , and white fingers (Burström et al. 2010) . In a Spanish study, 245 patients who visited the emergency room because of worsening pain due to RA were included (Abasolo et al. 2013) . The patients' symptoms and weather conditions were analyzed using a case-crossover design, and the results indicated that a low mean temperature might worsen symptoms in patients with RA (Abasolo et al. 2013) . Some study results have suggested that extreme changes in temperature might cause an inflammatory reaction that leads to adverse effect related to RA. Based on research results, some interesting ideas were put forward that extreme temperature changes may cause inflammatory reaction and adverse effect. For example, a decrease in temperature can increase C-reactive protein, especially in the winter and spring (Sung 2006) . The relationship between rapidly decreased temperature and increased IL-6 was significant with respect to the risk of cardiovascular disease (Schneider et al. 2008) . Furthermore, sudden changes in temperature could affect the homeostasis of B and T lymphocytes, resulting in autoimmune disorders (Maes et al. 1994) . These inflammatory factors play important roles not only in cardiovascular disease but also in RA development. In the Epidemiological Investigation of Rheumatoid Arthritis study, a case-control study based on Swedish population, increased risk of RA was found related to working in cold environments, with respect to both anticitrullinated protein antibody-positive and anticitrullinated protein antibodynegative RA (Zeng et al. 2017) . The subgroup analyses demonstrated that the lag effects of temperature changes were different between patients aged below 65 years and those aged ≥ 65 years, and the effect of changes in temperature was significant in RA patients aged below 65 years. A previous study showed that the association between the worsening of RA symptoms and exposure to low temperatures was significant only in patients aged between 50 and 65 years (Abasolo et al. 2013) , similar to our finding. A previous study reported that young people rated more pain (Lue et al. 2018) , which may promote young RA patients hospitalized when they got the exacerbation of joint pain. The other opposite direction is that RA patients aged ≥ 65 who have the need to go to the hospital, but actually don't go for various reasons. High medical costs may prevent the elderly who cannot afford treatment (Zeng et al. 2021) . They visit the hospital only when the pain becomes unbearable. This might be a possible reason for the lag effect of low temperature. Furthermore, elderly patients might have other severe conditions; hence, the aggravation of RA symptoms caused by exposure to low temperatures might not get due attention, leading to lower hospitalization rates. The subgroup analysis according to gender indicated that the risk of admission for RA was higher in females with RA than in males with RA. Compared to males, females were more sensitive to pain (Fillingim et al. 2009) , and when female patients feel pain due to temperature changes, they may go to the hospital for treatment. There are some strengths to this study. First, this study evaluated the effects of temperature changes on the risk of admission for RA, considering possible lag effect using the DLNM model. Second, the RH and air quality indicators were included as confounding variables to assess the effects of temperature changes, which made our results more convincing. Third, subgroup analyses with respect to age, gender, RF index, and admission route of RA patients were performed to identify the risk of admission for RA in patients with different characteristics. We hope that our results can emphasize the need for considering the impact of temperature change on the health of patients with RA and reduce both the health and economic burdens of RA. However, the study also has several limitations. As this was an ecological study, it was difficult to rule out the existence of ecological fallacy. Second, the data on admission for RA were collected from only one city, which cannot represent the data in other cities or regions. Third, the small samples of RA admission in Anqing City may affect the evaluation of the exposure-response lag association between temperature changes and RA admission. Furthermore, due to the lack of information on patient exposure, it was difficult to determine the patients' exposure to outdoor temperatures, which might bias the results. In this study, we found the association between exposure to extreme cold temperature, low and high DTRs, low TCN, and admission risk of RA in Anqing, China. Furthermore, subgroup analyses indicated that male patients and patients younger than 65 years old had higher risk of hospitalized admissions. To avoid RA worsening, RA patients should pay more attention to their health and keep warm according to the temperature changes. The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00484-021-02207-9. Author contribution Li-Juan Huang: writing-original draft, project administration, resources. Prevalence and incidence of upper respiratory tract infection events are elevated prior to the development of rheumatoid arthritis in first-degree relatives The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis Seasonal and weather effects on rheumatoid arthritis: myth or reality? Are weather conditions associated with chronic musculoskeletal pain? Review of Results and Methodologies Time series regression studies in environmental epidemiology White fingers, cold environment, and vibration-exposure among Swedish construction workers Back and neck pain due to working in a cold environment: a cross-sectional study of male construction workers Impacts of temperature change on ambulance dispatches and seasonal effect modification Interferon gamma gene polymorphism and susceptibility to, and severity of, rheumatoid arthritis The global burden of rheumatoid arthritis: estimates from the global burden of disease 2010 study Associations between weather conditions and clinical symptoms in patients with hip osteoarthritis: a 2-year cohort study Relationship between synoptic weather type and emergency department visits for different types of pain across the Triangle region of North Carolina Sex, gender, and pain: a review of recent clinical and experimental findings Modeling exposure-lag-response associations with distributed lag non-linear models Attributable risk from distributed lag models Distributed lag non-linear models Effect of weather conditions on rheumatic patients Genetic implications in the pathogenesis of rheumatoid arthritis; an updated review Impact of ambient temperature and atmospheric evaporation on the incidence of acute deep venous thrombosis in the northeast of China Shortterm exposure to ambient air pollution and daily atherosclerotic heart disease mortality in a cool climate Impact and timing of smoking cessation on reducing risk of rheumatoid arthritis among women in the nurses' health studies Thermal pain tolerance and pain rating in normal subjects: gender and age effects Applied mixed generalized additive model to assess the effect of temperature on the incidence of bacillary dysentery and its forecast Seasonal variation in peripheral blood leukocyte subsets and in serum interleukin-6, and soluble interleukin-2 and -6 receptor concentrations in normal volunteers Associations between ambient temperature and hepatobiliary and renal hospitalizations in California Nationwide epidemiological study for estimating the effect of extreme outdoor temperature on occupational injuries in Italy Burden of musculoskeletal disorders in the Eastern Mediterranean Region, 1990-2013: findings from the Global Burden of Disease Study Influence of seasonal changes on disease activity and distribution of affected joints in rheumatoid arthritis Relation of temperature and humidity to the risk of recurrent gout attacks Weather effects in rheumatoid arthritis: from controversy to consensus Global, regional and national burden of rheumatoid arthritis 1990-2017: a systematic analysis of the Global Burden of Disease study 2017 Air temperature and inflammatory responses in myocardial infarction survivors Seasonality and ambient temperature at time of conception in term-born individuals-influences on cardiovascular disease and obesity in adult life Influence of weather conditions on rheumatic pain Seasonal variation of C-reactive protein in apparently healthy Koreans Inverse association between air pressure and rheumatoid arthritis synovitis Mortality on extreme heat days using official thresholds in Spain: a multi-city time series analysis Weather conditions can influence rheumatic diseases Associations of inter-and intraday temperature change with mortality Impact of ambient temperature on cardiovascular disease hospital admissions in farmers in China's Western suburbs Working in cold environment and risk of developing rheumatoid arthritis: results from the Swedish EIRA case-control study Healthcare-seeking behavior among Chinese older adults: patterns and predictive factors The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment Association of diurnal temperature range with daily mortality in England and Wales: a nationwide time-series study Morbidity burden of respiratory diseases attributable to ambient temperature: a case study in a subtropical city in China