Characterizing functional pathways of the human olfactory system
Abstract
The central processing pathways of the human olfactory system are not fully understood. The olfactory bulb projects directly to a number of cortical brain structures, but the distinct networks formed by projections from each of these structures to the rest of the brain have not been well-defined. Here, we used functional magnetic resonance imaging and k-means clustering to parcellate human primary olfactory cortex into clusters based on whole-brain functional connectivity patterns. Resulting clusters accurately corresponded to anterior olfactory nucleus, olfactory tubercle, and frontal and temporal piriform cortices, suggesting dissociable whole-brain networks formed by the subregions of primary olfactory cortex. This result was replicated in an independent data set. We then characterized the unique functional connectivity profiles of each subregion, producing a map of the large-scale processing pathways of the human olfactory system. These results provide insight into the functional and anatomical organization of the human olfactory system.
https://doi.org/10.7554/eLife.47177.001Introduction
The human sense of smell serves a variety of important functions in everyday life (Bushdid et al., 2014; Devanand et al., 2015; McGann, 2017). It is used to monitor the safety of inhaled air (Pence et al., 2014) and edibility of food (Yeomans, 2006). It also strongly impacts our social and emotional lives (Durand et al., 2013; Endevelt-Shapira et al., 2018; Frumin et al., 2015; Gelstein et al., 2011; Krusemark et al., 2013; Walla, 2008; Walla et al., 2003). Thus, the brain must extract different types of information from odor stimuli, including information about the identities of objects and foods, environmental hazards, and social and emotional cues. These functions are likely carried out by distinct cortical networks within the olfactory system, yet the organization of these functional networks is not fully understood.
This incomplete understanding is due partly to ambiguity about the anatomical and functional properties of the cortical targets of human olfactory bulb projections. Collectively, these areas are commonly referred to as primary olfactory cortex (Carmichael et al., 1994; Feher and Feher, 2017; Gottfried, 2010; Mai and Paxinos, 2012; Price, 2009) (although see Wilson, 2009; Haberly, 2001; Chapuis and Wilson, 2011 and others for discussions of the accuracy of this definition of the pirmary olfactory cortex). In humans, this includes the anterior olfactory nucleus, the olfactory tubercle, the frontal and temporal piriform cortices, and subregions of both the amygdala and entorhinal cortex (Allison, 1954; Eslinger et al., 1982; Gonçalves Pereira et al., 2005; Insausti et al., 2002; Milardi et al., 2017). The fact that the olfactory bulb simultaneously projects directly to a number of structures suggests parallel functional pathways in the olfactory system (Haberly, 2001; Kauer, 1991), but the distinct roles of these primary olfactory areas and their functional pathways are not fully understood (Bensafi et al., 2007; Gottfried et al., 2006; Gottfried et al., 2004; Gottfried et al., 2002; Howard et al., 2009; Li et al., 2008; Li et al., 2006; Sobel et al., 2000; Sobel et al., 1999; Zelano et al., 2005). Additionally, the olfactory system is organized differently than other sensory systems, which contain pre-cortical thalamic relays, further suggesting a deeper understanding of the organization of olfactory networks in the human brain is warranted.
The vast majority of research on primary olfactory cortex has focused on piriform cortex, which is the largest recipient of bulbar projections. Most of this research has been conducted in rodents, where piriform cortex is divided into anatomically and functionally distinct anterior and posterior subdivisions (Calu et al., 2007; Grau-Perales et al., 2019; Haberly and Price, 1978; Stettler and Axel, 2009; Yang et al., 2017). In humans, the anatomy and functionality of piriform cortex is less understood. Although it can be divided into frontal- and temporal-lobe subregions (Mai et al., 2015; Vaughan and Jackson, 2014; Young et al., 2018; Allison, 1954), whether these correspond to rodent anterior and posterior subdivisions is unclear. While neuroimaging studies have pointed to functional heterogeneity within human piriform cortex (Bensafi, 2012; Fournel et al., 2016; Gottfried et al., 2002; Howard et al., 2009; Howard and Gottfried, 2014; Li et al., 2008; Porter et al., 2005; Seubert et al., 2013; Zelano et al., 2011; Zelano et al., 2005), its anatomical and functional distinctions are still not clearly defined.
While numerous rodent and human studies have focused on piriform cortex, far fewer have examined other primary olfactory structures, such as the anterior olfactory nucleus and the olfactory tubercle. These structures have been anatomically well-defined in rodents (Aqrabawi and Kim, 2018a; Haberly and Price, 1978; Shipley and Adamek, 1984), primates (Carmichael et al., 1994), and humans (Allison, 1954; Eslinger et al., 1982; Mai et al., 2015), but their roles in olfactory processing are not fully understood in any of these species (Gadziola et al., 2015; Wesson and Wilson, 2011). Recent rodent data suggest that the anterior olfactory nucleus may be involved in odor memory (Aqrabawi and Kim, 2018b; Oettl et al., 2016) and localization (Kikuta et al., 2010), and the olfactory tubercle may play an important role in multisensory integration and attention (Wesson and Wilson, 2010; Zelano et al., 2005), although a complete understanding of the functions of these areas is lacking.
Previous studies have used task-related and resting functional magnetic resonance imaging (fMRI) to examine olfactory networks, using primary olfactory and orbitofrontal cortices as seed regions (Banks et al., 2016; Cecchetto et al., 2019; Fjaeldstad et al., 2017; Karunanayaka et al., 2017; Karunanayaka et al., 2014; Kiparizoska and Ikuta, 2017; Kollndorfer et al., 2015; Krusemark and Li, 2012; Nigri et al., 2013; Sreenivasan et al., 2017; Sunwoo et al., 2015). These studies have contributed important broad knowledge of parallel olfactory networks (Karunanayaka et al., 2014), how they compare to trigeminal networks (Karunanayaka et al., 2017), and how they change with age (Wang et al., 2005) and disease (Caffo et al., 2010; Fjaeldstad et al., 2017; Killgore et al., 2013; Sunwoo et al., 2015; Wang et al., 2010; Wang et al., 2015). However, the functional connectivity profiles of the primary olfactory subregions have not been considered separately. This is important because these subregions, which receive direct and parallel input from the bulb, likely form the anatomical substrates of ethological, parallel olfactory networks. Therefore, a quantitative characterization of the distinct functional pathways of human primary olfactory subregions would be an important step toward understanding the large-scale networks that underlie the basic, parallel, purposes of olfactory processing. The discovery of unique whole-brain connectivity profiles for the different primary olfactory subregions could also provide insight into the nature of the distinct functions of these areas in olfactory perception. This information, in turn, could have clinical implications for diseases that impact particular primary olfactory subregions.
Thus, the goals of this study were two-fold: first, to test the hypothesis that primary olfactory subregions form distinct large-scale olfactory processing networks; and second, if so, to characterize these networks across the whole brain. For the first, we used well-established, unsupervised k-means clustering techniques (Eickhoff et al., 2018; Kahnt et al., 2012; Kahnt and Tobler, 2017; Wang et al., 2017), to parcellate primary olfactory cortex into distinct clusters based solely on whole-brain connectivity patterns. We reasoned that if whole-brain functional connectivity patterns alone could be used to accurately parcellate primary olfactory cortex into its established, anatomically-defined subregions, this would suggest that these subregions do in fact form distinct functional pathways. We also reasoned that if these parcellation results were robust, the results should replicate in an independent data set. For the second, we characterized the distinct functional connectivity patterns of each primary olfactory subregion in order to produce a whole-brain map of the networks formed by each area. Our results provide insight into the functional and anatomical organization of the human olfactory system and provide a basis for future investigation into the functions of the distinct cortical targets of the olfactory bulb.
Results
We used resting-state fMRI connectivity to examine the functional pathways of human primary olfactory cortex in two main steps. First, we used k-means clustering techniques to parcellate primary olfactory cortex into distinct clusters. These clusters were based on the group-level, whole-brain functional connectivity of all voxels within primary olfactory cortex, separately for each hemisphere. Second, using the results of the parcellation analysis, we characterized the distinct, large-scale networks of the human olfactory system. To do this, we first determined that there were no hemispheric differences in the connectivity profiles of primary olfactory subregions, suggesting we should combine corresponding clusters across the left and right hemispheres. We then quantified the whole-brain functional networks that were unique to each subdivision, and those that were common to all subdivisions. Finally, as a discussion point, we considered the functional properties of connected brain areas for each primary olfactory subdivision and attempted to form a speculative hypothetical model of human olfactory functional networks.
Parcellation of human primary olfactory cortex
To test the hypothesis that primary olfactory subregions form distinct, large-scale olfactory networks, we tested whether their anatomical boundaries could be accurately delineated based on whole-brain functional connectivity maps. To do this, we conducted a functional-connectivity-based parcellation of human primary olfactory cortex. Twenty-five subjects (average ± standard error age: 25.5 ± 1.2 years; 14 female) underwent a 10 min resting-state fMRI scan. We first outlined the entirety of primary olfactory cortex into a combined region-of-interest (ROI) on which to perform the k-means clustering analysis. This ROI was drawn for the left and right hemispheres separately, according to a human brain atlas which contains detailed demarcation of most primary olfactory areas (Mai et al., 2015; Ongür et al., 2003) (Figure 1A). The ROI included only those subregions with detailed boundaries in the atlas, and consisted of a combination of the anterior olfactory nucleus, olfactory tubercle, and frontal and temporal piriform cortices (Figure 1B), defined based on Mai et al. (2015). Note that additional primary olfactory areas, including amygdala and entorhinal cortex (Allison, 1954; Carmichael et al., 1994; Eslinger et al., 1982; Gonçalves Pereira et al., 2005; Zatorre et al., 1992), were not included in our ROI because the exact location of olfactory afferents into these areas is poorly understood (Gonçalves Pereira et al., 2005). This is an important topic for future investigation of human olfactory networks.

Region of interest.
(A) Panels show examples from the human brain atlas used to define the region of interest used in the parcellation analysis. Relevant areas include the anterior olfactory nucleus (AON), olfactory tubercle (TUB), and frontal (PirF) and temporal (PirT) piriform cortex (Mai et al., 2015). (B) The region of interest shown overlaid on the FSL’s MNI152_T1_1mm_brain. The coronal and axial slices correspond to the vertical and horizontal lines on the sagittal slice respectively. R, right hemisphere.
-
Figure 1—source data 1
Relates to Figure 1.
Atlas and region-of-interest of the human primary olfactory cortex.
- https://doi.org/10.7554/eLife.47177.003
We then estimated the whole-brain functional connectivity profile of each voxel within the ROI by computing the Pearson correlation coefficient between the resting-state fMRI time-series of a given voxel and that of every other voxel in the rest of the brain. This resulted in subject-wise connectivity matrices. We then performed a leave-one-out analysis (Kahnt et al., 2012) to estimate the stability of the connectivity profiles of individual primary olfactory cortex voxels across participants, as a prerequisite for averaging the connectivity matrices across subjects (Figure 2A). To examine the similarity between the individual functional connectivity matrices, we computed a histogram of correlation values between individual matrices and the group matrix. The histogram of correlation values showed that the similarity of connectivity patterns was above zero in all voxels in primary olfactory cortex (mean correlation coefficient: 0.19, standard error: 0.0041), justifying averaging across subjects. Note that because they are computed across rest-of-the-brain voxels, R values larger than 0.0088 are statistically significant at p<0.05 (Bonferroni corrected for the number of voxels in primary olfactory cortex). To parcellate within-ROI voxels into subdivisions based on their whole-brain functional connectivity profiles, we applied unsupervised k-means clustering methods to the average connectivity matrix. We used a priori k = 4 based on the fact that our ROI was comprised of four anatomically distinct brain regions. For both the left and right hemispheres, this analysis successfully parcellated the primary olfactory cortex ROI into four distinct brain regions that corresponded anatomically to the anterior olfactory nucleus, olfactory tubercle, and frontal and temporal piriform cortices (Figure 2B).

