Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Nutr., 07 January 2026

Sec. Nutrition, Psychology and Brain Health

Volume 12 - 2025 | https://doi.org/10.3389/fnut.2025.1699670

This article is part of the Research TopicEmerging Trends in Functional Food Formulation: From Concept to ConsumersView all articles

A low-calorie rebaudioside D–glucose–sodium sweetener matches sucrose in palatability and neural responses in a functional MRI study of healthy adults


Yuko Nakamura,*Yuko Nakamura1,2*Rika TakahashiRika Takahashi3Tadahiro OhkuriTadahiro Ohkuri3
  • 1Graduate School of Art and Sciences, Center for Evolutionary Cognitive Sciences, The University of Tokyo, Tokyo, Japan
  • 2University of Tokyo Institute for Diversity and Adaptation of Human Mind (UTIDAHM), Tokyo, Japan
  • 3Suntory Beverage and Food Limited, Kanagawa, Japan

Introduction: Non-nutritive sweeteners (NNSs) often taste less palatable than glucose-containing sugars, potentially because they differentially engage sweet taste pathways. In this study, we examined a low-calorie rebaudioside D–glucose–sodium (RebD mix) sweetener, formulated on the basis of previous evidence that its components can engage type 1 taste receptors 2 and 3 (T1r2/T1r3) and sodium–glucose cotransporter 1 (SGLT1)–related pathways. We tested whether the RebD mix matches sucrose in palatability and neural responses.

Methods: A within-subject fMRI study was performed on 28 healthy adults [mean age 22.36 ± standard deviation (SD) 5.20 years; male = 22: body mass index (BMI) kg/m2 = 20.34 ± 2.22 SD]. During scanning, participants received three gustatory conditions in randomized order, with tasteless rinses (9 s) between trials: sucrose (234 mM), RebD mix (29.2 mM sucrose + 194.3 mM glucose + 0.156 mM RebD + 5 mM sodium gluconate), and RebD mix without sodium (Na+) (29.2 mM sucrose + 194.3 mM glucose + 0.184 mM RebD), which is known to enhance SGLT1 activation. Liking and wanting ratings for gustatory stimuli were obtained at scanning. Imaging analyses focused on the orbitofrontal cortex (OFC), a key region associated with palatability, and OFC-centered functional connectivity, using within-participant contrasts across the three gustatory conditions.

Results: Liking and wanting did not differ across the three conditions. OFC response to the RebD mix did not differ from sucrose, whereas the RebD without Na+ condition showed greater OFC–postcentral gyrus connectivity than the other two conditions (psBonferroni–corrected < 0.001). Across participants, OFC–postcentral connectivity in the RebD without Na+ condition was negatively associated with fullness measured at scanning.

Conclusion: The RebD mix matched sucrose in subjective palatability and OFC responses, while removal of Na+ altered OFC– postcentral gyrus connectivity linked to satiety. These findings indicate that the RebD mix, a novel low-calorie formulation, may achieve sucrose-comparable palatability.

1 Introduction

The global prevalence of overweight and obesity has been rapidly increasing, with prevalence nearly doubling between 1990 and 2021 (1). The global increase in obesity and obese-related noncommunicable diseases, including metabolic disorders (e.g., type 2 diabetes mellitus), cardiovascular diseases, and some cancers, may be closely associated with the widespread consumption of sugar-sweetened beverages (SSBs), which are a major source of added sugars in the diet (2, 3). Given these public health concerns, non-nutritive sweeteners (NNSs) have emerged as popular alternatives to sugar, aiming to provide sweetness with minimal caloric impact. However, the low palatability of NNSs would impede their widespread utilization. For example, sucralose in Splenda® (Heartland Food Products Group, LLC., Indiana, USA), one of the most widely used NNSs, has a bitter, metallic, and chemical taste compared to sucrose (4). Furthermore, brain imaging studies have indicated that, compare to sugar-sweetened solutions, NNSs would elicit smaller responses in brain regions including the insula/operculum, orbitofrontal cortex (OFC) and striatum (5). Moreover, brain responses to a NNS solution that was matched to the sweetness level of sucrose differed significantly from brain responses to the sucrose solution (6). Therefore, even if the sweetness level of NNSs is matched to sucrose, these NNSs would elicit different brain responses in regions associated with palatability than sucrose, and thus would be less palatable.

More recently, in addition to artificial sweeteners such as the sucralose, naturally occurring NNSs like stevia has been introduced into the market (7). Steviol glycosides, which are the sweet-tasting compounds found in the plant stevia Rebaudiana, are perceived as 200–300 times sweeter than sucrose (8, 9). Steviol glycosides are often perceived by consumers as a healthier alternative due to its plant-derived origin (10). Moreover, in comparison with the 2-week glucose intake trial, sucralose or saccharin intake for 2 weeks impaired glucose tolerance, whereas stevia intake did not (11). Therefore, steviol glycosides are less likely than other artificial NNSs to disturb sugar metabolism.

Among the steviol glycosides, rebaudioside D (RebD) is more palatable due to its lower bitterness compared to other steviol glycosides, such as rebaudioside A (12, 13). Thus, RebD is a potential candidate for better substitute for sugar than other NNSs. Additionally, recent studies have revealed differences in sweet taste receptors between glucose-containing sugars, such as sucrose, and NNSs (14, 15). Therefore, adding components to RebD that stimulate glucose-specific receptors could improve the sweetness of the RebD mixture, making it comparable to sucrose. It is well established that sweet taste is mediated primarily by the type 1 taste receptors 2 and 3 (T1r2/T1r3) on taste cells, which responds to glucose-containing sugars such as sucrose and glucose as well as NNSs (16, 17). More recently, sodium–glucose cotransporter 1 (SGLT1), a metabolic sugar sensor in the intestine (18), has also been shown to mediate T1r2/T1r3-independent detection of glucose-containing sugars in the oral cavity (19). SGLT1-related sweet taste detection is enhanced by adding Na+ (20), and this beneficial effect is observed for glucose-containing sugars but not for NNSs (15). Thus, glucose-containing sugars can co-activate T1r2/T1r3 and SGLT1, whereas most NNSs primarily engage T1r2/T1r3 alone. These receptor-level differences may lead to distinct brain responses and lower palatability for NNSs, and suggest that SGLT1-related sweet taste signals from the oral cavity may contribute to the detection of nutritive sugars. These findings suggested that combining RebD with components that influence SGLT1-related signaling might improve its sweet taste to be more comparable to sucrose. We therefore previously invented a new sweetener by combining RebD, glucose, and Na+ (the RebD mix), which achieved moderate sweetness while keeping the energy content of the mixture low (13). RebD itself is a non-nutritive sweetener, but the RebD mix used in the present study combines RebD with sucrose and glucose and is therefore more appropriately described as a low-calorie sweetener rather than a strictly non-nutritive sweetener. Based on prior receptor assays showing that RebD can activate T1r2/T1r3 (12) and that glucose and Na+ modulate SGLT1-related signaling (12, 1517, 19, 20), the RebD mix was formulated with the aim of engaging both T1r2/T1r3- and SGLT1-related pathways, although receptor activation was not directly examined in the present study. However, the palatability of the RebD mix compared to sucrose has not been well examined. Additionally, it is unclear whether the brain response to the RebD mix would be similar to its response to glucose-containing sugars, such as sucrose.

