Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions

doi:10.1186/1471-2350-10-85

BMC Med Genet. 2009; 10: 85.

Published online 2009 September 3. doi: 10.1186/1471-2350-10-85

PMCID: PMC2751750

Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions

Margaretha L Casselbrant,^1,² Ellen M Mandel,^1,² Jeesun Jung,³ Robert E Ferrell,³ Kathleen Tekely,^1,² Jin P Szatkiewicz,³ Amrita Ray,³ and Daniel E Weeks^3,⁴

¹Division of Pediatric Otolaryngology, Children's Hospital of Pittsburgh of University of Pittsburgh Medical Center, Pittsburgh, PA, 15201, USA

²Department of Otolaryngology, University of Pittsburgh, School of Medicine, Pittsburgh, PA, 15213, USA

³Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, PA, 15261, USA

⁴Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, PA, 15261, USA

Corresponding author.

Margaretha L Casselbrant: margaretha.casselbrant/at/chp.edu; Ellen M Mandel: mandele/at/pitt.edu; Jeesun Jung: jeejung/at/iupui.edu; Robert E Ferrell: rferrell/at/pitt.edu; Kathleen Tekely: Kathy.Tekely/at/chp.edu; Jin P Szatkiewicz: ae4899jin/at/gmail.com; Amrita Ray: amrita.ray/at/gmail.com; Daniel E Weeks: weeks/at/pitt.edu

Received June 16, 2009; Accepted September 3, 2009.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Abstract

Background

Otitis media (OM) is a common worldwide pediatric health care problem that is known to be influenced by genetics. The objective of our study was to use linkage analysis to map possible OM susceptibility genes.

Methods

Results

In the Caucasian-only data set, our most significant linkage peak is on chromosome 17q12 at rs226088 with a p-value of 0.00007. Other peaks of potential interest are on 10q22.3 (0.00181 at rs1878001), 7q33 (0.00105 at rs958408), 6p25.1 (0.00261 at rs554653), and 4p15.2 (0.00301 at rs2133507). In the combined Caucasian and African American dataset, the 10q22.3 peak becomes more significant, with a minimal p-value of 0.00026 at rs719871. Family-based association testing reveals signals near previously implicated genes: 513 kb from SFTPA2 (10q22.3), 48 kb from IFNG (12q14), and 870 kb from TNF (6p21.3).

Conclusion

Our scan does not provide evidence for linkage in the previously reported regions of 10q26.3 and 19q13.43. Our best-supported linkage regions may contain susceptibility genes that influence the risk for recurrent/persistent OM. Plausible candidates in 17q12 include AP2B1, CCL5, and a cluster of other CCL genes, and in 10q22.3, SFTPA2.

Background

Otitis media (OM [MIM 166760]) is a worldwide pediatric health care problem. Nearly all children experience at least one episode of acute OM (AOM). Otitis media, particularly in children who experience recurrent/persistent disease, in addition to the short- and long-term effect on the child places considerable financial burden on the families (parental time, transportation, medications) and on the health-care system (clinic visits, medical and surgical therapies, complications). In 1996, Gates [1] estimated the costs for AOM and chronic otitis media with effusion (COME) in the U.S. at approximately 3 billion and 2 billion dollars, respectively. Today the annual cost can be expected to be much higher [2]. In addition to the financial burden, the psychosocial impact of OM on the child and the family is tremendous, but has not been well studied.

Epidemiologic studies using various methodologies have demonstrated evidence that genetics plays an important role in recurrent and persistent OM: population studies [3-10]; adoption study [11]; familial aggregation studies [12-16]; and twin and triplet studies. The estimated otitis media heritability in the retrospective twin study by Kvaerner et al [17] was 0.74 in females and 0.45 in males, and in the prospective twin and triplet study of children followed from birth to 2 years of age by Casselbrant et al [18] the heritability of time with middle ear effusion was 0.79 in females and 0.64 in males. Similarly, a prospective study of a community sample of twins estimated heritability of total OM symptom scores as 0.49 to 0.71, depending on age [19].

While there is substantial evidence that genetics plays a role in OM, there has been only one genome-wide linkage scan to date by Daly et al [20] that provided evidence of linkage of COME and recurrent OM (ROM) to 10q26.3 and to 19q13.43. Subsequent multipoint linkage analysis of both regions further strengthened the evidence of linkage [21].

Otitis media is a complex disease with contributions from immunocompetence, inflammatory regulation and craniofacial abnormalities that could lead to Eustachian tube dysfunction; therefore, the underlying genetic determinants are likely to be complex and involve several loci. Ilia et al [22] in a recent review describe how case-control studies and animal studies have implicated genes involved in non-specific immunity (epithelial barriers and mucins), innate immunity (toll-like receptors, mannose-binding lectin, surfactant proteins, cytokines), and adaptive immunity (Fc-γ receptors, immunoglobins).

Our present study aims to localize genes that contribute to susceptibility to recurrent/persistent middle ear disease by carrying out a genome-wide screen for linkage using families containing at least two affected siblings. The identification of susceptibility genes may enable us to better understand the pathogenesis of OM leading to development of better and more innovative methods for prevention and treatment, which is especially important in this era of multi-drug resistant bacteria.

Methods

Population, Enrollment and Assessment