Parcellation of human left and right primary olfactory cortex.
(A) Inter-subject stability of functional connectivity patterns. The correlation of the functional connectivity patterns between each subject and all other subjects was calculated for each voxel using a leave-one-out method. The coronal slices, corresponding to the vertical lines on the sagittal slice, show the average stability map. The bar plot shows the histogram of the correlation values. (B) k-means (k = 4) clustering results shown on the FSL’s MNI152_T1_1mm_brain. The right column shows one axial and one coronal slice of the Atlas (Mai et al., 2015). (C) Parcellation accuracy of each subregion. Left column: proportion of voxels from each parcellation subdivision located within each anatomical subregion. Right column: z score of the proportion maps. * indicates p<0.001 (false discovery rate corrected). R, right hemisphere; AON, anterior olfactory tubercle; TUB, olfactory tubercle; PirF, frontal piriform cortex; PirT, temporal piriform cortex.
-
Figure 2—source data 1
Relates to Figure 2, panel (A).
Correlation between individual and group-level (leave-one-out) connectivity profiles.
- https://doi.org/10.7554/eLife.47177.008
-
Figure 2—source data 2
Relates to Figure 2, panel (B) and (C).
Manual segmentation and parcellation of the human primary olfactory cortex.
- https://doi.org/10.7554/eLife.47177.009
To confirm the correspondence between our parcellation results and the anatomical delineation of primary olfactory subregions in the Atlas of the Human Brain (Mai et al., 2015), we computed the proportion of voxels from each parcellation cluster located within each of the atlas-derived subdivisions, drawn prior to performing the parcellation analysis (Figure 2B). The statistical significance of this proportion was tested using a permutation test. Specifically, for each permutation, we shuffled the labels of the anatomical subdivision and re-calculated the proportion. This procedure was repeated 10,000 times, resulting in a distribution of permuted proportions for each parcellation cluster. A z score of the actual proportion values was computed by subtracting the average and then dividing by the standard deviation, which was obtained by normal distribution fitting of the permuted data (Matlab’s normfit). We found that for each parcellated subdivision, there was one anatomical ROI that contained significantly more voxels than the other anatomical ROIs (Figure 2C, minimal z score = 7.18). Thus, we found that the location of voxels within each parcellated subdivision corresponded to a single anatomically-determined ROI. Specifically, the medial-rostral-most parcellated subdivision corresponded to the anterior olfactory nucleus. The adjacent caudal subdivision corresponded to the olfactory tubercle. Within the frontal lobe, the lateral-rostral subdivision corresponded to the frontal piriform cortex, and within the temporal lobe, the caudal-most subdivision corresponded to temporal piriform cortex.
Replication of parcellation results
To confirm the robustness of our parcellation results, we performed two control analyses aimed at replicating the initial findings. First, we performed the k-means clustering analysis on a different ROI of primary olfactory cortex, drawn independently by one of the co-authors of this paper. Second, we performed the k-means clustering analysis on an independent data set (N = 53), collected for a previous study on a different scanner, with different acquisition parameters and different subjects (Kahnt and Tobler, 2017).
In the first control analysis, performed on an independently drawn ROI, we found that k-means clustering still successfully parcellated primary olfactory cortex into four distinct regions, corresponding anatomically to the anterior olfactory nucleus, olfactory tubercle, and frontal and temporal piriform cortices (Figure 2—figure supplement 1A,B). In our second control analysis, performed on an independent data set, we found that, again, primary olfactory cortex successfully parcellated into the same four distinct regions (Figure 2—figure supplement 1C,D). Importantly, all analysis steps performed on this independent data set were identical to those performed in our initial analysis. These results suggest good reliability of our finding that human olfactory cortex can be accurately parcellated based on whole-brain functional connectivity patterns.
Parcellation results across hemispheres and k values
Thus far, we demonstrated that both the left and right primary olfactory areas can be accurately subdivided based on their functional connectivity profiles. To further examine differences between the left and right hemispheres and at different k values, we conducted additional parcellation analyses using a single primary olfactory ROI containing all subregions from both hemispheres, at a range of k values. We reasoned that if connectivity patterns were similar across hemispheres for each subregion, then parcellation analysis of this combined ROI should group left and right sides of each primary olfactory subregion, as opposed to grouping, for example, the neighboring subregions on the same hemisphere. We computed this analysis using k values ranging from 3 to 6 (Figure 3). We found that for a k value of 3, the parcellation analysis grouped left and right anterior olfactory nucleus and left tubercle as one cluster, left and right frontal piriform cortex and left temporal piriform cortex as a second cluster, and right temporal piriform cortex alone as the third cluster. For a k value of 4, the analysis successfully grouped left and right hemispheres for both piriform subregions, but it grouped left anterior olfactory nucleus with left olfactory tubercle and right anterior olfactory nucleus with right olfactory tubercle, suggesting these two primary olfactory areas have relatively more lateralized connectivity patterns. A k value of 5 grouped the left and right hemispheres of all subregions, with a fifth cluster consisting of only right temporal piriform cortex. Finally, a k value of 6 grouped all subregions across hemispheres except for the anterior olfactory nucleus. These results indicate a clear separation of frontal and temporal piriform cortex for a wide range of cluster solutions and even across hemispheres (Figure 3A). All parcellation results showed good agreement with the anatomical subregions (Figure 3B,C). Of note, the anterior olfactory nucleus and olfactory tubercle are classified as one subregion for a clustering solution of k = 3. The left and right anterior olfactory nucleus were separated into different subregions for k = 4, 6. These findings suggest stronger lateralization of connectivity patterns for the anterior olfactory nucleus compared to other primary olfactory areas.

Parcellation of primary olfactory cortex combined across left and right hemispheres.
(A) k-means clustering results shown on the FSL’s MNI152_T1_1mm_brain for k = 3 to 6. Each color represents one cluster. (B) Proportion of voxels of each parcellation subdivision within each anatomical subregion. (C) z score of the proportion maps in panel B. * indicates p<0.001 (false discovery rate corrected). L, left hemisphere; R, right hemisphere; AON, anterior olfactory tubercle; TUB, olfactory tubercle; PirF, frontal piriform cortex; PirT, temporal piriform cortex.
-
Figure 3—source data 1
Relates to Figure 3.
Parcellation result (k = 3 to 6) of left and right hemispheres combined ROI.
- https://doi.org/10.7554/eLife.47177.011
Primary olfactory cortical functional connectivity does not statistically differ between hemispheres
Our next step was to examine the whole-brain functional connectivity profiles of the different primary olfactory subregions. Prior to performing this analysis, we first determined whether to use hemispherically combined clusters, or hemispherically distinct clusters. We reasoned that if connectivity profiles of left and right primary olfactory areas did not statistically differ, then they should not be analyzed separately. We therefore conducted a lateralization-index analysis to directly statistically compare connectivity patterns across hemispheres. The lateralization index was defined as (Zleft – Zright)/(Zleft + Zright), where Zleft and Zright were the functional connectivity maps for the left and right seed regions, respectively. We found no statistically significant difference between whole-brain connectivity maps for the corresponding left and right primary olfactory subregions (Figure 4—figure supplement 1) (threshold-free cluster enhancement (TFCE) corrected p>0.001). Based on this result, all proceeding analyses were conducted using combined ROIs of corresponding left and right subregions.
Distinct whole-brain human olfactory networks
The fact that primary olfactory subregions were accurately anatomically parcellated based on their functional connectivity profiles suggests that they form distinct, parallel olfactory networks. To examine these functional networks, we produced whole-brain maps of the non-overlapping brain areas exhibiting functional connectivity with each subregion. To do this, we first applied a statistical threshold to the whole-brain functional connectivity map for each subregion (TFCE corrected p<0.001), and binarized them (assigned a value of 1 or 0). This resulted in four distinct whole-brain maps. We then further masked them according to whether each voxel exhibited statistically significant functional connectivity with a single subregion or with multiple subregions. This produced two maps of functional connectivity: one with the unique connectivity patterns for each subregion, and the other with connectivity patterns shared by multiple subregions. The complete list of areas showing subregion-specific connectivity is contained in Table 1.
Summary of functional connectivity results.
https://doi.org/10.7554/eLife.47177.012Volume (mm3) | |||||
---|---|---|---|---|---|
Label | Overlap | AON | TUB | PirF | PirT |
Frontal Orbital Cortex | 2592 | 3000 | 144 | - | 1048 |
Frontal Medial Cortex | 992 | 1120 | - | - | - |
Cingulate Gyrus | 2760 | 200 | 80 | 2304 | - |
Insular Cortex | 384 | 496 | - | 224 | 832 |
Subcallosal Cortex | 3632 | 616 | - | - | - |
Caudate | 136 | 120 | - | 2024 | - |
Paracingulate Gyrus | 1336 | - | 2600 | - | - |
Parahippocampal Gyrus | 296 | - | 464 | - | 2584 |
Temporal Pole | 328 | - | - | - | 9184 |
Putamen | 1368 | - | - | 3376 | 96 |
Hippocampus | 1176 | - | - | 136 | 1448 |
Amygdala | 2120 | - | - | - | - |
Accumbens | 336 | - | - | - | - |
Planum Polare | - | 248 | - | - | 480 |
Frontal Pole | - | 2792 | 1504 | - | 736 |
Temporal Fusiform Cortex | - | 688 | 352 | - | 1240 |
Inferior Frontal Gyrus | - | 208 | - | - | - |
Inferior Temporal Gyrus | - | 1224 | - | - | 248 |
Heschl's Gyrus (includes H1 and H2) | - | 80 | - | - | 208 |
Planum Temporale | - | 104 | - | - | 96 |
Brainstem | - | - | 592 | - | 6056 |
Thalamus | - | - | 120 | 1384 | - |
Pallidum | - | - | - | 504 | - |
Precentral Gyrus | - | - | - | 1616 | - |
Postcentral Gyrus | - | - | - | 216 | 296 |
Frontal Operculum Cortex | - | - | - | 128 | 336 |
Central Opercular Cortex | - | - | - | 224 | 376 |
Supramarginal Gyrus | - | - | - | 808 | - |
Juxtapositional Lobule Cortex | - | - | - | 1104 | - |
Superior Frontal Gyrus | - | - | - | - | 128 |
Temporal Occipital Fusiform Cortex | - | - | - | - | 344 |
Superior Temporal Gyrus | - | - | - | - | 1464 |
Middle Temporal Gyrus | - | - | - | - | 1368 |
Angular Gyrus | - | - | - | - | 256 |
Parietal Operculum Cortex | - | - | - | - | 640 |
-
The volumes of statistically significant voxels in each brain region are shown for overlapping and subregion-specific clusters. The Overlap column does not include subregion-specific regions. - indicates volume less than 80 mm3 (10 voxels). Atlas query was conducted with FSL's HarvardOxford-cort-maxprob-thr50-2mm and HarvardOxford-sub-maxprob-thr50-2mm atlases. AON, anterior olfactory nucleus; TUB, olfactory tubercle; PirF, frontal piriform cortex; PirT, temporal piriform cortex.
-
Table 1—source data 1
Functional connectivity profile of each subregion.
- https://doi.org/10.7554/eLife.47177.013
Below, we outline the unique connectivity patterns we found for each primary olfactory subregion.
Functional connectivity profiles of anterior olfactory nucleus and olfactory tubercle
The brain areas that showed connectivity unique to the anterior olfactory nucleus were largest in the orbitofrontal cortex, with smaller clusters in the left inferior temporal gyrus, bilateral anterior temporal gyri, the bilateral anterior insula and the mammillary bodies and retromammillary commissure (Figure 4A,E). Areas of connectivity in the orbitofrontal cortex were extensive, including the entire gyrus rectus and encompassing parts of the medial, anterior and lateral orbital gyri. Notably, there was a strong cluster of connectivity corresponding to bilateral areas along the medial orbital sulcus, close to its intersection with the transverse orbital sulcus—this is significant because this part of orbitofrontal cortex is sometimes referred to as human secondary olfactory cortex (Gottfried and Zald, 2005). Connectivity with the left inferior temporal gyrus was centered around a posterior region along the established object-recognition pathway. There were also large clusters in the bilateral anterior temporal gyri, and the bilateral anterior insula, in gustatory cortex. Finally, there was a medial cluster centered around the region between the mammillary bodies, the retromammillary commissure and the posterior hypothalamic nucleus.