In this study, we therefore aimed to compare the brain response to the RebD mix, the RebD mix without Na+, and sucrose, specifically in the OFC. The OFC is considered a secondary gustatory cortex (21), plays a broad role in processing gustatory pleasure—including hunger, attention, motivation, reward, and pleasantness (22)—and is strongly associated with the palatability of gustatory stimuli (23). It was hypothesized that the RebD mix would elicit OFC responses similar to those evoked by sucrose, whereas the RebD mix without Na+ might show different OFC responses because it lacks the Na+-dependent component that has been implicated in SGLT1-related signaling in previous studies (15, 20).

2 Materials and methods

2.1 Participants

A total of 28 participants were included in the current study [mean age 22.36 ± standard deviation (SD) 5.20 years; male = 22, female = 6; body mass index (BMI) kg/m2 = 20.34 ± 2.22 SD] (Table 1). Healthy individuals aged 18–45 years old are included. Individuals with a history of neurological injury; known genetic or medical disorders; food allergies; current or previous use of psychotropic medications; current or previous use of tobacco, cigarettes, or electronic cigarettes; or any exclusion criteria for MRI were excluded. The study procedures were approved by the Ethics Committee of the Department of Arts and Sciences at The University of Tokyo (Approval No. 955). Written informed consent was obtained from all participants.

TABLE 1
www.frontiersin.org

Table 1. Demographics of participants and ratings for internal state at the functional MRI (fMRI) session.

2.2 fMRI session

On the day of the fMRI session, participants were instructed to abstain from food or beverages except water for a minimum of 3 h prior to their visit to the laboratory. Subsequent to their arrival, the subjects were measured for height and weight. These measurements were used to calculate the subjects’ BMI. Approximately 30 min before the fMRI scanning, the participants were instructed to consume the predetermined snack (300 kcal), with the objective of standardizing their internal state (hunger and fullness sensations). Following consumption of the snack, participants undertook a brief version (two blocks) of the liquid consumption task for the fMRI scan. Subsequently, the participants were escorted to the MRI room. Immediately before the fMRI scan, participants were instructed to rate their internal state (hunger and fullness) using a visual analog scale (VAS) ranging from 0 to 100. The left end of the scale was labeled “not at all” and the right end was labeled “more than ever.” A paired t-test showed there was no significant difference between hunger and fullness ratings (p = 0.25). We thus assumed that participants were neither extremely hungry nor extremely full at the time of the scan. Additionally, after consuming a small amount of each gustatory solution for the liquid consumption task, they rated their liking (“How much do you like this juice?”), wanting (“How much do you want to drink this juice?”), and intensity (“How thick is this juice?”) using the VAS that used for their internal state ratings (Table 2). Subsequently, participants were instructed to complete four runs of the liquid consumption task. However, two participants completed only two runs. Thus, a total of 26 participants completed four runs. The average length of one run was 7 min 16.8 ± 6.5 s (mean ± SD) (Figure 1).

TABLE 2
www.frontiersin.org

Table 2. Ratings for liking, wanting, and intensity of each gustatory solution at the functional MRI (fMRI) session.

FIGURE 1
Functional MRI session timeline. It includes snack consumption, liquid task training, and ratings before four fMRI runs. Each run has ten blocks with breaks, featuring different liquid samples: sucrose, RebD mix, RebD mix without sodium, tasteless, and rinse.

Figure 1. The experimental procedures. Participants were instructed to abstain from food and beverages, except for water, for at least three hours before visiting the laboratory. Anthropometric assessments were performed. Then, approximately 30 min before the fMRI scan, participants were instructed to consume the predetermined 300-kcal snack. Subsequently, practice for the liquid consumption task was performed. Immediately before the scan, they were asked to rate their internal state, as well as their liking, wanting, and intensity for each stimulus. Participants were then instructed to complete four runs of the liquid consumption task while being scanned.

2.2.1 Solutions for the liquid consumption task at the fMRI session

For the liquid consumption task, sucrose solution (234 mM sucrose, 32 kcal/100 ml), the RebD mix solution (29.2 mM sucrose + 194.3 mM glucose + 0.156 mM RebD + 5 mM sodium gluconate, 18 kcal/100 ml), and the RebD mix without Na+ solution (29.2 mM sucrose + 194.3 mM glucose + 0.184 mM RebD, 18 kcal/100 ml) were used as gustatory solutions. During the fMRI task, the total energy delivered was 17.92 kcal for sucrose and 10.08 kcal for both RebD mix solutions. To match the intensity of sweetness, the RebD mix included a smaller amount of RebD compared to the RebD mix without Na+, as Na+ enhances sweetness when mixedwithglucose. Inadditiontothese solutions, a tasteless solution was used as a control. For the tasteless solution, first, four tasteless solutions were prepared, including the original tasteless solution (25 mM KCl + 2.5 mM NaHCO3) and solutions with 25%, 50%, and 75% of the original concentration. Participants tasted a small amount of each solution and selected the one that tasted “the most like nothing” as their tasteless solution.

2.2.2 Solution delivery

The timing and amount of solution delivered were controlled using programmable syringe pumps (NE-1000 One Channel Syringe Pump, NEW ERA PUMP SYSTEMS Inc., NY, USA). Details of the solution delivery system can be found elsewhere (24). In brief, through a clear polyvinyl chloride plastic tube, the syringe pumps delivered 0.80 mL of each solution over 2 s following a 1 s break.

2.2.3 The liquid consumption task

For the fMRI session, participants performed four runs of the liquid consumption task. In this task, gustatory and tasteless solutions were delivered to a participant through a mouthpiece placed between the participant’s upper and lower lips.

Each solution was delivered randomly in ten blocks over four runs. Two or three blocks of each solution were performed in random order within a run. Within a block, a solution was delivered six to eight times (18–24 s), followed by three rinses (9 s). A tasteless solution was used for rinse, and the next block commenced after a short break (14–16 s) (Figure 1).

2.2.4 MR image acquisition

