BACKGROUND Current breath tests for measurement of gastric emptying of solids are expensive, possibly inaccurate, and require cumbersome calculations.
AIMS We wished to validate a simplified solid gastric emptying test using a [13C]Spirulina platensisbreath test for accurate results relative to scintigraphy.
SUBJECTS Thirty healthy volunteers.
METHODS We measured gastric emptying of egg containing [13C]S platensis and99mTc sulphur colloid by breath13CO2 and scintigraphy over six hours. A generalised linear regression model was used to predict t1/2 and tLAG by scintigraphy from breath13CO2 data. The model was cross validated and normative data calculated for a prepacked [13C]meal.
RESULTS Regression models using all breath data over six hours, for the first three hours, and for samples at 75, 90, and 180 minutes (“reduced model”) predicted t1/2 and tLAG values similar to scintigraphy (tLAG 43 (SD 12) min; t1/2 100 (20) min). Standard deviations of differences in t1/2 and tLAG between scintigraphy and the “reduced model” were both 10 minutes. Gastric t1/2 for the prepacked [13C]meal was 91 (15) min (10–90% range: 74–118).
CONCLUSION The [13C]S platensis breath test and a simple formula using breath 13CO2 at baseline, 90, and 180 minutes measured gastric emptying t1/2 for solids with results that were comparable with scintigraphy.
- stable isotope
- breath test
- gastric emptying
Statistics from Altmetric.com
Gastric emptying of solids is used to screen for gastric motor disorders. Two parameters that are clinically useful are t1/2 and the proportion emptied at four hours1-3; these reflect overall emptying functions. The time for the first 10% to empty (or lag time1) is a marker for the ability of the stomach to triturate solid food to a particle size that can be emptied (<2 mm4). These measurements (tLAG, t1/2) are best estimated by scintigraphy and the state of the art analysis uses a power exponential or a nonlinear model.1-5
Stable isotope breath tests were recently introduced to measure gastric emptying of solids because of their potential advantages6-9 which include performance of the test at the bedside or in the community; remote laboratory automated analysis; application at sites where gamma camera facilities are not available; and avoidance of radiation, facilitating research studies in children and pregnant women.
The lack of validity for interindividual comparisons is a significant pitfall10; inaccuracies of the test may result partly from the nonlinear models used to fit the cumulative breath13CO2 recoveries. In the breath test, the total gastric residual is unknown and is estimated by a formula where time is infinite and 13CO2 excretion is assumed to have reached a steady state. Estimates of gastric half-emptying time (t1/2) and lag phase duration (tLAG ) using nonlinear models are determined by the shape of the curve and are independent of endogenous CO2 production. The results, however, are strongly associated with the estimated parameter “m” or estimated gastric residual.11 In contrast, scintigraphy can be used to follow complete emptying of the stomach, an unequivocal end point. In our previous studies10 12 using the conventional nonlinear model, we determined that accurate determination of the parameter “m” in healthy subjects (that is a steady 13CO2 excretion state) required at least six hours of breath collection. We anticipated that patients with impaired gastric emptying may require an even longer sampling time, suggesting the breath test may become impractical.
From a review of simultaneous scintigraphic and13CO2 cumulative excretion curves performed in nearly 70 healthy volunteers in our laboratory, we noted that the “tail” of the cumulative excretion curve often did not reach steady state conditions; moreover, this portion of the curve corresponded to a time when all of the 99mTc isotope had emptied from the stomach. Thus optimizing the curve fit might represent a futile exercise unrelated to the process of emptying.
In the present study we aimed to develop a new mathematical model for breath data analysis that was independent of the parameter “m” and would offer accuracy with shorter (three hours or less) sampling times. This new mathematical model was developed while assessing the accuracy of a breath test using a 13C substrate consisting predominantly of the protein [13C]Spirulina platensis(described below). Subsequently, the model was further validated using the leave-one-out method. Finally, we performed the [13C]S platensis breath test in a simulated office based practice using a standardised prepacked [13C]biscuit commercial meal and compared the results with those of the [13C]S platensis egg meal.
Materials and methods
Thirty healthy volunteers, aged 18–53 years (mean 34 (SE 1.6); 15 males, 15 females; median body mass index 23 kg/m2 (range 18.7–39.1)), were recruited by public advertisement. None had previously undergone abdominal surgery other than appendicectomy. None was receiving medications. Females of childbearing potential had a negative pregnancy test within 48 hours of beginning the study. A clinical interview and physical examination were performed and none of the subjects was receiving medications that could influence gastrointestinal motility or autonomic function. Written informed consent was obtained before participation in the study which was approved by the Institutional Review Board and the Radiation Control Committee of the Mayo Clinic.
SUBSTRATE FOR 13CO2 BREATH TEST ([13C]S PLATENSIS)
S platensis is an edible blue-green alga used as a protein source in many parts of the world and sold as a food supplement in health food stores.13 14 It contains 50–60% protein, 30% starch, and 10% lipid.15 TheS platensis used in our study was grown as a pure monoculture in a closed hydroponics chamber in a medium containing inorganic salts and purged with pure 13CO2. The resultant cells were uniformly labelled with 13C; the natural level of 13C in the environment (about 1%) was increased to 99% in the substrate used. When metabolised, the proteins, carbohydrates, and lipids of the S platensis give rise to respiratory CO2 that is enriched in 13C. Because the contents of the algal cells are not freely diffusible, incorporation of labelledS platensis into a solid phase meal occurs easily and represents the digestive state of the solid food matrix within which it is contained.
SCINTIGRAPHY AND [13C]S PLATENSIS BREATH TEST METHODS
Radioscintigraphic marker and test meal
The test meal comprised eggs dosed with 200 mg of [13C]S platensis (Meretek Diagnostics, Inc., Houston, Texas, USA) and 0.5 mCi 99mTc sulphur colloid. The egg whites were mixed with 0.5 mCi99mTc sulphur colloid. The egg meal was placed on a slice of whole wheat bread and given with a glass of skimmed milk for a total caloric value of 220 kcal with a nutrient composition of 35% protein, 40% carbohydrate, 25% fat, and 2.6 g of fibre.
Scintigraphic and breath test methods
After an overnight fast, scintigraphic scanning with anterior and posterior cameras was performed with the patient standing. Imaging began at the start of the test meal and scans were obtained every 15 minutes for two hours and then every 30 minutes for six hours. Breath samples were taken at baseline before the meal and followed the same time schedule as the scintigraphic technique. Breath samples were collected and stored in duplicate in glass vacutainers using a straw to blow into the bottom of the tube. After restoppering the tubes, the13C breath content was determined in a centralised laboratory (Dr P D Klein, Meretek Diagnostics) by isotope ratio mass spectrometry.16 17
Standard meal study using [13C]S platensis biscuit
About seven months (5–9 months) after the [13C]S platensis in egg meal studies, another study using a standardised meal including a [13C]S platensis biscuit (Meretek Diagnostics, Inc.) was performed in 27 of the 30 healthy subjects. After an overnight fast they ingested the standardised meal consisting of a 60 g rye roll (160 calories), 30 g cream cheese (90 calories), and 120 ml white grape juice (80 calories). The rye roll contained 200 mg of [13C]S platensis. Breath samples were collected at baseline (before meal) and at 75, 90, and 180 minutes, and gastric t1/2 and tLAG were estimated using the same mathematical method.
A region of interest was drawn around the stomach on the anterior and posterior images for each time frame. To correct for tissue attenuation, the counts in each region were multiplied together and the square root of the product was taken to obtain the geometric mean. Gastric emptying was summarised using a power exponential model1: where propt is the proportion remaining in the stomach at time t. The index k is an expression of the instantaneous slope of the curve and b is an index for the shape of the curve (b=1 implies a simple exponential emptying model). The parametersk and b were estimated using the nonlinear least squares (NLIN) procedure in the SAS software package.18 The duration of the lag phase was defined as the time when 10% of the test meal was emptied from the stomach. The gastric half-emptying time (t1/2S) and lag phase duration (tLAGS) were derived after estimatingk and b for each subject, solving the above formula for timet, and using 0.5 and 0.1, respectively, for propt.
[13C]S platensis breath test data
13C enrichment determined by isotope ratio mass spectrometry was expressed as the delta per mil difference between the13CO2 /12CO2 ratio of the sample and standard. 13CO2 in the breath samples was calculated as atom percent excess (percentage13CO2 in total CO2 above background13CO2 concentration in nature) as in previous studies.18-20 The results of the breath test were presented as a percentage of 13CO2 recovery per hour and as a cumulative value over six hours.13CO2 production was assumed to be 300 mmol per square metre of body surface area per hour as published by Ghooset al.6 Body surface area was calculated using a validated weight-height formula.
To calculate the quantity of 13C appearing in breath per unit time, delta over baseline (DOB) was used: 13C (μM/min)=DOB× 0.0112372×CO2 production, where 0.0112372 is the isotopic abundance of the limestone standard, Pee Dee Belemnite, and CO2 production was corrected for individual basal metabolic rates according to age (<3, 3–10, 10–18, 18–30, 30–60, >60 years) and gender as determined by Schofield21
Generalised linear regression model to predict scintigraphic estimates of gastric half-emptying time and lag phase using breath test values
A generalised linear multiple regression analysis19was used to develop a model for predicting the scintigraphic gastric emptying estimates for t1/2 and tLAG using13C values from breath samples obtained over the first three hours after the meal. In this generalised linear model approach, t1/2 and tLAG values were assumed to follow a gamma distribution (positively skewed distribution) to account for a common situation in clinical practice where data from patients with gastroparesis, or dumping syndrome, often lead to a skewed distribution of gastric emptying t1/2 and tLAG values. The link function used to relate the linear predictor (weighted combination of 13C values at the multiple time points) to scintigraphically observed t1/2 (or tLAG) values was the inverse (reciprocal) function. A preliminary analysis of the (reciprocally) transformed t1/2 and tLAGvalues using a simple stepwise multiple regression model (maximumr 2 method) provided candidate subsets of breath sample times to consider as “reduced models” (that is, fewer time points) in the generalised regression (gamma) model analysis.
As it is possible that evaluation of models to estimate tLAG and t1/2 may give an overly optimistic assessment when based on the same data,20 an internal cross validation analysis was used to assess the results of the “reduced models”. Thus a “leave-one-out” approach20 was used to obtain the breath sample prediction for each subject. In this analysis, each individual subject is left out of the calculations to estimate the regression parameters of the linear predictor for the gamma regression model; then the predicted value for the subject omitted from the model (or left out) is computed using the regression model results based on the remaining subjects included (or left in) the model. This is done repeatedly with each subject in turn being left out and used as a (quasi) external validation observation.
Comparisons of the adequacy of fit of the models between the full six hour data compared with the three hour data or the “reduced models” were based on the difference in log likelihood values for each model. Similarly, the adequacy of fit for the three hour data and the “reduced models” were compared.
In addition, estimates of t1/2 and tLAG by the full and “reduced models” were compared using a (non-parametric) Wilcoxon signed rank test of the differences resulting from different estimation methods (for example a paired comparison of the absolute values of (scintigraphic t1/2 minus breath sample estimated t1/2)). The correlation between scintigraphically determined values (observed values) and predicted values (for the various estimation methods) were computed, and the distribution of simple residuals (observed-predicted) summarised.
The distributions of t1/2 and tLAG estimates from the breath test using either the egg meal or the biscuit meal in 27 healthy subjects were compared.
COMPARISON OF GASTRIC EMPTYING INDICES BETWEEN SCINTIGRAPHY AND THE BREATH TEST USING A GENERALISED LINEAR REGRESSION MODEL PREDICTING SCINTIGRAPHIC RESULTS
Gastric emptying by scintigraphy yielded mean values of 43 (SD 12) minutes for t LAGS and 100 (20) minutes for t 1/2S. The generalised linear regression models predicting scintigraphic tLAG based on the entire six hour and the first three hour breath test data for 13C breath excretion did not differ significantly (r=0.79 and 0.72, respectively). Twice the difference in log likelihood values (which is approximated by a chi-square distribution) was 5.0 on 6 degrees of freedom (p=NS). The standard deviation of the difference between the observed tLAG and the predicted values computed using the three hour breath test data was eight minutes (range −22 to 19) (fig1A). The generalised linear models predicting scintigraphic t1/2 based on the entire six hour and the first three hour breath test data for 13C breath excretion did not differ significantly (r=0.91 and 0.86, respectively; difference in log likelihood, 7.5 on 6 degrees of freedom; NS). The standard deviation of the differences between t1/2 from scintigraphy and the predicted value computed using the three hour breath test data was 10 minutes (range −23 to 17) (fig 1B).
The generalised linear regression models identified three time points as being significant predictors of tLAG (75, 180 minutes) and t1/2 (90, 180 minutes). The model predicting scintigraphic tLAG using two data points (reduced model: 75, 180 minutes) was less accurate (r=0.61v full modelr=0.72) (fig 1A, C), although the difference in log likelihood was 4.6 on 8 degrees of freedom (NS). The standard deviation of differences in tLAG between the scintigraphy value and the predicted value computed using the “reduced model” was also 10 minutes (range −21 to 21). The generalised regression model predicting scintigraphic t1/2 based on just two data points (“reduced model”: 90, 180 minutes) for13CO2 breath excretion did not differ significantly from the full model (r=0.86v full modelr=0.80; log likelihood 5.8 on 8 degrees of freedom; NS) (fig 1B, D). The standard deviation of the differences in t1/2 from scintigraphy and the predicted value computed using the “reduced model” was 10 minutes (range −27 to 29).
A summary of the mathematical models using data from three breath test samples (75, 90, 180 minutes) is: tLAG=1/LPLAG, where LPLAG=0.0250+0.006313C75−0.003213C180; and for t 1/2: t1/2=1/LP1/2, where LP1/2=0.0097+ 0.002113C90−0.001213C180.
In these formulas, tLAG and t1/2 are predicted tLAG and t1/2 using breath test concentrations in μmol/min at the specific times 75, 90, and 180 minutes; LPLAG and LP1/2 are the linear predictors (that is, weighted sums of breath test concentrations); and13Ci is the excreted quantity of13C (μmol/min) for timet i=75, 90, or 180 minutes.
GASTRIC EMPTYING PARAMETERS IN HEALTHY SUBJECTS USING BREATH TEST AND SCINTIGRAPHY
Figure 2 summarises gastric emptying indices for 30 healthy volunteers from scintigraphy and various breath test results (a generalised linear regression model using the three hour data and a generalised linear regression “reduced model” using only two time points). The predicted tLAG and t1/2 values from the generalised linear regression (“reduced models”) were similar to the scintigraphic values (fig 2).
CROSS VALIDATION ANALYSIS OF THE GENERALISED LINEAR REGRESSION REDUCED MODEL
Figure 3 summarises the cross validation assessment of the “reduced model” for estimating the scintigraphic gastric emptying parameters (tLAG and t1/2). The standard deviation of the differences for tLAG was 11 minutes (range −23 to 23) and for t 1/2, 13 minutes (range −29 to 31). Figure 4 indicates a slightly broader distribution of residuals compared with those in fig 3C and D. However, the cross validated predicted values showed little change in the correlation coefficients.
STANDARD MEAL STUDY USING THE [13C]S PLATENSISBISCUIT
Twenty seven of the 30 original participants were recruited for a repeat gastric emptying study using [13C]S platensis biscuit, a prepacked standard meal. In fig 4, the distributions of tLAG and t1/2 estimates using the generalised linear regression model based on three breath13CO2 excretion data points are shown for the [13C]S platensis egg meal and the [13C]S platensis biscuit meal. There were small differences in t1/2 estimates within individuals (median three minutes, range −18 to 49) for the two meals and these differences are consistent with the day to day variation in normal gastric emptying (median 12–15%) previously reported.10 22 These differences are also summarised in table 1.
Our laboratory has further pursued the goal of developing a simple, non-radioactive office based, or field based, breath test using stable isotopes to measure gastric emptying of solids. Our data suggest that the methods proposed are ready for clinical application and for use in epidemiological studies. We showed that sufficiently accurate estimates of t1/2 and tLAG using the [13C]S platensis breath test could be obtained using the new method compared with scintigraphic methods. Mathematical analysis of breath test data appears to be important. In a subsequent application23 of the new approach, we showed that it can also be applied to measurement of gastric emptying of solids with [13C]octanoic acid as the substrate.
Our approach was to develop an analytical model which could overcome the essential difference between scintigraphy and breath tests, that is the unknown total gastric residual when a breath test is used. Following our observation that the early portion of the cumulative13CO2 breath curves was generally linear (and reciprocal to the scintigraphic emptying curve), we sought the time points during the three hour period of breath sampling which best reflected the scintigraphic estimates of gastric emptying in healthy subjects. Using this method we determined that three breath samples during the first three hours provided essentially equivalent information compared with estimates from all breath samples obtained during the three hour test. We have shown, using a generalised linear model approach, that the three hour data were essentially equivalent to the estimates from the six hour breath test, and that both gave good estimates of tLAG and t1/2 determined by scintigraphy. These time points (particularly at 90 and 180 minutes) most closely reflected the gastric half-emptying time. Interestingly, these time points are earlier than the time when one would anticipate errors from the redistribution of 13C from the body's bicarbonate pool which occurs typically more than three hours after a meal.
Using the generalised linear regression “reduced model”, breath samples were needed at only four times: baseline, 75, 90, and 180 minutes after the meal, with three data points each for tLAG (baseline, 75, and 180 minutes) and t1/2(baseline, 90, and 180 minutes). Moreover, the calculations were simplified relative to the nonlinear approach; predicted tLAG and t1/2 can be estimated using a pocket calculator. This simplified method can be applied in epidemiological or multicentre studies because of its safety, simplicity, accuracy, and low cost relative to scintigraphy.
As our model was developed based on data from 30 healthy volunteers, it is possible that our approach reflects only physiological gastric emptying. However, because the model incorporated a non-Gaussian data distribution, we anticipate it would be similarly accurate in disease states, particularly in delayed gastric emptying. We have validated the use of the model in diabetic patients using a [13C]octanoic acid meal as the source of13C.23 The relative accuracy of the test using13C spirulina substrates in disease states requires further study.
A final departure from previous studies in this analysis was the use of the actual quantity of 13C appearing in breath per minute. This was used to further simplify the calculation, obviating the need to collect all data points for the cumulative excretion plot. A delta over baseline value or a percentage of 13CO2recovery per hour can also be used, but we found that actual13C values reflected the most accurate overall results.
In summary, new calculations render the stable isotope [13C]S platensis breath test comparable with scintigraphy for estimating gastric emptying in healthy volunteers. The normal values for a standardised preprepared meal containing 13C were estimated. These approaches need to be applied in disease states and may ultimately enable the development of an office or field based non-radioactive test to measure gastric emptying of solids in humans.
This study was supported in part by a grant from Meretek Laboratories, Houston, Texas, by a grant to Meretek from the Janssen Research Foundation and by General Clinical Research Center grant (No RR00585) and RO1-DK54681–01 and K24-DK02638–01 (Dr M Camilleri) from the National Institutes of Health. The authors thank Mrs Cindy Stanislav for excellent secretarial assistance.
Potential conflict of interest: Dr P D Klein is an employee of Meretek Diagnostics, Inc. Two US patents have been awarded to Meretek Diagnostics, Inc. (Nos 5,707,602 and 5,785,949) covering the use of 13C spirulina (S platensis) in solid and liquid phase gastric emptying measurements.
- Abbreviations used in this paper:
- half-emptying time
- lag time
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.