Abstract

Acute gastrointestinal illness (AGI) is an important public health issue, with many pathogen sources and modes of transmission. A one-year telephone survey was conducted in Canada (2014-2015) to estimate the incidence of self-reported AGI in the previous 28 days and to describe health care seeking behaviour, using a symptom-based case definition. Excluding cases with respiratory symptoms, it is estimated that there are 0.57 self-reported AGI episodes per person-year, almost 19.5 million episodes in Canada each year. The proportion of cases seeking medical care was nearly 9%, of which 17% reported being requested to submit a sample for laboratory testing, and 49% of those requested complied and provided a sample. Results can be used to inform burden of illness and source attribution studies and indicate that AGI continues to be an important public health issue in Canada.

1. Introduction

Acute gastrointestinal illness (AGI) is an important public health issue, with substantial economic and human health impact [14]. Numerous countries in Europe, North and South America, Asia, and Australia and New Zealand have conducted population-based studies to estimate the incidence and burden of AGI in the community [521]. In Canada, the National Studies on Acute Gastrointestinal Illness (NSAGI) initiative was initiated in 1999 and developed population estimates in select regions of Canada [2224]. The NSAGI studies have been used to inform Canadian burden of enteric illness studies [2527].

In 2014-15, a population-based telephone survey, the Foodbook study, was conducted nationally to describe Canadians’ exposure to foods, animals, and water that may serve as sources for enteric illness pathogens and included questions related to AGI symptoms and care seeking behaviours [28]. The objective of this paper is to describe the prevalence, distribution, and symptoms of AGI and the health care seeking behaviours of individuals with AGI, across Canada for 2014-15 based on the Foodbook study.

2. Materials and Methods

2.1. Study Design and Data Collection

The Foodbook study was conducted in Canada’s 10 provinces and 3 territories over a one-year period (April 2014-April 2015) using a population-based telephone survey and included questions on AGI symptoms and related health care seeking behaviours. Households were randomly selected from a sampling frame of telephone numbers that consisted of land lines (70% listed and 10% random digit-dialing) and cell phones (20%). The survey was designed with a target total sample size of 11,016 surveys collected evenly over a 12-month period across four age groups (0–9, 10–19, 20–64, and 65+ years) and all 13 Canadian provinces and territories. In order to improve completion rates for younger age groups, when households contained children less than 18 years old, 50% of surveys were conducted with the child who would have the next birthday and 50% were conducted with the adult who would have the next birthday. If there were no children in the household, the survey was conducted with the adult who would have the next birthday. Complete details on participant selection and questionnaire administration can be found in the Foodbook Report [28]. The surveys were conducted by an independent research company contracted by the Public Health Agency of Canada (the Agency). Individuals were excluded if they could not speak the supported languages (English, French, Inuktitut, and on-demand verbal translation for other languages), if they did not have a listed land line or cellular telephone number or travelled outside their province or territory of residence during the seven days prior to interview. The Foodbook study was reviewed and approved by Health Canada and the Agency’s Research Ethics Board (REB 2013-0025) as well as the Newfoundland and Labrador Health Research Ethics Authority to meet a unique provincial legal requirement (HREB 13.238).

Weighted selection of survey participants to reflect the Canadian population was assigned using the following method. The forward sortation area (FSA, or first three digits of the postal code) collected for each respondent was converted to the most likely census metropolitan area (CMA). Using 2011 Census data, the CMA indicator along with age group, household type, province or territory, the number of people in the household, the number of land lines and cell phones in the household, and gender were used to calculate the individual-level survey expansion weight. To create the final weighting variable, a poststratification step used iterative proportional ranking with available control tables. The population reference year was 2011, representing a population of 33,400,000 [29].

The survey questions relating to AGI were developed to be consistent with the NSAGI studies previously conducted in Canada [2224]. Respondents were asked if they had experienced any vomiting or diarrhea in the 28 days prior to the interview. The module included questions about symptoms, their frequency and duration, existing chronic conditions and medication use, respiratory symptoms, and care seeking and stool submission behaviours. The survey was pilot-tested over a two-week period before full survey implementation.

