PYPE: A pipeline for phenome-wide association and Mendelian randomization in investigator-driven biobank scale analysis.

Phenome-wide association studies (PheWASs) serve as a way of documenting the relationship between genotypes and multiple phenotypes, helping to uncover unexplored genotype-phenotype associations (known as pleiotropy). Secondly, Mendelian randomization (MR) can be harnessed to make causal statements about a pair of phenotypes by comparing their genetic architecture. Thus, approaches that automate both PheWASs and MR can enhance biobank-scale analyses, circumventing the need for multiple tools by providing a comprehensive, end-to-end tool to drive scientific discovery. To this end, we present PYPE, a Python pipeline for running, visualizing, and interpreting PheWASs. PYPE utilizes input genotype or phenotype files to automatically estimate associations between the chosen independent variables and phenotypes. PYPE can also produce a variety of visualizations and can be used to identify nearby genes and functional consequences of significant associations. Finally, PYPE can identify possible causal relationships between phenotypes using MR under a variety of causal effect modeling scenarios.

Exposome-wide association study of cognition among older adults in the National Health and Nutrition Examination Survey.

Cognitive impairment among older adults is a growing public health challenge and environmental chemicals may be modifiable risk factors. A wide array of chemicals has not yet been tested for association with cognition in an environment-wide association framework. In the US National Health and Nutrition Examination Survey (NHANES) 1999-2000 and 2011-2014 cross-sectional cycles, cognition was assessed using the Digit Symbol Substitution Test (DSST, scores 0-117) among participants aged 60 years and older. Concentrations of environmental chemicals measured in blood or urine were log 2 transformed and standardized. Chemicals with at least 50% of measures above the lower limit of detection were included (n chemicals =147, n classes =14). We tested for associations between chemical concentrations and cognition using parallel survey-weighted multivariable linear regression models adjusted for age, sex, race/ethnicity, education, smoking status, fish consumption, cycle year, urinary creatinine, and cotinine. Participants with at least one chemical measurement (n=4,982) were mean age 69.8 years, 55.0% female, 78.2% non-Hispanic White, and 77.0% at least high school educated. The mean DSST score was 50.4 (standard deviation (SD)=17.4). In adjusted analyses, 5 of 147 exposures were associated with DSST at p-value<0.01. Notably, a SD increase in log 2 -scaled cotinine concentration was associated with 2.71 points lower DSST score (95% CI -3.69, -1.73). A SD increase in log 2 -scaled urinary tungsten concentration was associated with 1.34 points lower DSST score (95% CI -2.11, -0.56). Exposure to environmental chemicals, particularly heavy metals and tobacco smoke, may be modifiable factors for cognition among older adults.

Projected Changes in Statin and Antihypertensive Therapy Eligibility With the AHA PREVENT Cardiovascular Risk Equations.

Importance: Since 2013, the American College of Cardiology (ACC) and American Heart Association (AHA) have recommended the pooled cohort equations (PCEs) for estimating the 10-year risk of atherosclerotic cardiovascular disease (ASCVD). An AHA scientific advisory group recently developed the Predicting Risk of cardiovascular disease EVENTs (PREVENT) equations, which incorporated kidney measures, removed race as an input, and improved calibration in contemporary populations. PREVENT is known to produce ASCVD risk predictions that are lower than those produced by the PCEs, but the potential clinical implications have not been quantified. Objective: To estimate the number of US adults who would experience changes in risk categorization, treatment eligibility, or clinical outcomes when applying PREVENT equations to existing ACC and AHA guidelines. Design, Setting, and Participants: Nationally representative cross-sectional sample of 7765 US adults aged 30 to 79 years who participated in the National Health and Nutrition Examination Surveys of 2011 to March 2020, which had response rates ranging from 47% to 70%. Main Outcomes and Measures: Differences in predicted 10-year ASCVD risk, ACC and AHA risk categorization, eligibility for statin or antihypertensive therapy, and projected occurrences of myocardial infarction or stroke. Results: In a nationally representative sample of 7765 US adults aged 30 to 79 years (median age, 53 years; 51.3% women), it was estimated that using PREVENT equations would reclassify approximately half of US adults to lower ACC and AHA risk categories (53.0% [95% CI, 51.2%-54.8%]) and very few US adults to higher risk categories (0.41% [95% CI, 0.25%-0.62%]). The number of US adults receiving or recommended for preventive treatment would decrease by an estimated 14.3 million (95% CI, 12.6 million-15.9 million) for statin therapy and 2.62 million (95% CI, 2.02 million-3.21 million) for antihypertensive therapy. The study estimated that, over 10 years, these decreases in treatment eligibility could result in 107 000 additional occurrences of myocardial infarction or stroke. Eligibility changes would affect twice as many men as women and a greater proportion of Black adults than White adults. Conclusion and Relevance: By assigning lower ASCVD risk predictions, application of the PREVENT equations to existing treatment thresholds could reduce eligibility for statin and antihypertensive therapy among 15.8 million US adults.

