CARLA: Adjusted common average referencing for cortico-cortical evoked potential data.

Human brain connectivity can be mapped by single pulse electrical stimulation during intracranial EEG measurements. The raw cortico-cortical evoked potentials (CCEP) are often contaminated by noise. Common average referencing (CAR) removes common noise and preserves response shapes but can introduce bias from responsive channels. We address this issue with an adjusted, adaptive CAR algorithm termed "CAR by Least Anticorrelation (CARLA)". CARLA was tested on simulated CCEP data and real CCEP data collected from four human participants. In CARLA, the channels are ordered by increasing mean cross-trial covariance, and iteratively added to the common average until anticorrelation between any single channel and all re-referenced channels reaches a minimum, as a measure of shared noise. We simulated CCEP data with true responses in 0-45 of 50 total channels. We quantified CARLA's error and found that it erroneously included 0 (median) truly responsive channels in the common average with ≤42 responsive channels, and erroneously excluded ≤2.5 (median) unresponsive channels at all responsiveness levels. On real CCEP data, signal quality was quantified with the mean R2 between all pairs of channels, which represents inter-channel dependency and is low for well-referenced data. CARLA re-referencing produced significantly lower mean R2 than standard CAR, CAR using a fixed bottom quartile of channels by covariance, and no re-referencing. CARLA minimizes bias in re-referenced CCEP data by adaptively selecting the optimal subset of non-responsive channels. It showed high specificity and sensitivity on simulated CCEP data and lowered inter-channel dependency compared to CAR on real CCEP data.

Thalamic stimulation induced changes in effective connectivity.

Deep brain stimulation (DBS) is a viable treatment for a variety of neurological conditions, however, the mechanisms through which DBS modulates large-scale brain networks are unresolved. Clinical effects of DBS are observed over multiple timescales. In some conditions, such as Parkinson's disease and essential tremor, clinical improvement is observed within seconds. In many other conditions, such as epilepsy, central pain, dystonia, neuropsychiatric conditions or Tourette syndrome, the DBS related effects are believed to require neuroplasticity or reorganization and often take hours to months to observe. To optimize DBS parameters, it is therefore essential to develop electrophysiological biomarkers that characterize whether DBS settings are successfully engaging and modulating the network involved in the disease of interest. In this study, 10 individuals with drug resistant epilepsy undergoing intracranial stereotactic EEG including a thalamus electrode underwent a trial of repetitive thalamic stimulation. We evaluated thalamocortical effective connectivity using single pulse electrical stimulation, both at baseline and following a 145 Hz stimulation treatment trial. We found that when high frequency stimulation was delivered for >1.5 hours, the evoked potentials measured from remote regions were significantly reduced in amplitude and the degree of modulation was proportional to the strength of baseline connectivity. When stimulation was delivered for shorter time periods, results were more variable. These findings suggest that changes in effective connectivity in the network targeted with DBS accumulate over hours of DBS. Stimulation evoked potentials provide an electrophysiological biomarker that allows for efficient data-driven characterization of neuromodulation effects, which could enable new objective approaches for individualized DBS optimization.

Signatures of Electrical Stimulation Driven Network Interactions in the Human Limbic System.

Stimulation-evoked signals are starting to be used as biomarkers to indicate the state and health of brain networks. The human limbic network, often targeted for brain stimulation therapy, is involved in emotion and memory processing. Previous anatomic, neurophysiological, and functional studies suggest distinct subsystems within the limbic network (Rolls, 2015). Studies using intracranial electrical stimulation, however, have emphasized the similarities of the evoked waveforms across the limbic network. We test whether these subsystems have distinct stimulation-driven signatures. In eight patients (four male, four female) with drug-resistant epilepsy, we stimulated the limbic system with single-pulse electrical stimulation. Reliable corticocortical evoked potentials (CCEPs) were measured between hippocampus and the posterior cingulate cortex (PCC) and between the amygdala and the anterior cingulate cortex (ACC). However, the CCEP waveform in the PCC after hippocampal stimulation showed a unique and reliable morphology, which we term the "limbic Hippocampus-Anterior nucleus of the thalamus-Posterior cingulate, HAP-wave." This limbic HAP-wave was visually distinct and separately decoded from the CCEP waveform in ACC after amygdala stimulation. Diffusion MRI data show that the measured end points in the PCC overlap with the end points of the parolfactory cingulum bundle rather than the parahippocampal cingulum, suggesting that the limbic HAP-wave may travel through fornix, mammillary bodies, and the anterior nucleus of the thalamus (ANT). This was further confirmed by stimulating the ANT, which evoked the same limbic HAP-wave but with an earlier latency. Limbic subsystems have unique stimulation-evoked signatures that may be used in the future to help network pathology diagnosis.SIGNIFICANCE STATEMENT The limbic system is often compromised in diverse clinical conditions, such as epilepsy or Alzheimer's disease, and characterizing its typical circuit responses may provide diagnostic insight. Stimulation-evoked waveforms have been used in the motor system to diagnose circuit pathology. We translate this framework to limbic subsystems using human intracranial stereo EEG (sEEG) recordings that measure deeper brain areas. Our sEEG recordings describe a stimulation-evoked waveform characteristic to the memory and spatial subsystem of the limbic network that we term the "limbic HAP-wave." The limbic HAP-wave follows anatomic white matter pathways from hippocampus to thalamus to the posterior cingulum and shows promise as a distinct biomarker of signaling in the human brain memory and spatial limbic network.

