Evaluation of 15 Functional Candidate Genes for Association with Chronic Otitis Media with Effusion and/or Recurrent Otitis Media (COME/ROM)

Abstract

DNA sequence variants in genes involved in the innate immune response and secondary response to infection may confer susceptibility to chronic otitis media with effusion and/or recurrent otitis media (COME/ROM). We evaluated single nucleotide polymorphisms (SNPs) in 15 functional candidate genes. A total of 99 SNPs were successfully genotyped on the Sequenom platform in 142 families (618 subjects) from the Minnesota COME/ROM Family Study. Data were analyzed for association with COME/ROM using the Generalized Disequilibrium Test (GDT). Sex and age at exam were adjusted as covariates, relatedness was accounted for, and genotype differences from all phenotypically discordant relative pairs were utilized to measure the evidence of association between COME/ROM and each SNP. SNP rs2735733 in the region of the mucin 5, subtypes A/C gene (MUC5AC) exhibited nominal evidence for association with COME/ROM (P = 0.002). Two additional SNPs from this region had P values<0.05. Other variants exhibiting associations with COME/ROM at P<0.05 included the SCN1B SNP rs8100085 (P = 0.013), SFTPD SNP rs1051246 (P = 0.039) and TLR4 SNP rs2770146 (P = 0.038). However, none of these associations replicated in an independent sample of COME/ROM families. The candidate gene variants examined do not appear to make a major contribution to COME/ROM susceptibility, despite a priori evidence from functional or animal model studies for a role in COME/ROM pathology.

About the Authors:

Michèle M. Sale

* E-mail: [email protected]

Affiliations Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, United States of America, Department of Medicine, University of Virginia, Charlottesville, Virginia, United States of America, Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, Virginia, United States of America

Wei-Min Chen

Affiliations Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, United States of America, Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Daniel E. Weeks

Affiliations Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America, Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

Josyf C. Mychaleckyj

Affiliations Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, United States of America, Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Xuanlin Hou

Affiliations Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, United States of America, Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Miranda Marion

Affiliation: Department of Biostatistical Sciences, Wake Forest University School of Medicine, Winston-Salem, North Carolina, United States of America

Fernando Segade

Affiliation: Department of Anatomy and Cell Biology, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Margaretha L. Casselbrant

Affiliations Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America, Department of Otolaryngology, Children's Hospital of Pittsburgh of UPMC, Pittsburgh, Pennsylvania, United States of America

Ellen M. Mandel

Affiliation: Department of Otolaryngology, Children's Hospital of Pittsburgh of UPMC, Pittsburgh, Pennsylvania, United States of America

Robert E. Ferrell

Affiliation: Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

Stephen S. Rich

Affiliations Center for Public Health Genomics, University of Virginia, Charlottesville, Virginia, United States of America, Department of Medicine, University of Virginia, Charlottesville, Virginia, United States of America, Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Kathleen A. Daly

Affiliation: Department of Otolaryngology, University of Minnesota, Minneapolis, Minnesota, United States of America

Introduction

Incidence and prevalence of otitis media (OM) has declined since 2000 when infant immunization with pneumococcal conjugate vaccine was adopted in the United States [1], [2], [3]. However, OM is still a very common childhood disease and remains a major cause of morbidity and hearing loss in young children [4], [5]. Children with an affected sibling are at higher risk, and OM clusters in families across generations [6], [7]. Twin studies in the United States, Norway, England and Wales have reported a high degree of heritability for recurrent and chronic OM [8], [9], [10], defined as more than 3 episodes of acute OM in a year, or middle ear fluid persisting for 3–4 months [11], [12].

Our previous research on genetic contributions to chronic otitis media with effusion and/or recurrent otitis media (COME/ROM) include a complex segregation analysis of 173 families suggesting single gene effects [13], a genome-wide linkage scan in 121 families providing support for loci on chromosomes 10q and 19q [14], and an investigation of the human ortholog (FBXO11) of the deaf mouse mutant jeff gene, a single gene model of chronic/recurrent OM in our family-based study [15].

We selected 15 candidate genes involved in the innate immune response or the ability to clear infectious agents implicated in OM, and evaluated common variation across the coding regions of these genes for association with COME/ROM.

Methods

Ethics statement

This study was conducted under Institutional Review Board approval at the University of Minnesota, Wake Forest University, and University of Pittsburgh, and adhered to the tenets of the Declaration of Helsinki. Written consent was provided by all participants or, in the case of children, by their parent or legal guardian.

Subjects

Subjects who had tympanostomy tube surgery for COME/ROM (probands) and their families were recruited for the study, which has been described previously [14], [15], [16]. All probands were treated with tympanostomy tubes. Participating first and second degree relatives and parents of probands were classified as affected or unaffected based on 1) treatment with tympanostormy tubes, or 2) four data sources. Participants or parents completed a history form documenting history and risk factors for COME/ROM. The study otolaryngologist performed an ear examination to determine presence of OM sequelae without knowledge of the subject's prior OM history, an audiologist conducted tympanometric testing at three frequencies (226, 630 or 710, and 1400) to detect abnormal middle ear mechanics, and medical records were abstracted to document tympanic membrane and middle ear findings consistent with COME/ROM, and duration of OME. Hearing was screened at 20 dB for speech frequencies. Phenotype information comes from four sources: self or parent reported OM history, ear examination, tympanometric findings, and medical record abstract (see Material S1). Affected status of probands' relatives was based on having positive findings from at least two of these four sources (Table S1). A total of 142 families with at least two individuals with the diagnosis of COME/ROM were used in the candidate gene analyses.