2.2. Case Definitions

Illnesses reported to have started prior to the 28 days of the interview were excluded (3.6% weighted). Respondents who identified more than one episode of AGI during the 28 days prior to the interview were asked to respond only for their most recent episode. Respondents who did not report symptoms of AGI, as well as those identified as having self-reported that their diarrhea or vomiting was due to pregnancy, medical treatment (e.g., chemotherapy), or medical conditions (e.g., Crohn’s disease, colitis, irritable bowel syndrome, and alcoholism), were included in the noncase category.

Two case definitions of AGI were assessed: (a) a person reporting three or more loose stools in 24 hours or any vomiting in the past 28 days, according to the international AGI definition [24] and (b) a person reporting three or more loose stools in 24 hours or any vomiting in the past 28 days without concurrent respiratory symptoms (cough or sore throat) [25]. Removing cases with respiratory symptoms creates a more specific definition attempting to exclude respiratory infections that may cause gastrointestinal symptoms such as vomiting or diarrhea [30].

2.3. Missing Data and Resolution

Starting in November 2014 to the end of the survey time period (April 2015), due to an error in the Computer Assisted Telephone Interviewing (CATI) survey tool, 225 respondents who indicated AGI symptoms were not asked all the relevant survey questions. These participants responded “yes” to having any symptoms of vomiting or diarrhea but “no” to vomiting symptoms specifically; therefore, it is assumed that they only had diarrhea symptoms. Questions missed pertained to the subsection specific to diarrhea and included duration of diarrhea symptoms, number of stools, and related chronic conditions or medication use.

An adjustment for missing data was made by weighting completed interviews for respondents that reported symptoms of diarrhea only (April 2014 to October 2014; weighted n = 988,114.51) to account for the missing data (November 2014 to April 2015; weighted n = 855,500.56) by province (and respiratory symptoms for the more specific case definition), therefore increasing the weight assigned to the completed interviews (Appendix A). An assumption was made that there would be similar chronic disease and duration of symptoms between the diarrhea-only cases from April to October and November to April.

3. Analysis

Data analysis was performed in Stata 13.0 (StataCorp., Texas Station, TX) using the survey weight, and only weighted results are reported. Categorical variables were described using weighted percentages and the relative 95% confidence interval (CI). Individuals responding “don’t know” or “refused” to any question were excluded from the analysis of that question. Mean and median values were used to describe continuous variables.

For incidence rate calculations, respondents identifying multiple episodes were counted as a single episode. The primary outcome measure of monthly prevalence was defined as the number of respondents reporting AGI in the previous 28 days divided by the total number of respondents. Annual incidence rate and incidence proportion calculations were performed using formulas found in Appendix B [31].

The null hypothesis of no overall association between the prevalence of disease and province was tested using the Wald χ2 test, with a value cutoff of 0.05. The difference between the proportion of cases (i.e., the prevalence of AGI) in a specific province and the proportion of cases in all other provinces combined was tested using the χ2 test.

4. Results

The survey response rate was 19.9%, and a total of 10,798 residents responded to the survey (Table 1). There were 975 respondents (weighted n = 2,803,946) that indicated symptoms of diarrhea or vomiting in the past 28 days, reflecting a monthly prevalence of 8.5%. Of these respondents, 12% (weighted) reported that their diarrhea or vomiting in the past 4 weeks was caused by a medical condition, medication, or pregnancy and were counted in the noncase group. Using the international AGI case definition of three or more loose stools in 24 hours and any vomiting [24], the overall monthly prevalence was 5.7% (95% CI 4.6–7.2, weighted n = 1,887,588) corresponding to an annual incidence rate of 0.77 episodes/person-year (95% CI 0.61–0.97) and 26.3 million episodes of AGI per year in Canada (Table 2). After removal of cases with concurrent respiratory symptoms (25%), the monthly prevalence was 4.3% (95% CI 3.1–5.8, weighted n = 1,407,698), with an annual incidence rate of 0.57 episodes/person-year (95% CI 0.41–0.78). This reflects 19.4 million episodes of AGI per year in Canada.