Subregion-specific functional connectivity patterns.
(A–D) Brain regions that are uniquely positively connected to each of the primary olfactory subregions including the (A) anterior olfactory nucleus (AON, red), (B) olfactory tubercle (TUB, purple), (C) frontal piriform cortex (PirF, blue) and (D) temporal piriform cortex (PirT, green). The results are shown on the FSL’s MNI152_T1_1mm_brain. (E–H) Functional connectivity maps shown on Freesurfer's cvs_avg35_inMNI152 brain surface for AON, TUB, PirF and PirT. All functional connectivity maps were thresholded at threshold-free cluster enhancement corrected p<0.001. L, left hemisphere; PHN, posterior hypothalamic nucleus; ITG, inferior temporal gyrus; OFC, orbitofrontal cortex; INS, insula; PCC, paracingulate cortex; RN, retromammillary nucleus; ACCU, accumbens; RSC, retrosplenial area; FP, frontal pole; FFA, fusiform face areas; ACC, anterior cingulate cortex; THAL, thalamus; CAU, caudate; PUT, putamen; SUP, supramarginal gyrus; TP, temporal pole; HIPP, hippocampus; BA, Broca’s area.
-
Figure 4—source data 1
Relates to Figure 4.
Statistical images of the functional connectivity of human primary olfactory cortex subregions.
- https://doi.org/10.7554/eLife.47177.017
The brain areas that exhibited connectivity unique to the olfactory tubercle were largest in the medial prefrontal cortex, with smaller clusters in the perisplenial region, the left temporal fusiform cortex, the red nucleus of the brainstem and the accumbens (Figure 4B,F). Connectivity between the olfactory tubercle and the medial prefrontal cortex was mainly located in the paracingulate gyrus and frontal pole. In the paracingulate, it centered bilaterally around the anterior-most aspect of Brodman’s area 32. Connectivity with the frontal pole was stronger in the left hemisphere, and extended from the dorsal edge of the ventro-medial prefrontal cortex into the medial prefrontal cortex, with a smaller cluster extending even more dorsal, reaching the dorsomedial prefrontal cortex. There was also a cluster of connectivity located bilaterally in the ventral striatum, in the accumbens. There were bilateral clusters of connectivity in both the red nuclei of the brainstem and in the retrosplenial area.
Functional connectivity profiles of frontal and temporal piriform cortices
Interestingly, the frontal and temporal piriform cortices also showed distinct connectivity patterns, supporting the hypothesis that they are functionally distinct. The frontal piriform cortex showed a large cluster of connectivity with the dorsal striatum (both putamen and caudate), with smaller clusters in the precentral gyrus, the cingulate gyrus, mediodorsal thalamus and left supramarginal gyrus (Figure 4C,G). Connectivity with the caudate nucleus and putamen was extensive, covering the bilateral entirety of both regions. In the precentral gyrus, connectivity was evident along the primary motor strip, close to the lateral sulcus, corresponding to the face/lips/tongue/jaw area. There was also a large cluster of connectivity in the caudal anterior (anterior mid-cingulate) cortex, a cluster within the mediodorsal thalamic nucleus extending over the entire anterior-posterior axis, and a cluster unilaterally in the center of the left supramarginal gyrus.
The brain areas that exhibited unique connectivity to the temporal piriform cortex were largest in the bilateral brainstem (centered in the pons) and temporal pole, with smaller clusters in bilateral inferior frontal gyri, bilateral superior temporal gyrus, bilateral hippocampus and bilateral posterior insula (Figure 4D,H). In the brainstem, connectivity was maximal in the ventral aspect of the pons within the nucleus raphe magnus. There were also large clusters of connectivity in the anterior temporal pole. Within the inferior frontal gyrus, connectivity was stronger on the left side within Broca’s area, with a smaller cluster in the same area on the right hemisphere. There were clusters in the left superior temporal gyrus, in the language comprehension areas (Friederici et al., 2003; Leff et al., 2009) and clusters in the right medial temporal gyrus.
Functional connectivity common to all subregions
Although we have shown that primary olfactory subregions have distinct functional connectivity profiles, they likely also have some functional pathways in common, especially considering the strong reciprocal connectivity between the subregions. To identify the common primary olfactory connectivity network, the connectivity map of each subregion was binarized at a threshold of TFCE corrected p<0.001 to include only those clusters that exhibited connectivity with all subregions. This resulted in a whole-brain connectivity map of areas that are functionally connected to all primary olfactory cortical subregions (Figure 5, Table 1). We found that these areas included large clusters in the bilateral hippocampus, amygdala and subgenual area, with smaller clusters in the anterior insula and posterior orbitofrontal cortex (Figure 5A,B). Connectivity with the hippocampus covered the entire anterior-posterior extent. In the amygdala, connectivity was maximal in the medial subregions and in the central amygdala. Interestingly, for both the right and left primary olfactory cortical subregions, connectivity clusters with the rest of the brain were generally more extensive in the right hemisphere compared to the left hemisphere (Frasnelli et al., 2010; Herz et al., 1999; Royet and Plailly, 2004; Zatorre et al., 1992) (Figure 5A,B). The common connectivity of all primary olfactory subregions with the orbitofrontal cortex could explain why it is so reliably activated even in basic and passive olfactory tasks.

Functional connectivity common to all subregions.
(A) Brain regions that showed statistically significant positive functional connectivity with each of the primary olfactory subregions. Results are overlaid on the FSL’s sagittal and axial MNI152_T1_1mm_brain. (B) Same brain regions as in (A) shown on medial (top row) and lateral (bottom row) Freesurfer's cvs_avg35_inMNI152 brain surfaces. Red indicates the functional connectivity maps that were thresholded at threshold-free cluster enhancement corrected p<0.001. R, right hemisphere.
-
Figure 5—source data 1
Relates to Figure 5.
Functional connectivity patterns common to all human primary olfactory cortex subregions.
- https://doi.org/10.7554/eLife.47177.019
Discussion
In the olfactory system, information flows from the bulb to the cortex in a direct and parallel manner. Similarly, other sensory systems have parallel organization of processing pathways, but these typically occur at a later stage of processing. For example, information from other systems, including somatosensation, gustation, vision and audition, is processed in the thalamus prior to primary sensory cortex, and is parallelized downstream from there. In olfaction, thalamic processing occurs after direct parallel primary cortical processing, and, based on data from this study, may only occur for a subset of olfactory processing streams. From the olfactory bulb, mitral and tufted cells project to several different brain regions which are thought to play unique roles in olfactory processing (Ghosh et al., 2011; Haberly, 2001; Miyamichi et al., 2011; Sosulski et al., 2011). In rodents, the anatomical locations and properties of these regions have been well defined (Igarashi et al., 2012; Nagayama et al., 2010; Vassar et al., 1994). In humans, however, far fewer studies have attempted to identify and define these primary olfactory areas (Crosby and Humphrey, 1939; Eslinger et al., 1982; Allison, 1954; Shipley and Reyes, 1991). Similarly, whereas many studies have begun to outline unique functional properties of primary olfactory areas within the rodent olfactory system (Ikemoto, 2003; Illig, 2005; Lei et al., 2006; Myhrer et al., 2010; Wesson and Wilson, 2010), fewer have done the same in humans. In the current study, we used data-driven connectivity-based parcellation techniques to show that whole-brain functional connectivity patterns alone could be used to parcellate human primary olfactory cortical regions into subregions that anatomically corresponded to the anterior olfactory nucleus, olfactory tubercle and frontal and temporal piriform cortices. Our findings indicate a dissociation in whole-brain functional connectivity patterns across the subregions of human primary olfactory cortex. This suggests that the human olfactory system is comprised of distinct parallel processing pathways (Cecchetto et al., 2019; Savic et al., 2000) that may be related to the different recipients of projections from the olfactory bulb, in line with previous studies in rodents (Geramita et al., 2016; Haberly, 2001; Igarashi et al., 2012; Kauer, 1991; Payton et al., 2012). Our findings provide a detailed description of the particular brain areas that exhibit unique connectivity with each individual primary olfactory subregion. This could be used to gain insight into the specific role that each subregion plays in olfactory perception. We also demonstrate parallel organization of the olfactory system, in which olfactory networks reach a broader set of cortical targets at an earlier stage of processing compared to other sensory systems. This data might have implications for understanding the olfactory decline that appears in early stages of some neurological disease states. For example, olfactory structures that form networks with brain areas implicated in particular pathologies likely perform critical olfactory sensory functions (e.g. the olfactory tubercle is connected to areas involved in depression, see Croy and Hummel, 2017, and the temporal piriform cortex is connected to areas involved in primary progressive aphasia, see Olofsson et al., 2013).
Across species, the most frequently studied primary olfactory area is piriform cortex. In rodents, it has been consistently shown that projections from the bulb to piriform cortex are spatially distributed (Ghosh et al., 2011; Iurilli and Datta, 2017; Miyamichi et al., 2011) and that information about odor identity can be extracted from the spatiotemporal dynamics of these ensemble patterns (Haddad et al., 2013; Illig and Haberly, 2003; Poo and Isaacson, 2009; Rennaker et al., 2007; Stettler and Axel, 2009; Sugai et al., 2005; Zhan and Luo, 2010). At the same time, numerous studies suggest that the function of piriform cortex goes beyond simple odor-identity coding and is strongly impacted by its regional connectivity with other cortical areas (Cleland and Linster, 2003; Sadrian and Wilson, 2015). Numerous studies have implicated posterior piriform cortex in associative functions such as odor learning and memory (Calu et al., 2007; Chen et al., 2014; Choi et al., 2011; Gire et al., 2013; Gottfried and Dolan, 2003; Johnson et al., 2000; Karunanayaka et al., 2015; Martin et al., 2006; Roesch et al., 2007; Sacco and Sacchetti, 2010; Schoenbaum and Eichenbaum, 1995), and also suggest that the region may mediate learned olfactory responses and behaviors (Choi et al., 2011). Moreover, the strength and composition of piriform networks has also been shown to depend on experience and on the state of the organism (Chapuis et al., 2013; Cohen et al., 2015; Cohen et al., 2008; Hasselmo and Barkai, 1995; Kay and Freeman, 1998; Linster and Hasselmo, 2001; Wilson and Sullivan, 2011).
While human studies suggest a role for piriform cortex that goes beyond pure odor-object coding, (Bensafi, 2012; Gottfried, 2010; Porter et al., 2007; Schulze et al., 2017; Zelano et al., 2005), the anatomical and functional properties of human piriform cortex are not well understood compared to rodents. In fact, in some cases, it is not even clear which human bulbar recipients correspond to the well-defined rodent ones. While it has been suggested that the frontal and temporal subdivisions of human piriform may correspond to the anterior and posterior subdivisions of rodent piriform respectively, there is little anatomical or functional evidence to support this. For example, the rodent anterior piriform receives a greater density of bulb projections than the posterior piriform, while the same is not true for human frontal, relative to temporal, piriform. The lack of clear correspondence between human and rodent piriform subdivisions suggest that functional differences in primary olfactory areas across species could be substantial.
In our study, we found a clear differentiation between the functional connectivity patterns of human frontal and temporal piriform subdivisions, suggesting these two areas play different roles in olfactory processing in the human brain (Albrecht et al., 2010; Bao et al., 2016; Bensafi et al., 2007; Howard et al., 2009; Plailly et al., 2012; Porter et al., 2005; Zelano et al., 2005). This distinction between frontal and temporal piriform cortices was robust, surviving across k values and hemispheres. Interestingly, our findings suggest that frontal piriform cortex has strong functional connectivity with motor planning areas, including the caudate/putamen and the primary motor cortex, specifically at the face/nose/jaw section of the motor homunculus (perhaps the facial movement areas include sniffing). Intriguingly, frontal piriform cortex was also connected to the left supramarginal gyrus, an area that has been consistently implicated in tool-grasping in humans (Glover et al., 2012; Johnson-Frey et al., 2005), leading us to the tempting and highly speculative hypothesis that frontal piriform cortex may play a specific role in combining olfactory information with motor planning in order to guide food with the hand to the mouth. In contrast, we found that the temporal piriform cortex was connected to the brainstem raphe magnus and posterior insula, areas implicated in pain processing (Segerdahl et al., 2015; Woo et al., 2009) and respiratory modulations (Ackermann and Riecker, 2010; Evans et al., 2009), as well as the core language network (Ardila et al., 2014; Wible et al., 2005). Interestingly, in the context of olfaction, respiratory modulation and pain mediation are tightly linked, since many olfactory stimuli also activate trigeminal nerve endings that are located inside the nasal cavities. If a painful stimulus enters the nasal cavities (e.g. ammonia), a protective, fast respiratory reduction occurs, to minimize entry of dangerous chemicals into the body. Intriguingly, verbal communication also requires modulations in breathing, and thus, temporal piriform cortex may mediate both olfactory-related verbal communication and protective changes in respiration. These findings may be applicable to the hypothesized neurocognitive limitations of olfactory language (Cain, 1979; Engen and Ross, 1973; Olofsson and Gottfried, 2015), since our data may suggest that the same olfactory subregion exhibiting connectivity with language networks is involved in other critical olfactory functions.
Beyond piriform cortex, other cortical recipients of olfactory bulb output are less explored, including the anterior olfactory nucleus (Brunjes et al., 2005) and olfactory tubercle (Wesson and Wilson, 2011). In rodents, the anterior olfactory nucleus is a true cortical structure that can be divided into two main subdivisions, the pars externa and pars principalis (Pigache, 1970; Valverde et al., 1989). Studies suggest that some of the spatial properties of rodent glomerular activation are preserved in the pars externa, whereas the pars principalis exhibits a spatially distributed activation similar to other primary olfactory areas (Kay et al., 2011; Miyamichi et al., 2011). This suggests functional heterogeneity across the rodent anterior olfactory nucleus. In humans, the anterior olfactory nucleus is also a true cortical structure, but it does not appear to contain an analog to the pars externa (Crosby and Humphrey, 1941). It can be divided anatomically into retrobulbar and cortical anterior and posterior subdivisions (Mai et al., 2015; Ubeda-Bañon et al., 2017; Ongür et al., 2003). That said, only a small number of studies have directly explored the human anterior olfactory nucleus (Ubeda-Bañon et al., 2017), so its functional subdivisions remain unclear. Rodent studies suggest the pars externa is involved in odor localization (Esquivelzeta Rabell et al., 2017; Kikuta et al., 2010) and the pars medialis could be involved in top-down modulation of bulbar responses (Aqrabawi et al., 2016). More broadly, rodent studies suggest a role for the anterior olfactory nucleus in the initial formation of representations of odor objects (gestalts) (Haberly, 2001; Lei et al., 2006), allowing piriform cortex to perform more associative functions, relating information about odor objects to movement and behavior. In our study, we found that whole-brain connectivity of the human anterior olfactory nucleus was maximal in orbitofrontal areas and other regions associated with object recognition. Our findings are in line with rodent studies showing connectivity between anterior olfactory nucleus and orbitofrontal cortex (Illig, 2005), and with the hypothesized role for the anterior olfactory nucleus in the formation of odor-object representations (Haberly, 2001). We also found that connectivity profiles for the anterior olfactory nucleus were more lateralized than those of other primary olfactory subregions. This finding is in line with rodent studies suggesting that anterior olfactory nucleus neurons can distinguish between signals from the ipsilateral and contralateral nostrils, suggesting representation of lateralized inputs in this region (Kikuta et al., 2010).
Similar to the anterior olfactory nucleus, we lack a complete understanding of the function of the olfactory tubercle, especially in humans. The majority of research on the olfactory tubercle has been conducted in rodents and has focused on its relationship with the reward system (Ikemoto, 2007), with far fewer studies considering it as a primary olfactory cortical structure (Wesson and Wilson, 2011). Interestingly, the olfactory tubercle is the main recipient of tufted (as opposed to mitral) cell projections from the rodent olfactory bulb (Scott et al., 1980). Tufted cells show enhanced odor sensitivity, enhanced respiratory entrainment and broader receptive fields compared to mitral cells (Mori and Shepherd, 1994; Shepherd et al., 2004), suggesting that the olfactory tubercle is important for olfactory tasks requiring high sensitivity. Studies suggest potential roles for the rodent olfactory tubercle in odor discrimination (Murakami et al., 2005; Wesson and Wilson, 2010), olfactory multisensory integration (Wang et al., 2010), state-dependent modulation of olfactory bulb activity (Gervais, 1979) and odor reward value (Howard et al., 2016). Studies also suggest a role for the rodent olfactory tubercle in sensory hedonics and social behavior (Gervais, 1979; Hagamen et al., 1977; Hitt et al., 1973). Compared to rodents, much less is known about the human olfactory tubercle. In our study, we found connectivity between the olfactory tubercle and brain areas implicated in emotional processing, depression, and social cognition, including anterior paracingulate cortex and left frontal pole (Bludau et al., 2014; Eskenazi et al., 2015; Fettes et al., 2018; Jackson et al., 2003; Koch et al., 2018; Koechlin, 2011; Papmeyer et al., 2015; Veer et al., 2010). Interestingly, we also found connectivity between the olfactory tubercle and the left fusiform gyrus (fusiform face area), a brain region that is highly responsive to human faces (Çukur et al., 2013; Grill-Spector et al., 2004; Kanwisher et al., 1997; McCarthy et al., 1997). Our findings are consistent with rodent data suggesting a role for the tubercle in emotion and social cognition.
In an attempt to summarize and clarify the basic implications of our findings, we created a simplified, speculative illustrative diagram suggesting potential roles for each of the primary olfactory subdivisions that our study parcellated (Figure 6). One possible interpretation of our data suggests the following basic functions for human primary olfactory cortical subregions: Anterior olfactory nucleus, a cortical structure, may function as a first step in forming olfactory object representations; The olfactory tubercle may function in extracting social and emotional information from olfactory stimuli, mediating olfactory-related social and emotional responses; Frontal piriform cortex may be involved in planning olfactory-related movements, speculatively those involved in eating; Temporal piriform cortex may mediate pain-related (protective) respiratory modulations and verbal communication about olfactory stimuli. Importantly, the functions of primary olfactory areas cannot be determined from resting connectivity networks explored here, and there are numerous alternate interpretations of our data. Future experiments that combine psychophysics with the measurement of neural activity across primary olfactory areas are needed to gain a full understanding of parallel functional olfactory pathways in the human brain.

