Comparison for the effects of different components of temperature variability on mortality: A multi-country time-series study.

BACKGROUND: Temperature variability (TV) is associated with increased mortality risk. However, it is still unknown whether intra-day or inter-day TV has different effects. OBJECTIVES: We aimed to assess the association of intra-day TV and inter-day TV with all-cause, cardiovascular, and respiratory mortality. METHODS: We collected data on total, cardiovascular, and respiratory mortality and meteorology from 758 locations in 47 countries or regions from 1972 to 2020. We defined inter-day TV as the standard deviation (SD) of daily mean temperatures across the lag interval, and intra-day TV as the average SD of minimum and maximum temperatures on each day. In the first stage, inter-day and intra-day TVs were modelled simultaneously in the quasi-Poisson time-series model for each location. In the second stage, a multi-level analysis was used to pool the location-specific estimates. RESULTS: Overall, the mortality risk due to each interquartile range [IQR] increase was higher for intra-day TV than for inter-day TV. The risk increased by 0.59% (95% confidence interval [CI]: 0.53, 0.65) for all-cause mortality, 0.64% (95% CI: 0.56, 0.73) for cardiovascular mortality, and 0.65% (95% CI: 0.49, 0.80) for respiratory mortality per IQR increase in intra-day TV0-7 (0.9 °C). An IQR increase in inter-day TV0-7 (1.6 °C) was associated with 0.22% (95% CI: 0.18, 0.26) increase in all-cause mortality, 0.44% (95% CI: 0.37, 0.50) increase in cardiovascular mortality, and 0.31% (95% CI: 0.21, 0.41) increase in respiratory mortality. The proportion of all-cause deaths attributable to intra-day TV0-7 and inter-day TV0-7 was 1.45% and 0.35%, respectively. The mortality risks varied by lag interval, climate area, season, and climate type. CONCLUSIONS: Our results indicated that intra-day TV may explain the main part of the mortality risk related to TV and suggested that comprehensive evaluations should be proposed in more countries to help protect human health.

Global, regional, and national burden of heatwave-related mortality from 1990 to 2019: A three-stage modelling study.

BACKGROUND: The regional disparity of heatwave-related mortality over a long period has not been sufficiently assessed across the globe, impeding the localisation of adaptation planning and risk management towards climate change. We quantified the global mortality burden associated with heatwaves at a spatial resolution of 0.5°×0.5° and the temporal change from 1990 to 2019. METHODS AND FINDINGS: We collected data on daily deaths and temperature from 750 locations of 43 countries or regions, and 5 meta-predictors in 0.5°×0.5° resolution across the world. Heatwaves were defined as location-specific daily mean temperature ≥95th percentiles of year-round temperature range with duration ≥2 days. We first estimated the location-specific heatwave-mortality association. Secondly, a multivariate meta-regression was fitted between location-specific associations and 5 meta-predictors, which was in the third stage used with grid cell-specific meta-predictors to predict grid cell-specific association. Heatwave-related excess deaths were calculated for each grid and aggregated. During 1990 to 2019, 0.94% (95% CI: 0.68-1.19) of deaths [i.e., 153,078 cases (95% eCI: 109,950-194,227)] per warm season were estimated to be from heatwaves, accounting for 236 (95% eCI: 170-300) deaths per 10 million residents. The ratio between heatwave-related excess deaths and all premature deaths per warm season remained relatively unchanged over the 30 years, while the number of heatwave-related excess deaths per 10 million residents per warm season declined by 7.2% per decade in comparison to the 30-year average. Locations with the highest heatwave-related death ratio and rate were in Southern and Eastern Europe or areas had polar and alpine climates, and/or their residents had high incomes. The temporal change of heatwave-related mortality burden showed geographic disparities, such that locations with tropical climate or low incomes were observed with the greatest decline. The main limitation of this study was the lack of data from certain regions, e.g., Arabian Peninsula and South Asia. CONCLUSIONS: Heatwaves were associated with substantial mortality burden that varied spatiotemporally over the globe in the past 30 years. The findings indicate the potential benefit of governmental actions to enhance health sector adaptation and resilience, accounting for inequalities across communities.

