Skip to main content

ORIGINAL RESEARCH article

Front. Syst. Neurosci., 18 November 2021
This article is part of the Research Topic Cognitive Neurodynamics: Integrating Cognitive Science and Brain Dynamics View all 15 articles

Neuropeptide Y Variation Is Associated With Altered Static and Dynamic Functional Connectivity of the Salience Network

\r\nKatherine G. WarthenKatherine G. Warthen1Robert C. Welsh*Robert C. Welsh2*Benjamin SanfordBenjamin Sanford3Vincent KoppelmansVincent Koppelmans2Margit BurmeisterMargit Burmeister4Brian J. MickeyBrian J. Mickey2
  • 1Department of Biomedical Engineering, The University of Utah, Salt Lake City, UT, United States
  • 2Department of Psychiatry, The University of Utah, Salt Lake City, UT, United States
  • 3Department of Psychiatry, University of Michigan, Michigan, MI, United States
  • 4Michigan Neuroscience Institute and Departments of Computational Medicine & Bioinformatics, Human Genetics and Psychiatry, The University of Michigan, Michigan, MI, United States

Neuropeptide Y (NPY) is a neurotransmitter that has been implicated in the development of anxiety and mood disorders. Low levels of NPY have been associated with risk for these disorders, and high levels with resilience. Anxiety and depression are associated with altered intrinsic functional connectivity of brain networks, but the effect of NPY on functional connectivity is not known. Here, we test the hypothesis that individual differences in NPY expression affect resting functional connectivity of the default mode and salience networks. We evaluated static connectivity using graph theoretical techniques and dynamic connectivity with Leading Eigenvector Dynamics Analysis (LEiDA). To increase our power of detecting NPY effects, we genotyped 221 individuals and identified 29 healthy subjects at the extremes of genetically predicted NPY expression (12 high, 17 low). Static connectivity analysis revealed that lower levels of NPY were associated with shorter path lengths, higher global efficiency, higher clustering, higher small-worldness, and average higher node strength within the salience network, whereas subjects with high NPY expression displayed higher modularity and node eccentricity within the salience network. Dynamic connectivity analysis showed that the salience network of low-NPY subjects spent more time in a highly coordinated state relative to high-NPY subjects, and the salience network of high-NPY subjects switched between states more frequently. No group differences were found for static or dynamic connectivity of the default mode network. These findings suggest that genetically driven individual differences in NPY expression influence risk of mood and anxiety disorders by altering the intrinsic functional connectivity of the salience network.

Introduction

Neuropeptide Y (NPY) is among the most abundantly expressed peptides in the brain (Tatemoto et al., 1982; Adrian et al., 1983; Allen et al., 1983), and is known to affect the neural processing of emotion, appetite, and stress (Morgan et al., 2002; Van Den Heuvel et al., 2015; Reichmann and Holzer, 2016). Individuals vary in their expression of NPY, and low levels of NPY have been associated with the development of anxiety (Reichmann and Holzer, 2016), depression (Widdowson et al., 1992; Heilig et al., 2004; Mickey et al., 2011), and posttraumatic stress disorder (PTSD) (Yehuda et al., 2006). We recently showed that healthy individuals with genetically driven low NPY expression exhibit heightened activation of the nucleus accumbens in response to salient stimuli (Warthen et al., 2018). There is also evidence that high NPY levels predispose individuals to attention-deficit hyperactivity disorder (ADHD) (Oades et al., 1998; Lesch et al., 2011; Pjetri et al., 2012). While anxiety, depression, PTSD, and ADHD are associated with changes in the functional connectivity of intrinsic brain networks (Helm et al., 2018), how NPY affects functional connectivity has not been studied.

Functional connectivity techniques allow one to study correlational relationships between brain regions and reveal information about long distance communication in brain networks (van den Heuvel and Hulshoff Pol, 2010). Intrinsic brain networks may be analyzed using graph theory, which is the study of pairwise relationships or sets of relationships (Wang et al., 2010). Traditionally, functional connectivity techniques have assumed that the brain networks are static (Smith, 2012). More recent techniques have been developed to reveal dynamic connectivity, giving valuable information about how correlational relationships between brain regions change over time (Hutchison et al., 2013; Lurie et al., 2020). In this paper we investigate how NPY expression effects both the static and dynamic connectivity of the brain using graph theoretical methods and Leading Eigenvector Dynamics Analysis (LEiDA) (Cabral et al., 2017), respectively.

We focus here on connectivity within the salience network and within the default mode network, both of which have been implicated in anxiety and depression (Kaiser et al., 2015b), as well as ADHD (Norman et al., 2017). Depression has been associated with unstable connectivity (Wise et al., 2017) and hyperconnectivity (Liston et al., 2014; Kaiser et al., 2015a) in the default mode network. Anxiety has been associated with both lower (Geng et al., 2016; Xu et al., 2019) and higher intra-salience network functional connectivity (Sylvester et al., 2012), and an increase in salience network connectivity has been reported in children with greater behavioral inhibition (Taber-Thomas et al., 2016). ADHD has presented with altered connectivity to important nodes in the salience and default mode network (Fair et al., 2013) as well as increased temporal variability in default mode connectivity (Mowinckel et al., 2017).

Here we investigate differences in resting functional connectivity between two extreme groups of healthy young adults who were selected by NPY genotype (Warthen et al., 2018). Subjects in the high-NPY and low-NPY groups are genetically predisposed to express high and low levels of NPY, respectively. We use multiple techniques to evaluate functional connectivity within the salience network and default mode network. First, we used graph theory to evaluate static resting functional connectivity within the salience network and default mode network. Second, a dynamic connectivity technique [LEiDA, (Cabral et al., 2017)] was used to determine group differences in dynamic network patterns within the salience network and default mode network. Finally, we performed an exploratory seed-based analysis of connectivity strength of the nucleus accumbens (NAc) with the whole brain, based on previously discovered group differences in activation of this region during a monetary incentive delay task (Warthen et al., 2018). We hypothesized that low-NPY subjects would have a more tightly knit salience network due to their previously reported stronger NAc responses to salient stimuli (Warthen et al., 2018). Although increased NAc activity has not yet been tied to a more functionally connected salience network, early functional connectivity studies have indicated that higher activity in certain regions has been related to higher functional connectivity within their primary networks, possibly due to the plastic nature of the functional networks and their components (Greicius et al., 2003). We also hypothesized that high-NPY subjects would display faster switching rates between states, or defined connectivity patterns, within the salience network and the default mode network because of evidence that high levels of NPY predispose to hyperactivity or ADHD (Lesch et al., 2011; Warthen et al., 2018). Although intra-network variability in the salience network in hyperactivity or ADHD is not well characterized, inter-network variability has been shown to be increased in subjects with ADHD possibly indicating more volatile functional connectivity overall (Cai et al., 2018; Wang et al., 2018).

Materials and Methods

Participants and Design

Healthy adults (n = 222), aged 18–22 years, were genotyped and those with pre-specified NPY genotypes participated in an imaging visit (n = 53). This visit included questionnaires, drug and pregnancy screens, and task and resting state functional magnetic resonance imaging (fMRI). The characteristics of the subject sample and results of the task fMRI have been reported previously (Warthen et al., 2018).

Subject Selection by Genotyping

Subjects were healthy as determined by the Mini International Neuropsychiatric Interview, drug screen, and pregnancy screen. The genotyped groups were the same as we described previously (Warthen et al., 2018). In brief, 6 polymorphic markers were determined by PCR followed by Sanger sequencing (Sanger et al., 1977) to classify subjects with unambiguously low haplotypes, unambiguously high NPY haplotypes, or intermediate, rare or ambiguous combinations of haplotypes. Thirty-one low-NPY and 22 high-NPY subjects underwent MRI.