Estimates of monthly prevalence and incidence of AGI nationally and by province for both AGI case definitions are presented in Table 2. There were some regional differences identified: the monthly prevalence for the province of Quebec was significantly lower compared to the rest of the provinces/territories combined () when assessing the international AGI case definition, and the monthly prevalence for the province of Saskatchewan was significantly lower compared to the rest of the provinces/territories combined () when assessing the more specific case definition of three or more loose stools in 24 hours and any vomiting without concurrent respiratory symptoms. A higher monthly prevalence was identified in Manitoba, Newfoundland, and the Territories; however, these differences were not statistically significant. The Territories and New Brunswick showed the greatest difference between prevalence of AGI when comparing AGI case definitions with and without concurrent respiratory symptoms (Figure 1).

When considering predisposing factors, 5.6% of AGI cases took prescription antibiotics in the previous 28 days and 5.0% of AGI cases with no respiratory symptoms in the previous 28 days. There was no clear seasonal pattern: lower monthly prevalence in February, June, and October and higher monthly prevalence in December and April (Figure 2).

When assessing the most specific AGI case definition, 61.1% of respondents reported experiencing diarrhea symptoms only, while 24.3% reported both vomiting and diarrhea and 14.6% reported vomiting only (Table 3). Of the cases who experienced diarrhea, 9.9% reported bloody diarrhea (95% CI 3.4–25.8). Duration of symptoms was longest for those who experienced both vomiting and diarrhea compared to those experiencing only one symptom (Table 4). Cases of the more specific AGI definition reported a mean of 4.43 episodes of diarrhea and 3.97 episodes of vomiting in a 24-hour period.

When using the more specific AGI case definition to assess care seeking behaviour, overall, 8.8% (95% CI 4.9–15.1) of cases visited a physician (Table 5). Of these, 17.1% (95% CI 7.5–34.5) were requested to submit a stool sample by a physician, and 49.0% (95% CI 17.6–81.2) of these submitted a stool sample. Hospitalizations were reported by 0.68% of cases with a mean hospital stay of 3.38 days (median 2). Of all cases, 0.11% reported taking antibiotics to treat their illness.

5. Discussion

This is the first nationwide survey conducted in Canada to describe the magnitude and distribution of AGI in the general population. Based on the more specific definition of AGI, excluding cases with respiratory symptoms, it is estimated that there are 0.57 (95% CI 0.41–0.78) self-reported AGI episodes per person-year or almost 19.5 million episodes of AGI in Canada each year. This estimate is lower than the rate of 0.63 (95% CI 0.57–0.69) episodes per person-year that was estimated based on the combined previous NSAGI studies and used in the Canadian estimates of foodborne illnesses [25]; however, the 95% confidence intervals of the current and previous estimates overlap indicating a lack of statistical difference.

This lower estimated incidence in the current study year (2014/2015) compared to the previous NSAGI rates (2002–2006) may be related to the different approach in survey design, specifically the use of a weighted sampling technique. Other differences include the exclusion of respondents who travelled outside their province or territory of residence during the seven days prior to the interview in the current study that may have been experiencing symptoms. As well, a true lower incidence may be explained by epidemiological trends including variability in norovirus trends from year to year, the impact of rotavirus vaccine on illness associated with rotavirus [32, 33], or possibly other public health interventions.

The annual rate of AGI when using the specific definition and excluding respiratory symptoms is comparable to estimates from the United States (US) (individual population studies 0.49, 0.54, and 0.73 and overall 0.60 episodes per person-year) [34] and lower than Italy (0.76) [14]. When comparing the international AGI case definition, the Canadian annual incidence rate (0.77 episodes per person-year) is lower than Germany [12], Denmark [35], Italy [14], Chile [6], Australia, and the US [24], which ranged from 0.83 to 1.4 episodes per person-year, but is higher than Ireland (0.64) and Malta (0.37) [24]. The proportion of respondents excluded due to chronic conditions, medication use, or pregnancy as the cause of their symptoms in the present study (12%) was lower than that in the previous NSAGI studies (16–19%) [22, 23, 36].