Global, regional, and national burden of mortality associated with cold spells during 2000-19: a three-stage modelling study.

BACKGROUND: Exposure to cold spells is associated with mortality. However, little is known about the global mortality burden of cold spells. METHODS: A three-stage meta-analytical method was used to estimate the global mortality burden associated with cold spells by means of a time series dataset of 1960 locations across 59 countries (or regions). First, we fitted the location-specific, cold spell-related mortality associations using a quasi-Poisson regression with a distributed lag non-linear model with a lag period of up to 21 days. Second, we built a multivariate meta-regression model between location-specific associations and seven predictors. Finally, we predicted the global grid-specific cold spell-related mortality associations during 2000-19 using the fitted meta-regression model and the yearly grid-specific meta-predictors. We calculated the annual excess deaths, excess death ratio (excess deaths per 1000 deaths), and excess death rate (excess deaths per 100 000 population) due to cold spells for each grid across the world. FINDINGS: Globally, 205 932 (95% empirical CI [eCI] 162 692-250 337) excess deaths, representing 3·81 (95% eCI 2·93-4·71) excess deaths per 1000 deaths (excess death ratio), and 3·03 (2·33-3·75) excess deaths per 100 000 population (excess death rate) were associated with cold spells per year between 2000 and 2019. The annual average global excess death ratio in 2016-19 increased by 0·12 percentage points and the excess death rate in 2016-19 increased by 0·18 percentage points, compared with those in 2000-03. The mortality burden varied geographically. The excess death ratio and rate were highest in Europe, whereas these indicators were lowest in Africa. Temperate climates had higher excess death ratio and rate associated with cold spells than other climate zones. INTERPRETATION: Cold spells are associated with substantial mortality burden around the world with geographically varying patterns. Although the number of cold spells has on average been decreasing since year 2000, the public health threat of cold spells remains substantial. The findings indicate an urgency of taking local and regional measures to protect the public from the mortality burdens of cold spells. FUNDING: Australian Research Council, Australian National Health and Medical Research Council, EU's Horizon 2020 Project Exhaustion.

Ozone-related acute excess mortality projected to increase in the absence of climate and air quality controls consistent with the Paris Agreement.

Short-term exposure to ground-level ozone in cities is associated with increased mortality and is expected to worsen with climate and emission changes. However, no study has yet comprehensively assessed future ozone-related acute mortality across diverse geographic areas, various climate scenarios, and using CMIP6 multi-model ensembles, limiting our knowledge on future changes in global ozone-related acute mortality and our ability to design targeted health policies. Here, we combine CMIP6 simulations and epidemiological data from 406 cities in 20 countries or regions. We find that ozone-related deaths in 406 cities will increase by 45 to 6,200 deaths/year between 2010 and 2014 and between 2050 and 2054, with attributable fractions increasing in all climate scenarios (from 0.17% to 0.22% total deaths), except the single scenario consistent with the Paris Climate Agreement (declines from 0.17% to 0.15% total deaths). These findings stress the need for more stringent air quality regulations, as current standards in many countries are inadequate.

Environmental and social inequities in continental France: an analysis of exposure to heat, air pollution, and lack of vegetation.

