State of the art of brain function detection technologies in robot-assisted lower limb rehabilitation.

BACKGROUND: With an aging population, the prevalence of neurological disorders is increasing, leading to a rise in lower limb movement disorders and, in turn, a growing need for rehabilitation training. Previous neuroimaging studies have shown a growing scientific interest in the study of brain mechanisms in robot-assisted lower limb rehabilitation (RALLR). OBJECTIVE: This review aimed to determine differences in neural activity patterns during different RALLR tasks and the impact on neurofunctional plasticity. METHODS: Sixty-five articles in the field of RALLR were selected and tested using three brain function detection technologies (BFDT). RESULTS: Most studies have focused on changes in activity in various regions of the cerebral cortex during different lower limb rehabilitation tasks, but have also increasingly focused on functional changes in other cortical and deep subcortical structures. Our analysis also revealed a neglect of certain task types. CONCLUSION: We identify and discuss future research directions that may contribute to a clear understanding of neural functional plasticity under different RALLR tasks.

Data-driven exploratory method investigation on the effect of dyslexia education at brain connectivity in Turkish children: a preliminary study.

Dyslexia is a specific learning disability that is neurobiological in origin and is characterized by reading and/or spelling problems affecting the development of language-related skills. The aim of this study is to reveal functional markers based on dyslexia by examining the functions of brain regions in resting state and reading tasks and to analyze the effects of special education given during the treatment process of dyslexia. A total of 43 children, aged between 7 and 12, whose native language was Turkish, participated in the study in three groups including those diagnosed with dyslexia for the first time, those receiving special education for dyslexia, and healthy children. Independent component analysis method was employed to analyze functional connectivity variations among three groups both at rest and during the continuous reading task. A whole-brain scanning during task fulfillment and resting states revealed that there were significant differences in the regions including lateral visual, default mode, left frontoparietal, ventral attention, orbitofrontal and lateral motor network. Our results revealed the necessity of adding motor coordination exercises to the training of dyslexic participants and showed that training led to functional connectivity in some brain regions similar to the healthy group. Additionally, our findings confirmed that impulsivity is associated with motor coordination and visuality, and that the dyslexic group has weaknesses in brain connectivity related to these conditions. According to our preliminary results, the differences obtained between children with dyslexia, group of dyslexia with special education and healthy children has revealed the effect of education on brain functions as well as enabling a comprehensive examination of dyslexia.

A widespread electrical brain network encodes anxiety in health and depressive states.

In rodents, anxiety is charactered by heightened vigilance during low-threat and uncertain situations. Though activity in the frontal cortex and limbic system are fundamental to supporting this internal state, the underlying network architecture that integrates activity across brain regions to encode anxiety across animals and paradigms remains unclear. Here, we utilize parallel electrical recordings in freely behaving mice, translational paradigms known to induce anxiety, and machine learning to discover a multi-region network that encodes the anxious brain-state. The network is composed of circuits widely implicated in anxiety behavior, it generalizes across many behavioral contexts that induce anxiety, and it fails to encode multiple behavioral contexts that do not. Strikingly, the activity of this network is also principally altered in two mouse models of depression. Thus, we establish a network-level process whereby the brain encodes anxiety in health and disease.

Cerebral cortex functional reorganization in preschool children with congenital sensorineural hearing loss: a resting-state fMRI study.

