Rhythmicity of neuronal oscillations delineates their cortical and spectral architecture.

Neuronal oscillations are commonly analyzed with power spectral methods that quantify signal amplitude, but not rhythmicity or 'oscillatoriness' per se. Here we introduce a new approach, the phase-autocorrelation function (pACF), for the direct quantification of rhythmicity. We applied pACF to human intracerebral stereoelectroencephalography (SEEG) and magnetoencephalography (MEG) data and uncovered a spectrally and anatomically fine-grained cortical architecture in the rhythmicity of single- and multi-frequency neuronal oscillations. Evidencing the functional significance of rhythmicity, we found it to be a prerequisite for long-range synchronization in resting-state networks and to be dynamically modulated during event-related processing. We also extended the pACF approach to measure 'burstiness' of oscillatory processes and characterized regions with stable and bursty oscillations. These findings show that rhythmicity is double-dissociable from amplitude and constitutes a functionally relevant and dynamic characteristic of neuronal oscillations.

Neuronal synchrony and critical bistability: Mechanistic biomarkers for localizing the epileptogenic network.

OBJECTIVE: Postsurgical seizure freedom in drug-resistant epilepsy (DRE) patients varies from 30% to 80%, implying that in many cases the current approaches fail to fully map the epileptogenic zone (EZ). We aimed to advance a novel approach to better characterize epileptogenicity and investigate whether the EZ encompasses a broader epileptogenic network (EpiNet) beyond the seizure zone (SZ) that exhibits seizure activity. METHODS: We first used computational modeling to test putative complex systems-driven and systems neuroscience-driven mechanistic biomarkers for epileptogenicity. We then used these biomarkers to extract features from resting-state stereoelectroencephalograms recorded from DRE patients and trained supervised classifiers to localize the SZ against gold standard clinical localization. To further explore the prevalence of pathological features in an extended brain network outside of the clinically identified SZ, we also used unsupervised classification. RESULTS: Supervised SZ classification trained on individual features achieved accuracies of .6-.7 area under the receiver operating characteristic curve (AUC). Combining all criticality and synchrony features further improved the AUC to .85. Unsupervised classification discovered an EpiNet-like cluster of brain regions, in which 51% of brain regions were outside of the SZ. Brain regions in the EpiNet-like cluster engaged in interareal hypersynchrony and locally exhibited high-amplitude bistability and excessive inhibition, which was strikingly similar to the high seizure risk regime revealed by our computational modeling. SIGNIFICANCE: The finding that combining biomarkers improves SZ localization accuracy indicates that the novel mechanistic biomarkers for epileptogenicity employed here yield synergistic information. On the other hand, the discovery of SZ-like brain dynamics outside of the clinically defined SZ provides empirical evidence of an extended pathophysiological EpiNet.

Measuring digital intervention user experience with a novel ecological momentary assessment (EMA) method, CORTO.

Digital interventions often suffer from low usage, which may reflect insufficient attention to user experience. Moreover, the existing evaluation methods have limited applicability in the remote study of user experience of complex interventions that have expansive content and that are used over an extensive period of time. To alleviate these challenges, we describe here a novel qualitative Ecological Momentary Assessment (EMA) method: the CORTO method (Contextual, One-item, Repeated, Timely, Open-ended). We used it to gather digital intervention user experience data from Finnish adults (n = 184) who lived with interview-confirmed major depressive disorder (MDD) and took part in a randomized controlled trial (RCT) that studied the efficacy of a novel 12-week game-based digital intervention for depression. A second dataset on user experience was gathered with retrospective interviews (n = 22). We inductively coded the CORTO method and retrospective interview data, which led to four user experience categories: (1) contextual use, (2) interaction-elicited emotional experience, (3) usability, and (4) technical issues. Then, we used the created user experience categories and Template Analysis to analyze both datasets together, and reported the results qualitatively. Finally, we compared the two datasets with each other. We found that the data generated with the CORTO method offered more insights into usability and technical categories than the interview data that particularly illustrated the contextual use. The emotional valence of the interview data was more positive compared with the CORTO data. Both the CORTO and interview data detected 55 % of the micro-level categories; 20 % of micro-level categories were only detected by the CORTO data and 25 % only by the interview data. We found that the during-intervention user experience measurement with the CORTO method can provide intervention-specific insights, and thereby further the iterative user-centered intervention development. Overall, these findings highlight the impact of evaluation methods on the categories and qualities of insights acquired in intervention research.