Comparison of provincial/territorial results for Ontario, Quebec, and Nunavut using the international AGI definition showed that the estimates were lower than previous provincial/territorial illness estimates using the same definition [3638]. The large variation in incidence between provinces/territories, though not statistically significant for many of them, does speak to the apparent regional differences of AGI incidence in Canada and the importance of capturing national information that reflects all provinces and territories. Furthermore, having provincial and territorial specific estimates enables individual jurisdictions to assess their AGI burden more specifically. This could be used to generate regional estimates of enteric illness and specific transmission routes (e.g., foodborne illness estimates for a specific province or territory) that could be used to inform public health activities (e.g., education and prevention campaigns).

The proportion of cases with respiratory symptoms (25%) is at the low end of the range reported by other countries reporting between 19% and 48% of cases experiencing concurrent respiratory symptoms [13, 14, 24, 35]. Using the more specific case definition creates a more conservative estimate, attempting to account for cases whose AGI symptoms may be caused by respiratory infections [14, 30, 34].

Diarrhea only was the predominant symptom profile of cases (55.9% and 61.1%, resp., for the two case definitions); this result falls between other studies, reporting a higher proportion in Germany (78%) and Denmark (64%) and lower proportions in Sweden, Italy, and Chile (30–40%) [6, 12, 14, 35, 39]. The proportion of diarrheal cases with bloody diarrhea (9.9%) was higher than that in other countries (3-4%) [5, 12, 13, 35]. This may be due to the small number of cases reporting bloody diarrhea and a large assigned weight due to study design.

The proportion of persons with AGI varied somewhat by season, which is similar to higher rates of AGI in winter months as what has been reported by previous Canadian studies [22, 40] as well as internationally in the US, Denmark, Italy, Sweden, and Germany [5, 12, 14, 35, 39]. This pattern is likely driven by viruses circulating in the winter months, particularly norovirus which is the most common cause of AGI in Canada [41]. The higher monthly prevalence observed in April may be an artefact due to the lower number of survey respondents in April (4% of the survey) compared to other months (approximately 8% each).

The percentage of cases who reported seeing a physician was low with only 8.8% seeking care, and 17.1% of these were requested to submit a stool sample. These values are weighted and are lower than those in the previous NSAGI studies (11–23% and 26–54%) [22, 23, 36]. The previous NSAGI studies were not age adjusted; therefore, the more frequent care seeking among the elderly may contribute to the higher overall results in the previous studies. The exclusion of individuals who travelled in the past seven days in the Foodbook study may also influence the lower results. Recent travel is associated with seeking medical care and having a sample requested [42, 43]; thus, these individuals who travelled may have been more likely to seek care and be requested to submit a stool sample. Lower care seeking rates would influence pathogen-specific estimates as it indicates greater underdiagnosis of cases. This should be considered among future burden of illness activities and how survey respondent weighting may influence this phenomenon.

Possible limitations of this study include the retrospective study design as it may be subject to recall bias. Retrospective studies in the UK (IID2) gave higher estimates of disease burden than prospective studies [19]. However, retrospective studies with longer recall periods gave lower estimated rates than studies with shorter recall periods [44]. Extrapolation from a reported seven-day prevalence was almost twice the rate of illness estimated when extrapolating from the month recall period [6, 35, 45].

The study response rate of 20% is lower than that in previous NSAGI studies [22, 40] and may be a source of bias if those who did not respond had different symptom profiles compared to those who participated in the study. Furthermore, misclassification of cases may have occurred due to excluding cases with chronic conditions or respiratory symptoms that might have been true infectious AGI cases.

