Sensing in the Mouth: A Model for Filiform Papillae as Strain Amplifiers

1. Introduction

Food perception is based on the interaction between the human physiology and physical and chemical characteristics of food [1]. Molecular scale interactions between chemical compounds and receptors are usually responsible for taste, aroma and chemesthesis. Taste and chemesthesis pertain to the tongue [2–4] and the cells hosting taste receptors are located in onion-like structures called taste buds, located in circumvallate, foliate, and fungiform papillae. A fourth type of papilla, filiform papillae, are devoid of taste buds and are believed to be involved in texture perception, but the details of both the mechanical and molecular aspects of filiform papillae in encoding food texture are still unknown.

One important avenue of research consists in deciphering how the topological features of the tongue, mainly covered by filiform papillae, support the biological function of texture perception of foods—originally coined as psycho-rheology [5]. The particularly high sensitivity of the tongue in detecting small changes in stresses applied to the tissue remains to be conclusively linked to a physiological mechanism. Low yet perceptible stress changes can proceed from viscosity changes due to enzymatic [6], mechanical [7], or thermal breakdown of food structures, or due to the presence of micron-sized rigid particles in an otherwise homogeneous visco-elastic fluid such as glass beads in yogurt [8] or more universally a grain of sand in an oyster.

In this work, motivated by recent experiments reporting the presence of corpuscular endings in filiform papillae of mice tongues [9], we make the assumption that the origin of the high tactile acuity of the human tongue is the sum of mechanosensitive innervation [10], favorable transmission due to tissue mechanical properties and surface topology. Structures of similar aspect ratio to filiform papillae used to sense external fluid stresses have already been described in the animal kingdom. Theses include primary cilia at the cellular level [11] or superficial neuromasts in the fish lateral line [12, 13]. To function in a similar manner filiform papillae would need to bend significantly under typical in-mouth flows. Supporting the latter, we derive a coupled bio-fluid mechanical model of the mechanical advantage lying in the presence of filiform papillae.

This study is organized as follows. In Section 2 we introduce our modeling assumptions, including the geometry of papillae and mouth, their mechanical characteristics, the sensory function of the tongue and further physical assumptions. In Section 3 we then describe the flow as the tongue is moving and derive the scalings for the deformation strains of the papillae induced by the viscous fluid forces in Section 4. The scalings for the resulting elastic deformation of the tongue due to the bending of the clamped papillae are derived in Section 5 where we show that papillae are able to act as strain amplifiers by more than two orders of magnitude.

2. Model

2.1. Morphology

Over the last 60 years, a number of groups have investigated the morphology of filiform papillae (oftentimes motivated with associated diseases). Detailed studies exist on rats [14, 15], mice [16], porcupine [17], wild boar and pig [18], and humans [19, 20]. An overview of evolution in papilla morphology among vertebrates was detailed in Iwasaki [21]. An overview of the literature indicates the following morphological characteristics. In small vertebrates, individual papillae have a typical width ranging from 15 to 20 μm, and often taper conically at their tip. This width corresponds to the size of two to three epithelial cells, and is about the same order as the typical distance between papillae on the tongue. Their length can vary greatly between the anterior and posterior parts of the tongue and range from 50 to 100 μm. A precise set of measurements for humans were reported in Yamashita and OdDalkhsuren [22] where a distinction is made between the wide papillae body, which is anchored in the tongue, and the hairs on the papillae which do protrude upward from the tongue and deform under flow. Their size is reported to be 34 ± 16 μm in width and and 250 ± 62 μm long. Other measurements reports significantly larger sizes, with widths that can range 100–300 μm and length 200–500 μm [21, 23, 24]. Given the wide range of measured values among all the gathered papers, we will consider a half-width (radius) of 50 μm and length 250 μm in this paper as representative numbers.

2.2. Mechanics

No measurements of the Young's modulus of filiform papillae is available in the literature. A measure of E = 15 kPa of the pig's tongue modulus was obtain on fresh pig tongues [24]. That number in the range which would have been obtained by focusing solely on the elasticity of epithelial cells, for example E ≈ 1 − 10 kPa for monkey kidney epithelial cells [25] while E ≈ 0.1 − 300 kPa for human foreskin epithelial cells, with an average of 14 kPa for young cells and 33 kPa for old cells [26]. In light of these published data, we will use a value of E = 25 kPa in this study.

2.3. Mathematical Model

Based on the previous section, we can now put forward the model considered in this work, illustrated in Figure 1. We consider the fluid-solid interactions in a model mouth. The palate is modeled as a smooth, rigid, flat surface, at a distance H from the tongue. The tongue is assumed to be a smooth, rigid, flat surface on which filiform papillae are distributed. The total length of the flow region from the anterior to the posterior side of the mouth is denoted W. The focus of the work is on the deformation of the filiform papillae. Each papilla is modeled as a straight elastic rod of length L, radius a, clamped on the tongue. The Young's modulus of the rod is denoted E. As discussed above, we take E ≈ 25 kPa as a representative value. We also assume a = 50 μm, leading to a bending modulus B = πa⁴E/4 ≈ 1.2 × 10⁻¹³ J.m.

FIGURE 1

