Abstract

Background. Age-related hearing impairment (ARHI) is considered an unpreventable disorder. We aimed to detect specific genetic variants that are potentially related to ARHI via genome-wide association study (GWAS). Methods. A sample of 131 dizygotic twins was genotyped for single-nucleotide polymorphism- (SNP-) based GWAS. Gene-based test was performed using VEGAS2. Pathway enrichment analysis was conducted by PASCAL. Results. The twins are with a median age of 49 years, of which 128 were females and 134 were males. rs6633657 was the only SNP that reached the genome-wide significance level for better ear hearing level (BEHL) at 2.0 kHz (). Totally, 9, 10, 42, 7, 17, and 5 SNPs were suggestive evidence level for () BEHLs at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz and pure tone average (PTA), respectively. Several promising genetic regions in chromosomes (near the C20orf196, AQPEP, UBQLN3, OR51B5, OR51I2, OR52D1, GLTP, GIT2, and PARK2) nominally associated with ARHI were identified. Gene-based analysis revealed 165, 173, 77, 178, 170, and 145 genes nominally associated with BEHLs at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz and PTA, respectively (). For BEHLs at 0.5, 1.0, and 2.0 kHz, the main enriched pathways were phosphatidylinositol signaling system, regulation of ornithine decarboxylase, eukaryotic translation initiation factor (EIF) pathway, amine compound solute carrier (SLC) transporters, synthesis of phosphoinositides (PIPS) at the plasma membrane, and phosphatidylinositols (PI) metabolism. Conclusions. The genetic variations reported herein are significantly involved in functional genes and regulatory domains that mediate ARHI pathogenesis. These findings provide clues for the further unraveling of the molecular physiology of hearing functions and identifying novel diagnostic biomarkers and therapeutic targets of ARHI.

1. Introduction

Hearing impairment is the most prevalent sensory deficit, affecting over 50% of middle-aged people and the elderly in China [1, 2]. Age-related hearing impairment (ARHI) or presbycusis is the most common type of sensorineural hearing loss caused by the natural aging of the auditory system. It is considered an unpreventable and incurable disorder [3]. The typical characteristics of ARHI are gradual progression in later life and bilaterally symmetrical sensorineural hearing impairment, which starts at high frequencies in the early stages and then extends to medium and low frequencies over time. However, early-stage ARHI is often underrecognized. ARHI is a complex, multifactorial disease that is attributable to confounding genetic and environmental factors [4, 5]. The genetic predisposition to hearing impairment variations approximately have accounted for 25%-75% [511]. Genome-wide association study (GWAS) has had an enormous impact on our understanding of the molecular physiology of hearing impairment and has allowed the identification of several genetic loci located at or near the GRM7 [12, 13], DCLK1, PTPRD, GRM8, CMIP, ISG20, ACAN, and TRIOBP genes [14].

Only a small part of the genetic variants is explained by known genetic variation and many potential genes to be further discovered. As modern medicine cannot cure ARHI, the active prevention of it is particularly important. The twin study design relies on study twins raising in the same family environment. Simultaneously, on average, dizygotic (DZ) twins share 50% of the same genes, which can not only be regarded as ordinary sibling pairs but also have perfectly matched ages, prenatal intrauterine environment, and very similar life environment. Therefore, association analysis using DZ twins is more conducive for interpretation of the results. However, the molecular physiology of hearing impairment in middle- and old-aged Chinese population have not been investigated via GWAS yet. This undertaking is important because this population differs from other ethnic populations worldwide in terms of genetic constitutions and lifestyles.

Investigations into genetically related individuals, such as twins, will enhance genetic association studies, and the use of twin-based designs can efficiently identify both common and rare genetic variants underlying complex traits or diseases [15]. A previous study explored the magnitude of genetic impact on better ear hearing levels (BEHLs) and variations in pure tone average (PTA) via twin modelling analyses. Results indicated that heritability estimates range from 47.08% to 54.20% for BEHLs within 2.0-12.5 kHz [16]. Owing to the lack of studies on ARHI among middle- and old-aged Chinese twins via GWAS, we further conducted a GWAS to detect the specific genetic variants potentially associated with ARHI. We expect to be able to identify genetic mutations associated with ARHI and elucidate biological processes.