Critical-like Brain Dynamics in a Continuum from Second- to First-Order Phase Transition.

The classic brain criticality hypothesis postulates that the brain benefits from operating near a continuous second-order phase transition. Slow feedback regulation of neuronal activity could, however, lead to a discontinuous first-order transition and thereby bistable activity. Observations of bistability in awake brain activity have nonetheless remained scarce and its functional significance unclear. Moreover, there is no empirical evidence to support the hypothesis that the human brain could flexibly operate near either a first- or second-order phase transition despite such a continuum being common in models. Here, using computational modeling, we found bistable synchronization dynamics to emerge through elevated positive feedback and occur exclusively in a regimen of critical-like dynamics. We then assessed bistability in vivo with resting-state MEG in healthy adults (7 females, 11 males) and stereo-electroencephalography in epilepsy patients (28 females, 36 males). This analysis revealed that a large fraction of the neocortices exhibited varying degrees of bistability in neuronal oscillations from 3 to 200 Hz. In line with our modeling results, the neuronal bistability was positively correlated with classic assessment of brain criticality across narrow-band frequencies. Excessive bistability was predictive of epileptic pathophysiology in the patients, whereas moderate bistability was positively correlated with task performance in the healthy subjects. These empirical findings thus reveal the human brain as a one-of-a-kind complex system that exhibits critical-like dynamics in a continuum between continuous and discontinuous phase transitions.SIGNIFICANCE STATEMENT In the model, while synchrony per se was controlled by connectivity, increasing positive local feedback led to gradually emerging bistable synchrony with scale-free dynamics, suggesting a continuum between second- and first-order phase transitions in synchrony dynamics inside a critical-like regimen. In resting-state MEG and SEEG, bistability of ongoing neuronal oscillations was pervasive across brain areas and frequency bands and was observed only with concurring critical-like dynamics as the modeling predicted. As evidence for functional relevance, moderate bistability was positively correlated with executive functioning in the healthy subjects, and excessive bistability was associated with epileptic pathophysiology. These findings show that critical-like neuronal dynamics in vivo involves both continuous and discontinuous phase transitions in a frequency-, neuroanatomy-, and state-dependent manner.

The Development of Game-Based Digital Mental Health Interventions: Bridging the Paradigms of Health Care and Entertainment.

Game elements are increasingly used to improve user engagement in digital mental health interventions, and specific game mechanics may yield therapeutic effects per se and thereby contribute to digital mental health intervention efficacy. However, only a few commercial game-based interventions are available. We suggest that the key challenge in their development reflects the tension between the 2 underlying paradigms, health care and entertainment, which have disparate goals and processes in digital development. We describe 3 approaches currently used to negotiate the 2 paradigms: the gamification of health care software, designing serious games, and purpose shifting existing entertainment games. We advanced an integrative framework to focus attention on 4 key themes in intervention development: target audience, engagement, mechanisms of action, and health-related effectiveness. On each theme, we show how the 2 paradigms contrast and can complement each other. Finally, we consider the 4 interdependent themes through the new product development phases from concept to production. Our viewpoint provides an integrative synthesis that facilitates the research, design, and development of game-based digital mental health interventions.

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data.

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

Factors Affecting Digital Tool Use in Client Interaction According to Mental Health Professionals: Interview Study.

