ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article
Revised

Normalizing sleep quality disturbed by psychiatric polypharmacy: a single patient open trial (SPOT)

[version 2; peer review: 2 approved]
Previously titled: Normalizing sleep quality disturbed by psychiatric polypharmacy and sleep apnea: a comprehensive patient-centered N-of-1 study
PUBLISHED 05 Sep 2016
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

There is a growing interest in personalized and preventive medicine initiatives that leverage serious patient engagement, such as those initiated and pursued among participants in the quantified-self movement. However, many of the self-assessments that result are not rooted in good scientific practices, such as exploiting controls, dose escalation strategies, multiple endpoint monitoring, etc. Areas where individual monitoring and health assessments have great potential involve sleep and behavior, as there are a number of very problematic sleep and behavior-related conditions that are hard to treat without personalization. For example, winter depression or seasonal affective disorder (SAD) is a serious, recurrent, atypical depressive disorder impacting millions each year. In order to prevent yearly recurrence antidepressant drugs are used to prophylactically treat SAD. In turn, these antidepressant drugs can affect sleep patterns, further exacerbating the condition. Because of this, possibly unique combinatorial or ‘polypharmaceutical’ interventions involving sleep aids may be prescribed. However, little research into the effects of such polypharmacy on the long-term sleep quality of treated individuals has been pursued. Employing wireless monitoring in a patient-centered study we sought to gain insight into the influence of polypharmacy on sleep patterns and the optimal course of therapy for an individual being treated for SAD with duloxetine (Cymbalta) and temazepam. We analyzed continuous-time sleep data while dosages and combinations of these agents were varied. We found that the administration of Cymbalta led to an exacerbation of the subject’s symptoms in a statistically significant way. We argue that such analyses may be necessary to effectively treat individuals with similar overall clinical manifestations and diagnosis, despite their having a unique set of symptoms, genetic profiles and exposure histories. We also consider the limitations of our study and areas for further research.

Keywords

sleep, depression, polypharmacy, wireless monitoring, sleep apnea, personalized medicine, N-of-1 trials

Revised Amendments from Version 1

We have re-titled and extensively revised our manuscript to provide a concise focus on the medication trial aspect of the N-of-1 study, the type of which is properly defined as a single patient open trial (SPOT). 
 
Major sections were shortened, reorganized and amended as follows:

  • Superfluous patient history was removed from Introduction and Methods
     
  • Table 1 (drug trials description) was moved from Methods to Results
     
  • Text was added to Methods and Limitations to address and justify study design choices regarding medication blinding, washout periods and cross-overs
     
  • Suggested future directions were removed from Discussion
In order to simplify the manuscript, the auxiliary follow up study of mild sleep apnea was moved to the Supplementary Material. Rather than dispose of this work entirely, we chose to continue to make it available for two reasons:
  • It reflects the reproducibility of the subject’s basic no-drug sleep architecture 6 months after the medication trial.
     
  • It reveals the impact of a clinical intervention for mild sleep apnea for those who are interested

See the authors' detailed response to the review by Wilson D. Pace

Introduction

Winter depression or seasonal affective disorder (SAD) is an atypical depressive disorder that in most cases has onset in fall or winter with remission in spring or summer. It is estimated that approximately 5–10 percent of people in the U.S. (i.e., 10–20 million people) experience varying degrees of SAD in a given year1. While full syndromal SAD (frequently dependent on additional external negative stressors) is not reached every year, subsyndromal symptoms can be seen2. These symptoms are multiple, and include varying degrees of hypersomnia, carbohydrate-craving and jet-lagged physical and mental states (what is known as “brain fog”) resulting in fatigue and irritability. The annual shortening of the photoperiod is believed to be the main factor in SAD onset; however, responses to cold temperatures and epigenetic changes have been documented in seasonal mammals and exhibit evolutionary conservation down to lower forms of life36, suggesting that many very basic physiologic mechanisms could contribute to SAD. Ultimately, SAD is a complex disease with both chronobiological and neurobiological underpinnings711, which may include an etiology that for some could even begin in utero1216.

Treating SAD is far from trivial and will require tailoring the treatment to an individual and his or her circumstances, for a whole host of reasons, not the least of which concern both individual and societal expectations regarding work habits, lifestyle, communal conventions surrounding day vs. nighttime activities, and the use of pharmacotherapies to treat conditions affecting behavior. In addition, SAD, and depressive syndromes in general, are known to be accompanied by many co-morbidities and sequelae, including anxiety, detrimental body habitus, anhedonia, and, more importantly, sleep disturbances which may exacerbate any underlying depression as well as the additional associated conditions2. Tailored treatments for each and every condition possessed by an individual patient who also has SAD could adversely affect that patient’s sleep, thereby creating negative feedback for the SAD-related and other symptoms. Treatment of SAD includes a general recommendation for morning bright light therapy and/or antidepressant treatment which can be somewhat effective in managing symptoms, while melatonin, exercise and negative ion therapy are also suggested. However, a recent critical review of light therapy literature showed that most bright light therapy studies have methodological issues and evidence is not unequivocal17. Further, cognitive response to bright light therapy can vary based on genetics18. A proper prescription for light therapy requires knowing the dim light melatonin onset (DLMO) of SAD individuals (2/3 are phase-delayed) to determine circadian phase19. The same is true for using supplemental melatonin to advance sleep phase, as improper timing and dosing can exacerbate symptoms19. Because of the seasonal “on-off” nature of the disorder and difficulty in long-term compliance with bright light therapy (due to eyestrain and lack of individualized prescription), year-round prophylactic treatment with antidepressants may be prescribed.

Treatment for SAD and its sequelae are also compounded for peri- and post-menopausal females – a fact which may be under-appreciated in the primary care setting. The progression to menopause in normal women can result in circadian rhythm, vasomotor, and sleep disturbances and an increased risk for depression, possibly further exacerbating symptoms2022. Therefore, a clinician’s choice to potentially increase the dosage of, e.g., a previously effective SSRI antidepressant can in turn exacerbate side effects, such as sleep disturbances. Importantly, sleep apnea is one of the most under-diagnosed conditions in post-menopausal women and is a leading cause of cardiovascular morbidity and mortality2329. Prescribing sleep medications to aid in depression-related symptoms in peri- or post-menopausal women that may be susceptible, or have, sleep apnea is therefore highly problematic.

The fact that depression and sleep disturbances go hand in hand thus creates even more difficult treatment challenges. For example, ironically, it is known that many first-generation antidepressants exert their effects by, among other things, restoring sleep. Unfortunately, many second-generation antidepressants disrupt sleep. It is now accepted that SSRIs and SNRIs typically used to treat SAD can cause sleep disturbances, both in sleep quality (sleep initiation and maintenance) and sleep architecture (rapid eye movement (REM) and non-REM (NREM) sleep)3035. Further, these agents can induce or escalate parasomnias such as periodic leg movements (PLMs) and restless legs syndrome (RLS)36,37. These effects on sleep could further lead clinicians to routinely prescribe sleep medications to counter the stimulating effects of antidepressants, as was recommended for insomnia in patients taking fluoxetine3840. However, sleep medications can have their own negative impacts on sleep quality and architecture, and are not recommended for maintenance use. Thus, the resulting polypharmacy used to treat SAD is usually pursued without regard to the timing or dosage of the drugs or concern for drug-drug interactions. This fact, combined with unique patient characteristics such as age, gender, genetic and exposure profile, and co-morbid conditions, can further impact response to any prescribed drug or drug combination and may change over time.

In order to combat these issues, the management of SAD and related psychiatric disorders should, as noted, be pursued in a more patient-specific or ‘personalized’ manner – something that might not be accomplished at the level of a primary care provider. How such personalization can be achieved generally is an open question given the costs associated with the extra time a clinician might have to spend with a patient to determine an optimal course of therapy, but does suggest a greater number of empirical studies investigating the effects of polypharmacy and the utility of different treatment strategies are needed. In addition, patient-acceptance of the challenges surrounding treatment may motivate self-assessments of the type being pursued by members of the quantified self movement but perhaps in more objective ‘N-of-1’ clinical trial like settings41,42. We describe a study investigating the influence of polypharmacy involving a 58-year-old post-menopausal female who was diagnosed with SAD in 2001. The N-of-1 trial design utilized is known as a “single patient open trial” or SPOT41. The SPOT offers an alternative to the typical N-of-1 trial components. A SPOT requires no randomization, no placebo and no blinding and allows limited cross-overs of one or more. The ultimate goals of the study were two-fold: to determine if objective claims about the influence of her treatments on her psychological well-being could be made in a self-assessment-oriented but designed outcome measures study, and whether her medication use correlated with exacerbation of her various symptoms and conditions.

Ultimately, the study leveraged wireless monitoring devices and regression modeling to assess patient sleep quality (e.g., the Zeo Sleep Monitor43,44), and designed a drug removal and dose escalation study to determine drug effects. In the course of the study, a number of important insights were obtained. The study identified a number of statistically significant correlations between medication use and symptomology that led to a number of potential recommendations for future treatments. Although it is important to acknowledge the shortcomings of the study, we feel that such patient-engaged and initiated yet protocol-oriented and designed N-of-1 studies may be the best way to individualize treatments for individuals with multiple mood and sleep-related conditions for which polypharmaceutical interventions are common.

Methods

Participant

