Skip to main content
  • Research article
  • Open access
  • Published:

Genetic and functional evidence for a locus controlling otitis media at chromosome 10q26.3

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.

Peer Review reports

Background

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 [1]. 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 [2].

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. [3] 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.[4] reported heritability estimates of 0.73 for time with MEE in the first two years of life. Likewise, Rovers et al.[5] 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 [6]. The second identified two different regions of linkage on 17q12 and 10q22.3 [7]. Whilst the region of 19q13.43 has since been refined [8], 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) [9]. 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 [10]. Functional and bioinformatic studies were used to further clarify the putative etiological gene.

Methods

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 [9]. 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 [10]. 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 [11]. 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 [11]. 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.

DNA extraction

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.

Marker selection

For linkage analysis a minimum of four microsatellite markers [12] 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 [6] 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 [13]) from cleaned [14], imputed [13] Illumina 660 W-Quad Beadchip data for 831 individuals (256 cases and 575 controls).

Genotyping

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 [15]. 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 [16] 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 [17]. Briefly, QC checks were performed for individuals and SNPs using PLINK [18]. 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 [19] with genetic map distances obtained from the Rutgers Map v2 [12]. 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 [20] for 10,185 SNPs in the region 127 Mb to qter. Power approximations [21] 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 [22], and results presented as PLRT for the likelihood ratio test comparison. Regional plots of association were created using LocusZoom [23] 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 [24]. 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 [25]. cDNA was amplified using a touchdown PCR protocol [15] 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 [28]. 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 [29], and MatInspector v8.0.5 [30] 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 [31] following which CpG islands were identified using CpGIsland Searcher [32] and CpGPlot [33]. Pairwise LD between polymorphisms in the Raine Study was calculated using Haploview 4.1 [34].

Results

Linkage analysis

Results of linkage analysis using the 101 multi-case families available in WAFSOM are presented in Table 1. Multi-point (peak Zlr = 2.69; P = 0.0036 at D10S1770) and single-point (Zlr = 2.27; P = 0.012 at D10S1770) NPL analysis provided evidence for linkage between markers at chromosome 10q26.3 and OM, contributed to by both the rAOM and COME phenotypes (Table 2). Analysis of Caucasian families alone did not reveal a substantial difference in single-point results at any region (data not shown); therefore, all multi-case families were used in subsequent analyses. Linkage analysis conditional on the linkage at 10q26.3 provided no evidence of epistatic interactions with 19q13.43 or 3p25.3 in this study population (data not shown), as previously reported [6]. Borderline evidence (multi-point Zlr = 1.65; nominal P = 0.05) for replication of linkage was also observed at 10q22.3 whilst no evidence for linkage (Zlr < 1.65; nominal P > 0.05) was observed at 3p25.3, 17q12 or 19q13.43.

Table 1 Results of multi-point and single-point non-parametric linkage analysis in the WAFSOM across chromosomal regions previously[6, 7] linked to OM
Table 2 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

In our study, the peak of linkage at 10q26.3 is observed at 169 cM (Marshfield Genetic map)/132.57 Mb (NCBI Build 36), with evidence for linkage extending proximally to 130.85 Mb and distally to 135.05 Mb. Daly et al. [6] report clear evidence for linkage only at microsatellite marker D10S212, located at 171 cM/134.29 Mb, with no evidence for linkage ~126 Mb (Table 3). Although the disease locus in a linkage study might normally be found immediately under the peak of linkage [35], the true susceptibility gene can be displaced by up to 10 cM (~10 Mb) from the linkage peak, particularly in smaller samples [36, 37]. Therefore, we extended our search for genes through association mapping over the region 127 Mb to 135.35 Mb (=qter).

Table 3 Comparison of results from the WAFSOM linkage analysis for chromosome 10q26 with those published by Daly et al.[6]

For association mapping, genotype data from the region 127 Mb to 135.35 Mb on chromosome 10 was extracted from Illumina 660 W-Quad BeadChip data in 831 individuals (256 cases and 575 controls) from The Raine Study for whom full covariate data were available [11]. SNP locations on this Illumina chip are provided as NCBI Build 36 base pair locations. Data were available for 10,185 SNPs (2,270 genotyped, 7,915 by imputation). Logistic regression analysis under an additive model, controlling for population substructure using 2 PCs as covariates, provided evidence for several regions of association across the region 127 Mb to qter, with varying degrees of statistical support (Figure 1; Additional file 1: Table S3). In summary, individual top SNPs within genes/regions (P < 10-3) are: rs7902734 intronic in ADAM12 at 127.94 Mb (P = 8.04 × 10-4); rs9418832 intronic in DOCK1 at 128.76 Mb (P = 7.48 × 10-5); rs2996081 intronic in TCERG1L at 132.86 Mb (P = 9.15 × 10-4) and rs7922424 intergenic between TCERG1L and PPP2R2D at 133.21 Mb (P = 9.47 × 10-6). Association signals across these genes/regions were robust to adjustment for independent environmental/clinical covariates (Additional file 1: Table S4).

Figure 1
figure 1