2. Materials and Methods

Twin Samples CollectionSamples of twins were collected from the latest genetic epidemiology survey (2012-2013) on previously described aging phenotypes [1618]. In brief, information was collected via questionnaires and health examination, including anthropometric and laboratory measurements by well-trained clinicians. Participants were excluded if they were unconscious; unable or unwilling to participate; suffering from heart failure, kidney failure, cancer, or severe mental disorders; and currently pregnant or breast feeding; incomplete cotwin pairs were also dropped. Zygosity was determined using 16 multiple short-tandem sequence repeat DNA markers [1921]. Finally, the samples consisted of 131 complete DZ twins with a median age of 49 years (95% range: 41–67 years), of which 128 were females and 134 were males.

This study was approved by the Regional Ethics Committee of the Qingdao CDC Institutional Review Boards. Prior written informed consent was obtained from all participants. The ethical principles of Helsinki Declaration were followed.

2.1. Audiometric Examination

Audiometric examination was performed following the method described in a previous study [16]. In brief, the twins underwent otoscopy, and then, the pure-tone air-conducted hearing thresholds in each ear were separately measured at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz by using a diagnostic audiometer. BEHL was then calculated as the lower value of both ears at each frequency. Finally, the PTA at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz was separately calculated for the left and the right ear, and the better ear (i.e., the one with the lower value) was selected.

2.2. Genotyping and Quality Control

Genomic DNA was first extracted from the whole peripheral blood of the 131 DZ twins by using QIAamp DNA Blood Mini Kit (Qiagen GmbH, Hilden, Germany). Quantity and integrity of genomic DNA were then determined. Subsequently, DNA samples were genotyped on the Illumina’s Infinium Omni2.5Exome-8v1.2 Bead Chip platform (Illumina, San Diego, CA, USA). Autosome and chromosome X data were analyzed. Quality control was applied using the following criteria: , , , and according to genome-wide efficient mixed-model association (GEMMA) [22]. Linear-Mixed Models were used to test the genotype-phenotype association by using GEMMA. Genetic relationship matrix was included in the model analyses because of our twin pedigree data. Finally, a total of 1,365,315 single-nucleotide polymorphisms (SNPs) qualified for subsequent analyses.

2.3. Statistical Analysis
2.3.1. Basic Characteristics Analysis

Descriptive statistics were computed using SPSS version 22.0. Square-root transformation for BEHLs and rank transformation for PTA were performed for normality. We first performed a normality test for basic characteristics. For those that did not conform to the normal distribution, the Mann-Whitney test was used for comparison.

2.3.2. SNP-Based Analysis

The association between ARHI and SNP genotypes across the genome was tested using the GEMMA software [22]. Sex, age, educational level, and the first five principal components served as covariates in model fitting. SNPs that reached a suggestive evidence level () rather than the conventional genome-wide significance level () for the association were detected [23, 24]. The chromosome X-wide association study (XWAS) was used to find the possible trait association signals from chromosome X. Functional elaboration of the detected SNPs was further performed, and likely, cell types of action were predicted using the HaploReg v4.1 software [25, 26]. Enrichment results of cell type enhancers were reported (uncorrected ).

2.3.3. Gene-Based Analysis

Gene-based analysis was implemented using SNP-set association test via the versatile gene-based association study-2 (VEGAS2) approach, which incorporated information from a full set of GWAS summary data within one gene and accounts for linkage disequilibrium between them [27, 28]. SNPs from “1000G East Asian Population” were adopted. was considered as nominal significance level [29].

2.3.4. Pathway Enrichment Analysis

Pathway enrichment analysis was conducted using pathway scoring algorithm (PASCAL) [30, 31]. First, the location of genetic marker SNPs in the genes was determined, and the related scores of all genes in the pathway were calculated. Chi-squared or empirical scores were used to evaluate the pathway enrichment of high-scoring (possibly fused) genes, avoiding any standard binary enrichment test with inherent value threshold. The pathway and its corresponding genes were selected KEGG, Reactome, and Biocarta.