Image Acquisition and Preprocessing

Blood oxygenation-level-dependent (BOLD) responses were measured with T2 weighted imaging (TR = 2 s, TE = 28 ms, flip angle = 90°, 39 transverse slices, slice thickness = 3.5 mm, slice gap = 0 mm, FOV = 64 × 64 matrix, 3.75 × 3.75 mm in-plane resolution) on a 3.0-T Philips Ingenia scanner (Best, Netherlands) with a 15-channel head coil using single-shot echo-planar imaging. Participants were given a focus point and asked to relax with eyes open during a resting state acquisition, during which 300 volumes were collected over 10 min. Resting state acquisition for both groups took place after the monetary incentive delay task described in Warthen et al. (2018). The first 5 volumes of each session were discarded to allow for image stabilization. Slice-time correction was applied with SPM12 (v7219) using the middle slice as the reference. A rigid-body least-squares two-pass procedure was used for motion correction that first registered images to the first image of the first run, and then to the mean image that was calculated after the first pass. Images were interpolated with a 4th degree B-spline. A high-resolution T1-weighted image (turbo-field-echo, TR = 9.8 ms, TE = 4.6 ms, flip angle = 8°, 1 × 1 × 1 mm voxel size) was co-registered to the mean functional image with a rigid-body transform using a normalized mutual-information cost function. The registered T1 was segmented with unified segmentation (Ashburner and Friston, 2005) in SPM12 (v7219). The image was segmented into gray matter, white matter, and cerebrospinal fluid. Gray and white matter images were normalized to MNI space with DARTEL using the provided MNI-space template in the VBM8 toolbox. The estimated warp was then applied to the motion-corrected functional images, which were then resliced to 3 × 3 × 3 mm voxels. Smoothing was performed with an isotropic kernel (5 mm full width half maximum). Twenty-four subjects with excessive head motion (mean frame displacement > 0.20 mm) were excluded from analysis, leaving 29 subjects available for final analyses.

Connectivity Preprocessing

Resting-state time-series data were preprocessed with a custom toolbox in MATLAB (ConnTool, authored by Dr. Welsh, available upon request) (Jelsone-Swain et al., 2010; Jacobs et al., 2014). Time-series data were quadratically detrended, and regressed against six realignment parameters from motion correction along with their first derivatives and quadratic terms. Additionally, physiological artifacts were minimized through incorporation of COMPCOR by regressing out the signals from cerebral spinal fluid and white matter (Behzadi et al., 2007). Bandpass filtering was performed with a fast Fourier transform filter between 0.01 and 0.10 Hz. Preprocessing steps followed guidelines put forward by Lindquist et al. (2019). Pearson correlation coefficients were calculated between the average BOLD signal in the NAc ROI and all other voxels in the brain. To account for varying degrees-of-freedom in resting-state time-series data due to autocorrelative structures, we incorporated a recently developed technique (Afyouni et al., 2019) into ConnTool to produce the final z-score images. This technique corrects for inflated correlation coefficients due to autocorrelation as well as instantaneous and lagged cross-correlation (known as “xDF” z-scores). Each individual xDF corrected z-score image was included in a second-level SPM12 random effects analysis, to test for whole-brain group differences in connectivity strength in the NAc and dACC maps.

Seed-Based Functional Connectivity

Region of interest “seed-based” connectivity was evaluated for the nucleus accumbens (NAc) and dorsal anterior cingulate cortex (dACC) because of previously observed group differences (low-NPY vs. high-NPY) in the BOLD signal in the NAc during a monetary incentive delay task (Warthen et al., 2018) and previous associations of dACC connectivity with risk for depression (Wagner et al., 2017; Schwartz et al., 2019). The NAc region of interest was created from a standard atlas [wfu pickatlas, (Maldjian et al., 2003)] and dilated by 2 voxels to allow for shift in activation peaks (Sacchet and Knutson, 2013; Warthen et al., 2018), thresholded at 0.5 and then binarized. The dACC region of interest was created from the Desikan atlas (Desikan et al., 2006) and also thresholded at 0.5 and then binarized.

Network Definition

We defined salience network mask similarly to Vinod Menon (2015), from a publicly available functional connectivity map based on 1,000 healthy subjects (Yeo et al., 2011), with a 6 mm spherical seed at (36, 18, 4 mm) and its homolog (−36, 18, 4 mm), obtained from neurosynth.org. The salience network was thresholded to above 0.3 (Pearson’s r) and then binarized. We parcellated the resting state fMRI data from our subjects into 268 nodes as defined by Shen et al. (2013). Nodes were considered part of the salience network if they overlapped the neurosynth generated mask by at least 20% of their volume. This resulted in the inclusion of 20 nodes (Table 1). The default mode network was defined similarly, from a connectivity map from neurosynth.org with a seed at (4, −54, 28 mm) in the posterior cingulate cortex, as by Uddin et al. (2009). Nodes that overlapped by 20% were included, resulting in a 26-node network (Table 1). The objective way that we identified the networks did result in some asymmetry between hemispheres, we have chosen to leave this asymmetry in our network analyses so as not to bias the findings and remain true to the 20% overlap rule. Figure 1 shows the anatomical locations of these regions in Montreal Neurological Institute space. Here we focus only on the salience network and default mode network rather than all 268 nodes in the atlas because of a previously demonstrated effect of NPY on nodes of the salience network (Brown et al., 2000; Adewale et al., 2007; Quarta et al., 2011) in animal studies, as well as the relationship of network patterns in the default mode network and depression and anxiety (Sylvester et al., 2012; Qiao et al., 2017).

TABLE 1
www.frontiersin.org

Table 1. Shen atlas regions of interest.

FIGURE 1
www.frontiersin.org

Figure 1. The default mode network (green, top) and salience network (orange, bottom). Numbers indicate regions of the Shen atlas.

Static Connectivity Analysis

We used the BRAPH [BRain Analysis using graPH theory (1.0.0) (Mijalkov et al., 2017)] package for Matlab (R2015b) to calculate measures of network integration (characteristic path length, global efficiency), network segregation (clustering coefficient, modularity), and small-worldness (Rubinov and Sporns, 2010; Sporns, 2011). These measures allow us to compare the high and low-NPY groups in terms of global network cohesion and local communication of the network, as both may be important in functional differences. Characteristic path length is the average distance between a node and all other nodes in terms of edge length traversed. Efficiency is related to characteristic path length, but is the inverse of average shortest path length. The clustering coefficient is the fraction closed triangles to open triangles around an individual node; where triangles occur when three nodes are all connected to each other by edges. Clustering coefficient therefore gives an idea of how much grouping, or local coherent oscillation, is occurring around a given node. Modularity is used to measure how a network is divided into subunits, or communities, by comparing the number of edges inside a cluster with the number of edges that would be found in that cluster in a random graph. Small-worldness measures the balance between local grouping and global cooperation in a network. It consists of the ratio of the normalized clustering coefficient over the normalized characteristic path length. These measures (portrayed in Supplementary Figure 1) characterize how closely certain regions of interest are connected. We also investigated node-specific measures of strength (sum of all edge weights feeding into a node) as well as eccentricity (maximum length from one node to its farthest connection) to estimate how our nodes (brain regions) are connected in a local and global network sense. We used a weighted undirected graph with Pearson correlation, and a correlation threshold of zero.

