Skip to main content

Interactions of genetic variations in FAS, GJB2 and PTPRN2 are associated with noise-induced hearing loss: a case-control study in China

Abstract

Background

This study aimed to screen and validate noise-induced hearing loss (NIHL) associated single nucleotide polymorphisms (SNPs), construct genetic risk prediction models, and evaluate higher-order gene-gene, gene-environment interactions for NIHL in Chinese population.

Methods

First, 83 cases and 83 controls were recruited and 60 candidate SNPs were genotyped. Then SNPs with promising results were validated in another case-control study (153 cases and 252 controls). NIHL-associated SNPs were identified by logistic regression analysis, and a genetic risk model was constructed based on the genetic risk score (GRS), and classification and regression tree (CART) analysis was used to evaluate interactions among gene-gene and gene-environment.

Results

Six SNPs in five genes were significantly associated with NIHL risk (p < 0.05). A positive dose-response relationship was found between GRS values and NIHL risk. CART analysis indicated that strongest interaction was among subjects with age ≥ 45 years and cumulative noise exposure ≥ 95 [dB(A)·years], without personal protective equipment, and carried GJB2 rs3751385 (AA/AB) and FAS rs1468063 (AA/AB) (OR = 10.038, 95% CI = 2.770, 47.792), compared with the referent group. CDH23, FAS, GJB2, PTPRN2 and SIK3 may be NIHL susceptibility genes.

Conclusion

GRS values may be utilized in the evaluation of the cumulative effect of genetic risk for NIHL based on NIHL-associated SNPs. Gene-gene, gene-environment interaction patterns play an important role in the incidence of NIHL.

Peer Review reports

Introduction

The candidate-gene approach isdriven by hypothesis based on the prior knowledge of single nucleotide polymorphisms (SNPs) and gene functions, and has often produced informative but sometimes contraditory results in studies related to hearing loss. In many studies that reported significant correlations, the odds ratios (ORs) of individual variants were < 2 [1, 2]. The low risk induced by individual polymorphism is not surprising, because the occurrence of hearing loss is usually a multi-step, multigenic process, and any single genetic polymorphism is unlikely to have a dramatic impact on hearing loss risk. Thus, a single gene study may have limited power to predict risk. However, it is unkonwn whether genetic predisposition to hearing loss expansion is due to the interaction of multiple susceptibility genes. Multigenic approachs, which evaluated the joint effects of multiple polymorphisms, may enlarge the influence of individual polymorphisms and enhance the predictive ability. Several recent multigenic studies have demonstrated the promising potential of using such a multigenic approach in correlation studies [3, 4]. In this study, we use noise-induced hearing loss (NIHL) as the hearing loss prototype to illustrate our theme.

NIHL is one of the most widespread prevalent occupational disease, and the second most common sensorineural hearing loss following age-related hearing impairment [5]. NIHL is a polygenic disease induced by the interaction of environmental and genetic factors. The known environmental factors such as noise exposure, organic solvent exposure, smoking and drinking are responsible to NIHL [6, 7]. It was reported that individuals represent variable degrees of NIHL susceptibility, even when exposed to the same levels of noise intensity [8]. Human studies have revealed that the heritability of NIHL in twins is approximately 36% [9]. Knockout mice, such as CDH23+/− [10], GPX1−/− [11], PMCA2−/− [12], SOD1−/− [13], were more susceptible to noise than their wild-type mice. Previous studies have discovered that more than one hundred SNPs may be related to the NIHL susceptibility, including heat shock proteins genes, oxidative stress genes, cadherin proteins genes, potassium recycling channel genes, apoptosis signaling genes, calcium ions recycling channel genes, and inflammatory factor genes [8, 14]. However, there are countless examples in which correlation studies have failed to replicate positive candidate-gene findings, which may be due to small sample size, insufficient statistical methods, and failure to assess the joint effect of multiple pathophysiological associated genes. Therefore, these findings are not enough to explain the total heritability of NIHL.

In present study, we used a multigenic approach to estimate the relationships of gene polymorphisms with NIHL risk. We performed a two-stage case-control study to screen and validate NIHL-associated SNPs among the Chinese populations. As far as we know, this is the most comprehensive multigenic NIHL correlation study reported. We also built a genetic risk score (GRS) model based on these NIHL-related SNPs to better capture the joint effects of multiple genes. Besides, we applied classification and regression tree (CART) to investigate high-order gene-gene, gene-environment interactions in NIHL susceptibility. Our findings will expand our understanding of the genetic basis of NIHL.

Materials and methods

Study participants

We conducted a two-stage case-control study to screen and validate NIHL-associated SNPs among the Chinese populations. The flow chart of this study is presented in Fig. 1. All participants were enrolled from the people who underwent NIHL diagnosis and occupational physical examination from September 2019 to October 2021 in Guangdong Province Hospital for Occupational Disease Prevention and Treatment and the Sixth People’s Hospital of Dongguan. The noise exposed worker was defined as the worker working in the workplace that exceed the occupational noise exposure limit, that is, eight-hour continuous equivalent dB(A) weighted sound pressure levels (LAeq, 8 h) ≥ 85 dB(A). According to the Diagnostic Standards for Occupational Noise-induced Deafness of the People’s Republic of China (GBZ 49-2014) (China, 2014), participants diagnosed with occupational noise-induced deafness were included in the NIHL group. The hearing evaluation of subject was based on the results of air-conduction pure tone audiometry (PTA). NIHL group was diagnosed by the binaural high frequency threshold average (BHFTA) (3000 Hz, 4000 Hz, 6000 Hz), the monaural threshold of weighted value (MTWV) calculated by 90% low frequency (500 Hz, 1000 Hz, 2000 Hz) hearing threshold mean weighting and 10% high frequency (4000 Hz) hearing threshold. NIHL was defined as follows: the noise exposed workers with normal hearing before noise exposure, > 3 years of occupational noise exposure, the results of pure-tone air-conduction threshold pattern assessment were consistent with the trend of noise hearing impairment pattern characteristics, BHFTA were ≥ 40 dB, and MTWV were ≥ 26 dB. The inclusion criteria for the controls were as follows: > 3 years of occupational noise exposure, BHFTA were < 40 dB and MTWV were < 26 dB. All subjects had no hearing-related complications, ear trauma, otitis media, hereditary deafness and blast deafness, family history of hearing loss, craniocerebral injury and usage of certain drugs or toxins (such as aminoglycoside drugs, cisplatin, chloroquine). All the methods and procedures carried out in this study were in accordance with the Declaration of Helsinki. And this study was approved by the Science Ethics Committee of Guangdong Pharmaceutical University (20,190,212). All the participants signed the consent forms.