We studied a post-menopausal 58-year-old female (the ‘subject’, author VLM) treated for SAD since 2001. The subject was interested in self-monitoring and an N-of-1 study for her sleep disturbances given her lengthy dissatisfaction with available treatment options, lack of insights into her multiple conditions, and a very elaborate and complex treatment history. The subject had a long history of usage of benzodiazepine as a sleep medication while taking antidepressants. The subject loosely qualifies as evening prone or delayed sleep phase disorder according to Basic Language Morningness Scale (BALM) questionnaire, which uses a 6-item scale45. In summer 2012, she reported that under prolonged indoor low-light conditions she was susceptible to feeling fatigued, exhibiting seasonal symptomology even in summer months in San Diego. In fall 2012, the subject was taking 60 mg Cymbalta, 30 mg temazepam for sleep, and 100 mg sumatriptan as needed for morning headaches. An N-of-1 (SPOT design) study was pursued to explore how her medications affected her sleep in the context of her diagnosed winter depression (SAD), evening chronotype, delayed sleep phase, restless legs/PLMs and morning headaches.

Ethics

The present study was self-administered by one of the authors (VLM). Therefore, ethical approval from an Institutional Review Board was not sought because the Helsinki Declaration does not apply in this case.

Measures and wireless devices

Sleep and activity monitoring. To assess sleep patterns a Zeo Sleep Monitor (http://www.myzeo.com, model number ZEO 301) was used, which was worn nightly after entering bed per manufacturer instructions. The Zeo wirelessly tracks sleep stages at 5-minute intervals and has been validated against laboratory polysomnography43. The number of awakenings (after sleep onset), percent time in light, deep, REM and wake were recorded and assessed with an accompanying iPAD application (Zeo Sleep Manager v1.9.0). Until the manufacturer’s bankruptcy, the Zeo online application provided nightly tracking of sleep stages and tools for evaluating trends. In addition, educational materials reminding the user of good sleep hygiene practices and journaling and counseling options were also offered. The data obtained with the Zeo monitor was captured on an iPad and Zeo graphic image data obtained with the device is available from the authors. In addition to the Zeo monitor, an Actiwatch Spectrum (manufactured by Philips Respironics) was used to collect data at 15-second intervals and worn daily to track sleep and light exposure. It was synchronized to the Zeo monitor on the nights it was worn. Because Actiwatch relies on movement to score wake versus sleep, the Actiwatch tends to overestimate time in sleep and underestimate time resting in a quiet awake state (Actiware software version 04.00). Periodic leg movements were measured using the PAM-RL (also manufactured by Philips Respironics) right and left ankle sensors and scored using default settings in software (PAM-RL version 7.6.2). Finally, the Fitbit Ultra actigraphic monitor (http://www.fitbit.com) was worn daily to track walking or “step” activity. The subject wore the Fitbit on her waist from the start of her day through the evening. The Fitbit can be used to monitor sleep activity, but may overestimate sleep time since it keys off of movement (Fitbit app v1.8.2).

Procedures

Pharmacotherapy manipulation: effect on sleep. A schedule was developed for evaluating the effects of Cymbalta, temazepam and melatonin on the subject. Fourteen trials were conducted from 12-30-2012 through 07-05-2013. Description of the 14 trials and the number of nights with complete data are presented in the Results. Essentially, Cymbalta and temazepam were provided to the subject in pre-specified time periods with pre-specified doses initiated on weekends. Melatonin (Nature Made, 3 mg chocolate melts) was used to attempt to phase-shift the subject as needed to keep a work schedule, but several periods involving different combinations were pursued to explore the influence of melatonin on phase. Consistent with a SPOT design by definition and rationale, the study was pursued without randomization in a real-time, real-life setting, similar to a clinical practice drug de-escalation/withdrawal, and no medication blinding was utilized. In addition, because of the strong effects of the medications on our subject any placebo would have been detected. Similarly, a “no treatment washout period” between treatments was not employed or even feasible. There are several reasons for this, first, not wanting to destroy the continuity of the biological effects; but second, and more importantly, complete Cymbalta withdrawal causes undesirable side-effect symptoms such as “brain-zaps” for several months, the duration of which cannot be predicted. Hence in this case, washouts designed into this type of study would extend the timetable while causing further harms. We accept that this would add carryover and rebound effects at treatment boundaries. As an underlying goal of the study was to eliminate the benzodiazepine temazepam and to determine if any combination of Cymbalta and/or melatonin could normalize our subject’s sleep, we took an adaptive approach for which treatment cross-overs were only included in the latter portion of the study. It should also be noted that in designing a study like the one described there are a number of potential confounding variables that inevitably arise especially in any naturalistic, free-living setting assessing sleep quality: a) sleep consolidation could occur as sleep deprivation leads to sleep pressure as week progresses; b) sleeping in and changing sleep patterns on weekends could affect weekday trends; and c) percent time in wake after sleep onset can be increased by PLMs, sleep apnea or other sleep maintenance problems, which could be compounded by medication use.

Sleep analysis. Each night and morning, the subject manually entered start and stop times into the Zeo sleep monitor iPAD app. The time to REM sleep was manually calculated based on Zeo graphic histogram output showing first REM sleep bar. Percent wake, light, deep and REM sleep and number of awakenings were supplied by the Zeo device. We did not use the Zeo sleep latency parameter “Time to Z” due to the confounding presence of PLMs, which our subject has shown to exhibit upon sleep initiation (clinically validated via videotape). The subject also wore the Actiwatch Spectrum around the clock from April 2013 until August 2013 as well as the PAM-RL ankle sensors nightly from April 2013 to July 2013. Some missing sleep quality data occurred due to days for which the subject was traveling.

General statistical analysis