BACKGROUND: Digital tools and interventions are being increasingly developed in response to the growing mental health crisis, and mental health professionals (MHPs) considerably influence their adoption in client practice. However, how MHPs use digital tools in client interaction is yet to be sufficiently understood, which poses challenges to their design, development, and implementation. OBJECTIVE: This study aimed to create a contextual understanding of how MHPs use different digital tools in clinical client practice and what characterizes the use across tools. METHODS: A total of 19 Finnish MHPs participated in semistructured interviews, and the data were transcribed, coded, and inductively analyzed. RESULTS: We found that MHP digital tool use was characterized by 3 distinct functions: communication, diagnosis and evaluation, and facilitating therapeutic change. The functions were addressed using analog tools, digitized tools that mimic their analog counterparts, and digital tools that use the possibilities native to digital. The MHP-client communication included various media alongside face-to-face meetings, the MHPs increasingly used digitized tools in client evaluation, and the MHPs actively used digitized materials to facilitate therapeutic change. MHP tool use was generally characterized by adaptability-it was negotiated in client interactions. However, there was considerable variance in the breadth of MHPs' digital toolbox. The existing clinical practices emphasized MHP-client interaction and invited incremental rather than radical developments, which challenged the achievement of the scalability benefits expected from digital tools. CONCLUSIONS: MHPs use digitized and digital tools in client practice. Our results contribute to the user-centered research, development, and implementation of new digital solutions in mental health care by classifying them according to their function and medium and describing how MHPs use and do not use them.

Genetic polymorphisms in COMT and BDNF influence synchronization dynamics of human neuronal oscillations.

Neuronal oscillations, their inter-areal synchronization, and scale-free dynamics constitute fundamental mechanisms for cognition by regulating communication in neuronal networks. These oscillatory dynamics have large inter-individual variability that is partly heritable. We hypothesized that this variability could be partially explained by genetic polymorphisms in neuromodulatory genes. We recorded resting-state magnetoencephalography (MEG) from 82 healthy participants and investigated whether oscillation dynamics were influenced by genetic polymorphisms in catechol-O-methyltransferase (COMT) Val158Met and brain-derived neurotrophic factor (BDNF) Val66Met. Both COMT and BDNF polymorphisms influenced local oscillation amplitudes and their long-range temporal correlations (LRTCs), while only BDNF polymorphism affected the strength of large-scale synchronization. Our findings demonstrate that COMT and BDNF genetic polymorphisms contribute to inter-individual variability in neuronal oscillation dynamics. Comparison of these results to computational modeling of near-critical synchronization dynamics further suggested that COMT and BDNF polymorphisms influenced local oscillations by modulating the excitation-inhibition balance according to the brain criticality framework.

Phase-Based Cortical Synchrony Is Affected by Prematurity.

Inter-areal synchronization by phase-phase correlations (PPCs) of cortical oscillations mediates many higher neurocognitive functions, which are often affected by prematurity, a globally prominent neurodevelopmental risk factor. Here, we used electroencephalography to examine brain-wide cortical PPC networks at term-equivalent age, comparing human infants after early prematurity to a cohort of healthy controls. We found that prematurity affected these networks in a sleep state-specific manner, and the differences between groups were also frequency-selective, involving brain-wide connections. The strength of synchronization in these networks was predictive of clinical outcomes in the preterm infants. These findings show that prematurity affects PPC networks in a clinically significant manner, suggesting early functional biomarkers of later neurodevelopmental compromise that may be used in clinical or translational studies after early neonatal adversity.

Long range temporal correlations (LRTCs) in MEG-data during emerging psychosis: Relationship to symptoms, medication-status and clinical trajectory.

Long-Range Temporal Correlations (LRTCs) index the capacity of the brain to optimally process information. Previous research has shown that patients with chronic schizophrenia present altered LRTCs at alpha and beta oscillations. However, it is currently unclear at which stage of schizophrenia aberrant LRTCs emerge. To address this question, we investigated LRTCs in resting-state magnetoencephalographic (MEG) recordings obtained from patients with affective disorders and substance abuse (clinically at low-risk of psychosis, CHR-N), patients at clinical high-risk of psychosis (CHR-P) (n = 115), as well as patients with a first episode (FEP) (n = 25). Matched healthy controls (n = 47) served as comparison group. LRTCs were obtained for frequencies from 4 to 40 Hz and correlated with clinical and neuropsychological data. In addition, we examined the relationship between LRTCs and transition to psychosis in CHR-P participants, and the relationship between LRTC and antipsychotic medication in FEP participants. Our results show that participants from the clinical groups have similar LRTCs to controls. In addition, LRTCs did not correlate with clinical and neurocognitive variables across participants nor did LRTCs predict transition to psychosis. Therefore, impaired LRTCs do not reflect a feature in the clinical trajectory of psychosis. Nevertheless, reduced LRTCs in the beta-band over posterior sensors of medicated FEP participants indicate that altered LRTCs may appear at the onset of the illness. Future studies are needed to elucidate the role of anti-psychotic medication in altered LRTCs.

