Co-occurring ripple oscillations facilitate neuronal interactions between cortical locations in humans.

How the human cortex integrates ("binds") information encoded by spatially distributed neurons remains largely unknown. One hypothesis suggests that synchronous bursts of high-frequency oscillations ("ripples") contribute to binding by facilitating integration of neuronal firing across different cortical locations. While studies have demonstrated that ripples modulate local activity in the cortex, it is not known whether their co-occurrence coordinates neural firing across larger distances. We tested this hypothesis using local field-potentials and single-unit firing from four 96-channel microelectrode arrays in the supragranular cortex of 3 patients. Neurons in co-rippling locations showed increased short-latency co-firing, prediction of each other's firing, and co-participation in neural assemblies. Effects were similar for putative pyramidal and interneurons, during non-rapid eye movement sleep and waking, in temporal and Rolandic cortices, and at distances up to 16 mm (the longest tested). Increased co-prediction during co-ripples was maintained when firing-rate changes were equated, indicating that it was not secondary to non-oscillatory activation. Co-rippling enhanced prediction was strongly modulated by ripple phase, supporting the most common posited mechanism for binding-by-synchrony. Co-ripple enhanced prediction is reciprocal, synergistic with local upstates, and further enhanced when multiple sites co-ripple, supporting re-entrant facilitation. Together, these results support the hypothesis that trans-cortical co-occurring ripples increase the integration of neuronal firing of neurons in different cortical locations and do so in part through phase-modulation rather than unstructured activation.

Binding of cortical functional modules by synchronous high frequency oscillations.

High-frequency phase-locked oscillations have been hypothesized to facilitate integration ('binding') of information encoded across widespread cortical areas. Ripples (~100ms long ~90Hz oscillations) co-occur ('co-ripple') broadly in multiple states and locations, but have only been associated with memory replay. We tested whether cortico-cortical co-ripples subserve a general role in binding by recording intracranial EEG during reading. Co-rippling increased to words versus consonant-strings between visual, wordform and semantic cortical areas when letters are binding into words, and words to meaning. Similarly, co-ripples strongly increased before correct responses between executive, response, wordform and semantic areas when word meanings bind instructions and response. Task-selective co-rippling dissociated from non-oscillatory activation and memory reinstatement. Co-ripples were phase-locked at zero-lag, even at long distances (>12cm), supporting a general role in cognitive binding.

Co-occurring ripple oscillations facilitate neuronal interactions between cortical locations in humans.

Synchronous bursts of high frequency oscillations ('ripples') are hypothesized to contribute to binding by facilitating integration of neuronal firing across cortical locations. We tested this hypothesis using local field-potentials and single-unit firing from four 96-channel microelectrode arrays in supragranular cortex of 3 patients. Neurons in co-rippling locations showed increased short-latency co-firing, prediction of each-other's firing, and co-participation in neural assemblies. Effects were similar for putative pyramidal and interneurons, during NREM sleep and waking, in temporal and Rolandic cortices, and at distances up to 16mm. Increased co-prediction during co-ripples was maintained when firing-rate changes were equated, and were strongly modulated by ripple phase. Co-ripple enhanced prediction is reciprocal, synergistic with local upstates, and further enhanced when multiple sites co-ripple. Together, these results support the hypothesis that trans-cortical co-ripples increase the integration of neuronal firing of neurons in different cortical locations, and do so in part through phase-modulation rather than unstructured activation.

Cortical Ripples during NREM Sleep and Waking in Humans.