All images were collected using a MAGNETOM Prisma 3.0 Tesla scanner with a 64-channel head coil (Siemens Healthineers, Erlangen, Germany). All high-resolution structural images of the brain were collected using a T1-weighted 3D MPRAGE protocol (repetition time (TR) = 1,900 ms, echo time (TE) = 2.53 ms, flip angle = 9°, field-of-view (FOV) = 256 × 256 mm2, resolution 1.0 × 1.0 × 1.0 mm3). As functional images, T2*-weighted images reflecting blood-oxygen-level-dependent (BOLD) signals were acquired using 2D gradient-echo echo-planar imaging (EPI) with an isotropic resolution (3.0 × 3.0 × 3.0 mm3), parallel acquisition factor = 3, TR = 2,000 ms, TE = 25 ms, 39 slices, and flip angle = 80° in a 192 mm2 FOV, transverse slices with phase encoding in the P > > A direction. All the images were acquired in an interleaved manner.

2.3 Statistical analysis

All functional and anatomical data underwent preprocessing using Statistical Parametric Mapping 12 (SPM12) (25) software with the CONN functional connectivity toolbox (CONN, version 22a) (26, 27). Functional and anatomical data were preprocessed using a modular preprocessing pipeline (28) including slice timing correction, creation of voxel-displacement maps, realignment with susceptibility distortion correction using fieldmaps, outlier detection, indirect segmentation and MNI-space normalization, and smoothing. Temporal misalignment between different slices of the functional data (acquired in interleaved bottom-up order) was corrected following SPM slice-timing correction (STC) procedure (29, 30), using sinc temporal interpolation to resample each slice BOLD timeseries to a common mid-acquisition time. Functional data were realigned using SPM realign and unwarp procedure (31) integrating fieldmaps for susceptibility distortion correction, where all scans were coregistered to a reference image (first scan of the first session) using a least squares approach and a six parameter (rigid body) transformation, and resampled using b-spline interpolation (32) to simultaneously correct for motion, magnetic susceptibility geometric distortions, and their interaction. Potential outlier scans were identified using ART (33) as acquisitions with framewise displacement above 0.9 mm or global BOLD signal changes above five standard deviations (34), and a reference BOLD image was computed for each subject by averaging all scans, excluding outliers. Functional and anatomical data were coregistered and normalized into standard MNI space, segmented into gray matter, white matter, and cerebrospinal fluid (CSF) tissue classes, and resampled to 2 mm isotropic voxels following an indirect normalization procedure (35) using SPM unified segmentation and normalization algorithm (36, 37) with the default IXI-549 tissue probability map template. Last, functional data were smoothed using spatial convolution with a Gaussian kernel of 6 mm full width half maximum (FWHM).

Preprocessed functional images from a total of 28 participants were brought to the individual level (first-level). The condition-specific effects (gustatory, tasteless, and rinse) at each voxel were estimated using a general linear model (GLM). The canonical hemodynamic response function provided by SPM12 was used to model the responses to the events. In the time-series analysis, a high-pass filter (270 s) was included in the filtering matrix to remove low-frequency noise and slow drifts from the signal. Confounding regressors from the preprocessing stage were also included in the model as covariates of no interest. Then, the [sucrose > tasteless], [RebD mix > tasteless], and [RebD mix without Na+ > tasteless] contrast images were created for individual participants.

Each contrast image was used in the subsequent group level (second-level) analysis. To examine the effect of each gustatory condition, a repeated measure analysis of variance (ANOVA) was performed using a voxel-wise GLM including ratings for intensity as a covariate-of-no-interest since there was a difference in ratings for intensity at a trend revel (p = 0.063) across three gustatory conditions. The predicted effect of these analyses was tested using the region of interest (ROI) approach. The ROI was created based on the Neurosynth (38) meta-analytic functional map for the term “taste” (80 studies, downloaded 15 February 2025). First, the downloaded brain map was thresholded at 60 % of probability (greater than 60 % of the peak z-value), then masked with the anatomical mask of the OFC crated with the automatic anatomical labeling (AAL) atlas (39) using the WFU Pickatlas toolbox (40, 41). Then, the sphere with 6 mm radius centered on the peak voxel [(x, y, z) = (−21.8, 33.5, −16.7)] of the masked brain map was created as the ROI in the OFC. The unpredicted effect of these analyses was tested using whole brain analysis. The threshold was set at pfamily–wise error rate (FWE)–corrected < 0.05 for the ROI and whole brain approaches, respectively.

Subsequently, whole-brain searchlight multivoxel pattern analysis (MVPA) (42) was employed to identify the brain regions that exhibited varied responses to each solution using the Decoding Toolbox (43) for 28 participants. This analysis was executed in accordance with a linear support vector machine approach. The training and testing of the analysis were conducted on individual-level regression coefficient maps derived from the individual level analysis, employing a leave-one-run-out cross-validation strategy. For this analysis, fMRI data was preprocessed using a modular preprocessing pipeline without MNI-space normalization and smoothing. Then, to generate the regression coefficient maps used as input for the MVPA, a separate individual-level regression model was applied to the fMRI data, which modeled each of the five conditions (sucrose, RebD mix, RebD mix without Na+, tasteless, and rinse) including confounding regressors. For each participant, a separate searchlight analysis was performed using a searchlight with a radius of three voxels, evaluating all three pairwise gustatory contrasts (sucrose vs. RebD mix, sucrose vs. RebD mix without Na+, RebD mix vs. RebD mix without Na+). The output of this searchlight analysis was a voxelwise map of average classification accuracy minus chance (33% for the multiclass comparisons). To evaluate the classification results at the group level, the resulting classification maps were warped to the standard MNI space, and applied a spatial smoothing (6 mm FWHM) to normalize the distribution of scores. Then, a group-level voxel-wise GLM was applied using permutation-based non-parametric testing with 5,000 permutations (44). Threshold-free cluster enhancement (TFCE) was employed to assess cluster significance (45), with threshold of pFWE–corrected < 0.05.

Finally, for seed to whole brain functional connectivity analysis, the generalized psychophysiological interaction (gPPI) (46, 47) analysis was adapted using preprocessed functional images. Since two participants completed only two runs, these two participants were excluded from this analysis, and a total of 26 participants were included in the gPPI. For this, the ROI used for the ANOVA was adapted as a seed. In addition, a sphere with a radius of 6 mm centered on the peak voxel [(x, y, z) = (34, 18, 6)] of the results from the whole-brain searchlight MVPA was used as a seed. At the individual level, separately, for each pair of seed and target areas, a gPPI was defined with BOLD time-series extracted from the seed as physiological factors, boxcar signals characterizing each individual gustatory or taste condition convolved with an SPM canonical hemodynamic response function as psychological factors, and the product of the two as psychophysiological interaction terms. Functional connectivity changes across conditions were characterized by the multivariate regression coefficient of the psychophysiological interaction terms in each model. At the group level, a GLM was estimated with first-level connectivity maps as dependent variables and conditions (sucrose, RebD mix, RebD mix without Na+) as independent variables to test for differences in connectivity across conditions using F-test. Inferences were drawn at the level of individual clusters. Cluster-level inferences were based on parametric statistics from Gaussian Random Field theory (48). Results were thresholded using a combination of a cluster-forming p < 0.001 voxel-level threshold, and a pfalse discovery rate (FDR)–corrected < 0.05 cluster-size threshold (49). Two seeds, the OFC and the insula, were used to create seed-to-whole brain connectivity maps separately. Thus, the significance level was set at pFDR–corrected < 0.025 (0.05/2).