3. Results

3.1. Basic Characteristics

The basic characteristics of the 131 DZ twins were summarized in Table 1. The males showed a higher moderate and high BEHLs (4.0 and 8.0 kHz) and PTA than the females (), whereas no difference was found in terms of low BEHLs (0.5 and 1.0 kHz).

3.2. SNP-Based Analysis

A total of 1,365,315 SNPs genotyped from the current sample were included in the GWAS. The relationships between the observed and expected GWAS values for BEHLs and PTA were illustrated in quantile–quantile () plots (Figure 1). The values of -statistic were close to one (0.9906–1.0110), suggesting no evidence of bias from population stratification or genomic inflation of the test statistics. The slight deviation in the upper right tail from null distribution crudely suggested some form of associations.

As illustrated in Manhattan plots (Figure 2), rs6633657 was the only SNP that reached the genome-wide significance level (). This SNP was located in the intron region of PTCHD1-AS on chromosome 23 for BEHL at 2.0 kHz. Particular for the trait association signals from chromosome X, then we ascertained by using the XWAS. By analyzing the associations of rs6633657, we identified this SNP was associated with BEHL2.0 (Additional file 1). No other SNP reached the genome-wide significance level () for BEHLs at the other frequencies and PTA. However, 9, 10, 42, 7, and 17 SNPs were suggestive of association () for BEHLs at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz, respectively; by comparison, five SNPs were suggestive of association for PTA (Table 2). No consistent SNPs were observed for BEHLs at 0.5, 1.0, 2.0, 4.0, and 8.0 kHz frequencies.

As illustrated by the regional association plots (Figure 3), several chromosomal loci showed nominal association with ARHI. Among these top signals (Table 2), three SNPs () were located at or near the C20orf196 gene on chromosome 20p12.3 for BEHL at 0.5 kHz (Figure 3(a)); four SNPs () at or near the AQPEP gene on chromosome 5q23.1 for BEHL at 1.0 kHz (Figure 3(b)); nine SNPs () at or near the UBQLN3, OR51B5, OR51I2, and OR52D1 genes on chromosome 11p15.4 for BEHL at 2.0 kHz (Figure 3(c)); two SNPs ( and ) at or near the GLTP and GIT2 genes on chromosome 12q24.11 for BEHL at 4.0 kHz (Figure 3(d)); and six SNPs () at or near the PARK2 gene on chromosome 6q26 for BEHL at 8.0 kHz (Figure 3(e)). For PTA, two SNPs ( and ) were positioned within or closest to the GLTP and GIT2 genes on chromosome 12q24.11 (Figure 3(f)).

The primary T helper memory/regulatory cells from peripheral blood was identified for BEHL at 0.5 kHz by using the HaploReg v4.1 software (Additional file 2). The results were compared with meaningful ARHI-associated SNPs previously reported by other GWAS. No evidence of replication was found.

3.3. Gene-Based Analysis

A total of 165, 173, 77, 178, 170, and 145 genes were observed to be nominally associated with BEHLs at 0.5 kHz, 1.0 kHz, 2.0 kHz, 4.0 kHz, and 8.0 kHz and PTA, respectively () (Additional file 3). The top 20 genes for BEHL0.5 ranked by values are listed in Table 3 (BEHL1.0: Additional file 4, BEHL2.0: Additional file 5, BEHL4.0: Additional file 6, BEHL8.0: Additional file 7, PTA: Additional file 8). C20orf196 gene for 0.5 kHz, SLC16A9, UBQLN3, and OR51I2A genes for 2.0 kHz and FAM184A and TBC1D1genes for 8.0 kHz had already been shown in suggestive level SNP-based.

These results were compared with significant ARHI-associated genes previously reported by other GWAS. Two replicable genes were identified, namely, ACAN [14] for BEHL at 4.0 kHz () and CMIP [32] for BEHL at 8.0 kHz (). Although the well-known GRM7 gene [12, 13] was also identified for BEHL at 8.0 kHz and PTA, the values did not reach the nominal significance level.