We carried out a replication analysis in an independent study of otitis media [17] that consisted of 1,583 genotyped individuals from 441 Caucasian families. In order to assure a history of significant ear disease, two or more full siblings who both or all had undergone tympanostomy tube insertion were enrolled. The need for tympanostomy tube insertion established that the subject's history of middle ear disease was truly significant, resulting in the need for a surgical procedure. A subject was only considered “affected” if he/she had undergone tympanostomy tube insertion at least once for recurrent/persistent OM, while a subject was considered “unaffected” if he/she had never had tympanostomy tubes and had no known history of recurrent/persistent OM. The remaining subjects were considered as having “unknown” disease status. Otoscopic examinations and tympanograms were conducted at at entry if requested by the parent for patient information not study data since the condition of the ears at entry did not determine eligibility and the tubes may have been inserted many years prior to study entry.

SNP selection and genotyping

DNA was isolated from blood using the Gentra PureGene method (Qiagen, Valencia CA). Candidate genes were selected on the basis of a role in innate host defense and acute inflammatory response or chronic response and effusion, and prior evidence in the literature for involvement in OM. We selected the largest known isoform, added 5 kb upstream and downstream of the coding region, and used HapMap CEU (Centre d'Etude Polymorphisme Human) data (HapMap Data Rel#21/phaseII Jul06, on NBCI B35 assembly, dbSNP b125) to identify SNPs with minor allele frequency (MAF)≥0.05. For mucin 5, subtypes A/C (MUC5AC), we used the combined largest mRNA sequence (AF015521; AJ298317; AF043909), which encompasses MUC5AC NC_000011.8. We used the pairwise tagging option of Tagger [18] (r²>0.8) to select the miminal tagSNP set for each gene. Genotyping was conducted on the Sequenom platform using the iPlex assay [19]. A total of 99 SNPs were successfully genotyped in >95% of samples (see Table 1 for gene summary). Concordance rates were calculated on the basis of 21 replicates included for genotyping on the Sequenom platform at the same time that all samples were assayed. Fifty-four SNPs were 100% concordant; 19 SNPs had one discordant call; 25 SNPs had two discordant calls; and 1 SNP (rs778588 near CD14) had three discordant calls. The concordance rates for nominally associated SNPs were 100% (rs1051246, rs7396030, rs2075859, rs8100085) or 95% (rs2735733, rs2770146; equivalent to one discordant call out of 21 or 20 replicate pairs respectively).

[Figure omitted. See PDF.]

Table 1. Selected candidate genes and number of SNPs genotyped.

https://doi.org/10.1371/journal.pone.0022297.t001

Candidate genes selected on the basis of their function in innate host defense and acute inflammatory response included Toll-like receptor 4 (TLR4) [20], Toll-like receptor 2 (TLR2) [21], [22], Beta-defensin 1 (DEFB1) [23], Surfactant protein A1 (SFTPA1) [24], [25], [26], Surfactant protein D (SFTPD) [27], [28], Interleukin-8 (IL8) [29], [30], [31], Interleukin-1 β (IL1B) [31], Interleukin 10 (IL10) [20], and Lactotransferrin (LTF) [32]. Genes selected for their role in chronic response and effusion were mucin 2 (MUC2) [33], mucin 5, subtypes A/C (MUC5AC) [34], Mannose-binding lectin precursor (MBL2) [35], CD14 antigen precursor (CD14) [36], Sodium channel, voltage-gated, type I-β (SCN1B) [37], and Ecotropic viral integration site 1 (EVI1) [38].

Statistical analyses

Deviations from Hardy Weinberg equilibrium (HWE) in unaffected founder individuals were determined using the exact test [39]. Four SNPs with small HWE P values were excluded from the analysis and the remaining 99 SNPs have P>0.001. Data were checked for Mendelian errors using PedCheck [40] and pedstats [39] and poorly performing SNPs or a minimal dataset required to resolve errors were removed from analyses. Two monozygotic (MZ) twin pairs were detected and incorporated in the association analysis.

Data were analyzed for association with COME/ROM using the Generalized Disequilibrium Test (GDT) [41]. The GDT utilizes the genotype differences of all phenotypically discordant relative pairs in assessing association within families. This test has been shown to be consistently more powerful than PDT [42], FBAT [43] and several other family-based association tests for a common disease [41]. It allows modeling of covariates and IBD allele sharing, and information from extended pedigrees is efficiently used without breaking extended pedigrees into multiple nuclear families, as implemented in other methods. Compared to other association tests that also examine between-family association, the GDT method has the advantage of protecting from population stratification between families, and being less affected by genetic heterogeneity from multiple susceptibility genes. In GDT analyses, we modeled sex and age at exam (considered a proxy for generation/clinical practice) as covariates. Identity-by-descent (IBD) statistics were estimated using the Merlin package [44], and were incorporated in the GDT statistics.

Results

Population characteristics

The genotyped population included 618 individuals from 142 families (Table 2). One hundred thirty-two families were non-Hispanic and of European ancestry. The 10 non-European families consisted of seven families that described themselves as non-Hispanic and of mixed race, one mixed race Hispanic family, one Asian family, and one Native American family. As family-based tests of association using the GDT are robust to population stratification between families [41], all families were retained for analyses.

[Figure omitted. See PDF.]

Table 2. Participant characteristics.

https://doi.org/10.1371/journal.pone.0022297.t002

Associations with COME/ROM

A total of 99 SNPs encompassing 267 kb were successfully genotyped in subjects from the Minnesota COME/ROM Family Study. Age at exam was the factor most significantly associated with COME/ROM. The odds ratio of association between one year increment and COME/ROM was 0.94, with P = 1.2×10⁻¹⁶. This strong association suggests the importance of adjusting for the age effect in the candidate region association analysis. Although the sex effect is not statistically significant after adjusting for age, it was retained as a covariate in the candidate region association analyses.

SNP rs2735733 in the MUC5AC gene region exhibited nominal evidence (P = 0.002) for association with COME/ROM. At this SNP, the T allele (minor allele frequency (MAF) 0.461) is more frequent in unaffected individuals than in affected individuals, with the allele frequency difference being 0.106. The odds ratio (OR) of the minor allele was estimated to be 0.646, assuming a population prevalence of 0.1. Although mucin 5B (MUC5B) was not initially selected for study as an a priori candidate, the genomic region genotyped, including the associated region of MUC5AC, overlaps with MUC5B. In the region of chromosome 11 encompassing MUC5AC (MUC5B) and MUC2, there were 21 SNPs in total genotyped. Two other SNPs in this region had P values<0.05 (P = 0.049 at rs7396030 and P = 0.041 at rs2075859). The largest effect of a single SNP is OR = 1.57 at rs7396030. SNP rs7396030 is in weak LD with rs2075859 and rs2735733, with r² 0.036 and 0.013 respectively (D′ 0.293 and 0.216 respectively); rs2075859 and rs2735733 are in LD, with r² = 0.678 (and D′ = 0.987). Only founders were used in the LD calculation, conducted using PLINK [45]. Haplotype-based GDT analysis (as implemented in GDT [41]) identified association (although not significant) between haplotype CC at rs2075859 and rs2735733 and COME/ROM, with P value 0.0097, which is consistent with the single SNP association results. Further haplotype analysis based on three SNPs does not produce any significant association (smallest P = 0.43).

We identified several other nominal SNP-COME/ROM associations, including the SCN1B rs8100085 (P = 0.013), SFTPD rs1051246 (P = 0.039), and TLR4 rs2770146 (P = 0.038; Table 3).

[Figure omitted. See PDF.]

Table 3. GDT association results (P<0.05) for COME/ROM.

https://doi.org/10.1371/journal.pone.0022297.t003

We also carried out association analysis in families consisting of European American only (eight families were excluded) and observed the same associations. All six associated SNPs listed in Table 3 retain their significance with P values 0.027, 0.027, 0.045, 0.033, 0.002, and 0.022 respectively. An additional SNP (rs2672812) in the MUC5A/MUC5B region also had a P value <0.05 (P = 0.042).

We carried out a replication analysis of the six SNPs in Table 3 in an independent study [17]. The association results in the replication study are shown in Table 4. The allele frequencies in two studies are comparable. Only one SNP in the MUC2 region, rs7396030, reached a P-value <0.05 in the replication study (unadjusted P = 0.0075; P = 0.022 after adjusting for sex). However, the risk allele in the replication study is the major allele, opposite from our study. The allele flip could occur either by chance (with a small probability of 0.025), or as a genuine allele flip that is due to complex linkage disequilibrium with the causal variant in the region (Clarke and Cardon, 2010). The most likely explanation is that both results are false positives. In the replication sample, the LD structure across the MUC5AC/MUC5B/MUC2 region was similar to that observed in the University of Minnesota sample. LD between SNPs rs2075859 and rs2735733 was high (r² = 0.652, D′ = 0.988), but LD between these SNPs and MUC2 SNP rs7396030 was low (D′ = 0.022 and 0.034 respectively).

[Figure omitted. See PDF.]

Table 4. GDT association results for SNPs in Table 3 genotyped in the Caucasian subjects from the replication sample.

https://doi.org/10.1371/journal.pone.0022297.t004

Discussion

We performed an evaluation of 15 functional candidate genes in a well-characterized population of families from the Minnesota COME/ROM Family Study. We utilized available HapMap data and employed an LD-based tagging approach to survey common variation in these genes. In some cases, few tagging SNPs were available in the genic region at the time of SNP selection. The strongest observed associations were with SNPs in the mucin 5 genomic region, encompassing both MUC2 and MUC5AC/MUC5B.

Middle ear inflammation stimulates the production and release of inflammatory mediators, with subsequent upregulation of a suite of mucin genes [46]. Mucins are a family of glycosylated proteins that function as part of the mucociliary transport system within the middle ear, binding pathogens and helping clear cellular debris [46]. Muc5ac is upregulated in the OM rat model [47], and is one of the genes expressed in mucoid effusion from patients [48]. Similarly, muc5b is a major component of OME effusions in a rabbit model of OME [49], and is present in middle ear secretory cells of patients with COM [50]. Both MUC5AC and MUC5B are located on chromosome 11, with overlapping coding regions. One study found that a longer transcript of MUC5AC – but not variants of MUC5B or MUC2 – was associated with OME [51]. However, recent evidence showed MUC5B to be the predominant mucin identified in middle ear effusions from children with COME [52]. We plan to explore whether SNPs associated with COME/ROM in the present study (or variants in LD with them) influence isoform production or transcript stability.

Less robust associations were observed with SNPs in the SCN1B, SFTPD and TLR4 genes. Interestingly, in rats whose ears were inoculated with Streptococcus pneumoniae, one of the most common pathogens in otitis media, scn1b, muc2 and muc5 were among the genes upregulated [37], although polymorphisms of SCN1B do not appear to have been investigated in COME/ROM previously. Surfactant protein D (sftpd) knockout animals do not clear influenza A viral infections [53]. SFTPD polymorphisms in humans have been associated with surfactant protein D assembly, function, and concentration [54], as well as severe respiratory syncytial virus infections [27], a known precipitating factor for OM [55]. In a mouse model of spontaneous OM, an exon 3 missense mutation in Tlr4 resulted in lipid A/lipopolysaccharide (LPS) endotoxin insenstivity and an inability to clear Gram-negative bacteria such as Haemophilus influenzae [56], a bacterium frequent found in COME [57]. Similar defects in LPS response due to TLR4 mutations have been shown to exist in humans [58].

This study has several strengths and limitations. The primary strength is the family-based study design, permitting an analytical approach robust to population stratification. Our newly developed GDT approach also allowed us to take into account covariates age and sex. The study included a small sample of non-European families, consequently we are unable to draw any conclusions about whether the observed associations are specific to European American populations. Genotyping inaccuracy could potentially have led to false negative results however, as noted in the Methods, the concordance rates for the SNPs shown in Table 3 exceeded 95%. The candidate gene approach is necessarily limited by the need for evidence of a role in OM pathology, thus we are expanding our search for COME/ROM susceptibility alleles to a genome-wide association study.

Epidemiological studies suggest that six episodes of AOM by age six is not uncommon. From among individuals with medical record abstracts (65% of participants), only 8 subjects (7 probands and 1 sibling) had 6 episodes of OM by age 6, but did not have 3 episodes OM in any 12 month period. All 8 had tubes placement surgery and were also positive for a history of COME/ROM by either ear exam or tympanometry measures, so would have been classified as affected even without meeting more stringent ROM criteria. Although both studies used similar phenotype criteria to determine affected status based on reported history of COME/ROM, medical record review, ear examination and tympanograms, for the primary study, the sibling of the proband did not need to have undergone tympanostomy tube insertion. In contrast, in the replication sample, both siblings were required to have tympanostomy tube insertion. Failure to replicate nominal associations may be a consequence of diagnostic differences between samples, although it is probable that initial associations represent false positive results. Although we cannot rule out the possibility of type 2 error due to lack of power in the initial sample, we can conclude that there are no strong effects on COME/ROM risk by common variants genotyped in these genes.

Since we define COME as middle ear fluid that lasts 3 or more months, and ROM as 3 OM episodes in a year, or 6 episodes by age 6, there is considerable overlap in our sample between COME and ROM. It is not unreasonable to hypothesize that inadequate production of mucins may predispose to ROM whereas excessive production of mucins may predispose to COME, or that different loci contribute to one or the other. Often patients with middle ear fluid are prone to acute OM, so COME may predipose to ROM. In the UMN sample there were 143 participants with a history of both COME and ROM. Disentangling this relationship in future analyses would require a larger sample.

Our results suggest that common variants in the selected candidate genes do not appear to make a major contribution to COME/ROM susceptibility, despite strong a priori hypotheses for a role of these genes in disease pathology, suggesting agnostic genome-wide analyses of COME/ROM populations may be a more productive approach to identify susceptibility alleles.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Criteria for Classifying Family Members as Affected with COME/ROM for the University of Minnesota Study.

https://doi.org/10.1371/journal.pone.0022297.s001

(DOCX)

Material S1.

Criteria for Classifying Family Members as Affected with COME/ROM for the University of Minnesota Study.

https://doi.org/10.1371/journal.pone.0022297.s002

(DOCX)

Acknowledgments

We thank the families for their participation. The study would not have been possible without the recruitment team: Donna Brauer, Karen Savonna, Ruth Selvius, Ann Brauti, Cynthia Roh, Michele Bochert, Jennifer Palmer, Kendra Herrell, Judy Monroe, Norrita Rech, Kim Canfield, Heather Nelsonm and JoAnn Knox. Thanks also to Peter Perlegas, Dax Allred, and Pamela Hicks for technical assistance, as well as Donald Bowden, Ph.D., and Bronya Keats, Ph.D. for advice.

Author Contributions

Conceived and designed the experiments: MMS FS SSR KAD. Performed the experiments: MMS REF. Analyzed the data: W-MC XH DEW MM. Contributed reagents/materials/analysis tools: JCM XH MLC EMM KAD. Wrote the paper: MMS W-MC DEW SSR KAD.

Citation: Sale MM, Chen W-M, Weeks DE, Mychaleckyj JC, Hou X, Marion M, et al. (2011) Evaluation of 15 Functional Candidate Genes for Association with Chronic Otitis Media with Effusion and/or Recurrent Otitis Media (COME/ROM). PLoS ONE 6(8): e22297. https://doi.org/10.1371/journal.pone.0022297

References

1. Grijalva CG, Poehling KA, Nuorti JP, Zhu Y, Martin SW, et al. (2006) National impact of universal childhood immunization with pneumococcal conjugate vaccine on outpatient medical care visits in the United States. Pediatrics 118: 865–873.CG GrijalvaKA PoehlingJP NuortiY. ZhuSW Martin2006National impact of universal childhood immunization with pneumococcal conjugate vaccine on outpatient medical care visits in the United States.Pediatrics118865873

2. U.S. Department of Health and Human Services (2006) Chapter 28: Vision and Hearing. Healthy People 2010 Midcourse Review. Washington, D.C.: U.S. Government Printing Office. pp. 28–25.U.S. Department of Health and Human Services2006Chapter 28: Vision and Hearing. Healthy People 2010 Midcourse ReviewWashington, D.C.U.S. Government Printing Office2825

3. Zhou F, Shefer A, Kong Y, Nuorti JP (2008) Trends in acute otitis media-related health care utilization by privately insured young children in the United States, 1997–2004. Pediatrics 121: 253–260.F. ZhouA. SheferY. KongJP Nuorti2008Trends in acute otitis media-related health care utilization by privately insured young children in the United States, 1997–2004.Pediatrics121253260

4. Ho V, Daly KA, Hunter LL, Davey C (2002) Otoacoustic emissions and tympanometry screening among 0–5 year olds. Laryngoscope 112: 513–519.V. HoKA DalyLL HunterC. Davey2002Otoacoustic emissions and tympanometry screening among 0–5 year olds.Laryngoscope112513519

5. Hunter LL, Davey CS, Kohtz A, Daly KA (2007) Hearing screening and middle ear measures in American Indian infants and toddlers. Int J Pediatr Otorhinolaryngol 71: 1429–1438.LL HunterCS DaveyA. KohtzKA Daly2007Hearing screening and middle ear measures in American Indian infants and toddlers.Int J Pediatr Otorhinolaryngol7114291438

6. Teele DW, Klein JO, Rosner B (1989) Epidemiology of otitis media during the first seven years of life in children in greater Boston: a prospective, cohort study. J Infect Dis 160: 83–94.DW TeeleJO KleinB. Rosner1989Epidemiology of otitis media during the first seven years of life in children in greater Boston: a prospective, cohort study.J Infect Dis1608394

7. Rasmussen F (1993) Protracted secretory otitis media. The impact of familial factors and day-care center attendance. Int J Pediatr Otorhinolaryngol 26: 29–37.F. Rasmussen1993Protracted secretory otitis media. The impact of familial factors and day-care center attendance.Int J Pediatr Otorhinolaryngol262937

8. Kvaerner KJ, Tambs K, Harris JR, Magnus P (1997) Distribution and heritability of recurrent ear infections. Ann Otol Rhinol Laryngol 106: 624–632.KJ KvaernerK. TambsJR HarrisP. Magnus1997Distribution and heritability of recurrent ear infections.Ann Otol Rhinol Laryngol106624632

9. Casselbrant ML, Mandel EM, Fall PA, Rockette HE, Kurs-Lasky M, et al. (1999) The heritability of otitis media: a twin and triplet study. JAMA 282: 2125–2130.ML CasselbrantEM MandelPA FallHE RocketteM. Kurs-Lasky1999The heritability of otitis media: a twin and triplet study.JAMA28221252130

10. Rovers M, Haggard M, Gannon M, Koeppen-Schomerus G, Plomin R (2002) Heritability of symptom domains in otitis media: a longitudinal study of 1,373 twin pairs. Am J Epidemiol 155: 958–964.M. RoversM. HaggardM. GannonG. Koeppen-SchomerusR. Plomin2002Heritability of symptom domains in otitis media: a longitudinal study of 1,373 twin pairs.Am J Epidemiol155958964

11. Daly KA, Hunter LL, Giebink GS (1999) Chronic otitis media with effusion. Pediatr Rev 20: 85–93; quiz 94.KA DalyLL HunterGS Giebink1999Chronic otitis media with effusion.Pediatr Rev208593; quiz 94

12. Kvaerner KJ (2004) When is acute otitis media too often? In: Alper CM, Bluestone CD, Casselbrant ML, Dohar JE, Mandel EM, editors. Advanced Therapy of Otitis Media. Hamilton ON: BC Decker Inc. pp. 106–109.KJ Kvaerner2004When is acute otitis media too often?CM AlperCD BluestoneML CasselbrantJE DoharEM MandelAdvanced Therapy of Otitis MediaHamilton ONBC Decker Inc106109

13. Spray BJ, Daly KA, Sellers TA, Levine SC, Rich SS (2002) Evidence for single gene effects in the familial risk of chronic and recurrent otitis media. pp. 225–227.BJ SprayKA DalyTA SellersSC LevineSS Rich2002Evidence for single gene effects in the familial risk of chronic and recurrent otitis media.225227Proceedings of the Seventh International Symposium on Recent Advances in Otitis Media. Proceedings of the Seventh International Symposium on Recent Advances in Otitis Media.

14. Daly KA, Brown WM, Segade F, Bowden DW, Keats BJ, et al. (2004) Chronic and recurrent otitis media: a genome scan for susceptibility loci. Am J Hum Genet 75: 988–997.KA DalyWM BrownF. SegadeDW BowdenBJ Keats2004Chronic and recurrent otitis media: a genome scan for susceptibility loci.Am J Hum Genet75988997

15. Segade F, Daly KA, Allred D, Hicks PJ, Cox M, et al. (2006) Association of the FBXO11 gene with chronic otitis media with effusion and recurrent otitis media: the Minnesota COME/ROM Family Study. Arch Otolaryngol Head Neck Surg 132: 729–733.F. SegadeKA DalyD. AllredPJ HicksM. Cox2006Association of the FBXO11 gene with chronic otitis media with effusion and recurrent otitis media: the Minnesota COME/ROM Family Study.Arch Otolaryngol Head Neck Surg132729733

16. Daly KA, Rich SS, Levine S, Margolis RH, Le CT, et al. (1996) The family study of otitis media: design and disease and risk factor profiles. Genet Epidemiol 13: 451–468.KA DalySS RichS. LevineRH MargolisCT Le1996The family study of otitis media: design and disease and risk factor profiles.Genet Epidemiol13451468

17. Casselbrant ML, Mandel EM, Jung J, Ferrell RE, Tekely K, et al. (2009) Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions. BMC Med Genet 10: 85.ML CasselbrantEM MandelJ. JungRE FerrellK. Tekely2009Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions.BMC Med Genet1085

18. de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, et al. (2005) Efficiency and power in genetic association studies. Nat Genet 37: 1217–1223.PI de BakkerR. YelenskyI. Pe'erSB GabrielMJ Daly2005Efficiency and power in genetic association studies.Nat Genet3712171223

19. Jurinke C, van den Boom D, Cantor CR, Koster H (2002) Automated genotyping using the DNA MassArray technology. Methods Mol Biol 187: 179–192.C. JurinkeD. van den BoomCR CantorH. Koster2002Automated genotyping using the DNA MassArray technology.Methods Mol Biol187179192

20. Emonts M, Veenhoven RH, Wiertsema SP, Houwing-Duistermaat JJ, Walraven V, et al. (2007) Genetic polymorphisms in immunoresponse genes TNFA, IL6, IL10, and TLR4 are associated with recurrent acute otitis media. Pediatrics 120: 814–823.M. EmontsRH VeenhovenSP WiertsemaJJ Houwing-DuistermaatV. Walraven2007Genetic polymorphisms in immunoresponse genes TNFA, IL6, IL10, and TLR4 are associated with recurrent acute otitis media.Pediatrics120814823

21. Moon SK, Woo JI, Lee HY, Park R, Shimada J, et al. (2007) Toll-like receptor 2-dependent NF-kappaB activation is involved in nontypeable Haemophilus influenzae-induced monocyte chemotactic protein 1 up-regulation in the spiral ligament fibrocytes of the inner ear. Infect Immun 75: 3361–3372.SK MoonJI WooHY LeeR. ParkJ. Shimada2007Toll-like receptor 2-dependent NF-kappaB activation is involved in nontypeable Haemophilus influenzae-induced monocyte chemotactic protein 1 up-regulation in the spiral ligament fibrocytes of the inner ear.Infect Immun7533613372

22. Lee HY, Takeshita T, Shimada J, Akopyan A, Woo JI, et al. (2008) Induction of beta defensin 2 by NTHi requires TLR2 mediated MyD88 and IRAK-TRAF6-p38MAPK signaling pathway in human middle ear epithelial cells. BMC Infect Dis 8: 87.HY LeeT. TakeshitaJ. ShimadaA. AkopyanJI Woo2008Induction of beta defensin 2 by NTHi requires TLR2 mediated MyD88 and IRAK-TRAF6-p38MAPK signaling pathway in human middle ear epithelial cells.BMC Infect Dis887

23. Mason KM, Bruggeman ME, Munson RS, Bakaletz LO (2006) The non-typeable Haemophilus influenzae Sap transporter provides a mechanism of antimicrobial peptide resistance and SapD-dependent potassium acquisition. Mol Microbiol 62: 1357–1372.KM MasonME BruggemanRS MunsonLO Bakaletz2006The non-typeable Haemophilus influenzae Sap transporter provides a mechanism of antimicrobial peptide resistance and SapD-dependent potassium acquisition.Mol Microbiol6213571372

24. Tino MJ, Wright JR (1996) Surfactant protein A stimulates phagocytosis of specific pulmonary pathogens by alveolar macrophages. Am J Physiol 270: L677–688.MJ TinoJR Wright1996Surfactant protein A stimulates phagocytosis of specific pulmonary pathogens by alveolar macrophages.Am J Physiol270L677688

25. Ramet M, Lofgren J, Alho OP, Hallman M (2001) Surfactant protein-A gene locus associated with recurrent otitis media. J Pediatr 138: 266–268.M. RametJ. LofgrenOP AlhoM. Hallman2001Surfactant protein-A gene locus associated with recurrent otitis media.J Pediatr138266268

26. Lofgren J, Ramet M, Renko M, Marttila R, Hallman M (2002) Association between surfactant protein A gene locus and severe respiratory syncytial virus infection in infants. J Infect Dis 185: 283–289.J. LofgrenM. RametM. RenkoR. MarttilaM. Hallman2002Association between surfactant protein A gene locus and severe respiratory syncytial virus infection in infants.J Infect Dis185283289

27. Lahti M, Lofgren J, Marttila R, Renko M, Klaavuniemi T, et al. (2002) Surfactant protein D gene polymorphism associated with severe respiratory syncytial virus infection. Pediatr Res 51: 696–699.M. LahtiJ. LofgrenR. MarttilaM. RenkoT. Klaavuniemi2002Surfactant protein D gene polymorphism associated with severe respiratory syncytial virus infection.Pediatr Res51696699

28. LeVine AM, Whitsett JA, Gwozdz JA, Richardson TR, Fisher JH, et al. (2000) Distinct effects of surfactant protein A or D deficiency during bacterial infection on the lung. J Immunol 165: 3934–3940.AM LeVineJA WhitsettJA GwozdzTR RichardsonJH Fisher2000Distinct effects of surfactant protein A or D deficiency during bacterial infection on the lung.J Immunol16539343940

29. Takeuchi K, Maesako K, Yuta A, Sakakura Y (1994) Interleukin-8 gene expression in middle ear effusions. Ann Otol Rhinol Laryngol 103: 404–407.K. TakeuchiK. MaesakoA. YutaY. Sakakura1994Interleukin-8 gene expression in middle ear effusions.Ann Otol Rhinol Laryngol103404407

30. Hotomi M, Samukawa T, Yamanaka N (1994) Interleukin-8 in otitis media with effusion. Acta Otolaryngol 114: 406–409.M. HotomiT. SamukawaN. Yamanaka1994Interleukin-8 in otitis media with effusion.Acta Otolaryngol114406409

31. Smirnova MG, Kiselev SL, Gnuchev NV, Birchall JP, Pearson JP (2002) Role of the pro-inflammatory cytokines tumor necrosis factor-alpha, interleukin-1 beta, interleukin-6 and interleukin-8 in the pathogenesis of the otitis media with effusion. Eur Cytokine Netw 13: 161–172.MG SmirnovaSL KiselevNV GnuchevJP BirchallJP Pearson2002Role of the pro-inflammatory cytokines tumor necrosis factor-alpha, interleukin-1 beta, interleukin-6 and interleukin-8 in the pathogenesis of the otitis media with effusion.Eur Cytokine Netw13161172

32. Lee HY, Andalibi A, Webster P, Moon SK, Teufert K, et al. (2004) Antimicrobial activity of innate immune molecules against Streptococcus pneumoniae, Moraxella catarrhalis and nontypeable Haemophilus influenzae. BMC Infect Dis 4: 12.HY LeeA. AndalibiP. WebsterSK MoonK. Teufert2004Antimicrobial activity of innate immune molecules against Streptococcus pneumoniae, Moraxella catarrhalis and nontypeable Haemophilus influenzae.BMC Infect Dis412

33. Kerschner JE, Meyer TK, Burrows A (2004) Chinchilla middle ear epithelial mucin gene expression in response to inflammatory cytokines. Arch Otolaryngol Head Neck Surg 130: 1163–1167.JE KerschnerTK MeyerA. Burrows2004Chinchilla middle ear epithelial mucin gene expression in response to inflammatory cytokines.Arch Otolaryngol Head Neck Surg13011631167

34. Elsheikh MN, Mahfouz ME (2006) Up-regulation of MUC5AC and MUC5B mucin genes in nasopharyngeal respiratory mucosa and selective up-regulation of MUC5B in middle ear in pediatric otitis media with effusion. Laryngoscope 116: 365–369.MN ElsheikhME Mahfouz2006Up-regulation of MUC5AC and MUC5B mucin genes in nasopharyngeal respiratory mucosa and selective up-regulation of MUC5B in middle ear in pediatric otitis media with effusion.Laryngoscope116365369

35. Wiertsema SP, Herpers BL, Veenhoven RH, Salimans MM, Ruven HJ, et al. (2006) Functional polymorphisms in the mannan-binding lectin 2 gene: effect on MBL levels and otitis media. J Allergy Clin Immunol 117: 1344–1350.SP WiertsemaBL HerpersRH VeenhovenMM SalimansHJ Ruven2006Functional polymorphisms in the mannan-binding lectin 2 gene: effect on MBL levels and otitis media.J Allergy Clin Immunol11713441350

36. Wiertsema SP, Khoo SK, Baynam G, Veenhoven RH, Laing IA, et al. (2006) Association of CD14 promoter polymorphism with otitis media and pneumococcal vaccine responses. Clin Vaccine Immunol 13: 892–897.SP WiertsemaSK KhooG. BaynamRH VeenhovenIA Laing2006Association of CD14 promoter polymorphism with otitis media and pneumococcal vaccine responses.Clin Vaccine Immunol13892897

37. Lin J, Tsuboi Y, Pan W, Giebink GS, Adams GL, et al. (2002) Analysis by cDNA microarrays of altered gene expression in middle ears of rats following pneumococcal infection. Int J Pediatr Otorhinolaryngol 65: 203–211.J. LinY. TsuboiW. PanGS GiebinkGL Adams2002Analysis by cDNA microarrays of altered gene expression in middle ears of rats following pneumococcal infection.Int J Pediatr Otorhinolaryngol65203211

38. Parkinson N, Hardisty-Hughes RE, Tateossian H, Tsai HT, Brooker D, et al. (2006) Mutation at the Evi1 locus in Junbo mice causes susceptibility to otitis media. PLoS Genet 2: e149.N. ParkinsonRE Hardisty-HughesH. TateossianHT TsaiD. Brooker2006Mutation at the Evi1 locus in Junbo mice causes susceptibility to otitis media.PLoS Genet2e149

39. Wigginton JE, Abecasis GR (2005) PEDSTATS: descriptive statistics, graphics and quality assessment for gene mapping data. Bioinformatics 21: 3445–3447.JE WiggintonGR Abecasis2005PEDSTATS: descriptive statistics, graphics and quality assessment for gene mapping data.Bioinformatics2134453447

40. O'Connell JR, Weeks DE (1998) PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet 63: 259–266.JR O'ConnellDE Weeks1998PedCheck: a program for identification of genotype incompatibilities in linkage analysis.Am J Hum Genet63259266

41. Chen WM, Manichaikul A, Rich SS (2009) A generalized family-based association test for dichotomous traits. Am J Hum Genet 85: 364–376.WM ChenA. ManichaikulSS Rich2009A generalized family-based association test for dichotomous traits.Am J Hum Genet85364376

42. Martin ER, Monks SA, Warren LL, Kaplan NL (2000) A test for linkage and association in general pedigrees: the pedigree disequilibrium test. Am J Hum Genet 67: 146–154.ER MartinSA MonksLL WarrenNL Kaplan2000A test for linkage and association in general pedigrees: the pedigree disequilibrium test.Am J Hum Genet67146154

43. Horvath S, Xu X, Laird NM (2001) The family based association test method: strategies for studying general genotype–phenotype associations. Eur J Hum Genet 9: 301–306.S. HorvathX. XuNM Laird2001The family based association test method: strategies for studying general genotype–phenotype associations.Eur J Hum Genet9301306

44. Abecasis GR, Cherny SS, Cookson WO, Cardon LR (2002) Merlin–rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 30: 97–101.GR AbecasisSS ChernyWO CooksonLR Cardon2002Merlin–rapid analysis of dense genetic maps using sparse gene flow trees.Nat Genet3097101

45. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

46. Kubba H, Pearson JP, Birchall JP (2000) The aetiology of otitis media with effusion: a review. Clin Otolaryngol Allied Sci 25: 181–194.H. KubbaJP PearsonJP Birchall2000The aetiology of otitis media with effusion: a review.Clin Otolaryngol Allied Sci25181194

47. Kim YT, Jung HH, Ko TO, Kim SJ (2001) Up-regulation of MUC5AC mRNA expression in endotoxin-induced otitis media. Acta Otolaryngol 121: 364–370.YT KimHH JungTO KoSJ Kim2001Up-regulation of MUC5AC mRNA expression in endotoxin-induced otitis media.Acta Otolaryngol121364370

48. Takeuchi K, Yagawa M, Ishinaga H, Kishioka C, Harada T, et al. (2003) Mucin gene expression in the effusions of otitis media with effusion. Int J Pediatr Otorhinolaryngol 67: 53–58.K. TakeuchiM. YagawaH. IshinagaC. KishiokaT. Harada2003Mucin gene expression in the effusions of otitis media with effusion.Int J Pediatr Otorhinolaryngol675358

49. Schousboe LP, Rasmussen LM, Ovesen T (2001) Induction of mucin and adhesion molecules in middle ear mucosa. Acta Otolaryngol 121: 596–601.LP SchousboeLM RasmussenT. Ovesen2001Induction of mucin and adhesion molecules in middle ear mucosa.Acta Otolaryngol121596601

50. Kawano H, Paparella MM, Ho SB, Schachern PA, Morizono N, et al. (2000) Identification of MUC5B mucin gene in human middle ear with chronic otitis media. Laryngoscope 110: 668–673.H. KawanoMM PaparellaSB HoPA SchachernN. Morizono2000Identification of MUC5B mucin gene in human middle ear with chronic otitis media.Laryngoscope110668673

51. Ubell ML, Khampang P, Kerschner JE (2009) Mucin gene polymorphisms in otitis media patients. Laryngoscope. ML UbellP. KhampangJE Kerschner2009Mucin gene polymorphisms in otitis media patients.Laryngoscope

52. Preciado D, Goyal S, Rahimi M, Watson AM, Brown KJ, et al. (2010) MUC5B Is the predominant mucin glycoprotein in chronic otitis media fluid. Pediatr Res 68: 231–236.D. PreciadoS. GoyalM. RahimiAM WatsonKJ Brown2010MUC5B Is the predominant mucin glycoprotein in chronic otitis media fluid.Pediatr Res68231236

53. LeVine AM, Whitsett JA, Hartshorn KL, Crouch EC, Korfhagen TR (2001) Surfactant protein D enhances clearance of influenza A virus from the lung in vivo. J Immunol 167: 5868–5873.AM LeVineJA WhitsettKL HartshornEC CrouchTR Korfhagen2001Surfactant protein D enhances clearance of influenza A virus from the lung in vivo.J Immunol16758685873

54. Leth-Larsen R, Garred P, Jensenius H, Meschi J, Hartshorn K, et al. (2005) A common polymorphism in the SFTPD gene influences assembly, function, and concentration of surfactant protein D. J Immunol 174: 1532–1538.R. Leth-LarsenP. GarredH. JenseniusJ. MeschiK. Hartshorn2005A common polymorphism in the SFTPD gene influences assembly, function, and concentration of surfactant protein D.J Immunol17415321538

55. Heikkinen T (2000) The role of respiratory viruses in otitis media. Vaccine 19: Suppl 1S51–55.T. Heikkinen2000The role of respiratory viruses in otitis media.Vaccine19Suppl 1S5155

56. Poltorak A, He X, Smirnova I, Liu MY, Van Huffel C, et al. (1998) Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: mutations in Tlr4 gene. Science 282: 2085–2088.A. PoltorakX. HeI. SmirnovaMY LiuC. Van Huffel1998Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: mutations in Tlr4 gene.Science28220852088

57. Klein JO (1997) Role of nontypeable Haemophilus influenzae in pediatric respiratory tract infections. Pediatr Infect Dis J 16: S5–8.JO Klein1997Role of nontypeable Haemophilus influenzae in pediatric respiratory tract infections.Pediatr Infect Dis J16S58

58. Arbour NC, Lorenz E, Schutte BC, Zabner J, Kline JN, et al. (2000) TLR4 mutations are associated with endotoxin hyporesponsiveness in humans. Nat Genet 25: 187–191.NC ArbourE. LorenzBC SchutteJ. ZabnerJN Kline2000TLR4 mutations are associated with endotoxin hyporesponsiveness in humans.Nat Genet25187191