All analyses were performed using R version 3.1.3 (http://www.R-project.org). For the sleep analysis, the data used contained information for 188 consecutive nights from December 30, 2012 to July 5, 2013 with 21 nights having missing data attributable to lost records and was therefore treated as missing at random (MAR). The response variables focusing on sleep quality included the number of wakes, time to first REM sleep, percent time in REM sleep, percent time in deep sleep, percent time in light sleep, and percent time in wake. To accommodate the presence of serial correlation in the nightly data, linear models considering an autoregressive moving average (ARMA) serial correlation structure among the data were fit. Different assumptions about the degree of serial correlation were made and tested. Interestingly, little evidence for a strong serial correlation was found, and therefore simple univariate linear regressions were used for all response variables via the lm function in R, retaining predictor variables significant at p < 0.05. Analyses involving model residuals were pursued to assess goodness-of-fit and satisfaction of linear model criteria. These included a Durbin-Watson test (to detect serial correlation between residual values), Shapiro-Wilk normality check, Portmanteau test and ARCH test. In cases where residuals in final models did not satisfy normality, a Box-Cox procedure was performed on the model. The resulting optimal exponential transformation was applied to the response variable and the model refit. To determine best fit among similar models, linear regression model fit measures (Akaike information criteria (AIC), Bayesian information criteria (BIC) and log likelihood) were evaluated. Only the best final models meeting all linear model criteria including no serial correlation or autocorrelation are presented in the results. The univariate regression models for each dependent variable were pursued in very similar ways, as outlined in the following example. Let perstage tdenote series analysis response variables, where non-transformed variables are percent wake (perwake), percent light (perlight), percent deep (perdeep) and percent REM (perrem).

Mathematics

To be more specific, an example model for perstaget was created to follow the simple scheme below, with other variables leveraging similar models:

perstaget = μ0 + βcym30cym30 + βcym60cym60 + βmel3mel3 + βcym30mel3cym30mel3 + βcym30mel6cym30mel6 + βcym60mel3cym60mel3 + βcym60mel6cym60mel6 + βcym60tem15cym60tem15 + βcym60tem30cym60tem30 + t

where μ0 is a y-intercept term, the β terms are regression coefficients, t is an error term with 0 mean and variance σ2. The other terms in the model correspond to the drugs being evaluated and are denoted as follows: Cymbalta 30 mg (cym30); Cymbalta 60 mg (cym60); Melatonin 3 mg (mel3); Cymbalta 30 mg and Melatonin 3 mg (cym30mel3); Cymbalta 30 mg and Melatonin 6 mg (cym30mel6); Cymbalta 60 mg and Melatonin 3 mg (cym60mel3); Cymbalta 60 mg and Melatonin 6 mg (cym60mel6); Cymbalta 60 mg and Temazepam 15 mg (cym60tem15); Cymbalta 60 mg and Temazepam 30 mg (cym60tem30). Significant terms (i.e., p < 0.05 based on t-test of the coefficient value and its standard error) in the model were evaluated in an overall model fit as well as in a step-wise manner. Models were also fit to assess the impact of study design (night in time course) and days of the week (using Sunday as comparator per convention) by including these factors as independent variables in the model. The same analyses were performed for time to REM sleep.

Results

GroupTTOREMPERWAKEPERREMPERLIGHTPERDEEPNWAKESNIGHTDAYDAYCODEDAYSUNDAYDAYMONDAYDAYTUESDAYDAYWEDNESDAYDAYTHURSDAYDAYFRIDAYDAYSATURDAYCYMDOSEMELDOSETEMDOSECutgroupDateCYM30CYM60MEL3MEL6TEM15TEM30CYM30MEL3CYM30MEL6CYM60MEL3CYM60MEL6CYM60TEM15CYM60TEM30
A-CYM60TEM304.170.140.210.60.06111ASunday01000000202Cymbalta60mg12/30/12000000000001
A-CYM60TEM304.50.220.210.540.04122Monday10100000202Cymbalta60mg12/31/12000000000001
A-CYM60TEM305.170.160.170.630.04163Tuesday20010000202Cymbalta60mg1/1/13000000000001
A-CYM60TEM305.080.180.160.620.04194Wednesday30001000202Cymbalta60mg1/2/13000000000001
A-CYM60TEM305.580.260.180.490.06155Thursday40000100202Cymbalta60mg1/3/13000000000001
A-CYM60TEM300.920.110.250.610.03176Friday50000010202Cymbalta60mg1/4/13000000000001
A-CYM60TEM30NANANANANANA7Saturday60000001NANANACymbalta60mg1/5/13NANANANANANANANANANANANA
A-CYM60TEM300.170.040.260.660.0498ASunday01000000202Cymbalta60mg1/6/13000000000001
A-CYM60TEM302.920.20.160.610.03139Monday10100000202Cymbalta60mg1/7/13000000000001
A-CYM60TEM304.920.20.170.60.031210Tuesday20010000202Cymbalta60mg1/8/13000000000001
A-CYM60TEM305.420.140.210.60.041411Wednesday30001000202Cymbalta60mg1/9/13000000000001
A-CYM60TEM302.080.060.370.50.08812Thursday40000100202Cymbalta60mg1/10/13000000000001
B-CYM60TEM153.920.160.160.60.071113Friday50000010201Cymbalta60mg1/11/13000000000010
B-CYM60TEM15NANANANANANA14Saturday60000001NANANACymbalta60mg1/12/13NANANANANANANANANANANANA
B-CYM60TEM154.420.240.130.560.071515ASunday01000000201Cymbalta60mg1/13/13000000000010
B-CYM60TEM154.420.180.180.540.11416Monday10100000201Cymbalta60mg1/14/13000000000010
B-CYM60TEM154.580.20.110.620.071317Tuesday20010000201Cymbalta60mg1/15/13000000000010
B-CYM60TEM154.50.210.190.540.061418Wednesday30001000201Cymbalta60mg1/16/13000000000010
B-CYM60TEM153.330.10.220.590.11319Thursday40000100201Cymbalta60mg1/17/13000000000010
B-CYM60TEM153.830.130.220.580.071720Friday50000010201Cymbalta60mg1/18/13000000000010
C-CYM604.50.270.110.520.11821Saturday60000001200Cymbalta60mg1/19/13010000000000
C-CYM605.50.330.130.410.131122ASunday01000000200Cymbalta60mg1/20/13010000000000
C-CYM604.420.130.190.490.191423Monday10100000200Cymbalta60mg1/21/13010000000000
C-CYM603.330.190.180.450.171324Tuesday20010000200Cymbalta60mg1/22/13010000000000
C-CYM606.250.230.110.470.18925Wednesday30001000200Cymbalta60mg1/23/13010000000000
C-CYM603.080.120.240.460.181626Thursday40000100200Cymbalta60mg1/24/13010000000000
C-CYM604.250.180.230.430.161827Friday50000010200Cymbalta60mg1/25/13010000000000
C-CYM603.330.110.280.440.171928Saturday60000001200Cymbalta60mg1/26/13010000000000
C-CYM603.670.240.190.40.171529ASunday01000000200Cymbalta60mg1/27/13010000000000
C-CYM604.830.260.170.460.111430Monday10100000200Cymbalta60mg1/28/13010000000000
C-CYM605.920.160.20.510.141431Tuesday20010000200Cymbalta60mg1/29/13010000000000
C-CYM604.580.190.220.430.161532Wednesday30001000200Cymbalta60mg1/30/13010000000000
C-CYM6050.170.210.450.171733Thursday40000100200Cymbalta60mg1/31/13010000000000
C-CYM603.420.090.30.450.161534Friday50000010200Cymbalta60mg2/1/13010000000000
C-CYM604.330.330.230.320.121735Saturday60000001200Cymbalta60mg2/2/13010000000000
C-CYM6040.420.160.330.091536ASunday01000000200Cymbalta60mg2/3/13010000000000
C-CYM603.580.260.150.390.191137Monday10100000200Cymbalta60mg2/4/13010000000000
C-CYM606.50.30.220.340.151338Tuesday20010000200Cymbalta60mg2/5/13010000000000
C-CYM604.50.120.250.460.162039Wednesday30001000200Cymbalta60mg2/6/13010000000000
C-CYM604.670.350.190.340.121040Thursday40000100200Cymbalta60mg2/7/13010000000000
C-CYM602.750.110.280.470.142141Friday50000010200Cymbalta60mg2/8/13010000000000
C-CYM602.330.060.360.390.19842Saturday60000001200Cymbalta60mg2/9/13010000000000
C-CYM601.750.340.110.410.14943ASunday01000000200Cymbalta60mg2/10/13010000000000
C-CYM604.330.280.190.410.121444Monday10100000200Cymbalta60mg2/11/13010000000000
C-CYM603.50.170.190.490.151445Tuesday20010000200Cymbalta60mg2/12/13010000000000
D-CYM60MEL3-63.830.20.230.420.151546Wednesday30001000210Cymbalta60mg2/13/13000000001000
D-CYM60MEL3-64.250.190.240.440.131547Thursday40000100210Cymbalta60mg2/14/13000000001000
D-CYM60MEL3-64.170.240.280.370.121848Friday50000010210Cymbalta60mg2/15/13000000001000
D-CYM60MEL3-6NANANANANANA49Saturday60000001NANANACymbalta60mg2/16/13NANANANANANANANANANANANA
D-CYM60MEL3-640.350.250.290.11950ASunday01000000210Cymbalta60mg2/17/13000000001000
D-CYM60MEL3-65.750.30.140.440.121651Monday10100000220Cymbalta60mg2/18/13000000000100
D-CYM60MEL3-6NANANANANANA52Tuesday20010000NANANACymbalta60mg2/19/13NANANANANANANANANANANANA
D-CYM60MEL3-65.830.260.190.420.131053Wednesday30001000220Cymbalta60mg2/20/13000000000100
D-CYM60MEL3-64.330.230.250.380.131654Thursday40000100220Cymbalta60mg2/21/13000000000100
D-CYM60MEL3-6NANANANANANA55Friday50000010NANANACymbalta60mg2/22/13NANANANANANANANANANANANA
D-CYM60MEL3-64.330.260.230.390.131856Saturday60000001210Cymbalta60mg2/23/13000000001000
D-CYM60MEL3-64.830.290.170.420.131557ASunday01000000220Cymbalta60mg2/24/13000000000100
D-CYM60MEL3-6NANANANANANA58Monday10100000NANANACymbalta60mg2/25/13NANANANANANANANANANANANA
D-CYM60MEL3-6NANANANANANA59Tuesday20010000NANANACymbalta60mg2/26/13NANANANANANANANANANANANA
D-CYM60MEL3-6NANANANANANA60Wednesday30001000NANANACymbalta60mg2/27/13NANANANANANANANANANANANA
E-CYM30MEL3-61.920.120.250.450.181461Thursday40000100120Cymbalta30mg2/28/13000000010000
E-CYM30MEL3-62.080.260.280.320.151562Friday50000010110Cymbalta30mg3/1/13000000100000
E-CYM30MEL3-62.420.10.380.420.11463Saturday60000001110Cymbalta30mg3/2/13000000100000
E-CYM30MEL3-62.750.150.270.420.161064ASunday01000000110Cymbalta30mg3/3/13000000100000
E-CYM30MEL3-63.830.250.220.440.11265Monday10100000110Cymbalta30mg3/4/13000000100000
E-CYM30MEL3-630.060.330.490.121066Tuesday20010000120Cymbalta30mg3/5/13000000010000
E-CYM30MEL3-630.170.330.380.122267Wednesday30001000120Cymbalta30mg3/6/13000000010000
E-CYM30MEL3-64.080.20.20.480.12568Thursday40000100120Cymbalta30mg3/7/13000000010000
E-CYM30MEL3-63.250.180.280.40.141869Friday50000010110Cymbalta30mg3/8/13000000100000
E-CYM30MEL3-640.280.280.290.151470Saturday60000001110Cymbalta30mg3/9/13000000100000
E-CYM30MEL3-64.080.280.20.360.161771ASunday01000000110Cymbalta30mg3/10/13000000100000
E-CYM30MEL3-63.670.390.190.350.071072Monday10100000120Cymbalta30mg3/11/13000000010000
E-CYM30MEL3-630.20.30.380.121473Tuesday20010000120Cymbalta30mg3/12/13000000010000
E-CYM30MEL3-62.920.160.20.460.18974Wednesday30001000110Cymbalta30mg3/13/13000000100000
E-CYM30MEL3-62.920.090.250.490.171275Thursday40000100110Cymbalta30mg3/14/13000000100000
E-CYM30MEL3-630.270.30.30.131376Friday50000010110Cymbalta30mg3/15/13000000100000
E-CYM30MEL3-6NANANANANANA77Saturday60000001NANANACymbalta30mg3/16/13NANANANANANANANANANANANA
E-CYM30MEL3-64.250.320.190.370.121478ASunday01000000120Cymbalta30mg3/17/13000000010000
E-CYM30MEL3-63.250.240.210.410.151679Monday10100000120Cymbalta30mg3/18/13000000010000
E-CYM30MEL3-62.250.230.310.320.131680Tuesday20010000110Cymbalta30mg3/19/13000000100000
E-CYM30MEL3-63.170.140.240.440.181381Wednesday30001000110Cymbalta30mg3/20/13000000100000
E-CYM30MEL3-62.420.130.340.350.181782Thursday40000100110Cymbalta30mg3/21/13000000100000
E-CYM30MEL3-6NANANANANANA83Friday50000010NANANACymbalta30mg3/22/13NANANANANANANANANANANANA
E-CYM30MEL3-62.750.110.290.420.181984Saturday60000001110Cymbalta30mg3/23/13000000100000
E-CYM30MEL3-64.670.420.180.30.091385ASunday01000000120Cymbalta30mg3/24/13000000010000
E-CYM30MEL3-640.30.250.320.131886Monday10100000120Cymbalta30mg3/25/13000000010000
E-CYM30MEL3-64.250.370.210.30.12987Tuesday20010000110Cymbalta30mg3/26/13000000100000
F-CYM30MEL3T12.750.110.190.460.241488Wednesday30001000110Cymbalta30mg3/27/13000000100000
F-CYM30MEL3T13.750.180.330.360.142189Thursday40000100110Cymbalta30mg3/28/13000000100000
F-CYM30MEL3T14.920.320.220.340.121990Friday50000010110Cymbalta30mg3/29/13000000100000
F-CYM30MEL3T12.670.120.290.420.171791Saturday60000001110Cymbalta30mg3/30/13000000100000
F-CYM30MEL3T14.250.230.220.390.171392ASunday01000000110Cymbalta30mg3/31/13000000100000
F-CYM30MEL3T13.250.210.260.370.161393Monday10100000110Cymbalta30mg4/1/13000000100000
F-CYM30MEL3T12.670.150.210.450.181294Tuesday20010000110Cymbalta30mg4/2/13000000100000
F-CYM30MEL3T12.830.130.280.420.171895Wednesday30001000110Cymbalta30mg4/3/13000000100000
F-CYM30MEL3T12.670.150.240.450.151696Thursday40000100110Cymbalta30mg4/4/13000000100000
F-CYM30MEL3T14.170.280.270.30.151497Friday50000010110Cymbalta30mg4/5/13000000100000
G-CYM30T12.420.030.350.430.191098Saturday60000001100Cymbalta30mg4/6/13100000000000
G-CYM30T13.830.180.220.420.181399ASunday01000000100Cymbalta30mg4/7/13100000000000
G-CYM30T14.080.130.280.420.1719100Monday10100000100Cymbalta30mg4/8/13100000000000
G-CYM30T15.420.450.090.330.1310101Tuesday20010000100Cymbalta30mg4/9/13100000000000
G-CYM30T15.250.270.240.360.1317102Wednesday30001000100Cymbalta30mg4/10/13100000000000
G-CYM30T13.580.180.230.430.1623103Thursday40000100100Cymbalta30mg4/11/13100000000000
G-CYM30T12.830.140.280.330.259104Friday50000010100Cymbalta30mg4/12/13100000000000
G-CYM30T13.170.080.320.440.1515105Saturday60000001100Cymbalta30mg4/13/13100000000000
G-CYM30T1NANANANANANA106ASunday01000000NANANACymbalta30mg4/14/13NANANANANANANANANANANANA
G-CYM30T14.080.210.290.360.1418107Monday10100000100Cymbalta30mg4/15/13100000000000
G-CYM30T14.580.210.260.360.1713108Tuesday20010000100Cymbalta30mg4/16/13100000000000
G-CYM30T14.580.280.230.340.1518109Wednesday30001000100Cymbalta30mg4/17/13100000000000
G-CYM30T14.830.30.160.40.1514110Thursday40000100100Cymbalta30mg4/18/13100000000000
G-CYM30T12.830.160.310.370.1513111Friday50000010100Cymbalta30mg4/19/13100000000000
G-CYM30T13.420.160.320.370.1620112Saturday60000001100Cymbalta30mg4/20/13100000000000
G-CYM30T14.670.240.150.460.1510113ASunday01000000100Cymbalta30mg4/21/13100000000000
G-CYM30T12.750.280.250.320.1519114Monday10100000100Cymbalta30mg4/22/13100000000000
G-CYM30T13.420.160.230.440.1710115Tuesday20010000100Cymbalta30mg4/23/13100000000000
H-CYM30MEL3T23.170.210.290.360.1421116Wednesday30001000110Cymbalta30mg4/24/13000000100000
H-CYM30MEL3T24.750.470.180.240.1113117Thursday40000100110Cymbalta30mg4/25/13000000100000
H-CYM30MEL3T2NANANANANANA118Friday50000010NANANACymbalta30mg4/26/13NANANANANANANANANANANANA
H-CYM30MEL3T26.580.320.210.340.1225119Saturday60000001110Cymbalta30mg4/27/13000000100000
H-CYM30MEL3T23.670.460.20.230.1111120ASunday01000000110Cymbalta30mg4/28/13000000100000
H-CYM30MEL3T25.170.390.170.340.117121Monday10100000110Cymbalta30mg4/29/13000000100000
H-CYM30MEL3T24.580.260.230.370.1510122Tuesday20010000110Cymbalta30mg4/30/13000000100000
H-CYM30MEL3T23.170.220.290.350.1416123Wednesday30001000110Cymbalta30mg5/1/13000000100000
I-CYM30T2NANANANANANA124Thursday40000100NANANACymbalta30mg5/2/13NANANANANANANANANANANANA
I-CYM30T2NANANANANANA125Friday50000010NANANACymbalta30mg5/3/13NANANANANANANANANANANANA
I-CYM30T22.330.160.30.420.1218126Saturday60000001100Cymbalta30mg5/4/13100000000000
I-CYM30T24.50.260.220.420.115127ASunday01000000100Cymbalta30mg5/5/13100000000000
I-CYM30T2NANANANANANA128Monday10100000NANANACymbalta30mg5/6/13NANANANANANANANANANANANA
I-CYM30T2NANANANANANA129Tuesday20010000NANANACymbalta30mg5/7/13NANANANANANANANANANANANA
I-CYM30T23.330.160.280.380.1717130Wednesday30001000100Cymbalta30mg5/8/13100000000000
I-CYM30T23.330.160.330.390.1211131Thursday40000100100Cymbalta30mg5/9/13100000000000
I-CYM30T22.50.070.350.420.1613132Friday50000010100Cymbalta30mg5/10/13100000000000
I-CYM30T230.130.260.440.1616133Saturday60000001100Cymbalta30mg5/11/13100000000000
I-CYM30T23.580.180.280.410.1317134ASunday01000000100Cymbalta30mg5/12/13100000000000
I-CYM30T22.830.130.240.460.1716135Monday10100000100Cymbalta30mg5/13/13100000000000
I-CYM30T23.750.140.310.410.1421136Tuesday20010000100Cymbalta30mg5/14/13100000000000
I-CYM30T22.830.070.340.410.1712137Wednesday30001000100Cymbalta30mg5/15/13100000000000
I-CYM30T24.330.150.320.420.1123138Thursday40000100100Cymbalta30mg5/16/13100000000000
J-CYM30MEL3T34.170.260.20.440.1118139Friday50000010110Cymbalta30mg5/17/13000000100000
J-CYM30MEL3T32.750.110.340.420.1315140Saturday60000001110Cymbalta30mg5/18/13000000100000
J-CYM30MEL3T34.880.450.170.30.0714141ASunday01000000110Cymbalta30mg5/19/13000000100000
J-CYM30MEL3T33.40.190.270.370.1719142Monday10100000110Cymbalta30mg5/20/13000000100000
J-CYM30MEL3T31.820.180.280.370.1715143Tuesday20010000110Cymbalta30mg5/21/13000000100000
J-CYM30MEL3T33.720.20.220.410.1718144Wednesday30001000110Cymbalta30mg5/22/13000000100000
J-CYM30MEL3T36.030.260.230.370.1412145Thursday40000100110Cymbalta30mg5/23/13000000100000
J-CYM30MEL3T33.250.090.310.420.1814146Friday50000010110Cymbalta30mg5/24/13000000100000
J-CYM30MEL3T32.530.120.280.410.1810147Saturday60000001110Cymbalta30mg5/25/13000000100000
J-CYM30MEL3T33.330.150.30.380.1711148ASunday01000000110Cymbalta30mg5/26/13000000100000
J-CYM30MEL3T33.720.150.260.410.1813149Monday10100000110Cymbalta30mg5/27/13000000100000
J-CYM30MEL3T32.930.160.270.40.1713150Tuesday20010000110Cymbalta30mg5/29/13000000100000
J-CYM30MEL3T3NANANANANANA151Wednesday30001000NANANACymbalta30mg5/28/13NANANANANANANANANANANANA
J-CYM30MEL3T33.770.150.180.470.1911152Thursday40000100110Cymbalta30mg5/30/13000000100000
J-CYM30MEL3T32.50.160.280.420.1317153Friday50000010110Cymbalta30mg5/31/13000000100000
K-MEL3T12.650.040.370.410.188154Saturday60000001010Cymbalta0mg6/1/13001000000000
K-MEL3T11.950.050.410.350.26155ASunday01000000010Cymbalta0mg6/2/13001000000000
K-MEL3T11.180.270.380.210.147156Monday10100000010Cymbalta0mg6/3/13001000000000
L-NODRUGT11.30.190.380.210.225157Tuesday20010000000Cymbalta0mg6/4/13000000000000
L-NODRUGT11.30.080.440.240.247158Wednesday30001000000Cymbalta0mg6/5/13000000000000
L-NODRUGT11.070.10.430.260.228159Thursday40000100000Cymbalta0mg6/7/13000000000000
L-NODRUGT1NANANANANANA160Friday50000010NANANACymbalta0mg6/6/13NANANANANANANANANANANANA
L-NODRUGT10.480.060.440.290.217161Saturday60000001000Cymbalta0mg6/8/13000000000000
L-NODRUGT11.80.080.270.40.256162ASunday01000000000Cymbalta0mg6/9/13000000000000
L-NODRUGT10.880.280.330.210.187163Monday10100000000Cymbalta0mg6/10/13000000000000
L-NODRUGT11.10.030.310.370.295164Tuesday20010000000Cymbalta0mg6/11/13000000000000
M-MEL3T22.220.220.270.370.145165Wednesday30001000010Cymbalta0mg6/12/13001000000000
M-MEL3T20.830.040.410.320.236166Thursday40000100010Cymbalta0mg6/13/13001000000000
M-MEL3T21.180.250.30.280.179167Friday50000010010Cymbalta0mg6/14/13001000000000
M-MEL3T20.70.040.370.340.256168Saturday60000001010Cymbalta0mg6/15/13001000000000
N-NODRUGT21.720.050.320.340.296169ASunday01000000000Cymbalta0mg6/16/13000000000000
N-NODRUGT21.60.040.370.360.235170Monday10100000000Cymbalta0mg6/17/13000000000000
N-NODRUGT21.050.040.360.360.249171Tuesday20010000000Cymbalta0mg6/18/13000000000000
N-NODRUGT21.250.050.360.370.218172Wednesday30001000000Cymbalta0mg6/19/13000000000000
N-NODRUGT210.110.350.320.218173Thursday40000100000Cymbalta0mg6/20/13000000000000
N-NODRUGT20.930.210.360.250.184174Friday50000010000Cymbalta0mg6/21/13000000000000
N-NODRUGT20.970.020.350.340.286175Saturday60000001000Cymbalta0mg6/22/13000000000000
N-NODRUGT21.480.060.370.380.195176ASunday01000000000Cymbalta0mg6/23/13000000000000
N-NODRUGT22.080.040.30.440.227177Monday10100000000Cymbalta0mg6/24/13000000000000
N-NODRUGT20.650.020.40.340.242178Tuesday20010000000Cymbalta0mg6/25/13000000000000
N-NODRUGT20.870.070.410.360.1613179Wednesday30001000000Cymbalta0mg6/26/13000000000000
N-NODRUGT20.780.050.370.380.196180Thursday40000100000Cymbalta0mg6/27/13000000000000
N-NODRUGT21.070.080.390.360.1710181Friday50000010000Cymbalta0mg6/28/13000000000000
N-NODRUGT2NANANANANANA182Saturday60000001NANANACymbalta0mg6/29/13NANANANANANANANANANANANA
N-NODRUGT2NANANANANANA183ASunday01000000NANANACymbalta0mg6/30/13NANANANANANANANANANANANA
N-NODRUGT2NANANANANANA184Monday10100000NANANACymbalta0mg7/1/13NANANANANANANANANANANANA
N-NODRUGT21.230.040.270.420.279185Tuesday20010000000Cymbalta0mg7/2/13000000000000
N-NODRUGT21.920.090.360.370.1814186Wednesday30001000000Cymbalta0mg7/3/13000000000000
N-NODRUGT20.670.050.380.380.28187Thursday40000100000Cymbalta0mg7/4/13000000000000
Dataset 1.Drug dosage and sleep response data.
Group = Drug trials, TTOREM = Time to REM sleep in hours, PERWAKE = Percent time in Wake, PERREM = Percent time in REM sleep, PERLIGHT = Percent time in Light sleep, PERDEEP = Percent time in Deep sleep, NWAKES = Number of wakes per night, DAY = Day of the week, DAYCODE = numerical code for day of the week, DAYSUNDAY, DAYMONDAY, etc. = contrast data codes, CYMDOSE, MELDOSE, TEMDOSE = dosing codes for Cymbalta, Melatonin, Temazepam, respectively, Cut-group = Cymbalta dose groups for plotting, DATE = Date of medication doses and night of sleep data collection, CYM30 = Cymbalta 30 mg, CYM60 = Cymbalta 60 mg, MEL3 = Melatonin 3 mg, MEL6 = Melatonin 6 mg, TEM15 = Temazepam 15 mg, TEM30 = Temazepam 30 mg, CYM30MEL3 = Cymbalta 30 mg and Melatonin 3 mg, CYM30MEL6 = Cymbalta 30 mg and Melatonin 6 mg, CYM60MEL3 = Cymbalta 60 mg and Melatonin 3 mg, CYM60MEL6 = Cymbalta 60 mg and Melatonin 6 mg, CYM60TEM15 = Cymbalta 60 mg and Temazepam 15 mg, CYM60TEM30 = Cymbalta 60 mg and Temazepam 30 mg, NA = missing data.
Sleep DatePLMs per hour Right anklePLMs per hour Left ankle
4/19/1314.228.1
4/20/1320.916.2
4/21/1327.726.3
4/22/1311.715.8
4/23/1319.313.0
5/15/138.020.3
5/17/139.719.1
5/18/1310.716.7
5/19/1310.825.6
5/20/1311.330.2
5/21/1320.031.0
5/22/139.619.7
5/23/133.619.8
5/24/1312.130.4
5/25/1310.721.9
5/26/1323.029.9
5/27/134.015.3
5/28/1311.722.7
5/29/137.814.7
5/30/136.814.6
5/31/138.210.7
6/1/138.29.1
6/2/139.913.6
6/3/134.815.9
6/4/138.229.9
6/5/135.725.3
6/6/13716
6/7/1311.315.1
6/8/135.19.6
6/9/132.516
6/10/131.49.3
6/11/138.621.6
6/12/135.65.1
6/13/136.412.6
6/14/1310.519.2
6/15/13914.6
6/16/138.214.3
6/17/133.910.7
6/18/136.98
6/19/132.110.4
6/20/132.56.3
6/21/136.810.6
6/22/134.910.8
6/23/1311.217.6
6/24/133.716.6
6/25/134.411.4
6/26/138.711.3
6/27/1369.6
6/28/136.515.1
6/29/136.921.1
6/30/134.814
7/1/136.810
7/2/131.16.1
7/3/139.414
7/4/131.38.4
Dataset 2.PAM-RL Periodic Leg Movement Rates.
Sleep Date = Night of PLM collection, PLMs per hour Right ankle = Periodic leg movement counts per hour scored by PAM-RL software from right ankle sensor, PLMs per hour Left ankle = Periodic leg movement counts per hour scored by PAM-RL software from left ankle sensor.

Sleep quality analyses

Sleep data was collected for 188 consecutive nights from December 30, 2012 to July 5, 2013, with 21 nights having missing data (Dataset 1). A description of the 14 trials and the number of nights with complete data are listed in Table 1 (abbreviations: Cymbalta (CYM); temazepam (TEM); melatonin (MEL)). Table 2 gives a descriptive analysis of the sleep parameters used in the study. The mean and standard deviation (SD) for: the number of times per night the subject was awakened (wakes (N)); time to first REM sleep bout in hours (1st REM (h)); and percentage of time in each sleep stage (wake (%), light (%), deep (%), REM (%)) at each drug dose is shown. The number of days per dose and percent of the total nights are also shown (N days (%)). The dataset was not balanced in the sense that we had different numbers of observations while the subject was on different dosages of a drug.

Table 1. Drug trials and number of nights with complete data.

TrialStart dateNightsCYMTEMMEL
A12-30-121160 mg30 mg0 mg
B01-11-13760 mg15 mg0 mg
C01-19-132560 mg0 mg0 mg
D02-13-13960 mg0 mg3–6 mg
E02-28-132530 mg0 mg3–6 mg
F03-27-131030 mg0 mg3 mg
G04-06-131730 mg0 mg0 mg
H04-24-13730 mg0 mg3 mg
I05-02-131130 mg0 mg0 mg
J05-17-131430 mg0 mg3 mg
K06-01-1330 mg0 mg3 mg
L06-04-1370 mg0 mg0 mg
M06-12-1340 mg0 mg3 mg
N06-16-13170 mg0 mg0 mg

Table 2. Descriptive analysis of drug response variables.

Drug
dose
Response
variable
Cymbalta
mean (SD)
Melatonin
mean (SD)
Temazepam
mean (SD)
0Wakes (N)
1st REM (h)
Wake (%)
Light (%)
Deep (%)
REM (%)
7.03 (2.40)
1.26 (0.51)
9.00 (7.77)
33.68 (6.37)
21.48 (4.19)
35.81 (4.97)
12.67 (4.66)
3.20 (1.56)
16.08 (9.36)
42.66 (9.72)
15.20 (6.26)
25.98 (8.55)
13.11 (4.64)
3.23 (1.38)
18.68 (10.39)
38.48 (6.27)
16.11 (4.29)
26.68 (7.43)
N days (%)31 (18.56%)95 (56.89%)149 (89.22%)
1Wakes (N)
1st REM (h)
Wake (%)
Light (%)
Deep (%)
REM (%)
14.81 (3.80)
3.53 (0.93)
20.63 (9.74)
38.81 (5.46)
14.80 (3.15)
25.70 (5.49)
13.79 (4.20)
3.30 (1.15)
20.67 (10.08)
37.41 (6.20)
15.28 (3.44)
26.66 (5.89)
13.86 (1.86)
4.14 (0.46)
17.43 (4.83)
57.57 (3.05)
7.71 (1.60)
17.29 (4.23)
N days (%)84 (50.30%)58 (34.73%)7 (4.19%)
2Wakes (N)
1st REM (h)
Wake (%)
Light (%)
Deep (%)
REM (%)
14.13 (3.17)
4.15 (1.23)
20.50 (8.39)
47.65 (9.32)
11.42 (4.89)
20.44 (5.69)
13.79 (4.14)
3.97 (1.10)
25.00 (9.81)
39.93 (5.58)
12.36 (2.50)
22.71 (5.97)
13.27 (3.35)
4.14 (0.46)
15.55 (6.68)
58.73 (5.39)
4.45 (1.57)
21.36 (6.22)
N days (%)52 (31.14%)14 (8.38%)11 (6.59%)

Figure 1, Figure 2 and Figure S1Figure S4 (see Supplementary Material) graphically depict the impact of Cymbalta, melatonin and temazepam drug use on the subject’s sleep architecture. Figure 1 and Figure 2 show the percent of time per night that the subject was in deep sleep and light sleep, respectively, during 5-minute intervals detected by the Zeo Sleep Monitor throughout the entire study. Similar figures for the number of times the subject was awakened, time to REM sleep, percent time after sleep onset that the subject was awake and percent time in REM sleep during 5-minute intervals detected by the Zeo Sleep Monitor are presented in the Supplementary Material (Figure S1Figure S4, respectively).

522b7a32-6272-4eae-ad07-c4d826f0686b_figure1.gif

Figure 1. Percent deep sleep per night.

The percent time subject was in deep sleep during 5-minute intervals detected by the Zeo Sleep Monitor throughout the entire study. Dosages of Cymbalta (CYM60 = 60 mg, CYM30 = 30 mg), temazepam (TEM30 = 30 mg, TEM15 = 15 mg) and melatonin (MEL3 = 3 mg, MEL6 = 6 mg) were varied according to combinations A–N (T1, T2, T3 are trial replicates), including no drug trials (L, N).

522b7a32-6272-4eae-ad07-c4d826f0686b_figure2.gif

Figure 2. Percent light sleep per night.

The percent time subject was in light sleep during 5-minute intervals detected by the Zeo Sleep Monitor throughout the entire study. Dosages of Cymbalta (CYM60 = 60 mg, CYM30 = 30 mg), temazepam (TEM30 = 30 mg, TEM15 = 15 mg) and melatonin (MEL3 = 3 mg, MEL6 = 6 mg) were varied according to combinations A-N (T1, T2, T3 are trial replicates), including no drug trials (L, N).

A clear relationship can be seen between temazepam intake and reduced deep sleep in favor of light sleep (Figure 1 and Figure 2). However, Cymbalta had the strongest impact on the subject’s sleep architecture as shown in Figure 3, Figure 4, and Figure 5. Cymbalta intake increased the number of awakenings (Figure 3), time to first REM sleep (Figure 4), percent time after sleep onset that subject was awake (wake) (Figure 5A) and in light sleep (Figure 5B) at the expense of deep (Figure 5C) and REM (Figure 5D) sleep. Removal of Cymbalta decreased the number of awakenings, time to first REM sleep, percent time in wake and light sleep and increased percent time in deep and REM sleep (Figure 3Figure 5).

522b7a32-6272-4eae-ad07-c4d826f0686b_figure3.gif

Figure 3. Wakenings per night by Cymbalta dose.

The number of times per night subject was awake during 5-minute intervals detected by the Zeo Sleep Monitor. Doses of Cymbalta were decreased from 60 mg to 0 mg.

522b7a32-6272-4eae-ad07-c4d826f0686b_figure4.gif

Figure 4. Time to REM sleep per night by Cymbalta dose.

The number of hours (h) per night before subject achieved first REM sleep bout during 5-minute intervals detected by the Zeo Sleep Monitor. Doses of Cymbalta were decreased from 60 mg to 0 mg.

522b7a32-6272-4eae-ad07-c4d826f0686b_figure5.gif

Figure 5. Percent time in wake, light, deep and REM sleep per night by Cymbalta dose.

Percent time after sleep onset subject was awake (A); subject was in light sleep (B); subject was in deep sleep (C); or subject was in REM sleep (D) during 5-minute intervals detected by the Zeo Sleep Monitor. Doses of Cymbalta were decreased from 60 mg to 0 mg.

Because of the free-living nature of our study, the subject’s polypharmacy and struggle to counter sleep disturbances, a large variability in the data is seen. In addition, “normal” sleep staging typically follows a pattern wherein the first non-REM sleep (light plus deep sleep) and REM sleep cycle is completed in 70 to 100 minutes, followed by 90 to 120 minute cycles, with deep sleep bouts gradually disappearing and REM sleep bouts lengthening throughout the night46. Near the end of the night, usually only light and REM sleep periods make up the sleep cycles. As a result, we chose to analyze the percentage of time the subject was in each sleep/wake state, rather than total time. For the purposes of comparing Zeo monitored stages to classically defined sleep stages, we assumed the following to represent approximately “normal” sleep stage percentages: wake 5 percent; light 45–55 percent; deep 20–25 percent; REM 25 percent46.

Table 3 summarizes the results of our univariate analyses when the sleep stages, wake, light, deep and REM, were taken as dependent variables. The univariate linear regression models were performed as described (see Methods) and data is presented as mean percent for each sleep stage with treatment effects adjusted relative to the intercept. Analyses of percent wake and light sleep met Durbin-Watson test criteria once two outlier nights each were removed. Final model diagnosis showed that all linear regression assumption requirements were satisfied except for the normality condition for percent wake and percent light sleep. Therefore, the Box-Cox procedure and transformations were performed and the models refit. Final models satisfied all diagnostic tests and the transformed mean estimate values (denoted as ‘bc’) presented in Table 3 were adjusted and back-transformed to give mean percent wake and light sleep.

Table 3. Univariate regression analysis predicting percent wake, light, deep and REM sleep from drug and dose, as well as day of the week, information.

VariableMean
(%)
Wake
estimatebc
(SE)
p-valueMean
(%)
Light
estimatebc
(SE)
p-valueMean
(%)
Deep
estimate
(SE)
p-valueMean
(%)
REM
estimate
(SE)
p-value
μ010.50.39 (0.02)< 2e-1635.40.16 (0.01)< 2e-1622.30.22 (0.01)< 2e-1634.20.34 (0.01)< 2e-16
βcym30 20.80.13 (0.02)3.2e-0740.00.04 (0.01)0.001515.4-0.07 (0.01)7.7e-1424.8-0.09 (0.01)5.7e-11
βcym6025.30.17 (0.02)6.3e-1143.10.07 (0.01)1.8e-0715.0-0.07 (0.01)3.3e-1418.7-0.15 (0.01)< 2e-16
βmel318.7-0.04 (0.01)0.0068
βcym30mel325.00.17 (0.02)3.4e-1338.40.02 (0.01)0.020115.0-0.07 (0.01)< 2e-1623.6-0.11 (0.01)1.1e-15
βcym30mel625.80.18 (0.03)4.0e-0739.70.03 (0.02)0.035812.2-0.10 (0.01)2.0e-1523.8-0.10 (0.02)1.0e-07
βcym60mel331.00.22 (0.04)1.3e-0612.8-0.09 (0.01)2.1e-0922.2-0.12 (0.02)2.9e-06
βcym60mel630.60.22 (0.05)1.4e-0541.50.05 (0.02)0.035612.8-0.10 (0.02)3.3e-0818.1-0.16 (0.03)1.9e-08
βcym60tem1520.90.13 (0.04)0.001057.60.22 (0.02)< 2e-167.7-0.15 (0.01)< 2e-1615.8-0.18 (0.02)1.0e-14
βcym60tem3017.80.10 (0.03)0.003458.90.23 (0.02)< 2e-164.5-0.18 (0.01)< 2e-1620.5-0.14 (0.02)2.5e-12
βNight
βWeekday
βMonday
βTuesday7.1-0.06 (0.02)0.0123
βWednesday7.4-0.05 (0.02)0.0177
βThursday6.5-0.07 (0.02)0.001737.00.03 (0.01)0.0154
βFriday7.3-0.05 (0.02)0.023937.90.04 (0.01)0.0029
βSaturday4.8-0.11 (0.02)1.5e-0539.60.05 (0.01)3.1e-05
bc:exponent = 0.42bc:exponent = 1.79
Adjusted R2: 0.3900Adjusted R2: 0.6604Adjusted R2: 0.6629Adjusted R2: 0.5460
p-value: 1.3e-13p-value: < 2.2e-16p-value: < 2.2e-16p-value: < 2.2e-16

Adjusted mean in percent, mean estimate or transformed mean estimatebc, standard error (SE) and p-value (Pr > |t-value|). R2: R-squared; bc: Box-Cox transformed variable raised to exponent given in final model. Back-transformation to original units was performed (after adjustments relative to intercept) by taking the nth (exponent) root of estimate.

From Table 3 it is clear many of the drugs, doses and drug combinations have a highly significant and negative impact on deep and REM sleep (with the exception of melatonin at 3 mg). The estimate of the y-intercept (μ0) for the model with deep sleep as the dependent variable suggests that approximately 22.3 percent of the time the subject was in deep sleep without any drug effects (p < 2×10-16). The estimated coefficients for the drug and drug dosage independent variables in the model provide the effect on deep sleep of the drugs. The mean percent deep sleep ranged from 4.5 percent (-0.18 (SE: 0.01), p < 2×10-16) while the subject was taking 60 mg Cymbalta and 30 mg temazepam to 18.7 percent (-0.04 (SE: 0.01), p = 0.0068) while the subject was taking 3 mg melatonin. Although temazepam dosing in combination with Cymbalta had the greatest negative impact on deep sleep in favor of light sleep, Cymbalta alone continued to interfere with deep sleep.

Similarly, the impact of an antidepressant such as Cymbalta is expected to show a decrease in REM sleep, mainly through the delay in REM sleep onset (see Table 4). The estimate of the y-intercept (μ0) for the model with REM sleep as the dependent variable suggests that approximately 34.2 percent of the time the subject was in REM sleep without any drug effects (p < 2×10-16), which might be considered high compared to the usual 25 percent. The estimated coefficients for the drug and drug dosage independent variables in the model provide the effect on REM sleep of the drugs. The mean percent REM sleep ranged from 15.8 percent (-0.18 (SE: 0.02), p = 1×10-14) while the subject was taking 60 mg Cymbalta and 15 mg temazepam to 24.8 percent (-0.09 (SE: 0.01), p = 5.7×10-11) while the subject was taking 30 mg Cymbalta. Interestingly, there was an increase in REM sleep on Thursday, Friday and especially significant on Saturday (39.6 percent (0.05 (SE: 0.01), p = 3.1×10-5)).

Table 4. Univariate series analysis for predicting time to REM sleep.

VariableMean (h)Time to REM sleep
estimatebc (SE)
p-value
μ01.271.09 (0.03)< 2e-16
βcym303.730.54 (0.04)< 2e-16
βcym604.220.62 (0.04)< 2e-16
βmel3
βcym30mel33.550.51 (0.04)< 2e-16
βcym30mel63.430.49 (0.06)8.1e-15
βcym60mel34.330.64 (0.08)1.8e-14
βcym60mel65.160.76 (0.08)4.9e-16
βcym60tem154.270.63 (0.07)< 2e-16
βcym60tem304.430.64 (0.06)< 2e-16
βNight
βWeekday
βFriday1.05-0.08 (0.04)0.0458
βSaturday1.01-0.09 (0.04)0.0174
bc: exponent = 0.375
Adjusted R2: 0.6823
p-value: < 2.2e-16

Adjusted mean in hours (h), transformed mean estimatebc, standard error (SE) and p-value (Pr > |t-value|). R2: R-squared; bc: Box-Cox transformed variable raised to exponent given in final model. Back-transformation to original units was performed (after adjustments relative to intercept) by taking the nth (exponent) root of estimate.

Most drug combinations, except melatonin, significantly increased time in wake and light sleep. Of note, the Zeo monitor can detect micro-arousals as well as conscious wakes. Thus, some scores of the wakes at night may actually be classified as light sleep. However, from Table 2 and Table 3 the drug combinations increase both of these at the expense of deep and REM. The effect of increasing light sleep at the expense of deep sleep is most notably seen with temazepam use. The estimate of the y-intercept (μ0) for the model with light sleep as the dependent variable suggests that approximately 35.4 percent of the time the subject was in light sleep without any drug effects (p < 2×10-16). The estimated coefficients for the drug and drug dosage independent variables in the model provide the effect on light sleep of the drugs. The mean percent light sleep ranged from 38.4 percent (0.02 (SE: 0.01), p = 0.0201) while the subject was taking 30 mg Cymbalta and 3 mg melatonin to 58.9 percent (0.23 (SE: 0.02), p < 2×10-16) while the subject was taking 60 mg Cymbalta and 30 mg temazepam.

The major impact on wake after sleep onset occurred after the removal of temazepam and during Cymbalta use, indicating a possible sleep maintenance issue. The estimate of the y-intercept (μ0) for the model with wake as the dependent variable suggests that approximately 10.5 percent of the time the subject was in wake without any drug effects (p < 2×10-16). The estimated coefficients for the drug and drug dosage independent variables in the model provide the effect on wake of the drugs. The mean percent wake ranged from 17.8 percent (0.10 (SE: 0.03), p = 0.0034) while the subject was taking 60 mg Cymbalta and 30 mg temazepam to 31.0 percent (0.22 (SE: 0.04), p = 1.3×10-6) while the subject was taking 60 mg Cymbalta and 3 mg melatonin. Interestingly, there is evidence for decreased time classified as wake as the week progresses that might be attributed to a number of things such as increasing sleep pressure during the week, relaxed frame of mind and sleeping in on the weekend. In fact, the decrease in wake to 4.8 percent on Saturday seems to approximately parallel the increase in REM sleep on Saturday (approximately 5 percent) with similar p-values. There was no impact of the night of the study on any of the models.

Table 4 shows the univariate analyses of time to REM sleep in hours as a dependent variable. The univariate linear regression model exhibited no serial correlation based on the Durbin-Watson test once two Zeo technical outlier nights were removed (known REML error43). As above, models were also tested for the impact of study design (night in time course) and day of the week. An assessment of the normality and serial correlation among the residuals obtained from the model was performed by Portmanteau test, Durbin-Watson statistic, a standard normality check and ARCH test which showed that all linear regression assumption requirements were satisfied except normality. Therefore, the Box-Cox procedure and transformation was performed, and model refit as above. The mean estimates presented in Table 4 were adjusted and back-transformed to give the original unit of hours.

All drug combinations except for melatonin at the 3 mg dose caused large and highly significant increases in time to first REM sleep. Under normal circumstances the first REM bout is expected to occur before completing the first 70–100 minute full cycle of sleep (light + deep + REM), that is, in less than 2 hours. The estimate of the y-intercept (μ0) for the model with time to first REM sleep as the dependent variable suggests that time to first REM sleep for the subject was 1.27 hours (76.2 minutes) without any drug effects (p < 2×10-16), which is in the correct range for the first full sleep cycle. The estimated coefficients for the drug and drug dosage independent variables in the model provide the effect on time to REM sleep of the drugs. The drug effects ranging from most to least deleterious impact on mean percent time to REM sleep are: Cymbalta 60 mg and melatonin 6 mg, 5.16h (0.76 (SE: 0.08), (p = 4.9×10-16); Cymbalta 60 mg and temazepam 30 mg, 4.43h (0.64 (SE: 0.06), (p < 2×10-16); Cymbalta 60 mg and melatonin 3 mg, 4.33h (0.64 (SE: 0.08), (p = 1.8×10-14); Cymbalta 60 mg and temazepam 15 mg, 4.27h (0.63 (SE: 0.07), (p < 2×10-16); Cymbalta 60 mg, 4.22h (0.62 (SE: 0.04), (p < 2×10-16); Cymbalta 30 mg, 3.73h (0.54 (SE: 0.04), (p < 2×10-16); Cymbalta 30 mg and melatonin 3 mg, 3.55h (0.51 (SE: 0.04), (p < 2×10-16) and Cymbalta 30 mg and melatonin 6 mg, 3.43h (0.49 (SE: 0.06), (p = 8.1×10-15). There was no impact of the night of the study on the model. Of note, is the decrease in time to REM sleep on weekend nights.

The data shows an unequivocal Cymbalta dose-response, decreasing the time to REM sleep with decreasing Cymbalta dose as expected. Even under the least damaging drug regimen, time to first REM sleep was still delayed over 1.75 hours compared to the maximum in normal sleep architecture (3.43 hours versus 1.67 hours or 100 minutes). This delay in first REM sleep could possibly push normal REM sleep cycling into later parts of the night and interfere with the ability to naturally wake the next morning. Further, truncating REM sleep while keeping a daily work-week schedule might be expected to have additional functional and metabolic consequences.

We used the data to attempt to predict a lower Cymbalta drug dose which might not be expected to interfere with our subject’s sleep or perhaps normalize all of the percent sleep stages toward “normal” ranges (i.e., wake 5 percent; light 45–55 percent; deep 20–25 percent; REM 25 percent46) since our subject has increased REM (34 percent) and decreased light (35 percent) sleep if drug effects were accounted for. Table 5 provides the predicted values for 10 mg and 20 mg Cymbalta doses based on the fitted regression models.

Table 5. Predicted sleep architecture responses at intermediate Cymbalta doses.

Cymbalta
dose
Time to REM
hours
Percent
wake
Percent
light
Percent
deep
Percent
REM
  0 mg1.2710.535.422.334.2
10 mg2.0913.936.920.031.1
20 mg2.9117.438.517.727.9
30 mg3.7320.840.015.424.8

We note that even considering the removal of Cymbalta altogether, the percentage of the sleep time our subject was estimated to be in a ‘wake’ period as detected by the Zeo monitor is high. PLMs that tracked with Cymbalta use did decrease to less than 15 per hour during the study (see Dataset 2) which is considered to be normal and therefore not likely to be a source of confusion for the Zeo monitor since episodes of PLMs may confound time in the wake period. However, micro-arousals and unconscious wakes due to the possible presence of mild sleep apnea in our subject remained a concern and could be reflected in the sleep values we observed. A follow-up study performed to monitor a clinical intervention to correct mild sleep apnea is presented in the Supplementary Material.

Discussion

We have shown that monitoring an individual’s response to various drugs used to treat her severe sleep and sleep-related disturbances yielded important and actionable insights. For example, the subject’s sleep quality was highly compromised when taking Cymbalta at therapeutic (60 mg) and sub-therapeutic (30 mg) doses and was likely aggravated further by polypharmaceutical interventions she was prescribed. In addition, the subject’s other conditions, such as mild sleep apnea, may also have contributed to her sleep disturbances and general physical and psychological health. While sleep disruption is a common side effect of SSRIs and SNRIs, our finding that Cymbalta appears to have exacerbated the subject’s condition, is important for personalized care of patients with nuanced conditions. The problems associated with Cymbalta may have been due to the extended release formulation of the drug. It is known that Cymbalta is metabolized by CYP2D6, which has been recently shown to undergo a metabolizer phenotype conversion that cannot be assessed by genetic testing47. Drug-induced and particularly co-medication-induced phenoconversion is an increasing problem for personalized medicine48. Additionally, temazepam is not a short-acting benzodiazepine drug and can cause hangover effects in the course of a night that could contribute to the phase-delay our subject experienced. In fact, both temazepam and another highly used sleep aid, Ambien, were recently found to be associated with increased morbidity and mortality49. Despite the fact our subject was co-morbid for a number of circadian disruptors, her sleep architecture normalized when all drugs were removed. In addition, drug removal unmasked mild sleep apnea, manifesting mainly during an NREM sleep component. The temazepam-Cymbalta combination appears to have induced a removal of deep sleep that actually mimics the shallow sleep architecture seen in depressed patients50. Antidepressants are often touted as able to restore deep sleep and delay REM sleep in depression50. However, for the subject of focus here (and we suspect many others), the major destruction of her deep sleep occurred when a sleep aid was added to counteract the over-stimulation of the antidepressant.

A number of studies have shown that antidepressants can exacerbate symptoms associated with depression3037. Further, we found that our subject suffered from mild obstructive sleep apnea (OSA) and should probably never have been on sleep medication in the first place. Symptom clusters of poor sleep, migraines, and fatigue should motivate a physician to perform a sleep study. In fact, both in menopausal women and in psychiatric practice where mood and sleep disorders can show bi-directional causation, ordering sleep studies for patients has become the recommended course51,52.

Limitations

The drug withdrawal protocol for the subject discussed here ran from December to July. The days were getting longer across the time period (after winter solstice to after summer solstice) so changes in the subject’s responses to light and increased/decreased internal secretion of melatonin/serotonin could have had a beneficial influence on the direction of the changes in sleep parameters in parallel with drug removal. Alternatively, the hypersomnia expected in a SAD-susceptible individual during December-May could result in a more sound sleep (except for sleep latency issues expected from her phenotype/chronotype). However, we showed that the final (no-drug) sleep architecture in July 2013 was equivalent to that observed at the beginning of our sleep apnea intervention in December 2013 (see Supplementary Material). In the end, the subject demonstrated what is typical for SAD, normal sleep architecture, but tendency toward a delayed chronotype.

Due to the free-living nature of our study, attempts to follow/collect standardized food, exercise and sleep/wake behavior were not maintained, although, attempts to phase-shift to earlier sleep/wake regimens were documented. Applying a SPOT design, there was no randomization, drug placebo, blinding or washouts between trials, but we were able to compare our subject’s status to her status at times when no drug was provided in a crossover setting. Abrupt changes in treatment may have contributed some expected and some unexpected noise to the data. For example, temazepam dose decreases would be expected to result in delayed sleep onset, however, changes from Cymbalta 60 mg to Cymbalta 30 mg caused hot-flashes also impacting sleep initiation. For the most part, we collected enough data under each treatment studied (relative to drug or device on/off) to measure effects, including the capture of rebound and recovery effects, and the duration of our individual trial conditions were comparable to what is often seen in sleep literature. As stated in the Methods, our decision was to use a real-time/real-life dose withdrawal and not to use washout periods (the appropriate duration of “washout” would be hard to determine for Cymbalta) to avoid harms. As it was, our Cymbalta dose de-escalation was slower than what is used in clinical practice (Cymbalta 60 mg for 52 days, Cymbalta 30 mg for 84 days, Cymbalta 0 mg for 31 days). We also limited the number of times any one drug combination was provided. Given the number of drugs and the number of doses studied, it would be virtually impossible to accommodate multiple intervals with the same drugs and dosages given. Again, given the strong impact of the drugs used in this case, as evidenced by the variability in the data, only the lower dose trials included cross-overs. This resulted in trials at the beginning of the study having only one measurement, albeit covering periods from 7–25 days each. A similar range in days is seen when the duration of the individual same drug/dose regimen trial replicates are combined.

Conclusions

Many people suffering from circadian and sleep disturbances such as those found in SAD have very unique genetic determinants for their condition, different sets of sleep disturbance sequelae, secondary conditions, and nuanced lifestyles that make it hard to treat them exactly the same way. As a result, more focused attention on what intervention strategy makes the most sense to pursue is required. Such ‘personalized’ intervention strategies are not trivial to implement since they require an integrated, objective, and often-times completely empirical approach to identify and implement them. We describe our experience with, and the results of, a comprehensive investigation into the response of a single patient to designed manipulations of her sleep pharmacology. We find that the patient had underlying conditions (e.g., sleep apnea) that were confounded by the use of specific drugs to treat her SAD and that these drugs contributed to, or exacerbated, other issues in the subject’s life (e.g., alert time for work, attempts to make up for lack of quality sleep during the week on the weekends, etc.). Ultimately, our study and its results should set a precedent for patient-oriented, yet designed and objective, investigations into the impact of polypharmacy and general drug response in real-world settings.

Data availability

F1000Research: Dataset 1. Drug dosage and sleep response data, 10.5256/f1000research.7694.d11201653

F1000Research: Dataset 2. PAM-RL Periodic Leg Movement Rates, 10.5256/f1000research.7694.d13033854

Consent

Written informed consent for publication of their clinical details and/or clinical images was obtained from the patient.

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 03 Feb 2016
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Magnuson V, Wang Y and Schork N. Normalizing sleep quality disturbed by psychiatric polypharmacy: a single patient open trial (SPOT) [version 2; peer review: 2 approved]. F1000Research 2016, 5:132 (https://doi.org/10.12688/f1000research.7694.2)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 2
VERSION 2
PUBLISHED 05 Sep 2016
Revised
Views
7
Cite
Reviewer Report 25 Jul 2017
Catherine J. Nikles, Center for Clinical Research, The University of Queensland, Brisbane, QLD, Australia 
Thomas A. Harris, Center for Clinical Research, The University of Queensland, Brisbane, QLD, Australia 
Approved
VIEWS 7
Background
The area of the study is highly relevant to sleep medicine and addresses the compounding factors of polypharmacy in a single case study design. 

The pathology of SAD is well explained with all of the ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Nikles CJ and Harris TA. Reviewer Report For: Normalizing sleep quality disturbed by psychiatric polypharmacy: a single patient open trial (SPOT) [version 2; peer review: 2 approved]. F1000Research 2016, 5:132 (https://doi.org/10.5256/f1000research.9984.r24049)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
6
Cite
Reviewer Report 12 Dec 2016
Wilson D. Pace, DARTNet Institute, Aurora, CO, USA 
Approved
VIEWS 6
This version has been cleaned up considerably but is still long. Removing the sleep apnea sections and cleaning up most of the issues between background, methods, results and discussion helps improve the ability to read the document cleanly. Overall the ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Pace WD. Reviewer Report For: Normalizing sleep quality disturbed by psychiatric polypharmacy: a single patient open trial (SPOT) [version 2; peer review: 2 approved]. F1000Research 2016, 5:132 (https://doi.org/10.5256/f1000research.9984.r16094)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Version 1
VERSION 1
PUBLISHED 03 Feb 2016
Views
40
Cite
Reviewer Report 11 May 2016
Wilson D. Pace, DARTNet Institute, Aurora, CO, USA 
Not Approved
VIEWS 40
This manuscript purports to report on an N-of-1 trial to select therapy for a patient with seasonal affective disorder, sleep apnea and polypharmacy. The outcome of the activities was improved sleep quality. The condition, the treatments to be considered and ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Pace WD. Reviewer Report For: Normalizing sleep quality disturbed by psychiatric polypharmacy: a single patient open trial (SPOT) [version 2; peer review: 2 approved]. F1000Research 2016, 5:132 (https://doi.org/10.5256/f1000research.8286.r12973)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 05 Sep 2016
    Victoria Magnuson, Department of Human Biology, J Craig Venter Institute, La Jolla, USA
    05 Sep 2016
    Author Response
    We would like to thank Dr. Wilson Pace for his insightful review of our manuscript and his many constructive suggestions for improvement. We have amended our manuscript in response and ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 05 Sep 2016
    Victoria Magnuson, Department of Human Biology, J Craig Venter Institute, La Jolla, USA
    05 Sep 2016
    Author Response
    We would like to thank Dr. Wilson Pace for his insightful review of our manuscript and his many constructive suggestions for improvement. We have amended our manuscript in response and ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 03 Feb 2016
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.