Purpose: How cortical functional reorganization occurs after hearing loss in preschool children with congenital sensorineural hearing loss (CSNHL) is poorly understood. Therefore, we used resting-state functional MRI (rs-fMRI) to explore the characteristics of cortical reorganization in these patents. Methods: Sixty-three preschool children with CSNHL and 32 healthy controls (HCs) were recruited, and the Categories of Auditory Performance (CAP) scores were determined at the 6-month follow-up after cochlear implantation (CI). First, rs-fMRI data were preprocessed, and amplitude of low-frequency fluctuation (ALFF) and regional homogeneity (ReHo) were calculated. Second, whole-brain functional connectivity (FC) analysis was performed using bilateral primary auditory cortex as seed points. Finally, Spearman correlation analysis was performed between the differential ALFF, ReHo and FC values and the CAP score. Results: ALFF analysis showed that preschool children with CSNHL had lower ALFF values in the bilateral prefrontal cortex and superior temporal gyrus than HCs, but higher ALFF values in the bilateral thalamus and calcarine gyrus. And correlation analysis showed that some abnormal brain regions were weak negatively correlated with CAP score (p < 0.05). The ReHo values in the bilateral superior temporal gyrus, part of the prefrontal cortex and left insular gyrus were lower, whereas ReHo values in the bilateral thalamus, right caudate nucleus and right precentral gyrus were higher, in children with CSNHL than HCs. However, there was no correlation between ReHo values and the CAP scores (p < 0.05). Using primary auditory cortex (PAC) as seed-based FC further analysis revealed enhanced FC in the visual cortex, proprioceptive cortex and motor cortex. And there were weak negative correlations between the FC values in the bilateral superior temporal gyrus, occipital lobe, left postcentral gyrus and right thalamus were weakly negatively correlated and the CAP score (p < 0.05). Conclusion: After auditory deprivation in preschool children with CSNHL, the local functions of auditory cortex, visual cortex, prefrontal cortex and somatic motor cortex are changed, and the prefrontal cortex plays a regulatory role in this process. There is functional reorganization or compensation between children's hearing and these areas, which may not be conducive to auditory language recovery after CI in deaf children.

Individual-based morphological brain network changes in children with Rolandic epilepsy.

OBJECTIVE: To investigate the local cortical morphology and individual-based morphological brain networks (MBNs) changes in children with Rolandic epilepsy (RE). METHODS: Based on the structural MRI data of 56 children with RE and 56 healthy controls (HC), we constructed four types of individual-based MBNs using morphological indices (cortical thickness [CT], fractal dimension [FD], gyrification index [GI], and sulcal depth [SD]). The global and nodal properties of the brain networks were analyzed using graph theory. The between-group difference in local morphology and network topology was estimated, and partial correlation analysis was further analyzed. RESULTS: Compared with the HC, children with RE showed regional GI increases in the right posterior cingulate gyrus and SD increases in the right anterior cingulate gyrus and medial prefrontal cortex. Regarding the network level, RE exhibited increased characteristic path length in CT-based and FD-based networks, while decreased FD-based network node efficiency in the right inferior frontal gyrus. No significant correlation between altered morphological features and clinical variables was found in RE. CONCLUSIONS: These findings indicated that children with RE have disrupted morphological brain network organization beyond local morphology changes. SIGNIFICANCE: The present study could provide more theoretical basis for exploring the neuropathological mechanisms in RE.

Intentional binding for noninvasive BCI control.

Noninvasive brain-computer interfaces (BCIs) allow to interact with the external environment by naturally bypassing the musculoskeletal system. Making BCIs efficient and accurate is paramount to improve the reliability of real-life and clinical applications, from open-loop device control to closed-loop neurorehabilitation. By promoting sense of agency and embodiment, realistic setups including multimodal channels of communication, such as eye-gaze, and robotic prostheses aim to improve BCI performance. However, how the mental imagery command should be integrated in those hybrid systems so as to ensure the best interaction is still poorly understood. To address this question, we performed a hybrid EEG-based BCI training involving healthy volunteers enrolled in a reach-and-grasp action operated by a robotic arm. Main results showed that the hand grasping motor imagery timing significantly affects the BCI accuracy evolution as well as the spatiotemporal brain dynamics. Larger accuracy improvement was obtained when motor imagery is performed just after the robot reaching, as compared to before or during the movement. The proximity with the subsequent robot grasping favored intentional binding, led to stronger motor-related brain activity, and primed the ability of sensorimotor areas to integrate information from regions implicated in higher-order cognitive functions. Taken together, these findings provided fresh evidence about the effects of intentional binding on human behavior and cortical network dynamics that can be exploited to design a new generation of efficient brain-machine interfaces.

Serine racemase deletion alters adolescent social behavior and whole-brain cFos activation.

