Abstract

Background

Human rhinoviruses (HRVs), human enteroviruses (HEVs) and human parechoviruses (HPeVs) have been linked to acute otitis media (AOM). We evaluated this association in a prospective birth cohort setting.

Methods

A total of 324 healthy infants were followed up from birth to age 3 years. Nasal swab samples were collected at age 3, 6, 12, 18, 24, and 36 months and screened for HRV and HEV using real-time reverse-transcription quantitative polymerase chain reaction. Stool samples were collected monthly and analyzed for HRV, HEV, and HPeV. AOM episodes diagnosed by physicians were reported by parents in a diary. The association of viruses with AOM was analyzed using generalized estimation equations, and their relative contributions using population-attributable risk percentages.

Results

A clear association was found between AOM episodes and simultaneous detection of HEV (adjusted odds ratio for the detection of virus in stools, 2.04; 95% confidence interval, 1.06–3.91) and HRV (1.54; 1.04–2.30). HPeV showed a similar, yet nonsignificant trend (adjusted odds ratio, 1.44; 95% confidence interval, .81–2.56). HRV and HEV showed higher population-attributable risk percentages (25% and 20%) than HPeV (11%).

Conclusions

HEVs and HRVs may contribute to the development of AOM in a relatively large proportion of cases.

Globally, upper respiratory infections (URIs) continue to be the leading cause of acute morbidity in all age groups [1]. URIs, most of which have a viral cause, are often complicated by acute otitis media (AOM), especially in young children [2, 3]. In studies carried out among children, approximately one-third of all URIs have been complicated by AOM, which usually occurs within the first 7 days after URI onset [2–5]. Thus, owing to their high frequency, viral URIs constitute a significant risk factor for AOM, making it one of the most common diseases in childhood. Furthermore, children encountering AOM early in life have been shown to be at increased risk of concurrent AOM later on [6]. These episodes result in widespread use of antibiotics and in chronic cases may lead to otologic surgery, causing significant financial losses to families and society [7].

AOM has generally been considered a bacterial complication of URI. However, viruses alone or together with bacterial copathogens have also been shown to cause AOM [8–10]. Furthermore, according to current consensus, the mucosal damage caused by virus infection is a prerequisite for bacteria to be able to enter the middle ear [11]. Thus, elimination of these viruses could offer an efficient means of preventing AOM. Of all viruses causing URIs, picornaviruses (family Picornaviridae) have been found to be among the most common, with human rhinoviruses (HRVs; genus Enterovirus) accounting for up to 73% of all respiratory infections recorded during the first year of life [12, 13]. HRVs as well as human enteroviruses (HEVs; genus Enterovirus) and human parechoviruses (HPeVs; genus Parechovirus) are picornaviruses that have been linked to AOM in a number of studies. In Finnish and Australian studies, HRVs and HEVs were the predominant viruses detected in children with AOM and recurrent AOM [9, 14, 15]. HPeVs have been found to play a smaller role in the development of AOM [16–18].

The aforementioned studies have mostly relied on methods based on reverse-transcriptase (RT) polymerase chain reaction (PCR) assays, the high sensitivity of which has led to the questioning of the clinical relevance of the results. These viruses have also been found in children without concurrent respiratory symptoms [19]. Hence, the presence of HRVs, HEVs, or HPeVs in respiratory or stool samples alone does not establish causality of the concurrent illness. In this study we investigated the possible causal relationship by comparing the frequency of HRV, HEV, and HPeV infections during AOM episodes and AOM-free periods in a cohort of children who were followed up from birth to age 3 years. We also calculated the population-attributable risk (PAR) of these picornaviruses to the development of AOM.

METHODS

Children

This study is a secondary analysis of data collected during the prospective international DIABIMMUNE study carried out from September 2008 to October 2013. The birth cohort included in this study comprises 324 healthy infants (169 boys and 155 girls) born between September 2008 and May 2010, who were enrolled in the DIABIMMUNE study at Jorvi Hospital, Espoo, Finland. All newborns, whose families gave written informed consent to the study, were first screened for HLA DR/DQ alleles conferring susceptibility to type 1 diabetes, using cord blood. Those with an eligible HLA genotype were recruited to prospective follow-up; these constituted approximately 5% of all newborn infants [20]. The local ethics committee approved the DIABIMMUNE study protocol.

Recording of Otitis Media Episodes and Background Characteristics

During the study period, children visited the study clinics at age 3, 6, 12, 18, 24, and 36 months, making 1836 visits altogether to the study clinics. At the first visit, the parents were given a diary in which they were advised to record information related to the child’s diet, allergies, environmental exposures (eg, type of daycare, pets, living environment, and parental smoking), infections, use of medications, and vaccinations between the visits. The number of siblings was also documented and updated during the follow-up period. All data from the diaries were transferred to a central database during each visit. The number of otitis media (OM) episodes that the child had encountered since the last clinic visit was recorded in the database at each visit. The date, the codes from the International Statistical Classification of Diseases, Tenth Revision (ICD-10), and the type of OM diagnosis (parental report or diagnosed by a physician) were also recorded.

After data extraction, OM diagnoses were classified as acute or chronic based on the ICD-10 codes recorded. As customary in Finland for cases diagnosed by means of otoscopy and without myringotomy, >99% of the episodes were recorded as H66.9 (unspecified OM). These were assumed to be acute. All OM episodes reported by the parents were counted as separate episodes regardless of the interval between 2 consecutive episodes.

Detection of Viruses

To assess HRV, HEV, and HPeV infections, a study nurse collected nasal swab (NS) samples on each visit to the study clinic, using a pediatric flocked swab (Copan Diagnostics) After sampling, the swab was put into a tube containing UTM-RT medium (Copan Diagnostics). The samples were frozen immediately and stored at -80°C until analyzed. Parents were also asked to collect stool samples from the children at home every month starting at age 1 month. Parents stored the samples at home at -20°C and brought them to the study clinic at scheduled visits. The frozen samples were then stored at -80°C until analyzed.

Altogether, 1387 NS samples were collected at age 3, 6, 12, 18, 24, and 36 months (mean, 4.3 samples per child), and 2861 stool samples were collected monthly until age 1 year (mean, 8.8 samples per child). At least 1 NS sample was obtained and analyzed from 279 children, and at least 1 stool sample from 316 children.

The NS and stool samples were screened for the presence of viral RNA with real-time RT quantitative PCR (qPCR). First, a 10% stool suspension was prepared from each original stool sample using Hanks solution. Next, RNA was extracted from 140 µL of suspension using the modified Qiagen RNeasy96 kit (Qiagen) to obtain a final 60-uL volume of extracted RNA. Similarly, 140 uL of each NS sample was subjected to RNA extraction to obtain a final 60-uL volume of extracted RNA. RT-qPCR detection of HRVs [21], HEVs [22], and HPeVs [23] was performed using the QuantiTect Probe kit (Qiagen), as described in the respective references. All samples were analyzed in triplicate paralleling RT-qPCR runs, and a sample was considered positive if ≥1 of the triplicate runs yielded a positive result. The cycle threshold (CT) values for positivity were 44 for HRV and 42 for HEV and HPeV. Among the virus-positive NS samples, 77.9% and 73.3% of were positive for HRV and HEV in all 3 runs, respectively. Of the virus-positive stool samples, 63.3%, 68.8%, and 69.7% yielded a positive result in all 3 runs for HRV, HEV, and HPeV, respectively.

