gcarrasco-publications

PUBLICATIONS

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.

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

Objective To evaluate lag-response associations and effect modifications of exposure to floods with risks of all cause, cardiovascular, and respiratory mortality on a global scale.

Design Time series study.

Setting 761 communities in 35 countries or territories with at least one flood event during the study period.

Participants Multi-Country Multi-City Collaborative Research Network database, Australian Cause of Death Unit Record File, New Zealand Integrated Data Infrastructure, and the International Network for the Demographic Evaluation of Populations and their Health Network database.

Main outcome measures The main outcome was daily counts of deaths. An estimation for the lag-response association between flood and daily mortality risk was modelled, and the relative risks over the lag period were cumulated to calculate overall effects. Attributable fractions of mortality due to floods were further calculated. A quasi-Poisson model with a distributed lag non-linear function was used to examine how daily death risk was associated with flooded days in each community, and then the community specific associations were pooled using random effects multivariate meta-analyses. Flooded days were defined as days from the start date to the end date of flood events.

Results A total of 47.6 million all cause deaths, 11.1 million cardiovascular deaths, and 4.9 million respiratory deaths were analysed. Over the 761 communities, mortality risks increased and persisted for up to 60 days (50 days for cardiovascular mortality) after a flooded day. The cumulative relative risks for all cause, cardiovascular, and respiratory mortality were 1.021 (95% confidence interval 1.006 to 1.036), 1.026 (1.005 to 1.047), and 1.049 (1.008 to 1.092), respectively. The associations varied across countries or territories and regions. The flood-mortality associations appeared to be modified by climate type and were stronger in low income countries and in populations with a low human development index or high proportion of older people. In communities impacted by flood, up to 0.10% of all cause deaths, 0.18% of cardiovascular deaths, and 0.41% of respiratory deaths were attributed to floods.

Conclusions This study found that the risks of all cause, cardiovascular, and respiratory mortality increased for up to 60 days after exposure to flood and the associations could vary by local climate type, socioeconomic status, and older age.

Mortality risks associated with floods in 761 communities worldwide: time series study