Dynamic Connectivity Analysis

We were primarily interested in within-network state representation and switching activity within the salience and default mode networks, so we employed the Leading Eigenvector Dynamics Analysis (LEiDA) method1 described by Cabral et al. (2017) and Deco et al. (2019). LEiDA is a data-informed way of defining brain states and switching rates between brain states based on identification of large-scale functional connectivity patterns. A leading eigenvector is a vector containing information describing the largest variations in a matrix and is used to find the dominant connectivity pattern at a certain timepoint in vector form. The dominant connectivity pattern can then be reconstructed with the outer product of the vector, where each element of the vector is multiplied by every other element of the vector forming a matrix.

Dynamic connectivity patterns were determined as by Cabral et al. (2017). Leading eigenvectors were calculated at each fMRI time point for each person from the BOLD Phase Coherence Connectivity (Glerean et al., 2012), and the leading eigenvectors combined for all subjects were used for state calculation in each network. Repeating functional connectivity states were determined with k-means clustering (Lloyd, 1982), with a Dunn’s index (Dunn, 1974) determining the optimal number of states.

Statistical Analyses

Statistical comparisons of network measures were performed in R (version 3.6.1). Static and dynamic network measures were calculated per individual and group differences were determined with a standard two-sided t-test. Static and dynamic measures for both networks were normally distributed as determined by a Shapiro-Wilk test. Effect sizes were calculated with Cohen’s d (“cohen.d” function, “effsize” package, version 0.7.6). For tests of the primary hypothesis of NPY group differences in connectivity, Holm’s method was used for multiple-comparison correction (“p.adjust” function, Base R).

Several additional exploratory analyses were performed. Associations between static and dynamic network measures in the salience network as well as associations between traits (Table 2) and network measures were tested with a linear model (“lm” function, “lme4” package, version 1.1.21 in R) while controlling for the effects of NPY group. Sex differences were tested with a standard two-sample t-test. Salience condition BOLD contrasts (high-NPY vs. low-NPY) from a previously published monetary incentive delay task (Warthen et al., 2018) for the nucleus accumbens, dorsal anterior insula, and dorsal anterior cingulate cortex were also tested for associations with static and dynamic network measures with a linear model while controlling for NPY group because of the prominence of the nucleus accumbens, anterior insula, and dorsal anterior cingulate cortex in the salience network (Menon, 2015).

TABLE 2
www.frontiersin.org

Table 2. Demographic, physiological, and clinical characteristics of the sample.

Results

Subjects

Twenty-nine subjects were genotyped as high-NPY or low-NPY and provided resting functional MRI data that survived quality-control procedures. The final groups of included subjects did not vary in average frame displacement (p > 0.05). No group differences were found for any demographics, physiological, trait, or state measures (Table 2).

Static Connectivity: Salience Network

Low-NPY subjects showed higher values of salience network integration, through shorter characteristic path length (p = 0.0012 (uncorrected), d = 1.27) and higher global network efficiency [p = 0.0065 (uncorrected), d = 1.063; Figure 2]. High and low-NPY groups showed split results in measures of network segregation, with high-NPY subjects displaying higher modularity [p = 0.041 (uncorrected), d = 0.76], but low-NPY subjects having higher clustering [p = 0.0025 (uncorrected), d = 1.16]. Low-NPY subjects were found to have higher measures of small worldness [p = 0.0030 (uncorrected), d = 1.17] (Figure 2). These effects all survived Holm’s correction for multiple comparison (characteristic path length, p = 0.006; global efficiency, p = 0.013; modularity, p = 0.041; clustering, p = 0.0010; small worldness, p = 0.010).

FIGURE 2
www.frontiersin.org

Figure 2. Graph theory parameters for high and low NPY groups in the salience and default mode network. Error bars represent standard error of the mean. p < 0.05, ∗∗p < 0.001 as calculated by a two-sided t-test.

To further evaluate the underlying causes of these group differences, we tested graph-theoretic metrics at the level of individual nodes (eccentricity and strength). Although not surviving Holm’s correction for multiple comparison for twenty regions, higher node strength was found in the low-NPY group in the bilateral dACC (region 15: p = 0.24, 28: p = 0.094, 150: p = 0.088, 221: p = 0.36 after Holm’s correction for multiple comparisons), the right dorsolateral prefrontal cortex (region 9: p = 0.36, 11: p = 0.026, 21: p = 0.16), the right dorsal anterior insula (region 20: p = 0.24), the left anterior insula (region 169: p = 0.45), and the left lateral sulci (region 163: p = 0.36, 181: p = 0.45) (Figure 3). Higher eccentricity was found in the high-NPY group in the bilateral dACC (region 15: p = 0.011, 28: p = 0.041, 150: p = 0.024, 221: p = 0.17 after Holm’s correction for multiple comparisons), the right dorsolateral prefrontal cortex (region 9: p = 0.11, 11: p = 0.16, 19: p = 0.25), the right anterior (region 34: p = 0.17) and dorsal anterior insula (region 35: p = 0.25), the left anterior insula (region 169: p = 0.11), the left lateral sulci (region 46: p = 0.17, 163: p = 0.10, 181: p = 0.25) (not significant after correction) (Figure 3). Additionally average local efficiency across all nodes within the salience network was found to be significantly different (p < 2 × 10 –16). Further details can be found in Supplementary Material.

FIGURE 3
www.frontiersin.org

Figure 3. (A) Salience network node strength, where Low NPY > High NPY. (B) Salience network node eccentricity, where High NPY > Low nPY. *p < 0.05, **p < 0.001 as calculated by a two-sided t-test. Regions: 9, Right DLPFC; 11, Right DLPFC; 15, Right dACC; 34, Right Al; 35, Right Dorsal Al; 36, Right VLPFC; 45, Inf. Parietal Lobule; 46, Left lateral sulcus; 150, Left dACC; 155, Left Dosal Al; 163, Left lateral sulcus; 168, Left Al; 181, Left Lateral sulcus; 221, Left dACC.

Static Connectivity: Default Mode Network

No significant differences were found between high-NPY and low-NPY groups in network integration [global efficiency p = 0.43 (uncorrected), d = 0.32; characteristic path length p = 0.42 (uncorrected), d = 0.32], network segregation [clustering p = 0.41 (uncorrected), d = 0.32; modularity p = 0.28 (uncorrected), d = 0.42], or small-worldness (small worldness p = 0.21 (uncorrected), d = 0.45) in the default mode network (Figure 2). However, average local efficiency across all nodes within the default mode network was found to be significantly different (p = 0.026).

Dynamic Connectivity: Salience Network

The salience network was best described by two states according to k-clustering and Dunn’s index. Figure 4 shows each state along with its connectivity matrix representation. Positive and negative eigenvector values represent functional grouping of nodes within the state.

FIGURE 4
www.frontiersin.org

Figure 4. Salience network state 1 (A), and state 2 (B), where anatomical regions are colored by eigen vector value. Matrices show the outer product of the state defining eigen vector. (C) Fractional occupancy of high and low NPY groups in state 1 and state 2 of the salience network. (D) Switching probability between network states. Error bars represent standard error of the mean. *p < 0.05 as calculated by a two-sided t-test.

In state 1, the regions centered on dACC were more out of phase with the rest of the nodes in the salience network. In State 2, all nodes in the network were coherent.

High-NPY subjects spent more time in state 1 than low-NPY subjects, conversely low-NPY subjects spent more time in state 2 than high-NPY subjects [p = 0.038 (uncorrected), d = 0.89]. High-NPY subjects also displayed a higher switching probability between the two states [p = 0.047 (uncorrected), d = 0.75].