Previous studies have shown that SGLT1-related sweet taste signals from the oral cavity may play a role in detecting sugars (15, 20), and sweet taste rapidly inhibits hunger-promoting agouti-related peptide-expressing (AgRP) neurons in the brain (50). Therefore, differences in connectivity among the three gustatory conditions could be related to internal states. Thus, the correlations between the connectivity that showed a significant difference across the three conditions and the internal states were compared between conditions. First, correlations between hunger or fullness and connectivity, controlling for BMI as a potential confounder were calculated. For this, partial correlations with the ppcor package (51) in R (version 4.4.1, R Foundation for Statistical Computing, Vienna, Austria) was used. Each internal state was tested three times. Thus, statistical significance was set at α = 0.0167 (0.05/3). Then, to compare the strength of partial correlations, Steiger’s test (52) for dependent correlations using the cocor package (53) was employed. This test is specifically designed to compare correlations that share a common variable and accounts for the intercorrelations between predictor variables. Each correlation coefficient was tested twice. Thus, statistical significance was set at α = 0.025 (0.05/2).

3 Results

3.1 Participant characteristics and behavioral ratings

Participant characteristics and behavioral ratings (hunger, fullness, liking, wanting, and intensity for each gustatory condition) are summarized in Tables 1, 2. Exploratory comparisons between males (n = 22) and females (n = 6) revealed no sex differences in age (p = 0.69) or in ratings for internal states (hunger and fullness) (ps > 0.28). BMI was higher in males than in females (p = 0.038). Ratings of liking, wanting, or intensity for any of the gustatory stimuli did not significantly differ between males and females (ps > 0.26).

The Friedman test showed that liking, wanting, and intensity ratings for all gustatory stimuli were not significantly different (ps > 0.063) (Table 2).

3.2 Brain response to each gustatory condition

There was significant averaged brain response to each gustatory contrast in the ROI [(x, y, z) = (−20, 36, −16), z = 3.07, pFWE–corrected = 0.028, cluster size = 3 voxels], while there was no significant difference in brain response to each gustatory contrast (Figure 2). Since the sex distribution was not balanced well in this study, sex was included as a covariate in additional analyses. There was still significant averaged brain response to each gustatory contrast in the ROI [(x, y, z) = (−20, 36, −16), z = 3.12, pFWE–corrected = 0.024, cluster size = 3 voxels]. The whole-brain analysis showed that there were no unexpected effects of gustatory contrasts.

FIGURE 2
MRI brain scan image next to a scatter plot. The MRI highlights a specific region with coordinates \([-20, 36, -16]\), marked with a blue circle and red dot, indicating activation with an F-value scale from 3 to 9. The scatter plot compares contrast estimates for Sucrose, RebD mix, and RebD mix without Na\(^+\), showing data points spread across a vertical axis from -3 to 2. Significance is noted as z = 3.07 with \(P_{FWE\text{-}corrected} = 0.028\).

Figure 2. Averaged brain response to each gustatory solution. The scatter plot indicates the left orbitofrontal cortex response to each gustatory solution. The color bar indicates F-value.

3.3 Results of whole-brain searchlight MVPA

The whole-brain searchlight MVPA identified the brain region that exhibited significant, above chance classification accuracy for discriminating between gustatory conditions in the anterior insula at a trend level [(x, y, z) = (−34, 18, 6), pFWE–corrected = 0.072, t = 4.68, an averaged accuracy at the peak voxel = 40.71%) (Figure 3).

FIGURE 3
Brain scan image showing activation in the insula, highlighted by a blue circle and red square. T-values range from one to four, color-coded from black to red. Coordinates are [-34, 18, 6], with FWE-corrected p-value of 0.072 and t-value of 4.68.

Figure 3. The brain region identified using a multivariate searchlight analysis, trained to distinguish between sucrose, the RebD mix, and the RebD mix without Na+ solutions. The color bar indicates T-value.

3.4 Connectivity analysis with the gPPI

Firstly, an examination was conducted to detect any differences in functional connectivity with the seed areas (OFC or insula) across the three conditions. Functional connectivity between the left OFC and the left postcentral gyrus significantly differs across three conditions [(x, y, z) = (−38, −24, 54), pFDR–corrected = 0.010, F(2, 50) > 7.96, cluster size = 62 voxels] (Figure 4A). As additional analysis, sex was included as a covariate. There was still significant difference in connectivity between the left OFC and the left postcentral gyrus across three conditions [(x, y, z) = (−38, −22, 50), pFDR–corrected = 0.0024, F(2, 48) > 8.00, cluster size = 78 voxels]. There were no differences in insula-to-whole-brain connectivity maps across the three conditions.

FIGURE 4
Brain imaging results and scatter plots illustrate connectivity and sense of fullness with different sweetener mixes. Image A shows the OFC connectivity with postcentral gyrus with F-values from 8.26 to 14.53. Plot displays connectivity for sucrose, RebD mix, and RebD mix without sodium. Image B features three scatter plots correlating connectivity with fullness for sucrose (r = -0.034), RebD mix (r = -0.220), and RebD mix without sodium (r = -0.628), showing a stronger negative correlation with connectivity in the RebD mix without sodium condition.

Figure 4. Results of the generalized psychophysiological interaction analysis. (A) The brain region that showed significantly different connectivity with the left orbitofrontal cortex across three distinct gustatory conditions. The scatter plot represents connectivity between the left orbitofrontal cortex and the postcentral gyrus in each gustatory condition. The color bar indicates F-value. (B) Associations between ratings for fullness and the left orbitofrontal cortex-the postcentral gyrus connectivity in each condition. The y-axis indicates ratings for fullness, and the x-axis indicates connectivity.

Post hoc analysis using extracted connectivity values showed that the connectivity between the OFC and left postcentral gyrus in the RebD mix without Na+ condition was greater than that in the RebD mix or sucrose conditions (psBonferroni–corrected < 0.001). The connectivity in the RebD mix condition did not significantly differ from that in the sucrose condition.

Partial correlations between fullness and connectivity measures, controlling for BMI, were as follows: connectivity in the sucrose condition (r = −0.034, p = 0.871), in the RebD mix condition (r = −0.220, p = 0.291), and in the RebD mix without Na+ (r = −0.628, p < 0.001). Thus, there was a significant negative association between connectivity in the RebD without Na+ and fullness. Hunger was not significantly correlated with any connectivity (p > 0.157). Steiger’s test revealed no significant difference in correlation strength between the sucrose and RebD mix conditions (z = 0.679, p = 0.500). However, there was a significant difference between the sucrose and RebD mix without Na+ conditions (z = 2.591, p = 0.001). Additionally, there was a trend toward a significant difference in correlation strength between the RebD mix and the RebD mix without Na+ conditions (z = 2.008, p = 0.045) (Figure 4B).