BACKGROUND: Cumulative environmental exposures and social deprivation increase health vulnerability and limit the capacity of populations to adapt to climate change. OBJECTIVE: Our study aimed at providing a fine-scale characterization of exposure to heat, air pollution, and lack of vegetation in continental France between 2000 and 2018, describing spatiotemporal trends and environmental hotspots (i.e., areas that cumulate the highest levels of overexposure), and exploring any associations with social deprivation. METHODS: The European (EDI) and French (FDep) social deprivation indices, the normalized difference vegetation index, daily ambient temperatures, particulate matter (PM2.5 and PM10), nitrogen dioxide, and ozone (O3) concentrations were estimated for 48,185 French census districts. Reference values were chosen to characterize (over-)exposure. Hotspots were defined as the areas cumulating the highest overexposure to temperature, air pollution, and lack of vegetation. Associations between heat overexposure or hotspots and social deprivation were assessed using logistic regressions. RESULTS: Overexposure to heat was higher in 2015-2018 compared with 2000-2014. Exposure to all air pollutants except for O3 decreased during the study period. In 2018, more than 79% of the urban census districts exceeded the 2021 WHO air quality guidelines. The evolution of vegetation density between 2000 and 2018 was heterogeneous across continental France. In urban areas, the most deprived census districts were at a higher risk of being hotspots (odds ratio (OR): 10.86, 95% CI: 9.87-11.98 using EDI and OR: 1.07, 95% CI: 1.04-1.11 using FDep). IMPACT STATEMENT: We studied cumulative environmental exposures and social deprivation in French census districts. The 2015-2018 period showed the highest overexposure to heat between 2000 and 2018. In 2018, the air quality did not meet the 2021 WHO guidelines in most census districts and 8.6 million people lived in environmental hotspots. Highly socially deprived urban areas had a higher risk of being in a hotspot. This study proposes for the first time, a methodology to identify hotspots of exposure to heat, air pollution, and lack of vegetation and their associations with social deprivation at a national level.

Interactive effects of ambient fine particulate matter and ozone on daily mortality in 372 cities: two stage time series analysis.

OBJECTIVE: To investigate potential interactive effects of fine particulate matter (PM2.5) and ozone (O3) on daily mortality at global level. DESIGN: Two stage time series analysis. SETTING: 372 cities across 19 countries and regions. POPULATION: Daily counts of deaths from all causes, cardiovascular disease, and respiratory disease. MAIN OUTCOME MEASURE: Daily mortality data during 1994-2020. Stratified analyses by co-pollutant exposures and synergy index (>1 denotes the combined effect of pollutants is greater than individual effects) were applied to explore the interaction between PM2.5 and O3 in association with mortality. RESULTS: During the study period across the 372 cities, 19.3 million deaths were attributable to all causes, 5.3 million to cardiovascular disease, and 1.9 million to respiratory disease. The risk of total mortality for a 10 µg/m3 increment in PM2.5 (lag 0-1 days) ranged from 0.47% (95% confidence interval 0.26% to 0.67%) to 1.25% (1.02% to 1.48%) from the lowest to highest fourths of O3 concentration; and for a 10 µg/m3 increase in O3 ranged from 0.04% (-0.09% to 0.16%) to 0.29% (0.18% to 0.39%) from the lowest to highest fourths of PM2.5 concentration, with significant differences between strata (P for interaction <0.001). A significant synergistic interaction was also identified between PM2.5 and O3 for total mortality, with a synergy index of 1.93 (95% confidence interval 1.47 to 3.34). Subgroup analyses showed that interactions between PM2.5 and O3 on all three mortality endpoints were more prominent in high latitude regions and during cold seasons. CONCLUSION: The findings of this study suggest a synergistic effect of PM2.5 and O3 on total, cardiovascular, and respiratory mortality, indicating the benefit of coordinated control strategies for both pollutants.

Joint effect of heat and air pollution on mortality in 620 cities of 36 countries.