Schematic illustrative summary of the functional organization of human primary olfactory cortex.
(A) Summary of brain regions that are uniquely connected to each subregion, including the anterior olfactory nucleus (AON), olfactory tubercle (TUB), frontal piriform cortex (PirF), and temporal piriform cortex (PirT). (B) Summary of brain regions that are commonly connected to all subregions. mPFC, medial prefrontal cortex; CAU, caudate; PUT, putamen; AMC, anterior mid-cingulate cortex; THAL, thalamus; RN, red nucleus; RSC, retrosplenial cortex; ITG, inferior temporal gyrus; pINS, posterior insular cortex; MOT, motor area; TP, temporal pole; BA, Broca's area; aINS, anterior insular cortex; OFC, orbitofrontal cortex; ACC, anterior cingulate cortex; SC, subcallosal cortex; HIPP, hippocampus; AMY, amygdala.
-
Figure 6—source data 1
Relates to Figure 6.
Source files for constructing the schematic illustrative model.
- https://doi.org/10.7554/eLife.47177.021
Materials and methods
Participants
Request a detailed protocolTwenty-five healthy subjects (14 female) with normal or corrected-to-normal vision were recruited for this study. The average (standard error) age was 25.5 (1.2) years. All subjects were right-handed and had normal olfactory function by self-report. No participant reported a history of smell, taste or ear-nose-throat, psychiatry, or neurological disorder. This study was approved by Northwestern University's Institutional Review Board under Protocol #STU00201746. All participants gave their voluntary written consent before the experiment. All experiments were conducted according to the principles of the Declaration of Helsinki.
MRI data acquisition
Request a detailed protocolMagnetic resonance imaging (MRI) data were acquired on a 3T Siemens TIM Trio scanner equipped with a 64-channel head coil (Siemens Healthcare, Erlangen, Germany), at Northwestern University’s Center for Translational Imaging. A 10 min resting fMRI scan was acquired for each subject using a single-shot gradient-echo planar-imaging sequence with following parameters: repetition time (TR): 555 ms; echo time (TE): 22 ms; flip angle: 47°; MB-8 with Split-slice GRAPPA (Olman et al., 2009; Todd et al., 2016); field of view (FOV): 208 mm; voxel size: 2.0 × 2.0 × 2.0 mm3; 64 axial slices. To reduce the distortion and improve the signal-to-noise ratio in the primary olfactory and orbitofrontal areas, the slice orientation was set to approximately 30° from the AC-PC line (Deichmann et al., 2003). These acquisition parameters resulted in robust signals within orbitofrontal and olfactory areas in each subject (Figure 2—figure supplement 2). Further, our finding of strong connectivity between orbitofrontal cortex and all primary olfactory cortical areas suggests good orbitofrontal coverage. Subjects were instructed to look at a white fixation cross on a black background and to breathe in and out through their nose.
A high-resolution anatomical image was acquired for each subject using T1-weighted MPRAGE (TR: 2300 ms; TE: 2.94 ms; flip angle: 9°; FOV: 256 mm; voxel size: 1.0 × 1.0 × 1.0 mm3; 176 sagittal slices).
MRI data preprocessing
Request a detailed protocolThe structural images were skull-stripped and segmented into gray matter, white matter and cerebrospinal fluid using the BET (Smith, 2002) and FAST (Zhang et al., 2001) tools of FSL (FMRIB Software Library, www.fmrib.ox.ac.uk/fsl; RRID:SCR_002823) (Jenkinson et al., 2012; Smith et al., 2004; Woolrich et al., 2009). The resulting white matter and cerebrospinal fluid images were further eroded by one voxel (FSL's fslmaths).
Preprocessing of the resting fMRI data included removal of the first 10 volumes, motion correction and generating spatial registration matrices using FSL's FEAT. Each subject's functional images were normalized to their anatomical image using the brain-boundary registration method, and each individual anatomical image was registered to the Montreal Neurological Institute (MNI) standard brain (MNI152_T1_2mm_brain) using the non-linear registration method (12 degrees of freedom). Linear and quadratic trends were removed using Analysis of Functional NeuroImages (AFNI; RRID:SCR_005927) (Cox, 1996). Nuisance variables, including six head-movement parameters, and white matter and cerebrospinal signals, were regressed out using multiple linear regression methods (FSL's fsl_glm). Finally, the images were intensity normalized, band-pass filtered (0.008–0.01 Hz, AFNI's 3dFourier), registered to MNI space and spatially smoothed (Gaussian kernel, sigma = 3).
Functional connectivity-based parcellation
Request a detailed protocolTo perform the functional connectivity-based parcellation, we manually drew an ROI that included the anterior olfactory nucleus, olfactory tubercle, and piriform cortex on the MNI152 template brain, according to the Atlas of the Human Brain (Mai et al., 2015). The ROI was drawn onto coronal slices in the range of y = −3 to y = 15 on the MNI152_T1_1mm_brain, which was down-sampled to MNI152 2 mm space afterwards. For the comparison between parcellated clusters and anatomical subregions, we drew ROIs delineating the primary olfactory subregions prior to running the parcellation analysis. Notably, the human primary olfactory subregions differ in relative size and shape from their rodent analogs. We followed the definition of the olfactory tubercle and anterior olfactory nucleus used in the Atlas of the Human Brain (Mai et al., 2015; Ongür et al., 2003), and their locations agree with previous human studies (Allison, 1954; Crosby and Humphrey, 1941; Eslinger et al., 1982). Although we used an atlas that includes particularly detailed delineations of olfactory areas, based on rigorous techniques, including data from four human brains (Ongür et al., 2003), some researchers have postulated that the human anterior olfactory nucleus is part of the olfactory bulb (Daniel and Hawkes, 1992; Hyman et al., 1991) and that the human olfactory tubercle is actually the anterior perforated substance (Daniel and Hawkes, 1992; Hyman et al., 1991). However, these discrepencies are mainly terminology-based, as the different studies refer to the same anatomical structures.
We next parcellated the ROI into subregions based on their connectivity patterns with the rest of the brain. To do so, the Pearson correlation coefficient was computed between each voxel within the ROI and every other voxel in the rest of the whole brain, resulting a correlation matrix for each subject. The whole-brain mask was created using FSL's gray matter tissue prior image (avg152T1_gray.img, threshold of 100). The correlation coefficient was Fisher's z transformed and then averaged across subjects. The resulting matrix was transformed back into Pearson correlation coefficients, which were later used for parcellation analysis.
The parcellation analysis was performed using standard k-means methods, as implemented in the Matlab Statistics Toolbox (Matlab R2016b, The Mathworks Inc, Natick, MA USA; RRID:SCR_001622). The correlation between the connectivity pattern of the voxels within the ROI was used as the distance measure. Although this method is unsupervised, it requires an input of the number of clusters. Our hypothesis was that subregions of the primary olfactory areas within the ROI were separable based on their functional connectivity patterns. Thus, we chose a k value of four as the input for the algorithm. The clustering was performed on left and right hemisphere ROIs separately.
To evaluate the stability of the group-level connectivity patterns, we tested the stability of the correlation matrix using a leave-one-out method (Kahnt et al., 2012). This analysis tests whether the connectivity profile of individual voxels in primary olfactory cortex is similar across subjects, as required for averaging. For this, the correlation matrices were averaged across N−1 subjects, where N is the total number of subjects, and the correlation with the left-out subject was calculated for each voxel. We repeated this procedure N times and took the average of these repetitions as the final stability map.
To calculate the proportion of voxels from each parcellation cluster located within each anatomical subdivision of primary olfactory subregions, we first used the Atlas of the Human Brain (Mai et al., 2015 to outline anatomical ROIs of anterior olfactory nucleus, olfactory tubercle, frontal piriform cortex and temporal piriform cortex prior to performing the parcellation analysis. We then determined the proportion of voxels from each parcellated subdivision that were located within each ROI. The significance of the proportion was tested using a permutation method. In each permutation, we shuffled the labels of the voxels within the ROI and re-calculated the proportion number. This procedure was repeated 10,000 times, resulting in a null distribution of the proportion of voxels of each parcellated subdivision within each anatomical subregion. The mean and standard deviation of this distribution was computed by norm line fitting (Matlab's normfit). A z score of the real proportion was computed by subtracting the average and then dividing by the standard deviation.
To characterize the functional connectivity pattern of each subregion, the time series of all voxels were averaged for each ROI. The Pearson correlation between the average time series and every other voxel in the brain was calculated for each subject and Fisher's z transformed. Finally, random effects analysis of the functional connectivity maps was performed using one-sample t tests (FSL’s randomize, 10,000 permutations). Multiple comparisons were corrected using the TFCE method (Smith and Nichols, 2009). TFCE allows the identification of clusters in data sets without defining the clusters in a binary way. The output is a weighted sum of the local clustered signal. TCFE has been shown to produce results with better stability than other methods (Smith and Nichols, 2009). For functional connectivity network analysis, we focused on positive functional connectivity only, since the mechanisms of negative functional connectivity are less understood (Murphy et al., 2009).
To examine the functional connectivity difference between left and right subregions of the ROI, we computed a lateralization index of the functional connectivity. The lateralization index was calculated as (Zleft−Zright)/(Zleft + Zright), where Zleft and Zrightrepresent the whole-brain functional connectivity map of the left and right subregions respectively. The lateralization index was further spatially smoothed (Gaussian kernel, sigma = 3). A one-sample t test analysis of the lateralization index was conducted using FSL's randomize (10,000 permutations) and multiple comparisons were corrected using the TFCE method. Because we found no statistical difference between the maps across hemispheres, we combined the corresponding left and right primary olfactory cortex clusters, and performed all subsequent analyses on these hemispherically-combined seed regions.
For the replication of our parcellation analysis, we performed k-means clustering on subjects from an independent resting-state fMRI dataset that was published elsewhere (Kahnt and Tobler, 2017). In brief, 6 min of resting fMRI data were collected from fifty-three healthy subjects on a Philips Achieva 3T scanner (TR: 2000 ms; TE: 30 ms; voxel size: 2.75 × 2.75 × 3 mm3; flip angle: 90°). A high-resolution, T1-weighted MPRAGE anatomical image was acquired (TR: 8.2 ms; TE: 3.8 ms; FOV: 256 mm; voxel size: 1.0 × 1.0 × 1.0 mm3; 181 slices; flip angle: 8°) for each subject. All analyses performed on this data set were identical to the steps performed in our initial analysis. Full details on the acquisition parameters can be found in Kahnt and Tobler (2017), but importantly, their acquisition parameters differed from ours, highlighting the replicability of our findings.
Data availability
Request a detailed protocolSource data and code have been made available via GitHub: https://github.com/zelanolab/primaryolfactorycortexparcellation.git (Zelano Lab, 2019; copy archived at https://github.com/elifesciences-publications/primaryolfactorycortexparcellation).
Data availability
Source data have been provided for each figure and table. Source data and code have been made available via GitHub: https://github.com/zelanolab/primaryolfactorycortexparcellation.git (copy archived at https://github.com/elifesciences-publications/primaryolfactorycortexparcellation).
References
-
The contribution(s) of the insula to speech production: a review of the clinical and functional imaging literatureBrain Structure and Function 214:419–433.https://doi.org/10.1007/s00429-010-0257-x
-
The secondary olfactory areas in the human brainJournal of Anatomy 88:481–488.
-
Topographic organization of hippocampal inputs to the anterior olfactory nucleusFrontiers in Neuroanatomy 12:12:12.https://doi.org/10.3389/fnana.2018.00012
-
Participation of the insula in language revisited: a meta-analytic connectivity studyJournal of Neurolinguistics 29:31–41.https://doi.org/10.1016/j.jneuroling.2014.02.001
-
Structural and functional MRI differences in master sommeliers: a pilot study on expertise in the brainFrontiers in Human Neuroscience 10:414.https://doi.org/10.3389/fnhum.2016.00414
-
Hedonic-specific activity in piriform cortex during odor imagery mimics that during odor perceptionJournal of Neurophysiology 98:3254–3262.https://doi.org/10.1152/jn.00349.2007
-
A field guide to the anterior olfactory nucleus (cortex)Brain Research Reviews 50:305–335.https://doi.org/10.1016/j.brainresrev.2005.08.005
-
Central olfactory connections in the macaque monkeyThe Journal of Comparative Neurology 346:403–434.https://doi.org/10.1002/cne.903460306
-
Lateral entorhinal modulation of piriform cortical activity and fine odor discriminationJournal of Neuroscience 33:13449–13459.https://doi.org/10.1523/JNEUROSCI.1387-13.2013
-
Bidirectional plasticity of cortical pattern recognition and behavioral sensory acuityNature Neuroscience 15:155–161.https://doi.org/10.1038/nn.2966
-
BookCentral olfactory processingIn: Doty R. L, editors. Handbook of Olfaction and Gustation. New York: Marcel Dekker. pp. 165–180.https://doi.org/10.1201/9780203911457.ch8
-
Dynamic cortical lateralization during olfactory discrimination learningThe Journal of Physiology 593:1701–1714.https://doi.org/10.1113/jphysiol.2014.288381
-
AFNI: software for analysis and visualization of functional magnetic resonance neuroimagesComputers and Biomedical Research 29:162–173.https://doi.org/10.1006/cbmr.1996.0014
-
Studies of The vertebrate telencephalon. I. the nuclear configuration of the olfactory and accessory olfactory formations and of the nucleus olfactorius anterior of certain reptiles, birds, and mammalsThe Journal of Comparative Neurology 71:121–213.
-
Olfaction as a marker for depressionJournal of Neurology 264:631–638.https://doi.org/10.1007/s00415-016-8227-8
-
Functional subdomains within human FFAJournal of Neuroscience 33:16748–16766.https://doi.org/10.1523/JNEUROSCI.1259-13.2013
-
Olfactory identification deficits and increased mortality in the communityAnnals of Neurology 78:401–411.https://doi.org/10.1002/ana.24447
-
Imaging-based parcellations of the human brainNature Reviews Neuroscience 19:672–686.https://doi.org/10.1038/s41583-018-0071-7
-
Altered responses to social chemosignals in autism spectrum disorderNature Neuroscience 21:111–119.https://doi.org/10.1038/s41593-017-0024-x
-
Long-term memory of odors with and without verbal descriptionsJournal of Experimental Psychology 100:221–227.https://doi.org/10.1037/h0035492
-
Olfactory dysfunction in man: anatomical and behavioral aspectsBrain and Cognition 1:259–285.https://doi.org/10.1016/0278-2626(82)90028-8
-
Abnormal functional connectivity of frontopolar subregions in Treatment-Nonresponsive major depressive disorderBiological Psychiatry: Cognitive Neuroscience and Neuroimaging 3:337–347.https://doi.org/10.1016/j.bpsc.2017.12.003
-
Multidimensional representation of odors in the human olfactory cortexHuman Brain Mapping 37:2161–2172.https://doi.org/10.1002/hbm.23164
-
Neuroanatomical correlates of olfactory performanceExperimental Brain Research 201:1–11.https://doi.org/10.1007/s00221-009-1999-7
-
The olfactory tubercle encodes odor valence in behaving miceJournal of Neuroscience 35:4515–4527.https://doi.org/10.1523/JNEUROSCI.4750-14.2015
-
Unilateral lesions of the olfactory tubercle modifying general arousal effects in the rat olfactory bulbElectroencephalography and Clinical Neurophysiology 46:665–674.https://doi.org/10.1016/0013-4694(79)90104-4
-
Distinct cortical networks support the planning and online control of reaching-to-grasp in humansEuropean Journal of Neuroscience 35:909–915.https://doi.org/10.1111/j.1460-9568.2012.08018.x
-
MR volumetric analysis of the piriform cortex and cortical amygdala in drug-refractory temporal lobe epilepsyAJNR. American Journal of Neuroradiology 26:319–332.
-
Functional heterogeneity in human olfactory cortex: an event-related functional magnetic resonance imaging studyThe Journal of Neuroscience 22:10819–10828.https://doi.org/10.1523/JNEUROSCI.22-24-10819.2002
-
Central mechanisms of odour object perceptionNature Reviews Neuroscience 11:628–641.https://doi.org/10.1038/nrn2883
-
Flavor recognition memory related activity of the posterior piriform cortex in adult and aged ratsBehavioural Brain Research 360:196–201.https://doi.org/10.1016/j.bbr.2018.12.016
-
Association and commissural fiber systems of the olfactory cortex of the rat. I. systems originating in the piriform cortex and adjacent AreasThe Journal of Comparative Neurology 178:711–740.https://doi.org/10.1002/cne.901780408
-
Olfactory cortical neurons read out a relative time code in the olfactory bulbNature Neuroscience 16:949–957.https://doi.org/10.1038/nn.3407
-
Behavioral asymmetries following olfactory tubercle lesions in catsBrain, Behavior and Evolution 14:241–250.https://doi.org/10.1159/000125664
-
Effects of rostral medial forebrain bundle and olfactory tubercle lesions upon sexual behavior of male ratsJournal of Comparative and Physiological Psychology 82:30–36.https://doi.org/10.1037/h0033797
-
Odor quality coding and categorization in human posterior piriform cortexNature Neuroscience 12:932–938.https://doi.org/10.1038/nn.2324
-
Pathologic changes in the olfactory system in aging and Alzheimer's diseaseAnnals of the New York Academy of Sciences 640:14–19.https://doi.org/10.1111/j.1749-6632.1991.tb00184.x
-
Involvement of the olfactory tubercle in cocaine reward: intracranial self-administration studiesThe Journal of Neuroscience 23:9305–9311.https://doi.org/10.1523/JNEUROSCI.23-28-09305.2003
-
Projections from orbitofrontal cortex to anterior piriform cortex in the rat suggest a role in olfactory information processingThe Journal of Comparative Neurology 488:224–231.https://doi.org/10.1002/cne.20595
-
Odor-evoked activity is spatially distributed in piriform cortexThe Journal of Comparative Neurology 457:361–373.https://doi.org/10.1002/cne.10557
-
Connectivity-based parcellation of the human orbitofrontal cortexJournal of Neuroscience 32:6240–6250.https://doi.org/10.1523/JNEUROSCI.0257-12.2012
-
Dopamine modulates the functional organization of the orbitofrontal cortexThe Journal of Neuroscience 37:1493–1504.https://doi.org/10.1523/JNEUROSCI.2827-16.2016
-
The fusiform face area: a module in human extrastriate cortex specialized for face perceptionThe Journal of Neuroscience 17:4302–4311.https://doi.org/10.1523/JNEUROSCI.17-11-04302.1997
-
Rapidly acquired multisensory association in the olfactory cortexBrain and Behavior 5:e00390.https://doi.org/10.1002/brb3.390
-
Spatial distribution of neural activity in the anterior olfactory nucleus evoked by odor and electrical stimulationThe Journal of Comparative Neurology 519:277–289.https://doi.org/10.1002/cne.22519
-
Bidirectional processing in the olfactory-limbic Axis during olfactory behaviorBehavioral Neuroscience 112:541–553.https://doi.org/10.1037/0735-7044.112.3.541
-
Disrupted olfactory integration in schizophrenia: functional connectivity studyInternational Journal of Neuropsychopharmacology 20:740–746.https://doi.org/10.1093/ijnp/pyx045
-
Emotional control, reappraisedNeuroscience & Biobehavioral Reviews 95:528–534.https://doi.org/10.1016/j.neubiorev.2018.11.003
-
Frontal pole function: what is specifically human?Trends in Cognitive Sciences 15:241.https://doi.org/10.1016/j.tics.2011.04.005
-
When the sense of smell meets emotion: anxiety-state-dependent olfactory processing and neural circuitry adaptationThe Journal of Neuroscience 33:15324–15332.https://doi.org/10.1523/JNEUROSCI.1835-13.2013
-
Synaptic integration of olfactory information in mouse anterior olfactory nucleusJournal of Neuroscience 26:12023–12032.https://doi.org/10.1523/JNEUROSCI.2598-06.2006
-
Neuromodulation and the functional dynamics of piriform cortexChemical Senses 26:585–594.https://doi.org/10.1093/chemse/26.5.585
-
Learning-induced oscillatory activities correlated to odour recognition: a network activityEuropean Journal of Neuroscience 23:1801–1810.https://doi.org/10.1111/j.1460-9568.2006.04711.x
-
Face-specific processing in the human fusiform gyrusJournal of Cognitive Neuroscience 9:605–610.https://doi.org/10.1162/jocn.1997.9.5.605
-
The muted sense: neurocognitive limitations of olfactory languageTrends in Cognitive Sciences 19:314–321.https://doi.org/10.1016/j.tics.2015.04.007
-
Architectonic subdivision of the human orbital and medial prefrontal cortexThe Journal of Comparative Neurology 460:425–449.https://doi.org/10.1002/cne.10609
-
Risk factors for hazardous events in olfactory-impaired patientsJAMA Otolaryngology–Head & Neck Surgery 140:951–955.https://doi.org/10.1001/jamaoto.2014.1675
-
The anatomy of "paleocortex". A critical reviewErgebnisse Der Anatomie Und Entwicklungsgeschichte 43:3–62.https://doi.org/10.1007/978-3-642-65004-8
-
BookOlfactory higher centers anatomyIn: Squire L. arryR, editors. Encyclopedia Neuroscience. Academic Press. pp. 129–136.
-
Spatial and temporal distribution of odorant-evoked activity in the piriform cortexJournal of Neuroscience 27:1534–1542.https://doi.org/10.1523/JNEUROSCI.4072-06.2007
-
Lateralization of olfactory processesChemical Senses 29:731–745.https://doi.org/10.1093/chemse/bjh067
-
The organization of projections from the olfactory bulb to the piriform cortex and olfactory tubercle in the ratThe Journal of Comparative Neurology 194:519–534.https://doi.org/10.1002/cne.901940304
-
The dorsal posterior insula subserves a fundamental role in human painNature Neuroscience 18:499–500.https://doi.org/10.1038/nn.3969
-
BookOlfactory bulb.The Synaptic Organization of the BrainOxford University Press.https://doi.org/10.1093/acprof:oso/9780195159561.003.0005
-
Anatomy of the Human Olfactory Bulb and Central Olfactory PathwaysThe Human Sense of Smell, Anatomy of the Human Olfactory Bulb and Central Olfactory Pathways, 60, Berlin, Springer.
-
Fast robust automated brain extractionHuman Brain Mapping 17:143–155.https://doi.org/10.1002/hbm.10062
-
Time course of odorant-induced activation in the human primary olfactory cortexJournal of Neurophysiology 83:537–551.https://doi.org/10.1152/jn.2000.83.1.537
-
Structure of the nucleus olfactorius anterior of the hedgehog (Erinaceus europaeus)The Journal of Comparative Neurology 279:581–600.https://doi.org/10.1002/cne.902790407
-
The piriform cortex and human focal epilepsyFrontiers in Neurology 5:259.https://doi.org/10.3389/fneur.2014.00259
-
Whole brain resting-state analysis reveals decreased functional connectivity in major depressionFrontiers in Systems Neuroscience 4:.https://doi.org/10.3389/fnsys.2010.00041
-
Olfaction and face encoding in humans: a magnetoencephalographic studyCognitive Brain Research 15:105–115.https://doi.org/10.1016/S0926-6410(02)00144-1
-
Olfaction and its dynamic influence on word and face processing: cross-modal integrationProgress in Neurobiology 84:192–209.https://doi.org/10.1016/j.pneurobio.2007.10.005
-
Functional magnetic resonance imaging study of human olfaction and normal agingThe Journals of Gerontology Series A: Biological Sciences and Medical Sciences 60:510–514.https://doi.org/10.1093/gerona/60.4.510
-
Smelling sounds: olfactory-auditory sensory convergence in the olfactory tubercleJournal of Neuroscience 30:3013–3021.https://doi.org/10.1523/JNEUROSCI.6003-09.2010
-
Sniffing out the contributions of the olfactory tubercle to the sense of smell: hedonics, sensory integration, and more?Neuroscience & Biobehavioral Reviews 35:655–668.https://doi.org/10.1016/j.neubiorev.2010.08.004
-
BookOlfactory cortex physiologyIn: Squire L. R, editors. Encyclopedia of Neuroscience. Academic Press. pp. 95–100.https://doi.org/10.1016/B978-008045046-9.01690-9
-
Brain injury in autonomic, emotional, and cognitive regulatory areas in patients with heart failureJournal of Cardiac Failure 15:214–223.https://doi.org/10.1016/j.cardfail.2008.10.020
-
Olfactory influences on appetite and satiety in humansPhysiology & Behavior 87:800–804.https://doi.org/10.1016/j.physbeh.2006.01.029
-
Attentional modulation in human primary olfactory cortexNature Neuroscience 8:114–120.https://doi.org/10.1038/nn1368
-
Diverse patterns of odor representation by neurons in the anterior piriform cortex of awake miceJournal of Neuroscience 30:16662–16672.https://doi.org/10.1523/JNEUROSCI.4400-10.2010
-
Segmentation of brain MR images through a hidden markov random field model and the expectation-maximization algorithmIEEE Transactions on Medical Imaging 20:45–57.https://doi.org/10.1109/42.906424
Decision letter
-
Joel MainlandReviewing Editor; Monell Chemical Senses Center, United States
-
Joshua I GoldSenior Editor; University of Pennsylvania, United States
-
Jonas OlofssonReviewer
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Characterizing functional pathways of the human olfactory system" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a guest Reviewing Editor and Joshua Gold as the Senior Editor. The following individual involved in review of your submission has agreed to reveal his identity: Jonas Olofsson (Reviewer #1).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
The majority of our knowledge regarding primary olfactory cortex anatomy comes from rodents and the field needs more human data. This study uses brain imaging to add additional functional annotation to the human primary olfactory cortex.
Essential revisions:
Several reviewers asked for a more thorough review of the previous literature using fMRI to measure connectivity. Please revise the manuscript to place this study in context of previous work.
Additionally, the reviewers raised concerns and made a number of suggestions regarding technical aspects, lateralization, and addressing individual variability. Please see below for the individual points and address them as much as possible.
Reviewer #1:
The manuscript by Lane et al. describes an experiment whereby functional Magnetic Resonance Imaging connectivity analysis was used to parcellate the human primary olfactory cortex into separate connective pathways. Results suggest four distinct pathways emanating from primary areas (those receiving direct input from the olfactory bulb).
This manuscript is well written, especially considering the rodent neuroanatomical work which is covered in detail, and the results, which are clearly explained. Importantly, the study addresses a key issue in human olfaction, that of the presumed complementary roles of primary olfactory areas and how to define them. Prior work has used seed regions in the piriform, and sometimes even orbitofrontal, cortex to study the connectivity patterns of the olfactory system. The current study provides a more in-depth view, and will facilitate future comparisons within the olfactory network and across sensory systems. The independent validation dataset strengthens the authors' conclusions.
Overall, the study will provide a nice addition to the literature on human olfactory system connectivity, especially since the literature is currently very limited. I have a few issues for the authors to address, these issues are listed below.
The manuscript does not cite the small number of prior studies on fMRI connectivity of the olfactory cortex and it does not adequately describe how the present approach is different from that of prior work of the same kind.
Please motivate your definition of primary olfactory cortex, given the notion that the olfactory bulb is by some considered primary cortex.
Although the methods appear later in the text, there should be an early citation or mentioning of the dataset and its resting-state methodology; only the validation dataset is mentioned in the early stages of the main text.
Discuss the commonly observed problem of poor signal-to-noise ratio in olfactory regions and how the authors avoided the risk of biasing their results due to poor/variable SNR.
You found a mean correlation coefficient of 0.19 (SE 0.0041), how does that compare with prior applications of this method on similar data sets (e.g. Kahnt et al., 2012)? Please provide support for the claim that there is good similarity.
It is interesting to note that AON has a stronger lateralization pattern. Why do you think that is the case, and what are the possible implications? Does it imply a specific functional role or might it be an artifact, or evidence of less well-defined network memberships?
Please provide references for the olfactory interactions with the language regions, as is hypothesized in Discussion.
Regarding the framework, is might be noted that the shared connectivity with the OFC might explain why it is so reliably elicited even in basic/passive olfactory tasks; all primary regions project to it. Furthermore, regarding the connectivity patterns of the frontal piriform to regions responsible for facial movements, might this include sniffing?
It is a bit of a stretch to state that there is a "large body of evidence" suggesting olfaction is lateralized to the right side. The literature is in fact quite inconsistent. I suggest tempering this claim and also adding another more recent supporting citation.
Reviewer #2:
The conclusion appears somewhat weak: "… Our results provide important insight into the functional and anatomical organization of the human olfactory system…" Please be more detailed in terms of the significance of these findings and the consequences of it. Do we have to look at olfactory information processing differently from today on? Please speculate: Does this demonstration of parallel processing of olfactory information explain the vulnerability in aging? Or the decrease of olfactory functioning in depression? Other consequences?
Please compare the parallel processing n the olfactory system to other systems? How is that in taste? Touch?
I was always under the assumption that the AON in humans is part of the olfactory bulb (see Daniel and Hawkes, 1992; Hyman, Arriagada and Van Hoesen, 1991). Maybe you can discuss these discrepancies, too.
The authors talk about the olfactory tubercle. This is despite the anatomical impression that there is no bump in this vicinity. In fact the anterior perforated substance has been proposed as its human homologue (see Freiherr, Wiesmann, Witt in Welge-Luessen, Hummel (eds): 2014, Management of Smell and Taste Disorders; Theme, Stuttgart). Please comment.
It would be extremely valuable for the reader if the authors came up with a sketch that graphically summarizes the results in terms of the strength of the connectivity between the various regions.
Although difficult to do (and although the authors already pay a lot of attention to that), it would be excellent if the authors could try to separate more clearly between research that comes from animals and research on humans.
Reviewer #3:
Odor processing starts with transduction in the olfactory epithelium and continues with treatment in the olfactory bulb and cortical processing in primary and so-called secondary areas. In mammals, we have very good understanding of the anatomy and the functionality of these areas. In contrast, the functional neuroanatomy of human olfaction remains poorly understood. In this paper, Zhou and colleagues provide a significant contribution to this field by highlighting that the human primary olfactory cortex is clustered into sub-regions that anatomically corresponded to the anterior olfactory nucleus, olfactory tubercle and frontal and temporal piriform cortices. Furthermore, a dissociation in the whole-brain functional connectivity patterns was observed across the different sub-areas of the primary olfactory cortex suggesting the existence of distinct parallel processing pathways. The paper is very well written and the method used are innovative in the field of olfaction (data-driven connectivity-based parcellation techniques). My comments are below.
- Subsection “Parcellation results across hemispheres and k values”: Why it is important to combine data from both hemispheres? Asymmetry in human olfactory processing has been documented, especially by the Zatorre group in the early 90's and more recently in different studies. I would like to see how not merging the 2 hemispheres change the findings.
- The authors mentioned “Note that additional primary olfactory areas, including anterior cortical and peri- amygdala areas and entorhinal cortex (Allison, 1954; Carmichael et al., 1994; Eslinger et al., 1982; Gonçalves Pereira et al., 2005; Zatorre et al., 1992), were not included in our ROI because the exact location of olfactory afferents into these areas is poorly understood.”. I am wondering whether some areas included in the ROI (tubercule, AON), are less “poorly understood” than entorhinal cortex or amygdala in humans. Could the authors provide some more rationales for excluding these areas from the ROI?
In the same line, the authors wrote: “Interestingly, connectivity clusters were generally more extensive in the right hemisphere compared to the left hemisphere, supporting a large body of evidence suggesting that olfactory processing is lateralized to the right side (Zatorre et al., 1992).” So, why not considering the primary olfactory cortices separately from the right and the left hemispheres? There is much to gain in showing these data.
- The authors stated that “We then estimated the whole-brain functional connectivity profile of each voxel within the ROI by computing the Pearson correlation coefficient between the resting-state fMRI time-series of a given voxel and that of every other voxel in the rest of the brain. This resulted in subject-wise connectivity matrices, which were then averaged across all subjects.” Does this mean that individual variability is not considered in the rest of the analysis? If yes, how the authors could be sure that the values in the averaged matrix are not due to a sub-sample of subjects?
- In the same line, the authors wrote “To examine the similarity between the individual functional connectivity matrices, we computed a histogram of correlation values between individual matrixes and the groups matrix using a leave-one-out method (Kahnt et al., 2012). The results indicate a good similarity between the group-level and individual functional connectivity patterns (mean correlation coefficient: 0.19, standard error: 0.0041, Figure 2A).”. I can understand that the notion of similarity is relative, but here the averaged r-value is less than 0.20 which means that half of the data fall within a range of 0 to 0.20. I may not understand properly this analysis but to my view it seems that these r-values are low. Could you comment on that? Moreover, the authors must mention the significativity of these correlative analyses. Even if significant, there is a need to detail how many data-points were computed in these analyses and how such low r-values can be interpreted?
- Participants were “instructed to look at a white fixation cross on a black background and breath in-out through their nose.” As the authors know, and they are expert in the field, respiration can alter odor processing. Is there any effect of Nasal breathing on the brain imaging findings? If nasal respiration was measured, in addition to motion correction, the authors should correct for nasal breathing activity.
https://doi.org/10.7554/eLife.47177.024Author response
Reviewer #1:
[…] Overall, the study will provide a nice addition to the literature on human olfactory system connectivity, especially since the literature is currently very limited. I have a few issues for the authors to address, these issues are listed below.
The manuscript does not cite the small number of prior studies on fMRI connectivity of the olfactory cortex and it does not adequately describe how the present approach is different from that of prior work of the same kind.
This is a very important point, and we are grateful for being made aware of having missed important previous work. We have modified the Introduction to include a discussion of the following studies. We would be grateful if you could bring our attention to any additional studies we might have missed.
We now include a discussion of the following:
Banks et al., 2016; Caffo et al., 2010; Cecchetto et al., 2019; Fjaeldstad et al., 2017; Karunanayaka et al., 2014.; Karunanayaka, Tobia and Yang, 2017; Killgore et al., 2013; Kiparizoska and Ikuta2017; Kollndorfer et al., 2015; Krusemark and Li, 2012; Nigri et al., 2013; Sreenivasan et al., 2017; Sunwoo et al., 2015; Wang et al.,2010; Wang et al., 2005.; Wang et al., 2015.
Please motivate your definition of primary olfactory cortex, given the notion that the olfactory bulb is by some considered primary cortex.
We agree with you that the olfactory bulb is considered by some to be “primary olfactory cortex”. However, we feel that the “classic” and dominant terminology in the field still uses the term “primary olfactory cortex” to refer to structures that receive direct input from the olfactory bulb. For example, in Olfactory Cortex: Comparative Anatomy, Illig and Wilson discuss the concept that the olfactory bulb may perform primary sensory cortical functions, but still mention that piriform cortex is typically referred to as primary olfactory cortex. Most book chapters, reviews and studies still refer to primary olfactory cortex as brain areas that receive direct bulb projections. We feel it would be confusing to introduce new terminology in this manuscript, and might take away from the primary goal of our study, which was to provide a quantitative characterization of the parallel processing pathways of the human olfactory system. In humans, these parallel pathways occur downstream of the olfactory bulb. Thus, the important point here is not so much whether we are characterizing primary versus secondary olfactory cortex, but that we are characterizing the parallel pathways, all of which receive direct projections from the olfactory bulb. That said, your point is highly relevant for this paper, and reflects an important emerging movement in the field towards the idea that the classic terminology is flawed. Further, our initial phrasing of this in the Introduction may have mis-characterized the Wilson, 2009 reference. We have clarified these points in the second paragraph of the Introduction. We now initially refer to these areas as the cortical targets of olfactory bulb projections, and we refer readers to papers that discuss the accuracy of their being collectively termed “primary olfactory cortex”.
Although the methods appear later in the text, there should be an early citation or mentioning of the dataset and its resting-state methodology; only the validation dataset is mentioned in the early stages of the main text.
Very good point, thank you. This will make the paper read more clearly. We have added details about the dataset into the beginning of the Results section: “Twenty-five subjects (average ± standard error age: 25.5 ± 1.2 years; 14 female) underwent a 10-minute resting-state fMRI scan.”
Discuss the commonly observed problem of poor signal-to-noise ratio in olfactory regions and how the authors avoided the risk of biasing their results due to poor/variable SNR.
This is a helpful suggestion, and we have added a discussion of this to the Materials and methods as below. We have also added a supplementary figure showing the mean activation map for each subject, in order to show the data transparently.
“To reduce the distortion and improve the signal-to-noise ratio in the primary olfactory and orbitofrontal areas, the slice orientation was set to approximately 30° from AC-PC line (Deichmann et al., 2003). These acquisition parameters resulted in robust signals within orbitofrontal and olfactory areas in each subject (Figure 2—figure supplement 2). Further, our finding of strong connectivity between orbitofrontal cortex and all primary olfactory cortical areas suggests good orbitofrontal coverage.”
You found a mean correlation coefficient of 0.19 (SE 0.0041), how does that compare with prior applications of this method on similar data sets (e.g. Kahnt et al. 2012)? Please provide support for the claim that there is good similarity.
This is an important question, thank you for raising it. One of the other reviewers also asked about this R-value, further confirming that it’s important for us to clarify this analysis. The purpose of this leave-one-out analysis is only to estimate how similar (or stable) connectivity profiles of individual primary olfactory cortex voxels are across participants, as a prerequisite for averaging the connectivity matrices across participants. The overlap observed here was slightly lower than what Kahnt et al., 2012, observed in orbitofrontal cortex (r~0.3) using a different data set. However, this value cannot be reasonably compared across different studies because it depends on many study-specific factors including voxel size, the size of the regions of interest, and the number of data points used to compute the correlations. The important point to note about this histogram is that the similarity of connectivity patterns is above zero in all voxels in primary olfactory cortex. We have clarified this point in the second paragraph of the subsection “Parcellation of human primary olfactory cortex”.
It is interesting to note that AON has a stronger lateralization pattern. Why do you think that is the case, and what are the possible implications? Does it imply a specific functional role or might it be an artifact, or evidence of less well-defined network memberships?
We fully agree, this is a very interesting result and we should provide additional comments on this. We do not believe the finding is an artifact, and we now include, in the Results section, some discussion of this finding that ties it in with existing literature.
There are some interesting rodent findings that support the idea that the functional networks of the AON are more lateralized than other primary olfactory areas. For example, it has also been shown that the AON neurons can distinguish between signals from the ipsilateral and contralateral nostrils, suggesting representation of lateralized inputs in this region.
Please provide references for the olfactory interactions with the language regions, as is hypothesized in Discussion.
Thank you for this point. This comment made us realize that we also left out implications of our data on the neurocognitive limitations of olfactory language. We’ve added a discussion of this, along with references, in the Discussion.
Regarding the framework, is might be noted that the shared connectivity with the OFC might explain why it is so reliably elicited even in basic/passive olfactory tasks; all primary regions project to it. Furthermore, regarding the connectivity patterns of the frontal piriform to regions responsible for facial movements, might this include sniffing?
These are good and insightful points. We have added them to the Discussion.
It is a bit of a stretch to state that there is a "large body of evidence" suggesting olfaction is lateralized to the right side. The literature is in fact quite inconsistent. I suggest tempering this claim and also adding another more recent supporting citation.
We agree. We have toned down this language by removing
“supporting a large body of evidence suggesting that olfactory processing is lateralized to the right side”. We have also added more recent supporting citations.
Reviewer #2:
The conclusion appears somewhat weak: "…Our results provide important insight into the functional and anatomical organization of the human olfactory system..:" Please be more detailed in terms of the significance of these findings and the consequences of it. Do we have to look at olfactory information processing differently from today on?
Thank you for this suggestion, we agree. We have modified the Introduction and Discussion to better distinguish our findings from previous studies, and to clarify the reasons why we believe the study is important. Your suggestion really improved the Introduction, thank you. See the new additions in the fifth paragraph of the Introduction, and in the first paragraph of the Discussion.
Please speculate: Does this demonstration of parallel processing of olfactory information explain the vulnerability in aging? Or the decrease of olfactory functioning in depression? Other consequences?
Great points, thank you. We have added a discussion of these ideas to the Discussion, which has improved this part of the manuscript substantially:
“Our findings provide a detailed description of the particular brain areas which exhibit unique connectivity with each individual primary olfactory subregion. […] For example, olfactory structures that form networks with brain areas that are implicated in particular pathologies likely perform critical olfactory sensory functions (the olfactory tubercle is connected to areas involved in depression, see Croy and Hummel, 2017, and the temporal piriform cortex is connected to areas involved in primary progressive aphasia, see Olofsson et al., 2013).”
Please compare the parallel processing n the olfactory system to other systems? How is that in taste? Touch?
This is a good suggestion. We now highlight that the organization of the olfactory system differs from other sensory modalities. We attempted a detailed review of the organization of other sensory systems, however, in the end we felt that it disrupted the flow of the manuscript and we removed it. If the reviewer feels strongly that it should be included, we will add it back in.
I was always under the assumption that the AON in humans is part of the olfactory bulb (see Daniel and Hawkes, 1992; Hyman, Arriagada and Van Hoesen, 1991). Maybe you can discuss these discrepancies, too.
This is an important point, which we should have discussed more in the manuscript. The human AON certainly seems to differ anatomically from the rodent AON. In humans, it is often referred to as the “retrobulbar region”, but it has been demonstrated to be a true cortical structure, equivalent to paleocortical structures in rodents (Stephan H: Allocortex. In Bargmann W, editor: “Handbuch der mikroskopischen Anatomie des Menschen”, vol 4, Berlin and New York. Part 9, 1975, Springer Verlag, pp 1–998). In fact, the term anterior olfactory nucleus has been suggested to be misleading, given that it is a true cortical structure in humans (Zilles and Amunts, Architecture of the Cerebral Cortex, 3rd Edition, 2012). A recent study by Ongur et al. performed a detailed mapping of the human AON based on data from 5 post-mortem human brains (Ongur, 2003), which was updated in their Atlas of the Human Brain (Mai et al., 2016). Their study showed that the human AON is extensive, extending more than 11mm in the anterior-posterior direction, from MNI position y=14.22mm to MNI position y=2.53mm, as detailed in the 4th Edition of the Atlas of the Human Brain (Mai et al., 2016). Numerous studies (most of which focus on Parkinson’s Disease pathology) have characterized the human anterior olfactory nucleus into several subdivisions other than the bulbar part. Just as one example, a study by Ubeda-Banon et al. in 2010 divided the human AON into the following different parts: bulbar, intrapeduncular, retrobulbar, cortical anterior and cortical posterior subregions. The two cortical subregions could be further divided into medial and lateral components on either side of the olfactory tract. See Author response image 1, taken from this study, which shows the size and extent of the human AON. Some other studies including subdivisions of the AON beyond the olfactory bulb include:
- Pearce, Hawkes and Daniel, 1995, Movement Disorders;
- Price JL (1990) Olfactory system. In: Paxinos G (ed) The human nervous system. Academic Press, San Diego, pp 979–998;
- Hoogland PV, Huisman E (1999) Tyrosine hydroxylase immunoreactive structures in the aged human olfactory bulb and olfactory peduncle. J Chem Neuroanat 17:153–161; -Mai JK, Paxinos G, Voss T (2008) Atlas of the human brain. Elsevier, New York;
- Price JL (1990) Olfactory system. In: Paxinos G (ed) The human nervous system. Academic Press, San Diego, pp 979–998;
- Braak, H., Del Tredici, K., Rüb, U., de Vos, R. A., Jansen Steur, E. N., and Braak, E. (2003a). Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol. Aging 24, 197–211;
- Sengoku R, Saito Y, Ikemura M, Hatsuta H, Sakiyama Y, Kanemaru K, Arai T, Sawabe M, Tanaka N, Mochizuki H, Inoue K, Murayama S (2008) Incidence and extent of Lewy body-related α-synucleinopathy in aging human olfactory bulb. J Neuropathol Exp Neurol 67:1072–1083;
- Ubeda-Banon et al., 2017.
In our study, we did not distinguish the subdivisions of the human AON, though this could be of high interest for future studies.
We now include a brief description of these points in the Materials and methods subsection “Functional connectivity-based parcellation”.
The authors talk about the olfactory tubercle. This is despite the anatomical impression that there is no bump in this vicinity. In fact the anterior perforated substance has been proposed as its human homologue (see Freiherr, Wiesmann, Witt in Welge-Luessen, Hummel (eds): 2014, Management of Smell and Taste Disorders; Theme, Stuttgart). Please comment.
We agree that the anatomy of the rodent and human olfactory tubercle differs. Similar to the AON, the human tubercle differs in relative size and shape from the rodent homologue. We followed the definition of the olfactory tubercle used in the Atlas of the Human Brain (Mai et al., 2015; Öngür et al., 2003), as mentioned in our previous response above. The location of the OT in this atlas agrees with previous studies that we cite in the paper (Allison, 1954; Crosby and Humphrey, 1941; Eslinger et al., 1982). It does appear that the currently accepted location of the human olfactory tubercle overlaps with the location of the anterior perforated substance. We have added mention of this to our Materials and methods, and cite the book mentioned.
It would be extremely valuable for the reader if the authors came up with a sketch that graphically summarizes the results in terms of the strength of the connectivity between the various regions.
This is a great point. Our connectivity maps did not contain information about the relative strength of the different connections. We have modified our figures to include t-values rather than simple binarized values. These new maps, that contain information about the strength of the connectivity, have been inserted to replace Figure 4.
Although difficult to do (and although the authors already pay a lot of attention to that), it would be excellent if the authors could try to separate more clearly between research that comes from animals and research on humans.
We agree this is an important point, and a common issue in human olfactory studies/manuscripts, since the vast majority of studies are conducted in rodents. We have made efforts wherever possible to make this distinction more clear.
Reviewer #3:
[…] - Subsection “Parcellation results across hemispheres and k values”: Why it is important to combine data from both hemispheres? Asymmetry in human olfactory processing has been documented, especially by the Zatorre group in the early 90's and more recently in different studies. I would like to see how not merging the 2 hemispheres change the findings.
This is an excellent point. We will clarify our thinking here and in the text.
Not merging the 2 hemispheres does not change the findings in a statistically relevant way.
We merged hemispheres for three main reasons. First, there are no significant differences between hemispheres, based on a statistical comparison of the connectivity maps. Second, given this lack of statistical difference across hemispheres, to show each hemisphere for each subregion would add additional information without additional meaning, to the manuscript, while making it denser and more difficult to interpret for the reader. Third, and most important, to show such data would risk misleading the reader and could lead to misinterpretation of our results.
Because of these concerns, and particularly the concern that readers could mis-cite our paper as showing differences across hemispheres when we found none statistically, we feel strongly that we should not add supplementary figures showing the single hemisphere maps. However, if the reviewer feels very strongly that we should add this figure, we will create a supplementary figure showing these maps.
Based on this helpful comment, we have added clarification of this to the Materials and methods subsection “Functional connectivity-based parcellation”.
- The authors mentioned “Note that additional primary olfactory areas, including anterior cortical and peri- amygdala areas and entorhinal cortex (Allison, 1954; Carmichael et al., 1994; Eslinger et al., 1982; Gonçalves Pereira et al., 2005; Zatorre et al., 1992), were not included in our ROI because the exact location of olfactory afferents into these areas is poorly understood.”. I am wondering whether some areas included in the ROI (tubercule, AON), are less “poorly understood” than entorhinal cortex or amygdala in humans. Could the authors provide some more rationales for excluding these areas from the ROI?
Thank you for this insightful comment. It is a good point, and we did not adequately explain our rationale in the manuscript.
The AON and OT are better understood anatomically (though we agree, not functionally) than the olfactory inputs to the amygdala and entorhinal cortex. For example, a recent histological study by Goncalves et al. found that the borders of the subdivisions of the amygdala and entorhinal cortex that receive direct bulb projections are not discernable in T1 images (Gonçalves Pereira et al., 2005). Furthermore, we were unable to find a human brain atlas that clearly outlines olfactory amygdala and entorhinal areas receiving bulb input, whereas numerous studies and atlases do include delineation of the olfactory tubercle and AON. We have modified the manuscript to include this reference where we discuss this point in the first paragraph of the Results subsection “Parcellation of human primary olfactory cortex”.
In the same line, the authors wrote: “Interestingly, connectivity clusters were generally more extensive in the right hemisphere compared to the left hemisphere, supporting a large body of evidence suggesting that olfactory processing is lateralized to the right side (Zatorre et al., 1992).” So, why not considering the primary olfactory cortices separately from the right and the left hemispheres? There is much to gain in showing these data.
Thank you for this comment. Interestingly, we found that across the whole-brain, connectivity was stronger between all (both left and right) primary olfactory areas and clusters in the right hemisphere. Importantly, this whole-brain hemispheric effect did not differ between primary olfactory cortical hemispheres—that is, both the right and left primary olfactory cortices exhibited increased right-sided whole-brain connectivity. This is a subtle distinction, and is shown in the figure displaying the common connectivity pattern (Figure 5). We clarified this in the subsection “Functional connectivity common to all subregions”.
- The authors stated that “We then estimated the whole-brain functional connectivity profile of each voxel within the ROI by computing the Pearson correlation coefficient between the resting-state fMRI time-series of a given voxel and that of every other voxel in the rest of the brain. This resulted in subject-wise connectivity matrices, which were then averaged across all subjects.” Does this mean that individual variability is not considered in the rest of the analysis? If yes, how the authors could be sure that the values in the averaged matrix are not due to a sub-sample of subjects?
This is a very important point, thank you for bringing this up. We computed the connectivity maps on an individual basis, but averaged the resulting connectivity matrices to obtain a group-connectivity map, as done previously (Kahnt et al., 2012; Kahnt and Tobler, 2017). For this, we first determined whether the connectivity structure of individual primary olfactory cortex voxels was indeed comparable across subjects. This was done using a leave-one-out analysis, which showed that there was substantial overlap, as required for averaging across subjects. We have clarified this in the second paragraph of the subsection “Parcellation of human primary olfactory cortex”, where we mention averaging across subjects.
- In the same line, the authors wrote “To examine the similarity between the individual functional connectivity matrices, we computed a histogram of correlation values between individual matrixes and the groups matrix using a leave-one-out method (Kahnt et al., 2012). The results indicate a good similarity between the group-level and individual functional connectivity patterns (mean correlation coefficient: 0.19, standard error: 0.0041, Figure 2A).”. I can understand that the notion of similarity is relative, but here the averaged r-value is less than 0.20 which means that half of the data fall within a range of 0 to 0.20. I may not understand properly this analysis but to my view it seems that these r-values are low. Could you comment on that? Moreover, the authors must mention the significativity of these correlative analyses. Even if significant, there is a need to detail how many data-points were computed in these analyses and how such low r-values can be interpreted?
These are very important points for our study, thank you for these comments. The purpose of this leave-one-out analysis was solely to determine whether we could reasonably average across the connectivity maps of individual subjects. Similar to previous studies (e.g., Kahnt et al., 2012) we expect this R value to be relatively small. Moreover, this value will likely differ across studies because it depends on study-specific differences in ROI sizes and voxel size. However, as can be seen in the histogram, most of these correlations were above zero, and because they are computed across rest-of-the-brain voxels (N=170,660), also statistically significant (R values larger than 0.0088 are significant at p<0.05, Bonferroni corrected for the number of voxels in the primary olfactory cortex). We have clarified these points in the second paragraph of the subsection “Parcellation of human primary olfactory cortex”.
- Participants were “instructed to look at a white fixation cross on a black background and breath in-out through their nose.” As the authors know, and they are expert in the field, respiration can alter odor processing. Is there any effect of Nasal breathing on the brain imaging findings? If nasal respiration was measured, in addition to motion correction, the authors should correct for nasal breathing activity.
Excellent point. In this study, we focused on low-frequency (< 0.1 Hz) oscillations, which are below the human respiratory frequency (typically from 0.15 – 0.3 Hz). Although there’s increasing interest in studying higher frequency brain activity using fMRI, the resting functional connectivity profile of high frequency oscillations and their relationship with respiration is not well understood. Our previous study has shown that nasal breathing can modulate local brain activity in a way that is different from oral breathing, and it would be interesting to examine whether functional networks are distinct during those two different breathing routines in future studies.
https://doi.org/10.7554/eLife.47177.025Article and author information
Author details
Funding
National Institute on Deafness and Other Communication Disorders (R01-DC-015426)
- Thorsten Kahnt
National Institute on Deafness and Other Communication Disorders (R00-DC-0012803)
- Christina Zelano
National Institute on Deafness and Other Communication Disorders (R01-DC-016364)
- Christina Zelano
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This study was financially supported by the National Institutes on Deafness and Other Communications Disorders (NIDCD) grants R00-DC-012803 and R01-DC-016364 (to Christina Zelano) and R01-DC-015426 (to Thorsten Kahnt).
Ethics
Human subjects: This study was approved by Northwestern University's Institutional Review Board. All participants gave their written consent before the experiment. The experiment was conducted according to the principles of the Declaration of Helsinki. Protocol #STU00201746.
Senior Editor
- Joshua I Gold, University of Pennsylvania, United States
Reviewing Editor
- Joel Mainland, Monell Chemical Senses Center, United States
Reviewer
- Jonas Olofsson
Version history
- Received: March 27, 2019
- Accepted: July 9, 2019
- Version of Record published: July 24, 2019 (version 1)
Copyright
© 2019, Zhou et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 10,840
- Page views
-
- 1,034
- Downloads
-
- 87
- Citations
Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
Diabetes is caused by the inability of electrically coupled, functionally heterogeneous -cells within the pancreatic islet to provide adequate insulin secretion. Functional networks have been used to represent synchronized oscillatory [Ca2+] dynamics and to study -cell subpopulations, which play an important role in driving islet function. The mechanism by which highly synchronized -cell subpopulations drive islet function is unclear. We used experimental and computational techniques to investigate the relationship between functional networks, structural (gap-junction) networks, and intrinsic -cell dynamics in slow and fast oscillating islets. Highly synchronized subpopulations in the functional network were differentiated by intrinsic dynamics, including metabolic activity and KATP channel conductance, more than structural coupling. Consistent with this, intrinsic dynamics were more predictive of high synchronization in the islet functional network as compared to high levels of structural coupling. Finally, dysfunction of gap junctions, which can occur in diabetes, caused decreases in the efficiency and clustering of the functional network. These results indicate that intrinsic dynamics rather than structure drive connections in the functional network and highly synchronized subpopulations, but gap junctions are still essential for overall network efficiency. These findings deepen our interpretation of functional networks and the formation of functional sub-populations in dynamic tissues such as the islet.
-
- Computational and Systems Biology
- Immunology and Inflammation
T cells are required to clear infection, and T cell motion plays a role in how quickly a T cell finds its target, from initial naive T cell activation by a dendritic cell to interaction with target cells in infected tissue. To better understand how different tissue environments affect T cell motility, we compared multiple features of T cell motion including speed, persistence, turning angle, directionality, and confinement of T cells moving in multiple murine tissues using microscopy. We quantitatively analyzed naive T cell motility within the lymph node and compared motility parameters with activated CD8 T cells moving within the villi of small intestine and lung under different activation conditions. Our motility analysis found that while the speeds and the overall displacement of T cells vary within all tissues analyzed, T cells in all tissues tended to persist at the same speed. Interestingly, we found that T cells in the lung show a marked population of T cells turning at close to 180o, while T cells in lymph nodes and villi do not exhibit this “reversing” movement. T cells in the lung also showed significantly decreased meandering ratios and increased confinement compared to T cells in lymph nodes and villi. These differences in motility patterns led to a decrease in the total volume scanned by T cells in lung compared to T cells in lymph node and villi. These results suggest that the tissue environment in which T cells move can impact the type of motility and ultimately, the efficiency of T cell search for target cells within specialized tissues such as the lung.