Dynamic Connectivity: Default Mode Network

The default mode network was best described by 3 states, shown in Figure 5. In DMN state 1, all the regions of interest were closely correlated. In state 2, we see two distinct groupings. The first set of regions that includes the right anterior frontal cortex, the left dorsal anterior frontal cortex, left anterior and posterior temporal lobe as well as the right posterior temporal lobe, the right posterior cingulate cortex, and the left temporal parietal junction. The second set of regions includes the left anterior frontal cortex, the right dorsal anterior frontal cortex, the right anterior temporal lobe, the left posterior cingulate cortex and the left and right inferior parietal lobes. In the third state, part of the posterior cingulate cortex (−7, −42, 13) (region 227) is anti-correlated with the rest of the network.

FIGURE 5
www.frontiersin.org

Figure 5. Default mode network state 1 (A), state 2 (B), and state 3 (C), where anatomical regions are colored by eigen vector value. Matrices show the outer product of the state defining eigen vector. (D) Fractional occupancy of high and low NPY groups in state 1, state 2, and state 3 of the default mode network. (E) Switching probability between network states. Error bars represent standard error of the mean. No significant effects as measured by a two-sided t-test.

No significant differences between the high-NPY and low-NPY groups were found for fractional occupancy in any of the states, or the average switching probability (p > 0.05, two-sided t-test). Both groups spent a majority of time in state 1. Switching probabilities between each state in the DMN were tested for group differences; none were found that survived correction for multiple comparison. Details can be found in Supplementary Material.

Static Connectivity Association With Dynamic Connectivity in the Salience Network

Salience network switching probability was negatively associated with global efficiency [p = 6.1 × 10–5 (uncorrected), Pearson’s r = −0.75] and clustering [p = 5.4 × 10–4 (uncorrected), Pearson’s r = 0.68], but positively associated with characteristic path length [p = 2.1 × 10–4 (uncorrected), Pearson’s r = 0.72] while controlling for NPY group. Fractional occupancy was also negatively associated with global efficiency [p = 2.6 × 10–4 (uncorrected), Pearson’s r = 0.70] and clustering [p = 0.0017 (uncorrected), Pearson’s r = 0.63], but positively associated with characteristic path length [p = 9.3 × 10–4 (uncorrected), Pearson’s r = 0.66] and modularity [p = 0.030 (uncorrected), Pearson’s r = 0.37] while controlling for NPY group.

Functional Connectivity of the Nucleus Accumbens and Dorsal Anterior Cingulate Cortex

We compared seed-based functional connectivity from the NAc and dACC between groups in an exploratory analysis (Supplementary Figure 5 and Supplementary Table 4). No group difference was found between group connectivity maps after whole brain correction (p > 0.05, FWE).

Trait, Region of Interest Blood Oxygenation-Level-Dependent Contrast, and Sex Effects on Connectivity

Psychological traits were not found to be associated with connectivity while controlling for NPY group after correction for comparisons of 12 trait measures (p > 0.05, linear model). Exploratory correlations between network measures and trait and state questionnaires are shown in Supplementary Figure 6. No differences were found between men and women in any of the static or dynamic connectivity measures according to standard two-sided t-tests (p > 0.05, two-sided t-test). The MID task salience BOLD contrast (high vs. low) for the nucleus accumbens, dorsal anterior insula, and dorsal anterior cingulate cortex were not correlated with any of the static or dynamic connectivity measures (p > 0.05, linear model).

Discussion

NPY genotype has an impact on both static and dynamic functional connectivity in the salience network. Although no significant differences were found in functional connectivity maps from seeds in the NAc and dACC, in static network analyses of the salience network low-NPY subjects were found to have shorter path lengths, higher global efficiency, higher clustering, higher small-worldness, and higher node strength on average within the salience network. High-NPY subjects showed higher modularity and node eccentricity in the salience network. Salience network differences did not appear to be driven by any specific node, but rather by many nodes spread throughout the network. No differences were found between the groups in the default mode network. Dynamically, low-NPY subjects spent more time in a state displaying more coordination while high-NPY subjects spent more time in a state that had less coordinated connectivity in the salience network. High-NPY subjects also switched between states more often in this network. Again, no differences were found between groups in the default mode network.

We have previously shown that high and low-NPY subjects from the same sample differ in NAc activation during a monetary incentive delay task (Warthen et al., 2018). However, we did not find that functional connectivity maps based on the same NAc seed were statistically different between groups. These results indicate that, while NAc activation during reward behavior is different between groups, the two groups do not necessarily differ with respect to the information transfer to and from this region.

Because the salience network (Menon, 2015) and default mode network are implicated in many psychiatric disorders, we evaluated the impact of NPY expression on those networks. We found a hyperconnected salience network in low-NPY subjects. Similar findings have been observed in PTSD (Akiki et al., 2018; Fani et al., 2019) during rest as well as during eye-contact (Thome et al., 2014). In social and general anxiety, regions within the salience network (including the insula, anterior cingulate, and prefrontal cortex) have been observed to be hyperactive during rest and task functional connectivity (Etkin and Wager, 2007; Brühl et al., 2014). Structural hyperconnectivity between the anterior insula and amygdala, as measured with diffusion imaging, has also been shown to be associated with state and trait anxiety (Baur et al., 2013). Other studies on anxiety have reported heightened activity in the salience network (Sylvester et al., 2012), and point to increased error notification originating in the dACC (Hajcak et al., 2003; Paulus et al., 2004). Discordant results have also been observed, e.g., decreased connectivity in the salience network in anxiety (Geng et al., 2016; Xu et al., 2019). Increased coherence within the salience network may reflect increased vigilance in stimulus notification, which could lead to anxiety. Future work should examine the role of NPY in the activity of the salience network and how such network alterations impact risk for anxiety and depression.

Although low-NPY has been associated with depression, the hyperconnectivity within the salience network observed here has not been found in depression. A previous study showed a decrease in connectivity to the nucleus accumbens in subjects with depression compared to controls (Helm et al., 2018). Although we did not include this node in our salience network, we don’t see lower connectivity in low-NPY subjects in any node in the salience network when compared to high-NPY subjects. In ADHD, von Rhein et al. (2019) reported decreased functional connectivity between the salience network and executive control network, which may also point back to less coordinated functional connectivity within the salience network as we see in high-NPY subjects. Higher connectivity within the salience network may indicate a lower threshold for stimulus notification, which would fit symptoms seen in anxiety and PTSD. ADHD has been associated with lower functional connectivity within the salience network (von Rhein et al., 2019), as well as a possible higher necessary stimulus threshold for sustained attention (Tegelbeckers et al., 2015). The high-NPY group may share this parallel with subjects with ADHD.

The high-NPY group displayed a higher switching probability between states in the salience network. Subjects with ADHD have previously shown increased switching between networks (Cai et al., 2018; Scofield et al., 2019). Although here we examine within network switching probabilty, between and within network switching probabilities may be related. Depression has been shown to present with impaired salience network mediated switching into the central executive network (Wang et al., 2016), which may relate to the lower switching likelihood within the salience network shown here by the low-NPY group. Statically, we see lower long-distance global efficiency within the salience network in the high-NPY subjects, along with higher modularity. These findings are similar to network analyses in children with ADHD (Wang et al., 2009; Cai et al., 2018), which are of interest given that greater NPY function may be related to hyperactivity (Lesch et al., 2011), impulsivity (Bari et al., 2015), and possibly ADHD (Scassellati et al., 2012; Kourtesis et al., 2015).