Canonical Response Parameterization: Quantifying the structure of responses to single-pulse intracranial electrical brain stimulation.

Single-pulse electrical stimulation in the nervous system, often called cortico-cortical evoked potential (CCEP) measurement, is an important technique to understand how brain regions interact with one another. Voltages are measured from implanted electrodes in one brain area while stimulating another with brief current impulses separated by several seconds. Historically, researchers have tried to understand the significance of evoked voltage polyphasic deflections by visual inspection, but no general-purpose tool has emerged to understand their shapes or describe them mathematically. We describe and illustrate a new technique to parameterize brain stimulation data, where voltage response traces are projected into one another using a semi-normalized dot product. The length of timepoints from stimulation included in the dot product is varied to obtain a temporal profile of structural significance, and the peak of the profile uniquely identifies the duration of the response. Using linear kernel PCA, a canonical response shape is obtained over this duration, and then single-trial traces are parameterized as a projection of this canonical shape with a residual term. Such parameterization allows for dissimilar trace shapes from different brain areas to be directly compared by quantifying cross-projection magnitudes, response duration, canonical shape projection amplitudes, signal-to-noise ratios, explained variance, and statistical significance. Artifactual trials are automatically identified by outliers in sub-distributions of cross-projection magnitude, and rejected. This technique, which we call "Canonical Response Parameterization" (CRP) dramatically simplifies the study of CCEP shapes, and may also be applied in a wide range of other settings involving event-triggered data.

Electrical Stimulation of Temporal and Limbic Circuitry Produces Distinct Responses in Human Ventral Temporal Cortex.

The human ventral temporal cortex (VTC) is highly connected to integrate visual perceptual inputs with feedback from cognitive and emotional networks. In this study, we used electrical brain stimulation to understand how different inputs from multiple brain regions drive unique electrophysiological responses in the VTC. We recorded intracranial EEG data in 5 patients (3 female) implanted with intracranial electrodes for epilepsy surgery evaluation. Pairs of electrodes were stimulated with single-pulse electrical stimulation, and corticocortical evoked potential responses were measured at electrodes in the collateral sulcus and lateral occipitotemporal sulcus of the VTC. Using a novel unsupervised machine learning method, we uncovered 2-4 distinct response shapes, termed basis profile curves (BPCs), at each measurement electrode in the 11-500 ms after stimulation interval. Corticocortical evoked potentials of unique shape and high amplitude were elicited following stimulation of several regions and classified into a set of four consensus BPCs across subjects. One of the consensus BPCs was primarily elicited by stimulation of the hippocampus; another by stimulation of the amygdala; a third by stimulation of lateral cortical sites, such as the middle temporal gyrus; and the final one by stimulation of multiple distributed sites. Stimulation also produced sustained high-frequency power decreases and low-frequency power increases that spanned multiple BPC categories. Characterizing distinct shapes in stimulation responses provides a novel description of connectivity to the VTC and reveals significant differences in input from cortical and limbic structures.SIGNIFICANCE STATEMENT Disentangling the numerous input influences on highly connected areas in the brain is a critical step toward understanding how brain networks work together to coordinate human behavior. Single-pulse electrical stimulation is an effective tool to accomplish this goal because the shapes and amplitudes of signals recorded from electrodes are informative of the synaptic physiology of the stimulation-driven inputs. We focused on targets in the ventral temporal cortex, an area strongly implicated in visual object perception. By using a data-driven clustering algorithm, we identified anatomic regions with distinct input connectivity profiles to the ventral temporal cortex. Examining high-frequency power changes revealed possible modulation of excitability at the recording site induced by electrical stimulation of connected regions.

A motor association area in the depths of the central sulcus.