Fig. 1
figure 1

The subjects were recruited from the workers with occupational noise exposure. 83 NIHL subjects and 83 controls were genotyped. Then NIHL-associated SNPs were validated in another case-control study (153 cases and 252 controls), and GRS model and CART analysis were constructed

Data collection

All subjects were interviewed face to face by a well-trained public health practitioner using a structured questionnaire that included general characteristics (sex, age, and ethnicity), factors related to noise exposure (factory, work situation, daily noise exposure duration, work years with noise exposure, use of personal protective equipment (PPE), exposure to other occupational harmful factors), and lifestyle factors (family history of hearing loss, medical history, medication use, smoking and drinking status).

Other information, including the results of hearing test and environmental noise measurement, was obtained from personal occupational health examination reports and occupational environment monitoring files. Hearing test was conducted by well-trained physicians according to standard procedures. After being out of the noise environment for at least 48 h, all participants underwent a PTA in a sound-insulation room with background noise less than 25 dB(A). The hearing threshold for each ear was measured at six frequencies, and the deviation of hearing threshold of each frequency was ≤ 10dB. Cumulative noise exposure (CNE) was calculated by years of noise exposure and sound pressure level (SPL) to assess individual noise exposure dosage. Based on the technical requirements of Workplace Physical Factors Measurement Part 8: Noise (GBZ/T189.8-2007) (China, 2007), environmental noise exposure level was assessed using LAeq,8 h, the equivalent A-weighted sound pressure level of the contact noise intensity normalized to 8 h during the actual working time of a day. CNE is an evaluation parameter derived from the theory of equal energy exposure and equal biological effect, so that individual noise exposure levels of different noise exposure intensities and times have comparability.

Blood collection and DNA extraction

Whole blood samples were collected from each participant after an overnight fasting. DNA was extracted from blood leucocytes using the TIANamp Genomic DNA Kit (centrifugation column method) (Tiangen, Bejing, China) following the manufacturer’s protocol. The integrity of DNA was measured by agarose gel electrophoresis, and the content was measured by a microspectrophotometer (NanoDrop 5000, Thermo Scientific, USA).

SNPs selection and genotyping in the screening stage

83 NIHL cases and 83 controls were recruited. A total of 60 SNPs in 38 genes were selected as candidate SNPs according to the National Center for Biotechnology Information (NCBI) database and previous reports (Supplementary Table 1). The sequences of specific amplification primers and extension primers were designed on Primer Premier software (v6.0, Premier Biosoft Inc., CA). All the SNPs were genotyped on 3730XL gene sequencing instrument (Applied Biosystems Inc., USA) by the SNapShot assays at the Good Biotechnology Co., Ltd. (Guangzhou, China), and were analyzed with Gene Mapper software (v5.0, Premier Biosoft Inc., CA).

SNPs verification in the replication stage

A total of 153 NIHL and 252 controls were selected. All the SNPs screened from the first stage were genotyped using the Taqman-MGB method. The design methods of TaqMan-MGB probe and primer were the same as above. Genomic DNA was added into 2 x TaqMan qPCR Master Mix (Vazyme Biotech Co., Ltd, Nanjing, China) with genotyping probes and primers according to the instructions of manufacturer, and the amplified products were scanned by fluorescence in the whole wavelength range on Real-time fluorescent quantitative PCR instrument (CFX96, Bio-rad, USA). Analyses were performed using BIORAD CFX Maestro software (v3.1, Bio-rad, USA).

Statistical analysis

Continuous variables for the normal distribution were expressed as mean ± standard deviation (SD) and as median (P25, P75) for skewed distribution and analyzed by Student’s t-tests. Whether sample group was representative was decided by Hardy-Weinberg equilibrium (HWE) tests. Categorical variables were expressed as frequency (%) and analyzed by χ2 test. CNE = 10 × log (10SPL * years of noise exposure), where SPL is the sound pressure level [dB(A)] of noise exposure. HWE was tested using Pearson’s χ2 for each SNP in the controls. SNPs were eliminated when the call rates were < 90% and the controls deviated from HWE with p < 0.01. Multiple logistic regression was conducted to calculate the OR and 95% confidence interval (CI) under the dominant, super-dominant, recessive, additive and allele model, and adjusted by confounders, including age, years of occupational noise exposure, smoking and drinking status.