Hippocampal ripples index the reconstruction of spatiotemporal neuronal firing patterns essential for the consolidation of memories in the cortex during non-rapid eye movement sleep (NREM). Recently, cortical ripples in humans have been shown to enfold the replay of neuron firing patterns during cued recall. Here, using intracranial recordings from 18 patients (12 female), we show that cortical ripples also occur during NREM in humans, with similar density, oscillation frequency (∼90 Hz), duration, and amplitude to waking. Ripples occurred in all cortical regions with similar characteristics, unrelated to putative hippocampal connectivity, and were less dense and robust in higher association areas. Putative pyramidal and interneuron spiking phase-locked to cortical ripples during NREM, with phase delays consistent with ripple generation through pyramidal-interneuron feedback. Cortical ripples were smaller in amplitude than hippocampal ripples but were similar in density, frequency, and duration. Cortical ripples during NREM typically occurred just before the upstate peak, often during spindles. Upstates and spindles have previously been associated with memory consolidation, and we found that cortical ripples grouped cofiring between units within the window of spike timing-dependent plasticity. Thus, human NREM cortical ripples are as follows: ubiquitous and stereotyped with a tightly focused oscillation frequency; similar to hippocampal ripples; associated with upstates and spindles; and associated with unit cofiring. These properties are consistent with cortical ripples possibly contributing to memory consolidation and other functions during NREM in humans.SIGNIFICANCE STATEMENT In rodents, hippocampal ripples organize replay during sleep to promote memory consolidation in the cortex, where ripples also occur. However, evidence for cortical ripples in human sleep is limited, and their anatomic distribution and physiological properties are unexplored. Here, using human intracranial recordings, we demonstrate that ripples occur throughout the cortex during waking and sleep with highly stereotyped characteristics. During sleep, cortical ripples tend to occur during spindles on the down-to-upstate transition, and thus participate in a sequence of sleep waves that is important for consolidation. Furthermore, cortical ripples organize single-unit spiking with timing optimal to facilitate plasticity. Therefore, cortical ripples in humans possess essential physiological properties to support memory and other cognitive functions.

Widespread ripples synchronize human cortical activity during sleep, waking, and memory recall.

Declarative memory encoding, consolidation, and retrieval require the integration of elements encoded in widespread cortical locations. The mechanism whereby such "binding" of different components of mental events into unified representations occurs is unknown. The "binding-by-synchrony" theory proposes that distributed encoding areas are bound by synchronous oscillations enabling enhanced communication. However, evidence for such oscillations is sparse. Brief high-frequency oscillations ("ripples") occur in the hippocampus and cortex and help organize memory recall and consolidation. Here, using intracranial recordings in humans, we report that these ∼70-ms-duration, 90-Hz ripples often couple (within ±500 ms), co-occur (≥ 25-ms overlap), and, crucially, phase-lock (have consistent phase lags) between widely distributed focal cortical locations during both sleep and waking, even between hemispheres. Cortical ripple co-occurrence is facilitated through activation across multiple sites, and phase locking increases with more cortical sites corippling. Ripples in all cortical areas co-occur with hippocampal ripples but do not phase-lock with them, further suggesting that cortico-cortical synchrony is mediated by cortico-cortical connections. Ripple phase lags vary across sleep nights, consistent with participation in different networks. During waking, we show that hippocampo-cortical and cortico-cortical coripples increase preceding successful delayed memory recall, when binding between the cue and response is essential. Ripples increase and phase-modulate unit firing, and coripples increase high-frequency correlations between areas, suggesting synchronized unit spiking facilitating information exchange. co-occurrence, phase synchrony, and high-frequency correlation are maintained with little decrement over very long distances (25 cm). Hippocampo-cortico-cortical coripples appear to possess the essential properties necessary to support binding by synchrony during memory retrieval and perhaps generally in cognition.

Myofiber strain in healthy humans using DENSE and cDTI.