Detrended fluctuation analysis in the presurgical evaluation of parietal lobe epilepsy patients.

OBJECTIVE: To examine the usability of long-range temporal correlations (LRTCs) in non-invasive localization of the epileptogenic zone (EZ) in refractory parietal lobe epilepsy (RPLE) patients. METHODS: We analyzed 10 RPLE patients who had presurgical MEG and underwent epilepsy surgery. We quantified LRTCs with detrended fluctuation analysis (DFA) at four frequency bands for 200 cortical regions estimated using individual source models. We correlated individually the DFA maps to the distance from the resection area and from cortical locations of interictal epileptiform discharges (IEDs). Additionally, three clinical experts inspected the DFA maps to visually assess the most likely EZ locations. RESULTS: The DFA maps correlated with the distance to resection area in patients with type II focal cortical dysplasia (FCD) (p<0.05), but not in other etiologies. Similarly, the DFA maps correlated with the IED locations only in the FCD II patients. Visual analysis of the DFA maps showed high interobserver agreement and accuracy in FCD patients in assigning the affected hemisphere and lobe. CONCLUSIONS: Aberrant LRTCs correlate with the resection areas and IED locations. SIGNIFICANCE: This methodological pilot study demonstrates the feasibility of approximating cortical LRTCs from MEG that may aid in the EZ localization and provide new non-invasive insight into the presurgical evaluation of epilepsy.

Reduced evoked activity and cortical oscillations are correlated with anisometric amblyopia and impairment of visual acuity.

Amblyopia is a developmental disorder associated with abnormal visual experience during early childhood commonly arising from strabismus and/or anisometropia and leading to dysfunctions in visual cortex and to various visual deficits. The different forms of neuronal activity that are attenuated in amblyopia have been only partially characterized. In electrophysiological recordings of healthy human brain, the presentation of visual stimuli is associated with event-related activity and oscillatory responses. It has remained poorly understood whether these forms of activity are reduced in amblyopia and whether possible dysfunctions would arise from lower- or higher-order visual areas. We recorded neuronal activity with magnetoencephalography (MEG) from anisometropic amblyopic patients and control participants during two visual tasks presented separately for each eye and estimated neuronal activity from source-reconstructed MEG data. We investigated whether event-related and oscillatory responses would be reduced for amblyopia and localized their cortical sources. Oscillation amplitudes and evoked responses were reduced for stimuli presented to the amblyopic eye in higher-order visual areas and in parietal and prefrontal cortices. Importantly, the reduction of oscillation amplitudes but not that of evoked responses was correlated with decreased visual acuity in amblyopia. These results show that attenuated oscillatory responses are correlated with visual deficits in anisometric amblyopia.

Spectral and Anatomical Patterns of Large-Scale Synchronization Predict Human Attentional Capacity.

The capacity of visual attention determines how many visual objects may be perceived at any moment. This capacity can be investigated with multiple object tracking (MOT) tasks, which have shown that it varies greatly between individuals. The neuronal mechanisms underlying capacity limits have remained poorly understood. Phase synchronization of cortical oscillations coordinates neuronal communication within the fronto-parietal attention network and between the visual regions during endogenous visual attention. We tested a hypothesis that attentional capacity is predicted by the strength of pretarget synchronization within attention-related cortical regions. We recorded cortical activity with magneto- and electroencephalography (M/EEG) while measuring attentional capacity with MOT tasks and identified large-scale synchronized networks from source-reconstructed M/EEG data. Individual attentional capacity was correlated with load-dependent strengthening of theta (3-8 Hz), alpha (8-10 Hz), and gamma-band (30-120 Hz) synchronization that connected the visual cortex with posterior parietal and prefrontal cortices. Individual memory capacity was also preceded by crossfrequency phase-phase and phase-amplitude coupling of alpha oscillation phase with beta and gamma oscillations. Our results show that good attentional capacity is preceded by efficient dynamic functional coupling and decoupling within brain regions and across frequencies, which may enable efficient communication and routing of information between sensory and attentional systems.