4 Discussion

The current study examined the brain response to the newly developed sweetener RebD mix, RebD mix without Na+, and sucrose in the OFC, which processes the palatability of gustatory stimuli (23). Additionally, the brain response to these sweeteners was examined throughout the entire brain, and the anterior insula was found to be associated with the decoding of these sweeteners. Functional connectivity analysis revealed significantly greater connectivity between the OFC and postcentral gyrus in the RebD without Na+ condition than in the other two gustatory conditions. Furthermore, this connectivity in the RebD without Na+ condition was negatively correlated with fullness measured at the fMRI scan.

There were no significant differences in OFC responses to each sweet solution, which corresponded to the ratings comparison results for liking, wanting, and intensity. However, there was a significant averaged brain response to each sweet solution in the OFC. Since liking or wanting ratings for all sweet solutions were nearly or greater than half of the rating range, those solutions would be equally palatable and represented by the OFC response, which could have a role in processing palatability of gustatory stimuli (23).

Multivoxel pattern analysis indicated that the anterior insula is related to discriminating between sweet taste solutions. This region of the insula is involved in taste detection (54) and taste discrimination (55). Thus, although participants were not instructed to distinguish between three different sweet solutions, differences in the composition of the sweet substances would automatically be reflected in the anterior insula, even if the participants did not recognize the differences.

Connectivity analysis revealed that connectivity between the OFC and postcentral gyrus was greater in the RebD mix without Na+ condition compared to other two gustatory conditions. The postcentral gyrus has been reported to be connected with impulsivity (56). In addition, the postcentral gyrus showed a significant response to marijuana cues in regular marijuana users who had abstained from use for 3 days (57). Furthermore, working for the preferred food resulted in a stronger response in the postcentral gyrus (58). Thus, connectivity between the OFC and the postcentral gyrus could be involved in cravings or motivation for food. Since the connectivity in the RebD mix without Na+ was greater than in the other two conditions, the SGLT1-related signal would reduce this connectivity, as the RebD mix and sucrose are expected to stimulate the SGLT1-related signal (15, 19). Additionally, this connectivity in the RebD mix without Na+ condition was negatively associated with fullness measured at the fMRI scan, while the connectivity in other two conditions were not significantly associated with fullness. Thus, sweet taste associated with SGLT1-related signals might mediate the OFC-postcentral gyrus connectivity to reduce cravings for sweet tastes.

Glucose-induced intestinal SGLT1 activation modulates the secretion of intestinal hormones that regulate glucose homeostasis, including glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic peptide (GIP) (5961), which are well known as satiety hormones (62). Therefore, sweet taste associated with SGLT1-related signals may be linked to feelings of satiety. In addition, sweet taste rapidly inhibits hunger-promoting AgRP neurons in the brain, and the reduction in AgRP neuron activity was much larger following consumption of nutritive foods (e.g., sucrose) than their non-nutritive counterparts (e.g., sucralose) (50). In sum, SGLT1-related sweet taste signals from the oral cavity may play a role in detecting nutritive sugars (15, 20), and glucose-induced intestinal SGLT1 activation modulates satiety sensation (5962), and nutritive sweet taste significantly reduces hunger-promoting neural activation (50). Therefore, SGLT1-related sweet taste could be related to satiety sensation, although there is no conclusive evidence indicating that the sweet taste sensation elicited by SGLT1 plays a role in satiety for sweet tastes. Unlike sucrose or other glucose-containing sugars, most NNSs primarily activate the T1r2/T1r3 sweet taste receptors and do not substantially engage SGLT1 (1517). Consistent with this distinct sensory profile, a recent survey reported that 63.6% of participants perceived NNSs as sweet but not sugar-like and indicated that they would prefer them to taste more like sugar (63). Taken together with the above findings on SGLT1-related signaling, these observations suggest that targeting SGLT1-mediated sweet taste—for example, by enhancing SGLT1 activation through the addition of Na+—may provide a novel approach to developing lower-calorie sweeteners with a more desirable, sugar-like taste. Such lower-calorie, sugar-like sweeteners could support the reformulation of food products and may be beneficial for diabetes management and the prevention of excessive weight gain.

The current study has several limitations. First, the sample size would be relatively small. A previous study that performed a similar analysis included only 20 participants (55), and the median sample size for fMRI experiments with gustatory stimuli was 24 (64). Thus, the sample size of the current study would be acceptable. However, future studies should have a greater sample size. Second, the current study did not examine the brain response to the detailed characteristics of the sweet taste quality of the RebD mix, such as aftertaste or sense of viscosity. Future studies should examine differences in brain responses to the detailed characteristics of sweet taste quality between sucrose and the RebD mix to further understand the neural mechanisms of RebD mix perception. Third, this study included more males (n = 22) than females (n = 6). Therefore, the current results may be affected by sex, although exploratory analyses did not reveal sex differences in behavioral ratings, and including sex as a covariate did not significantly alter the imaging or connectivity results. Future studies with larger and more sex-balanced samples could further clarify the palatability of the RebD mix and the associated brain responses to this novel sweetener. Fourth, this study did not systematically assess alcohol consumption. Although participants were instructed to abstain from food and beverages except water for at least 3 h before scanning, possible effects of habitual alcohol use on palatability or neural responses cannot be entirely ruled out. Finally, another limitation is that we did not assess blood glucose or other metabolic responses to the gustatory solutions. Given the relatively small total energy delivered during scanning and the randomized presentation of conditions, large systematic effects of acute glycemic changes on the imaging results are unlikely; however, we did not directly measure blood glucose. Thus, we cannot make conclusions about glycemic effects.

The current study examined brain response to the RebD mix, the RebD mix without Na+, and sucrose. There were no significant differences in the ratings of liking, wanting, or intensity for those solutions or the OFC responses to them. However, the anterior insula response pattern differed across three solutions, reflecting differences in compounds of each sweetener. Additionally, connectivity analysis indicated that the RebD mix without Na+, rather than the RebD mix or sucrose, would be associated with reduced satiety for sweet tastes, as the OFC-postcentral gyrus connectivity in the RebD mix without Na+ condition was negatively correlated with fullness measured at the fMRI scan. This result would suggest that SGLT1 activation, enhanced by Na+, is important for satiety in response to sweet tastes. Overall, the RebD mix is a less nutritious sweetener, yet it would be comparable to sucrose with respect to palatability and satisfaction.

Data availability statement

The original contributions presented in this study are included in this article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by the Ethics Committee of the Department of Arts and Sciences at The University of Tokyo. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