LocusZoom plots of association for the region spanning 127 Mb to qter in the Raine study. Pairwise values of linkage disequilibrium with the top SNP (rs7922424; in purple) were calculated using the HapMap CEU population. The top SNP in other genomic regions have been labeled on the plot.

Testing for independent effects in genes across the chromosome 10q26 region

To evaluate whether the putative regions of association across the 127 Mb to qter region on chromosome 10q26 represented independent signals we repeated the single-point logistic regression analysis adjusting for the most significantly associated SNP (i.e. rs7922424; Figure 2). Upon adjusting for rs7922424 (Figure 2A), the significance at neighbouring SNPs within this intergenic region was lost indicating there is a single main effect across the TCERG1L_PPP2R2D intergenic region. However, improved significance was observed at DOCK1 (top SNP rs9418832; P = 4.77 × 10-5), suggesting that these genes/regions may be independently associated with OM. All other signals were at P > 10-3 after adjusting for rs7922424.

Figure 2
figure 2

LocusZoom plots of association adjusted for the top two SNPs, (A) rs7922424, (B) rs9418832 and (C) rs7922424 and rs9418832. Significance is retained at each SNP after adjusting the analysis for the other. Pairwise values of linkage disequilibrium with the top SNP (in purple) were calculated using the HapMap CEU population.

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.

To formally test for independent effects between ADAM12, DOCK1 and the TCERG1L_PPP2R2D intergenic region, we used forward stepwise logistic regression modelling. When comparing models in which the top SNPs are added to each of the other top SNPs (i.e. alternative 2-SNP model compared to the null 1-SNP model; Table 4) results show that they all add independent effects to each other (PLRT ≤ 0.001). Comparison of models containing all three SNPs compared to a null model containing only 2 of the top SNPs confirms these SNPs all add independent effects to each other (P ≤ 4.37 × 10-4). Overall, the data suggest that variants in the intergenic region between TCERG1L and PPP2R2D, which lies at 133.2 Mb directly under the peak of linkage (Table 3), as well as in DOCK1 (128.8 Mb) and ADAM12 (127.9 Mb), may contribute independently to susceptibility to OM. Therefore, we took these genes/regions forward in evaluating further evidence for the putative etiological gene(s) for OM under the peak of linkage on chromosome 10q26.

Table 4 Likelihood ratio tests (LRT) to determine independent effects across the 10q26.3 region

Expression analysis to support putative functional genes

Information on the expression pattern of several genes in the 10q26 region is currently limited. Therefore, another way to evaluate potential candidacy of disease associated genes is to determine whether they are expressed in relevant tissue sites ex vivo, or in model culture systems in vitro. To address the former, we looked at expression of genes in paired adenoid and tonsil tissue from OM cases (rAOM, OME) compared to controls diagnosed with obstructive sleep disorder (OSD) or recurrent acute tonsillitis (RAT) with no concurrent OM history. The role of adenoids in OM pathogenesis is supported by a systematic review showing that adenoidectomy significantly improves resolution of MEE in OM [38] whereas tonsillectomy alone does not [39]. Figure 3A summarises the RT-PCR expression data across adenoid and tonsil samples and shows that ADAM12, DOCK1, TCERG1L and PPP2R2D are, for the most part, expressed in all adenoid samples irrespective of OM or OSD diagnosis. In tonsil samples, ADAM12 and DOCK1 are again expressed in all samples irrespective of diagnosis. However, TCERG1L and PPP2R2D expression is less consistent in tonsils with at least two samples showing no expression.

Figure 3
figure 3

RT-PCR expression of various candidate genes in (A) adenoids and tonsils of differing phenotypes and (B) in epithelial and macrophage cell lines. OAZ1 is the housekeeping gene. A green diagonally striped box indicates expression was observed whilst a red dotted box indicates no expression was observed.

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).

Discussion

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 [40]. 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 [41], which is upregulated in human middle ear epithelial cells in response to tobacco smoke exposure [42]. 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 [43]. 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 [46]. 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 [47] and attention deficit disorder [48]. Hypermethylation of the TCERG1L promoter region leading to gene silencing has also been observed in colon cancer [49]. TCERG1L expression is documented [50] 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 [49], 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 [51]. This protein has a role as a modulator of the TGF-β/Activin/Nodal pathway [52], 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 [9]. 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 [11].

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 [40], 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 [53]. The role of this protein is widespread, ranging from roles in embryonic development [54] to promoters for a number of activity-induced genes within neuronal nuclei [55] and as transcription factors that regulate cell-type dependent cellular proliferation [56]. 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 [57], which is mutated in a mouse model of OM [58], 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.

Conclusions

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 [11]. 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.

