Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
    • Special Collections
  • EDITORIAL BOARD
    • Editorial Board
    • ECR Advisory Board
    • Journal Staff
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
    • Accessibility
  • SUBSCRIBE
PreviousNext
Research Articles, Behavioral/Cognitive

A Neural Mechanism in the Human Orbitofrontal Cortex for Preferring High-Fat Foods Based on Oral Texture

Putu A. Khorisantono, Fei-Yang Huang (黃飛揚), Michael P. F. Sutcliffe, Paul C. Fletcher, I. Sadaf Farooqi and Fabian Grabenhorst
Journal of Neuroscience 22 November 2023, 43 (47) 8000-8017; https://doi.org/10.1523/JNEUROSCI.1473-23.2023
Putu A. Khorisantono
1Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge CB2 3DY, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Putu A. Khorisantono
Fei-Yang Huang (黃飛揚)
1Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge CB2 3DY, United Kingdom
2Department of Experimental Psychology, University of Oxford, Oxford OX1 3TA, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael P. F. Sutcliffe
3Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul C. Fletcher
4Wellcome Trust-MRC Institute of Metabolic Science, Addenbrooke's Hospital, Cambridge CB2 0QQ, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
I. Sadaf Farooqi
4Wellcome Trust-MRC Institute of Metabolic Science, Addenbrooke's Hospital, Cambridge CB2 0QQ, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fabian Grabenhorst
1Department of Physiology, Development and Neuroscience, University of Cambridge, Cambridge CB2 3DY, United Kingdom
2Department of Experimental Psychology, University of Oxford, Oxford OX1 3TA, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Although overconsumption of high-fat foods is a major driver of weight gain, the neural mechanisms that link the oral sensory properties of dietary fat to reward valuation and eating behavior remain unclear. Here we combine novel food-engineering approaches with functional neuroimaging to show that the human orbitofrontal cortex (OFC) translates oral sensations evoked by high-fat foods into subjective economic valuations that guide eating behavior. Male and female volunteers sampled and evaluated nutrient-controlled liquid foods that varied in fat and sugar (“milkshakes”). During oral food processing, OFC activity encoded a specific oral-sensory parameter that mediated the influence of the foods' fat content on reward value: the coefficient of sliding friction. Specifically, OFC responses to foods in the mouth reflected the smooth, oily texture (i.e., mouthfeel) produced by fatty liquids on oral surfaces. Distinct activity patterns in OFC encoded the economic values associated with particular foods, which reflected the subjective integration of sliding friction with other food properties (sugar, fat, viscosity). Critically, neural sensitivity of OFC to oral texture predicted individuals' fat preferences in a naturalistic eating test: individuals whose OFC was more sensitive to fat-related oral texture consumed more fat during ad libitum eating. Our findings suggest that reward systems of the human brain sense dietary fat from oral sliding friction, a mechanical food parameter that likely governs our daily eating experiences by mediating interactions between foods and oral surfaces. These findings identify a specific role for the human OFC in evaluating oral food textures to mediate preference for high-fat foods.

SIGNIFICANCE STATEMENT Fat and sugar enhance the reward value of food by imparting a sweet taste and rich mouthfeel but also contribute to overeating and obesity. Here we used a novel food-engineering approach to realistically quantify the physical-mechanical properties of high-fat liquid foods on oral surfaces and used functional neuroimaging while volunteers sampled these foods and placed monetary bids to consume them. We found that a specific area of the brain's reward system, the orbitofrontal cortex, detects the smooth texture of fatty foods in the mouth and links these sensory inputs to economic valuations that guide eating behavior. These findings can inform the design of low-calorie fat-replacement foods that mimic the impact of dietary fat on oral surfaces and neural reward systems.

  • dietary fat
  • oral food processing
  • neuroeconomics
  • preference
  • reward value

Introduction

Oral food processing involves a chain of sensory and neural events in which a food's physical structure elicits oral sensations and subjective valuations that guide eating behavior. For example, foods high in sugar and fat typically produce a characteristic “mouthfeel” in the form of a sweet taste and thick, smooth texture. These sensory signals are thought to contribute to the reward value of food, to near-universal preferences for foods high in sugar and fat, and to the development of obesity (Alonso-Alonso et al., 2015; Carreiro et al., 2016; Rolls, 2020; Huang et al., 2021).

Specific brain areas, such as the orbitofrontal cortex (OFC), play key roles in assigning reward values to foods based on sensory inputs (Grabenhorst and Rolls, 2011; Rangel, 2013; Rolls, 2016; Morys et al., 2022). Studies in decision neuroscience identified neural responses to visual food stimuli that reflected valuation and related decision-making based on food attributes, such as nutrients, healthiness, and tastiness, specifically in the OFC and ventromedial PFC (Hare et al., 2009, 2011; Grabenhorst et al., 2013; Tang et al., 2014; Suzuki et al., 2017). Investigations in sensory neuroscience identified single-neuron activities in macaques and neuroimaging signals in humans that reflect oral food properties, including taste, odor, flavor, oral temperature, and texture, specifically in the OFC, insula, amygdala, and somatosensory cortex (de Araujo and Rolls, 2004; Small et al., 2004, 2008; Guest et al., 2007; Grabenhorst et al., 2008, 2010; Eldeghaidy et al., 2011; Grabenhorst and Rolls, 2014; Avery et al., 2020; Rolls, 2020). Despite these advances, the neural mechanisms that link the physical, material properties of foods to oral sensations, reward valuations, and economic preferences remain poorly understood.

Although sugar is sensed directly by taste receptors (Yarmolinsky et al., 2009), the mechanism for oral-fat sensing remains a topic of debate (Mattes, 2009). Substantial evidence points to a somatosensory, oral-texture mechanism for fat sensing (Kokini, 1987; Rolls, 2016). Fat-like textures elicit a smooth, creamy mouthfeel (Chojnicka-Paszun et al., 2012; Upadhyay et al., 2020) and activate neural sensory and reward systems (de Araujo and Rolls, 2004; Grabenhorst et al., 2010; Grabenhorst and Rolls, 2014; Rolls et al., 2018). Traditionally, viscosity has been considered a key oral-texture parameter important in fat sensing (Kokini, 1987; Mela, 1988; Stokes et al., 2013). Viscosity is a material, “bulk” food property that is measured as resistance to flow using rheological techniques (Bourne, 2002; Upadhyay et al., 2020). However, recent advances in food science identified the coefficient of sliding friction (CSF) as a distinct texture parameter important in oral lubrication. Unlike viscosity, sliding friction (i.e., the contact force measured between oral surfaces sliding against each other) results from dynamic interaction of food with oral surfaces, which is measured using tribological techniques that are sensitive to frictional rather than flow properties (Chen and Stokes, 2012; Stokes et al., 2013; Laguna et al., 2017). Separate neurons in primate OFC and insula encode CSF and viscosity (Rolls et al., 1999, 2003, 2018; Verhagen et al., 2003). An important question, therefore, is to what extent these two associated but distinct parameters — viscosity and sliding friction — contribute to oral-fat sensing and its rewarding properties.

Here we combined functional neuroimaging and behavioral experiments in human volunteers with food-engineering methods that quantified food textures in a realistic manner on biological oral surfaces (Huang et al., 2021). We investigated the neural mechanisms that link the oral processing of foods high in fat and sugar to subjective sensations and preferences. We hypothesized that CSF and viscosity would partly mediate the influence of fat content on subjective food valuations. We also hypothesized that neural reward systems process both the physical, oral-texture parameters of foods and the subjective sensations and valuations elicited by them. Finally, we sought to determine the predictive value of these measures by examining whether neural sensitivity to oral texture predicted subjects' eating behavior in a naturalistic context.

Materials and Methods

Participants

A total of 22 healthy volunteers participated in the full study (19-36 years old, mean age 25 years, SD = 4.8 years; 15 males, based on written self-report). Of 23 participants scanned, one dropped out during the first scanning session and was therefore excluded from the study. We initially aimed for a sample size of 24 participants, based on expected effect sizes from prior studies. However, it was not possible to complete data collection for all planned subjects because of Covid-related lockdown and restrictions in the United Kingdom at the time the study was conducted. Notably, each participant was scanned on 2 d to maximize the amount of data for decoding analyses. The sample size is comparable with previous papers in the field (Suzuki et al., 2017; DiFeliceantonio et al., 2018; Avery et al., 2020). Participants had a mean weight of 72.8 ± 14.2 kg and a mean BMI of 23.4 ± 2.56 kg/m2 (4 participants had a BMI between 25 and 29 kg/m2). Participants confirmed that they were not currently dieting, nonobese, and had no history of neurologic or psychiatric conditions. In post-study questionnaires, participants reported on average a liking for sweet and fatty foods (6.44 ± 1.92 and 6.28 ± 2.28, respectively; liking scales ranging from 0 to 10, labeled “not at all” and “very much,” respectively) and reported their milk preference for different fat levels (stated preferences: whole milk: 38.1% of participants; semi-skimmed milk: 33.3%; skimmed milk: 14.3%; nondairy milk: 14.3%). All participants provided written consent based on an information sheet. The study was approved by the Cambridge Local Research Ethics Committee. The different parts of the study were conducted at the Department of Physiology, Development and Neuroscience, University of Cambridge, at the Wolfson Brain Imaging Center, and at the Wellcome-MRC Institute of Metabolic Science Translational Research Facility.

Study design

To investigate behavioral and neural mechanisms underlying food preferences during oral food processing, we tested each participant (N = 22) in a series of experiments (see Fig. 1A). First, in a psychophysical test (day 1), subjects repeatedly sampled liquid food stimuli with well-controlled sensory and nutrient components and provided psychophysical ratings of food sensations; they also placed monetary bids to consume the foods in an economic auction-like task (Becker et al., 1964). Second, subjects participated in two fMRI scanning sessions (days 2 and 3) in which these tasks were repeated in the MRI scanner, including the oral sampling of the liquid foods (two scanning days were required to obtain sufficient data for multivoxel pattern fMRI analysis). Third, subjects participated in an ad libitum eating test (day 4) in which they made consumption choices for nutrient-controlled solid foods under naturalistic, life-like conditions. This design allowed us to investigate behavioral and neural responses to liquid foods under controlled experimental conditions and during realistic eating behavior involving solid foods.

Importantly, to orally process the liquid foods, subjects were trained to perform a standardized, left-right (or, in random alternation, right-left) tongue movement guided by a dynamic visual cue moving at a defined speed (replaced by a static left/right-pointing arrow for MRI scanning), and then hold their tongue still for several seconds before swallowing the liquids. This manipulation ensured (1) distribution of liquids around the oral cavity (Grabenhorst et al., 2010), (2) standardized oral mechanical and sensory stimulation across stimuli and subjects, and (3) oral sliding-friction stimulation at a defined speed that matched our tribology measurements (see below). This manipulation was important given our focus on oral texture and the importance of tongue movements in food-texture sensing (de Wijk et al., 2003, 2011; Koc et al., 2013).

Design of food stimuli