The functional connectivity differences observed here between low and high-NPY groups seem to be specific to the salience network, as we did not find group differences in the default mode network. It would be reasonable to expect group differences in default mode network connectivity, as low-NPY is associated with the development of depression (Mickey et al., 2011), and stronger default mode network connectivity has been reported in subjects with depression (Wang et al., 2016; Helm et al., 2018), and low-NPY status has been associated with greater activation of medial prefrontal cortex (Mickey et al., 2011). The functional connectivity of those at risk for the development of depression may not look like the functional connectivity of those with active depression. PTSD, also associated with low-NPY (Yehuda et al., 2006), has been characterized by a generally hypoactive default mode network (Akiki et al., 2018), which we do not see in low-NPY subjects. This may indicate that a state change is necessary from the at-risk state to the disease state (e.g., depressed or PTSD) before differences in default mode network connectivity are observed. Overall these findings suggest that NPY exerts effects on risk for psychiatric disorders primarily through the salience network, and not the default mode network.

We found a negative relationship between salience network switching probability and global efficiency and clustering, but a positive relationship with characteristic path length. This is not surprising because a higher network switching rate should result in lower correlations between regions in the network over time, as they spend less time in coherent oscillations. A lower network switching rate would allow for stronger long-distance correlations, resulting in a higher characteristic path length in a weighted graph. Fractional occupancy showed these same relationships, a negative correlation with global efficiency and clustering but a positive association with characteristic pathlength, possibly for the same reason. Fractional occupancy also showed a positive correlation with modularity. A longer time spent in a certain network state, or higher fractional occupancy, may allow for stronger local correlations, resulting in higher modularity.

Our study has several limitations. Although the imaged sub-sample was genetically selected from a relatively large sample (>200), the number of imaged subjects was reduced, and many had to be excluded to avoid artifacts related to head motion. The resulting sample size limited the power to detect moderate-sized between-group effects. A replication sample would boost the strength of these findings, and this should be investigated further in the future but the practical limitations inhibit our ability to recruit another 200 subjects to sample from at this time. Furthermore, to address the multiple-comparison problem, we had to focus our hypotheses on two brain networks of interest. We recognize that there is no hard definition of what is and is not the salience or default mode network. These findings should be validated in alternate definitions of these networks in future studies. The distribution of NPY haplotypes in the general population also led to an imbalance of NPY group numbers which could impact the k-clustering algorithm. As a reviewer stated we do not know the circulating, CSF, or cellular levels of NPY in these subjects in the genetically predicted groups. As we are studying healthy humans, CSF and cellular levels of NPY are not feasible to obtain. Additionally, although we do not have circulating or serum levels of NPY these measures likely would not represent CSF levels of NPY (Dötsch et al., 1997). We suggest animal studies to determine the relationship of these measures of the brain regions of interest. Here we only examined resting-state functional connectivity; task-based functional connectivity analysis of the salience network may also provide additional insight into the effects of NPY on intrinsic functional networks. Finally, the resting state scans in this study were collected after a monetary incentive delay task for both groups. Previous studies have shown that prior tasks may impact resting state (Waites et al., 2005), however, the effect would be similar for both groups.

Generally, we found increased coordination across the salience network among low-NPY subjects, and a salience network that was less efficient and coherent among high-NPY subjects. Differences in connectivity in this network may point to differences in emotional regulation and salience signaling. These differences in processing may provide high-NPY subjects with resilience against depression, anxiety, and PTSD, but which also may put them at higher risk for attentional or hyperactive disorders. Stronger, more coordinated connectivity in the salience network may indicate a lower threshold for stimulus detection. In disorders of anxiety and PTSD, this may result in hyper-vigilance and constant notification of realistically irrelevant stimuli. Lower connectivity, as generally seen in cases of ADHD, may represent a higher necessary stimulus threshold for coherent activity of the salience network, or paying attention to the given stimulus. This lower connectivity may also result in higher switching probabilities within the network, generally representing less coordinated activity in the network over time. Anxiety, PTSD, and ADHD have all been suggested to stem from a lack of regulation in brain networks (Etkin and Wager, 2007; Hoekzema et al., 2014). Our results suggest that NPY may play a role in network regulation, and this should be a focus of future studies on neurobiological risk factors for anxiety and mood disorders.

Conclusion

We have shown that low levels of NPY are associated with a more closely knit salience network, as determined by graph theory measures. Additionally, subjects with high levels of NPY displayed higher switching probabilities within the salience network, along with more time spent in a less strongly connected state. In the default mode network, no such differences were observed. NPY may exert effects on the risk of development of psychiatric disorders through subtly varied activity in the salience network.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the University of Michigan IRB. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

KW wrote the article, designed the analysis plan, and analyzed the data. RW developed some of the analytical tools, contributed to the analysis and analysis plan, and edited the article. BS collected the data. VK contributed to the analysis and analysis plan. MB contributed to the analysis plan and the initial data collection. BM edited the article, designed the initial data collection, collected the initial data, and contributed to the analysis and analysis plan. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by NIMH (K23 MH092648), NCRR (UL1 RR024986), and NCATS (2UL1 TR000433). ConnTool was developed in part with support from NINDS (R01 NS052514). We report no biomedical financial interests related to this work in the past or in the foreseeable future.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsys.2021.629488/full#supplementary-material

Footnotes

  1. ^ https://github.com/juanitacabral/LEiDA

References

Adewale, A. S., Macarthur, H., and Westfall, T. C. (2007). Neuropeptide Y-induced enhancement of the evoked release of newly synthesized dopamine in rat striatum: mediation by Y2 receptors. Neuropharmacology 52, 1396–1402. doi: 10.1016/j.neuropharm.2007.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Adrian, T. E. E., Allen, J. M. M., Bloom, S. R. R., Ghatei, M. A. A., Rossor, M. N. N., Roberts, G. W. W., et al. (1983). Neuropeptide Y distribution in human brain. Nature 306, 584–586. doi: 10.1038/306584a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Afyouni, S., Smith, S. M., and Nichols, T. E. (2019). Effective degrees of freedom of the Pearson’s correlation coefficient under autocorrelation. Neuroimage 199, 609–625. doi: 10.1016/j.neuroimage.2019.05.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Akiki, T. J., Averill, C. L., and Abdallah, C. G. (2018). A Network-Based Neurobiological Model of PTSD: Evidence from Structural and Functional Neuroimaging Studies. Curr. Psychiatry Rep. 19, 840–844. doi: 10.1007/s11920-017-0840-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen, Y. S., Adrian, T. E., Allen, J. M., Tatemoto, K., Crow, T. J., Bloom, S. R., et al. (1983). Neuropeptide Y distribution in the rat brain. Science 221, 877–879. doi: 10.1126/SCIENCE.6136091

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, J., and Friston, K. J. (2005). Unified segmentation. Neuroimage 26, 839–851. doi: 10.1016/j.neuroimage.2005.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Bari, A., Dec, A., Lee, A. W., Lee, J., Song, D., Dale, E., et al. (2015). Enhanced inhibitory control by neuropeptide Y Y5 receptor blockade in rats. Psychopharmacology 232, 959–973. doi: 10.1007/s00213-014-3730-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Baur, V., Hänggi, J., Langer, N., and Jäncke, L. (2013). Resting-state functional and structural connectivity within an insula-amygdala route specifically index state and trait anxiety. Biol. Psychiatry 73, 85–92. doi: 10.1016/j.biopsych.2012.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Behzadi, Y., Restom, K., Liau, J., and Liu, T. T. (2007). A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. Neuroimage 37, 90–101. doi: 10.1016/j.neuroimage.2007.04.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, C. M., Coscina, D. V., and Fletcher, P. J. (2000). The rewarding properties of neuropeptide Y in perifornical hypothalamus vs. nucleus accumbens. Peptides 21, 1279–1287. doi: 10.1016/S0196-9781(00)00270-9

