Genetic and functional evidence for a locus controlling otitis media at chromosome 10q26.3
© Rye et al.; licensee BioMed Central Ltd. 2014
Received: 18 October 2012
Accepted: 21 January 2014
Published: 6 February 2014
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.
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.
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.
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.
KeywordsAcute otitis media Otitis media with effusion Genetic polymorphisms Linkage Association Raine study WAFSOM Australia
Otitis media (OM) is a common childhood disease characterised by the presence of infection (acute otitis media; AOM) or fluid (otitis media with effusion; OME) in the middle ear cavity. Most children experience at least one episode of AOM by 1 year of age with up to 40% experiencing recurrent episodes of AOM (rAOM; ≥3 episodes of AOM in 6 months) or chronic episodes of OME (COME; middle ear effusion (MEE) ≥3 months) in childhood . Recurrent disease can result in perforation of the tympanic membrane and/or conductive hearing loss, leading to deficits in language development and poor educational outcomes. Treatment for recurrent disease may include insertion of tympanostomy tubes. High prevalence rates of OM result in substantial health care related costs and significant childhood morbidity in many countries .
The causal mechanisms that lead to recurrent disease are poorly understood but familial clustering and high heritability estimates point to a genetic component. Using self-reported AOM history from 2,570 Norwegian twin pairs, Kvaerner et al.  determined heritability estimates of 0.45 in males and 0.74 in females. Two prospective twin and triplet studies have confirmed these estimates. Casselbrant et al. reported heritability estimates of 0.73 for time with MEE in the first two years of life. Likewise, Rovers et al. reported heritability estimates of 0.49 at 2 years increasing to 0.71 at 4 years with a concomitant decrease in shared environment estimates from 0.41 at 2 years to 0.16 by 4 years of age. Overall these epidemiological studies confirm that susceptibility to OM has a substantial genetic component that increases with age.
To date, there have been few studies undertaken to pinpoint the genes involved. Two genome-wide linkage scans using multi-case families of Caucasian origin have identified specific regions of the genome harbouring putative susceptibility genes. The first identified two regions of linkage on chromosome 10q26.3 and 19q13.43, and a further region on 3p25.3 after conditioning on the linked regions, suggesting epistatic interactions . The second identified two different regions of linkage on 17q12 and 10q22.3 . Whilst the region of 19q13.43 has since been refined , there has been no replication of any of these regions in independent studies and the causal genes underlying them have yet to be identified.
To determine whether these genomic regions are important in recurrent/severe OM in a Western Australian population we have carried out linkage analysis using families who contain at least two individuals diagnosed with rAOM or chronic OME (COME) recruited to the Western Australian Family Study of Otitis Media (WAFSOM) . To identify the putative disease susceptibility locus in the only region that replicated in this population, association using SNPs was undertaken in the Western Australian Pregnancy Cohort (Raine) Study . Functional and bioinformatic studies were used to further clarify the putative etiological gene.
Sample collection and phenotype definition
Linkage studies utilised samples from the WAFSOM where probands with a history of tympanostomy tube insertion due to rAOM or COME were identified from the records of collaborating ear, nose and throat (ENT) specialists (HC and SV) as previously described . Parents and full siblings with a history of recurrent disease, defined as ≥3 physician diagnosed episodes of AOM or tympanostomy tube insertion for rAOM or COME, were also invited to participate. No exclusion was made on the basis of ethnicity or gender; 93.2% of the families in the WAFSOM self-identified as Caucasian. Recruitment to the WAFSOM was approved by the Human Ethics Committee at Princess Margaret Hospital for Children (PMH). Written, informed consents were obtained both for participation in the study and for DNA collection from all adults or from the parents of participants less than 18 years of age. For linkage analysis family members with no history of recurrent disease were classified as unaffected whilst all others were classified as unknown. A total of 101 multicase families (107 nuclear families; 208 affected individuals; 468 total individuals) with 2 to 4 affected individuals were included in the linkage analysis.
For association mapping, we utilised data from the Western Australian Pregnancy Cohort (Raine) Study (‘the Raine Study’), a longitudinal cohort of children whose mothers were recruited during early pregnancy . For the purposes of this study, data collected from clinical examinations and parental questionnaires completed each year at the (average) ages of one, two and three years was used to define a phenotype for OM. Children were defined as a case if clinical examination in the first three years of life indicated presence of inflamed, retracted or scarred TM, MEE or tympanostomy tubes in situ. Participants were also classified as a case where parents’ yearly reports indicated ≥3 episodes of AOM had occurred up to the age of 3 yrs; 35% of cases qualified on the basis of yearly questionnaire criteria alone. Children with no clinical or parental reported history of OM by the age of 3 yrs were classified as controls. Based on questionnaire data, 94% of the Raine Study participants self-identified their ethnicity as Caucasian. We used a subset of 831 Raine Study participants (256 cases and 575 controls) for whom both genome-wide data from an Illumina 660 W Quad Beadchip and complete epidemiological data for covariates were available . Covariates included day care attendance at <3 years of age, allergy diagnosed at <3 years of age, and non-exclusive breastfeeding from <6 months of age, all of which were significant risk factors for OM in the Raine Study cohort . Recruitment to the Raine Study and all follow-ups were approved by the Human Ethics Committee at King Edward Memorial Hospital and/or PMH, with specific adult re-consent for DNA for those individuals participating in this study.
High molecular weight genomic DNA was extracted from a 2 mL saliva sample using the Oragene technology (DNA Genotek) as per manufacturer’s instructions for the WAFSOM study. Extracted DNA was re-suspended in TE buffer, quantified by spectrophotometry, and stored at 50 ng/μl at -20C.
For the majority of Raine Study participants genomic DNA (gDNA) was extracted from whole blood collected at the 14- or 17-year follow-up via venipuncture utilising 4 mL K2EDTA vacuum tubes. Extraction of DNA from whole blood was performed utilising Qiagen PureGene chemistry. For a small subset of individuals (~5%), gDNA was extracted from saliva using the Oragene technology (DNA Genotek) as per manufacturer’s instructions. Samples were quantified using spectrophotometry, diluted to a normalised concentration with reduced EDTA TE buffer and stored at −80°C.
For linkage analysis a minimum of four microsatellite markers  spanning each reported [6, 7] region of linkage were identified. Marker names and primer sequences are provided in Additional file 1: Table S1. For 10q26.3, 19q13.43 and 3p25.3 the reported marker at the peak of linkage  was chosen in addition to four flanking markers. Association mapping across the 10q26.3 region was undertaken in the Raine Study cohort by extracting data for 10,185 SNPs (2,270 genotyped, 7,915 by imputation ) from cleaned , imputed  Illumina 660 W-Quad Beadchip data for 831 individuals (256 cases and 575 controls).
Microsatellite PCRs were performed on WAFSOM DNAs using 15 ng genomic DNA, 1X PCR buffer, 2.6 mM MgCl2, 0.45 mM dNTPs, 0.04 M betaine, 0.87 μM - 1.1 μM forward and reverse primers, and 0.25U of AmpliTaq Polymerase. Cycling was performed using a touchdown protocol . Amplification products were pooled, fragments resolved using an Applied Biosystems 3130xl Genetic Analyzer (Life Technologies, California) and alleles assigned via GeneMapper v4.1. All microsatellites were in Hardy-Weinberg equilibrium in genetically unrelated founders. Mendelian inconsistencies were identified using the Pedcheck software  and corrected or removed prior to analysis.
Genotype data for the Raine Study was generated at the Centre for Applied Genomics (Toronto, Canada) using an Illumina 660 W Quad Beadchip (Illumina, San Diego, California). Quality control (QC) and imputation for the Raine Study has previously been described . Briefly, QC checks were performed for individuals and SNPs using PLINK . Individuals were excluded based on gender mismatch, low genotyping rate (<97%), related to other participants or low level of heterozygosity. SNPs were excluded on the basis of deviation from Hardy-Weinberg Equilibrium (HWE P < 5.7 × 10-7), a genotype call rate <95%, or minor allele frequency <1%. For the chromosome 10q26.3 region, a total of 2,270 genotyped SNPs passed QC checks. Imputation was performed using MACH v1.0.16 with CEU samples from the HapMap Phase2 (Build 36, release 22) used as a reference panel and a default threshold of r2 > 0.3. After imputation, 10,185 SNPs (2,270 genotyped, 7,915 imputed) were available for analysis.
Linkage and association analysis
For microsatellite markers single-point and multi-point non-parametric linkage (NPL) analysis was carried out using all 101 multicase families with the Sall scoring function in the program ALLEGRO  with genetic map distances obtained from the Rutgers Map v2 . Conditional linkage analysis was also performed in ALLEGRO using only families showing evidence for linkage (Zlr >0) at a specified region.
Association mapping of the 10q26 region was performed using ProbABEL  for 10,185 SNPs in the region 127 Mb to qter. Power approximations  estimated for a disease prevalence of 0.4 show that the Raine Study has 86% power to detect associations at an alpha level of P = 1 × 10-5 with genotype relative risk (GRR) of 1.5 for SNPs with a MAF = 0.2 and 95% power at P = 1 × 10-4. Analysis was initially performed adjusting for population substructure (i.e. including the first two Principal Components; PCs), and then repeated, where appropriate, adjusting for PCs and independently associated covariates. If a strict Bonferroni correction was applied then the P-value required to correct for the number of SNPs analysed in the Raine Study data is P = 4.9 × 10-6 (P = 0.05/10,185). However, this is likely to be highly conservative due to the known presence of linkage disequilibrium (LD) in this region (i.e. 629 LD blocks across this 10q26 region, data not shown). Therefore, as this is a replicated linkage region [6, 7], a relaxed threshold of P < 10-3 was used to identify SNPs/genes of interest for follow-up. Logistic regression modelling for independent effects between pairs of SNPs was undertaken in R version 2.15.0 , and results presented as PLRT for the likelihood ratio test comparison. Regional plots of association were created using LocusZoom  in which -log10(P values) were graphed against their chromosomal location. Pairwise LD patterns between all regional SNPs and the respective top SNP were calculated using the HapMap CEU population.
Qualitative PCR expression
For qualitative RT-PCR, RNA was extracted from mononuclear cells isolated from 12 paired adenoid and tonsil samples collected with parent/guardian consent during adenoidectomy and/or tonsillectomy using TRI Reagent (Sigma-Aldrich) according to manufacturer’s instructions . Ethical approval for this part of the study was obtained from the Ethics Committee of PMH. RT-PCR was also performed on two respiratory epithelial and three macrophage cell lines with/without infection with otopathogens. All reagents were from Life Technologies, Australia unless otherwise stated. Adenocarcinomic alveolar basal epithelial cells (A549) were cultured in DMEM, 10% fetal bovine serum (FBS), 2 mM L-glutamine and 100U/ml penicillin/100 μg/ml streptomycin. A pharyngeal carcinoma epithelial cell line (Detroit 562) was grown in MEM with Earle’s salts, 10% FBS, 2 mM L-glutamine, 1X non-essential amino acids, 1 mM sodium pyruvate and 100U/ml penicillin/100 μg/ml streptomycin. Both the promonocytic U937 and monocytic THP-1 cell lines were grown in RPMI 1640, 10% FBS, 2 mM L-glutamine and 100U/ml penicillin/100 μg/ml streptomycin whilst the MonoMac 6 (MM6) mature monocytic cell line additionally required 1X non-essential amino acids and 9 μg/ml OPI bovine insulin (Sigma-Aldrich).
Cell lines were stimulated in vitro in 6-well plates with either non-typeable Haemophilus influenzae (NTHi) or Streptoccocus pneumoniae (SP) at a multiplicity of infection of 10:1. Bacterial cultures were kindly provided by Dr Lea-Ann Kirkham (School of Paediatrics and Child Health, University of Western Australia, Perth). Bacteria were streaked and incubated overnight at 37ºC/5% CO2 on agar plates (blood for SP and chocolate for NTHi). A viable count of bacteria was performed prior to infection using a Helber bacterial counting chamber with Thoma ruling (ProScitech). At the 0 hr time-point, 1 ml of bacteria (1 × 108 cfu/ml), or 1 ml of media for control samples, was added to 1 × 107 cells. After one hour media was replaced to remove non-adhered bacteria. Cells were harvested at 0, 1, 3, 6, 12 and 24 hours using TRI Reagent and RNA extracted accordingly.
A standardised 500 ng of extracted RNA was reverse transcribed using the High Capacity cDNA Reverse Transcription Kit (Life Technologies) as per manufacturer’s instructions. For each cell line an RT negative reaction was used as a negative control. Primers were designed across exon-exon boundaries for four candidate genes (Additional file 1: Table S2) under the linkage peak and for the OAZ1 housekeeping gene . cDNA was amplified using a touchdown PCR protocol  in a volume of 11.5 μl consisting of 20 ng cDNA, a PCR master mix containing 0.85X PCR buffer, 2.6 mM MgCl2, 0.45 mM dNTPs, 0.04 M betaine and 0.25 units of AmpliTaq Polymerase (Life Technologies) and 8.7 ng primer. Products were visualised on a 1.5% agarose gel stained with ethidium bromide.
In silico comparative genomics
The presence of conserved non-coding sequences (CNS), which may harbor important regulatory elements, were examined across the TCERG1L and PPP2R2D genes (± 10 kb flanking sequence) plus the intervening intergenic region. Genomic sequences and associated annotation for human, mouse and rat were exported from Ensembl (Genome Reference Consortium Release 37, Ensembl Release 67) in FASTA and General Feature File (GFF) format, respectively. Global alignment of genomic sequences was performed in Multi-LAGAN [26, 27] and the annotated alignment visualized in SynPlot . CNS regions were defined as regions with a nucleotide sequence conservation level of ≥0.6 that had no associated annotation. To search for putative transcription factor binding sites (TFBS) at SNP locations we used AliBaba v2.1 , and MatInspector v8.0.5  with a matrix similarity parameter >0.75. We also assessed selected polymorphisms (± 1000 bp flanking sequence) to determine whether they are located within a CpG island, defined using a criteria of C/G content >50%, ratio of observed to expected CpG >0.6 and length >200 bp. Repetitive elements were firstly masked using RepeatMasker  following which CpG islands were identified using CpGIsland Searcher  and CpGPlot . Pairwise LD between polymorphisms in the Raine Study was calculated using Haploview 4.1 .
Genetic map location (cM) Marshfield
Physical location (Start; Mb) NCBI build 36
Results of multi-point and single-point non-parametric linkage analysis in the WAFSOM for the 10q26.3 region stratified by rAOM (85 families; 174 affected individuals) and COME (53 families; 115 affected individuals)
Association mapping of the 10q26.3 region in the Raine Study cohort
Comparison of results from the WAFSOM linkage analysis for chromosome 10q26 with those published by Daly et al.
Genetic map position (cM) Marshfield
Physical location (Start; Mb) build 36
Daly LOD (Spt)
Daly LOD (Mpt)
Nearest genes Proximal_Distal
Testing for independent effects in genes across the chromosome 10q26 region
Given evidence of an independent effect at rs9418832 we also adjusted for this SNP (Figure 2B). When adjusting for rs9418832 significance was lost across the DOCK1 region, again indicating a single main effect within this region, whilst improved significance was observed at rs7922424 (P = 6.47 × 10-6), ADAM12 (top SNP rs7902734; P = 2.05 × 10-4) and TCERG1L (top SNP rs2996081; P = 4.26 × 10-4). The most significant association after adjusting for both rs7922424 and rs9418832 was observed at rs7902734 in ADAM12 (Figure 2C; P = 4.52 × 10-4), indicating that ADAM12 may also be independently associated with OM.
Likelihood ratio tests (LRT) to determine independent effects across the 10q26.3 region
Adding ADAM12 SNP
rs9418832 + rs7902734
1.94 × 10-4
rs7922424 + rs7902734
rs9418832 + rs7922424
rs9418832 + rs7922424 + rs7902734
4.37 × 10-4
Adding DOCK1 SNP
rs7902734 + rs9418832
2.07 × 10-5
rs7922424 + rs9418832
5.39 × 10-5
rs7902734 + rs7922424
rs7902734 + rs7922424 + rs9418832
1.45 × 10-5
Adding TCERG1L_PPP2R2D intergenic SNP
rs7902734 + rs7922424
2.31 × 10-5
rs9418832 + rs7922424
7.30 × 10-6
rs7902734 + rs9418832
rs7902734 + rs9418832 + rs7922424
1.61 × 10-5
Expression analysis to support putative functional genes
We also looked at epithelial (A549 & Detroit 562) and macrophage (U937, THP-1 & MM6) cell lines (Figure 3B). Again, ADAM12 and DOCK1 appear to be expressed in all five cell lines whilst PPP2R2D was seen in all cell lines except for MM6. Interestingly, TCERG1L was not expressed in epithelial or macrophage cell lines. TCERG1L expression was also not detected in any cell line following bacterial challenge with two common otopathogens (NTHi & SP; data not shown), suggesting that any role for this gene in OM would be independent of macrophage or epithelial cell inflammatory responses to otopathogens. It should be noted that this expression data is only qualitative. It is likely that quantitative RT-PCR data would add further information regarding expression levels at these candidates.
In silico comparative genomics analysis of the TCERG1L to PPP2R2D region
While statistical and expression analyses do not formally exclude any of the genes showing association under the chromosome 10q26 linkage peak with OM, the statistical analyses favour the interval TCERG1L to PPP2R2D (Additional file 2: Figure S1), which lies directly under the peak of linkage for OM at chromosome 10q26, as the most likely candidate region carrying the etiological variants associated with OM susceptibility. Since the peak of association lies in the intergenic region between these two genes, we carried out an in silico analysis to determine whether the associated SNPs in this region lie in highly conserved non-coding sequences that could contain regulatory elements (Additional file 2: Figure S2). Of the 9 associated SNPs (P < 0.001) within the intergenic TCERG1L to PPP2R2D interval, three are located within highly conserved, non-coding regions (nucleotide conservation level (NCL) ≥ 0.6); rs7037834 (NCL = 0.61), rs7914323 (NCL = 0.61) and rs11813611 (NCL = 0.67). However, rs11813611 lies within a long interspersed nucleotide element (LINEs) repeat region and is unlikely to be functionally relevant. The top SNP, rs7922424, does not lie in a conserved non-coding region but is in strong LD with rs7087384 (r2 = 0.97) and rs7914323 (r2 = 1.00), both of which do fall in conserved non-coding regions (Additional file 2: Figure S3).
We also looked for the presence of transcription factor binding sites (TFBS) that may be co-located with rs7922424 (top SNP), rs7087384 or rs7914232 (within CNS regions). Results (Additional file 1: Table S5) indicate that rs7922424 lies within the consensus sequence for the Upstream Stimulatory Factor-1 (USF1) transcription factor with the binding site abolished in the presence of the minor disease associated (A) allele. On the other hand, the minor disease associated allele (T) at rs7914323 potentially introduces a USF1 binding site. For rs7087384, the minor disease associated (A) allele disrupts a putative CREB/CREBP1 transcription factor binding site.
In addition to TFBS, we looked at the region containing rs7922424, rs7087384 and rs7914323 to determine whether these SNPs fall in regions identified as a CpG island and as such could disrupt biologically relevant DNA methylation at CpG motifs. Results from CpG Plot and CpG Island Searcher indicate that none of these polymorphisms falls at or near a CpG island (data not shown).
Genome-wide linkage analyses have highlighted five regions containing OM disease susceptibility loci on chromosomes 3p25.3, 10q22.3, 10q26.3, 17q12, and 19q13.43 [6–8]. Here, using a Western Australian cohort of children diagnosed with recurrent OM, we found evidence of replication of linkage at 10q26.3. We also found borderline evidence to support a region of linkage at 10q22.3 but did not find evidence for linkage at 3p25.3, 17q12, or 19q13.43, or for any region after conditioning on 10q26.3. Therefore, we focused attention on mapping the genes under the chromosome 10q26.3 linkage peak.
To identify the etiological gene/variant under the 10q26.3 peak, we performed fine-mapping using directed analysis of 2,270 genotyped SNPs (10,185 after imputation) spanning this region available on 256 cases and 575 controls within the Raine Study cohort. We focused our search in the region 127 Mb to qter, and found statistical support for variants in several genes/regions across this region contributing to the linkage peak. This is not an unusual phenomenon in complex diseases, where previous studies have also highlighted multiple genes contributing to peaks of linkage . In this case, for the region of linkage at 10q26.3, at least 49 genes lie within the interval, a number of which could be considered as potential functional candidates for OM susceptibility. Statistical support was observed for independent effects of SNPs at ADAM12, DOCK1 and the intergenic region between TCERG1L and PPP2R2D.
The ADAM12 (A Disintegrin and Metalloproteinase domain 12) gene is a member of the disintegrin and metalloproteinase (ADAM) family of proteins. This gene has been implicated in the epidermal growth factor receptor (EGFR) signalling pathway , which is upregulated in human middle ear epithelial cells in response to tobacco smoke exposure . The ADAM12 protein also interacts with the muscle-specific α-actinin-2 protein, with function centred on myoblast/muscle development showing increased expression during muscle generation . The DOCK1 (Dedicator of cytokinesis 1) protein has roles in phagocytosis of apoptotic cells in concert with ELMO1 during Rac signalling and cellular migration [44, 45]. The DOCK1 gene also has a role during embryogenesis in muscle development with knock-out mouse mutants having decreased skeletal and respiratory muscle tissues . This suggests the apoptotic function of DOCK1 could have a subtle role in the cell death of inflammatory factors and in apoptosis of mucous cells after an immune response to OM, whilst the ADAM12 protein interacts with pathways that could influence expression of inflammatory mediators within the middle ear.
The strongest statistical support for association was in the TCERG1L to PPP2R2D intergenic region, directly below the peak of linkage on chromosome 10q26.3. Very little is known about the function of the transcription elongation regulator like protein (TCERG1L) gene. In recent GWAS studies, variants at or near TCERG1L have been associated with fasting insulin, insulin resistance  and attention deficit disorder . Hypermethylation of the TCERG1L promoter region leading to gene silencing has also been observed in colon cancer . TCERG1L expression is documented  in a variety of tissues, including the brain, lung and eye. In this study, we have also demonstrated that TCERG1L is expressed in adenoids but not in macrophage or epithelial cell lines, either with or without otopathogen infection. Whilst expression in these cells may have been downregulated during immortalization, analogous to downregulation of this gene in cancer cells , the data suggest that any role TCERG1L may play in OM susceptibility is unlikely to occur through the innate inflammatory response to otopathogens. On balance, TCERG1L does not appear to be a strong candidate for OM susceptibility.
In contrast, PPP2R2D is a particularly interesting candidate gene for OM. PPP2R2D is a member of the B family of regulatory subunits of the protein phosphatase 2A (PP2A) and is widely expressed at the protein level in the brain, heart, placenta, skeletal muscle, testis and thymus . This protein has a role as a modulator of the TGF-β/Activin/Nodal pathway , where knockdown of the protein was shown to increase nuclear accumulation and phosphorylation of Smad2. This involvement with Smad2 is of specific interest, as this gene and others in the TGFb pathway have previously been highlighted as candidate susceptibility genes for OM within the WAFSOM cohort . Furthermore, a GWAS carried out in the Raine Study highlighted at least five other members of the TGFb pathway (BMP5, GALNT13, NELL1, TGFB3 and BPIFA1) as candidates for OM susceptibility .
The strongest signal for association in our study lay within the intergenic region between the TCERG1L and PPP2R2D genes. Many association signals in complex diseases have been found to lie within intergenic regions , leading to a search for potential regulatory functions within those regions. Whilst the top SNP (rs7922424) in the intergenic region does not itself lie within a conserved non-coding region, it is in strong LD with rs7087384 and rs7914323, which do lie in highly conserved non-coding regions. In addition, the minor disease associated alleles at rs7922424 and rs7914323 alter putative binding sites for the upstream stimulatory factor 1 (USF1). The USF1 transcription factor is a member of the helix-loop-helix leucine zipper family and is ubiquitously expressed in a variety of cells . The role of this protein is widespread, ranging from roles in embryonic development  to promoters for a number of activity-induced genes within neuronal nuclei  and as transcription factors that regulate cell-type dependent cellular proliferation . However, our data favours rs7914323 as the most likely regulatory polymorphism in this region. Not only does this SNP lie within a conserved, non-coding region that could harbor regulatory elements, but our analysis shows that the presence of the minor allele potentially eliminates a TFBS for the cAMP response element binding (CREB) transcription factor and its binding protein (CREBP or CBP). This is interesting in the context of OM as the function of CBP has been linked to the TGFb pathway via the recruitment of EVI1 , which is mutated in a mouse model of OM , although the EVI1 gene has not been associated with human susceptibility to OM to date [9, 59].
Overall, analysis of conserved, non-coding regions and putative TFBS sites indicate a number of regulatory elements that lie within this intergenic region that can potentially be influenced by polymorphisms associated with OM susceptibility. It is not possible to determine from these data which gene these regulatory elements may influence. However, the TFBS consensus sequences disrupted by these polymorphisms all lie on the forward strand, upstream of the PPP2R2D promoter. The TCERG1L gene on the other hand is encoded on the reverse strand. Taken together these observations strengthen the evidence for PPP2R2D as the likely gene contributing to linkage at 10q26.3.
Using the resources available within the family-based WAFSOM cohort the results of our study provide the first replication of linkage to OM susceptibility at chromosome 10q26.3. Subsequent association analysis using data available within the longitudinal Raine Study show that multiple genes could contribute to this linkage peak, however the weight of evidence supports PPP2R2D as the more likely functional candidate in this linkage region. It has previously been noted that there is phenotypic heterogeneity between WAFSOM and the Raine Study . The phenotype defined for the Raine Study is biased towards the milder end of the OM spectrum, being based on yearly clinical examinations and parental reports of ≥3 AOM episodes during the first three years of life. The WAFSOM cohort on the other hand is biased towards the more severe end of the OM spectrum, being based on children recommended for grommet insertion. This suggests that PPP2R2D plays a role in susceptibility to OM per se. In summary our results point to a role for PPP2R2D, a TGF-β/Activin/Nodal pathway modulator, as the more likely functional candidate for OM susceptibility at chromosome 10q26.3. This contributes to the growing evidence for a role for the TGFb pathway in susceptibility to this important childhood disease.
We would like to thank all the families who have participated in the Western Australian Family Study of Otitis Media. We would also like to thank the Raine Study participants for their contribution to the Raine Study as well as the Raine Study Team for cohort co-ordination and data collection. Thanks also go to Ms Nicole Warrington for her work with the QC analysis of the Raine Study GWAS data. The work was supported by iVEC through the use of advanced computing resources located at iVEC@Murdoch.
This work was supported by funds to SEJ from a Brightspark Foundation Fellowship, a Raine Medical Research Foundation Priming Grant and from core funds to JMB from The Stan Perron Foundation, the University of Western Australia (UWA) and the Western Australian State Government. MSR is supported by an Australian Post Graduate Scholarship. The Raine Study has been supported by the NH&MRC over the last 20 years with funding for Core Management provided by UWA, The Raine Medical Research Foundation at UWA, the UWA Faculty of Medicine, Dentistry and Health Sciences, the Telethon Institute for Child Health Research, the Women and Infants Research Foundation and Curtin University. Raine Study Illumina 660 W-Quad Beadchip Data was supported by the NH&MRC (ID 572613).
- Rovers MM, Schilder AGM, Zielhuis GA, Rosenfeld RM: Otitis media. Lancet. 2004, 363: 465-473. 10.1016/S0140-6736(04)15495-0.View ArticlePubMedGoogle Scholar
- Kong K, Coates HLC: Natural history, definitions, risk factors and burden of otitis media. Med J Aust. 2009, 191: S39-S43.PubMedGoogle Scholar
- Kvaerner KJ, Tambs K, Harris JR, Magnus P: Distribution and heritability of reccurent ear infections. Ann Otol Rhinol Laryngol. 1997, 106 (8): 624-632.View ArticlePubMedGoogle Scholar
- Casselbrant M, 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 (22): 2125-2130. 10.1001/jama.282.22.2125.View ArticlePubMedGoogle Scholar
- 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 (10): 958-964. 10.1093/aje/155.10.958.View ArticlePubMedGoogle Scholar
- 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. 10.1086/426061.View ArticlePubMedPubMed CentralGoogle Scholar
- Casselbrant ML, Mandel EM, Jung J, Ferrell RE, Tekely K, Szatkiewicz JP, Ray A, Weeks DE: Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions. BMC Med Genet. 2009, 10: 85-View ArticlePubMedPubMed CentralGoogle Scholar
- Chen W-M, Allen EK, Mychaleckyj JC, Chen F, Hou X, Rich SS, Daly KA, Sale MM: Significant linkage at chromosome 19q for otitis media with effusion and/or recurrent otitis media (COME/ROM). BMC Med Genet. 2011, 12 (1): 124-10.1186/1471-2350-12-124.View ArticlePubMedPubMed CentralGoogle Scholar
- Rye MS, Wiertsema SP, Scaman ESH, Oommen J, Sun W, Coates HL, Vijayasekaran S, Francis RW, Ang W, Pennell CE, et al: FBXO11, a regulator of the TGFβ pathway, is associated with severe otitis media in Western Australian children. Genes Immun. 2011, 12 (5): 352-359. 10.1038/gene.2011.2.View ArticlePubMedGoogle Scholar
- Newnham JP, Evans SF, Michael CA, Stanley FJ, Landau LI: Effects of frequent ultrasound during pregnancy: a randomised controlled trial. Lancet. 1993, 342: 887-891. 10.1016/0140-6736(93)91944-H.View ArticlePubMedGoogle Scholar
- Rye MS, Warrington NM, Scaman ESH, Vijayasekaran S, Coates HL, Anderson D, Pennell CE, Blackwell JM, Jamieson SE: Genome-wide association study to identify the genetic determinants of otitis media susceptibility in childhood. PLoS One. 2012, 7 (10): e48215-10.1371/journal.pone.0048215.View ArticlePubMedPubMed CentralGoogle Scholar
- Matise TC, Chen FCW, De La Vega FM, Hansen M, He C, Hyland FCL, Kennedy GC, Kong X, Murray SS, Ziegle JS, et al: A second-generation combined linkage-physical map of the human genome. Genome Res. 2007, 17 (12): 1783-1786. 10.1101/gr.7156307.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Y, Willer CJ, Sanna S, Abecasis GR: Genotype imputation. Annu Rev Genomics Hum Genet. 2009, 10: 387-406. 10.1146/annurev.genom.9.081307.164242.View ArticlePubMedPubMed CentralGoogle Scholar
- Teo YY: Common statistical issues in genome-wide association studies: a review on power, data quality control, genotype calling and population structure. Curr Opin Lipidol. 2008, 19 (2): 133-143. 10.1097/MOL.0b013e3282f5dd77.View ArticlePubMedGoogle Scholar
- Korbie DJ, Mattick JS: Touchdown PCR for increased specificity and sensitivity in PCR amplification. Nat Protoc. 2008, 3: 1452-1456. 10.1038/nprot.2008.133.View ArticlePubMedGoogle Scholar
- O’Connell JR, Weeks DE: PedCheck: a program for identifying marker typing incompatibilities in linkage analysis. Am J Hum Genet. 1998, 10: A288-Google Scholar
- Freathy RM, Mook-Kanamori DO, Sovio U, Prokopenko I, Timpson NJ, Berry DJ, Warrington NM, Widen E, Hottenga JJ, Kaakinen M, et al: Variants in ADCY5 and near CCNL1 are associated with fetal growth and birth weight. Nat Genet. 2010, 42 (5): 430-435. 10.1038/ng.567.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 toolset for whole-genome association and population-based linkage analysis. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.View ArticlePubMedPubMed CentralGoogle Scholar
- Gudbjartsson DF, Jonasson K, Frigge ML, Kong A: Allegro, a new computer program for multipoint linkage analysis. Nat Genet. 2000, 25 (1): 12-13. 10.1038/75514.View ArticlePubMedGoogle Scholar
- Aulchenko YS, Struchalin MV, van Duijn CM: ProbABEL package for genome-wide association analysis of imputed data. BMC Bioinformatics. 2010, 11: 134-10.1186/1471-2105-11-134.View ArticlePubMedPubMed CentralGoogle Scholar
- Skol AD, Scott LJ, Abecasis GR, Boehnke M: Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet. 2006, 38: 209-213. 10.1038/ng1706.View ArticlePubMedGoogle Scholar
- R: A Language and Environment for Statistical Computing. http://www.R-project.org,
- Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ: LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010, 26 (18): 2336-2337. 10.1093/bioinformatics/btq419.View ArticlePubMedPubMed CentralGoogle Scholar
- Chomczynski P: A reagent for the single-step simultaneous isolation of RNA, DNA and proteins from cell and tissue samples. Biotechniques. 1993, 15 (3): 532-537.PubMedGoogle Scholar
- de Jonge HJ, Fehrmann RS, Bont ES, Hofstra RM, Gerbens F, Kamps WA, de Vries EG, van der Zee AG, te Meerman GJ, ter Elst A: Evidence based selection of housekeeping genes. PLoS One. 2007, 2 (9): e898-10.1371/journal.pone.0000898.View ArticlePubMedPubMed CentralGoogle Scholar
- Brudno M, Chapman M, Gottgens B, Batzoglou S, Morgenstern B: Fast and sensitive multiple alignment of large genomic sequences. BMC Bioinformatics. 2003, 4: 66-10.1186/1471-2105-4-66.View ArticlePubMedPubMed CentralGoogle Scholar
- Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, Sidow A, Batzoglou S: LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 2003, 13 (4): 721-731. 10.1101/gr.926603.View ArticlePubMedPubMed CentralGoogle Scholar
- Gottgens B, Gilbert JG, Barton LM, Grafham D, Rogers J, Bentley DR, Green AR: Long-range comparison of human and mouse SCL loci: localized regions of sensitivity to restriction endonucleases correspond precisely with peaks of conserved noncoding sequences. Genome Res. 2001, 11 (1): 87-97. 10.1101/gr.153001.View ArticlePubMedPubMed CentralGoogle Scholar
- Grabe N: AliBaba2: context specific identification of transcription factor binding sites. In Silico Biol. 2002, 2 (1): S1-S15.PubMedGoogle Scholar
- Quandt K, Frech K, Karas H, Wingender E, Werner T: MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res. 1995, 23 (23): 4878-4884. 10.1093/nar/23.23.4878.View ArticlePubMedPubMed CentralGoogle Scholar
- RepeatMasker Open-3.0. http://www.repeatmasker.org,
- Takai D, Jones PA: Comprehensive analysis of CpG islands in human chromosomes 21 and 22. Proc Natl Acad Sci U S A. 2002, 99 (6): 3740-3745. 10.1073/pnas.052410099.View ArticlePubMedPubMed CentralGoogle Scholar
- Larsen F, Gundersen G, Lopez R, Prydz H: CpG islands as gene markers in the human genome. Genomics. 1992, 13 (4): 1095-1107. 10.1016/0888-7543(92)90024-M.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Weeks DE, Lathrop GM: Polygenic disease: methods for mapping complex disease traits. Trends Genet. 1995, 11 (12): 513-519. 10.1016/S0168-9525(00)89163-5.View ArticlePubMedGoogle Scholar
- Hovatta I, Lichtermann D, Juvonen H, Suvisaari J, Terwilliger JD, Arajarvi R, Kokko-Sahin ML, Ekelund J, Lonnqvist J, Peltonen L: Linkage analysis of putative schizophrenia gene candidate regions on chromosomes 3p, 5q, 6p, 8p, 20p and 22q in a population-based sampled Finnish family set. Mol Psychiatry. 1998, 3 (5): 452-457. 10.1038/sj.mp.4000443.View ArticlePubMedGoogle Scholar
- Roberts SB, MacLean CJ, Neale MC, Eaves LJ, Kendler KS: Replication of linkage studies of complex traits: an examination of variation in location estimates. Am J Hum Genet. 1999, 65 (3): 876-884. 10.1086/302528.View ArticlePubMedPubMed CentralGoogle Scholar
- van den Aardweg MT, Schilder AG, Herkert E, Boonacker CW, Rovers MM: Adenoidectomy for otitis media in children. Cochrane Database Syst Rev. 2010, 1: CD007810-PubMedGoogle Scholar
- Kadhim AL, Spilsbury K, Semmens JB, Coates HL, Lannigan FJ: Adenoidectomy for middle ear effusion: a study of 50,000 children over 24 years. Laryngoscope. 2007, 117 (3): 427-433. 10.1097/MLG.0b013e31802c938b.View ArticlePubMedGoogle Scholar
- Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA: Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci. 2009, 106 (23): 9362-9367. 10.1073/pnas.0903103106.View ArticlePubMedPubMed CentralGoogle Scholar
- Mori S, Tanaka M, Nanba D, Nishiwaki E, Ishiguro H, Higashiyama S, Matsuura N: PACSIN3 binds ADAM12/meltrin alpha and up-regulates ectodomain shedding of heparin-binding epidermal growth factor-like growth factor. J Biol Chem. 2003, 278 (46): 46029-46034. 10.1074/jbc.M306393200.View ArticlePubMedGoogle Scholar
- Cho J-G, Woo J-S, Lee H-M, Jung HH, Hwang S-J, Chae S: Effects of cigarette smoking on mucin production in human middle ear epithelial cells. Int J Pediatr Otorhinolaryngol. 2009, 73 (10): 1447-1451. 10.1016/j.ijporl.2009.07.016.View ArticlePubMedGoogle Scholar
- Galliano M-F, Huet C, Frygelius J, Polgren A, Wewer UM, Engvall E: Binding of ADAM12, a marker of skeletal muscle regeneration, to the muscle-specific actin-binding protein, alpha-actinin-2, is required for myoblast fusion. J Biol Chem. 2000, 275: 13933-13939. 10.1074/jbc.275.18.13933.View ArticlePubMedGoogle Scholar
- Hasegawa H, Kiyokawa E, Tanaka S, Nagashima K, Gotoh N, Shibuya M, Kurata T, Matsuda M: DOCK180, a major CRK-binding protein, alters cell morphology upon translocation to the cell membrane. Mol Cell Biol. 1996, 16 (4): 1770-1776.View ArticlePubMedPubMed CentralGoogle Scholar
- Gumienny TL, Brugnera E, Tosello-Trampont A-C, Kinchen JM, Haney LB, Nishiwaki K, Walk SF, Nemergut ME, Macara IG, Francis R, et al: CED-12/ELMO, a Novel Member of the CrkII/Dock180/Rac Pathway, Is Required for Phagocytosis andCell Migration. Cell. 2001, 107 (1): 27-41. 10.1016/S0092-8674(01)00520-7.View ArticlePubMedGoogle Scholar
- Laurin M, Fradet N, Blangy A, Hall A, Vuori K, Cote J-F: The atypical Rac activator Dock180 (Dock1) regulates myoblast fusion in vivo. Proc Natl Acad Sci. 2008, 105 (40): 15446-15451. 10.1073/pnas.0805546105.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen G, Bentley A, Adeyemo A, Shriner D, Zhou J, Doumatey A, Huang H, Ramos E, Erdos M, Gerry N, et al: Genome-wide association study identifies novel loci association with fasting insulin and insulin resistance in African Americans. Hum Mol Genet. 2012, 21: 4530-4536. 10.1093/hmg/dds282.View ArticlePubMedPubMed CentralGoogle Scholar
- Neale BM, Medland S, Ripke S, Anney RJL, Asherson P, Buitelaar J, Franke B, Gill M, Kent L, Holmans P, et al: Case–control genome-wide association study of attention-deficit/hyperactivity disorder. J Am Acad Child Adolesc Psychiatry. 2010, 49 (9): 906-920. 10.1016/j.jaac.2010.06.007.View ArticlePubMedPubMed CentralGoogle Scholar
- Yi J, Dhir M, Guzzetta A, Iacobuzio-Donahue C, Heo K, Yang K, Suzuki H, Toyota M, Kim H-M, Ahuja N: DNA methylation biomarker candidates for early detection of colon cancer. Tumor Biology. 2012, 33 (2): 363-372. 10.1007/s13277-011-0302-2.View ArticlePubMedPubMed CentralGoogle Scholar
- Thierry-Mieg D, Thierry-Mieg J: AceView: a comprehensive cDNA-supported gene and transcripts annotation. Genome Biol. 2006, 7 (Supp 1): S12-View ArticlePubMedPubMed CentralGoogle Scholar
- Strack S, Chang D, Zaucha JA, Colbran RJ, Wadzinski BE: Cloning and characterization of Bδ, a novel regulatory subunit of protein phosphatase 2A. FEBS Lett. 1999, 460 (3): 462-466. 10.1016/S0014-5793(99)01377-0.View ArticlePubMedGoogle Scholar
- Batut J, Schmierer B, Cao J, Raftery LA, Hill CS, Howell M: Two highly related regulatory subunits of PP2A exert opposite effects on TGF-β/Activin/Nodal signalling. Development. 2008, 135 (17): 2927-2937. 10.1242/dev.020842.View ArticlePubMedGoogle Scholar
- Atchley WR, Fitch WM: A natural classification of the basic helix-loop-helix class of transcription factors. Proc Natl Acad Sci U S A. 1997, 94: 5172-5176. 10.1073/pnas.94.10.5172.View ArticlePubMedPubMed CentralGoogle Scholar
- Sirito M, Lin Q, Deng JM, Behringer RR, Sawadogo M: Overlapping roles and asymmetrical cross-regulation of the USF proteins in mice. Proc Natl Acad Sci U S A. 1998, 95 (7): 3758-3763. 10.1073/pnas.95.7.3758.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen WG, West AE, Tao X, Corfas G, Szentirmay MN, Szentirmay MN, Sawadogo M, Vinson C, Greenberg C: Upstream stimulatory factors are mediators of ca2 + −responsive transcription in neurons. J Neurosci. 2003, 23 (7): 2572-2581.PubMedGoogle Scholar
- Qyang Y, Luo X, Lu T, Ismail PM, Krylov D, Vinson C, Sawadogo MI: Cell-type-dependent activity of the ubiquitous transcription factor usf in cellular proliferation and transcriptional activation. Mol Cell Biol. 1999, 19 (2): 1508-1517.View ArticlePubMedPubMed CentralGoogle Scholar
- Alliston T, Ko TC, Cao Y, Feng X-H, Chang C, Derynck R: Repression of bone morphogenetic protein and activin-inducible transcription by Evi-1. J Biol Chem. 2005, 280 (25): 24227-24237. 10.1074/jbc.M414305200.View ArticlePubMedGoogle Scholar
- Parkinson N, Hardisty-Hughes RE, Tateossian H, Tsai H-T, Brooker D, Morse S, Lalane Z, MacKenzie F, Fray M, Glenister P, et al: Mutation at the evi1 locus in junbo mice causes susceptibility to otitis media. PLoS Genet. 2006, 2 (10): e149-10.1371/journal.pgen.0020149.View ArticlePubMedPubMed CentralGoogle Scholar
- Sale MM, Chen W-M, Weeks DE, Mychaleckyj JC, Hou X, Marion M, Segade F, Casselbrant ML, Mandel EM, Ferrell RE, et al: Evaluation of 15 Functional Candidate Genes for Association with Chronic Otitis Media with Effusion and/or Recurrent Otitis Media (COME/ROM). PLoS One. 2011, 6 (8): e22297-10.1371/journal.pone.0022297.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/15/18/prepub
This article is published under license to BioMed Central Ltd. 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.