DECREASED VARIATION IN ESOPHAGEAL ACID CONCENTRATION IN DIFFERENT PHENOTYPES OF SYMPTOMATIC GASTROESOPHAGEAL REFLUX DISEASE =========================================================================================================================== * Jerry D. Gardner ## ABSTRACT **Background** Homeostatic physiologic systems are regulated to maintain a particular variable such as blood pressure within a restricted range with relatively low variation, and this variation is increased with disease. Fractal physiologic systems, on the other hand, are characterized by wide, complex variation of variables such as heart rate or the walking stride interval, and this variation is decreased with disease. The present report examines time-series recordings of esophageal pH from normal subjects and different GERD phenotypes to measure the distributions of esophageal pH values and the distributions of changes in esophageal acid concentrations. **Methods** Using Lyon consensus definitions of symptomatic GERD phenotypes, I analyzed 24-hour esophageal pH recordings from normal subjects (n=20), Functional Heartburn subjects (n=20), Reflux Hypersensitivity subjects (n=20), and nonerosive esophageal reflux disease (NERD) subjects (n=20). For each subject I calculated the distribution of pH values as well as the distribution of changes in esophageal acid concentrations. **Results** Esophageal pH values have a power law distribution in both normal and symptomatic GERD phenotypes, and esophageal acid concentrations vary over four orders of magnitude in each group. The variation in esophageal acid concentrations decreased progressively from normal subjects to Functional Heartburn subjects to Reflux Hypersensitivity subjects to NERD subjects. **Conclusions** The decreased variation in esophageal acid concentration in symptomatic GERD phenotypes represents changes associated with disease in a fractal physiologic system that shares characteristic features with other fractal systems where disease states have been associated with decreased variation of heart rate or walking stride interval. Keywords * homeostasis * fractal * power law * cumulative acidity ## INTRODUCTION Frequently, disease or ageing are accompanied by changes in one or more variables in physiologic systems (1-4). Homeostatic physiologic systems are regulated to maintain a particular variable such as blood pressure within a restricted range with relatively low variation, and this variation is increased with ageing or disease (1, 2). Fractal physiologic systems, on the other hand, are characterized by wide, complex variation of variables such as heart rate or the walking stride interval, and this variation is decreased with ageing or disease (3, 4). The Lyon Consensus Conference (5, 6) proposed criteria for the clinical diagnosis of three different phenotypes of gastroesophageal reflux disease (GERD): nonerosive gastroesophageal reflux disease (NERD), Reflux Hypersensitivity, and Functional Heartburn. Previous analyses, before recognition of different phenotypes of GERD, found that time-series of esophageal pH showed a fractal pattern in normal and GERD subjects (7). In other analyses esophageal pH values had a power law distribution in normal and GERD subjects (8). The present report examines time-series recordings of esophageal pH from normal subjects and different GERD phenotypes to measure the distributions of esophageal pH values as well as the distributions of changes in esophageal acid concentrations. ## SUBJECTS Patients were identified by exploring the electronic database at the Royal London Hospital GI Physiology Unit that contains clinically indicated impedance-pH recordings (Sandhill Scientific, Highlands Ranch, CO) from patients with typical symptoms of gastroesophageal reflux. Using Lyon consensus definitions of symptomatic GERD phenotypes,(5, 6), I selected 24-hour esophageal pH recordings from normal subjects (n=20), Functional Heartburn subjects (n=20), Reflux Hypersensitivity subjects (n=20), and nonerosive esophageal reflux disease (NERD) subjects (n=20). All subjects had a normal upper gastrointestinal endoscopy at the time of the impedance-pH study. The esophageal pH recordings from one normal subject were technically unsatisfactory and were omitted from the present analyses. For this retrospective analysis of clinically indicated tests with no identifiable patient data, the Stanford University Institutional Review Board determined that this research does not involve human subjects as defined in 45 CFR 46.102(f) or 21 CFR 50.3 (g) (9). Values from impedance-pH testing from subjects for the present analyses have been published previously (10-12). Normal subjects, Functional Heartburn subjects and Reflux Hypersensitivity subjects all had normal esophageal acid exposure time (AET) with esophageal pH <4 for less than 4% of the 24-hour esophageal pH recording. NERD subjects had increased esophageal AET of pH <4 for greater than 6% of the 24-hour esophageal pH recording. Reflux Hypersensitivity subjects had a positive association of symptoms with reflux episodes (13, 14), whereas Functional Heartburn subjects had no association of symptoms with reflux episodes. Normal subjects had no symptoms during the impedance-pH testing. ## METHODS The impedance-pH catheter was inserted with an esophageal pH sensor positioned 5cm above the upper border of the lower esophageal sphincter. Six impedance channels were positioned 3, 5, 7, 9, 15 and 17cm above the upper border of the lower esophageal sphincter, respectively. Software provided by Sandhill to process pH recordings automatically adjusts all pH values for the difference between the calibration temperature of 25C and the recording temperature of 37C. Software provided by Sandhill was also used to export pH data for every 4th second of the recording to an Excel file. All pH values below 0.5 were replaced by 0.5 and all pH values above 7.5 were replaced by 7.5. To determine cumulative acid concentration for a subject, all pH values were converted to acid concentration in mmol/L Cumulative acid concentration was calculated as the sum of all values for acid concentration for that subject, and each value of acid concentration was then expressed as a percentage of the cumulative acid concentration. The distribution of the change in sequential values of acid concentration was calculated for esophageal acidity for each subject. The frequency distributions of changes in acid concentration for each group of subjects were calculated as the means of the values from all subjects in the group for each bin of the distribution. Expressing acid concentration as a percentage of cumulative acid concentration for a given subject makes it possible to examine differences in the distributions of changes in acid concentration that do not depend on the magnitude of the acid concentration per se. Others (1) have calculated the change in values in a time-series as a percentage difference from the mean of all values. Curve fitting and statistical analyses were performed using GraphPad Prism 9.4.1 software. Because the present analyses were exploratory, P-values were not adjusted for multiple comparisons. ## RESULTS Figure 1 shows a linear relationship between the values of esophageal pH and the logarithm of the frequency of the values for Normal subjects and each GERD phenotype. Results in Table 1 show that each slope of the lines in Figure 1 is significantly different from zero and accompanying values of R2 are at least 0.90. The linear relationships illustrated in Figure 1 are described by a power law in that the frequency of what is being measured is a negative exponential function of the magnitude of what is being measured, and the value of the exponential is given by the slope. The term “power law” is used because the frequency of a particular pH value is a power of the pH value. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/04/03/2023.03.31.23288021/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2023/04/03/2023.03.31.23288021/F1) Figure 1. Distributions of values of esophageal pH in normal subjects and GERD phenotypes. Values given on the Y-axis are mean frequencies for the bin indicated on the X-axis. The solid lines are from linear, least-squares regression analyses. Abbreviations are NL – normal subjects; FX HB – Functional Heartburn subjects; RFX HYP – Reflux Hypersensitivity subjects; NERD – Nonerosive Reflux Disease subjects. The number of subjects in each group is given in parentheses after each abbreviation. View this table: [Table 1.](http://medrxiv.org/content/early/2023/04/03/2023.03.31.23288021/T1) Table 1. RESULTS FROM LINEAR REGRESSION ANALYSES OF THE DISTRIBUTIONS OF ESOPHAGEAL pH VALUES. Table 2 shows results from pairwise comparisons of the values of the slopes given in Table 1. The slope from normal subjects was significantly higher than that from Reflux Hypersensitivity subjects and NERD subjects. It seemed possible that the lack of statistical significance comparing the slope from Functional Heartburn subjects to that from normal subjects might be due to an under-powered sample size. The slope from NERD subjects was significantly lower than that from normal subjects, from Functional Heartburn subjects, and from Reflux Hypersensitivity subjects. The slope from Functional Heartburn subjects was not significantly different from that from Reflux Hypersensitivity subjects. View this table: [Table 2.](http://medrxiv.org/content/early/2023/04/03/2023.03.31.23288021/T2) Table 2. COMPARISON OF THE SLOPES OF DISTRIBUTIONS OF ESOPHAGEAL pH VALUES IN NORMAL SUBJECTS AND GERD PHENOTYPES. Figure 2 shows that the distribution of changes in cumulative esophageal acidity is biphasic for each GERD phenotype and normal subjects, and each distribution spans at least 4 orders of magnitude. The distribution of changes was significantly better fit by a 6th order polynomial than by a 5th order polynomial (P<0.0001 by F-test). A 7th order polynomial did not give a significantly better fit than a 6th order polynomial by an F-test. The polynomial model has no physiologic significance. It is simply a model that creates a curve that comes close to the data points and makes it possible to test the data for statistical differences. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/04/03/2023.03.31.23288021/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2023/04/03/2023.03.31.23288021/F2) Figure 2. Distributions of change in values for cumulative esophageal acidity for GERD phenotypes and normal subjects. Values given on the Y-axis are mean frequencies for the bin indicated on the X-axis. Abbreviations are CUM – cumulative; ESO – esophageal. Table 3 gives results from statistical comparisons of the data for distributions of changes in esophageal acid concentration in Figure 2 and indicates that all pairwise comparisons are significantly different at P<0.0001 by an F-test. Furthermore, the variation in esophageal acid concentrations is indicated by the range of values for change in esophageal acid concentration and is greatest in normal subjects, less in Functional Heartburn, still less in Reflux Hypersensitivity and least in NERD. View this table: [Table 3.](http://medrxiv.org/content/early/2023/04/03/2023.03.31.23288021/T3) Table 3. COMPARISON OF THE DISTRIBUTIONS OF CHANGES IN ESOPHAGEAL ACIDITY IN GERD PHENOTYPES USING FITS TO A 6TH ORDER POLYNOMIAL. ## DISCUSSION The present analyses of time-series of esophageal pH measurements show that esophageal pH values have a power law distribution in both normal and symptomatic GERD phenotypes, and that esophageal acid concentrations vary over four orders of magnitude in each group. Also, the variation in esophageal acid concentration decreased progressively from normal subjects to Functional Heartburn subjects to Reflux Hypersensitivity subjects to NERD subjects. Thus, the different GERD phenotypes based on clinical characteristics identified by the Lyon Consensus Conference (5, 6) may actually share a common pathophysiology of varying severity. The decreased variation in esophageal acid concentrations in symptomatic GERD phenotypes represents changes associated with disease in a fractal physiologic system that is characterized by wide, complex variation under normal circumstances. Thus, GERD phenotypes share characteristic features with other fractal systems where disease states have been associated with decreased variation of heart rate or walking stride interval (3, 4). All GERD phenotypes have symptoms of gastroesophageal reflux such as heartburn, regurgitation, or chest pain (5, 6) even though Functional Heartburn subjects and Reflux Hypersensitivity subjects have normal esophageal acid exposure times. It seems possible that regardless of the value of esophageal acid exposure time, the decreased variation in esophageal acid concentration, which can result in runs of self-similar esophageal acid concentrations, may represent a pathophysiologic signal to esophageal mucosa that triggers a symptom. Previous analyses of the same subjects used for the present study found that esophageal acid sensitivity appears to oscillate in each GERD phenotype, and for a given value of esophageal acid sensitivity, Reflux Hypersensitivity subjects have significantly more sequential symptoms associated with this sensitivity than do Functional Heartburn subjects (11). This difference between Functional Heartburn subjects and Reflux Hypersensitivity subjects might indicate that the lower variation of values of esophageal acid concentrations in Reflux Hypersensitivity subjects can cause an increase in GERD symptoms. ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## Acknowledgement I am grateful to Dr. Daniel Sifrim, Director of Upper GI Physiology Unit, Royal London Hospital for providing the impedance-pH records and for stimulating, helpful discussions. ## Footnotes * The work described in this manuscript has not been published and is not being considered for publication elsewhere. * **Conflict of Interest Statement:** Dr. Gardner is President of Science for Organizations ([www.scifororg.com](http://www.scifororg.com)), a company that provides consulting services to biotechnology and pharmaceutical companies. These activities have no conflict of interest with the analyses in the present paper. * **Ethics Approval:** For this retrospective analysis of clinically indicated tests with no identifiable patient data, the Stanford University Institutional Review Board determined that this research does not involve human subjects as defined in 45 CFR 46.102(f) or 21 CFR 50.3 (g). * All data produced in the present study are available upon reasonable request to the authors. * Received March 31, 2023. * Revision received March 31, 2023. * Accepted April 3, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. 1.Fossion R, River AL, Estanol B. A physicist’s view of homeostasis: how time series of continuous monitoring reflect the function of physiological variables in regulatory mechanisms. Physiol Meas. 2018, 39, 084007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/1361-6579/aad8db&link_type=DOI) 2. 2.Parati G, Ochoa JE, Lombardi c, Bilo G. Assessment and management of blood- pressure variability. Nature Reviews|Cardiology. 2013: 10:143–153. 3. 3.Goldberger AL, Amaral LAN, Hausdorff JM, Ivanov PC, Peng C-K, Stanley HE. Fractal dynamics in physiology: Alterations with disease and ageing. Proc Natl Acad Sci. 2002; 99 (Suppl 1): 2466–2472. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxNToiOTkvc3VwcGxfMS8yNDY2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDQvMDMvMjAyMy4wMy4zMS4yMzI4ODAyMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 4. 4.Guidelines. Heart rate variability. Eur Heart Journal 1996; 17:354–381. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/oxfordjournals.eurheartj.a014868&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8737210&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F04%2F03%2F2023.03.31.23288021.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UF54400011&link_type=ISI) 5. 5.Gyawali CP, Kahrilas PJ, Savarino E, Zerbib F, Mion F, Smout Ajpm, Vaezi M, Sifrim D, Fox MR, Vela MF, Tutuian R, Tack J, Bredenoord AJ, Pandofilno J, Roman S. Modern diagnosis of GERD: the Lyon Consensus. Gut 2018; 67:1351–1362. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjk6IjY3LzcvMTM1MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA0LzAzLzIwMjMuMDMuMzEuMjMyODgwMjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 6. 6.Ghisa M, Barberio B, Savarino V, Marabotto E, Ribolsi M, Bodini G, Zingone F, Frazzoni M, Savarino E. The Lyon Consensus: Does It Differ From the Previous Ones? J Neurogastroenterol Motil 2020;26:311–321. 7. 7.Gardner JD, Young W, Sloan S, Robinson M, Miner PB Jr. The fractal nature of human gastro-oesophageal reflux. Aliment Pharmacol Ther 2005; 22: 823–830. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2036.2005.02665.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16225491&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F04%2F03%2F2023.03.31.23288021.atom) 8. 8.Gardner JD, Sloan S, Robinson M, Miner PB. Oesophageal pH has a power-law distribution in control and gastro-oesophageal reflux disease subjects. Aliment Pharmacol Ther 2004; 20: 1373–1379. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2036.2004.02278.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15606400&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F04%2F03%2F2023.03.31.23288021.atom) 9. 9.Exceptions to Institutional Review Board requirements for research. United States Code of Federal Regulations, Title 45, Part 46, Subpart A, Criterion 4. 10. 10.Gardner JD. The relationship between esophageal acidity and symptom frequency in symptomatic nonerosive gastroesophageal reflux disease. Physiological Reports, 10, e15442. 11. 11.Gardner JD. (2022). Oscillating esophageal acid sensitivity in symptomatic reflux hypersensitivity and functional heartburn. Qeios. doi:10.32388/IJUE1J.2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.32388/IJUE1J.2&link_type=DOI) 12. 12.Sifrim D, Roman S, Savarino E, Bor S, Bredenoord AJ, Castell D, Cicala M, De Bortoli N, Frazzoni M, Gonlachanvit S, Iwakiri K, Kawamura O, Krarup A, Lee YY, Ngiu CS, Ndebia E, Patcharatraku T, Pauwels A, De la Serna JP, Ramos R, Remes-Troche Jm, Ribolsi M, Sammon A, Simren M,Tack J, Tutuian R, Valdovinos M, Xiao Zerbib F, Gyawali CP. Normal values and regional differences in oesophageal impedance-pH metrics: a consensus analysis of impedance-pH studies from around the world. Gut 2021; 70:1441–1449. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjk6IjcwLzgvMTQ0MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA0LzAzLzIwMjMuMDMuMzEuMjMyODgwMjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 13. 13.Weiner GJ, Richter JE, Copper JB, Wu WC, Castell DO. The symptom index: A clinically important parameter of 24-hour esophageal pH monitoring. Am J Gastroenterol 1988; 83:358–361. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3348191&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F04%2F03%2F2023.03.31.23288021.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1988M813300005&link_type=ISI) 14. 14.Weusten BL, Roelofs JM, Akkermans LM, Van Berge-Henegouwen GP, Smout AJ. The symptom-association probability: An improved method for symptom analysis of 24-hour esophageal pH data. Gastroenterology 1994; 107:1741–1745. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0016-5085(94)90815-X&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7958686&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F04%2F03%2F2023.03.31.23288021.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994PU54300021&link_type=ISI)