The missing data for diarrheal cases from November to April due to the survey interview error were adjusted for based on known diarrheal cases captured from April to October; this however may not have accurately reflected the symptoms and behaviours of the missing cases and may have impacted the results. An assumption was made that there would be similar patterns (e.g., chronic disease and medical causes of symptoms, duration of symptoms, and care seeking behaviours) between the diarrhea-only cases from April to October and November to April. However, there may have been seasonal differences due to different pathogens circulating (e.g., norovirus in the winter or bacterial pathogens in the summer) or behavioural patterns (e.g., international winter travel or domestic summer recreational water exposure). From previous NSAGI studies, the monthly prevalence of AGI fluctuates seasonally with peaks seen in winter/early spring and again in summer [2224]. Symptom-specific monthly variations were observed in Ontario where diarrhea only was the predominant symptom for most months, followed by both vomiting and diarrhea combined; however, the statistical significance of these variations was not reported, and the general relationship between symptom profiles does not vary much across seasons [23]. The adjustment was made based on province/territory only as there were insufficient data to allow for age-gender-province/territory-based adjustments. Differences due to age and gender would be inherently incorporated into the province-based adjustment. However, age- and gender-specific results could not be described due to this adjustment approach.

Generating an estimate of the total amount of AGI in Canada provides the foundation for pathogen- and transmission route-specific burden of illness estimates. The lower incidence of AGI reported here will inform future activities to refine estimates of food and waterborne illness in Canada. Additionally, provincial- and territorial-specific estimates will enable individual jurisdictions to assess their AGI burden and generate region-specific public health plans that could include, for example, focused education campaigns, public health policies, or resource allocations toward prevention of AGI strategies.

Appendices

A. Adjustment for Missing Data

A.1. Missing Data

Starting in November 2014 to the end of the survey time period (April 2015), due to an error in the CATI survey tool, 225 respondents (representing 30.5% of the respondents with any symptoms in the past 28 days) who indicated AGI symptoms were not asked all the relevant survey questions. These participants responded “yes” to having any symptoms of vomiting or diarrhea but “no” to vomiting symptoms specifically; therefore, it is assumed that they only had diarrhea symptoms. Questions missed pertained to the subsection specific to diarrhea and included duration of diarrhea symptoms, number of stools, and related chronic conditions or medication use. Information on respiratory symptoms was captured for all cases.

A.2. Approach to Adjustment

The responses from diarrhea-only cases from April to October (i.e., respondents where all questions were correctly asked) were adjusted to account for the diarrhea-only cases from November to April that had missing variables, so they could not be assessed if they met the case definition. Complete responses were given additional weight to account for those that were incomplete. The weighted value of the respondents with incomplete interviews (n = 855,500.56) was assigned to the weighted responses for the respondents with complete interviews (n = 988,114.51).

To estimate the more specific case definition of AGI without respiratory symptoms, a different adjustment was made as information on the respiratory symptoms was available for all respondents. The weighted value of the respondents with AGI and no respiratory symptoms but with incomplete interviews (n = 661,095.32) was assigned to the weighted responses for the respondents with AGI and no respiratory symptoms and complete interviews (n = 702,464.68). As this more specific definition necessitated cases not have concurrent respiratory symptoms and that information was available from all respondents, noncases could be identified directly, and thus, their true weight was incorporated into the noncase definition without need for adjustment.

Using this approach, province/territory-specific multipliers were developed for cases of AGI (regardless of respiratory symptoms) and for cases of AGI without respiratory symptoms (Table 6). Essentially, the weight of completed interviews from each province/territory was given additional weight to account for those with missing responses in that province. This approach assumed that the distribution of chronic disease, medical condition, or medication use as the cause of AGI symptoms would be the same over time and that the number of stools for cases occurring in April to October would have the same distribution as that for cases occurring in November to April. Similarly, care seeking behaviours, duration of illness, hospitalization, etc. would also have the same distribution.

A.3. Adjustment 1 Example for Cases of AGI (Regardless of Respiratory Symptoms)

where A = weighted value of complete interviews for diarrhea only; B = weighted value of incomplete interviews for diarrhea only; C = adjustment multiplier (province specific) for complete interviews for diarrhea only; BC multiplier: (172,666.37 + 85,951.20)/172,666.37 = 1.50.

Therefore, the individual complete diarrhea-only responses for BC are weighted an additional 1.5 times.