Cells in the precentral gyrus directly send signals to the periphery to generate movement and are principally organized as a topological map of the body. We find that movement-induced electrophysiological responses from depth electrodes extend this map three-dimensionally throughout the gyrus. Unexpectedly, this organization is interrupted by a previously undescribed motor association area in the depths of the midlateral aspect of the central sulcus. This 'Rolandic motor association' (RMA) area is active during movements of different body parts from both sides of the body and may be important for coordinating complex behaviors.

A canonical visualization tool for SEEG electrodes.

Stereoencephalographic (SEEG) electrodes are clinically implanted into the brains of patients with refractory epilepsy to locate foci of seizure onset. They are increasingly used in neurophysiology research to determine focal human brain activity in response to tasks or stimuli. Clear visualization of SEEG electrode location with respect to patient anatomy on magnetic resonance image (MRI) scan is vital to neuroscientific understanding. An intuitive way to accomplish this is to plot brain activity and labels at electrode locations on closest MRI slices along the canonical axial, coronal, and sagittal planes. Therefore, we've developed an open-source software tool in Matlab for visualizing SEEG electrode positions, determined from computed tomography (CT), onto canonical planes of resliced brain MRI. The code and graphical user interface are available at: https://github.com/MultimodalNeuroimagingLab/mnl_seegviewClinical Relevance- This tool enables precise communication of SEEG electrode activity and location by visualization on slices of MRI in canonical axial, coronal, and sagittal planes.

Gut microbial determinants of clinically important improvement in patients with rheumatoid arthritis.

BACKGROUND: Rapid advances in the past decade have shown that dysbiosis of the gut microbiome is a key hallmark of rheumatoid arthritis (RA). Yet, the relationship between the gut microbiome and clinical improvement in RA disease activity remains unclear. In this study, we explored the gut microbiome of patients with RA to identify features that are associated with, as well as predictive of, minimum clinically important improvement (MCII) in disease activity. METHODS: We conducted a retrospective, observational cohort study on patients diagnosed with RA between 1988 and 2014. Whole metagenome shotgun sequencing was performed on 64 stool samples, which were collected from 32 patients with RA at two separate time-points approximately 6-12 months apart. The Clinical Disease Activity Index (CDAI) of each patient was measured at both time-points to assess achievement of MCII; depending on this clinical status, patients were distinguished into two groups: MCII+ (who achieved MCII; n = 12) and MCII- (who did not achieve MCII; n = 20). Multiple linear regression models were used to identify microbial taxa and biochemical pathways associated with MCII while controlling for potentially confounding factors. Lastly, a deep-learning neural network was trained upon gut microbiome, clinical, and demographic data at baseline to classify patients according to MCII status, thereby enabling the prediction of whether a patient will achieve MCII at follow-up. RESULTS: We found age to be the largest determinant of the overall compositional variance in the gut microbiome (R2 = 7.7%, P = 0.001, PERMANOVA). Interestingly, the next factor identified to explain the most variance in the gut microbiome was MCII status (R2 = 3.8%, P = 0.005). Additionally, by looking at patients' baseline gut microbiome profiles, we observed significantly different microbiome traits between patients who eventually showed MCII and those who did not. Taxonomic features include alpha- and beta-diversity measures, as well as several microbial taxa, such as Coprococcus, Bilophila sp. 4_1_30, and Eubacterium sp. 3_1_31. Notably, patients who achieved clinical improvement had higher alpha-diversity in their gut microbiomes at both baseline and follow-up visits. Functional profiling identified fifteen biochemical pathways, most of which were involved in the biosynthesis of L-arginine, L-methionine, and tetrahydrofolate, to be differentially abundant between the MCII patient groups. Moreover, MCII+ and MCII- groups showed significantly different fold-changes (from baseline to follow-up) in eight microbial taxa and in seven biochemical pathways. These results could suggest that, depending on the clinical course, gut microbiomes not only start at different ecological states, but also are on separate trajectories. Finally, the neural network proved to be highly effective in predicting which patients will achieve MCII (balanced accuracy = 90.0%, leave-one-out cross-validation), demonstrating potential clinical utility of gut microbiome profiles. CONCLUSIONS: Our findings confirm the presence of taxonomic and functional signatures of the gut microbiome associated with MCII in RA patients. Ultimately, modifying the gut microbiome to enhance clinical outcome may hold promise as a future treatment for RA.

Plasma metabolomic profiling in patients with rheumatoid arthritis identifies biochemical features predictive of quantitative disease activity.