Figure 1. (A) Schematic representation of mouth model used in this work. The palate and the tongue are assumed to be smooth, rigid, flat surfaces. The tongue is located at a distance H from the palate and the total width of the fluid region in the anterior-posterior direction is denoted W; (B) Schematic representation of papilla model: a straight elastic rod of length L, radius a, clamped on the tongue.

2.4. Sensory Function of the Tongue

As stated in the introduction, the bovine tongue sensory innervation has already been partially described [10]. In addition to this, a recent study carried out in mice reported the presence of corpuscular endings innervating each filiform papillae [9], confirming a similar finding from previous work also describing innervation of filiform papillae in cats [27]. Since those endings do not project to the surface of the oral epithelium, we can hypothesize that they are mechanosensitive in nature and do not pertain to chemical stimulation. Our motivation to study the deformation mechanics of filiform papillae under flow lies in the fact that applied macroscopic strains will cascade into microscopic cell and transmembrane protein deformations and lead to somatosensory sensing as already known in the mamallian skin [28]. This should thus open a window for biophysicists regarding the levels of stresses that are likely to be applied to mechanosensory cells when mammals probe food for its texture.

3. Flow Induced by the Motion of the Tongue

3.1. Parameters

As we are interested in the interactions between deforming papillae and the flow of the food product, it is important to accurately model the driving of the flow between the palate and the tongue. During tasting and swallowing the palate is stationary and the tongue moves horizontally parallel to and vertically toward the palate driven by the action of lingual muscles. The typical setup is therefore that of (1) a squeeze flow driven by the vertical motion of the bottom surface (so-called transverse motion in the literature [29]) and (2) a shear component also present during tasting (referred to as longitudinal motion [29]). The shear component of the tongue movement can play an important role when acting alone, but as we see below the stresses it creates can safely be neglected when a transverse flow is present.

The study is motivated by the measurements of textures on complex fluids such as ice cream. Rheological measurements for 38% fat cream show a viscosity of about 20 mPa.s, hence approximately 20 times that of water [30]. Ice-creams typically have a viscosity in the range 20–300 times that of water [31], with strong variation with temperature and fat content. Recent measurements for ice cream with fat, fat replacers, and sweeteners indeed confirm this showing viscosities O(100 mPa.s) [32, 33]. In this study we pick a reference viscosity of μ = 100 mPa.s, equivalent to 100 times that of water at room tenperature.

It is also necessary to account for typical velocities in mouth during food consumption. Detailed measurements of movement in the mouth during swallowing in show velocities ranging from 1 to 15 cm/s, with an average around 10 cm/s Shawker et al. [29]. This is consistent with results from computational simulations [34]. For a tongue-palate distance of H ≈ 5 mm, a papilla size L ≈ 250 μm, and a viscosity μ = 100 mPa.s, this leads to Reynolds numbers O(1 − 10) in the mouth while the flows around the papillae are characterized by Re ~ 10⁻¹ or even smaller due to the no-slip boundary condition near the papillae. Given that the relevant Reynolds number in narrow geometries is the regular Reynolds number times the channel aspect ratio, we are in the low-Reynolds number regime and can neglect the influence of inertia on the dynamics of the fluid in the mouth [35, 36].

3.2. Newtonian Fluid Flow

For the flow in the mouth we consider the basic setup illustrated in Figure 2. The tongue is assumed to move both in the longitudinal (horizontal) direction, with a typical velocity U, and in the transverse (vertical) direct with typical velocity V. We use x and y to denote the directions along and perpendicular to the tongue (see Figure 2), with the tongue located at y = 0 and the palate at y = H.

FIGURE 2

Figure 2. Setup for the flow induced between the tongue and the palate. We denote by U the longitudinal velocity of the tongue (leading a shear flow) and V its vertical upward velocity (leading to a squeeze, pressure-driven flow).

The longitudinal motion leads to a shear flow in the fluid. The solution of Stokes equations in this case is simply u = u_∥e_x with a linear profile [37]

$\begin{array}{ll}{u}_{\parallel }=U\frac{H-y}{H},& (1)\end{array}$

where we use the subscript ∥ to indicate that the flow is driven by motion parallel to the tongue. The typical shear rate acting on the tongue in this case is on the order of ${\dot{\gamma }}_{\parallel }=O(U/H)$.

The transverse motion leads to a pressure-driven (squeeze) flow in the thin gap between the tongue and the palate. Assuming unidirectional flow, we have in this case u = u_⊥e_x with u_⊥ satisfying

$\begin{array}{ll}\mu \frac{{\text{d}}^2{u}_{\perp }}{\text{d}{y}^2}=\frac{\text{d}p}{\text{d}x},& (2)\end{array}$

and where we use the subscript ⊥ to indicate that the flow is driven by motion perpendicular to the tongue. The solution to Equation (2) is the standard parabolic flow

$\begin{array}{ll}{u}_{\perp }(y)=\frac{1}{2\mu }\frac{\text{d}p}{\text{d}x}y(y-H).& (3)\end{array}$

The pressure gradient is found using conservation of mass. Indeed, the flow rate along the tongue, Q, needs to the equal to the flow rate induced by the motion of the tongue, and so Q = VW. Integrating Equation (3) across the gap we obtain the scaling