Background: Neurodevelopmental disorders (NDDs) can cause debilitating impairments in social cognition and aberrant functional connectivity in large-scale brain networks, leading to social isolation and diminished everyday functioning. To facilitate the treatment of social impairments, animal models of NDDs that link N- methyl-D-aspartate receptor (NMDAR) hypofunction to social deficits in adulthood have been used. However, understanding the etiology of social impairments in NDDs requires investigating social changes during sensitive windows during development. Methods: We examine social behavior during adolescence using a translational mouse model of NMDAR hypofunction (SR-/-) caused by knocking out serine racemase (SR), the enzyme needed to make D-serine, a key NMDAR coagonist. Species-typical social interactions are maintained through brain-wide neural activation patterns; therefore, we employed whole-brain cFos activity mapping to examine network-level connectivity changes caused by SR deletion. Results: In adolescent SR-/- mice, we observed disinhibited social behavior toward a novel conspecific and rapid social habituation toward familiar social partners. SR-/- mice also spent more time in the open arm of the elevated plus maze which classically points to an anxiolytic behavioral phenotype. These behavioral findings point to a generalized reduction in anxiety-like behavior in both social and non-social contexts in SR-/- mice; importantly, these findings were not associated with diminished working memory. Inter-regional patterns of cFos activation revealed greater connectivity and network density in SR-/- mice compared to controls. Discussion: These results suggest that NMDAR hypofunction - a potential biomarker for NDDs - can lead to generalized behavioral disinhibition in adolescence, potentially arising from disrupted communication between and within salience and default mode networks.

Alterations in neural circuit dynamics between the limbic network and prefrontal/default mode network in patients with generalized anxiety disorder.

BACKGROUND: Widespread functional alterations have been implicated in patients with generalized anxiety disorder (GAD). However, most studies have primarily focused on static brain network features in patients with GAD. The current research focused on exploring the dynamics within functional brain networks among individuals diagnosed with GAD. METHODS: Seventy-five participants were divided into patients with GAD and healthy controls (HCs), and resting-state functional magnetic resonance imaging data were collected. The severity of symptoms was measured using the Hamilton Anxiety Scale and the Patient Health Questionnaire. Co-activation pattern (CAP) analysis, centered on the bed nucleus of the stria terminalis, was applied to explore network dynamics. The capability of these dynamic characteristics to distinguish between patients with GAD and HCs was evaluated using a support vector machine. RESULTS: Patients with GAD exhibited disruptions in the limbic-prefrontal and limbic-default-mode network circuits. Particularly noteworthy was the marked reduction in dynamic indicators such as occurrence, EntriesFromBaseline, ExitsToBaseline, in-degree, out-degree, and resilience. Moreover, these decreased dynamic features effectively distinguished the GAD group from the HC in this study. CONCLUSIONS: The current findings revealed the underlying brain networks associated with compromised emotion regulation in individuals with GAD. The dynamic reduction in connectivity between the limbic-default mode network and limbic-prefrontal networks could potentially act as a biomarker and therapeutic target for GAD in the future.

On null models for temporal small-worldness in brain dynamics.

Brain dynamics can be modeled as a temporal brain network starting from the activity of different brain regions in functional magnetic resonance imaging (fMRI) signals. When validating hypotheses about temporal networks, it is important to use an appropriate statistical null model that shares some features with the treated empirical data. The purpose of this work is to contribute to the theory of temporal null models for brain networks by introducing the random temporal hyperbolic (RTH) graph model, an extension of the random hyperbolic (RH) graph, known in the study of complex networks for its ability to reproduce crucial properties of real-world networks. We focus on temporal small-worldness which, in the static case, has been extensively studied in real-world complex networks and has been linked to the ability of brain networks to efficiently exchange information. We compare the RTH graph model with standard null models for temporal networks and show it is the null model that best reproduces the small-worldness of resting brain activity. This ability to reproduce fundamental features of real brain networks, while adding only a single parameter compared with classical models, suggests that the RTH graph model is a promising tool for validating hypotheses about temporal brain networks.

We show that the random temporal hyperbolic (RTH) graph is a suitable null model for testing hypotheses about brain dynamics, after comparing it with the current state of the art and two other geometric null models. The static version of this theoretical model captures properties of various real-world networks, and its temporal version exhibits the temporal small-world property, for which we propose a new proper temporal definition. In particular, we show that the model best reproduces the temporal small-worldness measured in the empirical temporal network extracted from fMRI signals.

Response inhibition in premotor cortex corresponds to a complex reshuffle of the mesoscopic information network.

