PT - JOURNAL ARTICLE AU - Indie C. Garwood AU - Sourish Chakravarty AU - Jacob Donoghue AU - Pegah Kahali AU - Shubham Chamadia AU - Oluwaseun Akeju AU - Earl K. Miller AU - Emery N. Brown TI - A hidden Markov model reliably characterizes ketamine-induced spectral dynamics in macaque LFP and human EEG AID - 10.1101/2020.11.12.20221366 DP - 2020 Jan 01 TA - medRxiv PG - 2020.11.12.20221366 4099 - http://medrxiv.org/content/early/2020/11/16/2020.11.12.20221366.short 4100 - http://medrxiv.org/content/early/2020/11/16/2020.11.12.20221366.full AB - Ketamine is an NMDA receptor antagonist commonly used to maintain general anesthesia. At anesthetic doses, ketamine causes bursts of 30-50 Hz oscillations alternating with 0.1 to 10 Hz oscillations. These dynamics are readily observed in local field potentials (LFPs) of non-human primates (NHPs) and electroencephalogram (EEG) recordings from human subjects. However, a detailed statistical analysis of these dynamics has not been reported. We characterize ketamine’s neural dynamics using a hidden Markov model (HMM). The HMM observations are sequences of spectral power in 10 Hz frequency bands between 0 to 50 Hz, where power is averaged within each band and scaled between 0 and 1. We model the observations as realizations of multivariate beta probability distributions that depend on a discrete-valued latent state process whose state transitions obey Markov dynamics. Using an expectation-maximization algorithm, we fit this beta-HMM to LFP recordings from 2 NHPs, and separately, to EEG recordings from 9 human subjects who received anesthetic doses of ketamine. Together, the estimated beta-HMM parameters and optimal state trajectory revealed an alternating pattern of states characterized primarily by gamma burst and slow oscillation activity, as well as intermediate states in between. The mean duration of the gamma burst state was 2.5s([1.9,3.4]s) and 1.2s([0.9,1.5]s) for the two NHPs, and 2.7s([1.9,3.8]s) for the human subjects. The mean duration of the slow oscillation state was 1.6s([1.1,2.5]s) and 0.7s([0.6,0.9]s) for the two NHPs, and 2.8s([1.9,4.3]s) for the human subjects. Our beta-HMM framework provides a useful tool for experimental data analysis. Our characterizations of the gamma-burst process offer detailed, quantitative constraints that can inform the development of rhythm-generating neuronal circuit models that give mechanistic insights into this phenomenon and how ketamine produces altered states of arousal.Competing Interest StatementENB has licensed IP to Masimo for the analysis of EEG under anesthesia. ENB is a cofounder of PASCALL, a company constructing control systems for anesthesia.Funding StatementThis work was partially supported by NIH Award P01-GM118629 (to E. N. Brown, E. K. Miller), NIH Award R01MH115592 (to E. K. Miller), NIH Award 5T32EB019940-05 (to I. C. Garwood), NSF GRFP (to I. C. Garwood), funds from Massachusetts General Hospital (to E. N. Brown), funds from the Picower Institute for Learning and Memory (to E. N. Brown), Picower Postdoctoral Fellowship (to S. Chakravarty).Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:The Partners Institutional Review Board approved the human retrospective observational study of the clinical EEG data. Regarding the non-human primate (NHP) data, all NHP procedures reported here followed the guidelines of the Massachusetts Institute of Technology's Committee on Animal Care and of the National Institutes of Health.All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesData and code will be made publicly available when the paper is published.