$\begin{array}{ll}Q~-\frac{H^3}{\mu }\frac{\text{d}p}{\text{d}x},& (4)\end{array}$

and therefore

$\begin{array}{ll}\frac{1}{\mu }\frac{\text{d}p}{\text{d}x}~-\frac{VW}{H^3},& (5)\end{array}$

leading to a flow profile scaling as

$\begin{array}{ll}{u}_{\perp }(y)~V\frac{Wy(H-y)}{H^3}·& (6)\end{array}$

The typical shear rate on the tongue in this case is therefore given by ${\dot{\gamma }}_{\perp }=O(VW/H^2)$ and the ratio between the typical longitudinal and transverse shear rates is thus given by

$\begin{array}{ll}\frac{{\dot{\gamma }}_{\parallel }}{{\dot{\gamma }}_{\perp }}~\frac{U/H}{VW/H^2}~\left(\frac{U}{V}\right)\left(\frac{H}{W}\right).& (7)\end{array}$

Detailed measurements on tongue movement during swallowing shows that U is usually inferior to V and, in the rare cases where it is above, at most a factor of 2 larger [29]. Given that we are in a very high aspect-ratio geometry, H ≪ W by about one order of magnitude, and thus we obtain the result that ${\dot{\gamma }}_{\parallel }\ll {\dot{\gamma }}_{\perp }$. The shear flow induced by the tongue in the anterior-posterior direction leads to stresses on the tongue which are smaller in comparison with the stresses induced by the tongue motion in the superior/inferior direction and the resulting squeeze flow. In summary, when a squeeze flow is present, the shear flow is small and can be safely ignored.

4. Deformation of Single Papilla into Newtonian Flow

4.1. Dynamical Regime

With the Newtonian flow in the tongue characterized in Equations (1) and (6), we now consider its effect on a single papilla. Before computing the shape, and the strain, of a deforming papilla, a number of questions need to be addressed in order to properly characterize their dynamic regime.

First, do papilla deform with unsteady dynamics (transient motion and relaxation) or are the time scales involved sufficiently long that we can model them as deforming in a quasi-steady fashion? This is a question in the realm of fluid-structure interactions of filaments in viscous fluids [38]. With a bending modulus B and a papilla of length L in a viscous fluid of viscosity μ, the typical time scale for the deformation to reach its steady-state shape is given by an elasto-hydrodynamics time scale, t_eh, [39]

$\begin{array}{ll}{t}_{eh}~\frac{\xi L^4}{B},& (8)\end{array}$

where ξ is the drag coefficient for flow near the papilla, ξ ≈ 4πμ/[ln(2L/a) + 1/2]. With μ = 100 mPa.s, L ~ 250 μm, and B = 1.2 × 10⁻¹³ J.m, we obtain a typical time scale of t_eh = 0(10 ms). How does this compare with the typical time scale, T, for the flow in the mouth? This timescale is given by H/V, it is the time scale over which the flow in the mouth is going to change. The vertical velocity is on the order of V ≈ 10 cm/s while the thickness is H ≈ 5 mm and thus T = H/V ≈ 50 ms. Since the high value for V can be considered an upper bound, the time scale we obtained of 50 ms should be viewed as a lower bound. Given that T ≳ t_eh we are thus in a regime where we expect unsteadiness to not play an important role and the deformation will be treated as quasi-static.

The second important question about the dynamics concerns the issue of small vs. large deformations. Do we expect the papillae to deform in the nonlinearly geometric regime where internal tension plays an important role and careful attention needs to be paid to the extensibility of the papilla, or do we remain safely in the linear regime? If we were to remain in the linear regime, the typical deflection of a linear beam would be scaling in the following manner. Let us call δ the typical magnitude of the papilla tip deflection. For a force of magnitude q per unit length, δ scales as δ ~ qL⁴/B [40]. Here the load is due to the fluid drag and thus q ~ ξu where we use u to denote the typical velocity of the fluid relative to that of the papilla, leading to $\delta ~\xi u{L}^4/B=t_{eh}u$.

We need to then consider separately the shear and squeeze flows. In the case of a shear flow, the typical relative velocity around the papilla is u_∥ ≈ UL/H leading to the scaling δ/L ~ t_ehU/H. With t_eh = 10 ms, H = 5 mm, and U ≲ 10 cm/s we obtain δ/L = 0(10⁻¹) and remain safely in the linear regime. In the case of a squeeze flow we have the scaling near the papilla u ~ VLW/H² and thus we obtain $\delta /L~t_{eh}VW/H^2$. Here again with t_eh = 10 ms, H = 5 mm, V ≲ 10 cm/s, and W = 5 cm, we obtain δ/L = 0(1), which is at the limit but is a regime in which the linearized equations of solid mechanics will be at least qualitatively accurate, and approximately quantitatively accurate.

4.2. Mathematical Model

In order to describe the deformation of a papilla in a flow, and the elastic strain that it applies to the tongue, we use the formulation for the dynamics of elastic filaments in viscous flows. Let us denote by r the location of a point along the filament and s the arc length along the filament (from s = 0 to s = L). We use u to denote the external flow field located at position r and τ the tension (force per unit length) enforcing inextensibility of the filament. The instantaneous balance of forces and moments on the filament leads then to [38]

