Abstract
BACKGROUND: Otitis media (OM) is a common childhood disease characterised by middle ear effusion and inflammation. Susceptibility to recurrent acute OM and chronic OM with effusion is 40-70% heritable. Linkage studies provide evidence for multiple putative OM susceptibility loci. This study attempts to replicate these linkages in a Western Australian (WA) population, and to identify the etiological gene(s) in a replicated region. METHODS: Microsatellites were genotyped in 468 individuals from 101 multicase families (208 OM cases) from the WA Family Study of OM (WAFSOM) and non-parametric linkage analysis carried out in ALLEGRO. Association mapping utilized dense single nucleotide polymorphism (SNP) data extracted from Illumina 660 W-Quad analysis of 256 OM cases and 575 controls from the WA Pregnancy Cohort (Raine) Study. Logistic regression analysis was undertaken in ProbABEL. RT-PCR was used to compare gene expression in paired adenoid and tonsil samples, and in epithelial and macrophage cell lines. Comparative genomics methods were used to identify putative regulatory elements and transcription factor binding sites potentially affected by associated SNPs. RESULTS: Evidence for linkage was observed at 10q26.3 (Zlr = 2.69; P = 0.0036; D10S1770) with borderline evidence for linkage at 10q22.3 (Zlr = 1.64; P = 0.05; D10S206). No evidence for linkage was seen at 3p25.3, 17q12, or 19q13.43. Peak association at 10q26.3 was in the intergenic region between TCERG1L and PPP2R2D (rs7922424; P = 9.47 × 10-6), immediately under the peak of linkage. Independent associations were observed at DOCK1 (rs9418832; P = 7.48 × 10-5) and ADAM12 (rs7902734; P = 8.04 × 10-4). RT-PCR analysis confirmed expression of all 4 genes in adenoid samples. ADAM12, DOCK1 and PPP2R2D, but not TCERG1L, were expressed in respiratory epithelial and macrophage cell lines. A significantly associated polymorphism (rs7087384) in strong LD with the top SNP (rs7922424; r2 = 0.97) alters a transcription factor binding site (CREB/CREBP) in the intergenic region between TCERG1L and PPP2R2D. CONCLUSIONS: OM linkage was replicated at 10q26.3. Whilst multiple genes could contribute to this linkage, the weight of evidence supports PPP2R2D, a TGF-β/Activin/Nodal pathway modulator, as the more likely functional candidate lying immediately under the linkage peak for OM susceptibility at chromosome 10q26.3.