BACKGROUND: The epidemiological evidence on the interaction between heat and ambient air pollution on mortality is still inconsistent. OBJECTIVES: To investigate the interaction between heat and ambient air pollution on daily mortality in a large dataset of 620 cities from 36 countries. METHODS: We used daily data on all-cause mortality, air temperature, particulate matter ≤ 10 µm (PM10), PM ≤ 2.5 µm (PM2.5), nitrogen dioxide (NO2), and ozone (O3) from 620 cities in 36 countries in the period 1995-2020. We restricted the analysis to the six consecutive warmest months in each city. City-specific data were analysed with over-dispersed Poisson regression models, followed by a multilevel random-effects meta-analysis. The joint association between air temperature and air pollutants was modelled with product terms between non-linear functions for air temperature and linear functions for air pollutants. RESULTS: We analyzed 22,630,598 deaths. An increase in mean temperature from the 75th to the 99th percentile of city-specific distributions was associated with an average 8.9 % (95 % confidence interval: 7.1 %, 10.7 %) mortality increment, ranging between 5.3 % (3.8 %, 6.9 %) and 12.8 % (8.7 %, 17.0 %), when daily PM10 was equal to 10 or 90 µg/m3, respectively. Corresponding estimates when daily O3 concentrations were 40 or 160 µg/m3 were 2.9 % (1.1 %, 4.7 %) and 12.5 % (6.9 %, 18.5 %), respectively. Similarly, a 10 µg/m3 increment in PM10 was associated with a 0.54 % (0.10 %, 0.98 %) and 1.21 % (0.69 %, 1.72 %) increase in mortality when daily air temperature was set to the 1st and 99th city-specific percentiles, respectively. Corresponding mortality estimate for O3 across these temperature percentiles were 0.00 % (-0.44 %, 0.44 %) and 0.53 % (0.38 %, 0.68 %). Similar effect modification results, although slightly weaker, were found for PM2.5 and NO2. CONCLUSIONS: Suggestive evidence of effect modification between air temperature and air pollutants on mortality during the warm period was found in a global dataset of 620 cities.

Global short-term mortality risk and burden associated with tropical cyclones from 1980 to 2019: a multi-country time-series study.

BACKGROUND: The global spatiotemporal pattern of mortality risk and burden attributable to tropical cyclones is unclear. We aimed to evaluate the global short-term mortality risk and burden associated with tropical cyclones from 1980 to 2019. METHODS: The wind speed associated with cyclones from 1980 to 2019 was estimated globally through a parametric wind field model at a grid resolution of 0·5°â€ˆ× 0·5°. A total of 341 locations with daily mortality and temperature data from 14 countries that experienced at least one tropical cyclone day (a day with maximum sustained wind speed associated with cyclones ≥17·5 m/s) during the study period were included. A conditional quasi-Poisson regression with distributed lag non-linear model was applied to assess the tropical cyclone-mortality association. A meta-regression model was fitted to evaluate potential contributing factors and estimate grid cell-specific tropical cyclone effects. FINDINGS: Tropical cyclone exposure was associated with an overall 6% (95% CI 4-8) increase in mortality in the first 2 weeks following exposure. Globally, an estimate of 97 430 excess deaths (95% empirical CI [eCI] 71 651-126 438) per decade were observed over the 2 weeks following exposure to tropical cyclones, accounting for 20·7 (95% eCI 15·2-26·9) excess deaths per 100 000 residents (excess death rate) and 3·3 (95% eCI 2·4-4·3) excess deaths per 1000 deaths (excess death ratio) over 1980-2019. The mortality burden exhibited substantial temporal and spatial variation. East Asia and south Asia had the highest number of excess deaths during 1980-2019: 28 744 (95% eCI 16 863-42 188) and 27 267 (21 157-34 058) excess deaths per decade, respectively. In contrast, the regions with the highest excess death ratios and rates were southeast Asia and Latin America and the Caribbean. From 1980-99 to 2000-19, marked increases in tropical cyclone-related excess death numbers were observed globally, especially for Latin America and the Caribbean and south Asia. Grid cell-level and country-level results revealed further heterogeneous spatiotemporal patterns such as the high and increasing tropical cyclone-related mortality burden in Caribbean countries or regions. INTERPRETATION: Globally, short-term exposure to tropical cyclones was associated with a significant mortality burden, with highly heterogeneous spatiotemporal patterns. In-depth exploration of tropical cyclone epidemiology for those countries and regions estimated to have the highest and increasing tropical cyclone-related mortality burdens is urgently needed to help inform the development of targeted actions against the increasing adverse health impacts of tropical cyclones under a changing climate. FUNDING: Australian Research Council and Australian National Health and Medical Research Council.

