Neuroscience meets behavior: A systematic literature review on magnetic resonance imaging of the brain combined with real-world digital phenotyping.

A primary goal of neuroscience is to understand the relationship between the brain and behavior. While magnetic resonance imaging (MRI) examines brain structure and function under controlled conditions, digital phenotyping via portable automatic devices (PAD) quantifies behavior in real-world settings. Combining these two technologies may bridge the gap between brain imaging, physiology, and real-time behavior, enhancing the generalizability of laboratory and clinical findings. However, the use of MRI and data from PADs outside the MRI scanner remains underexplored. Herein, we present a Preferred Reporting Items for Systematic Reviews and Meta-Analysis systematic literature review that identifies and analyzes the current state of research on the integration of brain MRI and PADs. PubMed and Scopus were automatically searched using keywords covering various MRI techniques and PADs. Abstracts were screened to only include articles that collected MRI brain data and PAD data outside the laboratory environment. Full-text screening was then conducted to ensure included articles combined quantitative data from MRI with data from PADs, yielding 94 selected papers for a total of N = 14,778 subjects. Results were reported as cross-frequency tables between brain imaging and behavior sampling methods and patterns were identified through network analysis. Furthermore, brain maps reported in the studies were synthesized according to the measurement modalities that were used. Results demonstrate the feasibility of integrating MRI and PADs across various study designs, patient and control populations, and age groups. The majority of published literature combines functional, T1-weighted, and diffusion weighted MRI with physical activity sensors, ecological momentary assessment via PADs, and sleep. The literature further highlights specific brain regions frequently correlated with distinct MRI-PAD combinations. These combinations enable in-depth studies on how physiology, brain function and behavior influence each other. Our review highlights the potential for constructing brain-behavior models that extend beyond the scanner and into real-world contexts.

FreeSurfer 7 quality control: Key problem areas and importance of manual corrections.

We have studied the effects of manual quality control of brain Magnetic Resonance Imaging (MRI) images processed with Freesurfer. T1 images of first episode psychosis patients (N = 60) and healthy controls (N = 41) were inspected for gray matter boundary errors. The errors were fixed, and the effects of error correction on brain volume, thickness, and surface area were measured. It is commonplace to apply quality control to Freesurfer MRI recordings to ensure that the edges of gray and white matter are detected properly, as incorrect edge detection leads to changes in variables such as volume, cortical thickness, and cortical surface area. We find that while Freesurfer v7.1.1. does regularly make mistakes in identifying the edges of cortical gray matter, correcting these errors yields limited changes in the commonly measured variables listed above. We further find that the software makes fewer gray matter boundary errors when processing female brains. The results suggest that manually correcting gray matter boundary errors may not be worthwhile due to its small effect on the measurements, with potential exceptions for studies that focus on the areas that are more commonly affected by errors: the areas around the cerebellar tentorium, paracentral lobule, and the optic nerves, specifically the horizontal segment of the middle cerebral artery.

Peripheral differentiation patterns of human T cells.

Long-term T-cell memory is dependent on the maintenance of memory T cells in the lymphoid tissues, and at the surface interfaces that provide entry routes for pathogens. However, much of the current information on human T-cell memory is based on analyzing circulating T cells. Here, we have studied the distribution and age-related changes of memory T-cell subsets in samples from blood, mesenteric LNs, spleen, and ileum, obtained from donors ranging in age from 5 days to 67 years of age. Our data show that the main reservoir of polyclonal naive cells is found in the LNs, and the resting memory subsets capable of self-renewal are also prominent there. In contrast, nondividing but functionally active memory subsets dominate the spleen, and especially the ileum. In general, the replacement of naive cells with memory subsets continues throughout our period of observation, with no apparent plateau. In conclusion, the analysis of lymphoid and nonlymphoid tissues reveals a dynamic pattern of changes distinct to each tissue, and with substantial differences between CD4+ and CD8+ compartments.

Critical Drift in a Neuro-Inspired Adaptive Network.

It has been postulated that the brain operates in a self-organized critical state that brings multiple benefits, such as optimal sensitivity to input. Thus far, self-organized criticality has typically been depicted as a one-dimensional process, where one parameter is tuned to a critical value. However, the number of adjustable parameters in the brain is vast, and hence critical states can be expected to occupy a high-dimensional manifold inside a high-dimensional parameter space. Here, we show that adaptation rules inspired by homeostatic plasticity drive a neuro-inspired network to drift on a critical manifold, where the system is poised between inactivity and persistent activity. During the drift, global network parameters continue to change while the system remains at criticality.