$\begin{array}{ll}\xi ({\text{r}}_t-\text{u})=-B{\text{r}}_{ssss}+{(\tau {\text{r}}_s)}_s.& (9)\end{array}$

In Equation (9) we have used the shorthand notation that (…)_s means a partial derivative along the s direction and (…)_t a partial time derivative (so that the first term on the right-hand side of Equation (9) has four derivatives along s).

To simplify the mathematical approach and the interpretation of the results we proceed to non-dimensionalize the problem. We use the length of the papilla, L, as the characteristic length, t_eh as the characteristic time (Equation 8), L/t_eh as the characteristic velocity and B/L² as characteristic tension. Doing so, Equation (9) becomes

$\begin{array}{ll}{\text{r}}_t=\text{u}-{\text{r}}_{ssss}+{(\tau {\text{r}}_s)}_s.& (10)\end{array}$

The equation for the tension τ is found by enforcing that r_s · r_s remains equal to one for all times, leading to

$\begin{array}{ll}{\tau }_{ss}-\tau ({\text{r}}_{ss}·{\text{r}}_{ss})+4{\text{r}}_{ss}·{\text{r}}_{ssss}+3{\text{r}}_{sss}·{\text{r}}_{sss}+{\text{r}}_s·{\text{u}}_s=0.& (11)\end{array}$

In Equation (12) the term including u_s is small and can be neglected since it quantifies gradients of the flow along the length of the papilla. In the linear regime, r_s = t = e_y is the tangent to the papilla while u_s = ∂u/∂y is the y derivative of the flow, which is nonzero along the x direction only, and therefore r_s · u_s = 0. We thus obtain the further simplification

$\begin{array}{ll}{\tau }_{ss}-\tau ({\text{r}}_{ss}·{\text{r}}_{ss})+4{\text{r}}_{ss}·{\text{r}}_{ssss}+3{\text{r}}_{sss}·{\text{r}}_{sss}=0.& (12)\end{array}$

Furthermore, using our scaling approach from the previous section we can neglect any nonlinear terms in the shape of the papillae. As can be seen in Equation (12), the tension τ is expected to be quadratic in the filament amplitude and therefore the last term in Equation (10) is cubic, and can be neglected. We thus end up with the linear beam equation

$\begin{array}{ll}\text{u}-{\text{r}}_t={\text{r}}_{ssss},& (13)\end{array}$

which, when simplified further due to the fast relaxation of unsteady modes (see Section 4.1), becomes the quasi-steady equation

$\begin{array}{ll}\text{u}={\text{r}}_{ssss},& (14)\end{array}$

Physically, Equation (14) indicates a balance between bending and the drag forces from the fluid.

4.3. Scaling Laws for Deflections

Calling δ the typical deflection of an individual papilla and coming back to dimensional equations, we thus obtain the linearized balance

$\begin{array}{ll}\xi u~B\frac{{\text{d}}^4\delta }{\text{d}{y}^4},& (15)\end{array}$

leading to a typical tip deflection on the order of

$\begin{array}{ll}\delta ~\frac{\xi L^{4}u}{B}·& (16)\end{array}$

We can then use the two velocity profiles, Equation (1) and Equation (6), to obtain the typical order of magnitude of the deflections.

4.3.1. Shear Flow

For a shear flow we have u = u_∥ ~ UL/H leding to a typical deflection

$\begin{array}{ll}{\delta }_{\parallel }~\frac{\xi U{L}^5}{BH}·& (17)\end{array}$

The longitudinal strain, ϵ, along the papilla is given by ϵ = aκ where a is the radius of the papilla and κ the curvature of the papilla, κ ≈ δ″(y) ~ δ/L². The typical longitudinal strain, ϵ_∥, scales thus as

$\begin{array}{ll}{ϵ}_{\parallel }~\frac{a\xi U{L}^3}{BH}·& (18)\end{array}$

Using the parameters of our model, μ = 100 mPa.s, ξ ≈ 4πμ/[ln(2L/a) + 1/2] ≈ 0.45 Pa.s, L ~ 250 μm, B = 1.2 × 10⁻¹³ J.m, U ≈ 10 cm/s, H ≈ 5 mm and a ≈ 50 μm, we obtain small strains of order |ϵ_∥| ≈ 6%.

4.3.2. Squeeze Flow

In the case of the squeeze flow the velocity scales as

$\begin{array}{ll}{u}_{\perp }(y)~V\frac{Wy(H-y)}{H^3}·& (19)\end{array}$

Furthermore, since the papillae are much smaller than the mouth, L ≪ H, the velocity profile close to the tongue is approximately a shear flow and the papilla is not able to feel the curvature of the velocity profile in the center of the mouth. The typical magnitude of the flow near a papilla is therefore given by

$\begin{array}{ll}{u}_{\perp }(y)~V\frac{WL}{H^2}·& (20)\end{array}$

Using Equation (16), we thus obtain a typical deflection

$\begin{array}{ll}{\delta }_{\perp }~\frac{\xi VW{L}^5}{B{H}^2},& (21)\end{array}$

and therefore a longitudinal strain in the papilla scaling as