Heat-related cardiorespiratory mortality: Effect modification by air pollution across 482 cities from 24 countries.

BACKGROUND: Evidence on the potential interactive effects of heat and ambient air pollution on cause-specific mortality is inconclusive and limited to selected locations. OBJECTIVES: We investigated the effects of heat on cardiovascular and respiratory mortality and its modification by air pollution during summer months (six consecutive hottest months) in 482 locations across 24 countries. METHODS: Location-specific daily death counts and exposure data (e.g., particulate matter with diameters ≤ 2.5 µm [PM2.5]) were obtained from 2000 to 2018. We used location-specific confounder-adjusted Quasi-Poisson regression with a tensor product between air temperature and the air pollutant. We extracted heat effects at low, medium, and high levels of pollutants, defined as the 5th, 50th, and 95th percentile of the location-specific pollutant concentrations. Country-specific and overall estimates were derived using a random-effects multilevel meta-analytical model. RESULTS: Heat was associated with increased cardiorespiratory mortality. Moreover, the heat effects were modified by elevated levels of all air pollutants in most locations, with stronger effects for respiratory than cardiovascular mortality. For example, the percent increase in respiratory mortality per increase in the 2-day average summer temperature from the 75th to the 99th percentile was 7.7% (95% Confidence Interval [CI] 7.6-7.7), 11.3% (95%CI 11.2-11.3), and 14.3% (95% CI 14.1-14.5) at low, medium, and high levels of PM2.5, respectively. Similarly, cardiovascular mortality increased by 1.6 (95%CI 1.5-1.6), 5.1 (95%CI 5.1-5.2), and 8.7 (95%CI 8.7-8.8) at low, medium, and high levels of O3, respectively. DISCUSSION: We observed considerable modification of the heat effects on cardiovascular and respiratory mortality by elevated levels of air pollutants. Therefore, mitigation measures following the new WHO Air Quality Guidelines are crucial to enhance better health and promote sustainable development.

Changes in the temperature-mortality relationship in France: Limited evidence of adaptation to a new climate.

CONTEXT: Documenting trends in the health impacts of ambient temperature is key to supporting adaptation strategies to climate change. This paper explores changes in the temperature-related mortality in 18 French urban centers between 1970 and 2015. METHOD: A multicentric time-series design with time-varying distributed lag nonlinear models was adopted to model the shape of the relationship and assess temporal changes in risks and impacts. RESULTS: The general shape of the temperature-mortality relationship did not change over time, except for an increasing risk at very low percentiles and a decreasing risk at very high percentiles. The relative risk at the 99.9th percentile compared to the 50th percentile of the 1970-2015 temperature distribution decreased from 2.33 [95% confidence interval (CI): 1.95:2.79] in 1975 to 1.33 [95% CI: 1.14:1.55] in 2015. Between 1970 and 2015, 302,456 [95% CI: 292,723:311,392] deaths were attributable to non-optimal temperatures, corresponding to 5.5% [95% CI: 5.3:5.6] of total mortality. This burden decreased progressively, representing 7.2% [95% CI: 6.7:7.7] of total mortality in the 1970s to 3.4% [95% CI: 3.2:3.6] in the 2000s. However, the contribution of hot temperatures to this burden (higher than the 90th percentile) increased. DISCUSSION: Despite the decreasing relative risk, the fraction of mortality attributable to extreme heat increased between 1970 and 2015, thus highlighting the need for proactive adaptation.

On the supra-linearity of the relationship between air pollution, mortality and hospital admission in 18 French cities.