Identifying the inheritable component of human thymic T cell repertoire generation in monozygous twins.

We have analyzed T cell receptor repertoires in a unique set of thymus samples from a pair of monozygotic twins. While genetics affect the V(D)J rearrangement and generation of junctional sequences, the thymic selections seem largely stochastic and import no detectable inheritable effect at clonal level.

Generation of self-reactive, shared T-cell receptor α chains in the human thymus.

The T-cell receptor (TCR) repertoire is generated in a semistochastic process of gene recombination and pairing of TCRα to TCRß chains with the estimated total TCR diversity of >108. Despite this high diversity, similar or identical TCR chains are found to recur in immune responses. Here, we analyzed the thymic generation of TCR sequences previously associated with recognition of self- and nonself-antigens, represented by sequences associated with autoimmune diabetes and HIV, respectively. Unexpectedly, in the CD4+ compartment TCRα chains associated with the recognition of self-antigens were generated in significantly higher numbers than TCRα chains associated with the recognition of nonself-antigens. The analysis of the circulating repertoire further showed that these chains are not lost in negative selection nor predominantly converted to the regulatory T-cell lineage. The high abundance of self-reactive TCRα chains in multiple individuals suggests that the human thymus has a predilection to generate self-reactive TCRα chains independently of the HLA-type and that the individual risk of autoimmunity may be modulated by the TCRß repertoire associated with these chains.

Effects of spatial smoothing on functional brain networks.

Graph-theoretical methods have rapidly become a standard tool in studies of the structure and function of the human brain. Whereas the structural connectome can be fairly straightforwardly mapped onto a complex network, there are more degrees of freedom in constructing networks that represent functional connections between brain areas. For functional magnetic resonance imaging (fMRI) data, such networks are typically built by aggregating the blood-oxygen-level dependent signal time series of voxels into larger entities (such as Regions of Interest in some brain atlas) and determining their connection strengths from some measure of time-series correlations. Although it is evident that the outcome must be affected by how the voxel-level time series are treated at the preprocessing stage, there is a lack of systematic studies of the effects of preprocessing on network structure. Here, we focus on the effects of spatial smoothing, a standard preprocessing method for fMRI. We apply various levels of spatial smoothing to resting-state fMRI data and measure the changes induced in functional networks. We show that the level of spatial smoothing clearly affects the degrees and other centrality measures of functional network nodes; these changes are non-uniform, systematic, and depend on the geometry of the brain. The composition of the largest connected network component is also affected in a way that artificially increases the similarity of the networks of different subjects. Our conclusion is that wherever possible, spatial smoothing should be avoided when preprocessing fMRI data for network analysis.

Persistence of social signatures in human communication.

The social network maintained by a focal individual, or ego, is intrinsically dynamic and typically exhibits some turnover in membership over time as personal circumstances change. However, the consequences of such changes on the distribution of an ego's network ties are not well understood. Here we use a unique 18-mo dataset that combines mobile phone calls and survey data to track changes in the ego networks and communication patterns of students making the transition from school to university or work. Our analysis reveals that individuals display a distinctive and robust social signature, captured by how interactions are distributed across different alters. Notably, for a given ego, these social signatures tend to persist over time, despite considerable turnover in the identity of alters in the ego network. Thus, as new network members are added, some old network members either are replaced or receive fewer calls, preserving the overall distribution of calls across network members. This is likely to reflect the consequences of finite resources such as the time available for communication, the cognitive and emotional effort required to sustain close relationships, and the ability to make emotional investments.

Graph coarse-graining reveals differences in the module-level structure of functional brain networks.

Networks have become a standard tool for analyzing functional magnetic resonance imaging (fMRI) data. In this approach, brain areas and their functional connections are mapped to the nodes and links of a network. Even though this mapping reduces the complexity of the underlying data, it remains challenging to understand the structure of the resulting networks due to the large number of nodes and links. One solution is to partition networks into modules and then investigate the modules' composition and relationship with brain functioning. While this approach works well for single networks, understanding differences between two networks by comparing their partitions is difficult and alternative approaches are thus necessary. To this end, we present a coarse-graining framework that uses a single set of data-driven modules as a frame of reference, enabling one to zoom out from the node- and link-level details. As a result, differences in the module-level connectivity can be understood in a transparent, statistically verifiable manner. We demonstrate the feasibility of the method by applying it to networks constructed from fMRI data recorded from 13 healthy subjects during rest and movie viewing. While independently partitioning the rest and movie networks is shown to yield little insight, the coarse-graining framework enables one to pinpoint differences in the module-level structure, such as the increased number of intra-module links within the visual cortex during movie viewing. In addition to quantifying differences due to external stimuli, the approach could also be applied in clinical settings, such as comparing patients with healthy controls.