3.4. Pathway Enrichment Analysis

The top 20 pathways for BEHL0.5 were sorted according to their empirical values in Table 4 (BEHL1.0: Additional file 9, BEHL2.0: Additional file 10, BEHL4.0: Additional file 11, BEHL8.0: Additional file 12, PTA: Additional file 13). For BEHLs at 0.5, 1.0, and 2.0 kHz, the main enriched pathways were phosphatidylinositol signaling system; regulation of ornithine decarboxylase; EIF pathway; amine compound SLC transporters; synthesis of PIPS at the plasma membrane; PI metabolism; O glycan biosynthesis; transport of glucose and other sugars; and transport of bile salts and organic acids, metal ions, and amine compounds. By comparison, the main enriched pathways for BEHLs at 4.0 and 8.0 kHz were cysteine and methionine metabolism and adherens junction.

4. Discussion

We explored the specific genetic variants in 131 DZ twin pairs that underlie ARHI. VEGAS2 analysis suggested that several genes were nominally associated with BEHLs and PTA. Five consistent genes, namely, C20orf196, GALNT9, INPP4B, SEMA7A, and ARID3B, were observed for BEHLs at 0.5, 1.0, and 2.0 kHz. The SEMA7A gene encodes a member of the semaphorin family of proteins that have been found in activated lymphocytes and erythrocytes and which may play a crucial role in immunomodulatory and neuronal processes [33]. Although their functions in ARHI are uncertain, the other genes can also serve as latent candidates for future work. Our comparison of the ARHI-related genes found herein with those reported by previous GWAS obtained two replicable genes, namely, ACAN [14] for BEHL at 4.0 kHz and CMIP [32] for BEHL at 8.0 kHz. Using the Shared Harvard Inner Ear Database, Hoffman et al. found that ACAN is expressed in the auditory tissues of mouse [14]. In several developmental phases of mouse, it is mainly expressed in the cochlea and cysts, inner and outer hair cells of the cochlea, and spiral and vestibular ganglia [3436]. By comparison, CMIP is expressed in the inner ear. Furthermore, Girotto et al. found that this gene is associated with hearing ability at 0.25, 1.0, and 2.0 kHz [32]. In contrast to our findings, a GWAS meta-analysis of ARHI using pure tone audiometry from multiple cohorts reported seven completely different associated loci. This may be explained by the different ethnic and genetic background [37].

Among the enriched ARHI-related pathways, amine compound SLC transporters [38]; phosphatidylinositol signaling system [39]; synthesis of PIPS at the plasma membrane [3942]; transport of glucose and other sugars bile salts and organic acids metal ions and amine compounds [4346]; cysteine and methionine metabolism [47, 48]; and adherens junction [49] have been previously reported to be associated with ARHI. Aside from these pathways, other pathways that may be related to ARHI were found, including EIF pathway, PI metabolism, O glycan biosynthesis, and regulation of ornithine decarboxylase. To the best of our knowledge, regulation of ornithine decarboxylase had not been reported as associated with ARHI. Ornithine decarboxylase is a key enzyme in the process of polyamine anabolism in the human body. Polyamines have various biological functions, such as antioxidation, free radical scavenging, and intracellular calcium regulation, all of which reportedly have an impact on hearing [50, 51]. Accumulating evidence shows that ornithine decarboxylase is associated with disordered cell growth regulation [5254]. Aside from this pathway, PI metabolism had not been reported to be associated with ARHI. With the discovery of the high expression of the TRPM7 gene in the organ of Corti and cochlea, as well as the detection of TRPM4 immunoreactivity in the inner ear, researchers gradually realized that the TRP channel plays an important role in auditory functions [5557]. However, TRP channels require PI metabolism to be activated [58]. Therefore, PI metabolism is also closely related to the production of hearing.