PURPOSE: Understanding the relationship between an environmental determinant and a given health outcome is key to inform public health policies. The short-term mortality and morbidity responses to outdoor air pollutants are traditionally assessed as a log-linear relationship, but few studies suggest a possible deviation from linearity. This paper investigates the shape of the relationship between ozone, NO2 and fine particulate matter (PM10 and PM2.5), mortality and hospital admissions in 18 French cities between 2000 and 2017. METHOD: A multi-centric time series design, using quasi-Poisson generalized additive models, was used. Four approaches were compared to model concentration-response curves (log-linear, piecewise-linear with a priori defined breakpoints, piecewise-linear with no a priori breakpoint and cubic spline). RESULTS: All the models indicated evidence of supra-linearity between PM10, PM2.5, NO2, mortality and hospital admissions. For instance, with a log-linear model, a 10 µg/m3 increase in PM2.5 was associated with a 0.4% [95% CI 0.2; 0.7] increase in non-accidental mortality. When using a piecewise model with a priori set breakpoint at 10 µg/m3, the mortality increase was 3.8% [4.4; 6.3] below 10 µg/m3, and 0.3% [0; 0.6] above. Non-significant impacts of ozone were found for concentrations below 90 µg/m3 to 120 µg/m3, with some variability in the identified threshold across the heath indicator studied. CONCLUSION: The supra-linearity of the relationship between PM10, PM2.5, NO2, mortality and hospital admissions supports the need to further reduce air pollution concentrations well below regulatory values.

The association of COVID-19 incidence with temperature, humidity, and UV radiation - A global multi-city analysis.

BACKGROUND AND AIM: The associations between COVID-19 transmission and meteorological factors are scientifically debated. Several studies have been conducted worldwide, with inconsistent findings. However, often these studies had methodological issues, e.g., did not exclude important confounding factors, or had limited geographic or temporal resolution. Our aim was to quantify associations between temporal variations in COVID-19 incidence and meteorological variables globally. METHODS: We analysed data from 455 cities across 20 countries from 3 February to 31 October 2020. We used a time-series analysis that assumes a quasi-Poisson distribution of the cases and incorporates distributed lag non-linear modelling for the exposure associations at the city-level while considering effects of autocorrelation, long-term trends, and day of the week. The confounding by governmental measures was accounted for by incorporating the Oxford Governmental Stringency Index. The effects of daily mean air temperature, relative and absolute humidity, and UV radiation were estimated by applying a meta-regression of local estimates with multi-level random effects for location, country, and climatic zone. RESULTS: We found that air temperature and absolute humidity influenced the spread of COVID-19 over a lag period of 15 days. Pooling the estimates globally showed that overall low temperatures (7.5 °C compared to 17.0 °C) and low absolute humidity (6.0 g/m3 compared to 11.0 g/m3) were associated with higher COVID-19 incidence (RR temp =1.33 with 95%CI: 1.08; 1.64 and RR AH =1.33 with 95%CI: 1.12; 1.57). RH revealed no significant trend and for UV some evidence of a positive association was found. These results were robust to sensitivity analysis. However, the study results also emphasise the heterogeneity of these associations in different countries. CONCLUSION: Globally, our results suggest that comparatively low temperatures and low absolute humidity were associated with increased risks of COVID-19 incidence. However, this study underlines regional heterogeneity of weather-related effects on COVID-19 transmission.

Big data, big consortia, and pain: UK Biobank, PAINSTORM, and DOLORisk.

Chronic pain (CP) is a common and often debilitating disorder that has major social and economic impacts. A subset of patients develop CP that significantly interferes with their activities of daily living and requires a high level of healthcare support. The challenge for treating physicians is in preventing the onset of refractory CP or effectively managing existing pain. To be able to do this, it is necessary to understand the risk factors, both genetic and environmental, for the onset of CP and response to treatment, as well as the pathogenesis of the disorder, which is highly heterogenous. However, studies of CP, particularly pain with neuropathic characteristics, have been hindered by a lack of consensus on phenotyping and data collection, making comparisons difficult. Furthermore, existing cohorts have suffered from small sample sizes meaning that analyses, especially genome-wide association studies, are insufficiently powered. The key to overcoming these issues is through the creation of large consortia such as DOLORisk and PAINSTORM and biorepositories, such as UK Biobank, where a common approach can be taken to CP phenotyping, which allows harmonisation across different cohorts and in turn increased study power. This review describes the approach that was used for studying neuropathic pain in DOLORisk and how this has informed current projects such as PAINSTORM, the rephenotyping of UK Biobank, and other endeavours. Moreover, an overview is provided of the outputs from these studies and the lessons learnt for future projects.