Full siblings, two or more, who had a history of tympanostomy tube insertion due to a significant history of OM, their parent(s) and other full sibling(s) with no history of tympanostomy tube insertion were eligible for the study. There was no exclusion of subjects based on gender or race. The study was explained in detail to the parent(s) and children, and informed consent was obtained prior to enrollment. This study was approved by the Human Rights Committee of the Children's Hospital of Pittsburgh (CHP) and by the University of Pittsburgh Institutional Review Board.

The study required only one visit to the Ear, Nose and Throat (ENT) Research Center at CHP. A detailed history regarding recurrent/persistent OM was obtained for each enrolled family member, as well as history regarding risk factors such as breast-feeding, day care attendance, siblings, and exposure to smoking for the enrolled children. Medical records of enrolled children were obtained for review whenever possible. When feasible, children and parent(s) had an ear examination using pneumatic otoscopy by the study physicians (MC and EM) who are validated otoscopists [23]. Tympanograms were obtained using a GSI-38 middle ear analyzer (Lucas-Grason-Stadler, INC). Children were excluded for the following reasons: major congenital malformations, medical conditions with a predisposition for OM (e.g. cleft palate, Down syndrome, or other craniofacial malformations), cared for in the Intensive Care Unit as a neonate; been on assisted ventilation, or known sensorineural hearing loss. The subjects were selected from patients who presented to the ENT Research Center from several sources: the ENT Clinic at the CHP and satellite clinics, subjects who were in other studies at the ENT Research Center, or from physician- or self-referral.

Phenotype definition

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. While still recording each subject's medical history, 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.

DNA isolation and genotyping

EDTA anticoagulated whole blood (5 cc if less than 5 years old; 5-10 cc if 5 years or older; 10-20 cc from parents) was collected from an arm vein in the Pediatric General Clinical Research Center or in the ENT Research Center or in satellite clinics and transferred to the Human Genetics Laboratory within 48 hours. High molecular weight DNA was isolated by the method of Miller et al [24]. Following re-suspension, the DNA was quantitated by fluorometry, diluted to 50 ng/μl, arrayed in a standard 96 well format and stored at -80°C.

Genotyping was carried out using Affymetrix's GeneChip^®Human Mapping 10K Array in two different versions: 1) the original GeneChip^®Human Mapping 10K Array Xba 131 with 11,560 SNPs with a mean intermarker distance of 210 kb and an average heterozygosity of 0.37, and 2) the newer GeneChip^®Human Mapping 10K Array Xba 142 2.0 with 10,204 SNPs with a mean intermarker distance of 258 kb and an average heterozygosity of 0.38. These two chips are very similar to each other: the newer SNP chip was created by removing 1,400 SNPs from the original SNP list and adding approximately 70 new SNPs. Samples were processed according to the GeneChip Mapping Assay Manual (Affymetrix) hybridized at 48°C for 16 hours in an Affymetrix GeneChip Hybridization Oven. After 16 hours the probe arrays were washed and stained according to the GeneChip Mapping Assay Manual using the DNAARRAY_WS2 protocol on the Affymetrix Fluidics Station 400. Arrays were scanned once with the Agilent GeneArray 2500 scanner and analyzed with Affymetrix GeneChip DNA Analysis Software (GDAS) to generate genotype calls for each of the SNP probes on the array. In this study we used the older Xba 131 SNP chip to genotype 515 people, and the newer Xba 142 2.0 SNP chip to genotype 1,216 people.

Genetic maps

Starting with the available annotation from Affymetrix, we carefully updated it, as a fair number of SNPs are actually bad, hybridizing to multiple positions in the genome. To detect these, we took the primer sequence for each SNP and used BLAST to find its position(s) in the genome, carefully identifying poor multi-hit SNPs, as well as verifying the current best physical position for each good SNP in the data set. After careful resolution of reference sequence (rs) numbers, current physical map positions were retrieved from Ensembl and used to linearly interpolate genetic map positions based on the Rutgers Combined Linkage-Physical Map [25].

Quality Control

To verify the stated pedigree structures, we carried out relationship testing using a subset of 1,534 well-typed informative autosomal SNPs (with MAF > = 0.40 and > = 99% genotyping rate) and then made a number of changes to the structures to minimize the number of relationship errors while keeping as many people as possible.

PREST identified 73 within-family pairs with p-values < 0.001 [26,27]. Similarly, RELPAIR examined 1,287,210 genotyped relative pairs, and identified 70 within-family pairs whose inferred relationship was > 1,000 times more likely than the putative relationship [28,29]. Checking for relationships between putative families with RELPAIR revealed two families that had entered the study twice.

Gender checking based on PLINK's algorithm led to the removal of genotypes for 6 people because of unresolved gender errors [30]. We replaced the stated gender with the inferred gender in 4 cases where such a change was symmetric in terms of relationship and OM affection status, and so would not alter the linkage results.

A number of conservative changes were made in response to the relationship testing results, including the following: To resolve the 10 putative monozygotic (MZ) twin pairs identified, from the larger families, we removed the least typed member of an MZ pair 6 times, and we completely removed 2 families that each contained only a single inferred MZ pair. To resolve putative half-sib (HS) families, we removed 6 families containing only a single inferred HS pair, and from larger families, we removed a single person to eliminate the putative HS pair 4 times. To resolve the two families that had participated twice in our study, in both cases, only the more completely typed family was kept.

After resolving the relationship errors, we then removed all genotypes for the 250 individuals with a per person genotyping success rate less than 90%.

Identity-by-state cluster analysis based on PLINK's hierarchical clustering algorithm was used to verify that the self-reported ethnicities were consistent with the clusters generated from the SNP marker data.