Genuine cross-frequency coupling networks in human resting-state electrophysiological recordings.

Phase synchronization of neuronal oscillations in specific frequency bands coordinates anatomically distributed neuronal processing and communication. Typically, oscillations and synchronization take place concurrently in many distinct frequencies, which serve separate computational roles in cognitive functions. While within-frequency phase synchronization has been studied extensively, less is known about the mechanisms that govern neuronal processing distributed across frequencies and brain regions. Such integration of processing between frequencies could be achieved via cross-frequency coupling (CFC), either by phase-amplitude coupling (PAC) or by n:m-cross-frequency phase synchrony (CFS). So far, studies have mostly focused on local CFC in individual brain regions, whereas the presence and functional organization of CFC between brain areas have remained largely unknown. We posit that interareal CFC may be essential for large-scale coordination of neuronal activity and investigate here whether genuine CFC networks are present in human resting-state (RS) brain activity. To assess the functional organization of CFC networks, we identified brain-wide CFC networks at mesoscale resolution from stereoelectroencephalography (SEEG) and at macroscale resolution from source-reconstructed magnetoencephalography (MEG) data. We developed a novel, to our knowledge, graph-theoretical method to distinguish genuine CFC from spurious CFC that may arise from nonsinusoidal signals ubiquitous in neuronal activity. We show that genuine interareal CFC is present in human RS activity in both SEEG and MEG data. Both CFS and PAC networks coupled theta and alpha oscillations with higher frequencies in large-scale networks connecting anterior and posterior brain regions. CFS and PAC networks had distinct spectral patterns and opposing distribution of low- and high-frequency network hubs, implying that they constitute distinct CFC mechanisms. The strength of CFS networks was also predictive of cognitive performance in a separate neuropsychological assessment. In conclusion, these results provide evidence for interareal CFS and PAC being 2 distinct mechanisms for coupling oscillations across frequencies in large-scale brain networks.

Influence of multiple hypothesis testing on reproducibility in neuroimaging research: A simulation study and Python-based software.

BACKGROUND: Reproducibility of research findings has been recently questioned in many fields of science, including psychology and neurosciences. One factor influencing reproducibility is the simultaneous testing of multiple hypotheses, which entails false positive findings unless the analyzed p-values are carefully corrected. While this multiple testing problem is well known and studied, it continues to be both a theoretical and practical problem. NEW METHOD: Here we assess reproducibility in simulated experiments in the context of multiple testing. We consider methods that control either the family-wise error rate (FWER) or false discovery rate (FDR), including techniques based on random field theory (RFT), cluster-mass based permutation testing, and adaptive FDR. Several classical methods are also considered. The performance of these methods is investigated under two different models. RESULTS: We found that permutation testing is the most powerful method among the considered approaches to multiple testing, and that grouping hypotheses based on prior knowledge can improve power. We also found that emphasizing primary and follow-up studies equally produced most reproducible outcomes. COMPARISON WITH EXISTING METHOD(S): We have extended the use of two-group and separate-classes models for analyzing reproducibility and provide a new open-source software "MultiPy" for multiple hypothesis testing. CONCLUSIONS: Our simulations suggest that performing strict corrections for multiple testing is not sufficient to improve reproducibility of neuroimaging experiments. The methods are freely available as a Python toolkit "MultiPy" and we aim this study to help in improving statistical data analysis practices and to assist in conducting power and reproducibility analyses for new experiments.

Neuronal correlates of full and partial visual conscious perception.