To combine the relatively small effects of individual genes and to better explore the complex relationship between genetics and NIHL disease, GRS model was performed. The SNPs, whose p values < 0.05 in logistic regression model during the replication stage, were selected due to the limited sample size and false negative signals. GRS refers to the sum of risk genotypes (0 for non-risk genotypes, 1 for other genotypes, and 2 for risk genotypes) across NIHL-associated SNPs, and was calculated as reported previously [3]. The participants were divided based on the GRS values and the ORs were calculated with the subjects in the lowest levels of GRS (≤ 5) as a reference. The dose-response relationship between GRS value and OR was assessed using the χ2 test.

CART analysis was applied to establish a decision tree through recursive partitioning to determine specific combinations of genetic and environmental factors related to disease risk [15]. Based on the splitting rule of information index, the data were layered into individual subsets, which are represented as nodes in this decision tree. This process terminated when the classification achieves the lowest cross-validation error in the terminal nodes. Taking the lowest percentage of termination nodes as a reference, logistic regression analysis was conducted to calculate the OR value and its 95% CI of different branches, and adjusted by confounders such as age, years of occupational noise exposure, smoking and drinking status. The learning set was applied to build the tree model, and the testing set was used to internally verify the generated tree model. The data were randomly divided into a learning set (65% of the data) and a testing set (35% of the data).

All statistical analyses were conducted using R software (v 4.1.2). Two-tailed p < 0.05 indicated statistical significance, and 0.05 < p < 0.10 was marginal significance.

Results

Eight candidate SNPs were replicated in the screening stage

A total of 83 NIHL subjects and 83 controls were recruited (Supplementary Table 2). Briefly, the mean (SD) age, years of noise exposure and CNE was 48.32 ± 5.69 years, 12.24 ± 5.75 years and 96.51 ± 5.49 dB(A)·years, respectively. There was no statistical significance between the NIHL cases and controls regarding age, sex, ethnicity, years of noise exposure, CNE, use of PPE, smoking and drinking status (p > 0.05). The binaural high frequency threshold average (BHFTA) in NIHL group was 55.23 ± 10.38 dB, and was significantly higher than the control group (28.02 ± 7.49 dB; p < 0.05).

There were 15 SNPs with marginal or statistical significance by χ2 test. Seven SNPs in six genes (CDH23 rs2394795, FAS rs1468063 and rs2862833, GJB2 rs3751385, HOTAIR rs87494, SIK3 rs6589574, STAT3 rs1053005) were found to be significantly different between two groups under at least one genetic model (p < 0.05; Supplementary Table 3), and eight SNPs in seven genes (FOXO3 rs2802292 and rs10457180, GAPDH rs6489721, HDAC2 rs10499080, hOGG1 rs1052133, PTPRN2 rs10081191, STAT3 rs1053023 and TAB2 rs521845) were marginally significant (0.05 < p < 0.10; Supplementary Table 3). Finally, eight SNPs in six genes (CDH23 rs2394795, FAS rs1468063 and rs2862833, GJB2 rs3751385, PTPRN2 rs10081191, SIK3 rs6589574, and STAT3 rs1053023 and rs1053005) were related to NIHL susceptibility in logistic regression model after adjusting for age, year of noise exposure, smoking and drinking status (p < 0.05; Table 1).

Table 1 Odd ratios (ORs) and 95% CI of NIHL-associated SNPs in the screening stage

Six candidate SNPs were verified in the replication stage

In this stage, 153 NIHL and 252 controls were enrolled. As shown in Supplementary Table 4, the mean age in NIHL group was significantly higher than in control group (p < 0.05). No statistical significance between the cases and controls in terms of sex, ethnicity, CNE and years of noise exposure were observed (p > 0.05). The use rate of PPE and drinking rate in the NIHL group were slightly lower than those in the control group (p < 0.05). The BHFTA in NIHL group was 55.86 ± 10.72 dB, and was significantly higher than the control group (24.92 ± 8.17 dB; p < 0.05).

We further verified the eight candidate SNPs screened from the first stage. Five SNPs in four genes (CDH23 rs2394795, FAS rs1468063 and rs2862833, GJB2 rs3751385, SIK3 rs6589574) showed significantly different between two groups under at least one genetic model (p < 0.05; Supplementary Table 5). Finally, the variant alleles of CDH23 rs2394795, FAS rs1468063 and rs2862833, GJB2 rs3751385, PTPRN2 rs10081191 and SIK3 rs6589574 showed significantly association with NIHL risk in logistic regression model after adjusting for multiple variables (Table 2).

Table 2 ORs and 95% CI of NIHL-associated SNPs in the replication stage

Genotyping and calculation of GRS

The remaining six SNPs mentioned above were used to build the GRS model. Different distributions of GRS values were presented between NIHL cases and the controls (Table 3). The median (P25, P75) of GRS in the NIHL group was 7.0 (5.0, 8.0), and was significantly higher than the control group [5.0 (4.0, 7.0)] (p < 0.001). The proportion of NIHL subjects was gradually higher than controls as the increase of GRS value. When the individuals were divided based on GRS values, the lowest levels of GRS (≤ 5) as the reference group, the positive trend and significant association of the GRS values with NIHL risk was observed (p < 0.05). The individuals were further divided into high and low GRS group according to a GRS value of 6, then the subjects with GRS ≥ 6 had increased risk for NIHL, with an OR of 2.734 (95%CI = 1.801, 4.195) compared to the subjects with GRS < 6.

Table 3 Associations between levels of GRS and risk for NIHL.

CART analysis