Consistent with the ethnicity distribution in our sampling area, the majority of our families, 403, were Caucasian, while there were only 26 African American families. For the purposes of analyses, the larger Caucasian data set was first analyzed by itself. Here, 147 SNPs were excluded due to a HWE p-value < = 0.001 (computed by PLINK using only founder genotypes); 1,936 SNPs were excluded because their genotyping success rate across all individuals was less than 90%; 999 SNPs were excluded due to an MAF < 0.05; leaving 8,802 SNPs that were used in the linkage analyses.

Allele frequencies

For the initial analyses using Merlin [31], allele frequencies were estimated from the data, ignoring family structure, by simple gene counting. On small nuclear families such as we have here, this approach results in unbiased estimates [32]. For the subsequent analyses using Mendel [33], allele frequencies were estimated by maximizing the pedigree likelihood while properly taking the family structures into account.

Linkage Analysis

We used PedCheck [34] to detect Mendelian inconsistencies; all genotypes at each problematic SNP were zeroed out within each family containing a Mendelian inconsistency. We used Merlin [31] to carry out non-parametric linkage analyses using the S_allstatistic using sex-averaged genetic maps and allele-counting allele frequency estimates; Minx was used to analyze the X-linked data. Summary statistics were computed using Pedstats [35]. File-format conversions were done using Mega2 [36].

To assay the effects of SNP-SNP linkage disequilibrium, we recomputed the LOD scores using Merlin's LD modeling (--rsq 0.1) option [37]. To examine robustness to other modeling assumptions, we also analyzed the data using the Mendel v.9.0.0 package [33] in a series of three steps: (a) We used Mendel's mistyping option to zero out all genotypes at each problematic SNP within each family containing a Mendelian inconsistency; (b) We used Mendel's allele frequency option to estimate the allele frequencies at each SNP while properly taking the family structures into account; (c) We used Mendel's non-parametric linkage (NPL) option to compute an 'additive all' linkage statistic, using the replicate pool method with 50 replicates and 100,000 samples to generate accurate empirical P-values. These analyses used sex-specific genetic maps for the autosomes, and the female map for the X chromosome.

Combined linkage analyses

After the initial analyses using the Caucasian-only data set, we carried out additional linkage analyses of the combined Caucasian and African American data set. For these analyses, we used Mendel's options for estimating and using ethnic-group specific allele frequencies. We did not analyze the African American families by themselves, as the sample size was too small by itself to have adequate power: only 26 African American families containing 27 genotyped affected sib pairs.

Candidate gene prediction

We used the GRAIL "Gene Relationships Across Implicated Loci" software to search our linkage peaks for possible candidate genes based on similarity to a set of previously implicated OM candidate genes [38]. Similarity is measured using automated statistical text mining of 250,000 PubMed abstracts. GRAIL can identify relatedness between genes even when there are no published co-citations or established pathways in common. GRAIL was run using a seed list of previously implicated candidate genes. While the GRAIL results do depend on the seed list, the GRAIL approach does avoid circular reasoning when assigning a score to each query region by not considering any seed regions that share genes with the query region.

Association analyses

Using the Caucasian-only data set, we tested for association of 8,585 autosomal SNPs (which met our quality control criteria) using the gamete competition statistic implemented in Mendel [33]. This statistic, which tests the null hypothesis of no linkage and no association, uses all the family information to test for over-transmission of a specific SNP allele to affecteds with complementary under-transmission to unaffecteds. These results were annotated using the WGA Viewer program [39].

Results

Using a stringent diagnostic model in which only those who underwent tympanostomy tube insertion at least once for recurrent/persistent OM are considered affected, we have carried out a genome-wide linkage scan using the 10K Affymetrix SNP panel. We genotyped a panel of 403 Caucasian families containing 1,431 genotyped individuals, and 26 African American families containing 75 genotyped individuals. For the purposes of analyses, the larger Caucasian data set was first analyzed by itself. After careful quality control, non-parametric linkage analysis was carried out using 8,802 SNPs; the cleaned Caucasian data contained 377 genotyped affected sib pairs (ASPs) (89 female-female ASPs, 110 male-male ASPs, and 178 male-female ASPs). Initial linkage analyses with Merlin found two regions with linkage peaks with linear S_allLOD scores > 2.0 with a maximum LOD of (a) 2.83 on chromosome 17q12 at rs226088 and (b) 2.25 on 6p25.1 at rs554653 (Additional file 1). When marker-to-marker linkage disequilibrium (LD) is modeled, we obtain a maximum LOD of 2.85 at rs226088, and of 2.16 at rs554653. Additionally, we observed a maximum LOD of 2.00 on chromosome 7q33 at rs1343697 (increased from 1.96 previously). Sensitivity to marker-to-marker LD should be modest, as 57.6% of the families have both founders genotyped, and 38.5% have one founder genotyped.

Mendel-based linkage results, which use better allele-frequency estimates, sex-specific maps, and empirical P-values, are summarized in Table Table11 and Figure Figure1.1. For the Caucasian-only data set, the most significant empirical p-value on chromosome 17q12 is 0.00007 at rs226088. There are four other peaks with suggestive p-values: chromosome 10q22.3 with a minimal p-value of 0.00181 at rs1878001; chromosome 7q33 with a minimal p-value of 0.00105 at rs958408; chromosome 6p25.1 with a minimal empirical p-value of 0.00261 at rs554653; and chromosome 4p15.2 with a minimal p-value of 0.00301 at rs2133507.