Reorganization of functionally connected brain subnetworks in high-functioning autism.

Previous functional connectivity studies have found both hypo- and hyper-connectivity in brains of individuals having autism spectrum disorder (ASD). Here we studied abnormalities in functional brain subnetworks in high-functioning individuals with ASD during free viewing of a movie containing social cues and interactions. Twenty-six subjects (13 with ASD) watched a 68-min movie during functional magnetic resonance imaging. For each subject, we computed Pearson's correlation between haemodynamic time-courses of each pair of 6-mm isotropic voxels. From the whole-brain functional networks, we derived individual and group-level subnetworks using graph theory. Scaled inclusivity was then calculated between all subject pairs to estimate intersubject similarity of connectivity structure of each subnetwork. Additional 54 individuals (27 with ASD) from the ABIDE resting-state database were included to test the reproducibility of the results. Between-group differences were observed in the composition of default-mode and ventro-temporal-limbic (VTL) subnetworks. The VTL subnetwork included amygdala, striatum, thalamus, parahippocampal, fusiform, and inferior temporal gyri. Further, VTL subnetwork similarity between subject pairs correlated significantly with similarity of symptom gravity measured with autism quotient. This correlation was observed also within the controls, and in the reproducibility dataset with ADI-R and ADOS scores. Our results highlight how the reorganization of functional subnetworks in individuals with ASD clarifies the mixture of hypo- and hyper-connectivity findings. Importantly, only the functional organization of the VTL subnetwork emerges as a marker of inter-individual similarities that co-vary with behavioral measures across all participants. These findings suggest a pivotal role of ventro-temporal and limbic systems in autism.

Genetic structure of native ant supercolonies varies in space and time.

Ant supercolonies are the largest cooperative units known in nature. They consist of networks of interconnected nests with hundreds of reproductive queens, where individuals move freely between nests, cooperate across nest boundaries and show little aggression towards non-nestmates. The combination of high queen numbers and free mixing of workers, queens and brood between nests results in extremely low nestmate relatedness. In such low-relatedness societies, cooperative worker behaviour appears maladaptive because it may aid random individuals instead of relatives. Here, we provide a comprehensive picture of genetic substructure in supercolonies of the native wood ant Formica aquilonia using traditional population genetic as well as network analysis methods. Specifically, we test for spatial and temporal variation in genetic structure of different classes of individuals within supercolonies and analyse the role of worker movement in determining supercolony genetic networks. We find that relatedness within supercolonies is low but positive when viewed on a population level, which may be due to limited dispersal of individuals and/or ecological factors such as nest site limitation and competition against conspecifics. Genetic structure of supercolonies varied with both sample class and sampling time point, which indicates that mobility of individuals varies according to both caste and season and suggests that generalizing has to be carried out with caution in studies of supercolonial species. Overall, our analysis provides novel evidence that native wood ant supercolonies exhibit fine-scale genetic substructure, which may explain the maintenance of cooperation in these low-relatedness societies.

Temporal motifs reveal homophily, gender-specific patterns, and group talk in call sequences.

Recent studies on electronic communication records have shown that human communication has complex temporal structure. We study how communication patterns that involve multiple individuals are affected by attributes such as sex and age. To this end, we represent the communication records as a colored temporal network where node color is used to represent individuals' attributes, and identify patterns known as temporal motifs. We then construct a null model for the occurrence of temporal motifs that takes into account the interaction frequencies and connectivity between nodes of different colors. This null model allows us to detect significant patterns in call sequences that cannot be observed in a static network that uses interaction frequencies as link weights. We find sex-related differences in communication patterns in a large dataset of mobile phone records and show the existence of temporal homophily, the tendency of similar individuals to participate in communication patterns beyond what would be expected on the basis of their average interaction frequencies. We also show that temporal patterns differ between dense and sparse neighborhoods in the network. Because also this result is independent of interaction frequencies, it can be seen as an extension of Granovetter's hypothesis to temporal networks.