$\begin{array}{ll}{ϵ}_{\perp }~\frac{a\xi VW{L}^3}{B{H}^2}·& (22)\end{array}$

Using the same numbers as above with in addition V ≈ 10 cm/s and W = 5 cm, leads to a much larger strain than in the longitudinal case, |ϵ_⊥| ≈ 60%, indicating significant strains. Although a more complete quantitative model would include all nonlinear terms in the beam equation, this results provides a qualitative estimate.

5. Papillae-Induced Substrate Deformation

5.1. Main Question

We are now in a position to answer the main question of this work: do the papillae, through their deformation and the fact that they are anchored to the tongue, lead to deformation of the tongue which is stronger than if the papillae were not present. In other words, if the mechano-sensors are distributed beneath the surface of the tongue, could papillae amplify deformation of the tongue mediated through fluid-structure interactions? We can quantitatively address this question within our mathematical framework. In addition, since sensory neurons have also been shown to be mechanosensitive themselves (and thus not solely serving conduction of surface mechanosensory cells [41]), it is important to ask whether such neurons could be impacted by papillae strains.

5.2. Fundamental Solution: Point-Force on an Elastic Substrate

The first step is to consider the Green's function for the deformation of a semi-infinite elastic substrate due to a point force. With the Green's function known, all other types of deformations can be tackled mathematically.

Consider an elastic medium located in the x₃ < 0 semi-infinite plane, as a model for the tongue. The Young's modulus of the tongue is denoted E and its Poisson's ratio ν. A constant force of magnitude F₃ is applied at the surface of the tongue, with the surface being parallel to the (x₁, x₂) plane. The force is being applied at the origin of the coordinate system. The solution to this Green's function problem is classical [40] and we summarize it here. The elastic displacements in the tongue along each direction, {u₁, u₂, u₃} at position {x₁, x₂, x₃} are given by

$\begin{array}{ll}{u}_1=\frac{1+\sigma }{2\pi E}\left[-\frac{x_1{x}_3}{r^3}-\frac{(1-2\sigma )x_1}{r(r-x_3)}\right]F_3,& (23\text{a})\end{array}$

$\begin{array}{ll}{u}_2=\frac{1+\sigma }{2\pi E}\left[-\frac{x_2{x}_3}{r^3}-\frac{(1-2\sigma )x_2}{r(r-x_3)}\right]F_3,& (23\text{b})\end{array}$

$\begin{array}{ll}{u}_3=\frac{1+\sigma }{2\pi E}\left[\frac{x_3^2}{r^3}+\frac{2(1-\sigma )}{r}\right]F_3,& (23\text{c})\end{array}$

where $r^2=\sum _i{x}_i^2$.

For simplicity we then assume that the tongue is incompressible, meaning its Poisson's ratio is given by σ = 1/2. This simplifies the results for the displacements as

$\begin{array}{ll}{u}_1=\frac{3}{4\pi E}\left(-\frac{x_1{x}_3}{r^3}\right)F_3,& (24\text{a})\end{array}$

$\begin{array}{ll}{u}_2=\frac{3}{4\pi E}\left(-\frac{x_2{x}_3}{r^3}\right)F_3,& (24\text{b})\end{array}$

$\begin{array}{ll}{u}_3=\frac{3}{4\pi E}\left(\frac{x_3^2}{r^3}+\frac{1}{r}\right)F_3.& (24\text{c})\end{array}$

5.3. Strain Field

With the elastic displacements known, we can then compute the strain field, and it is given by the symmetric part of the displacement gradients. As can be seen in Equation (24), the typical displacement inside the tongue scales as u ~ F/Er at a distance r from the point where the force F is being applied. The strains in the tongue, ε, are therefore expected to take the approximate values

$\begin{array}{ll}\epsilon ~\frac{u}{r}~\frac{F}{E{r}^2}·& (25)\end{array}$

5.4. Point moment on an Elastic Substrate

The deforming papilla is not applying a force on the tongue but instead it applies a moment at its base. If we denote by M the applied moment then we need to take one additional spatial derivative to obtain the strains in the tongue, leading to

$\begin{array}{ll}\epsilon ~\frac{M}{E{r}^3}·& (26)\end{array}$

The maximum tongue strain is found near the base of the papilla, r ~ a where a is the radius of the papilla leading to the relevant strain value of

$\begin{array}{ll}\epsilon ~\frac{M}{E{a}^3}·& (27)\end{array}$

5.5. Strains in Tongue from Papilla Bending

The result shown by Equation (27) needs now to be compared with what would happen in the absence of papillae-induced deformation.

5.5.1. Shear Flow

In the case of a shear flow, in the absence of papillae the pressure in the Newtonian fluid is constant and equal to its atmospheric value, meaning that the tongue would undergo no forcing normal to its surface. The shear flow would of course lead to shear stresses, but no shear strain in the case of an approximately incompressible tongue. So as far as strain receptors are concerned, they would not measure anything if there were no papillae in this situation.

What is the order of magnitude of the tongue strains, ε_∥, we would obtain as induced by the presence of papillae? The moment at the base of an individual papilla scales as M ~ aF with the force F ~ a²σ corresponding to typical bending stress σ ~ E_pϵ, where ϵ denotes the longitudinal strain in the papilla and E_p the Young's modulus of the papilla. We therefore finally obtain an applied moment $M~a^3{E}_pϵ$.