Recent studies have explored functional and effective neural networks in animal models; however, the dynamics of information propagation among functional modules under cognitive control remain largely unknown. Here, we addressed the issue using transfer entropy and graph theory methods on mesoscopic neural activities recorded in the dorsal premotor cortex of rhesus monkeys. We focused our study on the decision time of a Stop-signal task, looking for patterns in the network configuration that could influence motor plan maturation when the Stop signal is provided. When comparing trials with successful inhibition to those with generated movement, the nodes of the network resulted organized into four clusters, hierarchically arranged, and distinctly involved in information transfer. Interestingly, the hierarchies and the strength of information transmission between clusters varied throughout the task, distinguishing between generated movements and canceled ones and corresponding to measurable levels of network complexity. Our results suggest a putative mechanism for motor inhibition in premotor cortex: a topological reshuffle of the information exchanged among ensembles of neurons.

In this study, we investigated the dynamics of information transfer among functionally identified neural modules during cognitive motor control. Our focus was on mesoscopic neural activities in the dorsal premotor cortex of rhesus monkeys engaged in a Stop-signal task. Leveraging multivariate transfer entropy and graph theory, we uncovered insights on how behavioral control shapes the topology of information transmission in a local brain network. Task phases modulated the strength and hierarchy of information exchange between modules, revealing the nuanced interplay between neural populations during generated and canceled movements. Notably, during successful inhibition, the network displayed a distinctive configuration, unveiling a novel mechanism for motor inhibition in the premotor cortex: a topological reshuffle of information among neuronal ensembles.

Association of physical function with connectivity in the sensorimotor and dorsal attention networks: why examining specific components of physical function matters.

Declining physical function with aging is associated with structural and functional brain network organization. Gaining a greater understanding of network associations may be useful for targeting interventions that are designed to slow or prevent such decline. Our previous work demonstrated that the Short Physical Performance Battery (eSPPB) score and body mass index (BMI) exhibited a statistical interaction in their associations with connectivity in the sensorimotor cortex (SMN) and the dorsal attention network (DAN). The current study examined if components of the eSPPB have unique associations with these brain networks. Functional magnetic resonance imaging was performed on 192 participants in the BNET study, a longitudinal and observational trial of community-dwelling adults aged 70 or older. Functional brain networks were generated for resting state and during a motor imagery task. Regression analyses were performed between eSPPB component scores (gait speed, complex gait speed, static balance, and lower extremity strength) and BMI with SMN and DAN connectivity. Gait speed, complex gait speed, and lower extremity strength significantly interacted with BMI in their association with SMN at rest. Gait speed and complex gait speed were interacted with BMI in the DAN at rest while complex gait speed, static balance, and lower extremity strength interacted with BMI in the DAN during motor imagery. Results demonstrate that different components of physical function, such as balance or gait speed and BMI, are associated with unique aspects of brain network organization. Gaining a greater mechanistic understanding of the associations between low physical function, body mass, and brain physiology may lead to the development of treatments that not only target specific physical function limitations but also specific brain networks.

A synergistic workspace for human consciousness revealed by Integrated Information Decomposition.

How is the information-processing architecture of the human brain organised, and how does its organisation support consciousness? Here, we combine network science and a rigorous information-theoretic notion of synergy to delineate a 'synergistic global workspace', comprising gateway regions that gather synergistic information from specialised modules across the human brain. This information is then integrated within the workspace and widely distributed via broadcaster regions. Through functional MRI analysis, we show that gateway regions of the synergistic workspace correspond to the human brain's default mode network, whereas broadcasters coincide with the executive control network. We find that loss of consciousness due to general anaesthesia or disorders of consciousness corresponds to diminished ability of the synergistic workspace to integrate information, which is restored upon recovery. Thus, loss of consciousness coincides with a breakdown of information integration within the synergistic workspace of the human brain. This work contributes to conceptual and empirical reconciliation between two prominent scientific theories of consciousness, the Global Neuronal Workspace and Integrated Information Theory, while also advancing our understanding of how the human brain supports consciousness through the synergistic integration of information.

The human brain consists of billions of neurons which process sensory inputs, such as sight and sound, and combines them with information already stored in the brain. This integration of information guides our decisions, thoughts, and movements, and is hypothesized to be integral to consciousness. However, it is poorly understood how the brain regions responsible for processing this integration are organized in the brain. To investigate this question, Luppi et al. employed a mathematical framework called Partial Information Decomposition (PID) which can distinguish different types of information: redundancy (available from many regions) and synergy (which reflects genuine integration). The team applied the PID framework to the brain scans of 100 individuals. This allowed them to identify which brain regions combine information from across the brain (known as gateways), and which ones transmit it back to the rest of the brain (known as broadcasters). Next, Luppi et al. set out to find how these regions compared in unconscious and conscious individuals. To do this, they studied 15 healthy volunteers whose brains were scanned (using a technique called functional MRI) before, during, and after anaesthesia. This revealed that the brain integrated less information when unconscious, and that this reduction happens predominantly in gateway rather than broadcaster regions. The same effect was also observed in the brains of individuals who were permanently unconscious due to brain injuries. These findings provide a way of understanding how information is organised in the brain. They also suggest that loss of consciousness due to brain injuries and anaesthesia involve similar brain circuits. This means it may be possible to gain insights about disorders of consciousness from studying how people emerge from anaesthesia.