References

  1. Rovers MM, Schilder AGM, Zielhuis GA, Rosenfeld RM: Otitis media. Lancet. 2004, 363: 465-473. 10.1016/S0140-6736(04)15495-0.

    Article  PubMed  Google Scholar 

  2. Kong K, Coates HLC: Natural history, definitions, risk factors and burden of otitis media. Med J Aust. 2009, 191: S39-S43.

    PubMed  Google Scholar 

  3. Kvaerner KJ, Tambs K, Harris JR, Magnus P: Distribution and heritability of reccurent ear infections. Ann Otol Rhinol Laryngol. 1997, 106 (8): 624-632.

    Article  CAS  PubMed  Google Scholar 

  4. 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.

    Article  CAS  PubMed  Google Scholar 

  5. 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.

    Article  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 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-

    Article  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. 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.

    Article  CAS  PubMed  Google Scholar 

  10. 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.

    Article  CAS  PubMed  Google Scholar 

  11. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  CAS  PubMed  Google Scholar 

  15. 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.

    Article  CAS  PubMed  Google Scholar 

  16. 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 

  17. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. 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.

    Article  CAS  PubMed  Google Scholar 

  20. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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.

    Article  CAS  PubMed  Google Scholar 

  22. R: A Language and Environment for Statistical Computing. http://www.R-project.org,

  23. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. 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.

    CAS  PubMed  Google Scholar 

  25. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Grabe N: AliBaba2: context specific identification of transcription factor binding sites. In Silico Biol. 2002, 2 (1): S1-S15.

    PubMed  Google Scholar 

  30. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. RepeatMasker Open-3.0. http://www.repeatmasker.org,

  32. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 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.

    Article  CAS  PubMed  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. 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.

    Article  CAS  PubMed  Google Scholar 

  37. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. 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-

    PubMed  Google Scholar 

  39. 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.

    Article  PubMed  Google Scholar 

  40. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. 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.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  PubMed  Google Scholar 

  43. 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.

    Article  CAS  PubMed  Google Scholar 

  44. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. 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.

    Article  CAS  PubMed  Google Scholar 

  46. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Thierry-Mieg D, Thierry-Mieg J: AceView: a comprehensive cDNA-supported gene and transcripts annotation. Genome Biol. 2006, 7 (Supp 1): S12-

    Article  PubMed  PubMed Central  Google Scholar 

  51. 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.

    Article  CAS  PubMed  Google Scholar 

  52. 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.

    Article  CAS  PubMed  Google Scholar 

  53. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. 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.

    CAS  PubMed  Google Scholar 

  56. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. 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.

    Article  CAS  PubMed  Google Scholar 

  58. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  59. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgements

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).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Marie S Rye or Sarra E Jamieson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: SEJ JMB MSR. Clinical care and characterization for WAFSOM: SV HLC. WAFSOM recruitment and preparation: ESHS SEJ MSR. Raine Study management: CEP. Adenoid/Tonsil collection and preparation: RT SV HLC. WAFSOM genotyping, cell culture, RNA preparation and RT-PCR expression: MSR. Design and maintenance of in-house genetic database for WAFSOM and bioinformatics CNS pipeline: RWF. Analyzed the data: MSR SEJ JMB. Supervised the work: SEJ JMB CEP. Wrote the paper: MSR SEJ JMB. Reviewed the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12881_2012_1175_MOESM1_ESM.pdf

Additional file 1: Table S1: Microsatellite markers and their respective primer sequences used for replication of previously reported linkage regions. Table S2. Forward and reverse primer sequences for qualitative RT-PCR expression of four chromosome 10 candidate genes. Table S3. Results for SNPs P < 10-3 in the Raine Study across the region of chromosome 10 from 127 Mb to qter. Table S4. Raine Study results for the 10q26.3 region showing respective top SNPs for various genes/regions after adjusting for known environmental covariates using ProbABEL. Table S5. Details of transcription factor binding sites for top SNPs across the TCERG1L/PPP2R2D gene region. A transcription factor is shown only when a change in allele adds or removes a transcription factor binding site. (PDF 223 KB)

12881_2012_1175_MOESM2_ESM.pdf

Additional file 2: Figure S1: LocusZoom plot of association (Figure 1) focused on the TCERG1L/PPP2R2D intergenic region. Pairwise values of LD with the top SNP (rs7922424; in purple) were calculated using the HapMap CEU population. Figure S2. Conserved non-coding sequence plot of the TCERG1L/PPP2R2D intergenic region containing the 9 associated SNPs (P < 0.0001). The top SNP (rs7922424) is annotated in orange, with SNPs P < 0.0001 in green and additional SNPs (P < 0.001) annotated in purple. Colours indicate the following regions; repeats (blue), genes (brown) and coding sequence (red). The three SNPs of interest are highlighted by bold text. Species are human (H), mouse (M), and rat (R). Figure S3. Pairwise LD between the 9 associated SNPs (P < 0.0001) in the Raine Study cohort, where (A) is the D’ value and (B) is the r2 value. For both D’ and r2, LD measures are indicated at the matrix intercept between two markers and are indicative of a decimal value. A square with no value indicated at the intercept equates to a value of 1.00 (or complete LD). For D’ measures, red and shades of red indicate a higher degree of confidence (i.e. LOD ≥2.0). (PDF 349 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Rye, M.S., Scaman, E.S., Thornton, R.B. et al. Genetic and functional evidence for a locus controlling otitis media at chromosome 10q26.3. BMC Med Genet 15, 18 (2014). https://doi.org/10.1186/1471-2350-15-18

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2350-15-18

Keywords