We saw in Equation (18) that in the case of a shear flow we have a longitudinal strain inside the papilla

$\begin{array}{ll}{ϵ}_{\parallel }~\frac{a\xi U{L}^3}{BH},& (28)\end{array}$

and thus the moment applied at the base of the papilla on the tongue due to shear-induced bending is given by

$\begin{array}{ll}M~\frac{\xi a^4{E}_{p}U{L}^3}{BH}·& (29)\end{array}$

Using the scaling relationship between the Young's modulus of the papilla and the bending rigidity $B~a^4{E}_p$, we obtain a simpler formula for the torque magnitude as

$\begin{array}{ll}M~\frac{\xi U{L}^3}{H}·& (30)\end{array}$

With this value for the moment, we can then use Equation (27) to find the maximum elastic strains in the tongue

$\begin{array}{ll}{\epsilon }_{\parallel }~\frac{\xi U{L}^3}{a^{3}EH}·& (31)\end{array}$

Provided that the Young's modulus of the tongue is similar to that of individual papillae, which as a first approximation it should, then writing Equation (31), we recognize the elasto-hydrodynamics time scale, t_eh, from Equation (8), as a function of

$\begin{array}{ll}{\epsilon }_{\parallel }~t_{eh}\frac{aU}{LH}·& (32)\end{array}$

Putting in numbers, we recall that t_eh = 10 ms, U = 10 cm/s, L = 250 μm a = 50 μm, E = 25 kPa, and H = 5 mm. This leads to t_eh ≈ 3 ms and gives a maximum tongue strain of

$\begin{array}{ll}{\epsilon }_{\parallel }~4\%.& (33)\end{array}$

Thus in shear flow we propose the advantage of having papillae distributed on the tongue: they create a normal mechanical deformation of the tongue which would not take place in their absence.

5.5.2. Squeeze Flow

In the case of a shear flow, the scalings are different. As in previous cases discussed, we expect to get an even stronger effect than in the longitudinal situation.

From Equation (22), the magnitude of the elongational strain in the papilla due to its bending is given by

$\begin{array}{ll}{ϵ}_{\perp }~\frac{a\xi VW{L}^3}{B{H}^2},& (34)\end{array}$

and therefore the torque applied at the base of the papilla on the tongue has a typical magnitude

$\begin{array}{ll}M~\frac{\xi a^4{E}_{p}VW{L}^3}{B{H}^2}~\frac{\xi VW{L}^3}{H^2}·& (35)\end{array}$

With this, the typical strains in the tongue in the case of a squeeze flow, ε_⊥, are

$\begin{array}{ll}{\epsilon }_{\perp }~\frac{\xi VW{L}^3}{a^{3}E{H}^2}·& (36)\end{array}$

Comparing Equations (31) and (36) we find the ratio

$\begin{array}{ll}\frac{{\epsilon }_{\perp }}{{\epsilon }_{\parallel }}~\frac{VW}{UH}·& (37)\end{array}$

Using the same numbers as above except W = 5 cm, V = 10 cm/s, this leads to a significantly larger

$\begin{array}{ll}{\epsilon }_{\perp }~40\%.& (38)\end{array}$

This result needs to be compared with what would happen in the absence of papillae. In that case the elastic stress in the tongue would be uniform everywhere and equal to the fluid pressure. Approximately, this hydrodynamic pressure is given by

$\begin{array}{ll}p~\frac{\mu V{W}^2}{H^3},& (39)\end{array}$

and thus the strains in the tongue in the absence of papillae, ε_p, would be

$\begin{array}{ll}{\epsilon }_p~\frac{p}{E}~\frac{\mu V{W}^2}{E{H}^3},& (40)\end{array}$

leading to

$\begin{array}{ll}\frac{{\epsilon }_p}{{\epsilon }_{\perp }}~\frac{\mu }{\xi }\frac{W}{H}{\left(\frac{a}{L}\right)}^3·& (41)\end{array}$

With W = 5 cm, L = 250 μm, a = 50 μm, H = 5 mm this leads to

$\begin{array}{ll}\frac{{\epsilon }_p}{{\epsilon }_{\perp }}~0.02,& (42)\end{array}$

so the presence of papillae leads to an increase of two orders of magnitude of the strains in the tongue.

6. Conclusion

The analytical and numerical solutions derived in this article confirm that the topology of the tongue can offer mechanical advantage for texture perception. This effect is dependent on the presence and location of mechanosensitive cells and neurons in relation to such topological structures. This, to the best of our knowledge, remains to be studied, however recent work reporting corpuscular endings in mice filiform papillae indicates locations consistent with the mechanism described here [9].

Most interestingly, we prove that depending on the putative location of those mechanosensitive cells, filiform papillae may either serve sensing purposes in a direct fashion (papillae strains), similarly to superficial neuromasts or in an indirect fashion (tongue strains). In the latter, instead of sensing strains themselves, papillae might encode stresses by inducing strains in the tongue which then lead to bulk sensing. The underlying neurophysiology for those two cases would be dramatically different in terms of strain orientation and magnitude. The application from the developed model could then be used to define the levels of stress and strain that are experienced by different cells depending on the physiological location.