Epidemic outcomes following government responses to COVID-19: Insights from nearly 100,000 models.

Government responses to COVID-19 are among the most globally impactful events of the 21st century. The extent to which responses-such as school closures-were associated with changes in COVID-19 outcomes remains unsettled. Multiverse analyses offer a systematic approach to testing a large range of models. We used daily data on 16 government responses in 181 countries in 2020-2021, and 4 outcomes-cases, infections, COVID-19 deaths, and all-cause excess deaths-to construct 99,736 analytic models. Among those, 42% suggest outcomes improved following more stringent responses ("helpful"). No subanalysis (e.g. limited to cases as outcome) demonstrated a preponderance of helpful or unhelpful associations. Among the 14 associations with P values < 1 × 10-30, 5 were helpful and 9 unhelpful. In summary, we find no patterns in the overall set of models that suggests a clear relationship between COVID-19 government responses and outcomes. Strong claims about government responses' impacts on COVID-19 may lack empirical support.

Panel stacking is a threat to consensus statement validity.

Consensus statements can be very influential in medicine and public health. Some of these statements use systematic evidence synthesis but others fail on this front. Many consensus statements use panels of experts to deduce perceived consensus through Delphi processes. We argue that stacking of panel members toward one particular position or narrative is a major threat, especially in absence of systematic evidence review. Stacking may involve financial conflicts of interest, but nonfinancial conflicts of strong advocacy can also cause major bias. Given their emerging importance, we describe here how such consensus statements may be misleading, by analyzing in depth a recent high-impact Delphi consensus statement on COVID-19 recommendations as a case example. We demonstrate that many of the selected panel members and at least 35% of the core panel members had advocated toward COVID-19 elimination (Zero-COVID) during the pandemic and were leading members of aggressive advocacy groups. These advocacy conflicts were not declared in the Delphi consensus publication, with rare exceptions. Therefore, we propose that consensus statements should always require rigorous evidence synthesis and maximal transparency on potential biases toward advocacy or lobbyist groups to be valid. While advocacy can have many important functions, its biased impact on consensus panels should be carefully avoided.

Implications of Race Adjustment in Lung-Function Equations.

BACKGROUND: Adjustment for race is discouraged in lung-function testing, but the implications of adopting race-neutral equations have not been comprehensively quantified. METHODS: We obtained longitudinal data from 369,077 participants in the National Health and Nutrition Examination Survey, U.K. Biobank, the Multi-Ethnic Study of Atherosclerosis, and the Organ Procurement and Transplantation Network. Using these data, we compared the race-based 2012 Global Lung Function Initiative (GLI-2012) equations with race-neutral equations introduced in 2022 (GLI-Global). Evaluated outcomes included national projections of clinical, occupational, and financial reclassifications; individual lung-allocation scores for transplantation priority; and concordance statistics (C statistics) for clinical prediction tasks. RESULTS: Among the 249 million persons in the United States between 6 and 79 years of age who are able to produce high-quality spirometric results, the use of GLI-Global equations may reclassify ventilatory impairment for 12.5 million persons, medical impairment ratings for 8.16 million, occupational eligibility for 2.28 million, grading of chronic obstructive pulmonary disease for 2.05 million, and military disability compensation for 413,000. These potential changes differed according to race; for example, classifications of nonobstructive ventilatory impairment may change dramatically, increasing 141% (95% confidence interval [CI], 113 to 169) among Black persons and decreasing 69% (95% CI, 63 to 74) among White persons. Annual disability payments may increase by more than $1 billion among Black veterans and decrease by $0.5 billion among White veterans. GLI-2012 and GLI-Global equations had similar discriminative accuracy with regard to respiratory symptoms, health care utilization, new-onset disease, death from any cause, death related to respiratory disease, and death among persons on a transplant waiting list, with differences in C statistics ranging from -0.008 to 0.011. CONCLUSIONS: The use of race-based and race-neutral equations generated similarly accurate predictions of respiratory outcomes but assigned different disease classifications, occupational eligibility, and disability compensation for millions of persons, with effects diverging according to race. (Funded by the National Heart Lung and Blood Institute and the National Institute of Environmental Health Sciences.).