We also measured these variants by BEHLs and PTA via GWAS. Several SNPs were found to be suggestively associated with ARHI. We compared these SNPs with significant ARHI-associated SNPs previously reported by other GWAS [1214, 32, 59]. Nevertheless, we found several promising genetic regions on chromosomes that were nominally associated with ARHI. The association between the genes involved in these promising genetic regions and ARHI could serve as candidates for further research and validation. Furthermore, the enhancer of primary T helper memory/regulatory cells from the peripheral blood for BEHL at 0.5 kHz was found. Genes involved in immunity and apoptosis are probably related to ARHI [60, 61], and the maintenance of systemic immune functions can prevent accelerated ARHI [62]. Hence, T helper memory/regulatory cells may serve as candidate tissues for further investigation of gene expression in animal models.

Investigations into genetically related individuals, such as twins, will enhance genetic association studies, and the use of twin-based designs can efficiently identify both common and rare genetic variants underlying complex traits or diseases. We conducted this GWAS on ARHI in a sample of middle and old-aged Chinese twins, and the utilization of twin-based design will empower genetic association studies and efficiently identify genetic variants underlying ARHI. However, the present GWAS has several limitations. First, owing to the challenges in recruiting and confirming qualified twin participants, we obtained a relatively small sample size. Thus, a GWAS meta-analysis with a larger sample is warranted. In addition, we could not distinguish sensorineural hearing loss from conductive hearing loss because bone conduction test was not performed in this study. Finally, lack of replication of identified signals was performed.

5. Conclusions

We identified lists of SNPs reached the suggestive evidence level and found several promising genetic regions on chromosomes associated with ARHI measured by BEHLs and PTA. And sets of genes nominally associated with ARHI were involved in significant biological pathways potentially related to pathogenesis of auditory development and hearing impairment. Nevertheless, the potential candidate biomarkers of ARHI reported here should merit further verifications.

Data Availability

The SNPs datasets for this study have been deposited in the European Variation Archive (EVA) (Accession No. PRJEB23749).

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

We thank Gu Zhu for his help on installing and applicating the related software and performing the statistical analyses. This work was supported by the grants from the National Natural Science Foundation of China (grant numbers 81703292 and 31371024), the China Postdoctoral Science Foundation (2016M590622), the Qingdao Postdoctoral Application Research Project (2015150), the Qingdao Key Health Discipline Development Fund (201909A09), the Qingdao Outstanding Health Professional Development Fund (20190947), and the Shandong Medical and Health Science and Technology Development Program (2015wso329).

Supplementary Materials

Supplementary 1. Additional file 1: summarised results for chromosome X SNP rs6633657 with trait BEHL 2.0 kHz.

Supplementary 2. Additional file 2: query SNP enhancer summary for BEHLs and PTA.

Supplementary 3. Additional file 3: the genes nominally associated with BEHLs at each frequency and PTA from VEGAS2 gene-based analysis ().

Supplementary 4. Additional file 4: top 20 genes from VEGAS2 gene-based analysis showing the strongest association with BEHL1.0.

Supplementary 5. Additional file 5: top 20 genes from VEGAS2 gene-based analysis showing the strongest association with BEHL2.0.

Supplementary 6. Additional file 6: top 20 genes from VEGAS2 gene-based analysis showing the strongest association with BEHL4.0.

Supplementary 7. Additional file 7: top 20 genes from VEGAS2 gene-based analysis showing the strongest association with BEHL8.0.

Supplementary 8. Additional file 8: top 20 genes from VEGAS2 gene-based analysis showing the strongest association with PTA.

Supplementary 9. Additional file 9: top 20 KEGG, Reactome, and Biocarta pathway results for BEHL1.0 in the typed GWAS data.

Supplementary 10. Additional file 10: top 20 KEGG, Reactome, and Biocarta pathway results for BEHL2.0 in the typed GWAS data.

Supplementary 11. Additional file 11: top 20 KEGG, Reactome, and Biocarta pathway results for BEHL4.0 in the typed GWAS data.

Supplementary 12. Additional file 12: top 20 KEGG, Reactome, and Biocarta pathway results for BEHL8.0 in the typed GWAS data.

Supplementary 13. Additional file 13: top 20 KEGG, Reactome, and Biocarta pathway results for PTA in the typed GWAS data.