CART analysis was performed through incorporating both the genetic and other variables. Compared to individuals carried GJB2 rs3751385 (AA/BB) variant genotypes, individuals carried GJB2 rs3751385 (BB), FAS rs1468063 (AA/BB) and PTPRN2 rs10081191(AA) variant genotypes had higher risk of NIHL (OR = 4.250; 95% CI = 1.984, 9.539; p < 0.001), and individuals carried GJB2 rs3751385 (BB) and FAS rs1468063 (AA/BB) variant genotypes had the highest risk of NIHL (OR = 12.084; 95% CI = 3.923, 46.526; p < 0.001; Table 4). Table 5 depicted the resultant tree structure generated. There was an initial split on GJB2 rs3751385, suggesting that GJB2 rs3751385 is the most important risk factor for NIHL. Further analysis of the CART structure found distinct patterns for GJB2 rs3751385 (AA/BB) and GJB2 rs3751385 (BB). Compared to the reference group (individuals with age < 45 years and carried GJB2 rs3751385 (AA/BB)), individuals carried GJB2 rs3751385 (BB) and PTPRN2 rs10081191 (AA) variant genotypes had the higher risk of NIHL (OR = 7.585; 95%CI = 3.654, 16.443; p < 0.001; terminal node 8), and individuals carried GJB2 rs3751385 (BB) and PTPRN2 rs10081191 (BB/AB) had a 2.483-fold increased risk of NIHL (95% CI = 1.239, 4.998; p < 0.05; terminal node 7). For age ≥ 45, individuals carried GJB2 rs3751385 (AA/AB) and FAS rs1468063 (BB) in variant genotypes had a 3.012-fold increased risk of NIHL (95%CI = 1.389, 6.567; p < 0.05; terminal node 6), individuals with CNE < 95 (dB(A)·years), without PPE and carried GJB2 rs3751385 (AA/AB) and FAS rs1468063 (AA/AB) variant genotypes had a 10.038 -fold increased risk of NIHL (95%CI = 2.770, 47.792; p < 0.05; terminal node 5). For non-drinkers, individuals with age ≥ 45 years and CNE < 95 (dB(A)·years) carried GJB2 rs3751385 (AA/AB) and FAS rs1468063 (AA/AB) variant genotypes had a 2.677-fold increased risk of NIHL (95% CI = 1.103, 6.447; p < 0.05; terminal node 3).

Table 4 Risk estimate of gene-gene interaction in CART model
Table 5 Risk estimate of gene-environment interaction in CART model

Discussion

In present study, we have applied a multigenic approach to systematically explore the relationships between polymorphisms in susceptibility genes and NIHL risk. We screened and validated 6 NIHL-associated SNPs in five genes from 60 candidate-SNPs in a two-stage case-control study. The most important discovery in this study is that joint analyses of multiple SNPs in multiple susceptibility genes may find otherwise undetectable correlations between a single SNP and NIHL risk. Moreover, we contructed a genetic risk predictive model according to these SNPs and found a dose-response relationship between GRS values and NIHL risk. Our findings indicate that a more comprehensive multigenic approach combining multiple polymorphisms provides more accurate delineation of risk groups and may suggest the future direction of correlation studies.

Many studies have reported the relationships between genetic polymorphisms and susceptibility to complex diseases like NIHL. With the progress of science and technology, more and more genomic methods were applied to explore the genetic susceptibility of NIHL [16], such as MALDI-TOF-MS [17], TaqMan probe method [18], Sanger sequencing [19], SNaPshot sequencing [20], exome-wide association study (EWAS) [21] and genome wide association study (GWAS) [22, 23]. In this study, SNaPshot sequencing was used in the screening stage and Taqman-MGB probe method was conducted in the replication stage for genotyping. Most studies have conducted a candidate-gene approach to explore only one or several selected genes at a time. Nevertheless, a large number of positive candidate genes have not been replicated in different populations. Replication in independent sample sets plays a vital role in confirming susceptibility genes for complex diseases. It is now generally believed that replication of findings in several independent sample sets is more important than to obtain highly significant p values [24]. Therefore, only findings that caused significant associations in both sample sets were considered meaningfully. Applying this replication criterion may acquire more reliable association results and avoid the problem of multiple testing. In our study, six SNPs in five genes, including SIK3 rs6589574, FAS rs1468063 and rs2862833, GJB2 rs3751385, PTPRN2 rs10081191 and CDH23 rs2394795, showed significant (p < 0.05) association with NIHL risk.

NIHL is a polygenic disease involving variants of multiple SNPs in multiple genes, thus the impact of a single polymorphism locus is weak. GRS was used to evaluate the cumulative effect of the genetic markers mentioned above on individual genetic predisposition to NIHL, similar to that conducted on other specific complex diseases, such as cancer and diabetes [25, 26]. Our data indicated that NIHL cases have a higher genetic susceptibility than controls, and these SNPs have a cumulative effect on NIHL genetic risk. In addition, a positive dose–response relationship between GRS levels and NIHL risk was observed. These findings extend our understanding of the relationship between multiple genes and the susceptibility of NIHL. Using these genetic biomarkers, we could screen individuals susceptible to NIHL, and distinguish higher sensitivity to NIHL from the noise-exposed workers. It is relatively difficult to avoid noise exposure in most work environments. Thus, effective and efficient preventive measures for high-risk populations are very important. In this case, the susceptible individuals should be screened and identified, and appropriate measures can be taken to reduce noise exposure and strengthen protection (wearing earplugs or earmuffs) in the noise environment, so as to effectively reduce the risk of NIHL.