Exploring mechanisms of anhedonia in depression through neuroimaging and data-driven approaches.

BACKGROUND: Anhedonia is a core symptom of depression that is closely related to prognosis and treatment outcomes. However, accurate and efficient treatments for anhedonia are lacking, mandating a deeper understanding of the underlying mechanisms. METHODS: A total of 303 patients diagnosed with depression and anhedonia were assessed by the Snaith-Hamilton Pleasure Scale (SHAPS) and magnetic resonance imaging (MRI). The patients were categorized into a low-anhedonia group and a high-anhedonia group using the K-means algorithm. A data-driven approach was used to explore the differences in brain structure and function with different degrees of anhedonia based on MATLAB. A random forest model was used exploratorily to test the predictive ability of differences in brain structure and function on anhedonia in depression. RESULTS: Structural and functional differences were apparent in several brain regions of patients with depression and high-level anhedonia, including in the temporal lobe, paracingulate gyrus, superior frontal gyrus, inferior occipital gyrus, right insular gyrus, and superior parietal lobule. And changes in these brain regions were significantly correlated with scores of SHAPS. CONCLUSIONS: These brain regions may be useful as biomarkers that provide a more objective assessment of anhedonia in depression, laying the foundation for precision medicine in this treatment-resistant, relatively poor prognosis group.

Forensic Neurology and the Role of Neurologists in Forensic Evaluations.

Forensic psychiatrists may be asked to opine on neurological evidence or neurological diseases outside the scope of their expertise. This article discusses the value of involving experts trained in behavioral neurology in such cases. First, we describe the field of behavioral neurology and neuropsychiatry, the subspecialty available to both neurologists and psychiatrists focused on the behavioral, cognitive, and neuropsychiatric manifestations of neurological diseases. Next, we discuss the added value of behavioral neurologists in forensic cases, including assisting in the diagnostic evaluation for complex neuropsychiatric diseases, using expertise in localization to provide a strong scientific basis for linking neurodiagnostic testing to relevant neuropsychiatric symptoms, and assisting in relating these symptoms to the relevant legal question in cases where such symptoms may be less familiar to forensic psychiatrists, such as frontal lobe syndromes. We discuss approaches to integrating behavioral neurology with forensic psychiatry, highlighting the need for collaboration and mentorship between disciplines. Finally, we discuss several forensic cases highlighting the additional value of experts trained in behavioral neurology. We conclude that forensic psychiatrists should involve behavioral neurology experts when encountering neurological evidence that falls outside their scope of expertise, and the need for further cross-disciplinary collaboration and training.

Individual Variability in Brain Connectivity Patterns and Driving-Fatigue Dynamics.

Mental fatigue during driving poses significant risks to road safety, necessitating accurate assessment methods to mitigate potential hazards. This study explores the impact of individual variability in brain networks on driving fatigue assessment, hypothesizing that subject-specific connectivity patterns play a pivotal role in understanding fatigue dynamics. By conducting a linear regression analysis of subject-specific brain networks in different frequency bands, this research aims to elucidate the relationships between frequency-specific connectivity patterns and driving fatigue. As such, an EEG sustained driving simulation experiment was carried out, estimating individuals' brain networks using the Phase Lag Index (PLI) to capture shared connectivity patterns. The results unveiled notable variability in connectivity patterns across frequency bands, with the alpha band exhibiting heightened sensitivity to driving fatigue. Individualized connectivity analysis underscored the complexity of fatigue assessment and the potential for personalized approaches. These findings emphasize the importance of subject-specific brain networks in comprehending fatigue dynamics, while providing sensor space minimization, advocating for the development of efficient mobile sensor applications for real-time fatigue detection in driving scenarios.