One can also speculate that similarly to the skin, the tongue is populated by various types of cellular mechanosensors, encoding information from different locations and responding to different levels of stresses and frequencies. In this case both direct and indirect sensing may occur from the filiform papillae bending.

As a final perspective, since we have argued that the aspect ratio of filiform papillae is similar to that of other biological structures known to support external stress field encoding, we propose that models similar to ours will also be applicable to such structures. We hope that biologists and biophysicists will be able to apply these ideas to other biological systems.

Author Contributions

BL and CP designed the project; EL derived the analytical model; EL, CP, and BL wrote the article.

Conflict of Interest Statement

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

BL and CP are employees of the Nestlé Research Center. EL was consulting with Nestlé Research Center during the study.

References

1. Burbidge A, Le Révérend B. Biophysics of food perception. J Phys D Appl Phys. (2016) 49:114001. doi: 10.1088/0022-3727/49/11/114001

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Chandrashekar J, Hoon MA, Ryba NJP, Zuker CS. The receptors and cells for mammalian taste. Nature (2006) 444:288–94. doi: 10.1038/nature05401

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Ishida Y, Ugawa S, Ueda T, Murakami S, Shimada S. Vanilloid receptor subtype-1 (VR1) is specifically localized to taste papillae. Mol Brain Res. (2002) 107:17–22. doi: 10.1016/S0169-328X(02)00441-2

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Dhaka A, Earley TJ, Watson J, Patapoutian A. Visualizing cold spots: TRPM8-expressing sensory neurons and their projections. J Neurosci. (2008) 28:566–75. doi: 10.1523/JNEUROSCI.3976-07.2008

PubMed Abstract | CrossRef Full Text

5. Scott Blair GW. Measurements of Mind and Matter (1950).

6. de Wijk RA, Prinz JF, Engelen L, Weenen H. The role of α-amylase in the perception of oral texture and flavour in custards. Phys Behav. (2004) 83:81–91. doi: 10.1016/j.physbeh.2004.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Peyron MA, Mishellany A, Woda A. Particle size distribution of food boluses after mastication of six natural foods. J Dental Res. (2004) 83:578–82. doi: 10.1177/154405910408300713

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Engelen L, van der Bilt A, Schipper M, Bosman F. Oral size perception of particles: effect of size, type, viscosity and method. J Text Stud. (2005) 36:373–86. doi: 10.1111/j.1745-4603.2005.00022.x

CrossRef Full Text | Google Scholar

9. Moayedi Y, Duenas-Bianchi L, Michlig S, le Coutre J, Le Révérend B, Lumpkin E. Somatosensory contribution to food oral processing. In: Proceedings of the 4th Food Oral Processing Conference. Lausanne (2016).

10. Sato O, Maeda T, Kobayashi S, Iwanaga T, Fujita T. Filiform papillae as a sensory apparatus in the tongue: an immunohistochemical study of nervous elements by use of neurofilament protein (NFP) and S-100 protein antibodies. Cell Tissue Res. (1988) 252:231–8. doi: 10.1007/BF00214365

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Malone AMD, Anderson CT, Tummala P, Kwon RY, Johnston TR, Stearns T, et al. Primary cilia mediate mechanosensing in bone cells by a calcium-independent mechanism. Proc Natl Acad Sci USA. (2007) 104:13325–30. doi: 10.1073/pnas.0700636104

PubMed Abstract | CrossRef Full Text | Google Scholar

12. McHenry MJ, van Netten SM. The flexural stiffness of superficial neuromasts in the zebrafish (Danio rerio) lateral line. J Exp Biol. (2007) 210:4244–53. doi: 10.1242/jeb.009290

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Triantafyllou MS, Weymouth GD, Miao J. Biomimetic survival hydrodynamics and flow sensing. Annu Rev Fluid Mech. (2016) 48:1–24. doi: 10.1146/annurev-fluid-122414-034329

CrossRef Full Text | Google Scholar

14. Baratz RS, Farbman AI. Morphogenesis of rat lingual filiform papillae. Am J Anat. (1975) 143:283–301. doi: 10.1002/aja.1001430303

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Nagato T, Nagaki M, Murakami M, Tanioka H. Morphological studies of rat lingual filiform papillae. Okajimas Folia Anat Jpn. (1989) 66:195–209. doi: 10.2535/ofaj1936.66.4_195

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Toprak B. Light and scanning microscopic structure of filiform papillae in mice. Veterinarski Arhiv. (2007) 76:555. Available online at: http://hrcak.srce.hr/22261

17. Karan M, Yilmaz S, Aydin A. Morphology of the filiform lingual papillae in porcupine (Hystrix cristata). Anat Histol Embryol. (2011) 40:100–3. doi: 10.1111/j.1439-0264.2010.01045.x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Chamorro CA, Fernandez JG, Paz P, Pelaez B, Anel L. Scanning electron microscopy of the wild boar and pig lingual papillae. Histol Histopathol. (1994) 9:657–67.

Google Scholar