We designed a set of dairy-based liquid food stimuli (“milkshakes”) with controlled nutrient and energy content to produce variation in sensory food stimulation, subjective sensations and economic valuations (see Table 1). A 2 × 2 factorial design with sugar and fat as factors defined the following basic stimuli: (1) a low-fat, low-sugar liquid (LFLS), (2) a high-fat, low-sugar liquid (HFLS), (3) a low-fat, high-sugar liquid (LFHS), and (4) a high-fat, high-sugar liquid (HFHS). Energy content was lowest for the LFLS stimulus, intermediate for the HFLS and LFHS stimuli, and highest in the HFHS stimulus. In addition, we included (5) a low-fat, high-sugar and high-protein stimulus to assess fat-like texture produced by protein, (6) a HFHS based on soy rather than dairy milk to include both animal- and plant-based fat stimuli, and (7) a LFHS that was thickened by adding carboxymethyl cellulose (a food-thickening agent used in the food industry) to produce a fat-like texture. Following previous studies (de Araujo et al., 2003; Grabenhorst et al., 2008, 2010), “artificial saliva” was used as a rinse. This rinse stimulus contained the main ionic components of saliva (25 mm KCl + 2.5 mm NaHCO3). All stimuli except artificial saliva were vanilla-flavored (Sainsbury's Madagascan vanilla extract); temperature during the experiments was held constant by the use of ice-packs. In the first testing session, subjects were trained to associate a distinct, abstract visual cue with each of the eight liquid food stimuli by performing the psychophysical rating task in which visual cues were presented before liquid delivery on each trial. Subjects did not receive any feedback or explicit information regarding the nutrient content of the liquid food stimuli. The visual conditioned stimuli were presented before the liquid food stimuli were delivered in the psychophysical task and during MRI scanning to reduce participants' uncertainty about which liquid would be delivered a given trial. Visual cues were well separated from the liquid delivery by a jittered interstimulus interval of 4 ± 1 s. We used the same set of eight visual cues across subjects but changed the association between specific visual cues and liquids to avoid systematic visual confounds. We did not include an identification task to examine whether subjects could reliably recall the visual cues associated with the stimuli, although stimulus ratings were consistent between pretesting and scanning sessions and across the two scanning sessions (correlation across stimuli between pretest and scanning sessions: sweetness: R = 0.937, p = 0.0018; thickness: R = 0.978, p = 0.0001; correlation across stimuli between two scanning sessions: sweetness: R = 0.997, p = 9.1 × 10−7; thickness: R = 0.977, p = 0.0001).

View this table:
  • View inline
  • View popup
Table 1.

Composition of experimental liquid food stimulia

Delivery of liquid food stimuli

In order to deliver the seven liquid food stimuli and the rinse stimulus while the participant underwent MRI scanning, custom-made computer-controlled peristaltic pumps were used. The pumps were set up in the scanner control room, where the experimenter, the radiographer, and the testing laptop were located. Food-grade silicone tubing was connected to each pump output and threaded through the connecting aperture into the scanner room. Participants held a soft, custom-made mouthpiece (a modified infant food feeder; Losuya) between their lips that connected the eight tubes without mixing the stimuli. The pumps were the same as used in a previous human fMRI experiment involving similar liquid food stimuli (Zangemeister et al., 2016). The pumps were controlled by the laptop with the use of a National Instruments card. The rotational speed of the pump rotor was constant within each pump. We controlled the volume of the liquid delivered by specifying the length of time the pump was activated. From previous studies involving liquid reward delivery in the fMRI (Grabenhorst et al., 2010; Zangemeister et al., 2016), we aimed to deliver between 0.75 and 1 ml to participants as this was shown to be sufficient to activate neural systems processing taste, flavor, and texture. We carefully and regularly calibrated the pumps to ensure constant and matched delivery volumes across stimuli.

Instructed tongue movements

In order to ensure even dispersal of the stimuli and rinse, and to ensure that texture-related subjective ratings were not confounded by variability in tongue movement, participants were instructed in how to move their tongue during the tasting period. In the behavioral pretesting task, a cursor appeared in the center of the screen during the tasting and rinse period. This cursor would move either to the left, then to the right before coming back to the center or to the right, then the left before coming back to the center, and participants were asked to move their tongue across the palate mirroring the cursor's movement. The direction in which the cursor moved was randomized for every trial and taste period such that the participant could not predict the direction before the appearance of the cursor. This task trained the participants to move their tongues with a specific speed and reduced the likelihood of results arising from differences in tongue movement. During the fMRI scanning task, the tongue movement cursor was replaced with a static arrow. This replacement reduced the need for participants to actively track the cursor movement across the screen with their eyes to minimize visual effects.

Task structure for psychophysical testing and MRI scanning

We used a similar trial structure for psychophysical testing and MRI scanning. The event timing for the fMRI task is shown in Figure 1B. Each trial started with the display of one of the eight visual conditioned stimuli for 1.5 s, followed by an interstimulus interval during which a central fixation cross was displayed for 4 ± 1 s. The visual cue served to prepare the subjects for the subsequent delivery of a liquid food stimulus. Delivery of one of the food stimuli (chosen in a random permuted sequence) was cued by a colored fixation cross which stayed on the screen for 7 s. During this period, a cursor or arrowhead pointing to the left or right was presented next to the fixation spot for 1 s to guide the subject's tongue movement. As described above, subjects were pretrained to make one left-right tongue movement under visual guidance and then hold still until they were cued to swallow. After this 7 s tasting period, the letter S was presented for 1.5 s to cue subjects to swallow. After this period, two response periods followed in each of which a labeled visual analog rating scale was presented and subjects moved a bar to the appropriate point on the scale using a button box.

During psychophysical testing (day 1), subjects made five different rating responses by rating sweetness, thickness, oiliness, protein content, willingness to pay (WTP), and made a binary choice response whether or not the stimulus contained fat. For the rating scales, we used descriptors previously used in relation to taste and textural attributes, including sweetness, thickness, and oiliness (defined as an oil-like feel) (Szczesniak, 1963; Bourne, 2002; de Araujo and Rolls, 2004; Kadohisa et al., 2005). Subjects received written instructions with real-life examples in how to use these ratings scales to evaluate the food stimuli. Given our focus on oral-texture properties, the present manuscript focuses on analyzing ratings of sweetness, thickness, oiliness, and WTP. During MRI scanning, restrictions on scan time required a simplified design involving two types of trials in which different ratings were used in random alternation: subjects either rated the food stimuli on separate scales for sweetness and thickness (trial Type 1) or placed a WTP bid and made a binary choice about whether the food stimulus contained fat (trial Type 2). Each rating period was 4 s long. After the second rating period ended, a fixation cross was displayed for 4 ± 1 s, followed by the visual conditioned cue associated with the rinse stimulus, which was displayed for 1.5 s. The delivery and timing of the rinse period, including the tongue movement and swallowing periods, were identical to the liquid food delivery, except that no rating scales were used for the rinse stimulus. Thus, the trial ended with the swallowing of the rinse and was followed by an intertrial interval of 4 ± 1 s during which a fixation cross was displayed. A typical trial lasted ∼44 s (depending on the jittered periods).

MRI scanning was performed in 6 runs of ∼15 min each, to provide sufficient numbers of runs for training and testing for the multivariate fMRI data analysis. One subject was instead tested with 5 runs of ∼18 min. The 6 runs were performed on 2 separate days <1 week apart. Importantly, to create a balanced design for the multivariate analysis, each food stimulus was delivered 3 times within each run. Thus, in total each food stimulus was delivered 18 times. Occasional trials in which subjects reported that no stimulus was delivered (e.g., because of liquid backflow) were excluded from analysis. This general protocol and design have been used successfully in previous studies to investigate activations and their relation to subjective ratings in cortical areas (Grabenhorst et al., 2008, 2010; Grabenhorst and Rolls, 2014; Zangemeister et al., 2016).

WTP measure

We measured the subjective, economic value of the food stimuli in a Becker-DeGroot-Marschak auction-like task (Becker et al., 1964), that has previously been used in functional neuroimaging and behavioral experiments with visual and liquid food stimuli (Plassmann et al., 2007; Zangemeister et al., 2016; Suzuki et al., 2017). Participants provided a bid of 0-10 credits on the liquid food stimulus they received on a given trial. At the beginning of the behavioral pretest, as well as at the beginning of each MRI scanning run, the participant received 100 credits endowment to spend on their bids. Using this endowment, participants placed bids on the liquid food stimulus they had just sampled on the same trial. Participants were instructed that they would bid to consume 250 ml of the just-sampled food after the experiment, and that their bid would be successful if it exceeded a computer-generated bid on the same trial (which was generated from a uniform distribution with replacement). We provided written instructions to the participants which explained that the best strategy was to bid exactly the amount that reflected their subjective value of the food stimulus; in that way, they would ensure to obtain the preferred food for the price they were willing to pay and avoid obtaining nonpreferred foods. Participants placed bids by moving a cursor on a WTP scale. When participants lost their bid, they lost none of their budget, although when they won, the computer-generated random price was deducted from their budget.

Rheological measurements of viscosity

We measured the viscosity of the experimental liquid food stimuli at the Rheology Center of the Department of Chemical Engineering, University of Cambridge. Measurements were performed using a Rheometric Scientific ARES controlled strain rheometer (TA Instruments) with the following Couette geometry: cup diameter = 34.0 mm, bob diameter = 32.0 mm, bob length = 34.0 mm. Viscosity was measured by carrying out shear-rate sweeps from 100-0.1 s−1 or 1 s−1 (reverse sweep) after allowing samples to equilibrate at the experimental temperature for 300 s. We measured the samples at different temperatures including 18.5°C (the temperature at which stimuli were delivered into the subject's mouth), 25°C, and 37.5°C (body temperature). We used measurements taken at 18.5°C for our data analyses, based on the temperature of the testing liquids (17 ± 2°C) measured during the experiment, which were kept constant with ice packs. When the liquid sample reached the determined temperature, the cup rotated first clockwise and then counterclockwise. In the time-sweeping measurement, the cup was rotated in a clockwise direction and the shear rate was fixed at 50 s−1; rheological properties of liquid stimuli have been shown at these settings to be related to oral viscosity evaluations (Wood, 1968).

Tribological measurements of the CSF

We measured the CSF for our liquid rewards at the Department of Engineering, University of Cambridge. To reflect realistic lubrication conditions in the oral cavity, we devised a custom-designed tribometer using pig tongues as biological contacting surfaces (see Fig. 1C), as described previously (Huang et al., 2021). The design involved a flat aluminum platform that held a fixed base consisting of the (nominally flat) pig tongue. An upper moving tongue-tip tissue was attached with superglue onto a dome-shaped slider of radius 100 mm. Thus, the top and bottom tissues contacted with a nominal point contact, avoiding any issues with alignment during the sliding process. The anatomically upper surface of each tongue was used for each of the two contacting surfaces. We mounted the dome-shaped slider using low-friction bushings onto a track containing two rails, pivoted at one end. We used a counterweight to balance out the weight of the rail elements of the track. Therefore, the load through the contact between the tongues was determined by the weight of the dome-shaped slider and tongue tip specimen (2.58 ± 0.07 N); these did not vary as the slider moved along the track. The slider was attached via a pulley to an Instron 5544 Universal Testing machine (Instron). We confirmed in preliminary tests that pulley friction was negligible. During testing, we loaded the moving tissue against the fixed base pig tongue with the testing liquids interfaced as lubricating layers. The Instron machine then imposed a fixed velocity (v = 16 mm/s) to the slider, while measuring the traction force using a load cell attached to the Instron machine. Accordingly, this design measured the sliding friction between the test liquids and the oral tissues; it thereby approximated biological oral sensing conditions. As we maintained a constant velocity during the test, according to Newton's First Law of Motion, as follows: F−μN=0 μ=traction force (F)loading force (N)=traction force measured by the tensile machinetotal weight of the slider and the fixed tongue tip where F is the applied force (traction) by the Instron machine, N is the loading force perpendicular to the contact surface (normal force), and μ (the CSF) is the ratio of the traction force and the perpendicular loading force.

On the day before testing, we obtained fresh pig tongues from a local butcher (Leech & Sons) which we gently rinsed with water to remove residual blood and tissue fluids. We then placed the superficial 1 cm thick, anterior 18 cm of the tongues onto a flat contact surface that fitted onto the testing platform. We preserved the processed pig-tongue slices in isotonic saline buffer (PBS, 1×, pH 7.4) in a freezer under 4°C overnight. On the testing day, we prepared the contact surfaces by gluing one 18 cm tongue on the base platform and another tongue tip (5 cm) on the surface of the dome-shaped slider. We weighed the dome-shaped slider with the attached tongue tip to give the loading force for later calculation of the CSF. Before each measurement, both tongue surfaces were rinsed with 10 ml of PBS 3 times to remove residual testing liquids and hydrate the tongue surfaces. Next, we loaded 30 ml of the testing liquids between the pig tongues and pulled the slider from the posterior tongue bases forward to the anterior tongue tip (see Fig. 1C). Measurements were performed with testing liquids at room temperature. For practical reasons, we were unable to measure CSF at different temperatures because the measurements using fresh pig tongues were time-consuming and each pig tongue could only be used for a limited number of measurements before the surface degraded. The procedures were approved by the Departmental Safety Office, Department of Engineering, University of Cambridge, based on Control of Substances Hazardous to Health and biohazard risk assessments.

For the final testing of the stimuli, we measured each liquid with triplicate repeats using two pairs of pig tongues in opposite measuring orders to control for possible carryover effects (low-fat to high-fat liquids and reverse). We first averaged the triplicate measurements for each liquid and divided them by the corresponding loading force to obtain the CSF along the tongue surface. We selected the anterior 5-7 cm of the tongues as our analysis window because the mechanosensory receptors are located mainly within the anterior two-thirds of tongue (Ogawa et al., 1989; Toda and Taoka, 2002); more anterior tongue tips were too thin for stable measurements. Because inevitable variations in tongue conditions influenced the CSF measurement, we normalized for each pair of testing tongues the measured CSF with the CSF of water obtained from the same pair of tongues. This normalization adjusted the offset in absolute measurement because of variations of the tongues and provided comparable results between different tongue pairs. Finally, we averaged the two normalized CSF obtained from the different tongues for our analyses.

Naturalistic food-preference test

This test was performed at the Translational Research Facility of the Wellcome-MRC Institute for Metabolic Sciences at Addenbrooke's Hospital, Cambridge. It was possible to test 20 of the 22 subjects that participated in psychophysical testing and MRI scanning sessions also in this naturalistic food-preference test; 2 subjects did not return for this separate testing session. The stimulus design and testing procedure followed a previous study (van der Klaauw et al., 2016). We conducted an ad libitum eating test in which participants were allowed to choose freely and repeatedly between three curries that were approximately matched in terms of visual appearance and spicing but had varying nutrient compositions. The curries were based on the chicken korma from a previous study (van der Klaauw et al., 2016), except that the protein component was from a popular meat replacement (Quorn) to not exclude vegetarian participants. The nutrient compositions of the curries were designed such that they had three distinct fat and sugar levels, following a similar ratio to that of the liquid milkshakes used in the behavioral pretesting and fMRI component of the study (see Table 6). The curries were finely chopped and mixed to match their visual appearance. The test was conducted as follows.

After their second MRI scanning session, participants were invited to a fourth visit. Participants were informed that they would perform a behavioral rating task involving solid food in addition to completing questionnaires on their eating habits. They were also informed that, as the rating of solid foods may change depending on their satiety levels, they should adhere to specific breakfast suggestions on the testing day, which they received in writing and which would contain ∼250 kcal. They were then scheduled for a 1 h testing session around lunchtime (between 12:00 noon to 1:00 PM).

Upon arrival, participants completed a pretesting questionnaire asking about their hunger and thirst levels in addition to their adherence to the breakfast guidelines provided. Subsequently, the participant was given three small portions of the experimental curry dishes and asked to rate them on specific psychophysical scales. The positioning of the curry dishes was always randomized by the research chef who prepared the foods, and the researcher and participant were blind to the positioning until after the experiment. The participant was also provided with 100 ml of water to cleanse the palate between samples.

After the tasting test, the participant was moved to a different corner of the room to perform a computer-based Stroop Color and Word Test, which served as a distractor task, for 2 min. During this time, the researcher informed an assistant with a text message that the participant was ready for lunch. At the end of the Stroop Test, the assistant entered the room and offered them a lunch meal, informing them that they would need to be back in 30 min to complete a final questionnaire. The participant was then taken to a separate lounge and shown the three curries in separate large dishes (in randomized positions) and was informed that they should sample all of them and consume as much as they wanted. We encouraged to sample the three dishes without explicitly mentioning that the three meals were different and instructed them to eat until they were comfortably full. After 30 min, the participant was taken to the testing room again to complete a final questionnaire and for debriefing and remuneration. The consumed amount of each meal was covertly weighed.

Analysis of behavioral data

To analyze the rating data, we first performed one-way repeated-measures ANOVA with food stimulus as factor. Separate ANOVAs were performed for ratings of sweetness, thickness, oiliness, and WTP. Next, to test how ratings depended on the sugar, fat, and protein content of stimuli, and on the texture parameters CSF and viscosity, we used mixed-effects multilinear regression analysis within the GLM framework (implemented with the fitglme function, MATLAB) and specified the subject number as the categorical group variable to account for subject-specific effects (random effects). We adopted the global model in which we estimated both the main effects and random effects of all the relevant regressors. The response variable was the trial-by-trial rating (sweetness, thickness, oiliness, WTP). We implemented the regression shown in Tables 2 and 3. Based on previous studies (de Araujo et al., 2003; Grabenhorst et al., 2008, 2010; Zangemeister et al., 2016), we had no a priori hypotheses with respect to sex-based behavioral and neural effects and therefore did not perform sex-based analyses. We used Bonferroni-corrected p values to determine significance in these analyses, determined by the number of times a variable was used in statistical tests (e.g., sugar was used 10 times).

View this table:
  • View inline
  • View popup
Table 2.

Regression models of sweetness, thickness, and oiliness ratings

View this table:
  • View inline
  • View popup
Table 3.

Regression models of WTP bids

fMRI data acquisition

Functional imaging data were acquired using a 3T Skyra (Siemens) scanner at the Wolfson Brain Imaging Center (WBIC, University of Cambridge). Whole-brain T2*-weighted EPIs were acquired with a repetition time of 3000 ms, echo time of 30 ms, flip angle of 90°, and 51 axial oblique slices with 3 mm isotropic resolution. Subjects were scanned in six runs. In each run, 300 volumes were acquired in an imaging time of 15 min 10 s. Importantly, each liquid food stimulus was presented 3 times in each run, to achieve a balanced design for the leave-one-run-out analysis approach described below. A high-resolution structural MPRAGE scan for normalization purposes was also acquired (voxel size, 1 × 1 × 1 mm; TR, 2300 ms; TE, 2.98 ms; inversion time, 900 ms; flip angle, 9°; total scan time, 5 min 3 s).

fMRI data analysis

We used the Statistical Parametric Mapping package to analyze the neuroimaging data (SPM12; Wellcome Trust Center for Neuroimaging). We discarded the first six volumes of each run to allow for magnetic saturation effects. We preprocessed the data by applying slice-timing correction, realignment, and reslicing (separately for the fMRI data from day 1 and day 2). We then segmented the data to extract white matter, gray matter, and cerebrospinal fluid and followed by coregistering the day 1 and day 2 data using the T1-weighted structural scans from each day. The data preprocessed in this way were used for multivariate voxel pattern analysis (MVPA) in individual subjects without applying normalization or spatial smoothing to the data. For univariate ROI analyses, we normalized the scans to the standard MNI coordinate system, and then applied smoothing using a Gaussian kernel with the FWHM of 6 mm. We also removed global effects from the fMRI time series using voxel-level linear model of the global signal to minimize effects correlated with global fluctuations of the signal (Macey et al., 2004).

MVPA

We used linear support vector machine (SVM) regression to test for encoding of continuous variables, such as the CSF, viscosity, thickness, and oiliness ratings, and the WTP ratings in local activity patterns distributed across voxels, as implemented in the Decoding Toolbox (Hebart et al., 2014). Decoding of textural parameters (CSF and viscosity) was performed with SVM regression using the CSF and viscosity value for each liquid food stimulus, whereas decoding of subjective psychophysical ratings and WTP bids was performed on the within-subject mean rating given to each food stimulus. Using a leave-one-out cross-validation method, we trained the classifier on five scan runs within each subject to predict the value of the variable of interest and then tested on the remaining run to predict the variable from distributed multivoxel activity patterns; this approach was repeated 6 times to obtain an average statistical map of z-transformed correlation coefficients describing the relationship between voxel activity and the variable of interest. The whole-brain MVPA was performed using a searchlight method on unsmoothed realigned data in the subject's native space. We used a linear kernel and default regularization parameter of c = 1. A 9 mm sphere was defined as a searchlight and used for MVPA, after which the searchlight moves onto a neighboring voxel. Exploratory analyses with different sphere sizes suggested that a 9 mm searchlight radius provided robust and accurate results, while small changes in sphere size did not qualitatively alter the results. We report the specific results for our main findings also with a 4 mm sphere; no additional brain areas were identified with this 4 mm sphere analysis. The resultant accuracy maps were subsequently normalized into the MNI space, following the protocol used in similar MVPA studies (Kahnt et al., 2011; Howard and Kahnt, 2017). In the second-level (group, random-effects) stage, the normalized, smoothed statistical maps resulting from the first-level MVPA were entered into one-sample t tests, resulting in group-level statistical parametric maps.

Conjunction analyses

Conjunction analyses (Nichols et al., 2005) were used to test for common encoding of different effects within the same brain region. This was achieved by testing if all of the variables of interest were significantly encoded (p < 0.05), in which case the null hypothesis that at least one variable was not encoded could be rejected. Although primarily used to compare shared regions in contrast files from univariate GLMs (Friston et al., 1999, 2005), conjunction analysis has also been used with accuracy maps from MVPA to indicate shared regions in taste encoding (Avery et al., 2020).

Finite impulse response analyses

We used a time-resolved SVR regression as previously described (Kahnt et al., 2011) to examine the encoding of our variables of interest at different time points in a trial (see Fig. 4F). Because we were interested not only in the activations related to the delivery of liquid food stimulus but also in possible further encoding in the period following the liquid delivery, a finite impulse response analysis was performed as implemented in SPM12, which makes no assumption about the specific time course of the HRF (Henson, 2003; Grabenhorst and Rolls, 2008). We used δ functions in the finite impulse response analysis that were spaced at intervals of 2 s and starting at the onset of the liquid delivery and spanning a period of 10 s. This basis set considers each time bin after stimulus onset individually to model the fMRI signal and can capture any possible shape of response function up to a given frequency limit. In this model, the parameter estimate for each time bin represents the average fMRI signal at that time. The activity patterns in these 2 s time bins were then used for a searchlight MVPA as described above.

Statistical significance testing

For all fMRI analyses, we report effects that survive correction for multiple comparisons across the whole brain using a significance level of p < 0.05 (family-wise error) at cluster level; this threshold was applied to statistical maps that were displayed using a cluster-forming threshold of p < 0.001 (uncorrected) at voxel level (Eklund et al., 2016) with cluster extent of k = 15 or more contiguous voxels. In addition, we used small-volume correction (p < 0.05, cluster-level) in structures for which we had strong a priori hypotheses based on previous studies, including the mid-to-lateral OFC, oral somatosensory cortex (oSSC), and ventromedial PFC. Small-volume corrections were performed in spheres of 10 mm radius for cortical areas. The spheres were centred on specific coordinates reported in previous studies as follows: mid-to-lateral OFC [32, 34, −14] and oSSC [66, −18, 12], both taken from the main effects of previous studies on oral food processing using liquid stimuli with defined sugar and fat content (Grabenhorst et al., 2010, 2014), and ventromedial PFC [−2, 40, −4] taken from a meta-analysis of value-based decision-making (Clithero and Rangel, 2014). The coordinates for small-volume corrections were determined during study design but were not preregistered.

Univariate ROI analysis

To ensure that statistical inference in our ROI analyses was not circular, we followed approaches used in previous studies (Behrens et al., 2008; Zangemeister et al., 2016; Rosenthal-von der Putten et al., 2019). Specifically, we used a leave-one-subject-out method in which we re-estimated a second-level analysis 22 times, each time leaving out one subject to define the ROI coordinates for the left-out subject. We then extracted the signal from the subject-specific coordinates defined in this way. Thus, the data which we used for the ROI analysis were independent from those used to define the coordinates for extracting the signal. ROIs were determined during study design but were not preregistered. Following data extraction, we applied a high-pass filter with a cut-off period of 128 s. The data were then z-normalized, oversampled by a factor of 10 using sinc-interpolation, and separated into trials to produce a matrix of trials against time. We generated separate matrices for each event of interest (rating trials, choice trials). We then fitted GLMs to each oversampled time point across trials separately in each subject. We fitted a GLM containing as main regressors the two key oral-texture variables CSF and viscosity. In addition, the GLMs included fat, sugar, energy content, motion parameters, and session constants as covariates of no interest. To test the statistical significance, we entered individual-subject coefficients into one-sample t tests (random-effects analysis, p < 0.05) and calculated group averages and SEs for each time point across participants, yielding the across-subject effect size time courses shown in the figures. These mean effect size time courses are shown Figure 3B.

Testing relationships between neural measures and naturalistic eating behavior

To test for relationships between neural effect sizes and behavioral variables across subjects, we extracted neural effects sizes within each individual subject form peak coordinates defined using the leave-one-subject-out procedure described above. Specifically, we extracted the neural effects sizes for encoding of CSF and viscosity from the statistical decoding maps, based on group cluster coordinates, which we defined independently by re-estimating the second-level analysis 22 times, each time leaving out one subject to define the ROI coordinates for the left-out subject. We then fitted a multilinear regression model with the neural effect sizes of OFC to CSF and viscosity as regressors to the amount of fat eaten in the ad libitum test as the dependent variable. To produce Figure 5D, we plotted the relationship between the dependent variable and its predicted values as given by the coefficients of the fitted regression model.

Results

Study design

To investigate the behavioral and neural mechanisms that link oral-sensory food processing to economic preferences, we tested each participant in a series of experiments (Fig. 1A). First, in a psychophysical test (day 1), subjects (N = 22) repeatedly sampled liquid food stimuli with controlled sensory and nutrient components and provided psychophysical ratings of sweetness, thickness, and oiliness. As a measure of the subjective, “economic” reward value of the foods, they also placed WTP bids to consume the foods in an incentive-compatible auction-like task (Becker et al., 1964). Second, subjects participated in two fMRI scanning sessions (days 2 and 3) in which they performed these tasks in the MRI scanner while orally sampling the liquid foods (two scanning days were required to collect sufficient data for multivoxel-pattern fMRI analyses). Liquid foods were delivered orally via food-grade tubes and computer-controlled peristaltic pumps, following established protocols (de Araujo and Rolls, 2004; Grabenhorst et al., 2010; Grabenhorst and Rolls, 2014; Zangemeister et al., 2016). Third, subjects participated in an ad libitum eating test (day 4) in which they made consumption choices for nutrient-controlled solid foods under naturalistic, life-like conditions. This design allowed us to investigate behavioral and neural responses to nutrient-defined foods both under controlled experimental conditions and during realistic eating behavior.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Study design and oral-sensory measurements for experimental foods. A, Sequence of experiments in each participant. B, Trial structure for controlled sampling and oral processing of liquid foods in psychophysical and fMRI experiments. Liquid foods were delivered orally under computer control via peristaltic pumps. Subjects performed a pretrained left-right (or right-left, alternating randomly trial by trial) tongue movement guided by a visual cue to standardize oral processing. C, Sliding friction and viscosity as key texture variables in food oral processing. Left, Sliding friction results from oral-surface movements, with liquid foods acting as lubricant; viscosity is a “bulk” material food property. Middle, Custom tribometer measures the CSF of experimental liquid foods in a biologically realistic manner using sliding pig tongues. Right, Rotational rheometer measures the viscosity of liquid foods. D, Validation of tribological measures of sliding friction in a series of basic fat-containing liquids. CSF reflects liquids' fat content. E, Linear regression of CSF and viscosity on fat content in a large set of experimental liquid foods. Colored foods were used in the present experiments. Inset, Relationship between CSF and viscosity. HPHS, Low-fat, high-sugar, high-protein; Soy, soy-based HFHS; CMC, low-fat, high-sugar thickened with carboxymethyl cellulose. CMC represents a nonfat artificially thickened control stimulus that was excluded from the regression.

Importantly, for the psychophysical and fMRI experiments, subjects were trained to orally process the liquid foods by performing a standardized left-right tongue movement guided by a dynamic visual cue moving at a defined speed (replaced by a static cue for MRI scanning) and then hold their tongue still for several seconds before swallowing the liquids (Fig. 1B). This manipulation served to (1) distribute liquid foods around the oral cavity (de Araujo and Rolls, 2004; Grabenhorst et al., 2010), and (2) produce standardized oral-mechanical stimulation across stimuli and subjects, including sliding friction at a defined speed that matched our tribological measurements (see below). These factors were important given our focus on oral texture and the role of tongue movements in food-texture sensing (de Wijk et al., 2003, 2011; Koc et al., 2013).

Design of food stimuli

We designed a set of dairy-based liquid food stimuli (“milkshakes”) with controlled nutrient and energy content to produce meaningful variation in oral-sensory stimulation, subjective sensations, and reward valuations (Table 1). We focused on sugar and fat because of their relevance for human overeating and obesity, and their role in determining sensory food qualities (Alonso-Alonso et al., 2015; Carreiro et al., 2016; Rolls, 2020; Huang et al., 2021). A 2 × 2 factorial design with sugar and fat as factors defined the following basic stimuli: (1) a LFLS liquid, (2) a HFLS liquid, (3) a LFHS liquid, and (4) a HFHS liquid. Energy content was lowest for the LFLS stimulus, intermediate for the HFLS and LFHS stimuli, and highest for the HFHS stimulus. In addition, we included (5) a LFHS and high-protein stimulus to include fat-like texture produced by protein, (6) a HFHS liquid based on soy rather than dairy milk to include both animal- and plant-based fat stimuli, and (7) a LFHS liquid that was thickened by adding carboxymethyl cellulose (a food-thickening agent used in the food industry) to produce a fat-like texture. Following previous studies (de Araujo et al., 2003; de Araujo and Rolls, 2004; Grabenhorst et al., 2008, 2010), “artificial saliva,” a largely neutral stimulus approximating the ionic composition of human saliva, was used as a rinse on each trial, and in select control analyses as described below. All stimuli except artificial saliva were vanilla-flavored. Stimuli were served at a temperature of ∼17°C that was held constant across stimuli with ice-packs.

Measuring food texture: viscosity and CSF

We reasoned that the oral sensations produced by food textures influenced subjects' food valuations and related neural processing. Therefore, a main aim of our study was to relate subjects' behavioral and neural responses to foods to oral-texture parameters quantified with rheological and tribological food-engineering methods. Two specific texture parameters have been implicated in oral fat-sensing (Fig. 1C): viscosity and the CSF (Chen and Stokes, 2012; Chojnicka-Paszun et al., 2012; Stokes et al., 2013; Laguna et al., 2017).

Viscosity is a material, “bulk” food property that is measured as resistance to flow (i.e., the quotient of shear stress and shear rate) in a rotational rheometer (Fig. 1C); it has been shown to increase with increasing fat-droplet concentration, which is higher in fatty liquids (Bourne, 2002; de Wijk et al., 2011; Chen and Stokes, 2012; Laguna et al., 2017). By contrast, CSF is a systems property that emerges from the interaction of food with oral surfaces and is measured as the resistance to relative motion between contacting surfaces using tribological techniques (Fig. 1C); fatty liquids produce lower CSF, likely by coalescing fat droplets forming a load-bearing, lubricating film on oral surfaces (Dresselhuis et al., 2008b; Chojnicka-Paszun et al., 2012; Torres et al., 2018; Soltanahmadi et al., 2023). Thus, CSF and viscosity reflect physically distinct aspects of food texture (Chen and Stokes, 2012). We hypothesized that both CSF and viscosity would mediate the influence of fat content on oral sensations and food valuations, and related neural activity. Testing this hypothesis required measuring viscosity and CSF for our liquid food stimuli under biologically plausible conditions, as described next.

To approximate the friction conditions in the oral cavity, which involve the tongue sliding against the palate, we used our recently developed tribometer (Fig. 1C) (Huang et al., 2021) with parallel-sliding fresh pig tongues as biological surfaces that approximate the softness and surface properties of the human tongue. In using biological tissues, this approach differs from standard tribometers involving circular-rotating glass-metal surfaces. Validation tests using a series of basic fatty liquids confirmed that decreases in CSF reliably reflected increased stimulus fat content (Fig. 1D). In a large set of more complex liquid food stimuli that included the foods used in the present fMRI study, CSF correlated negatively with the fat content of the stimuli (R = −0.78, p = 4.6 × 10−5), confirming lower sliding friction for high-fat foods, whereas viscosity correlated positively with fat content (R = 0.83, p = 4.4 × 10−6), confirming that high-fat foods were more viscous (Fig. 1E). By contrast, sugar content in the stimuli showed weaker and nonsignificant correlation with CSF (R = −0.42, p = 0.063) and viscosity (R = 0.12, p = 0.59). Although CSF and viscosity themselves were correlated (R = −0.78, p = 3.4 × 10−5), the two parameters reflect fundamentally different aspects of food texture that become relevant during different stages of oral processing, with viscosity being a material food property and CSF resulting from food-oral surface interactions (Chen and Stokes, 2012; Stokes et al., 2013). Consistently, in the present stimulus set, CSF and viscosity discriminated between the liquid foods in different ways (Fig. 1E, colored stimuli): low CSF for the high-fat stimuli (HFLS, HFHS, and soy) and the fat-like carboxymethyl cellulose (CMC) differed substantially from the high CSF of the low-fat stimuli (LFLS, LFHS, and protein); in contrast, viscosity additionally distinguished the HFLS and HFHS stimuli from CMC and soy. As shown below, CSF and viscosity also had distinguishable relationships with behavioral and neural measures.

Thus, the oral-texture parameters viscosity and sliding friction tracked the fat content of the liquid foods in partly different ways and therefore likely influenced subjects' food sensations and valuations.

Psychophysical ratings of food stimuli

Subjects provided psychophysical ratings on visual scales with numerical labels while sampling the liquid foods under controlled conditions and visually guided tongue movements. In their ratings, subjects discriminated the food stimuli in terms of sweetness, thickness, and oiliness (ANOVAs; sweetness: F(6,147) = 32.81, p = 6.9 × 10−25; thickness: F(6,147) = 63.36, p = 2.5 × 10−38; oiliness: F(6,147) = 14.49, p = 6.0 × 10−13) (Fig. 2A). Notably, stimulus ratings were consistent between pretesting and scanning sessions and across the two scanning sessions (correlation across stimuli between pretest and scanning sessions: sweetness: R = 0.937, p = 0.0018; thickness: R = 0.978, p = 0.0001; correlation across stimuli between two scanning sessions: sweetness: R = 0.997, p = 9.1 × 10−7; thickness: R = 0.977, p = 0.0001).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Psychophysical ratings and economic WTP bids for experimental foods. A, Ratings of sweetness, thickness, and oiliness sensations for all stimuli. Error bars indicate mean ± SEM. Data points show individual subjects. HPHS, Low-fat, high-sugar, high-protein; SOY, soy-based HFHS; CMC, low-fat, high-sugar thickened with carboxymethyl cellulose. B, Thickness and oiliness ratings are explained by a combination of oral-texture variables viscosity and sliding friction (mixed-effects multilinear regression, subjects as random factor). **p < 0.005, Bonferroni-corrected. C, WTP bids, obtained in an incentive-compatible Becker-DeGroot-Marschak auction-like task as a measure of subjects' economic food valuations. D, Mediation analysis. Path diagram describing relationships between nutrient content, texture parameters, and WTP bids (mixed-effects multilinear regression). The influence of fat and sugar content on WTP bids was decomposed into indirect effects mediated by texture variables viscosity and sliding friction and direct effects. Significance of path coefficients derived from bootstrap (1000 iterations). Protein effects were included in the model but were not significant and are not shown. E, WTP depended on oral sensations of sweetness, thickness, and oiliness (mixed-effects multilinear regression; **p < 0.005).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Neural encoding of specific oral-sensory food properties. A, Activity pattern across voxels in OFC reflected sliding friction of liquid foods (MNI coordinates: [36, 42, −16]; z = 4.08; p = 0.012; small-volume corrected) and activity patterns across voxels in oSSC reflected viscosity (MNI coordinates: [54, −10, 22]; z = 4.34; p < 0.001, whole-brain corrected). Shown are decoding-accuracy maps derived from whole-brain searchlight analyses that used SVM regression to decode each texture variable from multivoxel activity patterns in a 9 mm sphere. Maps are shown at p = 0.001 with extent threshold of 15 voxels. B, ROI regression of OFC cluster activity on sliding friction (t(21) = 2.11, p = 0.046; one-sample t test across subject-specific regression betas) and viscosity (not significant). Activity timeseries was extracted from coordinates identified in a previous study. Yellow shaded region represents the period in which the regression was significant; absence of yellow shaded region represents that the regression was not significant. C, ROI regression of oSSC cluster activity on viscosity (t(21) = 2.30, p = 0.031; one-sample t test) and sliding friction (not significant). Activity timeseries was extracted from coordinates identified in a previous study. D, OFC activity patterns encoded both sliding friction and oiliness in a conjunction analysis (MNI coordinates: [36, 42, −18]; z = 3.68; p = 0.002, whole-brain corrected). E, oSSC activity patterns encoded both viscosity and thickness in a conjunction analysis (MNI coordinates: [58, −12, 16]; z = 3.76; p = 0.021, whole-brain corrected).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Neural encoding of economic food values and related oral-sensory properties. A, Activity pattern in OFC reflected subjective, economic food reward values, measured as WTP bids during oral processing of liquid foods (MNI coordinates: [38, 40, −14]; z = 3.83; p = 0.010, small-volume corrected). Decoding-accuracy maps derived from whole-brain searchlight analysis using SVM regression to decode WTP from multivoxel activity patterns in a 9 mm sphere. Maps are shown at p = 0.001 with extent threshold of 15 voxels. B, OFC activity patterns during oral processing encoded sliding friction, oiliness, and WTP in a three-way conjunction analysis (MNI coordinates: [36, 42, −14]; z = 3.74; p = 0.002, small-volume corrected). The OFC was the only brain area to show a significant three-way conjunction (right, transparent brain map). C, OFC encoding of sliding friction remained significant when performing the decoding on CSF residuals, after regressing out variance components because of WTP bids ([36, 42, −14]; z = 3.99; p = 0.010), small-volume corrected). D, Neural decoding accuracies for sliding friction and WTP in OFC were uncorrelated. E, Activity patterns in pgACC encoded WTP only after oral processing before subjects placed their bids (MNI coordinates: [−2, 44, 16]; z = 3.73; p = 0.004, whole-brain corrected). F, Time-resolved decoding of WTP bids in OFC and pgACC. Decoding accuracy betas obtained from a finite impulse response analysis performed in 2 s bins. The OFC β was significant in bin 3 (t(21) = 2.60, p = 0.016); the pgACC β was significant in bin 5 (t(21) = 3.05, p = 0.006).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Neural encoding of oral-sensory food properties in OFC predicts preferences in a naturalistic eating test. A, Design of food stimuli. Subjects could freely and repeatedly choose between three nutrient-controlled curry meals in an ad libitum eating test. B, Variation in preference (consumed amounts) for different foods across subjects. C, Modeling fat preference in the eating test from OFC neural texture sensitivity measured in the MRI scanner. We fitted a multilinear regression model with the neural sensitivity of OFC to CSF and viscosity, as well as subjects' BMI, to the amount of fat eaten in the ad libitum test. Neural betas were extracted from decoding-accuracy maps using independently determined OFC peak coordinates from a leave-one-subject-out cross-validation procedure (see Materials and Methods). The model provided a significant fit to the amount of fat eaten in the eating test (F(3,16) = 4.57, full model p = 0.017) with both neural texture sensitivities contributing significantly (CSF regressor: p = 0.026; viscosity regressor: p = 0.027; BMI regressor: p = 0.264). D, Relationship between the amount of fat eaten as predicted by the model based on OFC texture sensitivity and the actual amount of fat consumed by the subjects during the eating test.

Mixed-effects multilinear regressions (with subject as random factor) showed that sweetness ratings were influenced by the foods' sugar and fat content (sugar: p = 2.0 × 10−16; fat: p = 1.5 × 10−10; protein: p = 0.112), while thickness and oiliness ratings depended mostly on fat content (Table 2; thickness: sugar: p = 0.0007, fat: p = 4.9 × 10−115, protein: p = 0.114; oiliness: sugar: p = 0.0003, fat: p = 6.0 × 10−12, protein: p = 0.159). Thickness and oiliness ratings were of particular interest, as we expected them to reflect sensations related to oral-texture parameters, which in turn depended on fat content. We therefore substituted the fat content regressor with texture properties to test whether CSF and viscosity partly explained thickness and oiliness ratings. Both thickness and oiliness ratings were explained by combinations of the two texture parameters CSF and viscosity, in addition to a significant sugar effect (Fig. 2B; Table 2; thickness: sugar: p = 1.2 × 10−5, CSF: p = 1.2 × 10−33, viscosity: p = 7.0 × 10−16, protein: p = 0.017, oiliness: sugar: p = 0.0023, CSF: p = 0.0008, viscosity: p = 0.024, protein: p = 0.771). Specifically, higher reported thickness and oiliness were associated with lower CSF (consistent with the negative relationship between CSF and fat content, compare Fig. 1E). Viscosity also had a negative relationship with thickness and oiliness, which likely reflected relatively lower thickness and oiliness ratings for the non–fat-containing but viscous CMC stimulus. Indeed, when CMC was removed from the analysis, CSF coefficients remained significant and negative, but the viscosity coefficient on oiliness was no longer significant. Thickness and oiliness ratings were only partly correlated (shared variance, R2 = 0.249), indicating that they reflected partly distinct oral-texture sensations. In both thickness and oiliness regressions, the CSF regressor had higher t values and lower p values compared with viscosity (Table 2). CSF also explained more variance in both thickness and oiliness ratings: removing CSF from the regression model (with remaining regressors sugar, viscosity, protein) decreased the proportion of explained variance by 65.7% and 64.6% for thickness and oiliness ratings, respectively; removing viscosity from the model decreased the explained variance only by 34.9% and 50.0%. For comparison, removing sugar from the regression model while retaining both CSF and viscosity reduced the explained variance by 7.7% and 13.2% for thickness and oiliness ratings, respectively. Therefore, CSF was relatively more important in explaining rated oral sensations than viscosity.

Thus, the liquid food stimuli produced oral sensations that reflected the foods' nutrient composition and related texture qualities, with both CSF and viscosity contributing to sensations of oiliness and thickness. We next examined how these food properties and related oral sensations influenced subjective reward valuations, as measured by subjects' WTP to consume the foods.

Subjective reward valuations of foods depend on texture parameters and oral sensations

Subjects performed an incentive-compatible auction-like task (Becker et al., 1964) in which they placed WTP bids on each of the liquid foods to win the opportunity to consume 250 ml of the offered liquid after the testing session. Bids were consequential, in that subjects were informed that the bid from a randomly chosen trial would be implemented. WTP bids varied significantly across food stimuli (Fig. 2C, left; ANOVA: F(6,147) = 30.3, p = 1.7 × 10−23). On average, subjects placed the highest bids on the HFHS stimulus, which was highest in sugar, fat, and energy content (p < 0.001, Bonferroni-corrected post hoc t tests). Despite this general effect, there were also considerable interindividual differences in food valuation (Fig. 2C, gray data points), indicating that WTP bids reflected individual preferences. Notably, bids were consistent between pretesting and scanning and across scanning sessions (correlation between pretesting and scanning sessions: R = 0.976, p = 0.0001; correlation between scanning sessions: R = 0.967, p = 0.0003).

Mixed-effects regression confirmed that higher bids were associated with higher fat and sugar content of the food stimuli (Table 3; sugar: p = 4.8 × 10−5; fat: p = 1.2 × 10−30); effects of protein were limited by design (p = 0.043). A previous study found a fat–carbohydrate interaction effect on WTP-measured valuation in response to visually cued foods (DiFeliceantonio et al., 2018). We tested for a fat–sugar interaction effect on WTP in an additional analysis, but this interaction was not significant (p = 0.159, mixed-effect regression) while fat and sugar main effects remained significant (both p < 1.0 × 10−5).

We tested the effects of oral-texture parameters on WTP bids in a mediation analysis, using Structural Equation Modeling implemented by a series of mixed-effects regression models. This approach estimated both direct effects of fat and sugar on WTP and indirect effects that were mediated by the oral-texture variables CSF and viscosity. Regression-derived path coefficients (Fig. 2D) showed that CSF and viscosity had significant negative effects on WTP (CSF: p = 9.6 × 10−7; viscosity: p = 2.5 × 10−7), showing that subjects placed higher bids on foods that produced low sliding friction and that were less viscous. Direct effects of fat and sugar were significantly reduced by the inclusion of CSF and viscosity in the regression, although they remained significant (sugar: p = 2.5 × 10−15; fat: p = 2.9 × 10−11). Thus, CSF and viscosity partially mediated the effects of fat and sugar on WTP bids.

Notably, as for the oiliness and thickness ratings, CSF seemed relatively more important in explaining variance in WTP bids: removing CSF from the regression model (with remaining regressors sugar, viscosity, protein) reduced the proportion of explained variance by 72.6%, whereas removing viscosity reduced the explained variance by 58.7%. For comparison, removing sugar from the regression model while retaining both CSF and viscosity reduced the explained variance by 34.7%. Model comparisons showed that a simpler model based on the energy content of the food stimuli was insufficient for explaining WTP bids (AIC and BIC favored more complex models that included fat and sugar as well as oral-texture variables, Table 3).

Finally, we examined how WTP bids were related to psychophysically reported oral sensations using sweetness, thickness, and oiliness ratings as regressors. Both sweetness and thickness had significant positive influences on bids (sweetness: p = 1.6 × 10−9; thickness: p = 2.7 × 10−16); oiliness had a significant negative influence on bids (p = 0.00003) and therefore seemed to partially capture the negative effect of sliding friction on WTP bids reported above (Fig. 2E). The behavioral analyses described above were designed to test specific hypotheses regarding the relationships between particular variables. For descriptive purposes, we also report the full cross-correlation matrix between these variables (Table 4).

View this table:
  • View inline
  • View popup
Table 4.

Cross-correlation matrix of food components, psychophysical ratings, and WTP across all liquid food stimuli (mean Pearson correlation coefficients across participants)

Thus, subjective food valuations depended on the sugar and fat content of the stimuli. The effect of fat on valuations was partly mediated by the oral-texture parameters CSF and viscosity, and their subjective correlates oiliness and thickness. We next examined how oral-texture parameters and food valuations were related to neural activity patterns in different brain regions.

Neural encoding of food texture: sliding friction and viscosity

We used fMRI to measure neural responses when subjects orally processed the liquid foods and applied multivariate analysis approaches to identify neural encoding of our key oral-texture variables. Specifically, we used a whole-brain searchlight approach that decoded CSF and viscosity from local activity patterns (fMRI signal intensities in a local group of voxels) while the liquid food was in the subjects' mouth on each trial; the searchlight approach performed this decoding across the whole brain by moving a sphere of 9 mm radius that defined the voxels in the decoding sample systematically through the brain. To do so, we trained a SVM regression algorithm to decode oral-texture variables and subjective stimulus ratings from the fMRI signal intensity patterns distributed across a 9 mm sphere that defined the voxels for the current decoding sample; the sphere was systematically moved across the whole brain to produce decoding accuracies for each brain area (Hebart et al., 2014; Haynes, 2015). We performed cross-validation by systematically training the classifier on data from all but one scanning runs to test performance on the remaining run (leave-one-run-out).

Cross-validated decoding revealed that activity patterns in a region of the mid-to-lateral OFC reflected our key variable CSF (Fig. 3A; MNI coordinates: [36, 42, −16]; z = 4.08; p = 0.012; small-volume corrected based on coordinates taken from a previous study (Grabenhorst et al., 2010) (Table 5). This finding implied that during the oral processing of liquid foods, OFC activity patterns reflected the sliding friction produced by the interaction of a specific liquid food with oral surfaces (with high-fat foods producing lower sliding friction, see Fig. 1). In addition, CSF could be decoded from a region of the postcentral gyrus consistent with the location of the oSSC (Grabenhorst and Rolls, 2014), although this effect was weaker than in OFC (MNI coordinates: [58, −14, 14]; z = 3.20; p = 0.044; small-volume corrected based on coordinates from a previous study (Grabenhorst and Rolls, 2014) (Table 5). We found that overlapping areas of OFC and oSSC showed significant relationships to the material food property viscosity (Fig. 3B, Table 5; OFC: MNI coordinates: [40, 38, −12]; z = 4.72; p < 0.001, whole-brain corrected; oSSC: MNI coordinates: [54, −10, 22]; z = 4.34; p < 0.001, whole-brain corrected). Thus, activity patterns in OFC and oSSC also reflected the viscosity of the liquid foods, which differed based on the foods' fat content (with high-fat foods being more viscous, see Fig. 1). These results were confirmed when decoding was performed with a smaller sphere of 4 mm radius (CSF, OFC: z = 4.17, p = 0.042; CSF, oSSC: z = 3.39; p = 0.048; viscosity, OFC: z = 4.02; p = 0.030; viscosity, oSSC: z = 4.48; p = 0.034). The results were also not explained by energy content: in control analyses, decoding based on CSF and viscosity residuals after regressing out energy showed significant effects in OFC and oSSC (CSF, OFC: z = 4.75, p = 0.045; CSF, oSSC: z = 3.94; p = 0.005; viscosity, OFC: z = 4.92; p < 0.0001; viscosity, oSSC: z = 4.23; p < 0.0001). Decoding on CSF residuals after removing shared variance with viscosity (and vice versa) was unsuccessful (in either direction), indicating that regressing out shared variance removed important information. We performed control analyses by repeating the decoding procedure for CSF and viscosity with a delayed onset time that captured neural activity measured during the rinse period on the same trials, after subjects had swallowed the liquid foods. No significant effects were found in the rinse period even at a more lenient threshold (p < 0.005, uncorrected), which confirmed that neural encoding of CSF and viscosity was specific to the time when the liquid food was in the mouth, and that there was no lingering stimulus (e.g., because of residual mouthcoating) that could have enabled texture decoding from subsequent neural activity.

View this table:
  • View inline
  • View popup
Table 5.

fMRI resultsa

View this table:
  • View inline
  • View popup
Table 6.

Nutrient and caloric composition of the foods in the ad libitum eating testa

To further distinguish between the neural encoding of CSF and viscosity in OFC and oSSC, we used ROI analyses by extracting the cluster-based fMRI signal from predefined coordinates of OFC (Grabenhorst et al., 2010) and oSSC (Grabenhorst and Rolls, 2014). We then performed time-resolved multilinear regression analysis on the extracted neural signals in which both CSF and viscosity were entered as covariates and thus competed to explain variance in neural activity (other covariates were fat, sugar, and energy content; see Materials and Methods). This analysis indicated that activity in the identified OFC and oSSC clusters contributed in different ways to the encoding of oral texture: cluster-based OFC activity showed a significant relationship to CSF but not viscosity (Fig. 3B, CSF: t(21) = 2.11; p = 0.046; one-sample t test across subject-specific regression betas), whereas cluster-based oSSC activity showed the opposite pattern of a significant relationship to viscosity but not CSF (Fig. 3C, viscosity: t(21) = 2.30; p = 0.031). Although this result does not demonstrate a formal double dissociation, it shows that CSF and viscosity contributed differently to explaining OFC and oSSC activity when these variables competed in the same regression model. Simple fat and sugar effects were not significant in both OFC and oSSC. An additional analysis which included only fat and sugar regressors showed no significant effect in OFC and a significant sugar effect in oSSC (t(21) = 2.65, p = 0.015), indicating that objective fat and sugar content per se did not explain OFC and oSSC activity well. The observed partial neural separation of CSF and viscosity is consistent with concepts from the food-engineering literature whereby distinct aspects of oral food processing are thought to reflect the bulk (mass-related) food property viscosity and the thin-film (lubrication-related) property CSF (Chen and Stokes, 2012; Stokes et al., 2013; Shewan et al., 2020).

Next, we trained the SVM regression algorithm to decode subjective food ratings related to oral-texture variables from neural activity patterns. Multivariate decoding of oiliness and thickness ratings given to the different food stimuli showed that OFC activity patterns reflected subjective oiliness (Table 5; MNI coordinates: [34, 40, −8]; z = 3.58; p = 0.005, small-volume corrected), whereas oSSC activity patterns reflected subjective thickness (MNI coordinates: [60, −14, 18]; z = 3.54; p = 0.012, small-volume corrected). We performed conjunction analyses (Nichols et al., 2005) between the statistical decoding-accuracy maps for texture variables and ratings to test whether the same areas that encoded texture variables also encoded subjective oral sensations. Conjunction analyses showed that activity patterns in the mid-to-lateral OFC area encoded both CSF and subjective oiliness (Fig. 3D; MNI coordinates: [36, 42, −18]; z = 3.68; p = 0.002, whole-brain corrected), whereas oSSC activity patterns reflected both viscosity and subjective thickness (Fig. 3E; MNI coordinates: [58, −12, 16]; z = 3.76; p = 0.021, whole-brain corrected).

Thus, OFC and oSSC encoded partly distinct food-texture properties and related sensations during the oral processing of liquid food stimuli: OFC was particularly sensitive to the sliding friction produced by moving the liquid foods between the oral surfaces, whereas oSSC was particularly sensitive to the foods' viscosity. Given that these texture variables and related oral sensations contributed to subjective reward valuations (Fig. 2), we next investigated neural activity in relation to subjects' WTP bids.

Neural encoding of subjective food valuations

Multivariate decoding using SVM regression showed that activity patterns in the mid-to-lateral OFC region during oral food processing reflected subjects' economic food valuations, as measured by WTP bids (Fig. 4A, Table 5; MNI coordinates: [38, 40, −14]; z = 3.83; p = 0.010, small-volume corrected; decoding with 4 mm radius sphere: z = 3.89; p = 0.023). OFC activity therefore not only reflected the physical food property of sliding friction but also individuals' subjective food valuations, which depended on sliding friction. In a three-way conjunction analysis between CSF, oiliness, and WTP, the OFC was the only brain area to significantly encode all three variables (Fig. 4B, MNI coordinates: [36, 42, −14]; z = 3.74; p = 0.002, small-volume corrected). Thus, neural representations of CSF, oiliness, and food valuations overlapped in OFC.

We performed strict tests of whether the encoding of sliding friction in OFC was separate from its encoding of value, to rule out that the conjunction result described above was entirely driven by shared variance components between WTP bids and CSF. In these tests, we first regressed out the variance in CSF that was related to WTP bids in each subject, and then repeated the multivariate decoding on the CSF residuals. Decoding with these CSF residuals after correcting for WTP showed a significant effect of sliding friction in OFC (Fig. 4C; [36, 42, −14]; z = 3.99; p = 0.010, small-volume corrected). In a converse analysis, decoding with individual subjects' WTP bids after regressing out CSF also revealed a significant effect of WTP bids in OFC ([34, 46, −18]; z = 3.91; p = 0.026). Thus, both CSF and WTP bids explained distinct variance portions in OFC activity patterns. Further, neural decoding accuracies for CSF and WTP bids (extracted from predefined OFC coordinates: [32, 34, −14]) (Grabenhorst et al., 2010) were uncorrelated (Fig. 4D), indicating partly independent representations. Consistently, a model comparison on the univariate OFC signal (compare Fig. 3B) showed that a model based on subjective variables (sweetness, thickness, oiliness, WTP) did not provide a significantly better fit than a model based on objective variables (sugar, fat, viscosity, CSF, energy; t test on AIC, t(20) = 1.76, p = 0.093). The same results were found for the univariate oSSC signal (t(20) = 1.57, p = 0.131). Thus, the OFC appeared to be a key integration site during oral food processing that combined distinct information about physical food-texture properties and subjective food reward valuations.

Different from the OFC, encoding of WTP bids during oral processing was not found in the oSSC, which therefore encoded viscosity and thickness as shown above but not value. We also did not find evidence for WTP encoding during oral processing in the medial PFC, including the ventromedial PFC, which was previously shown to encode WTP for visually cued foods (Plassmann et al., 2007; Lebreton et al., 2009; Suzuki et al., 2017; DiFeliceantonio et al., 2018). Based on these previous studies, we reasoned that medial prefrontal activity might encode subjective value at a later trial stage when subjects were asked to report their bids in response to visual cues. Testing for encoding of WTP with a delayed onset time of seven seconds (just after the tasting period and before the rating period), we found a strong effect related to WTP in the pregenual anterior cingulate cortex (pgACC) (Fig. 4E, MNI coordinates: [−2, 44, 16]; z = 3.73; p = 0.004, whole-brain corrected; decoding with 4 mm radius sphere: z = 3.36; p = 0.027 at lower threshold of p = 0.005 uncorrected), a region that partly overlaps with ventromedial PFC and has previously been implicated in food valuation (Grabenhorst et al., 2010). The pgACC cluster extended into adjacent ventromedial PFC with a subpeak at [−10, 44, 4]. A time-resolved decoding analysis (see Materials and Methods) further confirmed that OFC encoded subjective value primarily during food oral processing, whereas pgACC encoded subjective value only after oral processing, before subjects were required to report their food evaluations (Fig. 4F). A conjunction analysis showed that neither OFC nor pgACC jointly encoded the oral-texture variables CSF and viscosity alongside WTP bids in this later period following oral processing.

Thus, during oral processing, activity patterns in OFC reflected subjects' food reward valuations alongside oral-texture variables. This initial integration of texture parameters with subjective values in the OFC during oral processing was followed by a “pure” value signal in pgACC shortly before subjects reported their valuations.

Naturalistic food-preference test and relationship to neural measures

We next asked whether the neural texture sensitivity of the OFC, measured during controlled experimental conditions in the MRI scanner, was related to subjects' food preferences under life-like, naturalistic conditions. The same subjects that participated in MRI scanning (N = 20) also participated in a three choice ad libitum eating test during which they freely and repeatedly chose between three solid foods, offered in different serving trays, to consume their lunch. The design of the foods and testing conditions followed a previous study which demonstrated differences in fat preference between healthy participants and participants with obesity related to mutations in the Melanocortin-4-receptor (van der Klaauw et al., 2016). The foods were matched in flavor (a vegetarian version of chicken korma) and visual appearance (van der Klaauw et al., 2016) but differed in fat and sugar composition and included a low-fat, high-sugar stimulus, a HFLS stimulus, and a mid-fat mid-sugar (MFMS) stimulus (Fig. 5A; Table 6). Before the ad libitum eating test, subjects sampled and rated a small quantity of each of the three curry foods (performed in a separate room under conditions designed to prevent carryover effects between ratings and the eating test). Ratings suggested that subjects discriminated the foods in terms of sweetness, oiliness, savoriness, and pleasantness (ANOVAs with post hoc t tests corrected for multiple comparisons, all p < 0.001), but not in terms of thickness and saltiness (p > 0.075).

During the ad libitum eating test, the offered foods elicited substantial interindividual variation in how much subjects consumed of each of the foods (Fig. 5B, ANOVA p < 1.8 × 10−5). On average, subjects consumed more of the HFLS and MFMS foods compared with the LFHS food (both p < 0.0001, t tests, Bonferroni-corrected). As our study focused on preferences for fat, the key variable of interest was the amount of fat consumed in the eating test. If oral texture sensations contributed to subjects' fat preferences, then neural sensitivity to texture variables in the OFC might be related to subjects' eating phenotypes measured under naturalistic conditions. To test this hypothesis, we fitted a multilinear regression model with the neural sensitivity of OFC to CSF and viscosity as regressors to the amount of fat eaten in the ad libitum test. Neural betas were extracted from decoding-accuracy maps using independently determined OFC peak coordinates from a leave-one-subject-out cross-validation procedure (see Materials and Methods).

Consistent with our prediction, neural sensitivity of OFC to CSF and viscosity explained individual differences in eating phenotypes: subjects whose OFC encoded CSF and viscosity more strongly in the scanner subsequently consumed more fat in the life-like eating test (Fig. 5C,D; F(3,16) = 4.57, full model: p = 0.017; CSF regressor: p = 0.026; viscosity regressor: p = 0.027). Similar results were found in single linear regressions (CSF: p = 0.046, viscosity: p = 0.031) and in a random-effects regression that treated subject as random factor (CSF: p = 3.3 × 10−31, viscosity: p = 0.046). The regression model based on neural sensitivity to oral texture explained a substantial portion of the variance in consumed fat between subjects (R2 = 0.461). The effect of neural texture sensitivity on fat preference was independent of subjects' BMI, which we included as a control covariate (BMI regressor: p = 0.264). In a model that also included sex and age covariates, the CSF regressor remained significant (p = 0.0028) but not the viscosity regressor (p = 0.081) with no effect of sex (p = 0.515) but a significant effect of age (p = 0.022). As a further control, no significant relationship was found between neural texture sensitivity and the total amount of food eaten (F(3,16) = 2.63, full model: p = 0.0855; CSF regressor: p = 0.122; viscosity regressor: p = 0.0521; BMI regressor: p = 0.499), indicating that the effect was specific to fat preference. To further test the robustness of the identified predictive relationship, we performed an out-of-sample prediction by predicting each subject's consumed amount of fat based on a neural texture model that did not include the predicted subject's data (leave-one-subject-out resampling). The model remained significant (p < 0.05) in 19 of 20 tests with the proportion of explained variance in fat consumption ranging from R2 = 0.377 to R2 = 0.595 (mean ± SEM, R2 = 0.466 ± 0.01). We also tested for relationships between fat consumption and neural texture sensitivity of oSSC, which encoded oral-texture variables but, unlike OFC, did not encode subjective reward valuations. Different from OFC, texture sensitivity of oSSC was not related to fat consumption in the eating test (F(3,16) = 0.556, full model: p = 0.649; CSF regressor: p = 0.424; viscosity regressor: p = 0.622; BMI regressor: p = 0.670). As described in the Introduction, CSF and viscosity reflect related but mechanically and conceptually distinct oral-texture parameters; accordingly, we included both variables in the analysis above. For completeness, we also compared models that regressed the amount of fat eaten in the ad libitum eating test on OFC sensitivity to both CSF and viscosity, only CSF, or only viscosity. Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) showed lowest values for the model based on both CSF and viscosity (AIC/BIC: −154.93/−151.05) compared with models based only on CSF (AIC/BIC: −151.96/−149.97) or only on viscosity (AIC/BIC: −152.54/−150.55).

Thus, neural sensitivity of the OFC to oral-texture parameters, measured under controlled experimental conditions, explained individual differences in subjects' preference for high-fat foods, measured on a separate testing day during naturalistic conditions.

Discussion

We measured behavioral and neural responses when healthy volunteers consumed liquid foods that varied in nutrient (fat, sugar) composition and sensory, oral-texture qualities. Our main finding is that a region of the human OFC responded to foods in the mouth by encoding the CSF, a mechanical oral-texture parameter that resulted from interactions between liquid foods and oral surfaces.

Specifically, OFC activity patterns reflected the low sliding friction and oiliness sensations produced by high-fat foods, as well as the subjective economic valuations that depended on sliding friction. The neural encoding of sliding friction in OFC was partly independent from its encoding of economic value, suggesting separate but overlapping neural representations. These effects were not explained by viscosity, a distinct texture parameter that reflects material food properties rather than food-oral surface interactions and separately influenced value and neural activity. The neural integration of food-texture parameters with economic value was specific to OFC and not found in oral SSC or pgACC, which separately encoded food textures and value, respectively. Importantly, oral-texture sensitivity of OFC predicted the amount of fat that subjects consumed in a naturalistic eating test. These results identify the human OFC as a key brain structure that detects oral fat from sliding friction and translates this sensory information into economics values that guide food preferences.

The conditions under which CSF was measured and related to behavior and neural activity are distinctive features of our study. Using tribological techniques that approximated oral surface conditions with pig tongues (Huang et al., 2021), we showed that the CSF reliably indicated the fat content of liquid foods. The physical mechanism by which fatty liquids lower oral sliding friction involves coalescence of fat droplets on oral surfaces which form an adhering fat layer that likely influences oral sensation (Dresselhuis et al., 2008b; Chojnicka-Paszun et al., 2012; Shewan et al., 2020). Our approach followed earlier studies that used pig tongues to link food fat content to CSF (Dresselhuis et al., 2008a; Huang et al., 2021) and applied this method for the first time to human fMRI and behavioral measures of economic food valuations. To match oral processing and tribological testing conditions, we trained subjects to perform controlled tongue movements with a speed approximating the measuring conditions of our tribometer.

Subjects' behavioral responses to the foods confirmed that CSF contributed to sensations of oiliness and thickness, and partly mediated the influence of fat on subjective valuations. Previous studies linked human food ratings to oral-texture attributes (Kokini, 1987; Kadohisa et al., 2005; Laguna et al., 2017) but did not study CSF valuations with incentive compatible economic auctions. We found that lower sliding friction typically elicited higher valuations, indicating that subjects placed higher bids for the opportunity to consume foods with a smooth, oily texture typical of high-fat stimuli. These oral-texture effects were not accounted for by sugar, sweetness, and caloric content, which were regression covariates. While our mediation model linked CSF and viscosity directly to subjective WTP bids, a further refinement could model subjective fat and sugar perception as an additional, intermediate step. Different from a previous study using visual food cues (DiFeliceantonio et al., 2018), we found no fat–sugar interaction effect on valuation, which may be because of our use of oral food stimuli or our specific stimulus set.

In previous fMRI studies, fat-containing and texture-controlled foods activated OFC, oSSC, and pgACC (de Araujo and Rolls, 2004; Grabenhorst et al., 2010; Eldeghaidy et al., 2011; Stice et al., 2011; Grabenhorst and Rolls, 2014; Eldeghaidy et al., 2016; Zangemeister et al., 2016). We extended these approaches by combining (1) novel food-engineering methods that quantified oral-texture parameters on biological tissues, (2) incentive-compatible measures of economic value, and (3) a naturalistic food-preference test. This approach allowed us to uncover a functionally and anatomically specific role of the human mid-to-lateral OFC in oral food processing. This OFC region parametrically encoded sliding friction shown using both multivoxel-pattern decoding, which reflects local distributed voxel-activity patterns, and univariate multilinear regression, which reflects the dominant cluster signal of these voxels (Haynes et al., 2007). Although oSSC also encoded oral texture and pgACC encoded economic value, the OFC was the only area in which these variables were jointly encoded. Previous studies identified roles of the insula in processing food attributes evoked either by oral stimulation or visual food pictures, including at 7 Tesla MRI field strengths (Avery et al., 2020, 2021). Our findings suggest that food texture and related value evoked by oral food stimuli are primarily represented in OFC and oSSC, although future studies could investigate possible insula contributions with targeted approaches.

The OFC has long been implicated in food processing and related functions in reward and emotion (Rolls et al., 1999, 2003; Tremblay and Schultz, 1999; Rolls and Grabenhorst, 2008; Grabenhorst et al., 2010; Grabenhorst and Rolls, 2011; Howard and Kahnt, 2017). It receives direct inputs from oSSC (Carmichael and Price, 1995) and anterior insular taste cortex (Baylis et al., 1995). The same region of human mid-to-lateral OFC was previously implicated in oral-fat processing (Grabenhorst et al., 2010; Grabenhorst and Rolls, 2014), but the underlying oral-mechanical parameters and their relation to economic value and eating behavior remained unclear. Single-neuron recordings in macaque OFC demonstrated detailed representations of sensory food properties (Rolls et al., 1999, 2003; Verhagen et al., 2003; Kadohisa et al., 2005), including neurons that encoded either CSF or viscosity, which were intermingled in OFC (Rolls et al., 2018). Although the spatiotemporal resolution of fMRI cannot recover such detailed representations (Logothetis, 2008), our multivoxel decoding results showed that human OFC was sensitive to both CSF and viscosity. Fully dissociating CSF and viscosity in a decoding analysis was not possible; nevertheless, the univariate fMRI signal in OFC was dominated by CSF rather than viscosity, indicating a robust sliding-friction representation in the human OFC (although a formal double dissociation was not demonstrated). We note that cleanly dissociating CSF and viscosity is a current unresolved challenge in the food-engineering field.

The present results inform our understanding of neural reward mechanisms that contribute to overconsumption of high-fat foods and obesity. We found that the sensitivity with which the OFC encoded the oral-texture variables CSF and viscosity predicted subjects' fat preferences in a naturalistic eating test. This finding is consistent with studies that linked structure and function of the human OFC to eating behavior (Morys et al., 2022). For example, abnormal eating in frontotemporal dementia has been associated with OFC pathology (Woolley et al., 2007), variation in OFC gray-matter volume has been linked to obesity (Shott et al., 2015; Hall et al., 2023), and resting-state OFC functional connectivity is associated with liking for sweet foods and BMI (Rolls et al., 2023). Our findings suggest a possible mechanism by which the OFC could contribute to eating behavior: the OFC may translate information about the “objective,” material properties of foods, including their nutrient content and oral texture, into subjective sensations and economic values that guide eating behavior. The mapping from nutrients and sensory food properties to subjective values likely depends on individual-specific synaptic weights stored on OFC neurons (Rolls and Grabenhorst, 2008). Indeed, our behavioral results showed that oral-texture parameters partly mediated the effect of fat on individuals' subjective valuations, matching our recent findings in macaques (Huang et al., 2021). We note that the relationship between fat preferences and OFC encoding of oral texture was found in a relatively small sample of N = 20 participants. Although convergent evidence points to a key role of the OFC in food-intake control (reviewed above) and our result was robust to several control analyses, future studies should replicate this finding in a larger sample.

Our findings have implications for a central question in decision neuroscience: how are subjective values neurally constructed from sensory inputs? Suzuki et al. (2017) showed that multivoxel activity in lateral OFC in response to visual food cues reflected nutrient components and both medial and lateral OFC regions reflected economic value. Other studies found that values evoked by visual food cues were encoded in ventromedial PFC (Lebreton et al., 2009; Tang et al., 2014). Here, by examining oral food processing, we provide evidence that subjective value may be initially constructed from oral-sensory inputs in the mid-to-lateral OFC before being encoded as sensory-independent value in pgACC. Specifically, whereas OFC activity encoded WTP while subjects orally sampled liquid foods, pgACC activity encoded WTP only after subjects had swallowed the foods and prepared to report their valuations. Notably, only OFC but not pgACC jointly represented value and oral-texture signals and oiliness sensations on which valuations were based. Thus, our findings point to a key role of the mid-to-lateral OFC in constructing food value from oral-sensory inputs.

Our findings have implications for food engineering and food design. Advances in food science identified sliding friction as a key factor in the design and acceptability of food (Chen and Stokes, 2012; Wang and Chen, 2017; Upadhyay et al., 2020). Here we showed that the frictional properties of liquid foods influence not only subjective perceptions but also economic valuations measured in incentive-compatible auctions, naturalistic eating behavior, and neural reward responses. Our approach using biologically plausible tribology tools to quantify food textures and relate them to human perceptions, economic valuations, and neural measures could help validate the design of foods that have both healthy nutrient compositions and attractive oral-texture properties.

In conclusion, our findings show that the human OFC responded to foods in the mouth by encoding the CSF, a physical oral-texture parameter that mediated the influence of fat content on subjects' economic food valuations. Activity patterns in OFC combined sliding-friction signals with representations of subjective oiliness and economic value and explained fat preferences during naturalistic eating. Future studies could investigate how the presently identified neural mechanism for linking oral-textural properties of high-fat foods to eating behavior could contribute to overeating and obesity.

Footnotes

  • This work was supported by the Wellcome Trust and the Royal Society Sir Henry Dale Fellowship Grants 206207/Z/17/Z and 206207/Z/17/A to F.G. Part of this study was conducted in the Wellcome-MRC Institute for Metabolic Sciences Translational Research Facility, which is supported by Wellcome Trust Major Award 208363/Z/17/Z. This work was supported in whole, or in part, by the Wellcome Trust. For the purpose of Open Access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission. Data and code underlying this research are available from the corresponding author upon reasonable request. We thank Wolfram Schultz and his group for support; Alejandro Casabuena-Rodriguez, Simon Marshall, and Stefan Savage for support with tribology; Carol Keel for support with the naturalistic food-preference test; and Jae-Chang Kim for discussions.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Fabian Grabenhorst at fabian.grabenhorst{at}psy.ox.ac.uk

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    1. Alonso-Alonso M,
    2. Woods SC,
    3. Pelchat M,
    4. Grigson PS,
    5. Stice E,
    6. Farooqi S,
    7. Khoo CS,
    8. Mattes RD,
    9. Beauchamp GK
    (2015) Food reward system: current perspectives and future research needs. Nutr Rev 73:296–307. https://doi.org/10.1093/nutrit/nuv002 pmid:26011903
    OpenUrlCrossRefPubMed
  2. ↵
    1. Avery JA,
    2. Liu AG,
    3. Ingeholm JE,
    4. Riddell CD,
    5. Gotts SJ,
    6. Martin A
    (2020) Taste quality representation in the human brain. J Neurosci 40:1042–1052. https://doi.org/10.1523/JNEUROSCI.1751-19.2019 pmid:31836661
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Avery JA,
    2. Liu AG,
    3. Ingeholm JE,
    4. Gotts SJ,
    5. Martin A
    (2021) Viewing images of foods evokes taste quality-specific activity in gustatory insular cortex. Proc Natl Acad Sci USA 118:e2010932118. https://doi.org/10.1073/pnas.2010932118
    OpenUrlCrossRef
  4. ↵
    1. Baylis LL,
    2. Rolls ET,
    3. Baylis GC
    (1995) Afferent connections of the orbitofrontal cortex taste area of the primate. Neuroscience 64:801–812. https://doi.org/10.1016/0306-4522(94)00449-f pmid:7715788
    OpenUrlCrossRefPubMed
  5. ↵
    1. Becker G,
    2. DeGroot M,
    3. Marschak J
    (1964) Measuring utility by a single-response sequential method. Behav Sci 9:226–232. https://doi.org/10.1002/bs.3830090304 pmid:5888778
    OpenUrlCrossRefPubMed
  6. ↵
    1. Behrens TE,
    2. Hunt LT,
    3. Woolrich MW,
    4. Rushworth MF
    (2008) Associative learning of social value. Nature 456:245–249. https://doi.org/10.1038/nature07538 pmid:19005555
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bourne MC
    (2002) Food texture and viscosity: concept and measurement, Ed 2. London: Academic.
  8. ↵
    1. Carmichael ST,
    2. Price JL
    (1995) Sensory and premotor connections of the orbital and medial prefrontal cortex of macaque monkeys. J Comp Neurol 363:642–664. https://doi.org/10.1002/cne.903630409 pmid:8847422
    OpenUrlCrossRefPubMed
  9. ↵
    1. Carreiro AL,
    2. Dhillon J,
    3. Gordon S,
    4. Higgins KA,
    5. Jacobs AG,
    6. McArthur BM,
    7. Redan BW,
    8. Rivera RL,
    9. Schmidt LR,
    10. Mattes RD
    (2016) The macronutrients, appetite, and energy intake. Annu Rev Nutr 36:73–103. https://doi.org/10.1146/annurev-nutr-121415-112624 pmid:27431364
    OpenUrlPubMed
  10. ↵
    1. Chen JS,
    2. Stokes JR
    (2012) Rheology and tribology: two distinctive regimes of food texture sensation. Trends Food Sci Technol 25:4–12. https://doi.org/10.1016/j.tifs.2011.11.006
    OpenUrlCrossRef
  11. ↵
    1. Chojnicka-Paszun A,
    2. de Jongh HH,
    3. de Kruif CG
    (2012) Sensory perception and lubrication properties of milk: influence of fat content. Int Dairy J 26:15–22. https://doi.org/10.1016/j.idairyj.2012.04.003
    OpenUrl
  12. ↵
    1. Clithero JA,
    2. Rangel A
    (2014) Informatic parcellation of the network involved in the computation of subjective value. Soc Cogn Affect Neurosci 9:1289–1302. https://doi.org/10.1093/scan/nst106 pmid:23887811
    OpenUrlCrossRefPubMed
  13. ↵
    1. de Araujo IE,
    2. Rolls ET
    (2004) The representation in the human brain of food texture and oral fat. J Neurosci 24:3086–3093. https://doi.org/10.1523/JNEUROSCI.0130-04.2004 pmid:15044548
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. de Araujo IE,
    2. Rolls ET,
    3. Kringelbach ML,
    4. McGlone F,
    5. Phillips N
    (2003) Taste-olfactory convergence, and the representation of the pleasantness of flavour, in the human brain. Eur J Neurosci 18:2059–2068. https://doi.org/10.1046/j.1460-9568.2003.02915.x pmid:14622239
    OpenUrlCrossRefPubMed
  15. ↵
    1. de Wijk RA,
    2. Engelen L,
    3. Prinz JF
    (2003) The role of intra-oral manipulation in the perception of sensory attributes. Appetite 40:1–7. https://doi.org/10.1016/s0195-6663(02)00172-1 pmid:12631499
    OpenUrlCrossRefPubMed
  16. ↵
    1. de Wijk RA,
    2. Janssen AM,
    3. Prinz JF
    (2011) Oral movements and the perception of semi-solid foods. Physiol Behav 104:423–428. https://doi.org/10.1016/j.physbeh.2011.04.037 pmid:21570419
    OpenUrlPubMed
  17. ↵
    1. DiFeliceantonio AG,
    2. Coppin G,
    3. Rigoux L,
    4. Thanarajah SE,
    5. Dagher A,
    6. Tittgemeyer M,
    7. Small DM
    (2018) Supra-additive effects of combining fat and carbohydrate on food reward. Cell Metab 28:33–44.e33. https://doi.org/10.1016/j.cmet.2018.05.018 pmid:29909968
    OpenUrlCrossRefPubMed
  18. ↵
    1. Dresselhuis DM,
    2. de Hoog EH,
    3. Stuart MA,
    4. van Aken GA
    (2008a) Application of oral tissue in tribological measurements in an emulsion perception context. Food Hydrocolloids 22:323–335. https://doi.org/10.1016/j.foodhyd.2006.12.008
    OpenUrl
  19. ↵
    1. Dresselhuis DM,
    2. de Hoog EH,
    3. Stuart MA,
    4. Vingerhoeds MH,
    5. van Aken GA
    (2008b) The occurrence of in-mouth coalescence of emulsion droplets in relation to perception of fat. Food Hydrocolloids 22:1170–1183. https://doi.org/10.1016/j.foodhyd.2007.06.013
    OpenUrlCrossRef
  20. ↵
    1. Eklund A,
    2. Nichols TE,
    3. Knutsson H
    (2016) Cluster failure: why fMRI inferences for spatial extent have inflated false-positive rates. Proc Natl Acad Sci USA 113:7900–7905. https://doi.org/10.1073/pnas.1602413113 pmid:27357684
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Eldeghaidy S,
    2. Marciani L,
    3. McGlone F,
    4. Hollowood T,
    5. Hort J,
    6. Head K,
    7. Taylor AJ,
    8. Busch J,
    9. Spiller RC,
    10. Gowland PA,
    11. Francis ST
    (2011) The cortical response to the oral perception of fat emulsions and the effect of taster status. J Neurophysiol 105:2572–2581. https://doi.org/10.1152/jn.00927.2010 pmid:21389303
    OpenUrlCrossRefPubMed
  22. ↵
    1. Eldeghaidy S,
    2. Marciani L,
    3. Hort J,
    4. Hollowood T,
    5. Singh G,
    6. Bush D,
    7. Foster T,
    8. Taylor AJ,
    9. Busch J,
    10. Spiller RC,
    11. Gowland PA,
    12. Francis ST
    (2016) Prior consumption of a fat meal in healthy adults modulates the brain's response to fat. J Nutr 146:2187–2198. https://doi.org/10.3945/jn.116.234104 pmid:27655761
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Friston KJ,
    2. Penny WD,
    3. Glaser DE
    (2005) Conjunction revisited. Neuroimage 25:661–667. https://doi.org/10.1016/j.neuroimage.2005.01.013 pmid:15808967
    OpenUrlCrossRefPubMed
  24. ↵
    1. Friston KJ,
    2. Holmes AP,
    3. Price CJ,
    4. Buchel C,
    5. Worsley KJ
    (1999) Multisubject fMRI studies and conjunction analyses. Neuroimage 10:385–396. https://doi.org/10.1006/nimg.1999.0484 pmid:10493897
    OpenUrlCrossRefPubMed
  25. ↵
    1. Grabenhorst F,
    2. Rolls ET
    (2008) Selective attention to affective value alters how the brain processes taste stimuli. Eur J Neurosci 27:723–729. https://doi.org/10.1111/j.1460-9568.2008.06033.x pmid:18279324
    OpenUrlCrossRefPubMed
  26. ↵
    1. Grabenhorst F,
    2. Rolls ET
    (2011) Value, pleasure, and choice in the ventral prefrontal cortex. Trends Cogn Sci 15:56–67. https://doi.org/10.1016/j.tics.2010.12.004 pmid:21216655
    OpenUrlCrossRefPubMed
  27. ↵
    1. Grabenhorst F,
    2. Rolls ET
    (2014) The representation of oral fat texture in the human somatosensory cortex. Hum Brain Mapp 35:2521–2530. https://doi.org/10.1002/hbm.22346 pmid:24038614
    OpenUrlCrossRefPubMed
  28. ↵
    1. Grabenhorst F,
    2. Rolls ET,
    3. Bilderbeck A
    (2008) How cognition modulates affective responses to taste and flavor: top down influences on the orbitofrontal and pregenual cingulate cortices. Cereb Cortex 18:1549–1559. https://doi.org/10.1093/cercor/bhm185 pmid:18056086
    OpenUrlCrossRefPubMed
  29. ↵
    1. Grabenhorst F,
    2. Rolls ET,
    3. Parris BA,
    4. D'Souza A
    (2010) How the brain represents the reward value of fat in the mouth. Cereb Cortex 20:1082–1091. https://doi.org/10.1093/cercor/bhp169 pmid:19684248
    OpenUrlCrossRefPubMed
  30. ↵
    1. Grabenhorst F,
    2. Schulte FP,
    3. Maderwald S,
    4. Brand M
    (2013) Food labels promote healthy choices by a decision bias in the amygdala. Neuroimage 74:152–163. https://doi.org/10.1016/j.neuroimage.2013.02.012 pmid:23428568
    OpenUrlCrossRefPubMed
  31. ↵
    1. Guest S,
    2. Grabenhorst F,
    3. Essick G,
    4. Chen Y,
    5. Young M,
    6. McGlone F,
    7. de Araujo I,
    8. Rolls ET
    (2007) Human cortical representation of oral temperature. Physiol Behav 92:975–984. https://doi.org/10.1016/j.physbeh.2007.07.004 pmid:17689575
    OpenUrlCrossRefPubMed
  32. ↵
    1. Hall PA,
    2. Best JR,
    3. Danckert J,
    4. Beaton EA,
    5. Lee JA
    (2023) Morphometry of the lateral orbitofrontal cortex is associated with eating dispositions in early adolescence: findings from a large population-based study. Soc Cogn Affect Neurosci 18:nsab084. https://doi.org/10.1093/scan/nsab084
    OpenUrl
  33. ↵
    1. Hare TA,
    2. Camerer CF,
    3. Rangel A
    (2009) Self-control in decision-making involves modulation of the vmPFC valuation system. Science 324:646–648. https://doi.org/10.1126/science.1168450 pmid:19407204
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Hare TA,
    2. Malmaud J,
    3. Rangel A
    (2011) Focusing attention on the health aspects of foods changes value signals in vmPFC and improves dietary choice. J Neurosci 31:11077–11087. https://doi.org/10.1523/JNEUROSCI.6383-10.2011 pmid:21795556
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Haynes JD
    (2015) A primer on pattern-based approaches to fMRI: principles, pitfalls, and perspectives. Neuron 87:257–270. https://doi.org/10.1016/j.neuron.2015.05.025 pmid:26182413
    OpenUrlCrossRefPubMed
  36. ↵
    1. Haynes JD,
    2. Sakai K,
    3. Rees G,
    4. Gilbert S,
    5. Frith C,
    6. Passingham RE
    (2007) Reading hidden intentions in the human brain. Curr Biol 17:323–328. https://doi.org/10.1016/j.cub.2006.11.072 pmid:17291759
    OpenUrlCrossRefPubMed
  37. ↵
    1. Hebart MN,
    2. Gorgen K,
    3. Haynes JD
    (2014) The Decoding Toolbox (TDT): a versatile software package for multivariate analyses of functional imaging data. Front Neuroinform 8:88. https://doi.org/10.3389/fninf.2014.00088 pmid:25610393
    OpenUrlCrossRefPubMed
  38. ↵
    1. Henson R
    (2003) Analysis of fMRI time series. In: Human brain function, Ed 2 (Ashburner J, Friston K, Penny W, eds). London: Academic.
  39. ↵
    1. Howard JD,
    2. Kahnt T
    (2017) Identity-specific reward representations in orbitofrontal cortex are modulated by selective devaluation. J Neurosci 37:2627–2638. https://doi.org/10.1523/JNEUROSCI.3473-16.2017 pmid:28159906
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Huang FY,
    2. Sutcliffe MP,
    3. Grabenhorst F
    (2021) Preferences for nutrients and sensory food qualities identify biological sources of economic values in monkeys. Proc Natl Acad Sci USA 118:e2101954118. https://doi.org/10.1073/pnas.2101954118
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Kadohisa M,
    2. Rolls ET,
    3. Verhagen JV
    (2005) Neuronal representations of stimuli in the mouth: the primate insular taste cortex, orbitofrontal cortex, and amygdala. Chem Senses 30:401–419. https://doi.org/10.1093/chemse/bji036 pmid:15829609
    OpenUrlCrossRefPubMed
  42. ↵
    1. Kahnt T,
    2. Heinzle J,
    3. Park SQ,
    4. Haynes JD
    (2011) Decoding different roles for vmPFC and dlPFC in multi-attribute decision making. Neuroimage 56:709–715. https://doi.org/10.1016/j.neuroimage.2010.05.058 pmid:20510371
    OpenUrlCrossRefPubMed
  43. ↵
    1. Koc H,
    2. Vinyard CJ,
    3. Essick GK,
    4. Foegeding EA
    (2013) Food oral processing: conversion of food structure to textural perception. Annu Rev Food Sci Technol 4:237–266. https://doi.org/10.1146/annurev-food-030212-182637 pmid:23244397
    OpenUrlPubMed
  44. ↵
    1. Kokini JL
    (1987) The physical basis of liquid food texture and texture-taste interactions. J Food Eng 6:51–81. https://doi.org/10.1016/0260-8774(87)90021-5
    OpenUrlCrossRef
  45. ↵
    1. Laguna L,
    2. Farrell G,
    3. Bryant M,
    4. Morina A,
    5. Sarkar A
    (2017) Relating rheology and tribology of commercial dairy colloids to sensory perception. Food Funct 8:563–573. https://doi.org/10.1039/c6fo01010e pmid:27722602
    OpenUrlPubMed
  46. ↵
    1. Lebreton M,
    2. Jorge S,
    3. Michel V,
    4. Thirion B,
    5. Pessiglione M
    (2009) An automatic valuation system in the human brain: evidence from functional neuroimaging. Neuron 64:431–439. https://doi.org/10.1016/j.neuron.2009.09.040 pmid:19914190
    OpenUrlCrossRefPubMed
  47. ↵
    1. Logothetis NK
    (2008) What we can do and what we cannot do with fMRI. Nature 453:869–878. https://doi.org/10.1038/nature06976 pmid:18548064
    OpenUrlCrossRefPubMed
  48. ↵
    1. Macey PM,
    2. Macey KE,
    3. Kumar R,
    4. Harper RM
    (2004) A method for removal of global effects from fMRI time series. Neuroimage 22:360–366. https://doi.org/10.1016/j.neuroimage.2003.12.042 pmid:15110027
    OpenUrlCrossRefPubMed
  49. ↵
    1. Mattes RD
    (2009) Is there a fatty acid taste? Annu Rev Nutr 29:305–327. https://doi.org/10.1146/annurev-nutr-080508-141108 pmid:19400700
    OpenUrlCrossRefPubMed
  50. ↵
    1. Mela DJ
    (1988) Sensory assessment of fat content in fluid dairy products. Appetite 10:37–44. https://doi.org/10.1016/s0195-6663(88)80031-x pmid:3355125
    OpenUrlCrossRefPubMed
  51. ↵
    1. Morys F,
    2. Shishikura M,
    3. Dagher A
    (2022) Population-based research in obesity: an overview of neuroimaging studies using big data approach. Curr Opin Endocr Metab Res 23:100323. https://doi.org/10.1016/j.coemr.2022.100323
    OpenUrl
  52. ↵
    1. Nichols T,
    2. Brett M,
    3. Andersson J,
    4. Wager T,
    5. Poline JB
    (2005) Valid conjunction inference with the minimum statistic. Neuroimage 25:653–660. https://doi.org/10.1016/j.neuroimage.2004.12.005 pmid:15808966
    OpenUrlCrossRefPubMed
  53. ↵
    1. Ogawa H,
    2. Ito S,
    3. Nomura T
    (1989) Oral cavity representation at the frontal operculum of macaque monkeys. Neurosci Res 6:283–298. https://doi.org/10.1016/0168-0102(89)90021-7 pmid:2725988
    OpenUrlCrossRefPubMed
  54. ↵
    1. Plassmann H,
    2. O'Doherty J,
    3. Rangel A
    (2007) Orbitofrontal cortex encodes willingness to pay in everyday economic transactions. J Neurosci 27:9984–9988. https://doi.org/10.1523/JNEUROSCI.2131-07.2007 pmid:17855612
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Rangel A
    (2013) Regulation of dietary choice by the decision-making circuitry. Nat Neurosci 16:1717–1724. https://doi.org/10.1038/nn.3561 pmid:24270272
    OpenUrlCrossRefPubMed
  56. ↵
    1. Rolls ET
    (2016) Reward systems in the brain and nutrition. Annu Rev Nutr 36:435–470. https://doi.org/10.1146/annurev-nutr-071715-050725 pmid:27146018
    OpenUrlCrossRefPubMed
  57. ↵
    1. Rolls ET
    (2020) The texture and taste of food in the brain. J Texture Stud 51:23–44. https://doi.org/10.1111/jtxs.12488 pmid:31598975
    OpenUrlPubMed
  58. ↵
    1. Rolls ET,
    2. Grabenhorst F
    (2008) The orbitofrontal cortex and beyond: from affect to decision-making. Prog Neurobiol 86:216–244. https://doi.org/10.1016/j.pneurobio.2008.09.001 pmid:18824074
    OpenUrlCrossRefPubMed
  59. ↵
    1. Rolls ET,
    2. Critchley HD,
    3. Browning AS,
    4. Hernadi A,
    5. Lenard L
    (1999) Responses to the sensory properties of fat of neurons in the primate orbitofrontal cortex. J Neurosci 19:1532–1540. https://doi.org/10.1523/JNEUROSCI.19-04-01532.1999 pmid:9952429
    OpenUrlAbstract/FREE Full Text
  60. ↵
    1. Rolls ET,
    2. Verhagen JV,
    3. Kadohisa M
    (2003) Representations of the texture of food in the primate orbitofrontal cortex: neurons responding to viscosity, grittiness and capsaicin. J Neurophysiol 90:3711–3724. https://doi.org/10.1152/jn.00515.2003 pmid:12917386
    OpenUrlCrossRefPubMed
  61. ↵
    1. Rolls ET,
    2. Mills T,
    3. Norton AB,
    4. Lazidis A,
    5. Norton IT
    (2018) The neuronal encoding of oral fat by the coefficient of sliding friction in the cerebral cortex and amygdala. Cereb Cortex 28:4080–4089. https://doi.org/10.1093/cercor/bhy213 pmid:30169795
    OpenUrlPubMed
  62. ↵
    1. Rolls ET,
    2. Feng R,
    3. Cheng W,
    4. Feng J
    (2023) Orbitofrontal cortex connectivity is associated with food reward and body weight in humans. Soc Cogn Affect Neurosci 18:nsab083. https://doi.org/10.1093/scan/nsab083
    OpenUrl
  63. ↵
    1. Rosenthal-von der Putten AM,
    2. Kramer NC,
    3. Maderwald S,
    4. Brand M,
    5. Grabenhorst F
    (2019) Neural mechanisms for accepting and rejecting artificial social partners in the uncanny valley. J Neurosci 39:6555–6570. https://doi.org/10.1523/JNEUROSCI.2956-18.2019 pmid:31263064
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Shewan HM,
    2. Pradal C,
    3. Stokes JR
    (2020) Tribology and its growing use toward the study of food oral processing and sensory perception. J Texture Stud 51:7–22. https://doi.org/10.1111/jtxs.12452 pmid:31149733
    OpenUrlPubMed
  65. ↵
    1. Shott ME,
    2. Cornier MA,
    3. Mittal VA,
    4. Pryor TL,
    5. Orr JM,
    6. Brown MS,
    7. Frank GK
    (2015) Orbitofrontal cortex volume and brain reward response in obesity. Int J Obes (Lond) 39:214–221. https://doi.org/10.1038/ijo.2014.121 pmid:25027223
    OpenUrlPubMed
  66. ↵
    1. Small DM,
    2. Voss J,
    3. Mak YE,
    4. Simmons KB,
    5. Parrish T,
    6. Gitelman D
    (2004) Experience-dependent neural integration of taste and smell in the human brain. J Neurophysiol 92:1892–1903. https://doi.org/10.1152/jn.00050.2004 pmid:15102894
    OpenUrlCrossRefPubMed
  67. ↵
    1. Small DM,
    2. Veldhuizen MG,
    3. Felsted J,
    4. Mak YE,
    5. McGlone F
    (2008) Separable substrates for anticipatory and consummatory food chemosensation. Neuron 57:786–797. https://doi.org/10.1016/j.neuron.2008.01.021 pmid:18341997
    OpenUrlCrossRefPubMed
  68. ↵
    1. Soltanahmadi S,
    2. Bryant M,
    3. Sarkar A
    (2023) Insights into the multiscale lubrication mechanism of edible phase change materials. ACS Appl Mater Interfaces 15:3699–3712. https://doi.org/10.1021/acsami.2c13017 pmid:36633252
    OpenUrlPubMed
  69. ↵
    1. Stice E,
    2. Yokum S,
    3. Burger KS,
    4. Epstein LH,
    5. Small DM
    (2011) Youth at risk for obesity show greater activation of striatal and somatosensory regions to food. J Neurosci 31:4360–4366. https://doi.org/10.1523/JNEUROSCI.6604-10.2011 pmid:21430137
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Stokes JR,
    2. Boehm MW,
    3. Baier SK
    (2013) Oral processing, texture and mouthfeel: from rheology to tribology and beyond. Curr Opin Colloid Interface Sci 18:349–359. https://doi.org/10.1016/j.cocis.2013.04.010
    OpenUrlCrossRef
  71. ↵
    1. Suzuki S,
    2. Cross L,
    3. O'Doherty JP
    (2017) Elucidating the underlying components of food valuation in the human orbitofrontal cortex. Nat Neurosci 20:1780–1786. https://doi.org/10.1038/s41593-017-0008-x pmid:29184201
    OpenUrlCrossRefPubMed
  72. ↵
    1. Szczesniak AS
    (1963) Classification of textural characteristics. J Food Sci 28:385–389. https://doi.org/10.1111/j.1365-2621.1963.tb00215.x
    OpenUrlCrossRef
  73. ↵
    1. Tang DW,
    2. Fellows LK,
    3. Dagher A
    (2014) Behavioral and neural valuation of foods is driven by implicit knowledge of caloric content. Psychol Sci 25:2168–2176. https://doi.org/10.1177/0956797614552081 pmid:25304885
    OpenUrlCrossRefPubMed
  74. ↵
    1. Toda T,
    2. Taoka M
    (2002) Hierarchical somesthetic processing of tongue inputs in the postcentral somatosensory cortex of conscious macaque monkeys. Exp Brain Res 147:243–251. https://doi.org/10.1007/s00221-002-1239-x pmid:12410339
    OpenUrlPubMed
  75. ↵
    1. Torres O,
    2. Andablo-Reyes E,
    3. Murray BS,
    4. Sarkar A
    (2018) Emulsion microgel particles as high-performance bio-lubricants. ACS Appl Mater Interfaces 10:26893–26905. https://doi.org/10.1021/acsami.8b07883 pmid:30036468
    OpenUrlPubMed
  76. ↵
    1. Tremblay L,
    2. Schultz W
    (1999) Relative reward preference in primate orbitofrontal cortex. Nature 398:704–708. https://doi.org/10.1038/19525 pmid:10227292
    OpenUrlCrossRefPubMed
  77. ↵
    1. Upadhyay R,
    2. Aktar T,
    3. Chen J
    (2020) Perception of creaminess in foods. J Texture Stud 51:375–388. https://doi.org/10.1111/jtxs.12509 pmid:32017109
    OpenUrlPubMed
  78. ↵
    1. van der Klaauw AA,
    2. Keogh JM,
    3. Henning E,
    4. Stephenson C,
    5. Kelway S,
    6. Trowse VM,
    7. Subramanian N,
    8. O'Rahilly S,
    9. Fletcher PC,
    10. Farooqi IS
    (2016) Divergent effects of central melanocortin signalling on fat and sucrose preference in humans. Nat Commun 7:13055. https://doi.org/10.1038/ncomms13055 pmid:27701398
    OpenUrlCrossRefPubMed
  79. ↵
    1. Verhagen JV,
    2. Rolls ET,
    3. Kadohisa M
    (2003) Neurons in the primate orbitofrontal cortex respond to fat texture independently of viscosity. J Neurophysiol 90:1514–1525. https://doi.org/10.1152/jn.00320.2003 pmid:12761278
    OpenUrlCrossRefPubMed
  80. ↵
    1. Wang XM,
    2. Chen JS
    (2017) Food oral processing: recent developments and challenges. Curr Opin Colloid Interface Sci 28:22–30. https://doi.org/10.1016/j.cocis.2017.01.001
    OpenUrl
  81. ↵
    1. Wood FW
    (1968) Psychophysical studies on the consistency of liquid foods. SCI Monogr 27:40–49.
    OpenUrl
  82. ↵
    1. Woolley JD,
    2. Gorno-Tempini ML,
    3. Seeley WW,
    4. Rankin K,
    5. Lee SS,
    6. Matthews BR,
    7. Miller BL
    (2007) Binge eating is associated with right orbitofrontal-insular-striatal atrophy in frontotemporal dementia. Neurology 69:1424–1433. https://doi.org/10.1212/01.wnl.0000277461.06713.23 pmid:17909155
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Yarmolinsky DA,
    2. Zuker CS,
    3. Ryba NJ
    (2009) Common sense about taste: from mammals to insects. Cell 139:234–244. https://doi.org/10.1016/j.cell.2009.10.001 pmid:19837029
    OpenUrlCrossRefPubMed
  84. ↵
    1. Zangemeister L,
    2. Grabenhorst F,
    3. Schultz W
    (2016) Neural basis for economic saving strategies in human amygdala-prefrontal reward circuits. Curr Biol 26:3004–3013. https://doi.org/10.1016/j.cub.2016.09.016 pmid:27773572
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 43 (47)
Journal of Neuroscience
Vol. 43, Issue 47
22 Nov 2023
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Masthead (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
A Neural Mechanism in the Human Orbitofrontal Cortex for Preferring High-Fat Foods Based on Oral Texture
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
A Neural Mechanism in the Human Orbitofrontal Cortex for Preferring High-Fat Foods Based on Oral Texture
Putu A. Khorisantono, Fei-Yang Huang (黃飛揚), Michael P. F. Sutcliffe, Paul C. Fletcher, I. Sadaf Farooqi, Fabian Grabenhorst
Journal of Neuroscience 22 November 2023, 43 (47) 8000-8017; DOI: 10.1523/JNEUROSCI.1473-23.2023

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
A Neural Mechanism in the Human Orbitofrontal Cortex for Preferring High-Fat Foods Based on Oral Texture
Putu A. Khorisantono, Fei-Yang Huang (黃飛揚), Michael P. F. Sutcliffe, Paul C. Fletcher, I. Sadaf Farooqi, Fabian Grabenhorst
Journal of Neuroscience 22 November 2023, 43 (47) 8000-8017; DOI: 10.1523/JNEUROSCI.1473-23.2023
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • dietary fat
  • oral food processing
  • neuroeconomics
  • preference
  • reward value

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • The Serotonergic Dorsal Raphe Promotes Emergence from Propofol Anesthesia in Zebrafish
  • Memory boost for recurring emotional events is driven by initial amygdala response promoting stable neocortical patterns across repetitions
  • Stimulus-dependent expression of Bdnf is mediated by ATF2, MYT1L, and EGR1 transcription factors
Show more Research Articles

Behavioral/Cognitive

  • Memory boost for recurring emotional events is driven by initial amygdala response promoting stable neocortical patterns across repetitions
  • Neurophysiology of perceptual decision-making and its alterations in attention-deficit hyperactivity disorder (ADHD)
  • Late Development of Sensory Thresholds for Horizontal Relative Disparity in Human Visual Cortex in the Face of Precocial Development of Thresholds for Absolute Disparity
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Notice
  • Contact
  • Accessibility
(JNeurosci logo)
(SfN logo)

Copyright © 2025 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.