YN: Conceptualization, Writing – review & editing, Formal analysis, Visualization, Writing – original draft, Funding acquisition, Project administration. RT: Writing – review & editing. TO: Writing – review & editing, Writing – original draft.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This research was supported by the Division of the Research & Development Department of Suntory Beverage and Food Limited and by the Grant-in-Aid for Transformative Research Areas (A) of Japan Society for the Promotion of Science (JP21H05172).

Acknowledgments

We would like to express our sincere gratitude to all the participants and staff who were involved in this study. We appreciate their valuable collaboration and support throughout the research process.

Conflict of interest

This study was partially funded by Suntory Beverage and Food Limited, Kanagawa, Japan. The funder had the following involvement in the study: providing the sweeteners (the RebD mix and the RebD mix without Na+) tested in this study. RT and TO were employed by the company Suntory Beverage and Food Limited.

The remaining author declares that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

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

References

1. Sørensen TIA. Forecasting the global obesity epidemic through 2050. Lancet. (2025) 405:756–7. doi: 10.1016/S0140-6736(25)00260-0

PubMed Abstract | Crossref Full Text | Google Scholar

2. Malik VS, Hu FB. The role of sugar-sweetened beverages in the global epidemics of obesity and chronic diseases. Nat Rev Endocrinol. (2022) 18:205–18. doi: 10.1038/s41574-021-00627-6

PubMed Abstract | Crossref Full Text | Google Scholar

3. Nguyen M, Jarvis SE, Tinajero MG, Yu J, Chiavaroli L, Mejia SB, et al. Sugar-Sweetened beverage consumption and weight gain in children and adults: a systematic review and meta-analysis of prospective cohort studies and randomized controlled trials. Am J Clin Nutrit. (2023) 117:160–74. doi: 10.1016/j.ajcnut.2022.11.008

PubMed Abstract | Crossref Full Text | Google Scholar

4. Tan VWK, Wee MSM, Tomic O, Forde CG. Temporal sweetness and side tastes profiles of 16 sweeteners using temporal check-all-that-apply (Tcata). Food Res Intern. (2019) 121:39–47. doi: 10.1016/j.foodres.2019.03.019

PubMed Abstract | Crossref Full Text | Google Scholar

5. Yeung AWK, Wong NSM. How does our brain process sugars and non-nutritive sweeteners differently: a systematic review on functional magnetic resonance imaging studies. Nutrients. (2020) 12:3010. doi: 10.3390/nu12103010

PubMed Abstract | Crossref Full Text | Google Scholar

6. Frank GKW, Oberndorfer TA, Simmons AN, Paulus MP, Fudge JL, Yang TT, et al. Sucrose activates human taste pathways differently from artificial sweetener. NeuroImage. (2008) 39:1559–69. doi: 10.1016/j.neuroimage.2007.10.061

PubMed Abstract | Crossref Full Text | Google Scholar

7. Lohner S, Toews I, Meerpohl JJ. Health outcomes of non-nutritive sweeteners: analysis of the research landscape. Nutrit J. (2017) 16:55. doi: 10.1186/s12937-017-0278-x

PubMed Abstract | Crossref Full Text | Google Scholar

8. Ceunen S, Geuns JMC. Steviol glycosides: chemical diversity, metabolism, and function. J Nat Prod. (2013) 76:1201–28. doi: 10.1021/np400203b

PubMed Abstract | Crossref Full Text | Google Scholar

9. Lemus-Mondaca R, Vega-Gálvez A, Zura-Bravo L, Ah-Hen K. Stevia rebaudiana bertoni, source of a high-potency natural sweetener: a comprehensive review on the biochemical, nutritional and functional aspects. Food Chem. (2012) 132:1121–32. doi: 10.1016/j.foodchem.2011.11.140

PubMed Abstract | Crossref Full Text | Google Scholar

10. Prakash I, Dubois GE, Clos JF, Wilkens KL, Fosdick LE. Development of rebiana, a natural, non-caloric sweetener. Food Chem Toxicol. (2008) 46(Suppl 7):S75–82. doi: 10.1016/j.fct.2008.05.004

PubMed Abstract | Crossref Full Text | Google Scholar

11. Suez J, Cohen Y, Valdés-Mas R, Mor U, Dori-Bachash M, Federici S, et al. Personalized microbiome-driven effects of non-nutritive sweeteners on human glucose tolerance. Cell. (2022) 185: 3307–28.e19. doi: 10.1016/j.cell.2022.07.016.

PubMed Abstract | Crossref Full Text | Google Scholar

12. Hellfritsch C, Brockhoff A, Stähler F, Meyerhof W, Hofmann T. Human psychometric and taste receptor responses to steviol glycosides. J Agric Food Chem. (2012) 60:6782–93. doi: 10.1021/jf301297n

PubMed Abstract | Crossref Full Text | Google Scholar

13. Yokoo Y, Ohkuri T, Fujie A. Suntory Product Development Center, Assignee. Japan: Food and Drink with Increased Sweetness patent (2020).

Google Scholar

14. Vandenbeuch A, Kinnamon SC. Why low concentrations of salt enhance sweet taste. Acta Physiol. (2020) 230:e13560. doi: 10.1111/apha.13560

PubMed Abstract | Crossref Full Text | Google Scholar

15. Breslin PAS, Izumi A, Tharp A, Ohkuri T, Yokoo Y, Flammer LJ, et al. Evidence that human oral glucose detection involves a sweet taste pathway and a glucose transporter pathway. PLoS One. (2021) 16:e0256989. doi: 10.1371/journal.pone.0256989

PubMed Abstract | Crossref Full Text | Google Scholar

16. Nelson G, Hoon MA, Chandrashekar J, Zhang Y, Ryba NJ, Zuker CS. Mammalian sweet taste receptors. Cell. (2001) 106:381–90. doi: 10.1016/s0092-8674(01)00451-2

PubMed Abstract | Crossref Full Text | Google Scholar

17. Zhao GQ, Zhang Y, Hoon MA, Chandrashekar J, Erlenbach I, Ryba NJ, et al. The receptors for mammalian sweet and umami taste. Cell. (2003) 115:255–66. doi: 10.1016/s0092-8674(03)00844-4

PubMed Abstract | Crossref Full Text | Google Scholar

18. Wright EM, Loo DDF, Hirayama BA. Biology of human sodium glucose transporters. Physiol Rev. (2011) 91:733–94. doi: 10.1152/physrev.00055.2009

PubMed Abstract | Crossref Full Text | Google Scholar

19. Yee KK, Sukumaran SK, Kotha R, Gilbertson TA, Margolskee RF. Glucose transporters and ATP-gated K+ (KATP) metabolic sensors are present in type 1 taste receptor 3 (T1r3)-expressing taste cells. Proc Natl Acad Sci U S A. (2011) 108:5431–6. doi: 10.1073/pnas.1100495108

PubMed Abstract | Crossref Full Text | Google Scholar