Specification curve analysis to identify heterogeneity in risk factors for dementia: findings from the UK Biobank.

BACKGROUND: In 2020, the Lancet Commission identified 12 risk factors as priorities for prevention of dementia, and other studies identified APOE e4/e4 genotype and family history of Alzheimer's disease strongly associated with dementia outcomes; however, it is unclear how robust these relationships are across dementia subtypes and analytic scenarios. Specification curve analysis (SCA) is a new tool to probe how plausible analytical scenarios influence outcomes. METHODS: We evaluated the heterogeneity of odds ratios for 12 risk factors reported from the Lancet 2020 report and two additional strong associated non-modifiable factors (APOE e4/e4 genotype and family history of Alzheimer's disease) with dementia outcomes across 450,707 UK Biobank participants using SCA with 5357 specifications across dementia subtypes (outcomes) and analytic models (e.g., standard demographic covariates such as age or sex and/or 14 correlated risk factors). RESULTS: SCA revealed variable dementia risks by subtype and age, with associations for TBI and APOE e4/e4 robust to model specification; in contrast, diabetes showed fluctuating links with dementia subtypes. We found that unattributed dementia participants had similar risk factor profiles to participants with defined subtypes. CONCLUSIONS: We observed heterogeneity in the risk of dementia, and estimates of risk were influenced by the inclusion of a combination of other modifiable risk factors; non-modifiable demographic factors had a minimal role in analytic heterogeneity. Future studies should report multiple plausible analytic scenarios to test the robustness of their association. Considering these combinations of risk factors could be advantageous for the clinical development and evaluation of novel screening models for different types of dementia.

Environmental chemical-wide associations with immune biomarkers in US adults: A cross-sectional analysis.

Environmental chemical exposures influence immune system functions, and humans are exposed to a wide range of chemicals, termed the chemical "exposome". A comprehensive, discovery analysis of the associations of multiple chemical families with immune biomarkers is needed. In this study, we tested the associations between environmental chemical concentrations and immune biomarkers. We analyzed the United States cross-sectional National Health and Nutrition Examination Survey (NHANES, 1999-2018). Chemical biomarker concentrations were measured in blood or urine (196 chemicals, 17 chemical families). Immune biomarkers included counts of lymphocytes, neutrophils, monocytes, basophils, eosinophils, red blood cells, white blood cells, and mean corpuscular volume. We conducted separate survey-weighted, multivariable linear regressions of each log2-transformed chemical and immune measure, adjusted for relevant covariates. We accounted for multiple comparisons using a false discovery rate (FDR). Among 45,528 adult participants, the mean age was 45.7 years, 51.4% were female, and 69.3% were Non-Hispanic White. 71 (36.2%) chemicals were associated with at least one of the eight immune biomarkers. The most chemical associations (FDR<0.05) were observed with mean corpuscular volume (36 chemicals) and red blood cell counts (35 chemicals). For example, a doubling in the concentration of cotinine was associated with 0.16 fL (95% CI: 0.15, 0.17; FDR<0.001) increased mean corpuscular volume, and a doubling in the concentration of blood lead was associated with 61,736 increased red blood cells per µL (95% CI: 54,335, 69,138; FDR<0.001). A wide variety of chemicals, such as metals and smoking-related compounds, were highly associated with immune system biomarkers. This environmental chemical-wide association study identified chemicals from multiple families for further toxicological, immunologic, and epidemiological investigation.

Heterogeneity in elevated glucose and A1C as predictors of the prediabetes to diabetes transition: Framingham Heart Study, Multi-Ethnic Study on Atherosclerosis, Jackson Heart Study, and Atherosclerosis Risk In Communities.