BACKGROUND: Rheumatoid arthritis (RA) is a chronic, autoimmune disorder characterized by joint inflammation and pain. In patients with RA, metabolomic approaches, i.e., high-throughput profiling of small-molecule metabolites, on plasma or serum has thus far enabled the discovery of biomarkers for clinical subgroups, risk factors, and predictors of treatment response. Despite these recent advancements, the identification of blood metabolites that reflect quantitative disease activity remains an important challenge in precision medicine for RA. Herein, we use global plasma metabolomic profiling analyses to detect metabolites associated with, and predictive of, quantitative disease activity in patients with RA. METHODS: Ultra-high-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) was performed on a discovery cohort consisting of 128 plasma samples from 64 RA patients and on a validation cohort of 12 samples from 12 patients. The resulting metabolomic profiles were analyzed with two different strategies to find metabolites associated with RA disease activity defined by the Disease Activity Score-28 using C-reactive protein (DAS28-CRP). More specifically, mixed-effects regression models were used to identify metabolites differentially abundant between two disease activity groups ("lower", DAS28-CRP ≤ 3.2; and "higher", DAS28-CRP > 3.2) and to identify metabolites significantly associated with DAS28-CRP scores. A generalized linear model (GLM) was then constructed for estimating DAS28-CRP using plasma metabolite abundances. Finally, for associating metabolites with CRP (an indicator of inflammation), metabolites differentially abundant between two patient groups ("low-CRP", CRP ≤ 3.0 mg/L; "high-CRP", CRP > 3.0 mg/L) were investigated. RESULTS: We identified 33 metabolites differentially abundant between the lower and higher disease activity groups (P < 0.05). Additionally, we identified 51 metabolites associated with DAS28-CRP (P < 0.05). A GLM based upon these 51 metabolites resulted in higher prediction accuracy (mean absolute error [MAE] ± SD: 1.51 ± 1.77) compared to a GLM without feature selection (MAE ± SD: 2.02 ± 2.21). The predictive value of this feature set was further demonstrated on a validation cohort of twelve plasma samples, wherein we observed a stronger correlation between predicted and actual DAS28-CRP (with feature selection: Spearman's ρ = 0.69, 95% CI: [0.18, 0.90]; without feature selection: Spearman's ρ = 0.18, 95% CI: [-0.44, 0.68]). Lastly, among all identified metabolites, the abundances of eight were significantly associated with the CRP patient groups while controlling for potential confounders (P < 0.05). CONCLUSIONS: We demonstrate for the first time the prediction of quantitative disease activity in RA using plasma metabolomes. The metabolites identified herein provide insight into circulating pro-/anti-inflammatory metabolic signatures that reflect disease activity and inflammatory status in RA patients.

A subset of topologically associating domains fold into mesoscale core-periphery networks.

Mammalian genomes are folded into a hierarchy of compartments, topologically associating domains (TADs), subTADs, and long-range looping interactions. The higher-order folding patterns of chromatin contacts within TADs and how they localize to disease-associated single nucleotide variants (daSNVs) remains an open area of investigation. Here, we analyze high-resolution Hi-C data with graph theory to understand possible mesoscale network architecture within chromatin domains. We identify a subset of TADs exhibiting strong core-periphery mesoscale structure in embryonic stem cells, neural progenitor cells, and cortical neurons. Hyper-connected core nodes co-localize with genomic segments engaged in multiple looping interactions and enriched for occupancy of the architectural protein CCCTC binding protein (CTCF). CTCF knockdown and in silico deletion of CTCF-bound core nodes disrupts core-periphery structure, whereas in silico mutation of cell type-specific enhancer or gene nodes has a negligible effect. Importantly, neuropsychiatric daSNVs are significantly more likely to localize with TADs folded into core-periphery networks compared to domains devoid of such structure. Together, our results reveal that a subset of TADs encompasses looping interactions connected into a core-periphery mesoscale network. We hypothesize that daSNVs in the periphery of genome folding networks might preserve global nuclear architecture but cause local topological and functional disruptions contributing to human disease. By contrast, daSNVs co-localized with hyper-connected core nodes might cause severe topological and functional disruptions. Overall, these findings shed new light into the mesoscale network structure of fine scale genome folding within chromatin domains and its link to common genetic variants in human disease.

Detecting hierarchical genome folding with network modularity.

Mammalian genomes are folded in a hierarchy of compartments, topologically associating domains (TADs), subTADs and looping interactions. Here, we describe 3DNetMod, a graph theory-based method for sensitive and accurate detection of chromatin domains across length scales in Hi-C data. We identify nested, partially overlapping TADs and subTADs genome wide by optimizing network modularity and varying a single resolution parameter. 3DNetMod can be applied broadly to understand genome reconfiguration in development and disease.