20. Yasumatsu K, Ohkuri T, Yoshida R, Iwata S, Margolskee RF, Ninomiya Y. Sodium-Glucose cotransporter 1 as a sugar taste sensor in mouse tongue. Acta Physiol. (2020) 230:e13529. doi: 10.1111/apha.13529

PubMed Abstract | Crossref Full Text | Google Scholar

21. Rolls ET. The functions of the orbitofrontal cortex. Brain Cogn. (2004) 55:11–29. doi: 10.1016/s0278-2626(03)00277-x

PubMed Abstract | Crossref Full Text | Google Scholar

22. Kringelbach ML. The pleasure of food: underlying brain mechanisms of eating and other pleasures. Flavour. (2015) 4:20. doi: 10.1186/s13411-014-0029-2

Crossref Full Text | Google Scholar

23. Small DM, Bender G, Veldhuizen MG, Rudenga K, Nachtigal D, Felsted J. The role of the human orbitofrontal cortex in taste and flavor processing. Ann N Y Acad Sci. (2007) 1121:136–51. doi: 10.1196/annals.1401.002

PubMed Abstract | Crossref Full Text | Google Scholar

24. Nakamura Y, Ishida T. The effect of multiband sequences on statistical outcome measures in functional magnetic resonance imaging using a gustatory stimulus. NeuroImage. (2024) 300:120867. doi: 10.1016/j.neuroimage.2024.120867

PubMed Abstract | Crossref Full Text | Google Scholar

25. Penny WD, Friston KJ, Ashburner JT, Kiebel SJ, Nichols TE. Statistical Parametric Mapping: The Analysis of Functional Brain Images. Cambridge, MA: Academic Press (2011).

Google Scholar

26. Nieto-Castanon A, Whitfield-Gabrieli S. Conn Functional Connectivity Toolbox: Rrid Scr_009550, Release 22. Boston, MA: Hilbert Press (2022).

Google Scholar

27. Whitfield-Gabrieli S, Nieto-Castanon A. Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. (2012) 2:125–41. doi: 10.1089/brain.2012.0073

PubMed Abstract | Crossref Full Text | Google Scholar

28. Nieto-Castanon A. fMRI minimal preprocessing pipeline. Handbook of Functional Connectivity Magnetic Resonance Imaging Methods in Conn. Boston, MA: Hilbert Press (2020).

Google Scholar

29. Henson R, Buechel C, Josephs O, Friston KJ. The slice-timing problem in event-related fMRI. NeuroImage. (1999) 9:125.

Google Scholar

30. Sladky R, Friston KJ, Tröstl J, Cunnington R, Moser E, Windischberger C. Slice-Timing effects and their correction in functional MRI. Neuroimage. (2011) 58:588–94. doi: 10.1016/j.neuroimage.2011.06.078

PubMed Abstract | Crossref Full Text | Google Scholar

31. Andersson JL, Hutton C, Ashburner J, Turner R, Friston K. Modeling geometric deformations in EPI time series. Neuroimage. (2001) 13:903–19. doi: 10.1006/nimg.2001.0746

PubMed Abstract | Crossref Full Text | Google Scholar

32. Friston KJ, Ashburner J, Frith CD, Poline JB, Heather JD, Frackowiak RSJ. Spatial registration and normalization of images. Hum Brain Mapp. (1995) 3:165–89. doi: 10.1002/hbm.460030303

Crossref Full Text | Google Scholar

33. Whitfield-Gabrieli S, Nieto-Castanon A, Ghosh S. Artifact Detection Tools (ART), Release version 7:11. Cambridge, MA: Artifact Detection Tools (2011).

Google Scholar

34. Power JD, Mitra A, Laumann TO, Snyder AZ, Schlaggar BL, Petersen SE. Methods to detect, characterize, and remove motion artifact in resting state fMRI. Neuroimage. (2014) 84:320–41. doi: 10.1016/j.neuroimage.2013.08.048

PubMed Abstract | Crossref Full Text | Google Scholar

35. Calhoun VD, Wager TD, Krishnan A, Rosch KS, Seymour KE, Nebel MB, et al. The impact of T1 versus epi spatial normalization templates for fMRI data analyses. Hum Brain Mapp. (2017) 38:5331–42. doi: 10.1002/hbm.23737

PubMed Abstract | Crossref Full Text | Google Scholar

36. Ashburner J, Friston KJ. Unified segmentation. Neuroimage. (2005) 26:839–51. doi: 10.1016/j.neuroimage.2005.02.018

PubMed Abstract | Crossref Full Text | Google Scholar

37. Ashburner JA. Fast diffeomorphic image registration algorithm. Neuroimage. (2007) 38:95–113. doi: 10.1016/j.neuroimage.2007.07.007

PubMed Abstract | Crossref Full Text | Google Scholar

38. Yarkoni T, Poldrack RA, Nichols TE, Van Essen DC, Wager TD. Large-Scale automated synthesis of human functional neuroimaging data. Nat Methods. (2011) 8:665–70. doi: 10.1038/nmeth.1635

PubMed Abstract | Crossref Full Text | Google Scholar

39. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in spm using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. NeuroImage. (2002) 15:273–89. doi: 10.1006/nimg.2001.0978

PubMed Abstract | Crossref Full Text | Google Scholar

40. Maldjian JA, Laurienti PJ, Kraft RA, Burdette JH. An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. Neuroimage. (2003) 19:1233–9. doi: 10.1016/s1053-8119(03)00169-1

PubMed Abstract | Crossref Full Text | Google Scholar

41. Maldjian JA, Laurienti PJ, Burdette JH. Precentral gyrus discrepancy in electronic versions of the talairach atlas. Neuroimage. (2004) 21:450–5. doi: 10.1016/j.neuroimage.2003.09.032

PubMed Abstract | Crossref Full Text | Google Scholar

42. Kriegeskorte N, Goebel R, Bandettini P. Information-Based functional brain mapping. Proc Natl Acad Sci U S A. (2006) 103:3863–8. doi: 10.1073/pnas.0600244103

PubMed Abstract | Crossref Full Text | Google Scholar

43. Hebart MN, Görgen K, Haynes J-D. The Decoding Toolbox (TDT): a versatile software package for multivariate analyses of functional imaging data. Front Neuroinform. (2015) 8:88. doi: 10.3389/fninf.2014.00088

PubMed Abstract | Crossref Full Text | Google Scholar

44. Winkler AM, Ridgway GR, Webster MA, Smith SM, Nichols TE. Permutation inference for the general linear model. Neuroimage. (2014) 92(100):381–97. doi: 10.1016/j.neuroimage.2014.01.060

PubMed Abstract | Crossref Full Text | Google Scholar

45. Smith SM, Nichols TE. Threshold-Free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage. (2009) 44:83–98. doi: 10.1016/j.neuroimage.2008.03.061