Introduction: There are a number of glycemic definitions for prediabetes; however, the heterogeneity in diabetes transition rates from prediabetes across different glycemic definitions in major US cohorts has been unexplored. We estimate the variability in risk and relative risk of adiposity based on diagnostic criteria like fasting glucose and hemoglobin A1C% (HA1C%). Research Design and Methods: We estimated transition rate from prediabetes, as defined by fasting glucose between 100-125 and/or 110-125 mg/dL, and HA1C% between 5.7-6.5% in participant data from the Framingham Heart Study, Multi-Ethnic Study on Atherosclerosis, Atherosclerosis Risk in Communities, and the Jackson Heart Study. We estimated the heterogeneity and prediction interval across cohorts, stratifying by age, sex, and body mass index. For individuals who were prediabetic, we estimated the relative risk for obesity, blood pressure, education, age, and sex for diabetes. Results: There is substantial heterogeneity in diabetes transition rates across cohorts and prediabetes definitions with large prediction intervals. We observed the highest range of rates in individuals with fasting glucose of 110-125 mg/dL ranging from 2-18 per 100 person-years. Across different cohorts, the association obesity or hypertension in the progression to diabetes was consistent, yet it varied in magnitude. We provide a database of transition rates across subgroups and cohorts for comparison in future studies. Conclusion: The absolute transition rate from prediabetes to diabetes significantly depends on cohort and prediabetes definitions.

Assessing the genetic contribution of cumulative behavioral factors associated with longitudinal type 2 diabetes risk highlights adiposity and the brain-metabolic axis.

While genetic factors, behavior, and environmental exposures form a complex web of interrelated associations in type 2 diabetes (T2D), their interaction is poorly understood. Here, using data from ~500K participants of the UK Biobank, we identify the genetic determinants of a "polyexposure risk score" (PXS) a new risk factor that consists of an accumulation of 25 associated individual-level behaviors and environmental risk factors that predict longitudinal T2D incidence. PXS-T2D had a non-zero heritability (h2 = 0.18) extensive shared genetic architecture with established clinical and biological determinants of T2D, most prominently with body mass index (genetic correlation [rg] = 0.57) and Homeostatic Model Assessment for Insulin Resistance (rg = 0.51). Genetic loci associated with PXS-T2D were enriched for expression in the brain. Biobank scale data with genetic information illuminates how complex and cumulative exposures and behaviors as a whole impact T2D risk but whose biology have been elusive in genome-wide studies of T2D.

Grilling the data: application of specification curve analysis to red meat and all-cause mortality.

OBJECTIVES: To present an application of specification curve analysis-a novel analytic method that involves defining and implementing all plausible and valid analytic approaches for addressing a research question-to nutritional epidemiology. STUDY DESIGN AND SETTING: We reviewed all observational studies addressing the effect of red meat on all-cause mortality, sourced from a published systematic review, and documented variations in analytic methods (eg, choice of model, covariates, etc.). We enumerated all defensible combinations of analytic choices to produce a comprehensive list of all the ways in which the data may reasonably be analyzed. We applied specification curve analysis to data from National Health and Nutrition Examination Survey 2007 to 2014 to investigate the effect of unprocessed red meat on all-cause mortality. The specification curve analysis used a random sample of all reasonable analytic specifications we sourced from primary studies. RESULTS: Among 15 publications reporting on 24 cohorts included in the systematic review on red meat and all-cause mortality, we identified 70 unique analytic methods, each including different analytic models, covariates, and operationalizations of red meat (eg, continuous vs quantiles). We applied specification curve analysis to National Health and Nutrition Examination Survey, including 10,661 participants. Our specification curve analysis included 1208 unique analytic specifications, of which 435 (36.0%) yielded a hazard ratio equal to or more than 1 for the effect of red meat on all-cause mortality and 773 (64.0%) less than 1. The specification curve analysis yielded a median hazard ratio of 0.94 (interquartile range: 0.83-1.05). Forty-eight specifications (3.97%) were statistically significant, 40 of which indicated unprocessed red meat to reduce all-cause mortality and eight of which indicated red meat to increase mortality. CONCLUSION: We show that the application of specification curve analysis to nutritional epidemiology is feasible and presents an innovative solution to analytic flexibility.

Decoding the exposome: data science methodologies and implications in exposome-wide association studies (ExWASs).