PURPOSE: Myofiber strain, Eff , is a mechanistically relevant metric of cardiac cell shortening and is expected to be spatially uniform in healthy populations, making it a prime candidate for the evaluation of local cardiomyocyte contractility. In this study, a new, efficient pipeline was proposed to combine microstructural cDTI and functional DENSE data in order to estimate Eff in vivo. METHODS: Thirty healthy volunteers were scanned with three long-axis (LA) and three short-axis (SA) DENSE slices using 2D displacement encoding and one SA slice of cDTI. The total acquisition time was 11 minutes ± 3 minutes across volunteers. The pipeline first generates 3D SA displacements from all DENSE slices which are then combined with cDTI data to generate a cine of myofiber orientations and compute Eff . The precision of the post-processing pipeline was assessed using a computational phantom study. Transmural myofiber strain was compared to circumferential strain, Ecc , in healthy volunteers using a Wilcoxon sign rank test. RESULTS: In vivo, computed Eff was found uniform transmurally compared to Ecc (-0.14[-0.15, -0.12] vs -0.18 [-0.20, -0.16], P < .001, -0.14 [-0.16, -0.12] vs -0.16 [-0.17, -0.13], P < .001 and -0.14 [-0.16, -0.12] vs Ecc_C = -0.14 [-0.15, -0.11], P = .002, Eff_C vs Ecc_C in the endo, mid, and epi layers, respectively). CONCLUSION: We demonstrate that it is possible to measure in vivo myofiber strain in a healthy human population in 10 minutes per subject. Myofiber strain was observed to be spatially uniform in healthy volunteers making it a potential biomarker for the evaluation of local cardiomyocyte contractility in assessing cardiovascular dysfunction.

Estimating cardiomyofiber strain in vivo by solving a computational model.

Since heart contraction results from the electrically activated contraction of millions of cardiomyocytes, a measure of cardiomyocyte shortening mechanistically underlies cardiac contraction. In this work we aim to measure preferential aggregate cardiomyocyte ("myofiber") strains based on Magnetic Resonance Imaging (MRI) data acquired to measure both voxel-wise displacements through systole and myofiber orientation. In order to reduce the effect of experimental noise on the computed myofiber strains, we recast the strains calculation as the solution of a boundary value problem (BVP). This approach does not require a calibrated material model, and consequently is independent of specific myocardial material properties. The solution to this auxiliary BVP is the displacement field corresponding to assigned values of myofiber strains. The actual myofiber strains are then determined by minimizing the difference between computed and measured displacements. The approach is validated using an analytical phantom, for which the ground-truth solution is known. The method is applied to compute myofiber strains using in vivo displacement and myofiber MRI data acquired in a mid-ventricular left ventricle section in N=8 swine subjects. The proposed method shows a more physiological distribution of myofiber strains compared to standard approaches that directly differentiate the displacement field.

Probing cardiomyocyte mobility with multi-phase cardiac diffusion tensor MRI.

PURPOSE: Cardiomyocyte organization and performance underlie cardiac function, but the in vivo mobility of these cells during contraction and filling remains difficult to probe. Herein, a novel trigger delay (TD) scout sequence was used to acquire high in-plane resolution (1.6 mm) Spin-Echo (SE) cardiac diffusion tensor imaging (cDTI) at three distinct cardiac phases. The objective was to characterize cardiomyocyte organization and mobility throughout the cardiac cycle in healthy volunteers. MATERIALS AND METHODS: Nine healthy volunteers were imaged with cDTI at three distinct cardiac phases (early systole, late systole, and diastasis). The sequence used a free-breathing Spin-Echo (SE) cDTI protocol (b-values = 350s/mm2, twelve diffusion encoding directions, eight repetitions) to acquire high-resolution images (1.6x1.6x8mm3) at 3T in ~7 minutes/cardiac phase. Helix Angle (HA), Helix Angle Range (HAR), E2 angle (E2A), Transverse Angle (TA), Mean Diffusivity (MD), diffusion tensor eigenvalues (λ1-2-3), and Fractional Anisotropy (FA) in the left ventricle (LV) were characterized. RESULTS: Images from the patient-specific TD scout sequence demonstrated that SE cDTI acquisition was possible at early systole, late systole, and diastasis in 78%, 100% and 67% of the cases, respectively. At the mid-ventricular level, mobility (reported as median [IQR]) was observed in HAR between early systole and late systole (76.9 [72.6, 80.5]° vs 96.6 [85.9, 100.3]°, p<0.001). E2A also changed significantly between early systole, late systole, and diastasis (27.7 [20.8, 35.1]° vs 45.2 [42.1, 49]° vs 20.7 [16.6, 26.4]°, p<0.001). CONCLUSION: We demonstrate that it is possible to probe cardiomyocyte mobility using multi-phase and high resolution cDTI. In healthy volunteers, aggregate cardiomyocytes re-orient themselves more longitudinally during contraction, while cardiomyocyte sheetlets tilt radially during wall thickening. These observations provide new insights into the three-dimensional mobility of myocardial microstructure during systolic contraction.