PubMed Abstract | Crossref Full Text | Google Scholar

46. Friston KJ, Buechel C, Fink GR, Morris J, Rolls E, Dolan RJ. Psychophysiological and modulatory interactions in neuroimaging. Neuroimage. (1997) 6:218–29. doi: 10.1006/nimg.1997.0291

PubMed Abstract | Crossref Full Text | Google Scholar

47. McLaren DG, Ries ML, Xu G, Johnson SC. A Generalized form of context-dependent psychophysiological interactions (Gppi): a comparison to standard approaches. Neuroimage. (2012) 61:1277–86. doi: 10.1016/j.neuroimage.2012.03.068

PubMed Abstract | Crossref Full Text | Google Scholar

48. Worsley KJ, Marrett S, Neelin P, Vandal AC, Friston KJ, Evans AC. A unified statistical approach for determining significant signals in images of cerebral activation. Hum Brain Mapp. (1996) 4:58–73. doi: 10.1002/(sici)1097-019319964:13.0.Co;2-o

Crossref Full Text | Google Scholar

49. Chumbley J, Worsley K, Flandin G, Friston K. Topological FDR for neuroimaging. NeuroImage. (2010) 49:3057–64. doi: 10.1016/j.neuroimage.2009.10.090

PubMed Abstract | Crossref Full Text | Google Scholar

50. Aitken TJ, Liu Z, Ly T, Shehata S, Sivakumar N, La Santa Medina N, et al. Negative feedback control of hypothalamic feeding circuits by the taste of food. Neuron. (2024) 112:3354–70.e5. doi: 10.1016/j.neuron.2024.07.017.

PubMed Abstract | Crossref Full Text | Google Scholar

51. Kim S. Ppcor: an R package for a fast calculation to semi-partial correlation coefficients. Commun Stat Appl Methods. (2015) 22:665–74. doi: 10.5351/csam.2015.22.6.665

PubMed Abstract | Crossref Full Text | Google Scholar

52. Steiger JH. Tests for comparing elements of a correlation matrix. Psychol Bull. (1980) 87:245–51. doi: 10.1037/0033-2909.87.2.245

Crossref Full Text | Google Scholar

53. Diedenhofen B, Musch J. Cocor: a comprehensive solution for the statistical comparison of correlations. PLoS One. (2015) 10:e0121945. doi: 10.1371/journal.pone.0121945

PubMed Abstract | Crossref Full Text | Google Scholar

54. Veldhuizen MG, Gitelman DR, Small DM. An fMRI study of the interactions between the attention and the gustatory networks. Chemosens Percept. (2012) 5:117–27. doi: 10.1007/s12078-012-9122-z

PubMed Abstract | Crossref Full Text | Google Scholar

55. Avery JA, Liu AG, Ingeholm JE, Gotts SJ, Martin A. Viewing images of foods evokes taste quality-specific activity in gustatory insular cortex. Proc Natl Acad Sci U S A. (2021) 118:e2010932118. doi: 10.1073/pnas.2010932118

PubMed Abstract | Crossref Full Text | Google Scholar

56. Pan N, Wang S, Zhao Y, Lai H, Qin K, Li J, et al. Brain gray matter structures associated with trait impulsivity: a systematic review and voxel-based meta-analysis. Hum Brain Mapp. (2021) 42:2214–35. doi: 10.1002/hbm.25361

PubMed Abstract | Crossref Full Text | Google Scholar

57. Filbey FM, Schacht JP, Myers US, Chavez RS, Hutchison KE. Marijuana craving in the brain. Proc Natl Acad Sci U S A. (2009) 106:13016–21. doi: 10.1073/pnas.0903863106

PubMed Abstract | Crossref Full Text | Google Scholar

58. Wegman J, van Loon I, Smeets PAM, Cools R, Aarts E. Top-Down expectation effects of food labels on motivation. Neuroimage. (2018) 173:13–24. doi: 10.1016/j.neuroimage.2018.02.011

PubMed Abstract | Crossref Full Text | Google Scholar

59. Song P, Onishi A, Koepsell H, Vallon V. Sodium glucose cotransporter Sglt1 as a therapeutic target in diabetes mellitus. Exp Opin Ther Targets. (2016) 20:1109–25. doi: 10.1517/14728222.2016.1168808

PubMed Abstract | Crossref Full Text | Google Scholar

60. Gorboulev V, Schürmann A, Vallon V, Kipp H, Jaschke A, Klessen D, et al. Na(+)-D-Glucose Cotransporter Sglt1 is pivotal for intestinal glucose absorption and glucose-dependent incretin secretion. Diabetes. (2012) 61:187–96. doi: 10.2337/db11-1029

PubMed Abstract | Crossref Full Text | Google Scholar

61. Moriya R, Shirakura T, Ito J, Mashiko S, Seo T. Activation of sodium-glucose cotransporter 1 ameliorates hyperglycemia by mediating incretin secretion in mice. Am J Physiol Endocrinol Metab. (2009) 297:E1358–65. doi: 10.1152/ajpendo.00412.2009

PubMed Abstract | Crossref Full Text | Google Scholar

62. Astrup A. Reflections on the discovery Glp-1 as a satiety hormone: implications for obesity therapy and future directions. Eur J Clin Nutrit. (2024) 78:551–6. doi: 10.1038/s41430-024-01460-6

PubMed Abstract | Crossref Full Text | Google Scholar

63. Navanukraw C, Daosodsai P, Chanaboon S, Kanjanarach T. Exploring alternative sweeteners: consumer perspectives and product development for healthier and affordable choices. J Health Sci Med Res. (2024) 43:e20241132. doi: 10.31584/jhsmr.20241132

Crossref Full Text | Google Scholar

64. Yeung AWK, Wong NSM, Eickhoff SB. Empirical assessment of changing sample-characteristics in task- fMRI over two decades: an example from gustatory and food studies. Hum Brain Mapp. (2020) 41:2460–73. doi: 10.1002/hbm.24957

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: non-nutritive sweeteners, functional MRI, stevia, rebaudioside D, palatability

Citation: Nakamura Y, Takahashi R and Ohkuri T (2026) A low-calorie rebaudioside D–glucose–sodium sweetener matches sucrose in palatability and neural responses in a functional MRI study of healthy adults. Front. Nutr. 12:1699670. doi: 10.3389/fnut.2025.1699670

Received: 05 September 2025; Revised: 01 December 2025; Accepted: 05 December 2025;
Published: 07 January 2026.

Edited by:

Smriti Gaur, Jaypee Institute of Information Technology, India

Reviewed by:

Yugal Kishore Mohanta, University of Science and Technology, Meghalaya, India
Galileo Escobedo, General Hospital of Mexico, Mexico

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

*Correspondence: Yuko Nakamura, bmFrYW11cmEueXVrbzA3MDdAbWFpbC51LXRva3lvLmFjLmpw

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