Longitudinal Physiological Data from a Wearable Device Identifies SARS-CoV-2 Infection and Symptoms and Predicts COVID-19 Diagnosis =================================================================================================================================== * Robert P. Hirten * Matteo Danieletto * Lewis Tomalin * Katie Hyewon Choi * Micol Zweig * Eddye Golden * Sparshdeep Kaur * Drew Helmus * Anthony Biello * Renata Pyzik * Ismail Nabeel * Alexander Charney * Benjamin Glicksberg * Matthew Levin * David Reich * Dennis Charney * Erwin P Bottinger * Laurie Keefer * Mayte Suarez-Farinas * Girish N. Nadkarni * Zahi A. Fayad ## ABSTRACT **Background** Changes in autonomic nervous system function, characterized by heart rate variability (HRV), have been associated with and observed prior to the clinical identification of infection. We performed an evaluation of this metric collected by wearable devices, to identify and predict Coronavirus disease 2019 (COVID-19) and its related symptoms. **Methods** Health care workers in the Mount Sinai Health System were prospectively followed in an ongoing observational study using the custom Warrior Watch Study App which was downloaded to their smartphones. Participants wore an Apple Watch for the duration of the study measuring HRV throughout the follow up period. Survey’s assessing infection and symptom related questions were obtained daily. **Findings** Using a mixed-effect COSINOR model the mean amplitude of the circadian pattern of the standard deviation of the interbeat interval of normal sinus beats (SDNN), a HRV metric, differed between subjects with and without COVID-19 (p=0.006). The mean amplitude of this circadian pattern differed between individuals during the 7 days before and the 7 days after a COVID-19 diagnosis compared to this metric during uninfected time periods (p=0.01). Significant changes in the mean MESOR and amplitude of the circadian pattern of the SDNN was observed between the first day of reporting a COVID-19 related symptom compared to all other symptom free days (p=0.01). **Interpretation** Longitudinally collected HRV metrics from a commonly worn commercial wearable device (Apple Watch) can identify the diagnosis of COVID-19 and COVID-19 related symptoms. Prior to the diagnosis of COVID-19 by nasal PCR, significant changes in HRV were observed demonstrating its predictive ability to identify COVID-19 infection. **Funding** Support was provided by the Ehrenkranz Lab For Human Resilience, the BioMedical Engineering and Imaging Institute, The Hasso Plattner Institute for Digital Health at Mount Sinai, The Mount Sinai Clinical Intelligence Center and The Dr. Henry D. Janowitz Division of Gastroenterology. ## INTRODUCTION Coronavirus disease 2019 (COVID-19) has resulted in over 41 million infections and more than 1.1 million deaths.1 A prolonged incubation period and variable symptomatology has facilitated disease spread, with approximately 30-45% of individuals having asymptomatic SARS-CoV-2 infections, and testing generally limited to only symptomatic individuals.2-4 Health care workers (HCWs), characterized as any type of worker in a health care system, represent a vulnerable population with a threefold increased risk of infection compared to the general population.5 This increased risk of transmission is important in healthcare settings, where asymptomatic or pre-symptomatic HCWs can shed the virus contributing to transmission within healthcare facilities and their households.6 Digital health technology offers an opportunity to address the limitations of traditional public health strategies aimed at curbing COVID-19 spread.7 Smart phone Apps are effective in using symptoms to identify those possibly infected with SARS-CoV-2, but they rely on ongoing participant compliance and self-reported symptoms.8 Wearable devices are commonly used for remote sensing and provide a means to objectively quantify physiological parameters including heart rate, sleep, activity and measures of autonomic nervous system (ANS) function (e.g., heart rate variability [HRV]).9 The addition of physiological data from wearable devices to symptom tracking Apps has been shown to increase the ability to identify those infected with SARS-CoV-2.10 HRV is a physiological metric providing insight into the interplay between the parasympathetic and sympathetic nervous system which modulate cardiac contractility and cause variability in the beat-to-beat intervals.11 It exhibits a 24 hour circadian pattern with relative sympathetic tone during the day and parasympathetic activity at night.12-14 Changes in this circadian pattern can be leveraged to identify different physiological states. Several studies have demonstrated that lower HRV, indicating increased sympathetic balance, is a reliable predictor of infection onset.15,16 However, HRV and its dynamic changes over time have not been evaluated as a marker or predictor of COVID-19. In response to the COVID-19 pandemic we launched The Warrior Watch Study™, employing a novel smartphone App to remotely enroll and monitor HCWs throughout the Mount Sinai Health System in New York City, a site of initial case surge. This digital platform enables remote survey delivery to Apple iPhones and passive collection of Apple Watch data, including HRV. The aim of this study is to determine if SARS-CoV-2 infections can be identified and predicted prior to a positive test result using the longitudinal changes in HRV metrics derived from the Apple Watch. ## METHODS ### Study Design The primary aim of the study was to determine whether changes in HRV can differentiate participants infected or not infected with SARS-CoV-2. The secondary aim was to see if changes in HRV can predict the development of a SARS-CoV-2 infection prior to diagnosis by a SARS-CoV-2 nasal PCR. Exploratory aims were (1) to determine whether changes in HRV can identify the presence of COVID-19 related symptoms; (2) to determine whether changes in HRV can predict the development of COVID-19 related symptoms; and (3) to evaluate how HRV changed throughout the infection and symptom period. HCWs in the Mount Sinai Health System were enrolled in an ongoing prospective observational cohort study. Eligible participants were ≥18 years of age, current employees in the Mount Sinai Health System, had an iPhone Series 6 or higher, and had or were willing to wear an Apple Watch Series 4 or higher. Participants were excluded if they had an underlying autoimmune disease or were on medications known to interfere with ANS function. A positive COVID-19 diagnosis was defined as a positive SARS-CoV-2 nasal PCR swab reported by the participant. Daily symptoms were collected including fevers/chills, tired/weak, body aches, dry cough, sneezing, runny nose, diarrhea, sore throat, headache, shortness of breath, loss of smell or taste, itchy eyes, none, or other. This study was approved by the Institutional Review Board at The Icahn School of Medicine at Mount Sinai. ### Study Procedures Participants downloaded the custom Warrior Watch App to complete eligibility questionnaires and sign an electronic consent form. Participants completed an App-based baseline assessment collecting demographic information, prior COVID-19 diagnosis history, occupation, and medical history and were then followed prospectively through the App. Daily survey questionnaires captured COVID-19 related symptoms, symptom severity, SARS-CoV-2 nasal PCR results, serum SARS-CoV-2 antibody test results, and daily patient care related exposure (**Supplementary Table 1**). Participants carried out their normal activities throughout the study and were instructed to wear the Apple Watch for a minimum duration of 8 hours per day. View this table: [Supplementary Table 1.](http://medrxiv.org/content/early/2020/11/07/2020.11.06.20226803/T7) Supplementary Table 1. Infection related survey questions. ### Wearable Monitoring Device and Autonomic Nervous System Assessment HRV was measured via the Apple Watch Series 4 or 5, which are commercially available wearable devices. Participants wore the device on the wrist and connected it via Bluetooth to their iPhone. The Watch is equipped with an enhanced photoplethysmogram (PPG) optical heart sensor that combines a green LED light paired with a light sensitive photodiode generating time series peaks that correlate with the magnitude of change in the green light generated from each heartbeat.17 Data are filtered for ectopic beats and artifact. The time difference between heartbeats is classified as the Interbeat Interval (IBI) from which HRV is calculated. The Apple Watch and the Apple Health app automatically calculate HRV using the standard deviation of the IBI of normal sinus beats (SDNN), measured in milliseconds (ms). This time domain index reflects both sympathetic and parasympathetic nervous system activity and is calculated by the Apple Watch during ultra-short-term recording periods of approximately 60 seconds.11 The Apple Watch generates several HRV measurements throughout a 24-hour period. HRV metrics are stored in a locally encrypted database accessible through the iPhone Health app which is retrieved through our custom Warrior Watch App. Data is transferred from the iPhone and Apple Watch upon completion of the e-consent and any survey in the App. Wearable data is stored locally allowing retrieval during the days when surveys are not completed by participants. ### Statistical Analysis #### Heart Rate Variability Modelling The HRV data collected through the Apple Watch was characterized by a circadian pattern, a sparse sampling over a 24-hour period, and a non-uniform timing across days and participants. These characteristics bias easily derived features including mean, maximum and minimum creating the need to derive methods that model the circadian rhythm of HRV. A COSINOR model was used to model daily circadian rhythm over a 24 hour period with the non-linear function Y(t) = M+*Acos*(2*π*t/*τ* + *ϕ*) + ei(t) [equation 1], where τ is the period (*τ* =24h), M is the Midline Statistic of Rhythm (MESOR), a rhythm-adjusted mean, A is the amplitude, a measure of half the extent of variation within a day and F is the Acrophase, a measure of the time of overall high values recurring in each day (**Supplementary Figure 1**). This non-linear model with 3 parameters has the advantage of being easily transformed into a linear model by recoding time (t) into two new variables x and z as *x* = sin(2*π*t/*τ*), *z* = sin(2*π*t/*τ*). HRV can then be written as Y(t)=M+*β*xt + *γ*zt + ei(t) [equation 2], where the linear coefficients *β, γ* of the linear model in equation 2 are related to the non-linear parameters of the non-linear model in equation 1 by *β* = *Acos*(*ϕ*) *γ* = −*Asin*(*ϕ*). One can estimate the linear parameters *β, γ* and then obtain the A and *ϕ* as: ![Formula][1] We took advantage of the longitudinal structure of the data to identify a participant specific daily pattern and then measured departures from this pattern as a function of COVID-19 diagnosis or other relevant covariates. In order to do so we used a mixed-effect COSINOR model, where the HRV measure of participant i at time t can be written as HRVit = (M+*β*.xit + *γ*.zit) + *W**it*. *θ**i*+ ei(t), ei(t)∼N(0,s), and where M, *β* and *γ* are the population parameters (fixed-effects) and *θ*i is a vector of random effects and assumed to follow a multivariate normal distribution *θ* i∼N(0,Σ). In this context the introduction of random effects intrinsically model the correlation due to the longitudinal sampling. To measure the impact of any covariate C on the participants’ daily curve, we can introduce such covariates as fixed-effects as its interactions with x and z: HRVit = M+*a*oCi+(*β* + *a*2Ci).xit + (*γ* + *a* 3Ci.)zit + *W**it*. *θ**i* + ei(t) [equation 3]. Model parameters and the standard errors of equation 3 can be estimated via maximum likelihood or reweighted least squares (REWL) and hypothesis testing can be carried out for any comparison that can be written as a linear function of *a*′*s, β and γ* parameters. ![Supplementary Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/11/07/2020.11.06.20226803/F4.medium.gif) [Supplementary Figure 1.](http://medrxiv.org/content/early/2020/11/07/2020.11.06.20226803/F4) Supplementary Figure 1. Locally estimated scatterplot smoothing (loess) curve showing a daily circadian pattern on HRV measures. Such pattern can be represented by the COSINOR model using 3 parameters: the rhythm-adjusted mean (MESOR), half the extent of variation within a day (Amplitude) and the time of overall high values recurring in each day (acrophase). Red and green dots represent hypothetical sampling times though the day from two subjects that have the same daily curve, showing that features like maximum, range, or CV will be easily biased by the sampling time. However, to test if the COSINOR curve, defined by the non-linear parameters M, A and *ϕ* in equation 1 differs between the populations defined by the covariate C, we proposed the following bootstrapping procedure where for each resampling iteration we: (1) Fit a linear mixed-effect model using REWL; (2) Estimated the marginal means obtaining the linear parameters for each group defined by covariate C; (3) Used the inverse relationship to estimate marginal means M, A and *ϕ* for each group defined by C; and (4) Defined the bootstrapping statistics as the pairwise differences of M, A and *ϕ* between groups defined by C. For such iterations, the confidence intervals for the non-linear parameter was defined using standard bootstrap techniques, as well deriving the p-values for the differences of each non-linear parameter between groups defined by Ci. Age and sex were included as a covariate in HRV analyses and admitted invariant and time-variant covariates. #### Association and Prediction of COVID-19 Diagnosis and Symptoms The relationship between a COVID-19 diagnosis and change in HRV curves were evaluated. To test this association, we defined the time variant covariate Cit for participant i at time t as: ![Formula][2] HRV metrics for the 14 days following the time of first positive SARS-CoV-2 nasal PCR test were used to define the positive SARS-CoV-2 infection window. To evaluate the predictive ability of changes in HRV prior to a COVID-19 diagnosis and to explore its changes during the infection period, the time variant covariate was used to characterize the following 4 groups: healthy uninfected individuals [t