CrossRef Full Text | Google Scholar

Brühl, A. B., Delsignore, A., Komossa, K., and Weidt, S. (2014). Neuroimaging in social anxiety disorder-A meta-analytic review resulting in a new neurofunctional model. Neurosci. Biobehav. Rev. 2014:003. doi: 10.1016/j.neubiorev.2014.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabral, J., Vidaurre, D., Marques, P., Magalhães, R., Moreira, P. S., Soares, J. M., et al. (2017). Cognitive performance in healthy older adults relates to spontaneous switching between states of functional connectivity during rest. Sci. Rep. 7, 5425–5427. doi: 10.1038/s41598-017-05425-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, W., Chen, T., Szegletes, L., Supekar, K., and Menon, V. (2018). Aberrant Time-Varying Cross-Network Interactions in Children With Attention-Deficit/Hyperactivity Disorder and the Relation to Attention Deficits. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 3, 263–273. doi: 10.1016/j.bpsc.2017.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Deco, G., Cruzat, J., Cabral, J., Tagliazucchi, E., Laufs, H., Logothetis, N. K., et al. (2019). Awakening: Predicting external stimulation to force transitions between different brain states. Proc. Natl. Acad. Sci. 2019:201905534. doi: 10.1073/pnas.1905534116

PubMed Abstract | CrossRef Full Text | Google Scholar

Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Dötsch, J., Adelmann, M., Englaro, P., Dötsch, A., Hänze, J., Blum, W., et al. (1997). Relation of leptin and neuropeptide Y in human blood and cerebrospinal fluid. J. Neurol. Sci. 151, 185–188. doi: 10.1016/S0022-510X(97)00116-0

CrossRef Full Text | Google Scholar

Dunn, J. C. (1974). A Fuzzy Relative of the ISODATA Process and Its Use in Detecting Compact Well-Separated Clusters A Fuzzy Relative of the ISODATA Process and Its Use in Detecting Compact Well-Separated Clusters. J. Cybern. 3, 32–57. doi: 10.1080/01969727308546046

CrossRef Full Text | Google Scholar

Etkin, A., and Wager, T. D. (2007). Functional neuroimaging of anxiety: A meta-ana lysis of emotional processing in PTSD, social anxiety disorder, and specific phobia. Am. J. Psychiatry 2007:07030504. doi: 10.1176/appi.ajp.2007.07030504

PubMed Abstract | CrossRef Full Text | Google Scholar

Fair, D. A., Nigg, J. T., Iyer, S., Bathula, D., Mills, K. L., Dosenbach, N. U. F., et al. (2013). Distinct neural signatures detected for ADHD subtypes after controlling for micro-movements in resting state functional connectivity MRI data. Front. Syst. Neurosci. 6:00080. doi: 10.3389/fnsys.2012.00080

PubMed Abstract | CrossRef Full Text | Google Scholar

Fani, N., King, T. Z., Clendinen, C., Hardy, R. A., Surapaneni, S., Blair, J. R., et al. (2019). Attentional control abnormalities in posttraumatic stress disorder: Functional, behavioral, and structural correlates. J. Affect. Disord. 253, 343–351. doi: 10.1016/j.jad.2019.04.098

PubMed Abstract | CrossRef Full Text | Google Scholar

Geng, H., Li, Xuebing, Chen, J., Li, Xinying, et al. (2016). Decreased Intra- and Inter-Salience Network Functional Connectivity is Related to Trait Anxiety in Adolescents. Front. Behav. Neurosci. 9:00350. doi: 10.3389/fnbeh.2015.00350

PubMed Abstract | CrossRef Full Text | Google Scholar

Glerean, E., Salmi, J., Lahnakoski, J. M., Jääskeläinen, I. P., and Sams, M. (2012). Functional Magnetic Resonance Imaging Phase Synchronization as a Measure of Dynamic Functional Connectivity. Brain Connect. 2, 91–101. doi: 10.1089/brain.2011.0068

PubMed Abstract | CrossRef Full Text | Google Scholar

Greicius, M. D., Krasnow, B., Reiss, A. L., Menon, V., and Raichle, M. E. (2003). Functional connectivity in the resting brain: A network analysis of the default mode hypothesis. PNAS 100, 253–258.

Google Scholar

Hajcak, G., Mcdonald, N., and Simons, R. F. (2003). Anxiety and error-related brain activity. Biol. Psychol. 64, 77–90. doi: 10.1016/S0301-0511(03)00103-0

CrossRef Full Text | Google Scholar

Heilig, M., Zachrisson, O., Thorsell, A., Ehnvall, A., Mottagui-Tabar, S., Sjögren, M., et al. (2004). Decreased cerebrospinal fluid neuropeptide Y (NPY) in patients with treatment refractory unipolar major depression: preliminary evidence for association with preproNPY gene polymorphism. J. Psychiatr. Res. 38, 113–121. doi: 10.1016/S0022-3956(03)00101-8

CrossRef Full Text | Google Scholar

Helm, K., Viol, K., Weiger, T. M., Tass, P. A., Grefkes, C., Del Monte, D., et al. (2018). Neuronal connectivity in major depressive disorder: a systematic review. Neuropsychiatr. Dis. Treat. 14, 2715–2737. doi: 10.2147/NDT.S170989

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoekzema, E., Carmona, S., Ramos-Quiroga, J. A., Richarte Fernández, V., Bosch, R., Soliva, J. C., et al. (2014). An independent components and functional connectivity analysis of resting state FMRI data points to neural network dysregulation in adult ADHD. Hum. Brain Mapp. 35, 1261–1272. doi: 10.1002/hbm.22250

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., et al. (2013). Dynamic functional connectivity: Promises, issues, and interpretations. Neuroimage 80, 360–378. doi: 10.1016/j.neuroimage.2013.05.079

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobs, R. H., Jenkins, L. M., Gabriel, L. B., Barba, A., and Ryan, K. A. (2014). Increased Coupling of Intrinsic Networks in Remitted Depressed Youth Predicts Rumination and Cognitive Control. PLoS One 9:104366. doi: 10.1371/journal.pone.0104366

PubMed Abstract | CrossRef Full Text | Google Scholar

Jelsone-Swain, L. M., Fling, B. W., Seidler, R. D., Hovatter, R., Gruis, K., and Welsh, R. C. (2010). Reduced interhemispheric funcitonal connectivity in the motor cortex during rest in limb-onset amyotrophic lateral sclerosis. Front. Syst. Neurosci. 4:00158. doi: 10.3389/fnsys.2010.00158

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaiser, R. H., Whitfield-Gabrieli, S., Dillon, D. G., Goer, F., Beltzer, M., Minkel, J., et al. (2015b). Dynamic Resting-State Functional Connectivity in Major Depression. Neuropsychopharmacology 41, 1822–1830. doi: 10.1038/npp.2015.352

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaiser, R. H., Andrews-Hanna, J. R., Wager, T. D., and Pizzagalli, D. A. (2015a). Large-scale network dysfunction in Major Depressive Disorder: Meta-analysis of resting-state functional connectivity HHS Public Access. JAMA Psychiatry 72, 603–611. doi: 10.1001/jamapsychiatry.2015.0071