Only one-fourth of the NS samples were screened for HPeV. These analyses were discontinued owing to the very low number of positive results (0.5%). Therefore, the HPeV findings for NS samples were excluded from statistical analyses. HRV and HEV PCR was repeatedly validated in regular external quality control rounds (Quality Control for Molecular Diagnostics), in which they got the best possible scores.

Statistical Methods

The background characteristics of children with or without AOM were compared by means of Mann-Whitney U test. The association between host and environmental factors and HRV, HEV, and HPeV infections was analyzed using the Fisher exact test. The association between virus detection and AOM episodes was assessed first by means of univariate logistic regression. The covariates used in the final multivariable models were identified using Akaike information criterion in a stepwise algorithm. Generalized linear models with binomial family and logit link-function were used to fit statistical models. Generalized estimation equations with binomial family and logit link-function were used to take into account intraindividual correlation. Adjusted odds ratios (ORs) and their 95% confidence intervals were presented as the measure of the associations analyzed.

Only the AOM episodes diagnosed by a physician were included in the statistical analyses. Differences were considered statistically significant at P < .05. The analyses were performed using R software, version 3.6.1 (2019-07-05; http://r-project.org), with the stats and gee packages. The contributions of HRV, HEV, and HPeV to AOM were quantified with PAR percentages, which were calculated using the following formula: PAR = P*(OR - 1)/[P*(OR - 1) + 1], where P represents the prevalence of infection in the population (here the proportion of children who experienced ≥1 infection by age 1 year).

RESULTS

The total duration of follow-up was 930 child-years, and the mean duration 34.4 months per child (standard deviation, 6.1 months). Altogether, 930 AOM episodes were recorded (median, 2; range, 0–18 per child), showing a peak at age 19–24 months (14 AOM episodes per 100 child-months) (Figure 1). Being firstborn offered protection from AOM, whereas having a pet living mostly indoors was significantly associated with a higher number of AOM episodes. Table 1 shows the associations between different background characteristics and AOM.

Table 1.

Background Characteristics of Children and Their Impact on the Mean Number of Acute Otitis Media Episodes During Follow-up

CharacteristicChildren, No. (%)AOM Episodes, Mean No.P Valuea
Sex
 Male169 (52.2)3.0.45
 Female155 (47.8)2.7
Firstborn child
 Yes170 (52.5)2.6.045
 No154 (47.5)3.2
No. of siblings
 ≥1200 (61.7)3.0.12
 0124 (38.2)2.6
Exclusive breastfeeding at age 3 mo
 Yes178 (66.9)3.1.19
 No88 (33.1)2.6
Duration of breastfeeding ≥6 mo
 Yes226 (73.6)3.1.053
 No81 (26.4)2.6
Use of pacifier
 Yes200 (73.3)3.1.80
 No73 (26.7)2.9
Parental smoking
 Yes61 (19.4)2.9.57
 No253 (80.6)2.9
Pet living mostly indoors
 Yes121 (38.1)3.3.04
 No197 (61.9)2.6
Age at starting in daycare
 <12 mo
  Yes37 (15.7)3.1.31
  No198 (84.3)3.4
 <18 mo
  Yes128 (54.5)3.6.17
  No107 (45.5)2.9
 <24 mo
  Yes181 (77.0)3.5.22
  No54 (23.0)2.7
Living environment
 Urban260 (82.8)2.9.74
 Rural35 (11.1)3.0
 Both19 (6.1)3.1
CharacteristicChildren, No. (%)AOM Episodes, Mean No.P Valuea
Sex
 Male169 (52.2)3.0.45
 Female155 (47.8)2.7
Firstborn child
 Yes170 (52.5)2.6.045
 No154 (47.5)3.2
No. of siblings
 ≥1200 (61.7)3.0.12
 0124 (38.2)2.6
Exclusive breastfeeding at age 3 mo
 Yes178 (66.9)3.1.19
 No88 (33.1)2.6
Duration of breastfeeding ≥6 mo
 Yes226 (73.6)3.1.053
 No81 (26.4)2.6
Use of pacifier
 Yes200 (73.3)3.1.80
 No73 (26.7)2.9
Parental smoking
 Yes61 (19.4)2.9.57
 No253 (80.6)2.9
Pet living mostly indoors
 Yes121 (38.1)3.3.04
 No197 (61.9)2.6
Age at starting in daycare
 <12 mo
  Yes37 (15.7)3.1.31
  No198 (84.3)3.4
 <18 mo
  Yes128 (54.5)3.6.17
  No107 (45.5)2.9
 <24 mo
  Yes181 (77.0)3.5.22
  No54 (23.0)2.7
Living environment
 Urban260 (82.8)2.9.74
 Rural35 (11.1)3.0
 Both19 (6.1)3.1

Abbreviation: AOM, acute otitis media.

aP values based on Wilcoxon rank sum test.

Table 1.

Background Characteristics of Children and Their Impact on the Mean Number of Acute Otitis Media Episodes During Follow-up

CharacteristicChildren, No. (%)AOM Episodes, Mean No.P Valuea
Sex
 Male169 (52.2)3.0.45
 Female155 (47.8)2.7
Firstborn child
 Yes170 (52.5)2.6.045
 No154 (47.5)3.2
No. of siblings
 ≥1200 (61.7)3.0.12
 0124 (38.2)2.6
Exclusive breastfeeding at age 3 mo
 Yes178 (66.9)3.1.19
 No88 (33.1)2.6
Duration of breastfeeding ≥6 mo
 Yes226 (73.6)3.1.053
 No81 (26.4)2.6
Use of pacifier
 Yes200 (73.3)3.1.80
 No73 (26.7)2.9
Parental smoking
 Yes61 (19.4)2.9.57
 No253 (80.6)2.9
Pet living mostly indoors
 Yes121 (38.1)3.3.04
 No197 (61.9)2.6
Age at starting in daycare
 <12 mo
  Yes37 (15.7)3.1.31
  No198 (84.3)3.4
 <18 mo
  Yes128 (54.5)3.6.17
  No107 (45.5)2.9
 <24 mo
  Yes181 (77.0)3.5.22
  No54 (23.0)2.7
Living environment
 Urban260 (82.8)2.9.74
 Rural35 (11.1)3.0
 Both19 (6.1)3.1
CharacteristicChildren, No. (%)AOM Episodes, Mean No.P Valuea
Sex
 Male169 (52.2)3.0.45
 Female155 (47.8)2.7
Firstborn child
 Yes170 (52.5)2.6.045
 No154 (47.5)3.2
No. of siblings
 ≥1200 (61.7)3.0.12
 0124 (38.2)2.6
Exclusive breastfeeding at age 3 mo
 Yes178 (66.9)3.1.19
 No88 (33.1)2.6
Duration of breastfeeding ≥6 mo
 Yes226 (73.6)3.1.053
 No81 (26.4)2.6
Use of pacifier
 Yes200 (73.3)3.1.80
 No73 (26.7)2.9
Parental smoking
 Yes61 (19.4)2.9.57
 No253 (80.6)2.9
Pet living mostly indoors
 Yes121 (38.1)3.3.04
 No197 (61.9)2.6
Age at starting in daycare
 <12 mo
  Yes37 (15.7)3.1.31
  No198 (84.3)3.4
 <18 mo
  Yes128 (54.5)3.6.17
  No107 (45.5)2.9
 <24 mo
  Yes181 (77.0)3.5.22
  No54 (23.0)2.7
Living environment
 Urban260 (82.8)2.9.74
 Rural35 (11.1)3.0
 Both19 (6.1)3.1

Abbreviation: AOM, acute otitis media.

aP values based on Wilcoxon rank sum test.

Proportion of virus-positive nasal swab (NS) and stool samples and numbers of acute otitis media (AOM) episodes per month during the indicated follow-up periods (by age). Abbreviations: HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus.
Figure 1.

Proportion of virus-positive nasal swab (NS) and stool samples and numbers of acute otitis media (AOM) episodes per month during the indicated follow-up periods (by age). Abbreviations: HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus.

Three hundred two NS samples (21.8%) and stool samples 594 (20.8%) were positive for HRV, and 25 (1.8%) and 83 (2.9%), respectively, for HEV. Furthermore, 145 stool samples (5.1%) tested positive for HPeV. The proportion of virus-positive NS samples was highest at age 19–24 months for both HRV and HEV. In stool samples obtained during the first year of life, the frequencies of HEV and HPeV were highest at age 7–12 months, whereas HRV infections were most common at age 4–6 months (Figure 1). A clear seasonal pattern of HRV, HEV, and HPeV infections was observed, with HRV infections peaking in late spring/early summer and in autumn, and HEV and HPeV infections peaking in late summer and autumn (Figure 2).

Proportions of acute otitis media (AOM) episodes and human rhinovirus (HRV), human enterovirus (HEV), and human parechovirus (HPeV) RNA findings in nasal swab and/or stool samples detected each month, of all AOM episodes and respective virus findings.
Figure 2.

Proportions of acute otitis media (AOM) episodes and human rhinovirus (HRV), human enterovirus (HEV), and human parechovirus (HPeV) RNA findings in nasal swab and/or stool samples detected each month, of all AOM episodes and respective virus findings.

Having siblings significantly increased the proportion of HRV-positive samples (Table 2). Conversely, having siblings and living in the countryside significantly reduced the proportion of HPeV-positive samples. The median starting age for daycare was 17 months (interquartile range, 13–23 months). Starting daycare before age 12 months significantly increased the number of HEV-positive samples.

Table 2.

Association Between Background Characteristics and Proportions of Virus-Positive Samples During Follow-upa

HRVHEVHPeV
CharacteristicSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P Valueb
Sex
 Male2210485 (21.9).16221051 (2.3).33153283 (5.4).39
 Female2038411 (20.2)203857 (2.8)132062 (4.7)
Firstborn child
 Yes2305374 (16.2)<.001230550 (2.2).10152695 (6.2).004
 No1943522 (26.9)194358 (3.0)132650 (3.8)
Siblings
 >12604628 (24.1)<.001260469 (2.6).618176872 (4.1).002
 None1644268 (16.3)164439 (2.4)108473 (6.7)
Exclusive breastfeeding at age 3 mo
 Yes2475539 (21.8).86247563 (2.5).49168490 (5.3).77
 No1040223 (21.4)104031 (3.0)69339 (5.6)
Duration of breastfeeding >6 mo
 Yes3077657 (21.4).55307779 (2.6).732074105 (5.1).54
 No980200 (20.4)98027 (2.8)65337 (5.7)
Use of pacifier
 Yes2690558 (20.7).85269068 (2.5)>.99179197 (5.4)>.99
 No976205 (21.0)97624 (2.5)66135 (5.3)
Parental smoking
 Yes822160 (19.5).2382223 (2.8).7154633 (6.0).33
 No3324713 (21.5)332485 (2.6)2236111 (5.0)
Pet living mostly indoors
 Yes1583322 (20.3).46158344 (2.8).55105053 (5.0)>.99
 No2619559 (21.3)261964 (2.4)177290 (5.1)
Age at starting in daycare
 <12 mo
  Yes664143 (21.5).5266430 (4.5)<.00144928 (6.2).49
  No2543518 (20.4)254350 (2.0)169591 (5.4)
 <18 mo
  Yes1869401 (21.5).17186955 (2.9).07126370 (5.5)>.99
  No1338260 (19.4)133825 (1.9)88149 (5.6)
 <24 mo
  Yes2488503 (20.2).32248867 (2.7).22167487 (5.2).17
  No719158 (22.0)71913 (1.8)47032 (6.8)
Living environment
 Urban3407705 (20.7).15340788 (2.6).142278124 (5.4).03
 Rural or both731169 (23.1)73112 (1.6)50816 (3.1)
HRVHEVHPeV
CharacteristicSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P Valueb
Sex
 Male2210485 (21.9).16221051 (2.3).33153283 (5.4).39
 Female2038411 (20.2)203857 (2.8)132062 (4.7)
Firstborn child
 Yes2305374 (16.2)<.001230550 (2.2).10152695 (6.2).004
 No1943522 (26.9)194358 (3.0)132650 (3.8)
Siblings
 >12604628 (24.1)<.001260469 (2.6).618176872 (4.1).002
 None1644268 (16.3)164439 (2.4)108473 (6.7)
Exclusive breastfeeding at age 3 mo
 Yes2475539 (21.8).86247563 (2.5).49168490 (5.3).77
 No1040223 (21.4)104031 (3.0)69339 (5.6)
Duration of breastfeeding >6 mo
 Yes3077657 (21.4).55307779 (2.6).732074105 (5.1).54
 No980200 (20.4)98027 (2.8)65337 (5.7)
Use of pacifier
 Yes2690558 (20.7).85269068 (2.5)>.99179197 (5.4)>.99
 No976205 (21.0)97624 (2.5)66135 (5.3)
Parental smoking
 Yes822160 (19.5).2382223 (2.8).7154633 (6.0).33
 No3324713 (21.5)332485 (2.6)2236111 (5.0)
Pet living mostly indoors
 Yes1583322 (20.3).46158344 (2.8).55105053 (5.0)>.99
 No2619559 (21.3)261964 (2.4)177290 (5.1)
Age at starting in daycare
 <12 mo
  Yes664143 (21.5).5266430 (4.5)<.00144928 (6.2).49
  No2543518 (20.4)254350 (2.0)169591 (5.4)
 <18 mo
  Yes1869401 (21.5).17186955 (2.9).07126370 (5.5)>.99
  No1338260 (19.4)133825 (1.9)88149 (5.6)
 <24 mo
  Yes2488503 (20.2).32248867 (2.7).22167487 (5.2).17
  No719158 (22.0)71913 (1.8)47032 (6.8)
Living environment
 Urban3407705 (20.7).15340788 (2.6).142278124 (5.4).03
 Rural or both731169 (23.1)73112 (1.6)50816 (3.1)

Abbreviations: HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus.

aFor HRV and HEV, results are shown for stool and nasal swab samples combined.

bP values based on Fisher exact test.

Table 2.

Association Between Background Characteristics and Proportions of Virus-Positive Samples During Follow-upa

HRVHEVHPeV
CharacteristicSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P Valueb
Sex
 Male2210485 (21.9).16221051 (2.3).33153283 (5.4).39
 Female2038411 (20.2)203857 (2.8)132062 (4.7)
Firstborn child
 Yes2305374 (16.2)<.001230550 (2.2).10152695 (6.2).004
 No1943522 (26.9)194358 (3.0)132650 (3.8)
Siblings
 >12604628 (24.1)<.001260469 (2.6).618176872 (4.1).002
 None1644268 (16.3)164439 (2.4)108473 (6.7)
Exclusive breastfeeding at age 3 mo
 Yes2475539 (21.8).86247563 (2.5).49168490 (5.3).77
 No1040223 (21.4)104031 (3.0)69339 (5.6)
Duration of breastfeeding >6 mo
 Yes3077657 (21.4).55307779 (2.6).732074105 (5.1).54
 No980200 (20.4)98027 (2.8)65337 (5.7)
Use of pacifier
 Yes2690558 (20.7).85269068 (2.5)>.99179197 (5.4)>.99
 No976205 (21.0)97624 (2.5)66135 (5.3)
Parental smoking
 Yes822160 (19.5).2382223 (2.8).7154633 (6.0).33
 No3324713 (21.5)332485 (2.6)2236111 (5.0)
Pet living mostly indoors
 Yes1583322 (20.3).46158344 (2.8).55105053 (5.0)>.99
 No2619559 (21.3)261964 (2.4)177290 (5.1)
Age at starting in daycare
 <12 mo
  Yes664143 (21.5).5266430 (4.5)<.00144928 (6.2).49
  No2543518 (20.4)254350 (2.0)169591 (5.4)
 <18 mo
  Yes1869401 (21.5).17186955 (2.9).07126370 (5.5)>.99
  No1338260 (19.4)133825 (1.9)88149 (5.6)
 <24 mo
  Yes2488503 (20.2).32248867 (2.7).22167487 (5.2).17
  No719158 (22.0)71913 (1.8)47032 (6.8)
Living environment
 Urban3407705 (20.7).15340788 (2.6).142278124 (5.4).03
 Rural or both731169 (23.1)73112 (1.6)50816 (3.1)
HRVHEVHPeV
CharacteristicSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P ValuebSamples, No.Positive Samples, No. (%)P Valueb
Sex
 Male2210485 (21.9).16221051 (2.3).33153283 (5.4).39
 Female2038411 (20.2)203857 (2.8)132062 (4.7)
Firstborn child
 Yes2305374 (16.2)<.001230550 (2.2).10152695 (6.2).004
 No1943522 (26.9)194358 (3.0)132650 (3.8)
Siblings
 >12604628 (24.1)<.001260469 (2.6).618176872 (4.1).002
 None1644268 (16.3)164439 (2.4)108473 (6.7)
Exclusive breastfeeding at age 3 mo
 Yes2475539 (21.8).86247563 (2.5).49168490 (5.3).77
 No1040223 (21.4)104031 (3.0)69339 (5.6)
Duration of breastfeeding >6 mo
 Yes3077657 (21.4).55307779 (2.6).732074105 (5.1).54
 No980200 (20.4)98027 (2.8)65337 (5.7)
Use of pacifier
 Yes2690558 (20.7).85269068 (2.5)>.99179197 (5.4)>.99
 No976205 (21.0)97624 (2.5)66135 (5.3)
Parental smoking
 Yes822160 (19.5).2382223 (2.8).7154633 (6.0).33
 No3324713 (21.5)332485 (2.6)2236111 (5.0)
Pet living mostly indoors
 Yes1583322 (20.3).46158344 (2.8).55105053 (5.0)>.99
 No2619559 (21.3)261964 (2.4)177290 (5.1)
Age at starting in daycare
 <12 mo
  Yes664143 (21.5).5266430 (4.5)<.00144928 (6.2).49
  No2543518 (20.4)254350 (2.0)169591 (5.4)
 <18 mo
  Yes1869401 (21.5).17186955 (2.9).07126370 (5.5)>.99
  No1338260 (19.4)133825 (1.9)88149 (5.6)
 <24 mo
  Yes2488503 (20.2).32248867 (2.7).22167487 (5.2).17
  No719158 (22.0)71913 (1.8)47032 (6.8)
Living environment
 Urban3407705 (20.7).15340788 (2.6).142278124 (5.4).03
 Rural or both731169 (23.1)73112 (1.6)50816 (3.1)

Abbreviations: HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus.

aFor HRV and HEV, results are shown for stool and nasal swab samples combined.

bP values based on Fisher exact test.

A majority of the participants (82.4%) experienced ≥1 HRV infection, whereas 18.2% and 29.2% had ≥1 HEV or HPeV infection during the first year of life, respectively. At least 1 AOM episode was registered in 40.7% of participants. Among infants who had ≥1 AOM episode diagnosed by age 1 year, the mean numbers of HRV- and HEV-positive stool samples were higher than among those who did not encounter any AOM episodes by that age (Table 3). Such an association was not seen in NS samples.

Table 3.

Mean Number of Virus-Positive Nasal Swab and Stool Samples Among Children With or Without ≥1 Acute Otitis Media Episode by Age 1 Year

Virus-Positive Samples, Mean No.P Valuea
Sample Type and Virus≥1 AOM EpisodeNo AOM Episodes
Nasal swabn = 95bn = 144b
 HRV0.50.5.50
 HEV<0.1<0.1.70
Stooln = 129n = 186
 HRV2.31.6<.001
 HEV0.40.2.04
 HPeV0.40.5.93
Nasal swab and/or stooln = 132n = 192
 HRV2.61.9.002
 HEV0.40.2.04
Virus-Positive Samples, Mean No.P Valuea
Sample Type and Virus≥1 AOM EpisodeNo AOM Episodes
Nasal swabn = 95bn = 144b
 HRV0.50.5.50
 HEV<0.1<0.1.70
Stooln = 129n = 186
 HRV2.31.6<.001
 HEV0.40.2.04
 HPeV0.40.5.93
Nasal swab and/or stooln = 132n = 192
 HRV2.61.9.002
 HEV0.40.2.04

Abbreviations: AOM, acute otitis media; HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus.

aP values base on Mann-Whitney U test.

bNumbers (n = 95, n = 144, etc) represent numbers of children.

Table 3.

Mean Number of Virus-Positive Nasal Swab and Stool Samples Among Children With or Without ≥1 Acute Otitis Media Episode by Age 1 Year

Virus-Positive Samples, Mean No.P Valuea
Sample Type and Virus≥1 AOM EpisodeNo AOM Episodes
Nasal swabn = 95bn = 144b
 HRV0.50.5.50
 HEV<0.1<0.1.70
Stooln = 129n = 186
 HRV2.31.6<.001
 HEV0.40.2.04
 HPeV0.40.5.93
Nasal swab and/or stooln = 132n = 192
 HRV2.61.9.002
 HEV0.40.2.04
Virus-Positive Samples, Mean No.P Valuea
Sample Type and Virus≥1 AOM EpisodeNo AOM Episodes
Nasal swabn = 95bn = 144b
 HRV0.50.5.50
 HEV<0.1<0.1.70
Stooln = 129n = 186
 HRV2.31.6<.001
 HEV0.40.2.04
 HPeV0.40.5.93
Nasal swab and/or stooln = 132n = 192
 HRV2.61.9.002
 HEV0.40.2.04

Abbreviations: AOM, acute otitis media; HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus.

aP values base on Mann-Whitney U test.

bNumbers (n = 95, n = 144, etc) represent numbers of children.

The time association between virus infections and AOM was studied by comparing the virus detection rate in samples that were collected either within 2 weeks before or after the onset of AOM (101 NS and 184 stool samples) with that in all other samples of the same type (Table 4). HEV detection in NS and stool samples was strongly associated with simultaneous AOM. Furthermore, detection of HRV in stool samples was significantly associated with simultaneous AOM. A similar trend was seen for HPeV, but the association was not statistically significant. The time associations between HRV, HEV, and HPeV detection and AOM episodes remained the same when a longer time window (±3 weeks) was used in these analyses (data not shown).

Table 4.

Proportions of Virus-Positive Samples Obtained From Participants With or Without Acute Otitis Media (AOM) and Association Between Virus Positivity and AOM

AOMa*No AOM OR (95% CI)
Sample Type and VirusSamples, No.Positive Samples, No. (%)Samples, No.Positive Samples, No. (%)UnadjustedAdjusted (GLM)Adjusted (GEEs)
Nasal swab
 HRV10128 (27.7)1286274 (21.3)1.42 (.90–2.23)1.38 (.88–2.18)b1.36 (.85–2.18)b
 HEV1015 (5.0)128620 (1.6)3.30 (1.21–8.98)3.05 (1.12–8.35)b3.26 (1.31–8.12)b
Stool
 HRV18450 (27.2)2677544 (20.3)1.46 (1.04–2.05)1.57 (1.10–2.24)c1.54 (1.04–2.30)c
 HEV18412 (6.5)267771 (2.7)2.56 (1.36–4.81)1.94 (1.01–3.72)c2.05 (1.10–3.80)c
 HPeV18415 (8.2)2668130 (4.9)1.73 (.99–3.02)1.44 (.81–2.56)d1.45 (.84–2.51)d
Nasal swab and/or stool e
 HRV22763 (27.8)3309654 (19.8)1.56 (1.15–2.11)1.53 (1.12–2.09)c1.40 (.96–2.04)c
 HEV22715 (6.6)330979 (2.4)2.89 (1.64–5.11)2.15 (1.19–3.88)c2.04 (1.06–3.91)c
AOMa*No AOM OR (95% CI)
Sample Type and VirusSamples, No.Positive Samples, No. (%)Samples, No.Positive Samples, No. (%)UnadjustedAdjusted (GLM)Adjusted (GEEs)
Nasal swab
 HRV10128 (27.7)1286274 (21.3)1.42 (.90–2.23)1.38 (.88–2.18)b1.36 (.85–2.18)b
 HEV1015 (5.0)128620 (1.6)3.30 (1.21–8.98)3.05 (1.12–8.35)b3.26 (1.31–8.12)b
Stool
 HRV18450 (27.2)2677544 (20.3)1.46 (1.04–2.05)1.57 (1.10–2.24)c1.54 (1.04–2.30)c
 HEV18412 (6.5)267771 (2.7)2.56 (1.36–4.81)1.94 (1.01–3.72)c2.05 (1.10–3.80)c
 HPeV18415 (8.2)2668130 (4.9)1.73 (.99–3.02)1.44 (.81–2.56)d1.45 (.84–2.51)d
Nasal swab and/or stool e
 HRV22763 (27.8)3309654 (19.8)1.56 (1.15–2.11)1.53 (1.12–2.09)c1.40 (.96–2.04)c
 HEV22715 (6.6)330979 (2.4)2.89 (1.64–5.11)2.15 (1.19–3.88)c2.04 (1.06–3.91)c

Abbreviations: AOM, acute otitis media; CI, confidence interval; GEEs, generalized estimation equations; GLM, generalized linear model; HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus; OR, odds ratio.

aAOM within 2 weeks before or after sampling.

bThe final model included virus positivity and having siblings.

cThe final model included virus positivity, being firstborn, and age at time of sampling.

dThe final model included being firstborn and age at time of sampling. Virus positivity was added as an additional variable.

eIncludes only samples obtained during the first 12 months of life.

Table 4.

Proportions of Virus-Positive Samples Obtained From Participants With or Without Acute Otitis Media (AOM) and Association Between Virus Positivity and AOM

AOMa*No AOM OR (95% CI)
Sample Type and VirusSamples, No.Positive Samples, No. (%)Samples, No.Positive Samples, No. (%)UnadjustedAdjusted (GLM)Adjusted (GEEs)
Nasal swab
 HRV10128 (27.7)1286274 (21.3)1.42 (.90–2.23)1.38 (.88–2.18)b1.36 (.85–2.18)b
 HEV1015 (5.0)128620 (1.6)3.30 (1.21–8.98)3.05 (1.12–8.35)b3.26 (1.31–8.12)b
Stool
 HRV18450 (27.2)2677544 (20.3)1.46 (1.04–2.05)1.57 (1.10–2.24)c1.54 (1.04–2.30)c
 HEV18412 (6.5)267771 (2.7)2.56 (1.36–4.81)1.94 (1.01–3.72)c2.05 (1.10–3.80)c
 HPeV18415 (8.2)2668130 (4.9)1.73 (.99–3.02)1.44 (.81–2.56)d1.45 (.84–2.51)d
Nasal swab and/or stool e
 HRV22763 (27.8)3309654 (19.8)1.56 (1.15–2.11)1.53 (1.12–2.09)c1.40 (.96–2.04)c
 HEV22715 (6.6)330979 (2.4)2.89 (1.64–5.11)2.15 (1.19–3.88)c2.04 (1.06–3.91)c
AOMa*No AOM OR (95% CI)
Sample Type and VirusSamples, No.Positive Samples, No. (%)Samples, No.Positive Samples, No. (%)UnadjustedAdjusted (GLM)Adjusted (GEEs)
Nasal swab
 HRV10128 (27.7)1286274 (21.3)1.42 (.90–2.23)1.38 (.88–2.18)b1.36 (.85–2.18)b
 HEV1015 (5.0)128620 (1.6)3.30 (1.21–8.98)3.05 (1.12–8.35)b3.26 (1.31–8.12)b
Stool
 HRV18450 (27.2)2677544 (20.3)1.46 (1.04–2.05)1.57 (1.10–2.24)c1.54 (1.04–2.30)c
 HEV18412 (6.5)267771 (2.7)2.56 (1.36–4.81)1.94 (1.01–3.72)c2.05 (1.10–3.80)c
 HPeV18415 (8.2)2668130 (4.9)1.73 (.99–3.02)1.44 (.81–2.56)d1.45 (.84–2.51)d
Nasal swab and/or stool e
 HRV22763 (27.8)3309654 (19.8)1.56 (1.15–2.11)1.53 (1.12–2.09)c1.40 (.96–2.04)c
 HEV22715 (6.6)330979 (2.4)2.89 (1.64–5.11)2.15 (1.19–3.88)c2.04 (1.06–3.91)c

Abbreviations: AOM, acute otitis media; CI, confidence interval; GEEs, generalized estimation equations; GLM, generalized linear model; HEV, human enterovirus; HPeV, human parechovirus; HRV, human rhinovirus; OR, odds ratio.

aAOM within 2 weeks before or after sampling.

bThe final model included virus positivity and having siblings.

cThe final model included virus positivity, being firstborn, and age at time of sampling.

dThe final model included being firstborn and age at time of sampling. Virus positivity was added as an additional variable.

eIncludes only samples obtained during the first 12 months of life.

The PAR percentage was calculated for each virus, using the adjusted ORs shown in Table 4 (combined stool and NS sample results). Based on these estimates, HRV and HEV could explain quite a significant proportion of AOM episodes, HRV slightly more (25%) than HEV (20%), whereas HPeV was less significant (11%).

Discussion

This study aimed to assess the frequency of HRV, HEV, and HPeV infections and their association with AOM in a cohort of Finnish children who were followed up from birth to age 3 years. The study setting allowed the role of these viruses in AOM to be evaluated in an unbiased manner in a child population that was not selected based on hospitalization, diagnosis of AOM, or respiratory infection. This offered an opportunity to evaluate the contribution of HRV, HEV, and HPeV infections to the development of AOM, taking into account the relatively high frequency of these infections in the background population, including children in whom AOM does not develop.

Our findings suggest that HRV and HEV are associated with the development of AOM. According to the PAR percentage calculations, HRVs could theoretically contribute to the development of AOM in up to 25% of AOM cases, and HEVs in up to 20%. HPeVs seem to be less important, because they showed a statistically nonsignificant trend that would equal 10% PAR. Even though the PAR prediction models are not accurate and should be interpreted with caution, the results indicate that the role of these viruses can be significant, particularly during the epidemic seasons. The results also demonstrate the high frequency of HRVs in young children generally, which may lead to overestimation of their role in AOM if a representative control group is not analyzed at the same time.

Overall, our findings regarding the incidence of studied virus infections and their role in the development of AOM are in line with previous reports [9, 14–18, 24, 25]. Of these 3 viruses, HRVs were detected most frequently, with 82.4% and 90.1% of the children testing positive for HRV at least once during the first year of life and by age 3 years, respectively. The incidences of HEV and HPeV infections were lower than that of HRV infections: by age 1 year HEV and HPeV was detected at least once in 18.2% and 29.2% of the participants, respectively, and 21.9% of the participants experienced ≥1 HEV infection by age 3 years. Even though HEV was less common and detected less frequently than HRV or HPeV in samples taken within 2 weeks of an AOM episode, HEV showed the strongest association with the development of AOM. This finding is supported by the results of a prospective, longitudinal cohort study using symptom-driven episodic sampling, which suggested that URI caused by HEV is complicated by AOM slightly more often than URI caused by HRV (34.4% vs 30.4%) [4]. On the other hand, a more recent prospective study among infants based on specimens collected monthly and during URI found HRV to be significantly associated with URI symptoms and AOM development, whereas HEV detection was associated with AOM development only marginally [2, 9].

The incidence of HEV infections in our study was somewhat lower than that found in some other studies [26]. This may partly be explained by our sampling scheme: from age 1 year onward, HEV detection was based solely on NS samples that were collected quite infrequently (every 6–12 months); hence, some HEV infections may have been missed. On the other hand, our finding is in line with those of previous studies suggesting that HEV infections are less common in Finland than in most other European countries [27].

Our findings regarding the incidence of HRV infections agree with those of previous studies, which have shown HRVs to be the most common cause of URIs in early childhood [28]. In the Finnish Otitis Media Cohort Study, 79%–91% of participating children had ≥1 HRV infection during the first 2 years of life, with >20% of children experiencing their first HRV episode by age 6 months [29]. In a recent birth cohort study conducted in Southwestern Finland, HRV was associated with 50% of AOM episodes diagnosed among children aged 0–2 years [30]. Similarly, in the previously mentioned study carried out among 0–12-month-old infants living in the United States, HRV was the most common virus detected during URI, and it was associated with the development of AOM, which complicated 27% of all URI episodes [2]. However, as described above, owing to the high frequency of HRVs in children, the frequency of HRVs during AOM episodes should be compared with those in a representative control group, and most previous studies have not made such comparisons.

One additional factor limiting the ability to draw conclusions from existing literature is the variability in the sensitivity and specificity of RT-qPCR assays used to detect HRV and HEV infections, particularly because both HRVs and HEVs are amplified by the same commonly used primer pairs targeting the conserved 5’ noncoding region of viral genomes. The RT-qPCR assays used in the present study were specific and sensitive for these viruses, as verified repeatedly by best possible scores in external HRV and HEV quality control rounds (Quality Control for Molecular Diagnostics). These assays used primers and probes that were unique to each virus, and therefore no cross-amplification occurred between the viruses.

Previous studies have found that HPeV may also contribute to the development of AOM in at least some cases [16–18]. In the current study, no significant association was seen between HPeV and AOM. This may be due to the relatively small size of our data set. Our study found HEV and HPeV infections to be equally common during the first year of life, because 29% of the participants tested positive for HPeV at least once by age 12 months. This frequency is well in line with the results of 2 birth cohort studies previously conducted in Finland [18, 31].

Our study has some limitations. Because all parental reports of AOM were counted as separate episodes regardless of the time period between consecutive reports, the incidence of AOM may have been overestimated. However, the incidence of AOM during the 3-year follow-up period (1 episode per person-year) agrees well with findings in previous studies in which the incidence of AOM in pediatric populations has varied depending on the study setting and geographic location, ranging from 161 and 256 events per 1000 person-years in Eastern and Western European countries, respectively, to the high incidence of 1470–1540 events per 1000 person-years observed in Finland [14, 32, 33].

The sampling scheme and diagnostic methods applied might be seen as another potential limitation of our study. On the other hand, regular sampling eliminated biases that could have been introduced by symptom-driven episodic sampling. In any case, as viruses were detected by RT-qPCR from samples collected at predefined 1–12-month intervals, some infections were probably missed owing to relatively short shedding of the virus. However, the excretion period of HPeV is relatively long, ranging from 2 to 24 weeks, and monthly collected stool samples have been suggested to be a fairly good method for detecting HPeV infections [31, 34]. Hence, we can assume that our study provides a good estimate for the incidence of HPeV infections among Finnish infants. The shedding periods for HRV and HEV, on the other hand, are thought to be shorter, lasting up to 5–6 weeks for HRV and 2–3 weeks for HEV in NS samples, although longer excretion periods for HEV in stool samples have also been reported [26, 35].

Stool specimens have been considered the classic sample type for HEVs and HPeVs, replicating primarily in the intestinal and oropharyngeal mucosa, but HRVs have also been detected frequently in stool samples [36–38]. A more frequent sampling scheme, including stool samples after 1 year of age or additional use of serologic methods, might have provided a more accurate estimate of the incidence of HRV and HEV infections during the second and third years of life.

Because viruses have been observed to play a significant role in the development of AOM, preventing viral URIs would be beneficial. A systematic review suggested that the influenza vaccine may reduce the incidence of AOM in infants and small children [39]. Our study and previous reports indicate that picornaviruses should also be targeted to prevent AOM in young children. Indeed, a previous study by Seppälä et al [40] implied that the oral polio vaccine, a live attenuated enterovirus vaccine, may reduce the incidence of AOM by providing cross-protection from nonpolio HEVs and other picornaviruses that have been associated with AOM. Although vaccines targeting certain HEVs have already been developed, and currently new HEV vaccines are under clinical development [41], attempts to produce an HRV vaccine have failed [42, 43]. However, recent studies from animal models suggest that it might be possible to develop a multivalent HRV vaccine that can induce immunity against several different HRV types [44]. Altogether, the development of vaccines against AOM-associated picornaviruses remains an important goal in future research, owing to the immense medical and health-economic burden caused by AOM.

In conclusion, our findings in this population-based birth cohort study support previous observations, suggesting that HRVs and HEVs contribute to the development of AOM in young children. HEV showed the strongest association with the development of AOM. These results suggest that antiviral strategies targeted at picornaviruses might be beneficial in preventing AOM in this age group.

Notes

Acknowledgments. The authors thank all the participating children and families in the DIABIMMUNE study, and the physicians, nurses, and technicians in the DIABIMMUNE study centers and laboratories.

Financial support. This work was supported by the European Union Seventh Framework Programme (grant 202063), the Academy of Finland (decision 292538), the Centre of Excellence in Molecular Systems Immunology and Physiology Research (decision 250114), the Liv och Hälsa Fund, Finska Läkaresällskapet, and the Sigrid Juselius Foundation.

Potential conflicts of interest. M. K. and H. H. are shareholders and board members of Vactech, which develops vaccines against picornaviruses. All other authors report no potential conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

References

1.

GBD 2016 Disease and Injury Incidence and Prevalence Collaborators
.
Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016
.
Lancet
2017
;
390
:
1211
59
.

2.

Chonmaitree
T
,
Alvarez-Fernandez
P
,
Jennings
K
, et al.
Symptomatic and asymptomatic respiratory viral infections in the first year of life: association with acute otitis media development
.
Clin Infect Dis
2015
;
60
:
1
9
.

3.

Winther
B
,
Alper
CM
,
Mandel
EM
,
Doyle
WJ
,
Hendley
JO
.
Temporal relationships between colds, upper respiratory viruses detected by polymerase chain reaction, and otitis media in young children followed through a typical cold season
.
Pediatrics
2007
;
119
:
1069
75
.

4.

Chonmaitree
T
,
Revai
K
,
Grady
JJ
, et al.
Viral upper respiratory tract infection and otitis media complication in young children
.
Clin Infect Dis
2008
;
46
:
815
23
.

5.

Kalu
SU
,
Ataya
RS
,
McCormick
DP
,
Patel
JA
,
Revai
K
,
Chonmaitree
T
.
Clinical spectrum of acute otitis media complicating upper respiratory tract viral infection
.
Pediatr Infect Dis J
2011
;
30
:
95
9
.

6.

Kvaerner
KJ
,
Nafstad
P
,
Hagen
JA
,
Mair
IW
,
Jaakkola
JJ
.
Recurrent acute otitis media: the significance of age at onset
.
Acta Otolaryngol
1997
;
117
:
578
84
.

7.

Niemelä
M
,
Uhari
M
,
Möttönen
M
,
Pokka
T
.
Costs arising from otitis media
.
Acta Paediatr
1999
;
88
:
553
6
.

8.

Rovers
MM
,
Schilder
AG
,
Zielhuis
GA
,
Rosenfeld
RM
.
Otitis media
.
Lancet
2004
;
363
:
465
73
.

9.

Ruohola
A
,
Meurman
O
,
Nikkari
S
, et al.
Microbiology of acute otitis media in children with tympanostomy tubes: prevalences of bacteria and viruses
.
Clin Infect Dis
2006
;
43
:
1417
22
.

10.

Ruohola
A
,
Pettigrew
MM
,
Lindholm
L
, et al.
Bacterial and viral interactions within the nasopharynx contribute to the risk of acute otitis media
.
J Infect
2013
;
66
:
247
54
.

11.

Nokso-Koivisto
J
,
Marom
T
,
Chonmaitree
T
.
Importance of viruses in acute otitis media
.
Curr Opin Pediatr
2015
;
27
:
110
5
.

12.

Kusel
MM
,
de Klerk
NH
,
Holt
PG
,
Kebadze
T
,
Johnston
SL
,
Sly
PD
.
Role of respiratory viruses in acute upper and lower respiratory tract illness in the first year of life: a birth cohort study
.
Pediatr Infect Dis J
2006
;
25
:
680
6
.

13.

van der Zalm
MM
,
Uiterwaal
CS
,
Wilbrink
B
, et al.
Respiratory pathogens in respiratory tract illnesses during the first year of life: a birth cohort study
.
Pediatr Infect Dis J
2009
;
28
:
472
6
.

14.

Nokso-Koivisto
J
,
Räty
R
,
Blomqvist
S
, et al.
Presence of specific viruses in the middle ear fluids and respiratory secretions of young children with acute otitis media
.
J Med Virol
2004
;
72
:
241
8
.

15.

Wiertsema
SP
,
Chidlow
GR
,
Kirkham
LA
, et al.
High detection rates of nucleic acids of a wide range of respiratory viruses in the nasopharynx and the middle ear of children with a history of recurrent acute otitis media
.
J Med Virol
2011
;
83
:
2008
17
.

16.

Pajkrt
D
,
Benschop
KS
,
Westerhuis
B
,
Molenkamp
R
,
Spanjerberg
L
,
Wolthers
KC
.
Clinical characteristics of human parechoviruses 4-6 infections in young children
.
Pediatr Infect Dis J
2009
;
28
:
1008
10
.

17.

Sillanpää
S
,
Oikarinen
S
,
Sipilä
M
, et al.
Human parechovirus as a minor cause of acute otitis media in children
.
J Clin Virol
2015
;
62
:
106
9
.

18.

Tauriainen
S
,
Oikarinen
S
,
Taimen
K
, et al.
Temporal relationship between human parechovirus 1 infection and otitis media in young children
.
J Infect Dis
2008
;
198
:
35
40
.

19.

Nokso-Koivisto
J
,
Kinnari
TJ
,
Lindahl
P
,
Hovi
T
,
Pitkäranta
A
.
Human picornavirus and coronavirus RNA in nasopharynx of children without concurrent respiratory symptoms
.
J Med Virol
2002
;
66
:
417
20
.

20.

Simre
K
,
Uibo
O
,
Peet
A
, et al. ;
DIABIMMUNE Study Group
.
Exploring the risk factors for differences in the cumulative incidence of coeliac disease in two neighboring countries: the prospective DIABIMMUNE study
.
Dig Liver Dis
2016
;
48
:
1296
301
.

21.

Honkanen
H
,
Oikarinen
S
,
Peltonen
P
, et al.
Human rhinoviruses including group C are common in stool samples of young Finnish children
.
J Clin Virol
2012
;
56
:
334
8
.

22.

Honkanen
H
,
Oikarinen
S
,
Pakkanen
O
, et al.
Human enterovirus 71 strains in the background population and in hospital patients in Finland
.
J Clin Virol
2013
;
56
:
348
53
.

23.

Nix
WA
,
Maher
K
,
Johansson
ES
, et al.
Detection of all known parechoviruses by real-time PCR
.
J Clin Microbiol
2008
;
46
:
2519
24
.

24.

Seppälä
E
,
Sillanpää
S
,
Nurminen
N
, et al.
Human enterovirus and rhinovirus infections are associated with otitis media in a prospective birth cohort study
.
J Clin Virol
2016
;
85
:
1
6
.

25.

Yano
H
,
Okitsu
N
,
Hori
T
, et al.
Detection of respiratory viruses in nasopharyngeal secretions and middle ear fluid from children with acute otitis media
.
Acta Otolaryngol
2009
;
129
:
19
24
.

26.

Witsø
E
,
Palacios
G
,
Cinek
O
, et al.
High prevalence of human enterovirus a infections in natural circulation of human enteroviruses
.
J Clin Microbiol
2006
;
44
:
4095
100
.

27.

Viskari
H
,
Ludvigsson
J
,
Uibo
R
, et al.
Relationship between the incidence of type 1 diabetes and enterovirus infections in different European populations: results from the EPIVIR project
.
J Med Virol
2004
;
72
:
610
7
.

28.

Kieninger
E
,
Fuchs
O
,
Latzin
P
,
Frey
U
,
Regamey
N
.
Rhinovirus infections in infancy and early childhood
.
Eur Respir J
2013
;
41
:
443
52
.

29.

Nokso‐Koivisto
J
,
Räty
R
,
Blomqvist
S
, et al.
Presence of specific viruses in the middle ear fluids and respiratory secretions of young children with acute otitis media
.
J Med Virol
2004
;
72
:
241
8
.

30.

Toivonen
L
,
Schuez-Havupalo
L
,
Karppinen
S
, et al.
Rhinovirus infections in the first 2 years of life
.
Pediatrics
2016
;
138
:
e20161309
.

31.

Kolehmainen
P
,
Oikarinen
S
,
Koskiniemi
M
, et al.
Human parechoviruses are frequently detected in stool of healthy Finnish children
.
J Clin Virol
2012
;
54
:
156
61
.

32.

Liese
JG
,
Silfverdal
SA
,
Giaquinto
C
, et al.
Incidence and clinical presentation of acute otitis media in children aged <6 years in European medical practices
.
Epidemiol Infect
2014
;
142
:
1778
88
.

33.

Usonis
V
,
Jackowska
T
,
Petraitiene
S
, et al.
Incidence of acute otitis media in children below 6 years of age seen in medical practices in five East European countries
.
BMC Pediatr
2016
;
16
:
108
.

34.

Wildenbeest
JG
,
Benschop
KS
,
Bouma-de Jongh
S
,
Wolthers
KC
,
Pajkrt
D
.
Prolonged shedding of human parechovirus in feces of young children after symptomatic infection
.
Pediatr Infect Dis J
2016
;
35
:
580
3
.

35.

Jartti
T
,
Lehtinen
P
,
Vuorinen
T
,
Koskenvuo
M
,
Ruuskanen
O
.
Persistence of rhinovirus and enterovirus RNA after acute respiratory illness in children
.
J Med Virol
2004
;
72
:
695
9
.

36.

de Crom
SC
,
Obihara
CC
,
de Moor
RA
,
Veldkamp
EJ
,
van Furth
AM
,
Rossen
JW
.
Prospective comparison of the detection rates of human enterovirus and parechovirus RT-qPCR and viral culture in different pediatric specimens
.
J Clin Virol
2013
;
58
:
449
54
.

37.

de Crom
SCM
,
Rossen
JWA
,
van Furth
AM
,
Obihara
CC
.
Enterovirus and parechovirus infection in children: a brief overview
.
Eur J Pediatr
2016
;
175
:
1023
9
.

38.

Savolainen-Kopra
C
,
Simonen-Tikka
ML
,
Klemola
P
, et al.
Human rhinoviruses in INDIS-study material-evidence for recovery of viable rhinovirus from fecal specimens
.
J Med Virol
2013
;
85
:
1466
72
.

39.

Norhayati
MN
,
Ho
JJ
,
Azman
MY
.
Influenza vaccines for preventing acute otitis media in infants and children
.
Cochrane Database Syst Rev
2017
;
10
:
CD010089
.

40.

Seppälä
E
,
Viskari
H
,
Hoppu
S
, et al.
Viral interference induced by live attenuated virus vaccine (OPV) can prevent otitis media
.
Vaccine
2011
;
29
:
8615
8
.

41.

Hyöty
H
,
Leon
F
,
Knip
M
.
Developing a vaccine for type 1 diabetes by targeting coxsackievirus B
.
Expert Rev Vaccines
2018
;
17
:
1071
83
.

42.

McLean
GR
.
Developing a vaccine for human rhinoviruses
.
J Vaccines Immun
2014
;
2
:
16
20
.

43.

Yi
EJ
,
Shin
YJ
,
Kim
JH
,
Kim
TG
,
Chang
SY
.
Enterovirus 71 infection and vaccines
.
Clin Exp Vaccine Res
2017
;
6
:
4
14
.

44.

Lee
S
,
Nguyen
MT
,
Currier
MG
, et al.
A polyvalent inactivated rhinovirus vaccine is broadly immunogenic in rhesus macaques
.
Nat Commun
2016
;
7
:
12838
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/journals/pages/open_access/funder_policies/chorus/standard_publication_model)