This paper explores the exposome concept and its role in elucidating the interplay between environmental exposures and human health. We introduce two key concepts critical for exposomics research. Firstly, we discuss the joint impact of genetics and environment on phenotypes, emphasizing the variance attributable to shared and nonshared environmental factors, underscoring the complexity of quantifying the exposome's influence on health outcomes. Secondly, we introduce the importance of advanced data-driven methods in large cohort studies for exposomic measurements. Here, we introduce the exposome-wide association study (ExWAS), an approach designed for systematic discovery of relationships between phenotypes and various exposures, identifying significant associations while controlling for multiple comparisons. We advocate for the standardized use of the term "exposome-wide association study, ExWAS," to facilitate clear communication and literature retrieval in this field. The paper aims to guide future health researchers in understanding and evaluating exposomic studies. Our discussion extends to emerging topics, such as FAIR Data Principles, biobanked healthcare datasets, and the functional exposome, outlining the future directions in exposomic research. This abstract provides a succinct overview of our comprehensive approach to understanding the complex dynamics of the exposome and its significant implications for human health.

Pediatric and Young Adult Household Transmission of the Initial Waves of SARS-CoV-2 in the United States: Administrative Claims Study.

BACKGROUND: The correlates responsible for the temporal changes of intrahousehold SARS-CoV-2 transmission in the United States have been understudied mainly due to a lack of available surveillance data. Specifically, early analyses of SARS-CoV-2 household secondary attack rates (SARs) were small in sample size and conducted cross-sectionally at single time points. From these limited data, it has been difficult to assess the role that different risk factors have had on intrahousehold disease transmission in different stages of the ongoing COVID-19 pandemic, particularly in children and youth. OBJECTIVE: This study aimed to estimate the transmission dynamic and infectivity of SARS-CoV-2 among pediatric and young adult index cases (age 0 to 25 years) in the United States through the initial waves of the pandemic. METHODS: Using administrative claims, we analyzed 19 million SARS-CoV-2 test records between January 2020 and February 2021. We identified 36,241 households with pediatric index cases and calculated household SARs utilizing complete case information. Using a retrospective cohort design, we estimated the household SARS-CoV-2 transmission between 4 index age groups (0 to 4 years, 5 to 11 years, 12 to 17 years, and 18 to 25 years) while adjusting for sex, family size, quarter of first SARS-CoV-2 positive record, and residential regions of the index cases. RESULTS: After filtering all household records for greater than one member in a household and missing information, only 36,241 (0.85%) of 4,270,130 households with a pediatric case remained in the analysis. Index cases aged between 0 and 17 years were a minority of the total index cases (n=11,484, 11%). The overall SAR of SARS-CoV-2 was 23.04% (95% CI 21.88-24.19). As a comparison, the SAR for all ages (0 to 65+ years) was 32.4% (95% CI 32.1-32.8), higher than the SAR for the population between 0 and 25 years of age. The highest SAR of 38.3% was observed in April 2020 (95% CI 31.6-45), while the lowest SAR of 15.6% was observed in September 2020 (95% CI 13.9-17.3). It consistently decreased from 32% to 21.1% as the age of index groups increased. In a multiple logistic regression analysis, we found that the youngest pediatric age group (0 to 4 years) had 1.69 times (95% CI 1.42-2.00) the odds of SARS-CoV-2 transmission to any family members when compared with the oldest group (18 to 25 years). Family size was significantly associated with household viral transmission (odds ratio 2.66, 95% CI 2.58-2.74). CONCLUSIONS: Using retrospective claims data, the pediatric index transmission of SARS-CoV-2 during the initial waves of the COVID-19 pandemic in the United States was associated with location and family characteristics. Pediatric SAR (0 to 25 years) was less than the SAR for all age other groups. Less than 1% (n=36,241) of all household data were retained in the retrospective study for complete case analysis, perhaps biasing our findings. We have provided measures of baseline household pediatric transmission for tracking and comparing the infectivity of later SARS-CoV-2 variants.

Prediction and stratification of longitudinal risk for chronic obstructive pulmonary disease across smoking behaviors.

Smoking is the leading risk factor for chronic obstructive pulmonary disease (COPD) worldwide, yet many people who never smoke develop COPD. We perform a longitudinal analysis of COPD in the UK Biobank to derive and validate the Socioeconomic and Environmental Risk Score which captures additive and cumulative environmental, behavioral, and socioeconomic exposure risks beyond tobacco smoking. The Socioeconomic and Environmental Risk Score is more predictive of COPD than smoking status and pack-years. Individuals in the highest decile of the risk score have a greater risk for incident COPD compared to the remaining population. Never smokers in the highest decile of exposure risk are more likely to develop COPD than previous and current smokers in the lowest decile. In general, the prediction accuracy of the Social and Environmental Risk Score is lower in non-European populations. While smoking status is often considered in screening COPD, our finding highlights the importance of other non-smoking environmental and socioeconomic variables.