Table 1

Additive ALL linkage statistics as computed by Mendel, grouped by peak, for the Caucasian-only data set and for the combined Caucasian and African American data set.

Figure 1

The empirical -log10(p-values) for the Additive ALL statistics for the five highest peaks. Solid line: Caucasian-only results; dashed line: combined Caucasian and African American results. The horizontal reference line indicates a nominal P-value of 0.003. (more ...)

When we carried out analyses of the combined Caucasian and African American sample, using group-specific allele frequency estimates, three of the peaks (on chromosomes 17, 6, 4) became a bit less significant, while the chromosome 7 peak significance level stayed essentially the same and the chromosome 10 peak became more significant, with a minimal p-value of 0.00026 at rs719871 (Table (Table11 and Figure Figure1).1). As mentioned above, the combined sample adds 26 African American families containing 75 genotyped individuals. This African American sample contains 27 genotyped affected sib pairs (ASPs) (5 female-female ASPs, 8 male-male ASPs, and 14 male-female ASPs).

We used the GRAIL "Gene Relationships Across Implicated Loci" software to search our linkage peaks for possible candidate genes based on similarity to a 'seed' set of previously implicated OM candidate genes [38]. To do this, we used a seed list of 36 candidate genes (ADAMTS13, BCL6, CD14, DNAH5, EVI1, EYA4, FBXO11, FCGR2A, GJB2, GJB6, HLA-A, IFNG, IL10, IL13, IL1A, IL1RN, IL4R, IL6, MAPK14, MAPK8, MBL2, MMP2, MMP9, MUC2, MUC4, MUC5AC, MUC5B, SERPINE1, SERPING1, SFTPA1B, SFTPA2B, SFTPD, TGFB1, TLR2, TLR4, TNF). For the chromosome 17 linkage peak, the most significant candidate genes with p-values < 10^-13are CCL5, CCL18, CCL3, and CCL4. For the chromosome 10 peak, the most significant candidate genes with p-values < 10^-13are SFTPA1B, SFTPA2, SFTPD, and PLAU. For the chromosome 7 peak, the most significant candidate gene EPHB6 had a p-value of 2.6 × 10^-6. For the chromosome 6 peak, the most significant candidate genes with p-values < 10^-11are LY86, SERPINB9, and F13A1. For the chromosome 4 peak, the most significant candidate gene with p-values < 10^-11was SOD3.

Using the Caucasian-only data set, we tested for association of 8,585 autosomal SNPs (which met our quality control criteria) using the gamete competition statistic implemented in Mendel [33]. The most significant association signals are listed in Additional file 2. While none of these are significant after Bonferroni or false discovery rate correction for multiple testing, there are a few highly ranked signals that are near previously implicated candidate genes. One SNP (rs1437803, P-value 0.0005, rank 5) is 513 kb from SFTPA2 on chromosome 10. Two SNPs (rs722749, P-value 0.0009, rank 11 and rs722748, P-value 0.0013, rank 15) are 48 kb from IFNG on chromosome 12. Two SNPs (rs2213584, P-value 0.0021, rank 22 and rs2227139, P-value 0.0021, rank 23) are 870 kb from TNF on chromosome 6. Note that the gamete competition test should detect association signals at much larger distances than a regular case/control association test in unrelateds because the gamete competition test relies on measuring transmission distortion within families.

Discussion

While there is fairly strong evidence that otitis media has a strong genetic component, to our knowledge there has been only one prior family-based linkage study of otitis media. This first study provided evidence of linkage of COME/ROM to chromosome 10q26.3 at marker D10S212 and to chromosome 19q13.43 at marker D19S254 [20]. We report here on our genome-wide linkage scan for risk genes for otitis media. In our study, we chose to focus on severe disease, collecting affected sib pair families where at least two full siblings had undergone tympanostomy tube insertion. By narrowing the disease definition, we hoped to identify a more severe form of OM that is perhaps due to a small number of loci. This approach has proven remarkably successful in several complex diseases, such as breast cancer and Alzheimer's disease, where subdivision according to age of onset contributed to the mapping and cloning disease genes [40,41]. Additionally, it has been shown that there are statistical advantages in using a "narrow-phenotype" approach [42], because the power of an affected sib pair approach increases as the population prevalence of the trait decreases.

The rate of tympanostomy tube insertion in second siblings may be increased due to the parents' wishes to have the procedure performed earlier or for less significant disease than the first child because of prior positive experiences in the first sibling. This could potentially make it difficult to separate the genetic effect of OM from the parents' influence on the treatment. However, we have reason to believe such an effect, if it exists, will be small in our sample. Firstly, the majority of our sample was treated at the Children's Hospital of Pittsburgh, where our criteria for insertion for tubes are consistent and very stringent. Secondly, at entry we obtained, for each subject, a history of middle-ear disease prior to tube insertion to assess the subject's eligibility.

Since the majority of our families were Caucasian, we present two linkage analyses, one where we analyzed only the 403 Caucasian families by themselves, and a combined analysis where we jointly analyzed both the 403 Caucasian and 26 African American families (Table (Table11 and Figure Figure1).1). The Caucasian-only analyses generate a strong linkage signal on chromosome 17q12, as well as four other potentially interesting peaks (10q22.3, 7q33, 6p25.1, 4p15.2). The combined analyses strengthen the evidence for the 10q22.3 peak. Our scan does not provide evidence for linkage in the previously reported regions of 10q26.3 and 19q13.43 [20,21].