Estimating Aggregate Cardiomyocyte Strain Using In Vivo Diffusion and Displacement Encoded MRI.

Changes in left ventricular (LV) aggregate cardiomyocyte orientation and deformation underlie cardiac function and dysfunction. As such, in vivo aggregate cardiomyocyte "myofiber" strain ( [Formula: see text]) has mechanistic significance, but currently there exists no established technique to measure in vivo [Formula: see text]. The objective of this work is to describe and validate a pipeline to compute in vivo [Formula: see text] from magnetic resonance imaging (MRI) data. Our pipeline integrates LV motion from multi-slice Displacement ENcoding with Stimulated Echoes (DENSE) MRI with in vivo LV microstructure from cardiac Diffusion Tensor Imaging (cDTI) data. The proposed pipeline is validated using an analytical deforming heart-like phantom. The phantom is used to evaluate 3D cardiac strains computed from a widely available, open-source DENSE Image Analysis Tool. Phantom evaluation showed that a DENSE MRI signal-to-noise ratio (SNR) ≥20 is required to compute [Formula: see text] with near-zero median strain bias and within a strain tolerance of 0.06. Circumferential and longitudinal strains are also accurately measured under the same SNR requirements, however, radial strain exhibits a median epicardial bias of -0.10 even in noise-free DENSE data. The validated framework is applied to experimental DENSE MRI and cDTI data acquired in eight ( N=8 ) healthy swine. The experimental study demonstrated that [Formula: see text] has decreased transmural variability compared to radial and circumferential strains. The spatial uniformity and mechanistic significance of in vivo [Formula: see text] make it a compelling candidate for characterization and early detection of cardiac dysfunction.

Concordance Between Electronic Clinical Documentation and Physicians' Observed Behavior.

Importance: Following the adoption of electronic health records into a regulatory environment designed for paper records, there has been little investigation into the accuracy of physician documentation. Objective: To quantify the percentage of emergency physician documentation of the review of systems (ROS) and physical examination (PE) that observers can confirm. Design, Setting, and Participants: This case series took place at emergency departments in 2 academic medical centers between 2016 and 2018. Participants' patient encounters were observed to compare real-time performance with clinical documentation. Exposures: Resident physicians were shadowed by trained observers for 20 encounters (10 encounters per physician per site) to obtain real-time observational data; associated electronic health record data were subsequently reviewed. Main Outcomes and Measures: Number of confirmed ROS systems (range, 0-14) divided by the number of documented ROS systems (range, 0-14), and number of confirmed PE systems (range, 0-14) divided by the number of documented PE systems (range, 0-14). Results: The final study cohort included 9 licensed emergency medicine residents who evaluated a total of 180 patients (mean [SD] age, 48.7 [20.0] years; 91 [50.5%] women). For ROS, physicians documented a median (interquartile range [IQR]) of 14 (8-14) systems, while audio recordings confirmed a median (IQR) of 5 (3-6) systems. Overall, 755 of 1961 documented ROS systems (38.5%) were confirmed by audio recording data. For PE, resident physicians documented a median (IQR) of 8 (7-9) verifiable systems, while observers confirmed a median (IQR) of 5.5 (3-6) systems. Overall, 760 of 1429 verifiable documented PE systems (53.2%) were confirmed by concurrent observation. Interrater reliability for rating of ROS and PE was more than 90% for all measures. Conclusions and Relevance: In this study of 9 licensed year emergency medicine residents, there were inconsistencies between the documentation of ROS and PE findings in the electronic health record and observational reports. These findings raise the possibility that some documentation may not accurately represent physician actions. Further studies should be undertaken to determine whether this occurrence is widespread. However, because such studies are unlikely to be performed owing to institution-level barriers that exist nationwide, payers should consider removing financial incentives to generate lengthy documentation.

High-Resolution Ex Vivo Microstructural MRI After Restoring Ventricular Geometry via 3D Printing.