Adding network structure onto the map of collective behavior.

We propose an extension to the map of Bentley et al. by incorporating an aspect of underlying network structure that is likely relevant for many modes of collective behavior. This dimension, which captures a feature of network community structure, is known both from theory and from experiments to be relevant for decision-making processes.

Distinguishing subsampled power laws from other heavy-tailed distributions.

Distinguishing power-law distributions from other heavy-tailed distributions is challenging, and this task is often further complicated by subsampling effects. In this work, we evaluate the performance of two commonly used methods for detecting power-law distributions-the maximum likelihood method of Clauset et al. and the extreme value method of Voitalov et al.-in distinguishing subsampled power laws from two other heavy-tailed distributions, the lognormal and the stretched exponential distributions. We focus on a random subsampling method commonly applied in network science and biological sciences. In this subsampling scheme, we are ultimately interested in the frequency distribution of elements with a certain number of constituent parts-for example, species with k individuals or nodes with k connections-and each part is selected to the subsample with an equal probability. We investigate how well the results obtained from low-subsampling-depth subsamples generalize to the original distribution. Our results show that the power-law exponent of the original distribution can be estimated fairly accurately from subsamples, but classifying the distribution correctly is more challenging. The maximum likelihood method falsely rejects the power-law hypothesis for a large fraction of subsamples from power-law distributions. While the extreme value method correctly recognizes subsampled power-law distributions with all tested subsampling depths, its capacity to distinguish power laws from the heavy-tailed alternatives is limited. However, these false positives tend to result not from the subsampling itself but from the estimators' inability to classify the original sample correctly. In fact, we show that the extreme value method can sometimes be expected to perform better on subsamples than on the original samples from the lognormal and the stretched exponential distributions, while the contrary is true for the main tests included in the maximum likelihood method.

Universal patterns in egocentric communication networks.

Tie strengths in social networks are heterogeneous, with strong and weak ties playing different roles at the network and individual levels. Egocentric networks, networks of relationships around an individual, exhibit few strong ties and more weaker ties, as evidenced by electronic communication records. Mobile phone data has also revealed persistent individual differences within this pattern. However, the generality and driving mechanisms of social tie strength heterogeneity remain unclear. Here, we study tie strengths in egocentric networks across multiple datasets of interactions between millions of people during months to years. We find universality in tie strength distributions and their individual-level variation across communication modes, even in channels not reflecting offline social relationships. Via a simple model of egocentric network evolution, we show that the observed universality arises from the competition between cumulative advantage and random choice, two tie reinforcement mechanisms whose balance determines the diversity of tie strengths. Our results provide insight into the driving mechanisms of tie strength heterogeneity in social networks and have implications for the understanding of social network structure and individual behavior.

Analysis of thymic generation of shared T-cell receptor α repertoire associated with recognition of tumor antigens shows no preference for neoantigens over wild-type antigens.

BACKGROUND: The number of mutations in cancer cells is an important predictor of a positive response to cancer immunotherapy. It has been suggested that the neoantigens produced by these mutations are more immunogenic than nonmutated tumor antigens, which are likely to be protected by immunological tolerance. However, the mechanisms of tolerance as regards tumor antigens are incompletely understood. METHODS: Here, we have analyzed the impact of thymic negative selection on shared T-cell receptor (TCR) repertoire associated with the recognition of either mutated or nonmutated tumor antigens by comparing previously known TCR-antigen-pairs to TCR repertoires of 21 immunologically healthy individuals. RESULTS: Our results show that TCRα chains associated with either type of tumor antigens are readily generated in the thymus, at a frequency similar to TCRα chains associated with nonself. In the peripheral repertoire, the relative clone size of nonself-associated chains is higher than that of the tumor antigens, but importantly, there is no difference between TCRα chains associated with mutated or nonmutated tumor antigens. CONCLUSION: This suggests that the tolerance mechanisms protecting nonmutated tumor antigens are non-deletional and therefore potentially reversible. As unmutated antigens are, unlike mutations, shared by a large number of patients, they may offer advantages in designing immunological approaches to cancer treatment.

Herd immunity and epidemic size in networks with vaccination homophily.