Evidence linking traffic noise to insulin resistance and diabetes is limited and unanswered questions remain regarding the potential effect modification by neighborhood socioeconomic status (nSES). We aimed to assess socioeconomic inequalities in noise exposure, whether road and aircraft noise exposures were associated with insulin resistance or diabetes, and whether nSES modified these relationships. Among the Community of Mine Study in San Diego County, road and aircraft noise exposure at enrollment was calculated based on the static (participant's administrative boundary, and circular buffer around participant homes), and dynamic (mobility data by global positioning system, GPS) spatio-temporal aggregation methods. Associations of noise with insulin resistance (HOMA-IR) or type 2 diabetes (T2DM) were quantified using generalized estimating equation models adjusted for sex, age, ethnicity, individual income, and air pollution (nitrogen dioxide) exposure. Additive interaction between noise and nSES was assessed. Among 573 participants (mean age 58.7 y), participants living in low nSES were exposed to higher levels of aircraft and road noise using noise level at the census tract, circular buffer, or Kernel Density Estimation (KDE) of GPS data. Participants exposed to road noise greater or equal to the median (53 dB(A)) at the census tract and living in low nSES had an increased level of insulin resistance (β = 0.15, 95%CI: −0.04, 0.34) and higher odds of T2DM (Odds Ratio = 2.34, 95%CI: 1.12, 4.90). A positive additive interaction was found as participants living in low nSES had higher odds of T2DM. The impact of noise exposure on insulin resistance and T2DM differs substantially by nSES. Public health benefits of reducing exposure to road or aircraft noise would be larger in individuals living in low nSES.

Aircraft and road traffic noise, insulin resistance, and diabetes: The role of neighborhood socioeconomic status in San Diego County

Climate change constitutes an unprecedented challenge for public health and one of its main direct effects are extreme temperatures. It varies between intra-urban areas and this difference is called surface urban heat island (SUHI) effect. We aimed to assess SUHI distribution among socioeconomic levels in Lima, Peru by conducting a cross-sectional study at the block-level. The mean land surface temperature (LST) from 2017 to 2021 were estimated using the TIRS sensor (Landsat-8 satellite [0.5 km scale]) and extracted to block level. SUHI was calculated based on the difference on mean LST values (2017–2021) per block and the lowest LST registered in a block. Socioeconomic data were obtained from the 2017 Peruvian census. A principal component analysis was performed to construct a socioeconomic index and a mixture analysis based on quantile g-computation was conducted to estimate the joint and specific effects of socioeconomic variables on SUHI. A total of 69 618 blocks were included in the analysis. In the Metropolitan Lima area, the mean SUHI estimation per block was 6.44 (SD = 1.44) Celsius degrees. We found that blocks with high socioeconomic status (SES) showed a decreased exposure to SUHI, compared to those blocks where the low SES were predominant (p-value < 0.001) and that there is a significant SUHI exposure variation (p-value < 0.001) between predominant ethnicities per block (Non-White, Afro-American, and White ethnicities). The mixture analysis showed that the overall mixture effect estimates on SUHI was −1.01 (effect on SUHI of increasing simultaneously every socioeconomic variable by one quantile). Our study highlighted that populations with low SES are more likely to be exposed to higher levels of SUHI compared to those who have a higher SES and illustrates the importance to consider SES inequalities when designing urban adaptation strategies aiming at reducing exposure to SUHI.

Disproportionate exposure to surface-urban heat islands across vulnerable populations in Lima city, Peru

Disease control programs are needed to identify the breeding sites of mosquitoes, which transmit malaria and other diseases, in order to target interventions and identify environmental risk factors. The increasing availability of very-high-resolution drone data provides new opportunities to find and characterize these vector breeding sites. Within this study, drone images from two malaria-endemic regions in Burkina Faso and Côte d’Ivoire were assembled and labeled using open-source tools. We developed and applied a workflow using region-of-interest-based and deep learning methods to identify land cover types associated with vector breeding sites from very- high-resolution natural color imagery. Analysis methods were assessed using cross-validation and achieved maximum Dice coefficients of 0.68 and 0.75 for vegetated and non-vegetated water bodies, respectively. This classifier consistently identified the presence of other land cover types associated with the breeding sites, obtaining Dice coefficients of 0.88 for tillage and crops, 0.87 for buildings and 0.71 for roads. This study establishes a framework for developing deep learning approaches to identify vector breeding sites and highlights the need to evaluate how results will be used by control programs.

Mapping Malaria Vector Habitats in West Africa: Drone Imagery and Deep Learning Analysis for Targeted Vector Surveillance

BACKGROUND: Epidemiological evidence on the health risks of sulfur dioxide (SO2) is more limited compared with other pollutants, and doubts remain on several aspects, such as the form of the exposure–response relationship, the potential role of copollutants, as well as the actual risk at low concentrations and possible temporal variation in risks.
OBJECTIVES: Our aim was to assess the short-term association between exposure to SO2 and daily mortality in a large multilocation data set, using advanced study designs and statistical techniques.
METHODS: The analysis included 43,729,018 deaths that occurred in 399 cities within 23 countries between 1980 and 2018. A two-stage design was applied to assess the association between the daily concentration of SO2 and mortality counts, including first-stage time-series regressions and second-stage multilevel random-effect meta-analyses. Secondary analyses assessed the exposure–response shape and the lag structure using spline terms and distributed lag models, respectively, and temporal variations in risk using a longitudinal meta-regression. Bi-pollutant models were applied to examine confounding effects of particulate matter with an aerodynamic diameter of ≤10 lm (PM10) and 2:5 lm (PM2:5), ozone, nitrogen dioxide, and carbon monoxide. Associations were reported as relative risks (RRs) and fractions of excess deaths.
RESULTS: The average daily concentration of SO2 across the 399 cities was 11:7lg=m3, with 4.7% of days above the World Health Organization (WHO) guideline limit (40lg=m3, 24-h average), although the exceedances occurred predominantly in specific locations. Exposure levels decreased considerably during the study period, from an average concentration of 19:0 lg=m3 in 1980–1989 to 6:3 lg=m3 in 2010–2018. For all locations com- bined, a 10-lg=m3 increase in daily SO2 was associated with an RR of mortality of 1.0045 [95% confidence interval (CI): 1.0019, 1.0070], with the risk being stable over time but with substantial between-country heterogeneity. Short-term exposure to SO2 was associated with an excess mortality fraction of 0.50% [95% empirical CI (eCI): 0.42%, 0.57%] in the 399 cities, although decreasing from 0.74% (0.61%, 0.85%) in 1980–1989 to 0.37% (0.27%, 0.47%) in 2010–2018. There was some evidence of nonlinearity, with a steep exposure–response relationship at low concentrations and the risk attenuating at higher levels. The relevant lag window was 0–3 d. Significant positive associations remained after controlling for other pollutants.
DISCUSSION: The analysis revealed independent mortality risks associated with short-term exposure to SO2, with no evidence of a threshold. Levels below the current WHO guidelines for 24-h averages were still associated with substantial excess mortality, indicating the potential benefits of stricter air quality standards.

Short-Term Association between Sulfur Dioxide and Mortality: A Multicountry Analysis in 399 Cities

Background
Globally, access to life-saving vaccines has improved considerably in the past 5 decades. However, progress has started to slow down and even reverse in recent years. Understanding subnational heterogeneities in essential child immunization will be critical for closing the global vaccination gap.

Methods and findings
We use vaccination information for over 220,000 children across 1,366 administrative regions in 43 low- and middle-income countries (LMICs) from the most recent Demographic and Health Surveys. We estimate essential immunization coverage at the national and subnational levels and quantify socioeconomic inequalities in such coverage using adjusted concentration indices. Within- and between-country variations are summarized via the Theil index. We use local indicator of spatial association (LISA) statistics to identify clusters of administrative regions with high or low values. Finally, we estimate the number of missed vaccinations among children aged 15 to 35 months across all 43 countries and the types of vaccines most often missed. We show that national-level vaccination rates can conceal wide subnational heterogeneities. Large gaps in child immunization are found across West and Central Africa and in South Asia, particularly in regions of Angola, Chad, Nigeria, Guinea, and Afghanistan, where less than 10% of children are fully immunized. Furthermore, children living in these countries consistently lack all 4 basic vaccines included in the WHO’s recommended schedule for young children. Across most countries, children from poorer households are less likely to be fully immunized. The main limitations include subnational estimates based on large administrative divisions for some countries and different periods of survey data collection.

Conclusions
The identified heterogeneities in essential childhood immunization, especially given that some regions consistently are underserved for all basic vaccines, can be used to inform the design and implementation of localized intervention programs aimed at eliminating child suffering and deaths from existing and novel vaccine-preventable diseases.

Essential childhood immunization in 43 low- and middle-income countries: Analysis of spatial trends and socioeconomic inequalities in vaccine coverage

In recent years, global health security has been threatened by the geographical expansion of vector-borne infectious diseases such as malaria, dengue, yellow fever, Zika and chikungunya. For a range of these vector-borne diseases, an increase in residual (exophagic) transmission together with ecological heterogeneity in everything from weather to local human migration and housing to mosquito species’ behaviours presents many challenges to effective mosquito control. The novel use of drones (or uncrewed aerial vehicles) may play a major role in the success of mosquito surveillance and control programmes in the coming decades since the global landscape of mosquito-borne diseases and disease dynamics fluctuates frequently and there could be serious public health consequences if the issues of insecticide resistance and outdoor transmission are not adequately addressed. For controlling both aquatic and adult stages, for several years now remote sensing data have been used together with predictive modelling for risk, incidence and detection of transmission hot spots and landscape profiles in relation to mosquito-borne pathogens. The field of drone-based remote sensing is under continuous change due to new technology development, operation regulations and innovative applications. In this review we outline the opportunities and challenges for integrating drones into vector surveillance (i.e. identification of breeding sites or mapping micro-environmental composition) and control strategies (i.e. applying larval source management activities or deploying genetically modified agents) across the mosquito life-cycle. We present a five-step systematic environmental mapping strategy that we recommend be undertaken in locations where a drone is expected to be used, outline the key considerations for incorporating drone or other Earth Observation data into vector surveillance and provide two case studies of the advantages of using drones equipped with multispectral cameras. In conclusion, recent developments mean that drones can be effective for accurately conducting surveillance, assessing habitat suitability for larval and/or adult mosquitoes and implementing interventions. In addition, we briefly discuss the need to consider permissions, costs, safety/privacy perceptions and community acceptance for deploying drone activities.

The use of drones for mosquito surveillance and control

Lactoferrin (LF) has in vitro antiviral activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study aimed to determine the effect of bovine lactoferrin (bLF) in the prevention of SARS-CoV-2 infection in health care personnel. A randomized, double-blinded, placebo-controlled clinical trial was conducted in two tertiary hospitals that provide care to patients with SARS-CoV-2 infection in Lima, Peru. Daily supplementation with 600 mg of enteral bLF versus placebo for 90 days was compared. Participants were weekly screened for symptoms suggestive of SARS-CoV-2 infection and molecular testing was performed on suspected episodes. A serological test was obtained from all participants at the end of the intervention. The main outcome included symptomatic and asymptomatic cases. A sub-analysis explored the time to symptomatic infection. Secondary outcomes were the severity, frequency, and duration of symptomatic infection. The study was prematurely cancelled due to the availability of vaccines against SARS-CoV-2 in Peru. 209 participants were enrolled and randomized, 104 received bLF and 105 placebo. SARS-CoV-2 infection occurred in 11 (10.6%) participants assigned to bLF and in 9 (8.6%) participants assigned to placebo without significant differences (Incidence Rate Ratio = 1.23, 95%CI 0.51–3.06, p-value = 0.64). There was no significant effect of bLF on time to symptomatic infection (Hazard Ratio = 1.61, 95%CI 0.62–4.19, p-value = 0.3). There were no significant differences in secondary outcomes. A significant effect of bLF in preventing SARS-CoV-2 infection was not proven. Further studies are needed to assess the effect of bLF supplementation on SARS-CoV-2 infection.

Bovine lactoferrin for the prevention of COVID-19 infection in health care personnel: a double-blinded randomized clinical trial (LF-COVID)

In the Amazon Region of Peru, occupational activities are important drivers of human mobility and may increase the individual risk of being infected while contributing to increasing malaria community-level transmission. Even though out-of-village working activities and other mobility patterns have been identified as determinants of malaria transmission, no studies have quantified the effect of out-of-village working activities on recent malaria exposure and proposed plausible intervention scenarios. Using two population-based cross-sectional studies in the Loreto Department in Peru, and the parametric g-formula method, we simulated various hypothetical scenarios intervening in out-of-village working activities to reflect their potential health benefits. This study estimated that the standardized mean outcome (malaria seroprevalence) in the unexposed population (no out-of-village workers) was 44.6% (95% CI: 41.7%–47.5%) and 66.7% (95% CI: 61.6%–71.8%) in the exposed population resulting in a risk difference of 22.1% (95% CI: 16.3%–27.9%). However, heterogeneous patterns in the effects of interest were observed between peri-urban and rural areas (Cochran’s Q test = 15.5, p < 0.001). Heterogeneous patterns were also observed in scenarios of increased prevalence of out-of-village working activities and restriction scenarios by gender (male vs. female) and age (18 and under vs. 19 and older) that inform possible occupational interventions targetting population subgroups. The findings of this study support the hypothesis that targeting out-of-village workers will considerably benefit current malaria elimination strategies in the Amazon Region. Particularly, males and adult populations that carried out out-of-village working activities in rural areas contribute the most to the malaria seropositivity (recent exposure to the parasite) in the Peruvian Amazon.

Effect of out-of-village working activities on recent malaria exposure in the Peruvian Amazon using parametric g-formula

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.

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

The impact of human population movement (HPM) on the epidemiology of vector-borne diseases, such as malaria, has been described. However, there are limited data on the use of new technologies for the study of HPM in endemic areas with difficult access such as the Amazon. In this study conducted in rural Peruvian Amazon, we used self-reported travel surveys and GPS trackers coupled with a Bayesian spatial model to quantify the role of HPM on malaria risk. By using a densely sampled population cohort, this study highlighted the elevated malaria transmission in a riverine community of the Peruvian Amazon. We also found that the high connectivity between Amazon communities for reasons such as work, trading or family plausibly sustains such transmission levels. Finally, by using multiple human mobility metrics including GPS trackers, and adapted causal inference methods we identified for the first time the effect of human mobility patterns on malaria risk in rural Peruvian Amazon. This study provides evidence of the causal effect of HPM on malaria that may help to adapt current malaria control programmes in the Amazon.

Quantifying the effect of human population mobility on malaria risk in the Peruvian Amazon

Background
Missed Opportunities for Vaccination (MOV) represent a major risk in the re-emergence of immunopreventable diseases. However, in the region, there are few published studies on MOVs using national databases such as demographic and health surveys (DHS). This study aims to describe the frequency and trends of MOVs for the first dose of vaccines against the leading causes of infant morbidity and mortality, their complete vaccination coverage, and trends in socioeconomic inequalities at the national and departmental levels for an 11-years period.
Methods
Using DHS data from an 11-year period (2010–2020), we calculated frequencies and trends in MOVs of vaccines for the leading causes of child morbidity and mortality, estimated inequalities in MOVs using the Slope Inequality Index (SII) and conducted a spatial autocorrelation test to identify clusters of higher or lower inequality in MOVs at the national level.
Findings
We found that, at the national level, greater inequality was concentrated in the wealthiest categories of each socioeconomic variable. We identified that departments with higher poverty rates concentrated higher levels of inequality in the MOVs in the lowest strata of the socioeconomic variables. In addition, we found that some departments with similar geographic and socioeconomic characteristics had spatially correlated levels of inequality on MOVs.
Interpretation
These findings can help to identify the heterogeneity that exists in the distribution of MOVs among departments and socioeconomic strata, which would help to prioritize specific areas and subpopulations for national immunization strategies.

Missed opportunities for vaccination in Peru 2010–2020: A study of socioeconomic inequalities

Understanding the reservoir and infectivity of Plasmodium gametocytes to vector mosquitoes is crucial to align strategies aimed at malaria transmission elimination. Yet, experimental information is scarce regarding the infectivity of Plasmodium vivax for mosquitoes in diverse epidemiological settings where the proportion of asymptomatically infected individuals varies at a microgeographic scale. We measured the transmissibility of clinical and subclinical P. vivax malaria parasite carriers to the major mosquito vector in the Amazon Basin, Nyssorhynchus darlingi (formerly Anopheles). A total of 105 participants with natural P. vivax malaria infection were recruited from a cohort study in Loreto Department, Peruvian Amazon. Four of 18 asymptomatic individuals with P. vivax positivity by blood smear infected colony-grown Ny. darlingi (22%), with 2.6% (19 of 728) mosquitoes infected. In contrast, 77% (44/57) of symptomatic participants were infectious to mosquitoes with 51% (890 of 1,753) mosquitoes infected. Infection intensity was greater in symptomatic infections (mean, 17.8 oocysts/mosquito) compared with asymptomatic infections (mean, 0.28 oocysts/mosquito), attributed to parasitemia/gametocytemia level. Paired experiments (N = 27) using direct skin-feeding assays and direct membrane mosquito-feeding assays showed that infectivity to mosquitoes was similar for both methods. Longitudinal studies with longer follow-up of symptomatic and asymptomatic parasite infections are needed to determine the natural variations of disease transmissibility.

Insights into Plasmodium vivax Asymptomatic Malaria Infections and Direct Skin-Feeding Assays to Assess Onward Malaria Transmission in the Amazon

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.

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

Background
The landscape of malaria transmission in the Peruvian Amazon is temporally and spatially heterogeneous, presenting different micro-geographies with particular epidemiologies. Most cases are asymptomatic and escape routine malaria surveillance based on light microscopy (LM). Following the implementation of control programs in this region, new approaches to stratify transmission and direct efforts at an individual and community level are needed. Antibody responses to serological exposure markers (SEM) to Plasmodium vivax have proven diagnostic performance to identify people exposed in the previous 9 months.

Methodology
We measured antibody responses against 8 SEM to identify recently exposed people and determine the transmission dynamics of P. vivax in peri-urban (Iquitos) and riverine (Mazán) communities of Loreto, communities that have seen significant recent reductions in malaria transmission. Socio-demographic, geo-reference, LM and qPCR diagnosis data were collected from two cross-sectional surveys. Spatial and multilevel analyses were implemented to describe the distribution of seropositive cases and the risk factors associated with exposure to P. vivax.

Principal findings
Low local transmission was detected by qPCR in both Iquitos (5.3%) and Mazán (2.7%); however, seroprevalence indicated a higher level of (past) exposure to P. vivax in Mazán (56.5%) than Iquitos (38.2%). Age and being male were factors associated with high odds of being seropositive in both sites. Higher antibody levels were found in individuals >15 years old. The persistence of long-lived antibodies in these individuals could overestimate the detection of recent exposure. Antibody levels in younger populations (<15 years old) could be a better indicator of recent exposure to P. vivax.

Conclusions
The large number of current and past infections detected by SEMs allows for detailed local epidemiological analyses, in contrast to data from qPCR prevalence surveys which did not produce statistically significant associations. Serological surveillance will be increasingly important in the Peruvian Amazon as malaria transmission is reduced by continued control and elimination efforts.

Malaria transmission structure in the Peruvian Amazon through antibody signatures to Plasmodium vivax

Background
In Brazil, malaria is concentrated in the Amazon Basin, where more than 99% of the annual cases are reported. The main goal of this study was to investigate the population structure and genetic association of the biting behavior of Nyssorhynchus (also known as Anopheles) darlingi, the major malaria vector in the Amazon region of Brazil, using low-coverage genomic sequencing data.

Methods
Samples were collected in the municipality of Mâncio Lima, Acre state, Brazil between 2016 and 2017. Different approaches using genotype imputation and no gene imputation for data treatment and low-coverage sequencing genotyping were performed. After the samples were genotyped, population stratification analysis was performed.

Results
Weak but statistically significant stratification signatures were identified between subpopulations separated by distances of approximately 2–3 km. Genome-wide association studies (GWAS) were performed to compare indoor/outdoor biting behavior and blood-seeking at dusk/dawn. A statistically significant association was observed between biting behavior and single nucleotide polymorphism (SNP) markers adjacent to the gene associated with cytochrome P450 (CYP) 4H14, which is associated with insecticide resistance. A statistically significant association between blood-seeking periodicity and SNP markers adjacent to genes associated with the circadian cycle was also observed.

Conclusion
The data presented here suggest that low-coverage whole-genome sequencing with adequate processing is a powerful tool to genetically characterize vector populations at a microgeographic scale in malaria transmission areas, as well as for use in GWAS. Female mosquitoes entering houses to take a blood meal may be related to a specific CYP4H14 allele, and female timing of blood-seeking is related to circadian rhythm genes.

Nyssorhynchus darlingi genome-wide studies related to microgeographic dispersion and blood-seeking behavior

Studies have investigated the effects of heat and temperature variability (TV) on mortality. However, few assessed whether TV modifies the heat-mortality association. Data on daily temperature and mortality in the warm season were collected from 717 locations across 36 countries. TV was calculated as the standard deviation of the average of the same and previous days’ minimum and maximum temperatures. We used location-specific quasi-Poisson regression models with an interaction term between the cross-basis term for mean temperature and quartiles of TV to obtain heat-mortality associations under each quartile of TV, and then pooled estimates at the country, regional, and global levels. Results show the increased risk in heat-related mortality with increments in TV, accounting for 0.70% (95% confidence interval [CI]: −0.33 to 1.69), 1.34% (95% CI: −0.14 to 2.73), 1.99% (95% CI: 0.29–3.57), and 2.73% (95% CI: 0.76–4.50) of total deaths for Q1–Q4 (first quartile–fourth quartile) of TV. The modification effects of TV varied geographically. Central Europe had the highest attributable fractions (AFs), corresponding to 7.68% (95% CI: 5.25–9.89) of total deaths for Q4 of TV, while the lowest AFs were observed in North America, with the values for Q4 of 1.74% (95% CI: −0.09 to 3.39). TV had a significant modification effect on the heat-mortality association, causing a higher heat-related mortality burden with increments of TV. Implementing targeted strategies against heat exposure and fluctuant temperatures simultaneously would benefit public health.

Fluctuating temperature modifies heat-mortality association around the globe

Exposure to air pollution disproportionately affects racial/ethnic minorities that could contribute to health inequalities including metabolic disorders. However, most existing studies used a static assessment of air pollution exposure (mostly using the residential address) and do not account for activity space when modelling exposure to air pollution. The aim of this study is to understand how exposure to air pollution impacts metabolic disorders biomarkers, how this effect differs according to ethnicity, and for the first time compare these findings with two methods of exposure assessment: dynamic and static measures.

Methods
Among the Community of Mine study, a cross-sectional study conducted in San Diego County, insulin resistance, diabetes, hypertension, obesity, dyslipidemia, and metabolic syndrome (MetS) were assessed. Exposure to air pollution (PM2.5, NO2, traffic) was calculated using static measures around the home, and dynamic measures of mobility derived from Global Positioning Systems (GPS) traces using kernel density estimators to account for exposure variability across space and time. Associations of air pollution with metabolic disorders were quantified using generalized estimating equation models to account for the clustered nature of the data.

Results
Among 552 participants (mean age 58.7 years, 42% Hispanic/Latino), Hispanics/Latinos had a higher exposure to PM2.5 compared to non-Hispanics using static measures. In contrast, Hispanics/Latinos had less exposure to PM2.5 using dynamic measures. For all participants, higher dynamic exposure to PM2.5 and NO2 was associated with increased insulin resistance and cholesterol levels, and increased risk of obesity, dyslipidemia and MetS (RR 1.17, 95% CI: 1.07–1.28; RR 1.21, 95% CI: 1.12–1.30, respectively). The association between dynamic PM2.5 exposure and MetS differed by Hispanic/Latino ethnicity.

Conclusion
These results highlight the importance of considering people's daily mobility in assessing the impact of air pollution on health.

Air pollution and metabolic disorders: Dynamic versus static measures of exposure among Hispanics/Latinos and non-Hispanics

Objectives. To determine the effect of heat waves on emergency department (ED) visits for individuals experiencing homelessness and explore vulnerability factors.

Methods. We used a unique highly detailed data set on sociodemographics of ED visits in San Diego, California, 2012 to 2019. We applied a time-stratified case–crossover design to study the association between various heat wave definitions and ED visits. We compared associations with a similar population not experiencing homelessness using coarsened exact matching.

Results. Of the 24 688 individuals identified as experiencing homelessness who visited an ED, most were younger than 65 years (94%) and of non-Hispanic ethnicity (84%), and 14% indicated the need for a psychiatric consultation. Results indicated a positive association, with the strongest risk of ED visits during daytime (e.g., 99th percentile, 2 days) heat waves (odds ratio = 1.29; 95% confidence interval = 1.02, 1.64). Patients experiencing homelessness who were younger or elderly and who required a psychiatric consultation were particularly vulnerable to heat waves. Odds of ED visits were higher for individuals experiencing homelessness after matching to nonhomeless individuals based on age, gender, and race/ethnicity.

Conclusions. It is important to prioritize individuals experiencing homelessness in heat action plans and consider vulnerability factors to reduce their burden.

Heat Waves and Emergency Department Visits Among the Homeless, San Diego, 2012–2019

Land-use practices such as agriculture can impact mosquito vector breeding ecology, resulting in changes in disease transmission. The typical breeding habitats of Africa’s second most important malaria vector Anopheles funestus are large, semipermanent water bodies, which make them potential candidates for targeted larval source management. This is a technical workflow for the integration of drone surveys and mosquito larval sampling, designed for a case study aiming to characterise An. funestus breeding sites near two villages in an agricultural setting in Côte d’Ivoire. Using satellite remote sensing data, we developed an environmentally and spatially representative sampling frame and conducted paired mosquito larvae and drone mapping surveys from June to August 2021. To categorise the drone imagery, we also developed a land cover classification scheme with classes relative to An. funestus breeding ecology. We sampled 189 potential breeding habitats, of which 119 (63%) were positive for the Anopheles genus and nine (4.8%) were positive for An. funestus. We mapped 30.42 km2 of the region of interest including all water bodies which were sampled for larvae. These data can be used to inform targeted vector control efforts, although its generalisability over a large region is limited by the fine-scale nature of this study area. This paper develops protocols for integrating drone surveys and statistically rigorous entomological sampling, which can be adjusted to collect data on vector breeding habitats in other ecological contexts. Further research using data collected in this study can enable the development of deep-learning algorithms for identifying An. funestus breeding habitats across rural agricultural landscapes in Côte d’Ivoire and the analysis of risk factors for these sites.

Technical Workflow Development for Integrating Drone Surveys and Entomological Sampling to Characterise Aquatic Larval Habitats of Anopheles funestus in Agricultural Landscapes in Côte d’Ivoire

There is conflicting evidence on the influence of weather on COVID-19 transmission. Our aim is to estimate weather-dependent signatures in the early phase of the pandemic, while controlling for socio-economic factors and non-pharmaceutical interventions. We identify a modest non-linear association between mean temperature and the effective reproduction number (Re) in 409 cities in 26 countries, with a decrease of 0.087 (95% CI: 0.025; 0.148) for a 10 °C increase. Early interventions have a greater effect on Re with a decrease of 0.285 (95% CI 0.223; 0.347) for a 5th - 95th percentile increase in the government response index. The variation in the effective reproduction number explained by government interventions is 6 times greater than for mean temperature. We find little evidence of meteorological conditions having influenced the early stages of local epidemics and conclude that population behaviour and government interventions are more important drivers of transmission.

A cross-sectional analysis of meteorological factors and SARS-CoV-2 transmission in 409 cities across 26 countries

Successful malaria control interventions, mostly based on the training of health workers, distribution of insecticide-treated nets, and spraying, decrease malaria incidence; however, when these interventions are interrupted, a resurgence may occur. In the Peruvian Amazon, after discontinuing the control activities implemented by the PAMAFRO project (2006–2010)-a Global Fund-sponsored project for the strengthening of malaria control and surveillance in multiple countries in Latin America– malaria cases re-emerged dramatically. In parallel, meteorological factors determine the conditions suitable for the development, reproduction, and survival of mosquito vectors and parasites. This study hypothesized that interruption of malaria interventions may have modified the meteorological-malaria relationships over time (i.e., temporal changes in the dose-response between meteorological variables and malaria incidence). In this panel data analysis, we assessed the extent that relationships between meteorological variables and malaria changed temporally using data of monthly malaria incidence due to Plasmodium vivax or P. falciparum in Loreto, Peru (2000–2017). Generalized additive models were used to explore how the effects of meteorological variables changed in magnitude before, during, and after the PAMAFRO intervention. We found that once the PAMAFRO intervention had been interrupted, the estimated effects (dose-response) of meteorological variables on incidence rates decreased for both malaria parasite species. However, these fitted effect estimates did not reach their baseline levels (before the PAMAFRO period); variations of time-varying slopes between 0.45 and 2.07 times were observed after the PAMAFRO intervention. We also reported significant heterogeneity in the geographical distributions of malaria, parasite species, and meteorological variables. High malaria transmission occurred consistently in the northwestern provinces of Loreto Department. Since the end of the PAMAFRO period, a higher effect of precipitation and actual evapotranspiration was described on P. falciparum compared to P. vivax. The effect of temperature on malaria was greater over a shorter time (1-month lag or less), compared with precipitation and actual evapotranspiration (12-month lag). These findings demonstrate the importance of sustained malaria control efforts since interruption may enhance the links between meteorological factors and malaria. Our results also emphasize the importance of considering the time-varying effect of meteorological factors on malaria incidence to tailor control interventions, especially to better manage the current and future climate change crisis.

Time-Varying Effects of Meteorological Variables on Malaria Epidemiology in the Context of Interrupted Control Efforts in the Amazon Rainforest, 2000–2017

Background:
Minimum mortality temperature (MMT) is an important indicator to assess the temperature-mortality association, indicating long-term adaptation to local climate. Limited evidence about the geographical variability of the MMT is available at a global scale.

Methods:
We collected data from 658 communities in 43 countries under different climates. We estimated temperature-mortality associations to derive the MMT for each community using Poisson regression with distributed lag nonlinear models. We investigated the variation in MMT by climatic zone using a mixed-effects meta-analysis and explored the association with climatic and socioeconomic indicators.

Results:
The geographical distribution of MMTs varied considerably by country between 14.2 and 31.1 °C decreasing by latitude. For climatic zones, the MMTs increased from alpine (13.0 °C) to continental (19.3 °C), temperate (21.7 °C), arid (24.5 °C), and tropical (26.5 °C). The MMT percentiles (MMTPs) corresponding to the MMTs decreased from temperate (79.5th) to continental (75.4th), arid (68.0th), tropical (58.5th), and alpine (41.4th). The MMTs indreased by 0.8 °C for a 1 °C rise in a community’s annual mean temperature, and by 1 °C for a 1 °C rise in its SD. While the MMTP decreased by 0.3 centile points for a 1 °C rise in a community’s annual mean temperature and by 1.3 for a 1 °C rise in its SD.

Conclusions:
The geographical distribution of the MMTs and MMTPs is driven mainly by the mean annual temperature, which seems to be a valuable indicator of overall adaptation across populations. Our results suggest that populations have adapted to the average temperature, although there is still more room for adaptation.

Geographical Variations of the Minimum Mortality Temperature at a Global Scale: A Multicountry Study

Background
Many regions of the world are now facing more frequent and unprecedentedly large wildfires. However, the association between wildfire-related PM2·5 and mortality has not been well characterised. We aimed to comprehensively assess the association between short-term exposure to wildfire-related PM2·5 and mortality across various regions of the world.
Methods
For this time series study, data on daily counts of deaths for all causes, cardiovascular causes, and respiratory causes were collected from 749 cities in 43 countries and regions during 2000–16. Daily concentrations of wildfire-related PM2·5 were estimated using the three-dimensional chemical transport model GEOS-Chem at a 0·25° × 0·25° resolution. The association between wildfire-related PM2·5 exposure and mortality was examined using a quasi-Poisson time series model in each city considering both the current-day and lag effects, and the effect estimates were then pooled using a random-effects meta-analysis. Based on these pooled effect estimates, the population attributable fraction and relative risk (RR) of annual mortality due to acute wildfire-related PM2·5 exposure was calculated.
Findings
65·6 million all-cause deaths, 15·1 million cardiovascular deaths, and 6·8 million respiratory deaths were included in our analyses. The pooled RRs of mortality associated with each 10 μg/m3 increase in the 3-day moving average (lag 0–2 days) of wildfire-related PM2·5 exposure were 1·019 (95% CI 1·016–1·022) for all-cause mortality, 1·017 (1·012–1·021) for cardiovascular mortality, and 1·019 (1·013–1·025) for respiratory mortality. Overall, 0·62% (95% CI 0·48–0·75) of all-cause deaths, 0·55% (0·43–0·67) of cardiovascular deaths, and 0·64% (0·50–0·78) of respiratory deaths were annually attributable to the acute impacts of wildfire-related PM2·5 exposure during the study period.
Interpretation
Short-term exposure to wildfire-related PM2·5 was associated with increased risk of mortality. Urgent action is needed to reduce health risks from the increasing wildfires.

Mortality risk attributable to wildfire-related PM2·5 pollution: a global time series study in 749 locations