The very tip of our 17q12 linkage peak occurs in AP2B1 (adaptor-related protein complex 2, beta 1 subunit), which plays a role in Nef-mediated CD8 down-regulation [43], and children with recurrent otitis media had low numbers of CD8+-producing IFN g cells in adenoids [44]. However, the chromosome 17 peak also contains a cluster of CCL (chemokine C-C motif ligand) genes, several of which were highlighted as possible candidates by the GRAIL analyses. CCL5, also known as RANTES, is 18 kb from the linkage peak, and has been previously associated with otitis media [45-52]. CCL5 is an eosinophil chemoattactrant that is thought to play a role in the accumulation of eosinophils often observed in middle ear effusions of OM with allergy.

The combined 10q22.3 linkage peak is 1.2 Mb from a previously implicated candidate gene SFTPA2. This peak also contains a strongly associated SNP rs1437803 (Additional file 2, P-value 0.0005, rank 5), which is 513 kb from SFTPA2. The human surfactant protein A (SP-A) is expressed in the Eustachian tube, plays a role in innate host defense, upregulates phagocytosis of many OM risk pathogens (including Streptococus pneumoniae, Haemophilus influenzae, and respiratory syncytial virus), and consists of two very similar functional genes SFTPA1 and SFTPA2 located 5 kb apart on chromosome 10. Ramet et al [53] reported that the frequency of specific SP-A haplotypes and genotypes differ between children who experience their first episode of acute otitis before age 6 months and the general population. Pettigrew et al [54] also found that polymorphisms within the SP-A loci were protective for otitis media among white infants at risk for asthma.

In the remaining linkage peaks, the majority of the genes highlighted by GRAIL as possible candidates do not have prior evidence of involvement with OM.

While none of the association results (Additional file 2) are significant after correction for multiple testing, the level of evidence required is less for previously implicated genes. Thus, in addition to the associated SNP 513 kb from SFTPA2 mentioned above, it is noteworthy that we have associated SNPs 48 kb from interferon-γ (IFNG) on chromosome 12. Genetic variants at IFNG are associated with risk for otitis media in infants infected with the respiratory syncytial virus (RSV) [55]. When RSV-induced and non-RSV-induced otitis media cases were compared, significantly higher levels of IFNG were found in the RSV-induced cases [56]. We also observed two associated SNPs 870 kb from tumor necrosis factor (TNF) on chromosome 6. A TNF -308 polymorphism was associated with both otitis media susceptibility and placement of tympanostomy tubes [57]. Variants in the promoter region of TNF were associated with being otitis-prone [58].

The gamete competition test of association used here relies on transmission distortion within families, and so should detect association signals at larger distances than the conventional case/control association test in unrelateds. Even so, the proximity of some association signals to previously implicated candidate genes mentioned above may not be that meaningful, as these regions contain many other plausible candidate genes.

Conclusion

Our linkage scan, the largest to date, has identified two strong linkage peaks, on 17q12 and 10q22.3 that are likely to contain genes which influence risk for severe OM. While both of these peaks contain intriguing and plausible candidate genes, further fine-mapping, replication, and functional studies are required to reach more firm conclusions. We also find association signals near previously implicated risk genes (SFTPA2, IFNG, and TNF). Our hope is that our genetic results will help contribute to an enhanced understanding of the etiology of otitis media that ultimately may help lead to improved treatment and prevention of this extremely common childhood disease.

Availability

The URLs for data presented herein are as follows:

PREST, http://fisher.utstat.toronto.edu/sun/Software/Prest/

RELPAIR, http://csg.sph.umich.edu/boehnke/relpair.php

PLINK, http://pngu.mgh.harvard.edu/~purcell/plink/

Ensembl, http://www.ensembl.org/index.html

Rutgers Combined Linkage-Physical Map, http://compgen.rutgers.edu/maps/

PedCheck, http://watson.hgen.pitt.edu/register/

Pedstats, http://www.sph.umich.edu/csg/abecasis/PedStats/

Merlin, http://www.sph.umich.edu/csg/abecasis/merlin/

Mega2, http://watson.hgen.pitt.edu/register/

Mendel, http://www.genetics.ucla.edu/software/

Online Mendelian Inheritance in Man (OMIM), http://www.ncbi.nlm.nih.gov/sites/entrez?db=omim

WGA Viewer, http://www.genome.duke.edu/centers/pg2/downloads/wgaviewer.php

GRAIL, http://www.broad.mit.edu/mpg/grail/

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MLC, EMM, REF, and DEW conceived of the study and obtained funding; KT, EMM, and MLC carried out the clinical work; REF supervised the genotyping; JPS, JJ, AR, and DEW carried out the statistical analyses; DEW, MLC, REF, and EMM wrote the initial draft of the manuscript; all authors read and approved the manuscript.

Authors' information

Current affiliations: JJ: Department of Medical and Molecular Genetics, Center for Computational Biology and Bioinformatics, Indiana University, School of Medicine, Indianapolis, IN, 46202, USA; JPS: Eating Disorders Program, Department of Psychiatry, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA; AR: Statistics, reporting, and analysis (Cardiac Rhythm Management division), St. Jude Medical, Sunnyvale, CA, 94086, USA.

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2350/10/85/prepub

Supplementary Material

Additional file 1

Supplemental Figure S1. This file provides multipoint linkage curves as computed by Merlin with and without linkage disequilibrium.

