Inferring cardiac phase response curve in vivo

Characterizing properties of biological oscillators with phase response cirves (PRC) is one of main theoretical tools in neuroscience, cardio-respiratory physiology, and chronobiology. We present a technique that allows the extraction of the PRC from a non-invasive observation of a system consisting of two interacting oscillators, in this case heartbeat and respiration, in its natural environment and under free-running conditions. We use this method to obtain the phase coupling functions describing cardio-respiratory interactions and the phase response curve of 17 healthy humans. We show at which phase the cardiac beat is susceptible to respiratory drive and extract the respiratory-related component of heart rate variability. This non-invasive method of bivariate data analysis for the determination of phase response curves of coupled oscillators may find application in other biological and physical systems.