Network Analysis in Multiple Sclerosis and Related Disorders.

Multiple sclerosis (MS) is a neuroinflammatory and neurodegenerative disease of the central nervous system, commonly featuring disability and cognitive impairment. The pathologic hallmark of MS lies in demyelination and hence impaired structural and functional neuronal pathways. Recent studies have shown that MS shows extensive structural disconnection of key network hub areas like the thalamus, combined with a functional network reorganization that can mostly be related to poorer clinical functioning. As MS can, therefore, be considered a network disorder, this review outlines recent innovations in the field of network neuroscience in MS.

Epilepsy alters brain networks in patients with insular glioma.

AIMS: We intend to elucidate the alterations of cerebral networks in patients with insular glioma-related epilepsy (GRE) based on resting-state functional magnetic resonance images. METHODS: We collected 62 insular glioma patients, who were subsequently categorized into glioma-related epilepsy (GRE) and glioma with no epilepsy (GnE) groups, and recruited 16 healthy individuals matched to the patient's age and gender to form the healthy control (HC) group. Graph theoretical analysis was applied to reveal differences in sensorimotor, default mode, visual, and executive networks among different subgroups. RESULTS: No significant alterations in functional connectivity were found in either hemisphere insular glioma. Using graph theoretical analysis, differences were found in visual, sensorimotor, and default mode networks (p < 0.05). When the glioma located in the left hemisphere, the degree centrality was reduced in the GE group compared to the GnE group. When the glioma located in the right insula, the degree centrality, nodal efficiency, nodal local efficiency, and nodal clustering coefficient of the GE group were lower than those of the GnE group. CONCLUSION: The impact of insular glioma itself and GRE on the brain network is widespread. The networks altered by insular GRE differ depending on the hemisphere location. GRE reduces the nodal properties of brain networks than that in insular glioma.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Age-related differences in resting-state, task-related, and structural brain connectivity: graph theoretical analyses and visual search performance.

Previous magnetic resonance imaging (MRI) research suggests that aging is associated with a decrease in the functional interconnections within and between groups of locally organized brain regions (modules). Further, this age-related decrease in the segregation of modules appears to be more pronounced for a task, relative to a resting state, reflecting the integration of functional modules and attentional allocation necessary to support task performance. Here, using graph-theoretical analyses, we investigated age-related differences in a whole-brain measure of module connectivity, system segregation, for 68 healthy, community-dwelling individuals 18-78 years of age. We obtained resting-state, task-related (visual search), and structural (diffusion-weighted) MRI data. Using a parcellation of modules derived from the participants' resting-state functional MRI data, we demonstrated that the decrease in system segregation from rest to task (i.e., reconfiguration) increased with age, suggesting an age-related increase in the integration of modules required by the attentional demands of visual search. Structural system segregation increased with age, reflecting weaker connectivity both within and between modules. Functional and structural system segregation had qualitatively different influences on age-related decline in visual search performance. Functional system segregation (and reconfiguration) influenced age-related decline in the rate of visual evidence accumulation (drift rate), whereas structural system segregation contributed to age-related slowing of encoding and response processes (nondecision time). The age-related differences in the functional system segregation measures, however, were relatively independent of those associated with structural connectivity.

A Bayesian multilevel model for populations of networks using exponential-family random graphs.

The collection of data on populations of networks is becoming increasingly common, where each data point can be seen as a realisation of a network-valued random variable. Moreover, each data point may be accompanied by some additional covariate information and one may be interested in assessing the effect of these covariates on network structure within the population. A canonical example is that of brain networks: a typical neuroimaging study collects one or more brain scans across multiple individuals, each of which can be modelled as a network with nodes corresponding to distinct brain regions and edges corresponding to structural or functional connections between these regions. Most statistical network models, however, were originally proposed to describe a single underlying relational structure, although recent years have seen a drive to extend these models to populations of networks. Here, we describe a model for when the outcome of interest is a network-valued random variable whose distribution is given by an exponential random graph model. To perform inference, we implement an exchange-within-Gibbs MCMC algorithm that generates samples from the doubly-intractable posterior. To illustrate this approach, we use it to assess population-level variations in networks derived from fMRI scans, enabling the inference of age- and intelligence-related differences in the topological structure of the brain's functional connectivity. Supplementary Information: The online version contains supplementary material available at 10.1007/s11222-024-10446-0.