key: cord-0727516-0wspp086 authors: Li, Adam Y.; Hannah, Theodore C; Durbin, John R; Dreher, Nickolas; McAuley, Fiona M; Marayati, Naoum Fares; Spiera, Zachary; Ali, Muhammad; Gometz, Alex; Kostman, JT; Choudhri, Tanvir F. title: Multivariate Analysis of Black Race and Environmental Temperature on COVID-19 In the US date: 2020-06-20 journal: Am J Med Sci DOI: 10.1016/j.amjms.2020.06.015 sha: 4d38d1fdbd9b5ab145dcb4d76d0655196a2b333a doc_id: 727516 cord_uid: 0wspp086 BACKGROUND: There has been much interest in environmental temperature and race as modulators of Coronavirus disease-19 (COVID-19) infection and mortality. However, in the United States race and temperature correlate with various other social determinants of health, comorbidities, and environmental influences that could be responsible for noted effects. This study investigates the independent effects of race and environmental temperature on COVID-19 incidence and mortality in United States counties. METHODS: Data on COVID-19 and risk factors in all United States counties was collected. 661 counties with at least 50 COVID-19 cases and 217 with at least 10 deaths were included in analyses. Upper and lower quartiles for cases/100,000 people and halves for deaths/100,000 people were compared with t-tests. Adjusted linear and logistic regression analyses were performed to evaluate the independent effects of race and environmental temperature. RESULTS: Multivariate regression analyses demonstrated Black race is a risk factor for increased COVID-19 cases (OR=1.22, 95% CI: 1.09-1.40, P=0.001) and deaths independent of comorbidities, poverty, access to health care, and other risk factors. Higher environmental temperature independently reduced caseload (OR=0.81, 95% CI: 0.71-0.91, P=0.0009), but not deaths. CONCLUSIONS: Higher environmental temperatures correlated with reduced COVID-19 cases, but this benefit does not yet appear in mortality models. Black race was an independent risk factor for increased COVID-19 cases and deaths. Thus, many proposed mechanisms through which Black race might increase risk for COVID-19, such as socioeconomic and healthcare-related predispositions, are inadequate in explaining the full magnitude of this health disparity. Since the first case of coronavirus disease 2019 in December 2019, COVID-19 has spread rapidly around the world. 1 including patient demographics, 2 social determinants of health, 3 environmental variables, [4] [5] [6] housing, 7 and underlying health conditions. [8] [9] [10] In particular, there have been reports of an environmental temperature effect on COVID-19 cases, 4-6 however the intersection of environmental temperature and socioeconomic variables has not been fully elucidated in the United States. Additionally, numerous headlines have brought race to the forefront as a potentially significant risk factor for increased COVID-19 incidence. 11 For example, African Americans have accounted for all but three COVID-19 related deaths in St. Louis, Missouri and three quarters of those in Milwaukee, Wisconsin. 12 However, the effects of race are often intermixed with the effects of other socioeconomic factors. In this study, using multivariate regression analyses controlling for demographics, underlying health conditions, social determinants of health, environmental variables, and social distancing adherence, we seek to illuminate the independent effects of both race and environmental temperature on the incidence of COVID-19 cases and mortality at the US county level. COVID-19 confirmed case number and death number through April 14, 2020 were obtained for each US county from the Center for Systems Science and Engineering (CSSE) Coronavirus Resource Center at Johns Hopkins University. This data source contains COVID-19 case number and deaths data from all 3,143 US counties in 50 states and the District of Columbia (D.C.). Cases/100,000 people and deaths/100,000 people for each US county were calculated using CSSE and census data. Counties were excluded from the analyses if they had fewer than 50 cases or their first case occurred fewer than 3 weeks prior to the end of the study. With these constraints, 661 counties from 50 states and Washington D.C were included in case analysis. Counties were excluded from death analyses if they reported fewer than 10 deaths or the first death had occurred fewer than 2 weeks prior to the end of the study. The deaths analyses include 217 counties from 37 states and Washington D.C. (Table 1 ). Race demographics for counties was obtained from the County Health Rankings and Roadmaps Program database. Daily environmental temperature data for counties was obtained from the National Oceanic and Atmospheric Administration. County temperature was calculated using mean temperature for a period starting 10 days before the first confirmed county case and through the most current date (April 14, 2020). Unacast social distancing data was obtained through a research agreement with the company. Data from mobile phones was used to track average changes in distance travelled within counties after the first reported COVID-19 case. Change in average distance travelled was used as an indicator of social distancing. Data for other potential risk factors used in analysis was obtained through publicly available data sources. (Table 2) . T-tests were used to compare differences between covariates for the highest quartile and lowest quartile cases/100,000 and the top 50% versus the bottom 50% deaths/100,000. Since modeling complex phenomenon without controlling for a variety of confounding variables can lead to skewed results, sequential regression modelling was used to demonstrate the effect of race and temperature on cases and deaths/100,000. The stepwise approach to regression helps mitigate excessive confounding risk through sensitivity analysis. Logistic regression was used to compare counties in the highest quartile for cases/100,000 to those in the lowest. There were an insufficient number of counties meeting the deaths/100,000 requirements to run logistic regression, thus linear regression was used for deaths/100,000 analyses. Results of linear regression for cases per/100,000 including all 661 counties are provided for direct comparison between cases and deaths. The sequential model had four parts. Each model included all variables from the previous models. Model 1 was a univariate analysis. Model 2 added the following county macroeconomic and COVID-19-specific variables: population density, GDP/capita, COVID-19 tests/100,000 people (state level data), COVID-19 cases/100,000 people (deaths analysis only), a marker for which State the county is in, average percent reduction in cellphone movement from day of first confirmed case to end of study, the percent of population living in overcrowded housing, and the number of days since the first confirmed COVID-19 case. Model population over 65, proportion of Black residents, percent of the population that is female, percent of population living in rural areas, the Food Environment Index (a measure of accessibility and affordability of healthy food), the rate of violent crime/100,000 people, the average environmental temperature from 10 days before the first case to the end of the study, air quality measured as the average annual ambient concentrations of atmospheric particulate matter under 2.5 micrometers (PM2.5), percent of the population considered to be in fair or poor health, and the percent of the population living in poverty. Model 4 was the full model and included medical comorbidities and access to health care variables in addition to all variables in Model 3. Specifically these new variables were: diabetes, obesity, physical inactivity, excessive drinking, and smoking were all reported as percent of the population; liver disease, hypertension, coronary heart disease, and chronic respiratory disease were reported as mortality/100,000 people; and the patient to primary care physician ratio, the percent of people sleeping fewer than seven hours/night, the percent of the population without health insurance and the percent of the population who received the flu vaccine were also included. All statistical analyses were performed with Prism 8.0 (GraphPad Software, San Diego, CA). For all analyses, α = 0.05. Our study contained data from 3,143 US counties across all 50 states and the District of Columbia (D.C). 661 US counties (21.0%) from 50 states and D.C. were included in case analysis ( Figure 1A) . 217 US counties (6.9%) from 37 states and D.C. were included in death analysis ( Figure 1B ). There were significant differences in demographics between the two cohorts for cases/100,000 Table 4 ). The United States has become a major epicenter of the coronavirus pandemic, now reporting the greatest number of total cases and deaths worldwide. By conducting sequential multivariate analyses of variables that span social, structural, and environmental spheres, our study aimed to more completely characterize the effects of Black race and environmental temperature on the incidence and mortality of COVID-19 in the U.S. Our results suggest that higher environmental temperatures may lead to decreased incidence, however a similar effect was not seen for mortality. Black race was an independent risk factor for increased COVID-19 cases and mortality. As noted previously, Black Americans have accounted for all but three COVID-19 related deaths in St. Louis, Missouri and three quarters of those in Milwaukee, Wisconsin. 12 Further, recent analyses of patients hospitalized for COVID-19 in 99 United States counties spanning 14 states have demonstrated significantly increased hospitalization rates amongst Black Americans. 11 Our study corroborates and expands on these findings to suggest that Black Americans may be at significantly higher risk of COVID-19 infection and mortality nationwide. Studies on COVID-19 incidence amongst Black Americans have pointed to the fact that, compared to the general population of the United States, Black Americans have been shown to have increased rates of smoking tobacco use and chronic medical comorbidities such as cardiovascular disease, diabetes mellitus, hypertension, obesity, and chronic respiratory disease, all of which convey increased risk of adverse COVID-19 outcomes. [8] [9] [10] [11] Recent work has also cited the many inequities in structural variables, which now manifest in the disproportionate number of Black Americans suffering from poverty, environmental hazards, overcrowded housing, and decreased access to healthcare. 11 However, after controlling for these factors in this study, the effect of Black race on number of COVID-19 cases and deaths persisted, indicating other modulating factors must be considered. One potential domain includes additional sociocultural variables that have become relevant within the context of the pandemic. For example, African Americans make up a large percentage of the healthcare, transportation, government, and food supply industries, job sectors that have been deemed "essential" services in light of the SARS-CoV-2 pandemic. 13 Despite the majority of Americans currently living in counties with some form of social distancing mandate, fewer than one in five Black Americans have a job with flexibility to work from home, compared with more than a third of white and Asian American workers. 14 Thus, racial differences in social distancing may exist and could be contributing to increased risk of COVID-19 for Black Americans. Another potential consideration is that non-Hispanic Black populations in the United States suffer from vitamin D deficiency at higher rates than other racial groups. In a recent comprehensive analysis of serum 25-hydroxyvitamin D (25(OH)D) levels collected from more than 26,000 adults in the United States, 71.9% of non-Hispanic Black individuals were found to suffer from vitamin D deficiency and 22.6% had vitamin D insufficiency. Of the many sociodemographic, behavioral, and clinical variables analyzed, being a racial or ethnic minority was the strongest predictor of vitamin D deficiency, with non-Hispanic Black populations suffering from vitamin D deficiency at rates 24 times more than those seen in the white population. 15 Vitamin D has been suggested as a potentially mitigative factor in the COVID-19 pandemic due to its important modulatory effect on immune response. 16, 17 Furthermore, vitamin D deficiency has been implicated in numerous adverse health conditions such as acute respiratory syndromes, tuberculosis, and cardiovascular disease. 18 However, no clinical evidence has shown vitamin D supplements are beneficial in prevention or treatment of COVID-19, as no randomized trials on vitamin D supplementations have been performed. 19 While some registered open label trials are being performed, no results have been reported, and it is still unclear whether vitamin D is beneficial. 19 Environmental variables such as temperature have been studied significantly around the world, with varying results. Studies have positive correlations, negative correlations, or no correlations between rates of COVID-19 and increased temperature. [4] [5] [6] Here, we demonstrate an independent effect of environmental temperature that results in reduced COVID-19 cases, but not mortality, across US counties. If warmer temperatures do, in fact, play a role in mitigating disease spread, it is reasonable to expect a potential seasonal trend in global cases. As environmental temperatures begin to warm in the Northern hemisphere during the summer months, we may begin to see a decrease in disease burden; however, we should prepare to see a resurgence in COVID-19 incidence when environmental temperatures again decrease in the fall and winter months, as has been exhibited in previous pandemics like the Spanish flu. 20 Likewise, countries in the Southern hemisphere should be prepared to experience seasonal trends in the opposite direction. The mechanism through which environmental temperature may execute its mitigative effects remains unclear, though previous studies have shown other viruses including SARS-CoV and influenza have similar associations with temperature. [21] [22] [23] [24] Warmer temperatures have been associated with decreased transmission for these viruses. [21] [22] [23] [24] Some suggested mechanisms include increased stability of viral particles and smaller droplets at colder climates. For droplets, colder climates are associated with lower humidity, causing smaller droplets to remain airborne for extended periods of time and increase spread. 24 Lastly, while warmer temperature may or may not affect rates of COVID-19 cases and deaths, temperature alone will likely not be sufficient in halting further spread of COVID-19. 25 In univariate analysis, chronic respiratory disease did not significantly affect COVID-19 case and death numbers, and increased air pollution was associated with decreased COVID-19 cases. These results were unexpected given that both variables have been shown to increase risk for COVID-19. 9, 26 However, univariate analysis may be affected by confounding variables and may not represent true effects. For example, environmental temperature in this study was not significant in univariate analysis but was significant in multivariate sequential modeling. Additionally, air pollution data in this study is from 2014 and may not reflect current air pollution. Given increased social distancing, many counties with historically high air pollution levels currently have significantly decreased pollution levels, which may have decreased pollution related COVID-19 cases. 27, 28 Limitations This study is limited by excluding counties with insufficient COVID-19 case or death data. This limitation is unavoidable and waiting for sufficient data may deter public health response. Also, the COVID-19 pandemic is still in progress and counties are at different stages. Thus, case and mortality statistics for the counties provide a snapshot of the current state, but may not, and likely will not, be reflective of the ultimate case and death tolls in these counties. Moreover, most counties are still in the early stages of accumulating deaths which axiomatically lags behind the number of cases. Since we relied heavily on publicly available sources to create our database, we used sources that may be outdated or inconsistent with the statistics during the past few months. Social distancing was measured through changes in distance travelled within US counties. While this correlates with COVID-19 cases and is an indicator of social distancing, mobility data does not capture whether people are near each other and truly social distancing. we report as granularly as practical to account for this. This study evaluated the independent effects of Black race and environmental temperature on the incidence of cases and mortality of COVID-19 at the US county level. In multivariate regression analyses that controlled for county demographics, socioeconomic factors, and medical comorbidities, counties with higher average environmental temperatures were associated with decreased COVID-19 cases, but not deaths. Black race was significantly associated with both increased cases and increased deaths. This suggests that many of the proposed mechanisms through which Black race might increase risk for COVID-19, such as socioeconomic and healthcare-related predispositions, are inadequate in explaining the full magnitude of this health disparity. Florida 35 11 Ohio 27 10 Georgia 44 12 Oklahoma 9 3 Hawaii 2 0 Oregon 6 1 Idaho 4 0 Pennsylvania 31 13 Illinois 15 5 Rhode Island 4 0 Indiana 30 8 South Carolina 17 2 Iowa 8 1 South Dakota 2 0 Kansas 5 2 Tennessee 14 4 Kentucky 6 1 Texas 32 10 Louisiana 34 19 Utah 6 0 Maine 3 1 Vermont 4 1 Maryland 15 8 Virginia 23 2 Massachusetts 12 11 Washington 14 8 Michigan 23 6 West Virginia 4 0 Minnesota 6 1 Wisconsin 10 2 Mississippi 16 0 Wyoming 2 0 Missouri 10 3 Novel coronavirus: where we are and what we know Evolving epidemiology and transmission dynamics of coronavirus disease 2019 outside Hubei province, China: a descriptive and modelling study Knowledge, attitudes, and practices towards COVID-19 among Chinese residents during the rapid rise period of the COVID-19 outbreak: a quick online crosssectional survey The correlation between the spread of COVID-19 infections and weather variables in 30 Chinese provinces and the impact of Chinese government mitigation plans Temperature significantly changes COVID-19 transmission in (sub)tropical cities of Brazil Effect of temperature on the infectivity of COVID-19 COVID-19: a potential public health problem for homeless populations Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection? Covid-19: risk factors for severe disease and death Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study COVID-19 and African Americans COVID-19 exacerbating inequalities in the US How can we save black and brown lives during a pandemic? Data from past studies can point the way. Science Magazine Website Not everybody can work from home: Black and Hispanic workers are much less likely to be able to telework. Economic Policy Institute Website Vitamin D deficiency and insufficiency among US adults: prevalence, predictors and clinical implications Latitude Dependence of the COVID-19 Mortality Rate-A Possible Relationship to Vitamin D Deficiency? Evidence that Vitamin D Supplementation Could Reduce Risk of Influenza and COVID-19 Infections and Deaths Vitamin D deficiency: a worldwide problem with health consequences Vitamin D: A rapid review of the evidence for treatment or prevention in COVID-19. Centre for Evidence-Based Medicine Website The seasonality of pandemic and non-pandemic influenzas: the roles of solar radiation and vitamin D The Effects of Temperature and Relative Humidity on the Viability of the SARS Coronavirus Temperature Decreases Spread Parameters of the New Covid-19 Case Dynamics Environmental factors on the SARS epidemic: air temperature, passage of time and multiplicative effect of hospital infection Decline in temperature and humidity increases the occurrence of influenza in cold climate Seasonality of SARS-CoV-2: Will COVID-19 go away on its own in warmer weather? Center for Communicable Disease Dynamics Website Exposure to air pollution and COVID-19 mortality in the United States COVID-19 as a factor influencing air pollution? But right now it has some of the cleanest of any major city. CNN Website We thank Lauren Spinazze from the Unacast Group for providing us with data from their Social Distancing Scoreboard. We would also like to thank Aman M. Choudhri for his help with data acquisition.