Computational modeling of the heart requires accurately incorporating both gross anatomical detail and local microstructural information. Together, these provide the necessary data to build 3D meshes for simulation of cardiac mechanics and electrophysiology. Recent MRI advances make it possible to measure detailed heart motion in vivo, but in vivo microstructural imaging of the heart remains challenging. Consequently, the most detailed measurements of microstructural organization and microanatomical infarct details are obtained ex vivo. The objective of this work was to develop and evaluate a new method for restoring ex vivo ventricular geometry to match the in vivo configuration. This approach aids the integration of high-resolution ex vivo microstructural information with in vivo motion measurements. The method uses in vivo cine imaging to generate surface meshes, then creates a 3D printed left ventricular (LV) blood pool cast and a pericardial mold to restore the ex vivo cardiac geometry to a mid-diastasis reference configuration. The method was evaluated in healthy (N = 7) and infarcted (N = 3) swine. Dice similarity coefficients were calculated between in vivo and ex vivo images for the LV cavity (0.93 ± 0.01), right ventricle (RV) cavity (0.80 ± 0.05), and the myocardium (0.72 ± 0.04). The R 2 coefficient between in vivo and ex vivo LV and RV cavity volumes were 0.95 and 0.91, respectively. These results suggest that this method adequately restores ex vivo geometry to match in vivo geometry. This approach permits a more precise incorporation of high-resolution ex vivo anatomical and microstructural data into computational models that use in vivo data for simulation of cardiac mechanics and electrophysiology.

Model of Left Ventricular Contraction: Validation Criteria and Boundary Conditions.

Computational models of cardiac contraction can provide critical insight into cardiac function and dysfunction. A necessary step before employing these computational models is their validation. Here we propose a series of validation criteria based on left ventricular (LV) global (ejection fraction and twist) and local (strains in a cylindrical coordinate system, aggregate cardiomyocyte shortening, and low myocardial compressibility) MRI measures to characterize LV motion and deformation during contraction. These validation criteria are used to evaluate an LV finite element model built from subject-specific anatomy and aggregate cardiomyocyte orientations reconstructed from diffusion tensor MRI. We emphasize the key role of the simulation boundary conditions in approaching the physiologically correct motion and strains during contraction. We conclude by comparing the global and local validation criteria measures obtained using two different boundary conditions: the first constraining the LV base and the second taking into account the presence of the pericardium, which leads to greatly improved motion and deformation.

TIME RESOLVED DISPLACEMENT-BASED REGISTRATION OF IN VIVO CDTI CARDIOMYOCYTE ORIENTATIONS.

In vivo cardiac microstructure acquired using cardiac diffusion tensor imaging (cDTI) is a critical component of patient-specific models of cardiac electrophysiology and mechanics. In order to limit bulk motion artifacts and acquisition time, cDTI microstructural data is acquired at a single cardiac phase necessitating registration to the reference configuration on which the patient-specific computational models are based. Herein, we propose a method to register subject-specific microstructural data to an arbitrary cardiac phase using measured cardiac displacements. We validate our approach using a subject-specific computational phantom based on data from human subjects. Compared to a geometry-based non-rigid registration method, the displacement-based registration leads to improved accuracy (less than 1° versus 10° average median error in cardiomyocyte angular differences) and tighter confidence interval (3° versus 65° average upper limit of the 95% confidence interval).

Microstructurally Anchored Cardiac Kinematics by Combining In Vivo DENSE MRI and cDTI.

Metrics of regional myocardial function can detect the onset of cardiovascular disease, evaluate the response to therapy, and provide mechanistic insight into cardiac dysfunction. Knowledge of local myocardial microstructure is necessary to distinguish between isotropic and anisotropic contributions of local deformation and to quantify myofiber kinematics, a microstructurally anchored measure of cardiac function. Using a computational model we combine in vivo cardiac displacement and diffusion tensor data to evaluate pointwise the deformation gradient tensor and isotropic and anisotropic deformation invariants. In discussing the imaging methods and the model construction, we identify potential improvements to increase measurement accuracy. We conclude by demonstrating the applicability of our method to compute myofiber strain in five healthy volunteers.