- 1Department of Mechanical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, United States
- 2Theoretical and Applied Mechanics Program, McCormick School of Engineering, Northwestern University, Evanston, IL, United States
- 3Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States
Introduction: Plotting the pressure-cross-sectional area (P-CSA) hysteresis loops within the esophagus during a contraction cycle can provide mechanistic insights into esophageal motor function. Pressure and cross-sectional area during secondary peristalsis can be obtained from the functional lumen imaging probe (FLIP). The pressure-cross-sectional area plots at a location within the esophageal body (but away from the sphincter) reveal a horizontal loop shape. The horizontal loop shape has phases that appear similar to those in cardiovascular analyses, whichinclude isometric and isotonic contractions followed by isometric and isotonic relaxations. The aim of this study is to explain the various phases of the pressurecross-sectional area hysteresis loops within the esophageal body.
Materials and Methods: We simulate flow inside a FLIP device placed inside the esophagus lumen. We focus on three scenarios: long functional lumen imaging probe bag placed insidethe esophagus but not passing through the lower esophageal sphincter, long functional lumen imaging probe bag that crosses the lower esophageal sphincter, and a short functional lumen imaging probe bag placed in the esophagus body that does not pass through the lower esophageal sphincter.
Results and Discussion: Horizontal P-CSA area loop pattern is robust and is reproduced in all three cases with only small differences. The results indicate that the horizontal loop pattern is primarily a product of mechanical conditions rather than any inherently different function of the muscle itself. Thus, the distinct phases of the loop can be explained solely based on mechanics.
1 Introduction
A healthy functioning esophagus transports swallowed material from the mouth towards the stomach through peristaltic contraction (Regan et al., 2012; Jain et al., 2019). Patients with systemic sclerosis, achalasia, dysphagia, and other motor disorders have impaired esophageal peristalsis (Kahrilas, 2008; Shaheen et al., 2010; Gregersen et al., 2011; Aziz et al., 2016). Thus, understanding esophageal motor function and investigating the contractile behavior of esophageal smooth muscles can help in identifying esophageal disease progression and develop diagnostic tools. Extensive computational studies have been performed in recent years on the mechanics of esophageal transport (Kou et al., 2017a; Kou et al., 2017b; Kou et al., 2018; Kou et al., 2017c; Acharya et al., 2022; Halder et al., 2021; Elisha et al., 2022a; Acharya et al., 2021c; Acharya et al., 2021a; Elisha et al., 2021; Liao et al., 2020) and emptying as well as the state of various esophageal diseases (Acharya et al., 2021b; Halder et al., 2022a; Halder et al., 2022c; Halder et al., 2022b).
In recent years, several studies used pressure-cross-sectional area (P-CSA) relationship during esophageal contraction cycle to investigate esophageal wall function (Villadsen et al., 2001; Gregersen et al., 2011; Liao et al., 2013; Leibbrandt et al., 2016; Elisha et al., 2021). This approach is borrowed from cardiovascular analyses, which uses pressure-volume (P-V) plots to study dynamic motor activity and explain the pumping mechanism of the left ventricle and its dysfunction (Gregersen et al., 2007; Gregersen et al., 2011; Liao et al., 2013; Gregersen and Lo, 2018).
Tracking the pressure and CSA during a distention-induced contraction cycle of the esophagus was made possible by the functional lumen imaging probe (FLIP) (Liao et al., 2014; Moonen and Boeckxstaens, 2014; Lottrup et al., 2015). The FLIP device, illustrated in Figure 1, consists of a catheter surrounded by a fluid filled bag, placed inside the esophagus. The filled bag results in mechanical distention, which induces secondary peristalsis in the esophagus (Clarke et al., 2020). The response of the esophagus wall to the distention is captured by the FLIP in the form of CSA readings at 16 locations and one distal pressure measurement as a function of time (Carlson et al., 2015; Carlson et al., 2016). The FLIP in the esophagus can be placed such that it passes through the lower esophageal sphincter (LES), as in the middle illustration in Figure 1, or not, as in the left illustration in Figure 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1066351/fphys-13-1066351-HTML-r1/image_m/fphys-13-1066351-g001.jpg)
FIGURE 1. Diagram showing the FLIP bag and catheter assembly placed within the esophagus, passing through the LES or not. The bag used in this study is 16 cm long and 22 mm diameter. The catheter measures CSA at 16 location and one distal pressure measurement at each time instant.
As reported by Gregersen et al. (2011), the P-CSA hysteresis loop, at any location within the esophagus body (i.e., away from the lower esophageal sphincter or the LES), during distention-induced secondary peristalsis reveals a horizontal loop shape, presented in Figure 2A. This loop pattern looks similar to the one obtained from the left ventricle, and therefore, the loop was explained analogously. The esophageal P-CSA hysteresis loop was divided into four phases, each named analogously after the corresponding phase on the P-V loop plotted for the left ventricle. These phases are isometric contraction, isotonic contraction, isometric relaxation, and isotonic relaxation (Gregersen et al., 2011). Similar terminology was used by others Leibbrandt et al. (2016); Omari et al. (2015) when analyzing P-CSA hysteresis loop within the esophagus body or at the upper esophageal sphincter using a different measuring device.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1066351/fphys-13-1066351-HTML-r1/image_m/fphys-13-1066351-g002.jpg)
FIGURE 2. P-CSA relationships in the esophagus body for four different contraction cycles with different experimental setups. (A) P-CSA extracted from a 4-cm long FLIP obtained by Gregersen et al. (2011). (B) and (C) P-CSA extracted from a 16-cm long FLIP at 60 ml, where the FLIP bag passes through the LES. These hystereses are plotted from data recorded at a sensor located away from the LES. (D) P-CSA extracted from a 16-cm long FLIP at 60 ml, where the FLIP bag does not cross the LES.
While the resulting P-CSA loop within the esophagus body during the FLIP experiment and the P-V loop in the left ventricle have a similar pattern, the mechanical functioning is different. For example, the volume inside the ventricle changes through the cardiac cycle whereas the volume inside the esophagus during secondary peristalsis in a FLIP device is constant. Here, we aim to provide a mechanics-based explanation for the P-CSA loop within the esophagus body. Specifically, we show that the various phases in the loop shape are a result of the mechanical conditions (muscle contraction, boundary condition, and geometric configuration) rather than any inherently different function of the muscle itself.
1.1 Pressure-CSA loops within the esophagus body
Figure 2 presents four plots of the pressure within the esophageal body as a function of wall CSA during contraction cycles from four different experiments. Figure 2A displays the results obtained by Gregersen et al. (2011). In their work, the FLIP bag was 4 cm long and the maximum bag diameter was 5 cm. The middle of the bag was placed 7 cm above the LES and did not pass through it. Figures 2B, C display two examples of P-CSA loops obtained by a 16 cm FLIP bag at 60 ml bag fill, where the bag passes through the LES. Lastly, the loop plotted in Figure 2D is extracted by a 16 cm FLIP bag at 60 ml bag fill, where the bag does not cross the LES.
The data for Figures 2B–D was collected at the Esophageal Center of Northwestern between November 2012 and October 2018, using EndoFLIP® (EF-322N; Medtronic, Inc., Shoreview, MN, United States) (Acharya et al., 2021a; Carlson et al., 2021). We examined a total of 24 randomly selected, asymptomatic volunteers for the cases where the FLIP passes through the LES, and two cases where it does not. Informed consent was obtained for subject participation; they were paid for their participation. The study protocol was approved by the Northwestern University Institutional Review Board. Additional details on the data collection process and subject cohort selection is available in (Lin et al., 2013; Carlson et al., 2016; Carlson et al., 2021; Acharya et al., 2021a). Since FLIP provides only one distal pressure reading, pressure at different locations along the esophagus length was calculated using the distal pressure, CSA measurements, the mass conservation equation, and the momentum equation (discussed in Section 2.1), as proposed by Halder et al. (2021).
As Figure 2 reveals, the P-CSA hysteresis plotted at a location within the esophagus body displays a loop of average horizontal orientation. This is opposed to the loops of positive or negative slope orientation which have been reported at sphincter locations of different organs (Zifan et al., 2019; Elisha et al., 2022b). A loop in considered horizontal when the magnitude of the slope of its best fit line for the point cloud representing the hysteresis is such that |slope| < .1. The point cloud is obtained by discretizing the loop into equal CSA segments (∆A) and finding the corresponding pressure values. In the positive and negative slope loops at sphincter locations of different organs, the slopes are clearly greater in magnitude than the .1 threshold used above to define a horizontal loop.
On each plot in Figure 2, the four phases identified by Gregersen et al. (2011) are marked. Gregersen et al. (2011) refer to phases 1 and 2 as contraction, while phases 3 and 4 as relaxation. To explain these phases, we reproduce esophageal secondary peristalsis using simulations, discussed next.
2 Materials and methods
In our previous studies, we developed a one dimensional (1D) model of a flow inside an elastic tube closed on both ends to imitate a flow inside a FLIP device (Elisha et al., 2022a; Elisha et al., 2021; Acharya et al., 2021c; Halder et al., 2022c). Our goal was to study the relation between tube properties, fluid properties, muscle activation pattern, and their effect on pressure and CSA of the tube. In this section, we summarize the mathematical model.
2.1 Governing equations in 1D
The 1D mass and momentum conservation equations are
and
respectively, which were derived by Ottesen (2003). Here, A(x, t), u(x, t), P(x, t), ρ and μ are the tube CSA, fluid velocity (averaged at each cross-sectional area), pressure inside the tube, fluid density, and fluid viscosity, respectively. To complete the system, we relate pressure and CSA by a linear constitutive equation of the form
This relation, also known as the “Tube Law,” was derived by Whittaker et al. (2010) and validated experimentally by Kwiatek et al. (2011). In the equation above, Po is the outside pressure, Ke is tube stiffness, and Ao is the undeformed reference area (CSA of the tube when P = Po) (Elisha et al., 2022a; Elisha et al., 2021; Elisha et al., 2022b; Acharya et al., 2021c). Lastly, to mimic muscle peristaltic contraction, we vary the reference area sinusoidally by multiplying Ao by an activation function θ(x, t) (Ottesen, 2003; Manopoulos et al., 2006; Bringley et al., 2008; Acharya et al., 2021c; Elisha et al., 2022a). As Eq. 3 suggests, contraction implies reduction of the tube’s reference CSA (θ < 1), and relaxation is associated with increase in the reference CSA (θ > 1). When θ = 1, neither contraction nor relaxation are present. Note that in the simulations described in this report, θ is either less than or equal to 1. The expression for θ(x, t) is discussed in Section 2.3.
2.2 Non-dimensional dynamic equations
We non-dimensionalize the dynamic system using
to reduce the number of independent variables. The terms α, τ, U, p, and χ are the non-dimensional variables of area, time, velocity, pressure, and position, respectively (Acharya et al., 2021c; Elisha et al., 2022a). The dimensional wave speed and tube length are
Using the above, the mass conservation, momentum conservation, and tube law equations can be written as
respectively. The terms
2.3 Peristaltic wave input and active relaxation
The activation function (θ) can take different forms, and implementing it as described in Section 2.1 is a common approach when the external activation pressure at a specific location varies sinusoidally with time. As discussed in Elisha et al. (2022a), the θ pattern chosen for the model in this study aims to mimic clinical observations of esophageal secondary peristalsis reported by Goyal and Chaudhury (2008) and Crist et al. (1984). When the FLIP bag does not pass through the sphincter, we define the activation function as simply a traveling sinusoidal wave with constant amplitude, wavelength and wave speed of the form
where w is the non-dimensional width of the peristaltic wave (dimensional width/L) and θc is the maximum contraction strength (wave amplitude = θc) (Acharya et al., 2021c; Elisha et al., 2022a). Figure 3A presents θp(χ, τ) at five consecutive time instants. Given that an average esophageal peristaltic wave is about 3–4 cm wide, w = .25 for a 16 cm bag and w = .75 for a 4 cm bag. In the activation function presented in Figure 3A, w = .25.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1066351/fphys-13-1066351-HTML-r1/image_m/fphys-13-1066351-g003.jpg)
FIGURE 3. Activation functions (A) θp and (B) θsup (B) plotted as a function of χ at five consecutive time instants.
For the scenario where the FLIP bag passes through the LES, the activation function is defined as a superposition of θp(χ, τ) with θLES(χ, τ), where θLES(χ, τ) is responsible for the initial tone and relaxation of the LES, such that
The function θLES(χ, τ) is defined by
where θr is the LES tone, wLES is the non-dimensional width of the LES, χ1 and χ2 are the boundaries of the LES segment (χ2 > χ1 and wLES = χ2−χ1), and
In the equation above, τi and τm are time instants, where τi marks the beginning of LES relaxation and τm marks the end of the relaxation process (Elisha et al., 2021).
2.4 Initial and boundary conditions
At τ = 0, the fluid inside the tube is at rest, thus, its velocity (U) is equal to zero for all χ (Acharya et al., 2021c; Elisha et al., 2022a). Additionally, it implies that the pressure in the tube is uniform (typically positive), which according to the tube law indicates that α/θ (at τ = 0) is a constant (= SIC). Consequently, the initial condition can be written as
The velocity boundary condition is simply
since the tube is closed at both ends. The boundary condition for α is not as straightforward, but can be derived by plugging Eq. 7 and the velocity condition in Eq. 13 into Eq. 6, which yields
This condition is constructed by assuming that η ≈ 0. However, as explained by Acharya et al. (2021c), the effects due to this damping at the boundary are negligible for small η.
2.5 Numerical implementation
The system is solved using the MATLAB
In previous studies, a large range of possible values for the different parameters were examined to conduct a parametric study. For example, although typical clinical values correspond to small β (order 1) and large ψ (order 103), analyzing simulations only within this range means only capturing scenarios in which viscous effects are negligible. Therefore, larger β and smaller ψ were also considered (Elisha et al., 2022a; Elisha et al., 2021; Acharya et al., 2021c). However, based of the findings summarized in Elisha et al. (2022a) and Elisha et al. (2021), we choose to focus on simulations with parameter values listed in Table 1.
3 Results and discussion
Clinical data reveal that the P-CSA hysteresis loop during a single contraction cycle at a location within the esophageal body recorded by a FLIP device results a horizontal loop shape (Figure 2). To explain the underlying mechanisms of this contraction pattern by analyzing the different phases of these loops, we simulated different contraction cycles, as described in the previous section. In this section, we display and discuss the simulation results by plotting the P-CSA hysteresis for different scenarios. We focus on the horizontal loop orientation obtained at the esophageal body using a 16 cm FLIP.
Figure 4 presents the results of a simulation of an esophageal contraction cycle using a 16 cm long FLIP which does not pass through the LES. The activation function is θp with the input parameters listed in Table 1, w = .25, and ψ = 2400. The resulting tube shape at seven consecutive time instants are presented on the right and the location where the esophagus body loop is plotted is marked χ = χo. The P-CSA loop at that location is plotted on the top left, and the pressure and CSA as a function of time are plotted on the bottom left. As the figure shows, a horizontal loop emerges.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1066351/fphys-13-1066351-HTML-r1/image_m/fphys-13-1066351-g004.jpg)
FIGURE 4. Simulation results for a single contraction cycle with muscle activity defined by a traveling peristaltic wave (θp). Top left: P-CSA hysteresis loop at χo. Bottom left: P and CSA vs. time at χ = χo. Right: tube shape at seven consecutive instants, ordered chronologically.
As in the prior clinical analysis by Gregersen et al. (2011), we divide the loop into phases (a–e) (Figure 4), which will help us understand the underlying mechanics causing the horizontal loop shape. To explain the different phases, we use the diagram in Figure 5. The diagram presents the shape of the elastic tube at two time instants during a contraction cycle. In the simulation, the wave travels from left to right (proximal to distal). The P-CSA loop is plotted (Figure 4) at the point χ = χo, which is static. Point χ = c is where the maximum contraction of the advancing wave takes place (θ = θc), and it moves with time. Regions 1 and 2 refer to the regions distal and proximal to the contraction, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1066351/fphys-13-1066351-HTML-r1/image_m/fphys-13-1066351-g005.jpg)
FIGURE 5. Two instants showing the corresponding tube shapes from the same simulation. The locations of point χ = c, regions 1 and 2, and how they change over time are depicted. P-CSA is plotted at χ = χo and its location does not change with time.
As the peristalsis wave travels down the deformable tube, it pushes the fluid forward, such that the CSA, A1, at region 1 is greater than the CSA, A2, at region 2 (Acharya et al., 2021c; Elisha et al., 2022a). Therefore, the pressure, p1, at region 1 is also greater than the pressure, p2, at region 2. The minimum CSA, Ac, is at χ = c. Thus, when χo is distal to the contraction,
The tube shape and the relation between A1, A2, and Ac were described by in our prior work (Acharya et al., 2021c; Elisha et al., 2022a). As the advancing contraction wave travels down the deformable tube, a competition between the elastic forces generated due to deformation of the tube wall and resistance to flow through the narrowest part of the contraction is created. When tube stiffness is high and fluid viscosity is relatively low, the tube wall resists deformation and forces the fluid that was displaced by the peristaltic wave to flow back through the contraction. In this scenario, the energy that is required to expand the tube walls is greater than the energy required to overcome viscous resistance across the contraction. Yet, A1 > A2, even if the difference between the two is small. When tube stiffness is relatively low and fluid viscosity is relatively high, an opposite phenomenon occurs. In this scenario, the resistance to flow is high, and therefore it is favorable for the tube wall at region 1 to expand and accommodate for the displaced fluid. In this second scenario, A1 is much greater than A2. With the above understanding, the different phases, as marked in Figure 4 are discussed next.
(a) As the contraction cycle starts, the peristaltic wave pushes the fluid forward, increasing the pressure distal to the contraction, which is also where point χo is located (Elisha et al., 2022a). Because fluid is displaced, CSA also increases at χ = χo, alongside pressure. Due to the high tube stiffness (ψ = 2400), the increase in pressure and CSA is small. When viscosity increases and/or tube stiffness decreases, it is favorable for the tube to expand and accommodate for the displaced fluid. Thus, the line in phase (a) is longer as the change in pressure and CSA is more prominent, as seen in Figure 6A. Figure 6B displays a visual representation of the affect wall stiffness and fluid viscosity have on the loop, and specifically stage (a). Phase (a) in esophageal FLIP experiment is not perfectly vertical or isometric (Figures 2, 4).
(b) When the contraction wave approaches χo, the CSA at this location starts to decrease due to the contraction (i.e., θ(χo, τ) decreases). However, during phase (b) point χo is still distal to the minimum area at c, and θ(χo, τ) > θc. Thus, χo is still part of region 1 (Figure 5) where the pressure is approximately uniform at high wall stiffness and low fluid viscosity. This is because, in this case, the dynamic pressure gradients (pressure changes due to fluid motion) are not dominant in the momentum equation. Consequently, although CSA decreases, pressure remains approximately constant
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1066351/fphys-13-1066351-HTML-r1/image_m/fphys-13-1066351-g006.jpg)
FIGURE 6. P-CSA plots obtained from three different contraction cycle simulations. Input parameters are as listed in Table 1 and ψ = 100. (A, B) The FLIP bag does not pass through the LES, w = .25, (C) the FLIP bag passes through the LES, w = .25, and (D) simulation results for a short FLIP bag, w = .75.
Dynamic pressure gradients do become relevant for low tube stiffness or high viscosity (to overcome fluid friction). In this case, pressure gradients develop within region 1 (Figure 5) and one can see that, although minor compared to phase (a),
(c) In this phase, point χo transitions from being proximal to being distal to the contraction (point c). At the middle of phase (c), θ(χo, τ) = θc and CSA at χo is minimum. Around the minimum area the viscous resistance for the fluid to flow from region 1 to region 2 is maximum. The pressure gradient across the minimum area provides the force necessary to overcome the viscous resistance. Thus,
(d) As the contraction wave exits location χo, θ(χo, τ) increases, which implies that CSA increases. However, χ = χo is proximal to the contraction and in region 2 (Figure 5). For reasons similar to those discussed in phase (b), the pressure is almost uniform in region 2. Consequently,
(e) The constant pressure during phase (d) is lower than the initial pressure. In phase (e), the pressure at χ = χo increases back to the initial pressure. As the contraction wave reaches the closed distal end, and eventually exits the domain, it allows more back-flow, causing the displaced fluid to flow upstream. This enhanced back-flow increases the pressure at point χo, as the tube returns to its initial state. Note that in clinical data, the contraction is a repetitive process, meaning that as the advancing contraction wave exits at the distal end of the tube, another peristaltic wave enters the domain at the proximal end (Carlson et al., 2020). Given this set up, it is hard to separate phases (a) and (e), and the loop looks like one continuous process.
An additional complication is introduced when the FLIP bag is placed such that it crosses the LES. As the peristaltic contraction cycle starts, (phase (a)), the LES muscle begins to relax, as illustrated by the activation function in Figure 3B. Thus, the pressure in the tube is controlled by LES tone and the peristaltic contraction, each mechanism affecting the pressure at χ = χo differently. Peristaltic wave increases pressure distal of the contraction. Sphincter relaxation leads to pressure drop. Figure 6C presents the resulting P-CSA loop of such a simulation. As it turns out, although both mechanisms are present and cause different outcomes in pressure, the resulting loop has the same horizontal pattern. Most importantly, the pressure at point χ = χo during phase (a) increases. The reason for that was explained in a prior work (Elisha et al., 2021), where the P-CSA hysteresis loop at the LES was plotted. Their study revealed that pressure and CSA at the LES increase and decrease together. That is, as the LES opens, the pressure at the LES increases. While originally counter intuitive, this P-CSA relation was explained by Elisha et al. (2022b). The LES opening is a superposition of two mechanisms; neurogenic mediated relaxation of the LES muscle, and peristaltic contraction that causes pressurized distention. Thus, although the LES opens, which in an isolated case would reduce the pressure in region 1, the pressure increase from the peristaltic contraction dominates.
With this understanding in mind, the parallel between the horizontal loops in Figures 6A, C is clear. In phase (a), as the peristaltic wave travels down the deformable tube, the pressure and therefore the CSA increases distal to the contraction. Although the LES opens, the pressure rise from the traveling contraction dominates such that the pressure increases in phase (a). Note that this pressure increase is not as high as in the case where the FLIP does not pass through the LES (in Figure 6A).
By the time the contraction wave reaches point χo, marking the beginning of phase (b), the LES is fully open and remains open for the rest of the cycle. Hence, phases (b), (c), (d) are analogous to the ones described for Figure 4.
Lastly, note that the P-CSA loop plotted by Gregersen et al. (2011) and presented in Figure 2A has a slightly different shape than the remaining clinical cases presented in Figure 2. This loop does not have a well defined phase (a). Since the FLIP length is the main parameter that differentiates the case in Figure 2A from the rest, we hypothesize that the FLIP length causes this difference. To test this hypothesis, we simulate flow inside a short FLIP placed in the esophagus by changing the contraction width (w) from .25 (as in the case in Figure 6A) to w = .75. The P-CSA loop from this simulation is displayed in Figure 6D. It is seen that decreasing the tube length conserves the general horizontal orientation of the loop but eliminates the original phase (a). The resulting loop has a more drop-like shape, similar to the corresponding clinical loop by Gregersen et al. (2011). The reason for that is simple: the wave hardly travels outside point χo because the width of the contraction is almost equal to the width of the FLIP bag.
4 Conclusion
The goal of this study was to provide a mechanistic understanding of the P-CSA loops within the esophageal body. We examined the P-CSA hysteresis loop within the esophageal body obtained from clinical FLIP testing of various configurations. The considered cases were 1) 16-cm long FLIP bag placed inside the esophagus but did not pass through the LES, 2) 16-cm long FLIP bag which crossed the LES, and 3) an experimental setup reported by Gregersen et al. (2011), in which a 4-cm long FLIP bag was placed in the esophageal body, above the LES. The clinical P-CSA plots revealed a horizontal loop shape.
Similar to cardiovascular analyses, we divided the P-CSA loop into phases. However, we aimed to offer a new mechanics-based understanding of the different phases in the P-CSA loop within the esophageal body. To obtain this understanding, we simulated flow inside a FLIP device placed inside the esophagus lumen. This gave us a closer look into the esophageal mechanism at each time instant. This information allowed us to develop a clear view of the phases of peristaltic contraction inside the FLIP. We concluded that the horizontal loop pattern within the esophagus body and the phases are a product of different mechanical conditions (muscle contraction, boundary condition, and geometric configuration) rather than any inherently different function of the muscle itself.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors upon request, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Northwestern University Institutional Review Board. The patients/participants provided their written informed consent to participate in this study and were paid for their participation.
Author contributions
GE wrote the main manuscript text and prepared all the figures; GE and SH wrote the code; NP, GE, and SH developed the mathematical formulation; DC contributed to data acquisition; GE, SH, DC, PK, JP, and NP interpreted results of calculations; PK, JP, and NP contributed to obtaining funding, critical revision of the manuscript and final approval.
Funding
This work was funded by the National Institute of Diabetes and Digestive and Kidney Diseases R01- DK079902 and P01-DK117824, and National Science Foundation Grants OAC 1450374 and OAC 1931372. Their support is greatly appreciated.
Conflict of interest
PK and JP hold shared intellectual property rights and ownership surrounding FLIP panometry systems, methods, and apparatus with Medtronic Inc. DC: Medtronic (Speaking, Consulting) PK: Ironwood (Consulting), Reckitt (Consulting), Johnson \& Johnson (Consulting) JP: Crospon, Inc. (stock options), Given Imaging (Consultant, Grant, Speaking), Sandhill Scientific (Consulting, Speaking), Takeda (Speaking), Astra Zeneca (Speaking), Medtronic (Speaking, Consulting), Torax (Speaking, Consulting), Ironwood (Consulting), Impleo (Grant).
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Acharya S., Halder S., Carlson D. A., Kou W., Kahrilas P. J., Pandolfino J. E., et al. (2021a). Assessment of esophageal body peristaltic work using functional lumen imaging probe panometry. Am. J. Physiology-Gastrointestinal Liver Physiology 320, G217–G226. doi:10.1152/ajpgi.00324.2020
Acharya S., Halder S., Carlson D. A., Kou W., Kahrilas P. J., Pandolfino J. E., et al. (2021b). Estimation of mechanical work done to open the esophagogastric junction using functional lumen imaging probe panometry. Am. J. Physiology-Gastrointestinal Liver Physiology 320, G780–G790. doi:10.1152/ajpgi.00032.2021
Acharya S., Halder S., Kou W., Kahrilas P. J., Pandolfino J. E., Patankar N. A. (2022). A fully resolved multiphysics model of gastric peristalsis and bolus emptying in the upper gastrointestinal tract. Comput. Biol. Med. 143, 104948. doi:10.1016/j.compbiomed.2021.104948
Acharya S., Kou W., Halder S., Carlson D. A., Kahrilas P. J., Pandolfino J. E., et al. (2021c). Pumping patterns and work done during peristalsis in finite-length elastic tubes. J. Biomechanical Eng. 143, 071001. doi:10.1115/1.4050284
Aziz Q., Fass R., Gyawali C. P., Miwa H., Pandolfino J. E., Zerbib F. (2016). Esophageal disorders. Gastroenterology 150, 1368–1379. doi:10.1053/j.gastro.2016.02.012
Bringley T. T., Childress S., Vandenberghe N., Zhang J. (2008). An experimental investigation and a simple model of a valveless pump. Phys. Fluids 20, 033602. doi:10.1063/1.2890790
Carlson D. A., Gyawali C. P., Khan A., YadlapatiChen J., Chokshi R. V., et al. (2021). Classifying esophageal motility by flip panometry: A study of 722 subjects with manometry. Am. J. Gastroenterology 116, 2357–2366. doi:10.14309/ajg.0000000000001532
Carlson D. A., Kahrilas P. J., Lin Z., Hirano I., Gonsalves N., Listernick Z., et al. (2016). Evaluation of esophageal motility utilizing the functional lumen imaging probe. Am. J. Gastroenterology 111, 1726–1735. doi:10.1038/ajg.2016.454
Carlson D. A., Kou W., Masihi M., Acharya S., Baumann A. J., Donnan E. N., et al. (2020). Repetitive antegrade contraction: A novel response to sustained esophageal distension is modulated by cholinergic influence. Am. J. Physiology-Gastrointestinal Liver Physiology 319, G696–G702. doi:10.1152/ajpgi.00305.2020
Carlson D. A., Lin Z., Kahrilas P. J., Sternbach J., Donnan E. N., Friesen L., et al. (2015). The functional lumen imaging probe detects esophageal contractility not observed with manometry in patients with achalasia. Gastroenterology 149, 1742–1751. doi:10.1053/j.gastro.2015.08.005
Clarke J. O., Ahuja N. K., Fernandez-Becker N. Q., Gregersen H., Kamal A. N., Khan A., et al. (2020). The functional lumen imaging probe in gastrointestinal disorders: The past, present, and future. Ann. N. Y. Acad. Sci. 1482, 16–25. doi:10.1111/nyas.14463
Crist J., Gidda J. S., Goyal R. K. (1984). Intramural mechanism of esophageal peristalsis: Roles of cholinergic and noncholinergic nerves. Proc. Natl. Acad. Sci. U. S. A. 81, 3595–3599. doi:10.1073/pnas.81.11.3595
Elisha G., Acharya S., Halder S., Carlson D. A., Kou W., Kahrilas P. J., et al. (2022a). Peristaltic regimes in esophageal transport. Biomechanics Model. Mechanobiol. doi:10.1007/s10237-022-01625-x
Elisha G., Halder S., Acharya S., Carlson D. A., Kou W., Kahrilas P. J., et al. (2021). A mechanics-based perspective on the function of the esophagogastric junction during functional luminal imaging probe manometry. (preprint).
Elisha G., Halder S., Carlson D. A., Kou W., Kahrilas P. J., Pandolfino J. E., et al. (2022b). A mechanics–based perspective on the function of human sphincters during functional luminal imaging probe manometry. (preprint). doi:10.48550/arXiv.2209.08419
Goyal R. K., Chaudhury A. (2008). Physiology of normal esophageal motility. J. Clin. gastroenterology 42, 610–619. doi:10.1097/MCG.0b013e31816b444d
Gregersen H., Liao D., Pedersen J., Drewes A. (2007). A new method for evaluation of intestinal muscle contraction properties: Studies in normal subjects and in patients with systemic sclerosis. Neurogastroenterol. Motil. 19, 11–19. doi:10.1111/j.1365-2982.2006.00837.x
Gregersen H., Lo K. M. (2018). What is the future of impedance planimetry in gastroenterology? Neurogastroenterol. Motil. 24, 166–181. doi:10.5056/jnm18013
Gregersen H., Villadsen G. E., Liao D. (2011). Mechanical characteristics of distension-evoked peristaltic contractions in the esophagus of systemic sclerosis patients. Dig. Dis. Sci. 56, 3559–3568. doi:10.1007/s10620-011-1777-9
Halder S., Acharya S., Kou W., Campagna R. A. J., Triggs J. R., Carlson D. A., et al. (2022a). Myotomy technique and esophageal contractility impact blown-out myotomy formation in achalasia: An in silico investigation. Am. J. Physiology-Gastrointestinal Liver Physiology 322, G500–G512. doi:10.1152/ajpgi.00281.2021
Halder S., Acharya S., Kou W., Kahrilas P. J., Pandolfino J. E., Patankar N. A. (2021). Mechanics informed fluoroscopy of esophageal transport. Biomechanics Model. Mechanobiol. 20, 925–940. doi:10.1007/s10237-021-01420-0
Halder S., Johnson E. M., Yamasaki J., Kahrilas P. J., Markl M., Pandolfino J. E., et al. (2022b). Mri-mech: Mechanics-informed mri to estimate esophageal health. (preprint). doi:10.48550/arXiv.2209.07492
Halder S., Yamasaki J., Acharya S., Kou W., Elisha G., Carlson D. A., et al. (2022c). Virtual disease landscape using mechanics-informed machine learning: Application to esophageal disorders. Artif. Intell. Med. 134, 102435. doi:10.1016/j.artmed.2022.102435
Jain A. S., Carlson D. A., Triggs J., Tye M., Kou W., Campagna R., et al. (2019). Esophagogastric junction distensibility on functional lumen imaging probe topography predicts treatment response in achalasia—Anatomy matters. Am. J. Gastroenterology 114, 1455–1463. doi:10.14309/ajg.0000000000000137
Kahrilas P. J. (2008). Clinical practice. Gastroesophageal reflux disease. N. Engl. J. Med. 359, 1700–1707. doi:10.1056/NEJMcp0804684
Kou W., Griffith B. E., Pandolfino J. E., Kahrilas P. J., Patankar N. A. (2017a). A continuum mechanics-based musculo-mechanical model for esophageal transport. J. Comput. Phys. 348, 433–459. doi:10.1016/j.jcp.2017.07.025
Kou W., Pandolfino J. E., Kahrilas P. J., Patankar N. A. (2017b). Could the peristaltic transition zone be caused by non-uniform esophageal muscle fiber architecture? A simulation study. Neurogastroenterol. Motil. 29, e13022. doi:10.1111/nmo.13022
Kou W., Pandolfino J. E., Kahrilas P. J., Patankar N. A. (2017c). Simulation studies of the role of esophageal mucosa in bolus transport. Biomechanics Model. Mechanobiol. 16, 1001–1009. doi:10.1007/s10237-016-0867-1
Kou W., Pandolfino J. E., Kahrilas P. J., Patankar N. A. (2018). Studies of abnormalities of the lower esophageal sphincter during esophageal emptying based on a fully coupled bolus–esophageal–gastric model. Biomechanics Model. Mechanobiol. 17, 1069–1082. doi:10.1007/s10237-018-1014-y
Kwiatek M. A., Hirano I., Kahrilas P. J., Rothe J., Luger D., Pandolfino J. E. (2011). Mechanical properties of the esophagus in eosinophilic esophagitis. Gastroenterology 140, 82–90. doi:10.1053/j.gastro.2010.09.037
Leibbrandt R. E., Dinning P. G., Costa M., Cock C., Wiklendt L., Wang G., et al. (2016). Characterization of esophageal physiology using mechanical state analysis. Front. Syst. Neurosci. 10, 10. doi:10.3389/fnsys.2016.00010
Liao D., Krarup A. L., Lundager F. H., Drewes A. M., Gregersen H. (2014). Quantitative differences between primary and secondary peristaltic contractions of the esophagus. Dig. Dis. Sci. 59, 1810–1816. doi:10.1007/s10620-014-3070-1
Liao D., Mark E. B., Zhao J., Drewes A. M., Brock C. (2020). Modeling and measurement of the mechano-physiological function of the gastrointestinal organs. Physiol. Meas. 41, 114004. doi:10.1088/1361-6579/abc40a
Liao D., Villadsen G., Gregersen H. (2013). Distension-evoked motility analysis in human esophagus. Neurogastroenterol. Motil. 25, 407–412. doi:10.1111/nmo.12081
Lin Z., Kahrilas P. J., Xiao Y., Nicodeme F., Gonsalves N., Hirano I., et al. (2013). Functional luminal imaging probe topography: An improved method for characterizing esophageal distensibility in eosinophilic esophagitis. Ther. Adv. Gastroenterology 6, 97–107. doi:10.1177/1756283X12470017
Lottrup C., Gregersen H., Liao D., Fynne L., Frøkjær J. B., Krogh K., et al. (2015). Functional lumen imaging of the gastrointestinal tract. J. gastroenterology 50, 1005–1016. doi:10.1007/s00535-015-1087-7
Manopoulos C. G., Mathioulakis D. S., Tsangaris S. G. (2006). One-dimensional model of valveless pumping in a closed loop and a numerical solution. Phys. Fluids 18, 017106. doi:10.1063/1.2165780
Moonen A., Boeckxstaens G. (2014). Measuring mechanical properties of the esophageal wall using impedance planimetry. Gastrointest. Endosc. Clin. 24, 607–618. doi:10.1016/j.giec.2014.06.001
Omari T. I., Wiklendt L., Dinning P., Costa M., Rommel N., Cock C. (2015). Upper esophageal sphincter mechanical states analysis: A novel methodology to describe ues relaxation and opening. Front. Syst. Neurosci. 8, 241. doi:10.3389/fnsys.2014.00241
Ottesen J. (2003). Valveless pumping in a fluid-filled closed elastic tube-system: One-dimensional theory with experimental validation. J. Math. Biol. 46, 309–332. doi:10.1007/s00285-002-0179-1
Peskin C. S. (2002). The immersed boundary method. Acta Numer. 11, 479–517. doi:10.1017/S0962492902000077
Regan J., Walshe M., Rommel N., Tack J., McMahon B. P. (2012). New measures of upper esophageal sphincter distensibility and opening patterns during swallowing in healthy subjects using EndoFLIP. Neurogastroenterol. Motil. 25, e25–e34. doi:10.1111/nmo.12041
Shaheen N. J., Hansen R. A., Morgan D. R., Gangarosa L. M., Ringel Y., Thiny M. T., et al. (2010). The burden of gastrointestinal and liver diseases. Am. J. Gastroenterology 101, 2128–2138. doi:10.1111/j.1572-0241.2006.00723.x
Villadsen G., Storkholm J., Zachariae H., Hendel L., Bendtsen F., Gregersen H. (2001). Oesophageal pressure–cross-sectional area distributions and secondary peristalsis in relation to subclassification of systemic sclerosis. Neurogastroenterol. Motil. 13, 199–210. doi:10.1046/j.1365-2982.2001.00259.x
Whittaker R. J., Heil M., Jensen O. E., Waters S. L. (2010). A rational derivation of a tube law from shell theory. Q. J. Mech. Appl. Math. 63, 465–496. doi:10.1093/qjmam/hbq020
Keywords: sphincter, esophagus, peristalsis, pressure-area hysteresis loop, functional luminal imaging probe, mechanical states
Citation: Elisha G, Halder S, Carlson DA, Kahrilas PJ, Pandolfino JE and Patankar NA (2023) A mechanics-based perspective on the pressure-cross-sectional area loop within the esophageal body. Front. Physiol. 13:1066351. doi: 10.3389/fphys.2022.1066351
Received: 10 October 2022; Accepted: 12 December 2022;
Published: 09 January 2023.
Edited by:
Luca Dede’, Politecnico di Milano, ItalyReviewed by:
Francesco Regazzoni, Politecnico di Milano, ItalyGreger Lindberg, Karolinska Institutet (KI), Sweden
Copyright © 2023 Elisha, Halder, Carlson, Kahrilas, Pandolfino and Patankar. 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: Neelesh A. Patankar, n-patankar@northwestern.edu