Stimuli may induce only partial consciousness-an intermediate between null and full consciousness-where the presence but not identity of an object can be reported. The differences in the neuronal basis of full and partial consciousness are poorly understood. We investigated if evoked and oscillatory activity could dissociate full from partial conscious perception. We recorded human cortical activity with magnetoencephalography (MEG) during a visual perception task in which stimulus could be either partially or fully perceived. Partial consciousness was associated with an early increase in evoked activity and theta/low-alpha-band oscillations while full consciousness was also associated with late evoked activity and beta-band oscillations. Full from partial consciousness was dissociated by stronger evoked activity and late increase in theta oscillations that were localized to higher-order visual regions and posterior parietal and prefrontal cortices. Our results reveal both evoked activity and theta oscillations dissociate partial and full consciousness.

Preterm Birth Changes Networks of Newborn Cortical Activity.

Preterm birth is the greatest risk factor for lifelong neurocognitive deficits, globally. The effect of prematurity on early cortical network function has, however, remained poorly understood. Here, we developed a novel methodology that allows reliable assessment of functional connectivity in neonatal brain activity at millisecond and multisecond scales in terms of cortical phase and amplitude correlations, respectively. We measured scalp electroencephalography at term-equivalent age in infants exposed to very early prematurity as well as in healthy controls. We found that newborn cortical activity organizes into multiplex networks that differ significantly between vigilance states. As compared with healthy control infants, prematurity was found to cause frequency-specific patterns of dysconnectivity in cortical network, changes that were distinct for networks of phase and amplitude correlations. Neuroanatomically, the most prominent markers of prematurity were found in connections involving the frontal regions. Phase synchrony in frontally connected networks was correlated with newborn neurological performance, suggesting the first measure of cortical functional coupling that correlates with neurological performance in human infant.

Comparison of Methods to Identify Modules in Noisy or Incomplete Brain Networks.

Community structure, or "modularity," is a fundamentally important aspect in the organization of structural and functional brain networks, but their identification with community detection methods is confounded by noisy or missing connections. Although several methods have been used to account for missing data, the performance of these methods has not been compared quantitatively so far. In this study, we compared four different approaches to account for missing connections when identifying modules in binary and weighted networks using both Louvain and Infomap community detection algorithms. The four methods are "zeros," "row-column mean," "common neighbors," and "consensus clustering." Using Lancichinetti-Fortunato-Radicchi benchmark-simulated binary and weighted networks, we find that "zeros," "row-column mean," and "common neighbors" approaches perform well with both Louvain and Infomap, whereas "consensus clustering" performs well with Louvain but not Infomap. A similar pattern of results was observed with empirical networks from stereotactical electroencephalography data, except that "consensus clustering" outperforms other approaches on weighted networks with Louvain. Based on these results, we recommend any of the four methods when using Louvain on binary networks, whereas "consensus clustering" is superior with Louvain clustering of weighted networks. When using Infomap, "zeros" or "common neighbors" should be used for both binary and weighted networks. These findings provide a basis to accounting for noisy or missing connections when identifying modules in brain networks.

Dynamic large-scale network synchronization from perception to action.

Sensory-guided actions entail the processing of sensory information, generation of perceptual decisions, and the generation of appropriate actions. Neuronal activity underlying these processes is distributed into sensory, fronto-parietal, and motor brain areas, respectively. How the neuronal processing is coordinated across these brain areas to support functions from perception to action remains unknown. We investigated whether phase synchronization in large-scale networks coordinate these processes. We recorded human cortical activity with magnetoencephalography (MEG) during a task in which weak somatosensory stimuli remained unperceived or were perceived. We then assessed dynamic evolution of phase synchronization in large-scale networks from source-reconstructed MEG data by using advanced analysis approaches combined with graph theory. Here we show that perceiving and reporting of weak somatosensory stimuli is correlated with sustained strengthening of large-scale synchrony concurrently in delta/theta (3-7 Hz) and gamma (40-60 Hz) frequency bands. In a data-driven network localization, we found this synchronization to dynamically connect the task-relevant, that is, the fronto-parietal, sensory, and motor systems. The strength and temporal pattern of interareal synchronization were also correlated with the response times. These data thus show that key brain areas underlying perception, decision-making, and actions are transiently connected by large-scale dynamic phase synchronization in the delta/theta and gamma bands.