19. Toyoda M, Sakita S, Kagoura M, Morohashi M. Electron microscopic characterization of filiform papillae in the normal human tongue. Arch Hist Cytol. (1998) 61:253–68. doi: 10.1679/aohc.61.253

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Manabe M, Lim HW, Winzer M, Loomis CA. Architectural organization of filiform papillae in normal and black hairy tongue epithelium: dissection of differentiation pathways in a complex human epithelium according to their patterns of keratin expression. Arch Dermat. (1999) 135:177–81. doi: 10.1001/archderm.135.2.177

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Iwasaki S. Evolution of the structure and function of the vertebrate tongue. J Anat. (2002) 201:1–13. doi: 10.1046/j.1469-7580.2002.00073.x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yamashita K, OdDalkhsuren S. The papillae of the human tongue. In: Kat H, Shimizu T, editors. Tongue: Anatomy, Kinematics and Diseases. Hauppauge, NY: Nova Science Publishers, Inc. (2012). pp. 143–54.

23. Ranc H, Servais C, Chauvy PF, Debaud S, Mischler S. Effect of surface structure on frictional behaviour of a tongue/palate tribological system. Trib Int. (2006) 39:1518–26. doi: 10.1016/j.triboint.2006.01.017

CrossRef Full Text | Google Scholar

24. Dresselhuis DM, De Hoog EHA, Cohen Stuart MA, Van Aken GA. Application of oral tissue in tribological measurements in an emulsion perception context. Food Hydrocoll. (2008) 22:323–35. doi: 10.1016/j.foodhyd.2006.12.008

CrossRef Full Text | Google Scholar

25. Kamm R, Lammerding J, Mofrad M. Cellular nanomechanics. In: Bhushan B, editor. Springer Handbook of Nanotechnology. New York, NY: Springer (2010). pp. 1171–200.

26. Berdyyeva TK, Woodworth CD, Sokolov I. Human epithelial cells increase their rigidity with ageing in vitro: direct measurements. Phys Med Biol. (2005) 50:81. doi: 10.1088/0031-9155/50/1/007

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Spassova I. Ultrastructure of the simple encapsulated nerve endings (simple end-bulbs of Krause) in the tongue of the cat. J Anat. (1974) 118(Pt 1):1–9.

PubMed Abstract | Google Scholar

28. Lumpkin EA, Caterina MJ. Mechanisms of sensory transduction in the skin. Nature (2007) 445:858–65. doi: 10.1038/nature05662

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Shawker TH, Sonies B, Stone M, Baum BJ. Real-time ultrasound visualization of tongue movement during swallowing. J Clin Ultrasound (1983) 11:485–90. doi: 10.1002/jcu.1870110906

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kristensen D, Jensen PY, Madsen F, Birdi KS. Rheology and surface tension of selected processed dairy fluids: influence of temperature. J Dairy Sci. (1997) 80:2282–90. doi: 10.3168/jds.S0022-0302(97)76177-0

CrossRef Full Text | Google Scholar

31. Leighton A, Williams OE. The effect of temperature on the basic viscosity of ice-cream mixes. J Phys Chem. (1927) 31:1663–8. doi: 10.1021/j150281a004

CrossRef Full Text | Google Scholar

32. Adapa S, Dingeldein H, Schmidt KA, Herald TJ. Rheological properties of ice cream mixes and frozen ice creams containing fat and fat replacers. J Dairy Sci. (2000) 83:2224–9. doi: 10.3168/jds.S0022-0302(00)75106-X

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Muse MR, Hartel RW. Ice cream structural elements that affect melting rate and hardness. J Dairy Sci. (2004) 87:1–10. doi: 10.3168/jds.S0022-0302(04)73135-5

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Nicosia MA, Robbins J. The fluid mechanics of bolus ejection from the oral cavity. J Biomech. (2001) 34:1537–44. doi: 10.1016/S0021-9290(01)00147-6

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Kim S, Karilla JS. Microhydrodynamics: Principles and Selected Applications. Boston, MA: Butterworth-Heinemann (1991).

Google Scholar

36. Happel J, Brenner H. Low Reynolds Number Hydrodynamics. Englewood Cliffs, NJ: Prentice Hall (1965).

Google Scholar

37. Batchelor GK. An Introduction to Fluid Dynamics. Cambridge, UK: Cambridge University Press (1967).

Google Scholar

38. Powers TR. Dynamics of filaments and membranes in a viscous fluid. Rev Mod Phys. (2010) 82:1607–31. doi: 10.1103/RevModPhys.82.1607

CrossRef Full Text | Google Scholar

39. Wiggins CH, Goldstein RE. Flexive and propulsive dynamics of elastica at low Reynolds number. Phys Rev Lett. (1998) 80:3879–82. doi: 10.1103/PhysRevLett.80.3879

CrossRef Full Text | Google Scholar

40. Landau LD, Lifshitz EM. Theory of Elasticity. 3rd Edn. Oxford: Pergamon Press (1986).

41. Woo SH, Ranade S, Weyer AD, Dubin AE, Baba Y, Qiu Z, et al. Piezo2 is required for Merkel-cell mechanotransduction. Nature (2014) 509:622–6. doi: 10.1038/nature13251

PubMed Abstract | CrossRef Full Text | Google Scholar