A.4. Adjustment 2 Example for Cases of AGI without Respiratory Symptoms

where X = weighted value of complete interviews for diarrhea only and no respiratory symptoms; Y = weighted value of incomplete interviews for diarrhea only and no respiratory symptoms; Z = adjustment multiplier (province/territory specific) for complete interviews for diarrhea only and no respiratory symptoms; BC multiplier: (159,460.83 + 44,916.74)/159,460.83 = 1.28.

Therefore, the individual complete diarrhea-only and no respiratory symptom responses for BC are weighted an additional 1.28 times.

A.5. Adjustment for Seasonality Comparison

To report on seasonality by symptom, a different approach was needed to account for the missing data. The proportion of cases that experienced diarrhea only from April to October was calculated (60.04% and 46.72% for the international AGI case definition and the more specific case definition where cases with concurrent respiratory symptoms were removed, resp.). This proportion was then multiplied by the weighted total of respondents reporting diarrhea only that had missing information for each month to generate the estimate of the monthly number of cases that would have experienced diarrhea only. This value was combined with the reported values for the other symptom profiles (i.e., vomiting only and both diarrhea and vomiting) to estimate the seasonality of each symptom profile by month.where E = weighted value for month with incomplete interviews for diarrhea only; F = weighted value of completed interviews for diarrhea only that met case definition; G = weighted value of completed interviews for diarrhea only; H = estimated weighted value of respondents from month with incomplete data that met case definition.

A.6. Example

For January and the international AGI case definition, there were 47 respondents with a weighted value of 98,976, indicating that they experienced diarrhea symptoms, but with missing data, the weighted value of these respondents was multiplied by the proportion of completed weighted interviews that were diarrhea only (60.04% or 593,306.07/988,144.51) to estimate that 59,430 weighted respondents would have experienced diarrhea only in January. These values were then used to estimate the monthly prevalence by the symptom profile reported in Figure 2.

The difference of cases (98,967–59,430) would have been considered noncases that did not meet the case definition of having three or more loose stools in 24 hours or that their symptoms were related to a chronic condition or medication use.

A.7. Limitations

As the survey was a weighted study design, it required a weighted analysis, and thus, it was necessary to devise a weighted adjustment to address the CATI survey error. Ideally, we would have liked to perform this adjustment based on province/territory + age + gender weights so that we could comment on the differences among these demographics. However, due to data limitations in the province/territory + age + gender combinations where certain combinations were missing, this was not possible. Of the interviews with missing data from November 2014 to April 2015, about 8% (19/225) of AGI cases and 12% (18/144) of AGI cases without respiratory symptoms did not have completed interview data, affecting 11 and 14 of the province/territory + age + gender combinations, respectively. Therefore, province/territory + age + gender combinations were not used in the adjustment, and results on demographics could not be reported. The differences between age + gender combinations are inherently captured in the province-specific adjustment but are not able to be described explicitly. Available age and gender data were explored, and AGI results indicated little demographic difference to previous NSAGI studies (e.g., higher rates in children and lower rates in 65 years+ age group). Table 7 illustrates the demographic distribution of AGI cases with known responses and without adjustment for the missing 225 respondents for comparison with the adjusted values reported in the manuscript. This information however underestimates the true burden of AGI and thus cannot be used as the result for this research and is included only to demonstrate the general demographic conclusions that age and gender estimates do not differ greatly from what has been seen in previous NSAGI studies.

B. Prevalence and Annual Incidence Rate Calculation

Formulas for calculating prevalence and annual incidence rate are as follows [31]:

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors thank the Centre for Foodborne, Environmental and Zoonotic Infectious Diseases (Public Health Agency of Canada), Outbreak Management Division and Enteric Surveillance and Population Studies Division; the Foodbook team; Jennifer Cutler and Matt Hurst of the Foodbook team for their advice on adjustment for missing data; the staff at R.A. Malatest & Associates Ltd. for their expert interviewing; and the survey respondents for their participation. The authors thank the Public Health Agency of Canada for their funding support.