Click here for file^{(5.7M, pdf)}

Additional file 2

Supplemental Table S1. This file provides the gamete competition results in the Caucasian families for those SNPs with p-values < 0.01.

Click here for file^{(103K, pdf)}

Acknowledgements

We would like to thank all of our families for participating in our research study, and James Seroky, MA, for help in data management. We would like to thank the University of Pittsburgh Genomics and Proteomics Core Laboratories for carrying out the genotyping. Our research was supported primarily by the National Institutes of Health grant R01DC005630 [to M.L.C.], as well as by the National Institutes of Health General Clinical Research Center grant M01RR000084, National Institutes of Health Clinical and Translational Science Award UL1RR024153, Children's Hospital of Pittsburgh of UPMC, and the University of Pittsburgh. Amrita Ray was supported by the Fogarty International Center (FIC) grant D43TW006180 "India-US Research Training Program in Genetics" [to D.E.W.]. Jin Szatkiewicz was supported by the National Institutes of Health National Institute of Mental Health training grant "Discovering Genes for Mental Health" [T32MH020053 to D.E.W.]. Jeesun Jung was supported by the National Institutes of Health grant R01MH064205 [to D.E.W.].

References

Gates GA. Cost-effectiveness considerations in otitis media treatment. Otolaryngol Head Neck Surg. 1996;114:525–530. doi: 10.1016/S0194-5998(96)70243-7. [PubMed] [Cross Ref]
Rovers MM. The burden of otitis media. Vaccine. 2008;26:G2–4. doi: 10.1016/j.vaccine.2008.11.005. [PubMed] [Cross Ref]
Brody JA, Overfield T, McAlister R. Draining ears and deafness among Alaskan Eskimos. Arch Otolaryngol. 1965;81:29–33. [PubMed]
Maynard JE. Otitis media in Alaskan Eskimo children: an epidemiologic review with observations on control. Alaska Med. 1969;11:93–98. [PubMed]
Wiet RJ. Patterns of ear disease in the southwestern American Indian. Arch Otolaryngol. 1979;105:381–385. [PubMed]
Moran DJ, Waterford JE, Hollows F, Jones DL. Ear disease in rural Australia. Med J Aust. 1979;2:210–212. [PubMed]
Beery QC, Doyle WJ, Cantekin EI, Bluestone CD, Wiet RJ. Eustachian tube function in an American Indian population. Ann Otol Rhinol Laryngol Suppl. 1980;89:28–33. [PubMed]
Casselbrant ML, Mandel EM, Kurs-Lasky M, Rockette HE, Bluestone CD. Otitis media in a population of black American and white American infants, 0-2 years of age. Int J Pediatr Otorhinolaryngol. 1995;33:1–16. doi: 10.1016/0165-5876(95)01184-D. [PubMed] [Cross Ref]
Paradise JL, Rockette HE, Colborn DK, Bernard BS, Smith CG, Kurs-Lasky M, Janosky JE. Otitis media in 2253 Pittsburgh-area infants: prevalence and risk factors during the first two years of life. Pediatrics. 1997;99:318–333. doi: 10.1542/peds.99.3.318. [PubMed] [Cross Ref]
Doyle WJ. dissertation. University of Pittsburgh; 1977. A functional-anatomic description of eustachian tube vector relations in four ethinic populations: an osteologic study.
Spivey GH, Hirschhorn N. A migrant study of adopted Apache children. Johns Hopkins Med J. 1977;140:43–46. [PubMed]
Daly KA, Rich SS, Levine S, Margolis RH, Le CT, Lindgren B, Giebink GS. The family study of otitis media: design and disease and risk factor profiles. Genet Epidemiol. 1996;13:451–468. doi: 10.1002/(SICI)1098-2272(1996)13:5<451::AID-GEPI2>3.0.CO;2-5. [PubMed] [Cross Ref]
Sipila M, Karma P, Pukander J, Timonen M, Kataja M. The Bayesian approach to the evaluation of risk factors in acute and recurrent acute otitis media. Acta Otolaryngol. 1988;106:94–101. doi: 10.3109/00016488809107375. [PubMed] [Cross Ref]
Todd NW. Familial predisposition for otitis media in Apache Indians at Canyon Day, Arizona. Genet Epidemiol. 1987;4:25–31. doi: 10.1002/gepi.1370040104. [PubMed] [Cross Ref]
Stenstrom C, Ingvarsson L. Otitis-prone children and controls: a study of possible predisposing factors. 1. Heredity, family background and perinatal period. Acta Otolaryngol. 1997;117:87–93. doi: 10.3109/00016489709117997. [PubMed] [Cross Ref]
Rovers MM, Zielhuis GA, Straatman H, Ingels K, Wilt GJ van der, Broek P van den. Prognostic factors for persistent otitis media with effusion in infants. Arch Otolaryngol Head Neck Surg. 1999;125:1203–1207. [PubMed]
Kvaerner KJ, Tambs K, Harris JR, Magnus P. Distribution and heritability of recurrent ear infections. Ann Otol Rhinol Laryngol. 1997;106:624–632. [PubMed]
Casselbrant ML, Mandel EM, Fall PA, Rockette HE, Kurs-Lasky M, Bluestone CD, Ferrell RE. The heritability of otitis media: a twin and triplet study. JAMA. 1999;282:2125–2130. doi: 10.1001/jama.282.22.2125. [PubMed] [Cross Ref]
Rovers M, Haggard M, Gannon M, Koeppen-Schomerus G, Plomin R. Heritability of symptom domains in otitis media: a longitudinal study of 1,373 twin pairs. Am J Epidemiol. 2002;155:958–964. doi: 10.1093/aje/155.10.958. [PubMed] [Cross Ref]
Daly KA, Brown WM, Segade F, Bowden DW, Keats BJ, Lindgren BR, Levine SC, Rich SS. Chronic and recurrent otitis media: a genome scan for susceptibility loci. Am J Hum Genet. 2004;75:988–997. doi: 10.1086/426061. [PubMed] [Cross Ref]
Sale M, Perlegas P, Marion M, Hicks P, Rich S, Daly K. Fine mapping of two loci influencing chronic otitis media with effusion and/or recurrent otitis media (COME/ROM) 9th International Symposium on Recent Advances in Otitis Media: 2007; St Pete Beach, Florida: 9th International Symposium on Recent Advances in Otitis Media. 2007. p. 87.
Ilia S, Goulielmos GN, Samonis G, Galanakis E. Host's response in otitis media: understanding genetic susceptibility. Pediatr Infect Dis J. 2008;27:929–933. doi: 10.1097/INF.0b013e318174e0b0. [PubMed] [Cross Ref]
Stool SE, Flaherty MR. Validation of diagnosis of otitis media with effusion. Ann Otol Rhinol Laryngol. 1983;92:5–6.
Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16:1215. doi: 10.1093/nar/16.3.1215. [PMC free article] [PubMed] [Cross Ref]
Matise TC, Chen F, Chen W, De La Vega FM, Hansen M, He C, Hyland FC, Kennedy GC, Kong X, Murray SS, et al. A second-generation combined linkage physical map of the human genome. Genome Res. 2007;17:1783–1786. doi: 10.1101/gr.7156307. [PubMed] [Cross Ref]
McPeek MS, Sun L. Statistical tests for detection of misspecified relationships by use of genome-screen data. Am J Hum Genet. 2000;66:1076–1094. doi: 10.1086/302800. [PubMed] [Cross Ref]
Sun L, Wilder K, McPeek MS. Enhanced pedigree error detection. Hum Hered. 2002;54:99–110. doi: 10.1159/000067666. [PubMed] [Cross Ref]
Boehnke M, Cox NJ. Accurate inference of relationships in sib-pair linkage studies. Am J Hum Genet. 1997;61:423–429. doi: 10.1086/514862. [PubMed] [Cross Ref]
Epstein MP, Duren WL, Boehnke M. Improved inference of relationship for pairs of individuals. Am J Hum Genet. 2000;67:1219–1231. [PubMed]
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. doi: 10.1086/519795. [PubMed] [Cross Ref]
Abecasis GR, Cherny SS, Cookson WO, Cardon LR. Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002;30:97–101. doi: 10.1038/ng786. [PubMed] [Cross Ref]
Broman KW. Estimation of allele frequencies with data on sibships. Genet Epidemiol. 2001;20:307–315. doi: 10.1002/gepi.2. [PubMed] [Cross Ref]
Lange K, Cantor R, Horvath S, Perola M, Sabatti C, Sinsheimer J, Sobel E. MENDEL version 4.0: A complete package for the exact genetic analysis of discrete traits in pedigree and population data sets. Am J Hum Genet. 2001;69:504. [PubMed]
O'Connell JR, Weeks DE. PedCheck: A program for identifying genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998;63:259–266. doi: 10.1086/301904. [PubMed] [Cross Ref]
Wigginton JE, Abecasis GR. PEDSTATS: descriptive statistics, graphics and quality assessment for gene mapping data. Bioinformatics. 2005;21:3445–3447. doi: 10.1093/bioinformatics/bti529. [PubMed] [Cross Ref]
Mukhopadhyay N, Almasy L, Schroeder M, Mulvihill WP, Weeks DE. Mega2: data-handling for facilitating genetic linkage and association analyses. Bioinformatics. 2005;21:2556–2557. doi: 10.1093/bioinformatics/bti364. [PubMed] [Cross Ref]
Abecasis GR, Wigginton JE. Handling marker-marker linkage disequilibrium: pedigree analysis with clustered markers. Am J Hum Genet. 2005;77:754–767. doi: 10.1086/497345. [PubMed] [Cross Ref]
Raychaudhuri S, Plenge RM, Rossin EJ, Ng AC, Purcell SM, International Schizophrenia Consortium. Sklar P, Scolnick EM, Xavier RJ, Altshuler D, et al. Identifying relationships among genomic disease regions: Predicting genes at pathogenic SNP associations and rare deletions. PLoS Genet. 2009;5:e1000534. doi: 10.1371/journal.pgen.1000534. [PMC free article] [PubMed] [Cross Ref]
Ge D, Zhang K, Need AC, Martin O, Fellay J, Urban TJ, Telenti A, Goldstein DB. WGAViewer: software for genomic annotation of whole genome association studies. Genome Res. 2008;18:640–643. doi: 10.1101/gr.071571.107. [PubMed] [Cross Ref]
Nathanson KN, Wooster R, Weber BL. Breast cancer genetics: what we know and what we need. Nat Med. 2001;7:552–556. doi: 10.1038/87876. [PubMed] [Cross Ref]
Schellenberg GD, D'Souza I, Poorkaj P. The genetics of Alzheimer's disease. Curr Psychiatry Rep. 2000;2:158–164. doi: 10.1007/s11920-000-0061-z. [PubMed] [Cross Ref]
Rice JP, Saccone NL, Rasmussen E. Definition of the phenotype. Adv Genet. 2001;42:69–76. full_text. [PubMed]
Stove V, Walle I Van de, Naessens E, Coene E, Stove C, Plum J, Verhasselt B. Human immunodeficiency virus Nef induces rapid internalization of the T-cell coreceptor CD8alphabeta. J Virol. 2005;79:11422–11433. doi: 10.1128/JVI.79.17.11422-11433.2005. [PMC free article] [PubMed] [Cross Ref]
Avanzini AM, Castellazzi AM, Marconi M, Valsecchi C, Marseglia A, Ciprandi G, De Silvestri A, Marseglia GL. Children with recurrent otitis show defective IFN gamma-producing cells in adenoids. Pediatr Allergy Immunol. 2008;19:523–526. doi: 10.1111/j.1399-3038.2007.00682.x. [PubMed] [Cross Ref]
Iino Y, Kakizaki K, Katano H, Saigusa H, Kanegasaki S. Eosinophil chemoattractants in the middle ear of patients with eosinophilic otitis media. Clin Exp Allergy. 2005;35:1370–1376. doi: 10.1111/j.1365-2222.2005.02330.x. [PubMed] [Cross Ref]
Das S, Palmer OP, Leight WD, Surowitz JB, Pickles RJ, Randell SH, Buchman CA. Cytokine amplification by respiratory syncytial virus infection in human nasal epithelial cells. Laryngoscope. 2005;115:764–768. doi: 10.1097/01.MLG.0000159527.76949.93. [PubMed] [Cross Ref]
Liang G, Cui Y, Wei M, Liu Z, Wang C. [The expression of RANTES interleukin 1beta and endotoxin in the effusion of secretory otitis media] Lin Chuang Er Bi Yan Hou Ke Za Zhi. 2004;18:411–413. [PubMed]
Kariya S, Okano M, Aoji K, Kosaka M, Chikumoto E, Hattori H, Yuen K, Nishioka S, Nishioka K, Nishizaki K. Role of macrophage migration inhibitory factor in otitis media with effusion in adults. Clin Diagn Lab Immunol. 2003;10:417–422. [PMC free article] [PubMed]
Jang CH, Kim YH. Demonstration of RANTES and eosinophilic cataionic protein in otitis media with effusion with allergy. Int J Pediatr Otorhinolaryngol. 2003;67:531–533. doi: 10.1016/S0165-5876(03)00015-6. [PubMed] [Cross Ref]
Schousboe LP, Rasmussen LM, Ovesen T. Induction of mucin and adhesion molecules in middle ear mucosa. Acta Otolaryngol. 2001;121:596–601. doi: 10.1080/000164801316878881. [PubMed] [Cross Ref]
Schousboe LP, Rasmussen LM, Ovesen T. RANTES in otitis media with effusion: presence, role and correlation with cytokines and microbiology. APMIS. 2001;109:441–446. doi: 10.1034/j.1600-0463.2001.090606.x. [PubMed] [Cross Ref]
Moyse E, Lyon M, Cordier G, Mornex JF, Collet L, Froehlich P. Viral RNA in middle ear mucosa and exudates in patients with chronic otitis media with effusion. Arch Otolaryngol Head Neck Surg. 2000;126:1105–1110. [PubMed]
Ramet M, Lofgren J, Alho OP, Hallman M. Surfactant protein-A gene locus associated with recurrent otitis media. J Pediatr. 2001;138:266–268. doi: 10.1067/mpd.2001.110133. [PubMed] [Cross Ref]
Pettigrew MM, Gent JF, Zhu Y, Triche EW, Belanger KD, Holford TR, Bracken MB, Leaderer BP. Association of surfactant protein A polymorphisms with otitis media in infants at risk for asthma. BMC Med Genet. 2006;7:68. doi: 10.1186/1471-2350-7-68. [PMC free article] [PubMed] [Cross Ref]
Gentile DA, Doyle WJ, Zeevi A, Howe-Adams J, Kapadia S, Trecki J, Skoner DP. Cytokine gene polymorphisms moderate illness severity in infants with respiratory syncytial virus infection. Hum Immunol. 2003;64:338–344. doi: 10.1016/S0198-8859(02)00827-3. [PubMed] [Cross Ref]
Patel JA, Nair S, Grady J, Revai K, Victor S, Brasier AR, Chonmaitree T. Systemic Cytokine Response Profiles Associated With Respiratory Virus-Induced Acute Otitis Media. Pediatr Infect Dis J. 2009 [PMC free article] [PubMed]
Patel JA, Nair S, Revai K, Grady J, Saeed K, Matalon R, Block S, Chonmaitree T. Association of proinflammatory cytokine gene polymorphisms with susceptibility to otitis media. Pediatrics. 2006;118:2273–2279. doi: 10.1542/peds.2006-0764. [PubMed] [Cross Ref]
Emonts M, Veenhoven RH, Wiertsema SP, Houwing-Duistermaat JJ, Walraven V, de Groot R, Hermans PW, Sanders EA. Genetic polymorphisms in immunoresponse genes TNFA, IL6, IL10, and TLR4 are associated with recurrent acute otitis media. Pediatrics. 2007;120:814–823. doi: 10.1542/peds.2007-0524. [PubMed] [Cross Ref]

Articles from BMC Medical Genetics are provided here courtesy of
BioMed Central