PubMed Abstract | CrossRef Full Text | Google Scholar

Kourtesis, I., Kasparov, S., Verkade, P., and Teschemacher, A. G. (2015). Ultrastructural correlates of enhanced norepinephrine and neuropeptide y cotransmission in the spontaneously hypertensive rat brain. ASN Neuro 7:1759091415610115. doi: 10.1177/1759091415610115

PubMed Abstract | CrossRef Full Text | Google Scholar

Lesch, K.-P., Selch, S., Renner, T. J., Jacob, C., Nguyen, T. T., Hahn, T., et al. (2011). Genome-wide copy number variation analysis in attention-deficit/hyperactivity disorder: association with neuropeptide Y gene dosage in an extended pedigree. Mol. Psychiatry 16, 491–503. doi: 10.1038/mp.2010.29

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindquist, M. A., Geuter, S., Wager, T. D., and Caffo, B. S. (2019). Modular preprocessing pipelines can reintroduce artifacts into fMRI data. Hum. Brain Mapp. 40, 2358–2376. doi: 10.1002/hbm.24528

PubMed Abstract | CrossRef Full Text | Google Scholar

Liston, C., Chen, A. C., Zebley, B. D., Drysdale, A. T., Gordon, R., Leuchter, B., et al. (2014). Default mode network mechanisms of transcranial magnetic stimulation in depression. Biol. Psychiatry 76, 517–526. doi: 10.1016/j.biopsych.2014.01.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Lloyd, S. P. (1982). Least squares quantization in {PCM}. Special issue on quantization. IEEE Trans. Inf. Theory 28, 129–137.

Google Scholar

Lurie, D. J., Kessler, D., Bassett, D. S., Betzel, R. F., Breakspear, M., Kheilholz, S., et al. (2020). Questions and controversies in the study of time-varying functional connectivity in resting fMRI. Netw. Neurosci. 4, 30–69. doi: 10.1162/netn_a_00116

CrossRef Full Text | Google Scholar

Maldjian, J. A., Laurienti, P. J., Kraft, R. A., and Burdette, J. H. (2003). An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. Neuroimage 19, 1233–1239.

Google Scholar

Menon, V. (2015). Salience Network. Brain Mapp. Encycl. Ref. 2, 597–611. doi: 10.1016/B978-0-12-397025-1.00052-X

CrossRef Full Text | Google Scholar

Mickey, B. J., Zhou, Z., Heitzeg, M. M., Heinz, E., Hodgkinson, C. A., Hsu, D. T., et al. (2011). Emotion Processing, Major Depression, and Functional Genetic Variation of Neuropeptide Y. Arch. Gen. Psychiatry 68:158. doi: 10.1001/archgenpsychiatry.2010.197

PubMed Abstract | CrossRef Full Text | Google Scholar

Mijalkov, M., Kakaei, E., Pereira, J. B., Westman, E., and Volpe, G. (2017). BRAPH: A graph theory software for the analysis of brain connectivity for the Alzheimer’s Disease Neuroimaging Initiative. PLoS One 12:1–23. doi: 10.1371/journal.pone.0178798

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan, C. A., Rasmusson, A. M., Wang, S., Hoyt, G., Hauger, R. L., and Hazlett, G. (2002). Neuropeptide-Y, cortisol, and subjective distress in humans exposed to acute stress: replication and extension of previous report. Biol. Psychiatry 52, 136–142. doi: 10.1016/S0006-3223(02)01319-7

CrossRef Full Text | Google Scholar

Mowinckel, A. M., Alnaes, D., Pedersen, M. L., Ziegler, S., Fredriksen, M., Kaufmann, T., et al. (2017). Increased default-mode variability is related to reduced task-performance and is evident in adults with ADHD. Neuroimage Clin. 16, 369–382. doi: 10.1016/j.nicl.2017.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Norman, L. J., Carlisi, C. O., Christakou, A., Cubillo, A., Murphy, C. M., Chantiluke, K., et al. (2017). Shared and disorder-specific task-positive and default mode network dysfunctions during sustained attention in paediatric Attention-Deficit/Hyperactivity Disorder and obsessive/compulsive disorder. NeuroImage Clin. 15, 181–193. doi: 10.1016/j.nicl.2017.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Oades, R. D., Daniels, R., and Rascher, W. (1998). Plasma neuropeptide-Y levels, monoamine metabolism, electrolyte excretion and drinking behavior in children with attention-deficit hyperactivity disorder. Psychiatry Res. 80, 177–186. doi: 10.1016/S0165-1781(98)00064-X

CrossRef Full Text | Google Scholar

Paulus, M. P., Feinstein, J. S., Simmons, A., and Stein, M. B. (2004). Anterior cingulate activation in high trait anxious subjects is related to altered error processing during decision making. Biol. Psychiatry 55, 1179–1187. doi: 10.1016/j.biopsych.2004.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Pjetri, E., Adan, R. A., Herzog, H., de Haas, R., Oppelaar, H., Spierenburg, H. A., et al. (2012). NPY receptor subtype specification for behavioral adaptive strategies during limited food access. Genes Brain Behav. 11, 105–112. doi: 10.1111/j.1601-183X.2011.00732.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiao, J., Li, A., Cao, C., Wang, Z., Sun, J., and Xu, G. (2017). Aberrant Functional Network Connectivity as a Biomarker of Generalized Anxiety Disorder. Front. Hum. Neurosci. 11:00626. doi: 10.3389/fnhum.2017.00626

PubMed Abstract | CrossRef Full Text | Google Scholar

Quarta, D., Leslie, C. P. P., Carletti, R., Valerio, E., and Caberlotto, L. (2011). Central administration of NPY or an NPY-Y5 selective agonist increase in vivo extracellular monoamine levels in mesocorticolimbic projecting areas. Neuropharmacology 60, 328–335. doi: 10.1016/j.neuropharm.2010.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Raichle, M. E. (2015). The Brain’s Default Mode Network. Annu. Rev. Neurosci. 38, 433–447. doi: 10.1146/annurev-neuro-071013-014030

PubMed Abstract | CrossRef Full Text | Google Scholar

Reichmann, F., and Holzer, P. (2016). Neuropeptide Y: A stressful review. Neuropeptides 55, 99–109. doi: 10.1016/j.npep.2015.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: Uses and interpretations. Neuroimage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sacchet, M. D., and Knutson, B. (2013). Spatial smoothing systematically biases the localization of reward-related brain activity. Neuroimage 66, 270–277. doi: 10.1016/j.neuroimage.2012.10.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanger, F., Nicklen, S., and Coulson, A. R. (1977). DNA sequencing with chain-terminating inhibitors. Proc. Natl. Acad. Sci. U S A. 74, 5463–5467. doi: 10.1073/pnas.74.12.5463

PubMed Abstract | CrossRef Full Text | Google Scholar

Scassellati, C., Bonvicini, C., Faraone, S. V., and Gennarelli, M. (2012). Biomarkers and attention-deficit/hyperactivity disorder: A systematic review and meta-analyses. J. Am. Acad. Child Adolesc. Psychiatry 2012:015. doi: 10.1016/j.jaac.2012.08.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, J., Ordaz, S. J., Kircanski, K., Ho, T. C., Davis, E. G., Camacho, M. C., et al. (2019). Resting-state functional connectivity and inflexibility of daily emotions in major depression. J. Affect. Disord. 2019:040. doi: 10.1016/j.jad.2019.01.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Scofield, I. D., Johnson, J. D., Wood, P. K., and Geary, D. C. (2019). Latent resting-state network dynamics in boys and girls with attention-deficit/hyperactivity disorder. PLoS One 2019:0218891. doi: 10.1371/journal.pone.0218891

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, X., Tokoglu, F., Papademetris, X., and Constable, R. T. (2013). Groupwise whole-brain parcellation from resting-state fMRI data for network node identification. Neuroimage 82, 403–415. doi: 10.1016/j.neuroimage.2013.05.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M. (2012). The future of FMRI connectivity. Neuroimage 2012:022. doi: 10.1016/j.neuroimage.2012.01.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Sporns, O. (2011). Networks of the Brain. Massachusetts, MA: Massachusetts Institute of Technology.