Effect modification of greenness on the association between heat and mortality: A multi-city multi-country study.

BACKGROUND: Identifying how greenspace impacts the temperature-mortality relationship in urban environments is crucial, especially given climate change and rapid urbanization. However, the effect modification of greenspace on heat-related mortality has been typically focused on a localized area or single country. This study examined the heat-mortality relationship among different greenspace levels in a global setting. METHODS: We collected daily ambient temperature and mortality data for 452 locations in 24 countries and used Enhanced Vegetation Index (EVI) as the greenspace measurement. We used distributed lag non-linear model to estimate the heat-mortality relationship in each city and the estimates were pooled adjusting for city-specific average temperature, city-specific temperature range, city-specific population density, and gross domestic product (GDP). The effect modification of greenspace was evaluated by comparing the heat-related mortality risk for different greenspace groups (low, medium, and high), which were divided into terciles among 452 locations. FINDINGS: Cities with high greenspace value had the lowest heat-mortality relative risk of 1·19 (95% CI: 1·13, 1·25), while the heat-related relative risk was 1·46 (95% CI: 1·31, 1·62) for cities with low greenspace when comparing the 99th temperature and the minimum mortality temperature. A 20% increase of greenspace is associated with a 9·02% (95% CI: 8·88, 9·16) decrease in the heat-related attributable fraction, and if this association is causal (which is not within the scope of this study to assess), such a reduction could save approximately 933 excess deaths per year in 24 countries. INTERPRETATION: Our findings can inform communities on the potential health benefits of greenspaces in the urban environment and mitigation measures regarding the impacts of climate change. FUNDING: This publication was developed under Assistance Agreement No. RD83587101 awarded by the U.S. Environmental Protection Agency to Yale University. It has not been formally reviewed by EPA. The views expressed in this document are solely those of the authors and do not necessarily reflect those of the Agency. EPA does not endorse any products or commercial services mentioned in this publication. Research reported in this publication was also supported by the National Institute on Minority Health and Health Disparities of the National Institutes of Health under Award Number R01MD012769. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Also, this work has been supported by the National Research Foundation of Korea (2021R1A6A3A03038675), Medical Research Council-UK (MR/V034162/1 and MR/R013349/1), Natural Environment Research Council UK (Grant ID: NE/R009384/1), Academy of Finland (Grant ID: 310372), European Union's Horizon 2020 Project Exhaustion (Grant ID: 820655 and 874990), Czech Science Foundation (22-24920S), Emory University's NIEHS-funded HERCULES Center (Grant ID: P30ES019776), and Grant CEX2018-000794-S funded by MCIN/AEI/ 10.13039/501100011033 The funders had no role in the design, data collection, analysis, interpretation of results, manuscript writing, or decision to publication.

Classification of painful or painless diabetic peripheral neuropathy and identification of the most powerful predictors using machine learning models in large cross-sectional cohorts.