DNA characterization reveals potential operon-unit packaging of extracellular vesicle cargo from a gut bacterial symbiont.

Extracellular vesicles (EVs) are lipid bilayer-bound entities secreted by cells across all domains of life, known to contain a range of components, including protein complexes, RNA, and DNA. Recent studies on microbial extracellular vesicles indicate that these virus-sized nanoparticles, 40-90nm in diameter, readily cross the epithelial barrier and reach systemic circulation, can be detected in tissues throughout the body in mice and that 1mL of plasma from healthy humans contains up to one million bacterial EVs. They have been recently recognized for their biologically functional roles, including modulation of bacterial physiology and host-microbe interactions, hence their gain in the microbiome research community's attention. However, the exact understanding of their functionality is still a subject of active research and debate. Here, we employ long-read DNA sequencing on purified extracellular vesicles from a common mammalian gut symbiont, Parabacteroides goldsteinii, to characterize the genomic component within EV cargos. Our findings challenge the notion of DNA packaging into EVs as a stochastic event. Instead, our data demonstrate that the DNA packaging is non-random. Here, we suggest a novel hypothesis of selective EV-DNA packaging, potentially arranged in operon units, hence providing new insights into our understanding of its genetic makeup and its potential role, underlining the importance of our findings in microbial community dynamics.

Association between timing and consistency of physical activity and type 2 diabetes: a cohort study on participants of the UK Biobank.

AIMS/HYPOTHESIS: We sought to quantify the relationship between morning, afternoon or evening physical activity and consistency (e.g. routine) and risk of type 2 diabetes. METHODS: A cohort of 93,095 UK Biobank participants (mean age 62 years) without a history of type 2 diabetes wore a wrist-worn accelerometer for 1 week. We converted accelerometer information to estimate metabolic equivalent of task (MET), summing MET h of total physical activity completed within three intra-day time segments (morning, afternoon and evening). We quantified physical activity consistency as the SD of participants' daily total physical activity. We ultimately associated each of the following with incident type 2 diabetes: (1) morning, afternoon or evening 'time-segmented' MET h per week; and (2) consistency. We also considered moderate-to-vigorous physical activity (MVPA) and vigorous physical activity (VPA) in association with type 2 diabetes incidence. RESULTS: When considering MET as the physical activity measure, we observed protective associations of morning (HR 0.90 [95% CI 0.86, 0.93], p=7×10-8) and afternoon (HR 0.91 [95% CI 0.87, 0.95], p=1×10-5) but did not have evidence for evening physical activity (HR 0.95 [95% CI 0.90, 1.00], p=0.07) with type 2 diabetes. There was no difference between MET-measured morning and afternoon physical activity. Our substitution model highlighted the importance of adjusting for lifestyle factors (e.g. sleep time and diet); the effect of a substitution between afternoon and evening physical activity was attenuated after adjustment for lifestyle variables. Consistency of MET-measured physical activity was not associated with type 2 diabetes (p=0.07). MVPA and VPA were associated with decreased risk for type 2 diabetes at all times of the day. CONCLUSIONS/INTERPRETATION: Total metabolic equivalents of physical activity in the morning and afternoon had a protective effect on diabetes risk and evening activity was not associated with diabetes. Consistency of physical activity did not play a role in decreasing risk for diabetes. Vigorous activity is associated with lower risk no matter the time of day of activity.

Spatio-temporal interpolation and delineation of extreme heat events in California between 2017 and 2021.

Robust spatio-temporal delineation of extreme climate events and accurate identification of areas that are impacted by an event is a prerequisite for identifying population-level and health-related risks. In prior research, attributes such as temperature and humidity have often been linearly assigned to the population of the study unit from the closest weather station. This could result in inaccurate event delineation and biased assessment of extreme heat exposure. We have developed a spatio-temporal model to dynamically delineate boundaries for Extreme Heat Events (EHE) across space and over time, using a relative measure of Apparent Temperature (AT). Our surface interpolation approach offers a higher spatio-temporal resolution compared to the standard nearest-station (NS) assignment method. We show that the proposed approach can provide at least 80.8 percent improvement in identification of areas and populations impacted by EHEs. This improvement in average adjusts the misclassification of about one million Californians per day of an extreme event, who would be either unidentified or misidentified under EHEs between 2017 and 2021.