Google Scholar

Sylvester, C. M., Corbetta, M., Raichle, M. E., Rodebaugh, T. L., Schlaggar, B. L., Sheline, Y. I., et al. (2012). Functional network dysfunction in anxiety and anxiety disorders. Trends Neurosci. 2012:012. doi: 10.1016/j.tins.2012.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Taber-Thomas, B. C., Morales, S., Hillary, F. G., and Pérez-Edgar, K. E. (2016). Altered topography of intrinsic functional connectivity in childhood risk for social anxiety. Depress. Anxiety 33, 995–1004. doi: 10.1002/da.22508

PubMed Abstract | CrossRef Full Text | Google Scholar

Tatemoto, K., Carlquist, M., and Mutt, V. (1982). Neuropeptide Y—a novel brain peptide with structural similarities to peptide YY and pancreatic polypeptide. Nature 296, 659–660. doi: 10.1038/296659a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Tegelbeckers, J., Bunzeck, N., Duzel, E., Bonath, B., Flechtner, H. H., and Krauel, K. (2015). Altered salience processing in attention deficit hyperactivity disorder. Hum. Brain Mapp. 36, 2049–2060. doi: 10.1002/hbm.22755

PubMed Abstract | CrossRef Full Text | Google Scholar

Thome, J., Frewen, P., Daniels, J. K., Densmore, M., and Lanius, R. A. (2014). Altered connectivity within the salience network during direct eye gaze in PTSD. Borderline Personal. Disord. Emot. Dysregulation 1, 17. doi: 10.1186/2051-6673-1-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Uddin, L. Q., Kelly, A. M., Biswal, B. B., Castellanos, F. X., and Milham, M. P. (2009). Functional connectivity of default mode network components: correlation, anticorrelation, and causality. Hum. Brain Mapp. 30, 625–637. doi: 10.1002/hbm.20531

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Den Heuvel, J. K., Furman, K., Gumbs, M. C. R. R., Eggels, L., Opland, D. M., Land, B. B., et al. (2015). Neuropeptide Y Activity in the Nucleus Accumbens Modulates Feeding Behavior and Neuronal Activity. Biol. Psychiatry 77, 633–641. doi: 10.1016/j.biopsych.2014.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Heuvel, M. P., and Hulshoff Pol, H. E. (2010). Exploring the brain network: A review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 20, 519–534. doi: 10.1016/j.euroneuro.2010.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

von Rhein, D., Beckmann, C. F., Franke, B., Oosterlaan, J., Heslenfeld, D. J., Hoekstra, P. J., et al. (2019). Network-level assessment of reward-related activation in patients with ADHD and healthy individuals. Hum. Brain Mapp. 38, 2359–2369. doi: 10.1002/hbm.23522

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, G., De La Cruz, F., Köhler, S., and Bär, K.-J. (2017). Treatment Associated Changes of Functional Connectivity of Midbrain/Brainstem Nuclei in Major Depressive Disorder OPEN. Nature 7, 9077–9075. doi: 10.1038/s41598-017-09077-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Waites, A. B., Stanislavsky, A., Abbott, D. F., and Jackson, G. D. (2005). Effect of prior cognitive state on resting state networks measured with functional connectivity. Hum. Brain Mapp. 24, 59–68. doi: 10.1002/hbm.20069

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Zuo, X., and He, Y. (2010). Graph-based network analysis of resting-state functional MRI. Front. Syst. Neurosci. 4:16. doi: 10.3389/fnsys.2010.00016

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhu, C., He, Y., Zang, Y., Cao, Q., Zhang, H., et al. (2009). Altered small-world brain functional networks in children with attention-deficit/hyperactivity disorder. Hum. Brain Mapp. 30, 638–649. doi: 10.1002/hbm.20530

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. H., Jiao, Y., and Li, L. (2018). Identifying individuals with attention deficit hyperactivity disorder based on temporal variability of dynamic functional connectivity. Sci. Rep. 81, 1–12. doi: 10.1038/s41598-018-30308-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Öngür, D., Auerbach, R. P., and Yao, S. (2016). Cognitive vulnerability to major depression: View from the intrinsic network and cross-network interactions. Harv. Rev. Psychiatry 2016:081. doi: 10.1097/HRP.0000000000000081

PubMed Abstract | CrossRef Full Text | Google Scholar

Warthen, K. G., Sanford, B., Walker, K., Jones, K. G., Angstadt, M., Sripada, C., et al. (2018). Neuropeptide Y and representation of salience in human nucleus accumbens. Neuropsychopharmacology 44, 495–502. doi: 10.1038/s41386-018-0230-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Widdowson, P. S., Ordway, G. A., and Halaris, A. E. (1992). Reduced Neuropeptide Y Concentrations in Suicide Brain. J. Neurochem. 59, 73–80. doi: 10.1111/j.1471-4159.1992.tb08877.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wise, T., Marwood, L., Perkins, A. M., Herane-Vives, A., Joules, R., Lythgoe, D. J., et al. (2017). Instability of default mode network connectivity in major depression: A two-sample confirmation study. Transl. Psychiatry 7:40. doi: 10.1038/tp.2017.40

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Van Dam, N. T., Feng, C., Luo, Y., Ai, H., Gu, R., et al. (2019). Anxious brain networks: A coordinate-based activation likelihood estimation meta-analysis of resting-state functional connectivity studies in anxiety. Neurosci. Biobehav. Rev. 2019:005. doi: 10.1016/j.neubiorev.2018.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Yehuda, R., Brand, S., and Yang, R.-K. (2006). Plasma Neuropeptide Y Concentrations in Combat Exposed Veterans: Relationship to Trauma Exposure, Recovery from PTSD, and Coping. Biol. Psychiatry 59, 660–663. doi: 10.1016/J.BIOPSYCH.2005.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeo, B. T. T., Krienen, F. M., Sepulcre, J., Sabuncu, M. R., Lashkari, D., Hollinshead, M., et al. (2011). The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. 106, 1125–1165. doi: 10.1152/jn.00338.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: neuropeptide Y, dynamic connectivity, LEiDA, graph theory, fMRI

Citation: Warthen KG, Welsh RC, Sanford B, Koppelmans V, Burmeister M and Mickey BJ (2021) Neuropeptide Y Variation Is Associated With Altered Static and Dynamic Functional Connectivity of the Salience Network. Front. Syst. Neurosci. 15:629488. doi: 10.3389/fnsys.2021.629488

Received: 15 November 2020; Accepted: 25 October 2021;
Published: 18 November 2021.

Edited by:

Alessandro E. P. Villa, Neuro-Heuristic Research Group (NHRG), Switzerland

Reviewed by:

Behnaz Jarrahi, Stanford University, United States
Laurel Morris, Icahn School of Medicine at Mount Sinai, United States

Copyright © 2021 Warthen, Welsh, Sanford, Koppelmans, Burmeister and Mickey. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Robert C. Welsh, Robert.c.welsh@utah.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.