- Research article
- Open Access
- Open Peer Review
Use of laboratory and administrative data to understand the potential impact of human parainfluenza virus 4 on cases of bronchiolitis, croup, and pneumonia in Alberta, Canada
BMC Infectious Diseases volume 16, Article number: 402 (2016)
Human Parainfluenza Virus (hPIV) causes severe respiratory illness in infants and adults. Our study describes the association of hPIV1–4 with bronchiolitis, croup, and pneumonia using retrospective laboratory, administrative and public health data. Due to issues including the historic lack of hPIV4 in some commercial respiratory virus panels, the description of the impact of hPIV4 on croup, bronchiolitis, and pneumonia at population levels has often been limited. This study will use routine clinical laboratory data, and administrative data to provide a preliminary description of the impact of hPIV4 on these diseases in our population.
A three year cohort of patients positive for hPIV was linked with data from physician visits and hospital admissions to define cases and hospitalization status. International Classification of Disease (ICD-9) codes were used to determine if cases had croup, bronchiolitis, and pneumonia. We also looked at differences in hospitalization status, age and gender among hPIV1–4. All statistical analysis was done using SPSS (Version 19.0.0, IBM Corp© 2010) and Graphpad Prism V6 (GraphPad Software, Inc., 2012).
Only hPIV1 and hPIV4 specimens had positivity rates greater than 5 % of all specimens sent for respiratory virus panel testing. hPIV1 exhibited a biennial pattern while the pattern for hPIV3 was less interpretable due to lower positivity rates. Circulation patterns for hPIV2 and hPIV4 were not assessed due to the low positivity rates of theses specimens. From 2010 to 2013, there were 2300 hPIV cases with hPIV3 (46 %) being the most common, followed by hPIV1 (27 %), hPIV4 (16 %) and hPIV2 (11 %). The median age was 2 years for all hPIV types. Males were slightly greater than females for hPIV1 and hPIV2, with an equal distribution for hPIV3 and slightly more females than males for hPIV4. hPIV1 and hPIV2 had the highest or proportion of croup while hPIV3 and hPIV4 had the highest proportion of pneumonia. Within hPIV4 cases, distributions of diseases were; pneumonia (21 %, 95 % CI 17.1–25.7), bronchiolitis (18 %, 95 % CI 14.3–22.5), croup (2 %, 95 % CI 0.8–3.9), mixed illness of any of pneumonia, bronchiolitis or croup (4 %, 95 % CI 2.5–7.0) or other respiratory diseases (54 %, 95 % CI 49.1–59.6).
We used laboratory and administrative data to undertake a descriptive analysis of the association of hPIV1–4 with croup, bronchiolitis and pneumonia. hPIV4 appears to be more associated more with bronchiolitis and pneumonia and less with croup in our population.
Respiratory illness due to viral infections represents a significant burden on the healthcare system and our jurisdiction is highly impacted by respiratory viruses from later summer until early spring . Human Parainfluenzavirus (hPIV) is a single stranded RNA virus belonging to Paramyxoviridae family, and all four types are a significant cause of respiratory illness, in infants and elderly . hPIVs cause upper respiratory tract illness (URTI) such as colds, otitis media, and pharyngitis or lower respiratory tract illness (LRTI) such as croup, bronchiolitis and pneumonia [3–5]. In adults, the infection tends to stay in the upper respiratory tract; however morbidity in children under the age of five can be severe, with hospitalizations due to pneumonia, croup, and bronchiolitis . The four types of parainfluenza differ in their clinical outcomes. hPIV1 and 2 have been associated with croup in children (acute laryngotracheobronchitis) , while hPIV3 has been described in bronchiolitis and pneumonia .
Although hPIV4 has been described in studies using laboratory developed tests, it has not often been included in some commercial respiratory viral panels and so large scale systematic surveillance may not have included hPIV4 due to regulatory issues . We note that there is already some work suggesting an association of hPIV4 infections with pneumonia  and some speculation that hPIV4 infections might resemble hPIV3 infections [10, 11]. However, studies on the impact of HPIV4 on specific diseases have been limited to locations that able to overcome this limitation and integrate testing into diagnostic algorithms and that studies may not have focused on whole populations. In contrast, our laboratory provides centralized testing for our Province (Population 4.1 Million) and was one of the early adopters of a diagnostic that could differentiate hPIV4 from hPHIV1-3 on a population level over a period of several years [10, 12, 13]. Due to a number of historical factors, we also undertake extremely high levels of respiratory virus testing a year with this technology (approximately 30,000 specimens a year), for a variety of patient populations which allows us extensive coverage of a large population from a variety of settings. Similarly, our publically funded health care system allows us to link billing codes indicating suspect diagnosis to patient laboratory information.
We also wanted to determine whether hPIV4 infections exhibited similar patterns of prevalence to other hPIV types. Variability of temporal and seasonal patterns of prevalence (annual, biannual, biennial) for hPIV types have been described in the literature but the data for hPIV4 has been relatively less described . hPIV4 infections in other populations have been described as occurring in a year round distribution, with biennial peaks in odd-numbered years and a spectrum of disease . In contrast prevalence patterns for the other hPIV types 1–3 may vary. Some studies in the United States have also shown that hPIV1 occurs in the fall of odd numbered years and hPIV2 occurs in even-numbered years . While hPIV3 is observed year round, it tends to peak primarily in the spring and early summer  or might have a biannual pattern of prevalence . HPIV2 has shown to have varied seasonal distribution, often occurring in multiple seasons and characterization of seasonality for HPIV2 is further confounded by the low prevalence of this virus [17, 18].
Therefore, the purpose of the study is to use routine clinical laboratory data, and administrative data collected by public health to provide a preliminary description of the impact of hPIV4 on croup, bronchiolitis and pneumonia in our population. We intended to do this by comparing the distribution of ICD9 codes associated with these diseases between different hPIV types. We also wanted to provide a better understanding of the circulation patterns of PIV4 over time in our province.
Alberta is the fourth largest province in Canada with a population of 4.1 million . Alberta has universal health care; as such, all inpatient and outpatient visits to physicians, emergency departments, and hospitals for all Alberta residents are billed to the Alberta Ministry of Health. All Alberta residents are issued a unique Public Health Number (PHN) at birth or upon moving to the province. Physicians are required to include this identifier and the diagnosis of the patient when billing for inpatient or outpatient services.
The Provincial Laboratory for Public Health (ProvLab) is a major diagnostic lab in the province of Alberta. All inpatient and outpatient respiratory samples throughout the province are labeled with the patient’s PHN and sent to ProvLab for laboratory testing and the resulting test results are exported and stored in a database. For the time period in this analysis, the determination to order respiratory virus testing was left to the discretion of the ordering physician and there were not recommendations provided on when a respiratory virus panel should be ordered. In combining the ProvLab and Ministry of Health data sources, it is possible to determine the organism an Alberta resident was infected with, whether they were seen as an inpatient or an outpatient, and whether they were hospitalized. The study utilized this capability to describe the epidemiology of hPIV1–4 in the population of Alberta.
During the study period, all respiratory specimens were routinely tested for influenza A and influenza B using a real-time reverse transcriptase-real time polymerase chain reaction (RT-PCR) assay as previously described ; if found negative, they were tested by the XTAG® Respiratory Viral Panel (RVP, Luminex, Austin, TX, USA) for several other common respiratory pathogens including, hPIV1–4 . Clinical respiratory samples from Alberta residents received at the ProvLab between November 1, 2010 and December 31, 2013 that tested positive for hPIV1–4 were eligible for the study. Specimen positivity rates were calculated by determining the number of positive hPIV1–4 specimens detected versus the number of specimens that were tested by the RVP panel.
As patients could have multiple samples sent for testing, only the first positive sample for an individual within a 365 day period for a specific type of hPIV was allowed in the study. Subsequent positive tests were excluded within the 365 day period, unless they were positive for a different type of hPIV, at which time they were considered a different case of disease. This period of time was chosen as we wanted to ensure that the values we presented were conservative estimates of the number of cases diagnosed during the study period. Cases that were co-infected with more than one type of hPIV or co-infected with another virus (such as rhino-enterovirus or respiratory syncytial virus [RSV].) were also excluded from the analysis to prevent the confounding effect of other viruses that may also be associated with croup, bronchiolitis and pneumonia .
Positive cases of hPIV were deterministically linked to administrative data using the PHN. The physician claims data in the Supplemental Enhanced Service Event (SESE) database was used to provide the World Health Organization (WHO) International Classification of Diseases 9 Codes (ICD 9)  associated with the physician visit and associated hIPV laboratory result. Three ICD 9 codes were extracted from the SESE database along with patient demographic characteristics. There were often multiple ICD9 codes, for the non-croup, non-bronchiolitis, and non-pneumonia case group. For this group frequencies of other disease codes included; “infectious diseases” (001–139), “ill defined condition” (780–799), acute respiratory infections (460–466); other diseases of the respiratory system (467–519), diseases of the ear and mastoid process (380–389), other diseases of the central nervous system and sense organs (320–379).
Physician claims data is housed in the Supplemental Enhanced Service Event System (SESE database); inpatient hospitalizations, emergency department visits, and outpatient clinic data is housed in the Morbidity and Ambulatory Care Abstracting Reporting (MACAR) system, which feeds into the Canadian Institute for Health Information’s (CIHI) Discharge Abstract Database (DAD) and National Ambulatory Care Reporting System (NACRS). . The hospitalization data provided information on whether the case was hospitalized and the duration of hospitalization.
These databases were utilized to determine demographic information, and to determine if a hPIV case was diagnosed with croup, pneumonia, or bronchiolitis within 14 days of their positive parainfluenza specimen. Table 1 lists the ICD-9 codes used to define each diagnosis. “Single diagnosis” was defined as hPIV cases that had one event of croup, bronchiolitis or pneumonia. “Mixed diagnoses” were defined as hPIV cases diagnosed with two or more events (i.e. croup and pneumonia).
Descriptive statistics regarding demographic characteristics and clinical diagnosis of hPIV1–4 cases were performed. We also described the distribution of hPIV cases over time. Proportions, means, medians were performed using SPSS (Version 19.0.0, IBM Corp© 2010). Graphpad Prism (V6, (GraphPad Software, Inc., 2012) was used to obtain 95 % confidence interval (CI) of proportions. Fisher’s exact and Chi- square tests were utilized to test for statistical significance.
Identification of positive hPIV cases was carried out as follows. A total of 55,112 specimens were tested for hPIV and included in the study. 63.5 % (51,888/81,587) specimens tested negative for hPIV, leaving 3224 positive specimens. Negatives were not included in Table 2 because the point of this analysis was to compare the distribution of ICD9 codes for croup, bronchiolitis, and pneumonia with specific hPIV types. Of those 3224 specimens, 316 were considered duplicate (i.e. based on case definition), leaving 2908 unique cases of hPIV. Twenty five cases were unable to be linked to the SESE and MACAR databases due to invalid PHN’s or other personal identifiers and were therefore excluded from study at this point. We were left with 2883 cases of hPIV that were linked to SESE and MACAR. Furthermore we excluded 19 cases that had mixed hPIVs, and 564 cases that were co-infected with other respiratory viruses leaving a total of 2300 unique hPIV cases over a three year period.
As in Table 2, hPIV4 was the third most common hPIV type making up 16 % of all hPIV cases. Out of all the cases, hPIV3 was the most common (46 %) type followed by hPIV1 (27 %). hPIV2 was the least common of all the hPIV cases (11 %). The median age was 2 years (range: 0.01–103 years) for all hPIV types. Males were slightly greater than females for hPIV1 and hPIV2, with an equal distribution for hPIV3 and slightly more females than males for hPIV4. Greater than 50 % of cases were non-hospitalized (Table 2). Over 50 % of cases were associated with ICD-9 codes related to other respiratory diseases that were not croup, bronchiolitis and pneumonia. As shown in Table 3, 95 % confidence intervals of the hPIV types for hospitalized cases all overlapped which indicated to use that were no able to comfortably infer a dominance of any one type in this population. The same inability to infer a dominant type of hPIV was also found with the non-hospitalized group. hPIV1 and hPIV2 had the highest percentages of croup while hPIV3 and hPIV4 had the highest percentages of pneumonia. Within hPIV4 cases, distributions of diseases were; pneumonia (21 %, 95 % CI 17.1–25.7), bronchiolitis (18 %, 95 % CI 14.3–22.5), croup (2 %, 95 % CI 0.8–3.9), mixed illness of any of pneumonia, bronchiolitis or croup (4 %, 95 % CI 2.5–7.0) or other respiratory diseases (54 %, 95 % CI 49.1–59.6).
The association of PIV1-4 was also assessed for other disease codes for the non-croup, non-bronchiolitis, and non-pneumonia cases (n = 1267). As there were multiple ICD9 codes, for each non-case in this group, frequencies of other disease codes were described. Many of these were relatively non-specific. Frequencies for “infectious diseases” (001–139) for each PIV type were; PIV1 [52/313, 16.61 %), PIV2 (30/132, 22.72 %), PIV3 (78/624, 12.50 %), and PIV4 (23/198. 11.62 %). Frequencies for “ill defined condition” (780–799) for PIV types were; PIV1 (137/313, 43.77 %), PIV2 (50/132, 37.88 %), PIV3 (252/624, 40.38 %), and PIV4 (87/198, 43.94 %). Frequencies for acute respiratory infections (460–466) for PIV types were; PIV1 (114/313, 36.42 %), PIV2 (38/132, 28.79 %), PIV3 (162/624, 25.96 %), and PIV4 (54/198, 27.27 %). Frequencies for other diseases of the respiratory system (467–519) for PIV types were; PIV1 (58/313, 18.53 %), PIV2 (26/132, 19.70 %), PIV3 (159/625, 25.44 %), an PIV4 (40/198, 20.20 %). Frequencies of, diseases of the ear and mastoid process (380–389) for each PIV type were; PIV1 (13/313, 4.15 %), PIV2 (4/132, 3.03 %), PIV3 (27/624, 4.32 %), and PIV4 (7/198, 3.53 %). Frequencies for other diseases of the central nervous system and sense organs (320–379) for each PIV type were; PIV1 (17/313, 5.43 %), PIV2 (5/132, 3.79 %), PIV3 (56/624, 8.97 %), and PIV4 (14/198, 7.07 %). No significant difference was seen between PIV type and diseases of the “ear and mastoid” (Chi-square, p = 0.89) nor for other diseases of the “nervous system and sense organs (Chi-square, p = 0.0724).
As in Fig. 1, only hPIV1 and hPIV3 specimens had peak positivity rates greater than 5 %. The peak patterns for hPIV1 were more defined and there appeared to be two winter peaks each separated by a one year period with no hPIV1 peak in what can be described as a biennial pattern (Fig. 1a). The peak patterns for hPIV3 (Fig. 1c) were less defined with one major peak in the winter period of 2010–2011, which occurred out of cycle with the hPIV1 peaks. Latter peaks for hPIV3 were less easy to interpret due to lower positivity rates that hovered around the five percent positivity level. Temporal peak patterns for hPIV2 (Fig. 1b) and hPIV4 (Fig. 1d) were not assessed as these did not peak above 5 % positivity.
This study uses administrative and hospitalization codes to link clinical laboratory data from hPIV1–4 cases to a clinical diagnosis (e.g. croup, bronchiolitis, pneumonia) and to determine patient setting (hospital or non-hospitalized). ICD-9 codes have been shown to be valid markers of a true history of pneumonia and can be used to successfully identify infection-related conditions in epidemiologic studies . Previous epidemiologic studies in Canada have used ICD-9 codes to understand patterns of croup. These codes have also been used to analyze bronchiolitis in epidemiologic studies in a variety of settings and over time [25–28]. ICD-9 codes have also been utilized to characterize croup cases on a population basis in other studies [29, 30]. The use of these data sets allows for analysis of general trends of illness at population levels in a manner that could not be accomplished with later more detailed chart reviews and surveys of electronic health records. We have also previously used hospitalization diagnosis codes in our jurisdiction to understand patient location trends on a population level . However, we hope that linking multiple databases will lead to additional studies that can provide a greater understanding of the role that that these hPIV types play in acute respiratory infections. Using administrative and clinical data we were able to study a largely pediatric patient population with nearly three quarters of all hPIV-positive patients age 12 or younger, and 42 % of hPIV cases as being hospitalized.
This work sheds light on some basic characteristics of hPIV4 infection at a larger population level. Of all the hPIV types, the clinical and epidemiologic presentation of hPIV4 on a large population level have been the least studied , although hPIV4 has been described as a cause of acute influenza-like-illness (ILI) in both adults and children in a variety of settings [13, 31–34]. This increasing body of work contrasts to earlier studies which identified hPIV4 as being associated with mild disease  or predominantly upper respiratory tract infections . Instead, our descriptive analysis aligns more with recent studies describing hPIV4 in acute lower respiratory tract infections in pediatric patients [36, 37], as an etiologic agent in hospitalized cases of community acquired pneumonia , or being associated with bronchiolitis. This difference in disease patterns may be due to the broader sampling from our study (e.g. the whole jurisdiction) including acute care and tertiary care facilities instead of community settings. Thus, we may be identifying a physician bias towards testing patients who are more ill and we maybe cases who have a greater severity of illness and are more likely to be evaluated.
One recent study by Frost et al., recently did an excellent job at describing the epidemiology of hPIV4 in hospitalized pediatric patients in Colorado, USA. That study did not identify any cases of hPIV4 associated with croup . In contrast, our study found that seven percent of hPIV4 cases were diagnosed with croup. This may be accounted for by the fact that our study was a large population-based study over an extended time period and large geographic area and included both pediatric and adult, and hospitalized and non-hospitalized cases.
Although some differences in type-associate disease profiles are seen in the study against previous work, the seasonality of hPIV4 was not determined in this study due to low peak positivity rates of hPIV4 specimens. Peak hPIV4 activity has been described in the late autumn and early winter . Depending on the study both hPIV4 and hPIV2, have either annual or biennial circulation patterns . In contrast a biennial circulation trend was observed for the hPIV1 cases and peak hPIV3 activity was offset from peak hPIV1 activity by one year [6, 15, 39–43]. Published and unpublished differences in temporal circulation patterns between studies is not unexpected, as multiple host and environmental factors impact on the seasonal distributions of these viruses .
There are several limitations to this study. First, the analysis does not describe the differences between the two subtypes of hPIV4 , and the work excludes co-infections of hPIV types or co-infections with other respiratory viruses. The study also does not attempt to understand the deeper clinical histories of patients, including factors such as immunocompromised status or transplant that might be associated with acute respiratory infections by specific hPIV types. The lack of a standardized case definition, reinfection rates, or maximum shedding periods for parainfluenza impact comparisons between studies [43, 46–48]. While 1 year may be high, we felt this was justified as we wanted to ensure that the values we presented were conservative estimates of the number of cases diagnosed during the study period, and that the estimates of hospitalization were not overestimated by the potential inclusion of very ill patients being tested multiple times for the same illness. However, we are reassured that only 316 of 3224 specimens were excluded from the study as duplicates. Finally, this is an analysis of large sets of administrative and laboratory data which may yield different conclusions than clinical studies using chart reviews. Although there is the possibility of errors in diagnostic codes, the authors feel that the diagnostic codes for emergency department and inpatient data are valid and reliable due to the rigorous training and review that the Canadian Institute for Health Information (CIHI) imposes [49, 50].
This study suggests that the distribution clinical illnesses may be different between hPIV4 and hPIV1. While the proportion of hPIV4 cases diagnosed with croup was statistically significantly lower than hPIV1 cases, hPIV1 cases were statistically significantly more likely to be diagnosed with pneumonia and bronchiolitis. While the results did not attain statistical significance, a similar trend was seen with hPIV2. The difference in the distribution of clinical diagnoses in hPIV types suggests that in the future it may be worthwhile to support hPIV type specific patient management strategies.
CIHI, Canadian Institute for Health Information; DAD, Discharge Abstract Database; hPIV, human parainfluenza virus; ICD-9, international classification of diseases-9; ILI, influenza-like illness; LRTI, lower respiratory tract infection; MACAR, morbidity and ambulatory care reporting database; NACRS, National Ambulatory Care Reporting System; PHN, public health number; RSV, respiratory syncytial virus; RT-PCR, real-time reverse transcriptase-real time polymerase chain reaction; SESE, Supplemental Enhanced Service Event System; URTI, upper respiratory tract infection; WHO, World Health Organization
Lee BE, Mukhi SN, May-Hadford J, Plitt S, Louie M, Drews SJ. Determination of the relative economic impact of different molecular-based laboratory algorithms for respiratory viral pathogen detection, including Pandemic (H1N1), using a secure web based platform. Virol J. 2011;8:277.
Counihan ME, Shay DK, Holman RC, Lowther SA, Anderson LJ. Human parainfluenza virus-associated hospitalizations among children less than five years of age in the United States. Pediatr Infect Dis J. 2001;20(7):646–53.
Stockmann C, Ampofo K, Hersh AL, Carleton ST, Korgenski K, Sheng X, Pavia AT, Byington CL. Seasonality of acute otitis media and the role of respiratory viral activity in children. Pediatr Infect Dis J. 2013;32(4):314–9.
Cantais A, Mory O, Pillet S, Verhoeven PO, Bonneau J, Patural H, Pozzetto B. Epidemiology and microbiological investigations of community-acquired pneumonia in children admitted at the emergency department of a university hospital. J Clin Virol. 2014;60(4):402–7.
McCuskee S, Kirlew M, Kelly L, Fewer S, Kovesi T. Bronchiolitis and pneumonia requiring hospitalization in young first nations children in Northern Ontario, Canada. Pediatr Infect Dis J. 2014;33(10):1023–6.
Marx A, Torok TJ, Holman RC, Clarke MJ, Anderson LJ. Pediatric hospitalizations for croup (laryngotracheobronchitis): biennial increases associated with human parainfluenza virus 1 epidemics. J Infect Dis. 1997;176(6):1423–7.
Schomacker H, Schaap-Nutt A, Collins PL, Schmidt AC. Pathogenesis of acute respiratory illness caused by human parainfluenza viruses. Curr Opin Virol. 2012;2(3):294–9.
McFarlane HJ, MacDonald J, Collins TC, Molyneaux PJ, Carman WF. Severe pneumonia after cardiac surgery as a result of infection with parainfluenza virus type 4. J Cardiothorac Vasc Anesth. 2009;23(1):84–6.
Frost HM, Robinson CC, Dominguez SR. Epidemiology and clinical presentation of parainfluenza type 4 in children: a 3-year comparative study to parainfluenza types 1–3. J Infect Dis. 2014;209(5):695–702.
Lau SK, Li KS, Chau KY, So LY, Lee RA, Lau YL, Chan KH, Lim WW, Woo PC, Yuen KY. Clinical and molecular epidemiology of human parainfluenza virus 4 infections in hong kong: subtype 4B as common as subtype 4A. J Clin Microbiol. 2009;47(5):1549–52.
Qu JX, Gu L, Pu ZH, Yu XM, Liu YM, Li R, Wang YM, Cao B, Wang C. Viral etiology of community-acquired pneumonia among adolescents and adults with mild or moderate severity and its relation to age and severity. BMC Infect Dis. 2015;15:89.
Lekana-Douki SE, Nkoghe D, Drosten C, Ngoungou EB, Drexler JF, Leroy EM. Viral etiology and seasonality of influenza-like illness in Gabon, March 2010 to June 2011. BMC Infect Dis. 2014;14:373.
Reed G, Jewett PH, Thompson J, Tollefson S, Wright PF. Epidemiology and clinical impact of parainfluenza virus infections in otherwise healthy infants and young children <5 years old. J Infect Dis. 1997;175(4):807–13.
Fry AM, Curns AT, Harbour K, Hutwagner L, Holman RC, Anderson LJ. Seasonal trends of human parainfluenza viral infections: United States, 1990–2004. Clin Infect Dis. 2006;43(8):1016–22.
Henrickson KJ. Parainfluenza viruses. Clin Microbiol Rev. 2003;16(2):242–64.
Ge X, Han Z, Chen H, Cheng J, Gao M, Sun H. Characterization of acute respiratory infections among 340 infants in Wuxi, Jiangsu Province. Ann Transl Med. 2015;3(18):264.
Gulen F, Yildiz B, Cicek C, Demir E, Tanac R. Ten year retrospective evaluation of the seasonal distribution of agent viruses in childhood respiratory tract infections. Turk Pediatr Ars. 2014;49(1):42–6.
http://www.finance.alberta.ca/aboutalberta/population-projections/index.html. Accessed 11 Jun 2015.
Drews SJ, Simmonds K, Usman HR, Yee K, Fathima S, Tipples G, Tellier R, Pabbaraju K, Wong S, Talbot J. Characterization of Enterovirus activity, including Enterovirus D68 in pediatric patients in Alberta, Canada: 2014. J Clin Microbiol. 2015;53(3):1042-5. doi: 10.1128/JCM.02982-14. Epub 2015 Jan 14.
Miller EK, Gebretsadik T, Carroll KN, Dupont WD, Mohamed YA, Morin LL, Heil L, Minton PA, Woodward K, Liu Z, Hartert TV, Williams JV. Viral etiologies of infant bronchiolitis, croup and upper respiratory illness during 4 consecutive years. Pediatr Infect Dis J. 2013;32(9):950–5.
http://www.cdc.gov/nchs/icd/icd9cm.htm. Accessed 2 Nov 2015. 31 Oct 2015.
https://www.cihi.ca/en/data-and-standards/data-holdings. Accessed 2 Nov 2015. 31 Oct 2015.
Drahos J, Vanwormer JJ, Greenlee RT, Landgren O, Koshiol J. Accuracy of ICD-9-CM codes in identifying infections of pneumonia and herpes simplex virus in administrative data. Ann Epidemiol. 2013;23(5):291–3.
Parikh K, Hall M, Mittal V, Montalbano A, Mussman GM, Morse RB, Hain P, Wilson KM, Shah SS. Establishing benchmarks for the hospitalized care of children with asthma, bronchiolitis, and pneumonia. Pediatrics. 2014;134(3):555–62.
Riese J, McCulloh RJ, Koehn KL, Alverson BK. Demographic factors associated with bronchiolitis readmission. Hosp Pediatr. 2014;4(3):147–52.
Sharifi M, Krishanswami S, McPheeters ML. A systematic review of validated methods to capture acute bronchospasm using administrative or claims data. Vaccine. 2013;31 Suppl 10:K12–20.
Mansbach JM, Emond JA, Camargo Jr CA. Bronchiolitis in US emergency departments 1992 to 2000: epidemiology and practice variation. Pediatr Emerg Care. 2005;21(4):242–7.
Roxburgh CS, Youngson GG, Townend JA, Turner SW. Trends in pneumonia and empyema in Scottish children in the past 25 years. Arch Dis Child. 2008;93(4):316–8.
Nath S, Thomas M, Spencer D, Turner S. Has the incidence of empyema in Scottish children continued to increase beyond 2005? Arch Dis Child. 2015;100(3):255-8.
Roy MT, Chanda S, Mullick S, De P, Dey-Sarkar M, Chawla-Sarkar M. Spectrum of respiratory viruses circulating in eastern India: prospective surveillance among patients with influenza-like illness during 2010–2011. J Med Virol. 2013;85(8):1459–65.
Abiko C, Mizuta K, Aoki Y, Ikeda T, Itagaki T, Noda M, Kimura H, Ahiko T. An outbreak of parainfluenza virus type 4 infections among children with acute respiratory infections during the 2011–2012 winter season in Yamagata, Japan. Jpn J Infect Dis. 2013;66(1):76–8.
Hasman H, Pachucki CT, Unal A, Nguyen D, Devlin T, Peeples ME, Kwilas SA. Aetiology of influenza-like illness in adults includes parainfluenzavirus type 4. J Med Microbiol. 2009;58(Pt 4):408–13.
Bialasiewicz S, McVernon J, Nolan T, Lambert SB, Zhao G, Wang D, Nissen MD, Sloots TP. Detection of a divergent Parainfluenza 4 virus in an adult patient with influenza like illness using next-generation sequencing. BMC Infect Dis. 2014;14:275.
Billaud G, Morfin F, Vabret A, Boucher A, Gillet Y, Crassard N, Galambrun C, Ferraris O, Legrand L, Aymard M, Lina B, Freymuth F, Thouvenot D. Human parainfluenza virus type 4 infections: a report of 20 cases from 1998 to 2002. J Clin Virol. 2005;34(1):48–51.
Lau SK, To WK, Tse PW, Chan AK, Woo PC, Tsoi HW, Leung AF, Li KS, Chan PK, Lim WW, Yung RW, Chan KH, Yuen KY. Human parainfluenza virus 4 outbreak and the role of diagnostic tests. J Clin Microbiol. 2005;43(9):4515–21.
Ren L, Gonzalez R, Xie Z, Xiong Z, Liu C, Xiang Z, Xiao Y, Li Y, Zhou H, Li J, Yang Q, Zhang J, Chen L, Wang W, Vernet G, Paranhos-Baccala G, Shen K, Wang J. Human parainfluenza virus type 4 infection in Chinese children with lower respiratory tract infections: a comparison study. J Clin Virol. 2011;51(3):209–12.
Caglayan SD, Pullukcu H, Cicek C, Sipahi OR, Tasbakan S, Atalay S. Bacterial and viral etiology in hospitalized community acquired pneumonia with molecular methods and clinical evaluation. J Infect Dev Ctries. 2014;8(4):510–8.
de Silva LM, Cloonan MJ. Brief report: parainfluenza virus type 3 infections: findings in Sydney and some observations on variations in seasonality world-wide. J Med Virol. 1991;35(1):19–21.
Rosychuk RJ, Klassen TP, Voaklander DC, Senthilselvan A, Rowe BH. Seasonality patterns in croup presentations to emergency departments in Alberta, Canada: a time series analysis. Pediatr Emerg Care. 2011;27(4):256–60.
Segal AO, Crighton EJ, Moineddin R, Mamdani M, Upshur RE. Croup hospitalizations in Ontario: a 14-year time-series analysis. Pediatrics. 2005;116(1):51–5.
Liu WK, Liu Q, Chen DH, Liang HX, Chen XK, Huang WB, Qin S, Yang ZF, Zhou R. Epidemiology and clinical presentation of the four human parainfluenza virus types. BMC Infect Dis. 2013;13:28.
Murphy B, Phelan PD, Jack I, Uren E. Seasonal pattern in childhood viral lower respiratory tract infections in Melbourne. Med J Aust. 1980;1(1):22–4.
Nascimento-Carvalho CM, Cardoso MR, Barral A, Araujo-Neto CA, Oliveira JR, Sobral LS, Saukkoriipi A, Paldanius M, Vainionpaa R, Leinonen M, Ruuskanen O. Seasonal patterns of viral and bacterial infections among children hospitalized with community-acquired pneumonia in a tropical region. Scand J Infect Dis. 2010;42(11–12):839–44.
Wang CY, Arden KE, Greer R, Sloots TP, Mackay IM. A novel duplex real-time PCR for HPIV-4 detects co-circulation of both viral subtypes among ill children during 2008. J Clin Virol. 2012;54(1):83–5.
Burke CW, Bridges O, Brown S, Rahija R, Russell CJ. Mode of parainfluenza virus transmission determines the dynamics of primary infection and protection from reinfection. PLoS Pathog. 2013;9(11):e1003786.
Welliver R, Wong DT, Choi TS, Ogra PL. Natural history of parainfluenza virus infection in childhood. J Pediatr. 1982;101(2):180–7.
Henderson FW. Pulmonary infections with respiratory syncytial virus and the parainfluenza viruses. Semin Respir Infect. 1987;2(2):112–21.
Canadian Institute for Health Information. Quality assurance processes applied to the discharge abstract and hospital mortality databases. Ottawa: Government of Canada; 2008.
Canadian Institute for Health Information. Canadian coding standards for ICD-10-CA and CCI for 2008. Ottawa: Government of Canada; 2007.
The Provincial Laboratory for Public Health, Alberta, Canada for testing all the specimens and providing the line level data.
The Ministry of Health, Alberta, Canada for extracting data from administrative database and merging laboratory data.
No funding was obtained for this study.
Availability of data and materials
Privacy legislation prohibits the sharing of individual level data to non-custodians of the data. Aggregate data can be shared upon request. Most data is presented in the analysis.
SF: Data acquisition, data analysis, case classification, writing. KS: Data analysis, Linking cases to public Health databases, statistical analysis. JI: Data analysis and statistical analysis. ANS: Data analysis, Linking cases to public Health databases, statistical analysis. SJD: Conceptualization, data analysis, writing. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was assessed by the Chair of the Human Research Ethics Board (HREB) at the University of Alberta. It was determined that the work described would not require review as per the Tri-Council Policy Statement 2 (TCPS2) Article 2.5.