Since environment factors such as sex, age, ethnicity and lifestyle may interact with genetic factors to promote the occurrence of NIHL, the impact of each individual SNP is unlikely to be substantial [27, 28]. Therefore, CART was conducted to analysis higher-order interactions among gene-gene, gene-environment. Compared to individuals carried GJB2 rs3751385 (AA/BB) variant genotypes, individuals carried GJB2 rs3751385 (BB), FAS rs1468063 (AA/BB) and PTPRN2 rs10081191 (AA) variant genotypes had a 4.250-fold increased NIHL risk (95% CI = 1.984, 9.539), and individuals carried GJB2 rs3751385 (BB) and FAS rs1468063 (AA/BB) variant genotypes had a 12.084-fold increased NIHL risk (95% CI = 3.923, 46.526). We also observed a significant multiplicative interaction between age and three SNPs. All genetic effects were apparent only in “age ≥ 45” (subjects aged ≥ 45 years ) and not in “age < 45”. More interestingly, subgroups of individuals with higher NIHL risks were found according to simple combinations of age and genotypes. Our data presented an interaction between rs3751385, rs1468063, rs10081191, age, CNE, PPE and drinking status. Aging can aggravate the apoptosis of hair cells in cochlear, especially for the elderly exposed to noise, which may increase the prevalence of NIHL. It was reported that noise exposed workers aged 45–59 have a higher risk of NIHL than those aged 30–44 [29]. In addition, an obvious dose-response relationship between CNE and hearing loss has been confirmed [30]. It has been confirmed that drinking can increase the risk of NIHL, which may be related to vasospasm and contraction of inner ear terminals caused by excessive drinking, as well as ischemia and hypoxia of cochlear cells [31, 32]. It was well known that the easiest method to prevent NIHL is to wear earmuffs or earplugs. A study showed that in the logistic regression analysis, there was a significant difference in the rates of hearing loss between wearing earplugs sometimes (OR = 1.48, 95% CI = 1.07, 2.05) and never wearing earplugs (OR = 1.53, 95% CI = 1.12, 2.10) [33]. CART algorithm may rapidly identify the potential genetic and environmental interactions when dealing with numerous variables in complex diseases. However, the CART analysis is a postdata-mining tool and the results should be interpreted more cautious.

Briefly, six SNPs in five genes showed significant association with risk of NIHL in our study. Interestingly, all of them were non-coding SNPs, which may cause single base changes, resulting in 3’-UTR or intron variants, ultimately altering gene expression levels rather than proteins. The current studies on disease-related SNPs mainly focus on the regulatory regions or coding regions of the genome. However, in the human genome, SNPs in non-coding regions are more than in coding regions. The 3’ untranslated region (3’UTR) is crucial for the genes transcriptional regulation [34]. Gene variation in 3’UTR does not directly change the coding sequence of the gene, but can cause the change of gene function. MiRNAs are required to regulate gene expression by binding to the 3’UTR of target mRNAs. SNPs in 3’UTR can affect the recognition of miRNAs and target genes and the expression level of mature miRNAs, so as to enhance the expression of target gene by weakening the translation inhibitory effect of miRNA on target gene, or inhibit the expression of target gene by forming new miRNA binding sites. A recent study based on a large sample noise exposure cohort has reported that a variation in non-coding regions of gene adaptor-associated kinase 1 (AAK1) associated with the susceptibility of NIHL. Individuals with T/T genotype of AAK1 rs1396793 had stronger resistance to NIHL than those with G/T or G/G genotype, which may related to the upregulation of expression of AAK1. Compared with the AAK1G/G and AAK1T/G mice, the AAK1T/T mice also showed significantly less auditory brainstem response (ABR) threshold elevation and wave I amplitude decrement at higher frequencies after noise exposure, which may provide a new target for the prevention and treatment of NIHL [21]. Our data suggest an interaction between GJB2, FAS and PTPRN2. We speculate that one or more transcriptional regulators with similar functions and structures may simultaneously regulate these four genes expression levels, which together act to enhance the risk of NIHL. Therefore, these findings suggested that the detection of non-coding SNPs is also of great significance for the study of disease-related SNPs.

GJB2 protein encoded by GJB2 (Cx26) has about 220 mutation SNPs associated with hearing loss. The mutations of GJB2 were closely related to late-onset progressive hearing loss, especially among East Asia populations [35]. GJB2, as a potassium circulation channel gene, is closely related to NIHL and has been validated in various populations [3, 36]. GJB2 knock-in mice were more prone to NIHL, which may be due to the reduction of cochlear amplifier caused by lowered endocochlear potential, or the excitotoxicity of inner hair cells induced by potassium accumulation around hair cells [37].

The FAS gene belongs to the tumor necrosis factor receptor superfamily, and plays an essential role in the physiological regulation of programmed cell death [38]. FAS gene polymorphisms have been mainly studied in association with malignancies and immune system diseases [39, 40]. A previous study confirmed that the genetic polymorphisms in the FAS gene are associated with NIHL risk [41]. Noise can induce large amounts of reactive oxygen species (ROS) in cochlear tissue and blood, leading to oxidative stress, lipid peroxidation, and DNA damage, and further activating the FAS gene. The mutation of 3’-UTR of FAS resulted in the failure to suppress the high expression of FAS under noise stress, which eventually led to the apoptosis of hair cells.

PTPRN2 is a member of the Protein Tyrosine Phosphatases family (PTPs). PTPRN2 is an autoantigen in insulin-dependent diabetes and has been mainly studies in relation to metabolic diseases such as obesity, diabetes, and cancer [42, 43]. PTPRN2 has been shown to be linked to severe bilateral hearing loss [44, 45]. A Chinese population study showed an increased prevalence of NIHL in individuals with PTPRN2-rs10081191 allele A [23]. Recently, an increasing number of PTPs is identified as clinically relevant targets [45]. Thus larger population samples and more in-depth experimental studies are needed to confirm the mechanism of PTPRN2 in hearing loss.