BACKGROUND: To improve the treatment of painful Diabetic Peripheral Neuropathy (DPN) and associated co-morbidities, a better understanding of the pathophysiology and risk factors for painful DPN is required. Using harmonised cohorts (N = 1230) we have built models that classify painful versus painless DPN using quality of life (EQ5D), lifestyle (smoking, alcohol consumption), demographics (age, gender), personality and psychology traits (anxiety, depression, personality traits), biochemical (HbA1c) and clinical variables (BMI, hospital stay and trauma at young age) as predictors. METHODS: The Random Forest, Adaptive Regression Splines and Naive Bayes machine learning models were trained for classifying painful/painless DPN. Their performance was estimated using cross-validation in large cross-sectional cohorts (N = 935) and externally validated in a large population-based cohort (N = 295). Variables were ranked for importance using model specific metrics and marginal effects of predictors were aggregated and assessed at the global level. Model selection was carried out using the Mathews Correlation Coefficient (MCC) and model performance was quantified in the validation set using MCC, the area under the precision/recall curve (AUPRC) and accuracy. RESULTS: Random Forest (MCC = 0.28, AUPRC = 0.76) and Adaptive Regression Splines (MCC = 0.29, AUPRC = 0.77) were the best performing models and showed the smallest reduction in performance between the training and validation dataset. EQ5D index, the 10-item personality dimensions, HbA1c, Depression and Anxiety t-scores, age and Body Mass Index were consistently amongst the most powerful predictors in classifying painful vs painless DPN. CONCLUSIONS: Machine learning models trained on large cross-sectional cohorts were able to accurately classify painful or painless DPN on an independent population-based dataset. Painful DPN is associated with more depression, anxiety and certain personality traits. It is also associated with poorer self-reported quality of life, younger age, poor glucose control and high Body Mass Index (BMI). The models showed good performance in realistic conditions in the presence of missing values and noisy datasets. These models can be used either in the clinical context to assist patient stratification based on the risk of painful DPN or return broad risk categories based on user input. Model's performance and calibration suggest that in both cases they could potentially improve diagnosis and outcomes by changing modifiable factors like BMI and HbA1c control and institute earlier preventive or supportive measures like psychological interventions.

Global, regional, and national burden of mortality associated with short-term temperature variability from 2000-19: a three-stage modelling study.

BACKGROUND: Increased mortality risk is associated with short-term temperature variability. However, to our knowledge, there has been no comprehensive assessment of the temperature variability-related mortality burden worldwide. In this study, using data from the MCC Collaborative Research Network, we first explored the association between temperature variability and mortality across 43 countries or regions. Then, to provide a more comprehensive picture of the global burden of mortality associated with temperature variability, global gridded temperature data with a resolution of 0·5°â€ˆ× 0·5° were used to assess the temperature variability-related mortality burden at the global, regional, and national levels. Furthermore, temporal trends in temperature variability-related mortality burden were also explored from 2000-19. METHODS: In this modelling study, we applied a three-stage meta-analytical approach to assess the global temperature variability-related mortality burden at a spatial resolution of 0·5°â€ˆ× 0·5° from 2000-19. Temperature variability was calculated as the SD of the average of the same and previous days' minimum and maximum temperatures. We first obtained location-specific temperature variability related-mortality associations based on a daily time series of 750 locations from the Multi-country Multi-city Collaborative Research Network. We subsequently constructed a multivariable meta-regression model with five predictors to estimate grid-specific temperature variability related-mortality associations across the globe. Finally, percentage excess in mortality and excess mortality rate were calculated to quantify the temperature variability-related mortality burden and to further explore its temporal trend over two decades. FINDINGS: An increasing trend in temperature variability was identified at the global level from 2000 to 2019. Globally, 1 753 392 deaths (95% CI 1 159 901-2 357 718) were associated with temperature variability per year, accounting for 3·4% (2·2-4·6) of all deaths. Most of Asia, Australia, and New Zealand were observed to have a higher percentage excess in mortality than the global mean. Globally, the percentage excess in mortality increased by about 4·6% (3·7-5·3) per decade. The largest increase occurred in Australia and New Zealand (7·3%, 95% CI 4·3-10·4), followed by Europe (4·4%, 2·2-5·6) and Africa (3·3, 1·9-4·6). INTERPRETATION: Globally, a substantial mortality burden was associated with temperature variability, showing geographical heterogeneity and a slightly increasing temporal trend. Our findings could assist in raising public awareness and improving the understanding of the health impacts of temperature variability. FUNDING: Australian Research Council, Australian National Health & Medical Research Council.