We study how the herd immunity threshold and the expected epidemic size depend on homophily with respect to vaccine adoption. We find that the presence of homophily considerably increases the critical vaccine coverage needed for herd immunity and that strong homophily can push the threshold entirely out of reach. The epidemic size monotonically increases as a function of homophily strength for a perfect vaccine, while it is maximized at a nontrivial level of homophily when the vaccine efficacy is limited. Our results highlight the importance of vaccination homophily in epidemic modeling.

Quantifying daily rhythms with non-negative matrix factorization applied to mobile phone data.

Human activities follow daily, weekly, and seasonal rhythms. The emergence of these rhythms is related to physiology and natural cycles as well as social constructs. The human body and its biological functions undergo near 24-h rhythms (circadian rhythms). While their frequencies are similar across people, their phases differ. In the chronobiology literature, people are categorized into morning-type, evening-type, and intermediate-type groups called chronotypes based on their tendency to sleep at different times of day. Typically, this typology builds on carefully designed questionnaires or manually crafted features of time series data on people's activity. Here, we introduce a method where time-stamped data from smartphones are decomposed into components using non-negative matrix factorization. The method does not require any predetermined assumptions about the typical times of sleep or activity: the results are fully context-dependent and determined by the most prominent features of the activity data. We demonstrate our method by applying it to a dataset of mobile phone screen usage logs of 400 university students, collected over a year. We find four emergent temporal components: morning activity, night activity, evening activity and activity at noon. Individual behavior can be reduced to weights on these four components. We do not observe any clear categories of people based on the weights, but individuals are rather placed on a continuous spectrum according to the timings of their phone activities. High weights for the morning and night components strongly correlate with sleep and wake-up times. Our work points towards a data-driven way of characterizing people based on their full daily and weekly rhythms of activity and behavior, instead of only focusing on the timing of their sleeping periods.

Mobility Signatures: A Tool for Characterizing Cities Using Intercity Mobility Flows.

Understanding the patterns of human mobility between cities has various applications from transport engineering to spatial modeling of the spreading of contagious diseases. We adopt a city-centric, data-driven perspective to quantify such patterns and introduce the mobility signature as a tool for understanding how a city (or a region) is embedded in the wider mobility network. We demonstrate the potential of the mobility signature approach through two applications that build on mobile-phone-based data from Finland. First, we use mobility signatures to show that the well-known radiation model is more accurate for mobility flows associated with larger Finnish cities, while the traditional gravity model appears a better fit for less populated areas. Second, we illustrate how the SARS-CoV-2 pandemic disrupted the mobility patterns in Finland in the spring of 2020. These two cases demonstrate the ability of the mobility signatures to quickly capture features of mobility flows that are harder to extract using more traditional methods.

Characterization of human T cell receptor repertoire data in eight thymus samples and four related blood samples.

T cell receptor (TCR) is a heterodimer consisting of TCRα and TCRß chains that are generated by somatic recombination of multiple gene segments. Nascent TCR repertoire undergoes thymic selections where non-functional and potentially autoreactive receptors are removed. During the last years, the development of high-throughput sequencing technology has allowed a large scale assessment of TCR repertoire and multiple analysis tools are now also available. In our recent manuscript, Human thymic T cell repertoire is imprinted with strong convergence to shared sequences[1], we show highly overlapping thymic TCR repertoires in unrelated individuals. In the current Data in Brief article, we provide a more detailed characterization of the basic features of these thymic and related peripheral blood TCR repertoires. The thymus samples were collected from eight infants undergoing corrective cardiac surgery, two of whom were monozygous twins [2]. In parallel with the surgery, a small aliquot of peripheral blood was drawn from four of the donors. Genomic DNA was extracted from mechanically released thymocytes and circulating leukocytes. The sequencing of TCRα and TCRß repertoires was performed at ImmunoSEQ platform (Adaptive Biotechnologies). The obtained repertoire data were analysed applying relevant features from immunoSEQ® 3.0 Analyzer (Adaptive Biotechnologies) and a freely available VDJTools software package for programming language R [3]. The current data analysis displays the basic features of the sequenced repertoires including observed TCR diversity, various descriptive TCR diversity measures, and V and J gene usage. In addition, multiple methods to calculate repertoire overlap between two individuals are applied. The raw sequence data provide a large database of reference TCRs in healthy individuals at an early developmental stage. The data can be exploited to improve existing computational models on TCR repertoire behaviour as well as in the generation of new models.