CDH23 encodes a member of the cadherin super-family and has been validated associations with NIHL risk [1, 3]. The CDH23 protein is strongly associated with the structure and function of the inner ear and is essential for cellular adhesion. The mutation of CDH23 changes the cellular adhesion in inner hair-cells, reduces the stability of cadherin and increases potassium (K+) influx and frequent depolarization, which may cause the opening of the calcium ion channel on the basolateral membrane and the entery of excessive calcium ions, resulting in hair-cells death and hearing loss [46, 47].

SIK3 is a subfamily of serine/threonine protein kinase and is widely involved in glucose, cholesterol and lipid metabolism. SIK3 plays an important role in the early development of hair cells or sterecilia, and plays a maintenance function of non-neural cells throughout development and adulthood of the spiral ganglion [26, 41]. A meta-analysis of the whole genome related to hearing function from G-EAR Association and Twins UK found that rs681524 locus of SIK3 gene is closely related to hearing function of European population [48]. Another case-control study suggested that 3 SNPs (rs493134, rs6589574, and rs7121898) of SIK3 may be an important part of NIHL susceptibility [49].

This study had several advantages that could improve the robustness of our findings. First, two independent sample sets were used to confirm the true association between NIHL susceptibility and genetic polymorphisms, which may improve the reproducibility and power of the study. Second, GRS was employed as an risk prediction model to screen subjects sensitive to NIHL, so as to decrease the risk of NIHL. In addition, a multigenic approach such as CART model was used to estimate the gene-gene, gene-enviromnent interaction to more comprehensively understand the risk of NIHL. However, the shortcoming of this study is that the sample size remains limited. Further studies are warranted to confirm our findings.

Data availability

The data reported in this paper have been deposited in the OMIX, China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences (https://ngdc.cncb.ac.cn/omix), and the OMIX ID (shared link) of the data of screening stage and replication stage is OMIX005285 (https://ngdc.cncb.ac.cn/omix/preview/QtCW2nhS) and OMIX005284 (https://ngdc.cncb.ac.cn/omix/preview/GvGTuAln), respectively.

References

  1. Kowalski TJ, Pawelczyk M, Rajkowska E, Dudarewicz A, Sliwinska-Kowalska M. Genetic variants of CDH23 associated with noise-induced hearing loss. Otol Neurotol. 2014;35(2):358–65.

    Article  PubMed  Google Scholar 

  2. Zhang X, Liu Y, Zhang L, Yang Z, Shao Y, Jiang C, et al. Genetic variations in protocadherin 15 and their interactions with noise exposure associated with noise-induced hearing loss in Chinese population. Environ Res. 2014;135:247–52.

    Article  CAS  PubMed  Google Scholar 

  3. Zhang X, Ni Y, Liu Y, Zhang L, Zhang M, Fang X, et al. Screening of noise-induced hearing loss (NIHL)-associated SNPs and the assessment of its genetic susceptibility. Environ Health. 2019;18(1):30.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Sharma KL, Rai R, Srivastava A, Sharma A, Misra S, Kumar A, et al. A multigenic approach to evaluate genetic variants of PLCE1, LXRs, MMPs, TIMP, and CYP genes in gallbladder cancer predisposition. Tumour Biol. 2014;35(9):8597–606.

    Article  CAS  PubMed  Google Scholar 

  5. Śliwińska-Kowalska M, Zaborowski K. WHO Environmental noise guidelines for the European Region: a systematic review on environmental noise and permanent hearing loss and Tinnitus. Int J Environ Res Public Health. 2017;14(10).

  6. Rosati R, Jamesdaniel S. Environmental exposures and hearing loss. Int J Environ Res Public Health. 2020;17(13).

  7. Rabinowitz PM, Galusha D, Slade MD, Dixon-Ernst C, O’Neill A, Fiellin M, et al. Organic solvent exposure and hearing loss in a cohort of aluminium workers. Occup Environ Med. 2008;65(4):230–5.

    Article  CAS  PubMed  Google Scholar 

  8. Sliwinska-Kowalska M, Pawelczyk M. Contribution of genetic factors to noise-induced hearing loss: a human studies review. Mutat Res. 2013;752(1):61–5.

    Article  CAS  PubMed  Google Scholar 

  9. Heinonen-Guzejev M, Vuorinen HS, Mussalo-Rauhamaa H, Heikkilä K, Koskenvuo M, Kaprio J. Genetic component of noise sensitivity. Twin Res Hum Genet. 2005;8(3):245–9.

    Article  PubMed  Google Scholar 

  10. Holme RH, Steel KP. Progressive hearing loss and increased susceptibility to noise-induced hearing loss in mice carrying a Cdh23 but not a Myo7a mutation. J Assoc Res Otolaryngol. 2004;5(1):66–79.

    Article  PubMed  Google Scholar 

  11. Ohlemiller KK, McFadden SL, Ding DL, Lear PM, Ho YS. Targeted mutation of the gene for cellular glutathione peroxidase (Gpx1) increases noise-induced hearing loss in mice. J Assoc Res Otolaryngol. 2000;1(3):243–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Kozel PJ, Davis RR, Krieg EF, Shull GE, Erway LC. Deficiency in plasma membrane calcium ATPase isoform 2 increases susceptibility to noise-induced hearing loss in mice. Hear Res. 2002;164(1–2):231–9.

    Article  CAS  PubMed  Google Scholar 

  13. Ohlemiller KK, McFadden SL, Ding DL, Flood DG, Reaume AG, Hoffman EK, et al. Targeted deletion of the cytosolic Cu/Zn-superoxide dismutase gene (Sod1) increases susceptibility to noise-induced hearing loss. Audiol Neurootol. 1999;4(5):237–46.

    Article  CAS  PubMed  Google Scholar 

  14. Miao L, Ji J, Wan L, Zhang J, Yin L, Pu Y. An overview of research trends and genetic polymorphisms for noise-induced hearing loss from 2009 to 2018. Environ Sci Pollut Res Int. 2019;26(34):34754–74.

    Article  PubMed  Google Scholar 

  15. Wu X, Gu J, Grossman HB, Amos CI, Etzel C, Huang M, et al. Bladder cancer predisposition: a multigenic approach to DNA-repair and cell-cycle-control genes. Am J Hum Genet. 2006;78(3):464–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Tsuchihashi Z, Dracopoli NC. Progress in high throughput snp genotyping methods. Pharmacogenomics J. 2002;2(2):103–10.

    Article  CAS  PubMed  Google Scholar 

  17. Tumane RG, Thakkar L, Pingle SK, Jain RK, Jawade AA, Raje DV. Expression of serum proteins in noise induced hearing loss workers of mining based industry. J Proteom. 2021;240:104185.

    Article  CAS  Google Scholar 

  18. Shen H, Dou J, Han L, Bai Y, Li Q, Hong Z, et al. Genetic variation in ape1 gene promoter is associated with noise-induced hearing loss in a Chinese population. Int Arch Occup Environ Health. 2016;89:621–28.

    Article  CAS  PubMed  Google Scholar 

  19. Zhou H, Zhou J, Li H, Hui C, Bi J. Paraoxonase 3 gene polymorphisms are associated with occupational noise-induced deafness: a matched case-control study from China. PLoS ONE. 2020;15(10):e240615.

    Article  Google Scholar 

  20. Wang B, Wan L, Sun P, Zhang L, Han L, Zhang H, et al. Associations of genetic variation in e3 sumo-protein ligase cbx4 with noise-induced hearing loss. Hum Mol Genet. 2022;31(13):2109–20.

    Article  CAS  PubMed  Google Scholar 

  21. Wang Q, Wang X, Yang T, et al. A common variant in AAK1 reduces risk of noise-induced hearing loss. Natl Sci Rev. 2023;10(7):nwad080.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Lavinsky J, Crow AL, Pan C, Wang J, Aaron KA, Ho MK, et al. Genome-wide association study identifies nox3 as a critical gene for susceptibility to noise-induced hearing loss. Plos Genet. 2015;11(4):e1005094.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Niu Y, Xie C, Du Z, Zeng J, Chen H, Jin L, et al. Genome-wide association study identifies 7q11.22 and 7q36.3 associated with noise-induced hearing loss among Chinese population. J Cell Mol Med. 2021;25(1):411–20. https://doi.org/10.1111/jcmm.16094.

    Article  CAS  PubMed  Google Scholar 

  24. Neale BM, Sham PC. The future of association studies: gene-based analysis and replication. Am J Hum Genet. 2004;75(3):353–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Ye Z, Austin E, Schaid DJ, Kullo IJ. A multi-locus genetic risk score for abdominal aortic aneurysm. Atherosclerosis. 2016;246:274–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Go MJ, Lee Y, Park S, Kwak SH, Kim BJ, Lee J. Genetic-risk assessment of GWAS-derived susceptibility loci for type 2 diabetes in a 10 year follow-up of a population-based cohort study. J Hum Genet. 2016;61(12):1009–12.

    Article  CAS  PubMed  Google Scholar 

  27. Zhou L, Ruan X, Wang T, Xie H, Hu Y, Shi Z, Xin J, Zhou J, Xue P, Wei F, Zhang Y, Zhang M, Zou H. Epidemiological characteristics of hearing loss associated with noise temporal structure among manufacturing workers. Front Integr Neurosci. 2022;16:978213.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Lie A, Skogstad M, Johannessen HA, Tynes T, Mehlum IS, Nordby KC, Engdahl B, Tambs K. Occupational noise exposure and hearing: a systematic review. Int Arch Occup Environ Health. 2016;89(3):351–72.

    Article  PubMed  Google Scholar 

  29. Grobler LM, Swanepoel W, Strauss S, Becker P, Eloff Z. Occupational noise and age: a longitudinal study of hearing sensitivity as a function of noise exposure and age in South African gold mine workers. S Afr J Commun Disord. 2020;67(2):e1–e7.

    PubMed  Google Scholar 

  30. Xie HW, Zhou JN, Zhou LF, Gao XJ, Quan CJ, Luan YQ, et al. Dose-response relationship between noise kurtosis and noise-induced hearing loss. Zhonghua Lao Dong Wei Sheng Zhi Ye Bing Za Zhi. 2021;39(7):493–7.

    CAS  PubMed  Google Scholar 

  31. Gopinath B, Flood VM, McMahon CM, Burlutsky G, Smith W, Mitchell P. The effects of smoking and alcohol consumption on age-related hearing loss: the Blue mountains hearing study. Ear Hear. 2010;31(2):277–82.

    Article  CAS  PubMed  Google Scholar 

  32. Lao XQ, Yu IT, Au DK, Chiu YL, Wong CC, Wong TW. Noise exposure and hearing impairment among Chinese restaurant workers and entertainment employees in Hong Kong. PLoS ONE. 2013;8(8):e70674.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Kim HJ, Oh SY, Won SY, Kim HJ, Kim TK, Ko BC, et al. Associations between earplug use and hearing loss in ROK military personnel. BMJ Mil Health. 2021;167(6):398–401.

    Article  PubMed  Google Scholar 

  34. Marzluff WF, Wagner EJ, Duronio RJ. Metabolism and regulation of canonical histone mRNAs: life without a poly(A) tail. Nat Rev Genet. 2008;9(11):843–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Xiao ZA, Xie DH. GJB2 (Cx26) gene mutations in Chinese patients with congenital sensorineural deafness and a report of one novel mutation. Chin Med J (Engl). 2004;117(12):1797–801.

    CAS  PubMed  Google Scholar 

  36. Pawelczyk M, Van Laer L, Fransen E, Rajkowska E, Konings A, Carlsson PI, et al. Analysis of gene polymorphisms associated with K ion circulation in the inner ear of patients susceptible and resistant to noise-induced hearing loss. Ann Hum Genet. 2009;73(Pt 4):411–21.

    Article  CAS  PubMed  Google Scholar 

  37. Lin X, Li G, Zhang Y, Zhao J, Lu J, Gao Y, et al. Hearing consequences in Gjb2 knock-in mice: implications for human p.V37I mutation. Aging. 2019;11(18):7416–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Op de Beeck K, Schacht J, Van Camp G. Apoptosis in acquired and genetic hearing impairment: the programmed death of the hair cell. Hear Res. 2011;281(1–2):18–27.

    Article  PubMed  Google Scholar 

  39. Jalali SA, Homaei Shandiz F, Tavakol Afshari J, Davarpanah Tanha Ghochan M, Nikpoor AR, Mohammadi M. Status of FAS and FAS Ligand Gene Polymorphisms in patients with breast Cancer in Northeastern IRAN. Rep Biochem Mol Biol. 2018;7(1):23–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Lu MM, Ye QL, Feng CC, Yang J, Zhang T, Li J, et al. Association of FAS gene polymorphisms with systemic lupus erythematosus: a case-control study and meta-analysis. Exp Ther Med. 2012;4(3):497–502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Xu S, Wang B, Han L, Pu Y, Zhu B, Zhang J. Polymorphisms in the FAS gene are associated with susceptibility to noise-induced hearing loss. Environ Sci Pollut Res Int. 2021;28(17):21754–65.

    Article  CAS  PubMed  Google Scholar 

  42. van den Maagdenberg AM, Schepens JT, Schepens MT, Pepers B, Wieringa B, van Kessel AG, et al. Assignment of Ptprn2, the gene encoding receptor-type protein tyrosine phosphatase IA-2beta, a major autoantigen in insulin-dependent diabetes mellitus, to mouse chromosome region 12F. Cytogenet Cell Genet. 1998;82(3–4):153–5.

    Article  PubMed  Google Scholar 

  43. Lee S. The association of genetically controlled CpG methylation (cg158269415) of protein tyrosine phosphatase, receptor type N2 (PTPRN2) with childhood obesity. Sci Rep. 2019;9(1):4855.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Van Laer L, Huyghe JR, Hannula S, Van Eyken E, Stephan DA, Mäki-Torkko E, et al. A genome-wide association study for age-related hearing impairment in the Saami. Eur J Hum Genet. 2010;18(6):685–93.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Stanford SM, Bottini N. Targeting tyrosine phosphatases: time to end the Stigma. Trends Pharmacol Sci. 2017;38(6):524–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Jaiganesh A, Narui Y, Araya-Secchi R, Sotomayor M. Beyond cell-cell adhesion: sensational cadherins for hearing and balance. Cold Spring Harb Perspect Biol. 2018;10(9).

  47. Elledge HM, Kazmierczak P, Clark P, Joseph JS, Kolatkar A, Kuhn P, et al. Structure of the N terminus of cadherin 23 reveals a new adhesion mechanism for a subset of cadherin superfamily members. Proc Natl Acad Sci U S A. 2010;107(23):10708–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Wolber LE, Girotto G, Buniello A, Vuckovic D, Pirastu N, Lorente-Cánovas B, et al. Salt-inducible kinase 3, SIK3, is a new gene associated with hearing. Hum Mol Genet. 2014;23(23):6407–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Yin H, Guo J, Ding E, Zhang H, Han L, Zhu B. Salt-Inducible kinase 3 haplotypes Associated with noise-Induced hearing loss in Chinese workers. Audiol Neurootol. 2020;25(4):200–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors are grateful to all the subjects and all the investigators in this study.

Funding

This work was supported by the Natural Foundation of Guangdong Province (2020A1515010582), Dongguan Science and Technology of Development key Program (202050715005209), Guangdong Basic and Applied Basic Research Foundation (2019A1515110656), Guangdong Medical Scientific Research Foundation (A2020141, C2022097) and the National Innovation Training Program for College Students (202210573019).

Author information

Authors and Affiliations

Authors

Contributions

Shan Wu, Zhidan Wu, and Manlian Chen, performed the laboratorial determination and the statistical analysis for this study, and wrote the manuscript. Xiangbin Zhong, Haoyan Gu, Wenjing Du and Weidong Liu, participated in the epidemiological investigation and laboratorial determination. Li Lang and Junyi Wang, designed concept for the study, revised the manuscriptand and obtained funding. All authors reviewed the manuscript and had final approval of the submitted version of the paper.

Corresponding authors

Correspondence to Li Lang or Junyi Wang.

Ethics declarations

Ethics approval and consent to participate

All the methods and procedures carried out in this study were in accordance with the Declaration of Helsinki. The study was approved by the Science Ethics Committee of Guangdong Pharmaceutical University (20190212). Informed consent was obtained form all the participants.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, S., Wu, Z., Chen, M. et al. Interactions of genetic variations in FAS, GJB2 and PTPRN2 are associated with noise-induced hearing loss: a case-control study in China. BMC Med Genomics 17, 18 (2024). https://doi.org/10.1186/s12920-023-01790-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-023-01790-7

Keywords