Skip to main content

Overinterpretation of high throughput sequencing data in medical genetics: first evidence against TMPRSS3/GJB2 digenic inheritance of hearing loss

Abstract

Background

Hearing loss (HL) is the most common disability of human senses characterized by a great allelic heterogeneity. GJB2 and TMPRSS3 are two well-known HL genes typically underlying its monogenic form. Recently, TMPRSS3/GJB2 digenic inheritance has been proposed. As results of genetic testing can be easily overinterpreted, we aimed to verify the hypothesis.

Methods

From genetic database of HL patients with at least one TMPRSS3 pathogenic variants we have selected individuals with additional GJB2 pathogenic variants. All of the available family members were recruited for the study. Segregation analysis of the respective TMPRSS3 and GJB2 pathogenic variants was performed within the families.

Results

The strategy has allowed to identify four individuals who were double heterozygous for known pathogenic TMPRSS3 and GJB2 variants. Two individuals from different families had GJB2 c.35delG and TMPRSS3 c.208delC and in two other individuals from one family GJB2 c.35delG together with TMPRSS3 c.1343T>C variants were found. None of these subjects has ever reported hearing problems and their hearing status was normal.

Conclusions

Our data provide evidence against TMPRSS3/GJB2 digenic inheritance of HL. As high throughput sequencing is increasingly used for genetic testing, particular caution should be taken to provide the patients with accurate genetic counseling.

Background

Digenic inheritance can be defined as a mechanism, which requires an interaction of two loci for expression of a phenotype. Each of the loci can exert a different level of influence on the phenotype, which means that (i) one locus may represent a primary locus or (ii) both loci may be roughly equal in importance. The first example of digenic inheritance of human diseases is retinitis pigmentosa attributed to recessive variants in the PRPH2 and ROM1 genes [1]. Their causative role was provided based on convincing data from different family studies and confirmed interaction between the two gene products.

The growing number of human disorders with a digenic pattern of inheritance is now being collected in digenic diseases database (DIDA), a comprehensive repository with records on genes and genetic variants involved in digenic diseases (http://dida.ibsquare.be/). Hereditary hearing loss (HL) is a genetically heterogenous condition with over 100 genes involved in its development, which makes HL a good candidate for digenic inheritance. Congenital hereditary hearing impairment is a common disorder which affects approx. 1:2000 neonates. The major genetic cause of HL are biallelic GJB2/GJB6 (DFNB1 locus) pathogenic variants [2]. Another important contributor to autosomal recessive non-syndromic HL is the TMPRSS3 gene [3].

In the study we have aimed to verify the hypothesis of putative TMPRSS3/GJB2 digenic inheritance of HL. The idea was precipitated by recent publications suggesting the novel gene combination involved in HL development and reinforced by inquires from health care professionals and patients on how to interpret the co-occurrence of heterozygous recessive variants in both genes. According to previous reports up to one-third of “solved cases” may be attributable to overinterpretation due to incorrect assessment of variant pathogenicity or poor-evidenced digenic or oligogenic inheritance [4,5,6].

Methods

Patients

The study was approved by the local ethics committee (IFPS:/KB/03/2012). From patients tested for TMPRSS3 pathogenic variants (n = 2277) we have first selected individuals with at least one TMPRSS3 pathogenic variant (n = 42) and next patients with additional pathogenic variants in the GJB2 gene (n = 4). The probands and their family members (n = 18) gave written informed consent for participation in the study. Detailed medical history was evaluated. Hearing status was determined by self-assessment or pure tone-audiometry with hearing thresholds measured at selected frequencies 125, 250, 500, 1000, 2000, 4000 and 6000 Hz using the AC40 clinical audiometer (Interacoustics, Middelfart, Denmark). Hearing thresholds up to 20 dB are considered to be normal hearing. DNA was isolated from blood samples or buccal swabs using a standard protocol. Presence of the respective TMPRSS3 (NM_024022.2, NP_076927.1) and GJB2 (NM_004004.5, NP_003995.2) pathogenic variants was verified by Sanger sequencing [3] and family segregation study was performed.

Databases

DIDA was used to obtain information on genetic variants involved in digenic HL (http://dida.ibsquare.be/; accessed 10/2018). Interaction networks were generated using STRING v.10.5 (https://string-db.org/; accessed 12/2018), a database of known and predicted protein–protein interactions. Expression pattern of GJB2 and TMPRSS3 in the inner ear was collected from previous studies (Table 1).

Table 1 Expression pattern of GJB2 and TMPRSS3 genes

Results

In Family 1 prelingual profound bilateral sensorineural HL was diagnosed in the proband and his brother. They received cochlear implants (CI) at the age of 4 y and 2 y, respectively. Genetic testing revealed that HL in the family was a consequence of homozygous pathogenic c.208delC (p.His70Thrfs*19) variant in the TMPRSS3 gene. Carrier status of the TMPRSS3 variant was confirmed in both parents. In addition to that both HL siblings and their normal hearing father (hearing status confirmed by pure tone audiometry) were carriers of a pathogenic heterozygous c.35delG (p.Gly12Valfs*2) variant in the GJB2 gene. Neither paternal grandparents nor his siblings were available for genetic testing but none of them have complained of hearing impairment (Fig. 1a).

Fig. 1
figure 1

Coexistence of GJB2 and TMPRSS3 pathogenic variants in normal hearing individuals supports lack of interactions between the two genes. ac Pedigrees of the analyzed families. Probands are marked with an arrow. Black symbols indicate individuals affected with HL and open symbols indicate unaffected individuals. Red rectangles mark normal hearing individuals being double heterozygous for GJB2 and TMPRSS3 pathogenic variants. d Map of GJB2 and TMPRSS3 interactions according to STRING v.10.5 database. Round symbols indicate different genes. Interactions between genes are represented by color lines described in the legend

The proband in Family 2 had progressive bilateral sensorineural HL affecting mainly high frequencies from the age of about 20 years. She received hearing aids at the age of 30 y and CI of the right ear was performed at the age of 40 y. The patient was compound heterozygous for two pathogenic variants, i.e. c.1276G>A (p.Ala426Thr) and c.1343T>C (p.Met448Thr) in the TMPRSS3 gene that were inherited from her normal hearing parents. The TMPRSS3 c.1343T>C variant was detected also in the proband’s father, brother and her 23-year old normal hearing son. Except for the TMPRSS3 variants the proband, her father and her son were also heterozygous for the pathogenic GJB2 c.35delG variant (Fig. 1b).

Family 3 was characterized by prelingual profound bilateral sensorineural HL diagnosed in the proband and both of his parents. The proband received CI at the end of the first year of life. He is compound heterozygous for GJB2 c.235delC (p.Leu79Cysfs*3) and c.313_326del (p.Lys105Glyfs*5) pathogenic variants. HL in the proband’s father resulted from two GJB2 pathogenic variants c.35delG and c.235delC in a trans configuration. Proband’s mother was compound heterozygous for GJB2 c.35delG and c.313_326del pathogenic variants. In addition to the GJB2 variants, the proband and his mother were also carriers of a TMPRSS3 c.208delC pathogenic variant. A combination of one TMPRSS3 (c.208delC) and one GJB2 (c.35delG) variant was detected in the proband’s maternal uncle, who does not have hearing impairment. Both, proband’s mother and her brother, have most probably inherited the combination of GJB2 and TMPRSS3 variants from their normal hearing father, who was not available for genetic testing (Fig. 1c).

In DIDA there are two records on possible TMPRSS3/GJB2 digenic inheritance of HL. The first one refers to a combination of GJB2 c.35delG together with TMPRSS3 c.208delC for which only familial evidence (based on one family) has been provided by Battelino et al. [7]. The second combination is GJB2 c.487A>G (p.Met163Val) and TMPRSS3 c.1276G>A (p.Ala426Thr) that is supposed to be evidenced by familial and functional data [8].

Detailed expression pattern of GJB2 and TMPRSS3 has been presented in Table 1. It overlaps partially and is restricted mainly to supporting cells. GJB2 encodes connexin 26, a gap junction protein, which forms cell-to-cell channels permeable in the inner ear to K+ ions and other small molecules [9]. TMPRSS3 is a transmembrane protein with protease activity. Recent studies indicate that it is involved in actin-related hair cell mechanics and support of cell motility and integrity [10]. Currently, there is no evidence of direct relationships between the two genes. In STRING neither known (from curated databases or experimentally determined) nor predicted direct interactions between TMPRSS3 and GJB2 genes and proteins was reported (Fig. 1d).

Discussion

Based on family studies we have identified as many as four individuals from three different HL families who were double heterozygous for pathogenic variants in the GJB2 and TMPRSS3 genes (Fig. 1a–c). Interestingly, none of the subjects was diagnosed with HL. Our data provide strong evidence against TMPRSS3/GJB2 digenic inheritance of HL that has been proposed by two recent reports [7, 8], which in our opinion are far from being convincing. Both studies rely exclusively on data from only two families with either one [7] or two [8] affected individuals carrying heterozygous GJB2 and TMPRSS3 variants.

In the study by Battelino et al. one patient has been proposed to have HL due to digenic inheritance of GJB2 c.35delG and TMPRSS3 c.208delC. However, a second variant c.579dupA in TMPRSS3 was also detected but it was erroneously interpreted as occurring in the non-coding sequence and being non-pathogenic. As noticed by our group [11] congenital profound HL in the patient is not a consequence of TMPRSS3/GJB2 digenic inheritance but results from compound heterozygous c.208delC (p.His70Thrfs*19) and c.579dupA (p.Cys194Metfs*17) TMPRSS3 pathogenic variants.

The second paper on the presumed TMPRSS3/GJB2 digenic inheritance of HL presents two siblings with moderate-to-severe HL at mid and high frequencies. Sequencing of 71 known HL genes revealed a heterozygous GJB2 c.487A>G (p.Met163Val) variant and a heterozygous TMPRSS3 c.1276G>A (p.Ala426Thr) variant. Hearing status of the patients' parents was normal and each of them carried one of the detected variants [8]. While the TMPRSS3 variant was found as causative of postlingual progressive HL with recessive mode of inheritance [3] and its pathogenic potential was confirmed by in vitro studies [12], the pathogenic role of the GJB2 p.Met163Val variant is still intriguing. Up to now it has been detected in approximately 18 HL patients. In the majority of them p.Met163Val was identified in a simple heterozygous state and its pathogenicity was described as unknown. In vitro studies have shown that p.Met163Val affects the formation of homotypic connexin 26 junctional channels [13]. According to ACMG/AMP standards and guidelines for the interpretation of sequence variants, p.Met163Val should be classified as likely pathogenic [14].

There are at least three different possible explanations for the genetic cause of HL (other than TMPRSS3/GJB2 digenic inheritance) in the family studied by Leone et al. [8]. The first option is a dominant character of the identified GJB2 p.Met163Val variant with a possible incomplete penetrance in the patients’ mother. Another alternative is the omission of the second pathogenic DFNB1 variant due to e.g. poor quality sequencing, incomplete coverage or presence of undetected copy number variants. The third option is the involvement of genes other than GJB2 or TMPRSS3. Testing of a limited number of HL genes does not exclude the possibility that the genetic cause of HL is located in the untargeted regions of these genes such as alternative exons, introns or regulatory regions as well as in other known (currently more than 100 HL genes have been discovered) or still unknown HL genes or regulatory RNAs.

Among the relationships determining a presumed digenic inheritance between a pair of genes are their direct or indirect interactions, involvement in a common pathway, co-expression or similar function [15]. This should be taken into account if a digenic combination segregates with phenotype in a family. Based on what is currently known about the expression pattern, biological function or interaction between GJB2 and TMPRSS3, not only in the inner ear but also in other organs, it is difficult to find a common denominator for GJB2 and TMPRSS3 that could justify a phenomenon of their digenic inheritance. To experimentally test the hypothesis of TMPRSS3/GJB2 digenic inheritance, introduction of animal models and performing extensive cell biology experiments would be relevant [16,17,18,19]. In our opinion detection of double heterozygous TMPRSS3/GJB2 variants in HL patients is a purely coincidental finding. The carrier frequency of GJB2 variants in the general population has been estimated at approximately 3% [2] while TMPRSS3 at 0.4% (Lechowicz U. unpublished data). Considering the frequencies, one should expect a higher prevalence of HL due to digenic TMPRSS3/GJB2 inheritance than due to TMPRSS3 pathogenic variants alone. Since the publication by Leone et al. [8] there are no other reports confirming the presumed new phenomenon of TMPRSS3/GJB2 digenic inheritance. GJB2 is the most common cause of hereditary HL and one should keep in mind that it can be frequently find in tested individuals. Moreover, there is a significant enrichment of simple heterozygous GJB2 pathogenic variants in HL patients (~ 5% vs 2-3% in the general population) [20]. This observation may be a consequence of a higher complexity of the DFNB1 locus than currently accepted or related to an increased number of pathogenic variants in HL genes found in HL patients, referred to as mutational load [21].

Conclusion

High throughput sequencing data in HL patients may suggest a more complex oligogenic inheritance of HL but our results provide evidence that simple co-occurrence of heterozygous GJB2 and TMPRSS3 recessive variants is not related to the development of disabling HL. The number of detected, probably pathogenic variants increases with the use of advanced sequencing techniques, which may lead to overinterpretation of genotyping results. Particular caution and use of specific guidelines is required by data analysis and selection of causative variants in order to provide accurate diagnosis and counseling for the patients and reliable data for the medical genetic community [5, 14].

Availability of data and materials

Please contact author for data requests.

Abbreviations

HL:

hearing loss

ACMG/AMP:

the American College of Medical Genetics and Genomics and the Association for Molecular Pathology

DIDA:

digenic diseases database

STRING:

Search Tool for the Retrieval of Interacting Genes/Proteins

References

  1. Deltas C. Digenic inheritance and genetic modifiers. Clin Genet. 2018;93:429–38.

    Article  CAS  Google Scholar 

  2. Snoeckx RL, Huygen PL, Feldmann D, Marlin S, Denoyelle F, Waligora J, Mueller-Malesinska M, Pollak A, Ploski R, Murgia A, et al. GJB2 mutations and degree of hearing loss: a multicenter study. Am J Hum Genet. 2005;77:945–57.

    Article  CAS  Google Scholar 

  3. Lechowicz U, Gambin T, Pollak A, Podgorska A, Stawinski P, Franke A, Petersen BS, Firczuk M, Oldak M, Skarzynski H, Ploski R. Iterative Sequencing and Variant Screening (ISVS) as a novel pathogenic mutations search strategy—application for TMPRSS3 mutations screen. Sci Rep. 2017;7:2543.

    Article  Google Scholar 

  4. Farwell KD, Shahmirzadi L, El-Khechen D, Powis Z, Chao EC, Tippin Davis B, Baxter RM, Zeng W, Mroske C, Parra MC, et al. Enhanced utility of family-centered diagnostic exome sequencing with inheritance model-based analysis: results from 500 unselected families with undiagnosed genetic conditions. Genet Med. 2015;17:578–86.

    Article  CAS  Google Scholar 

  5. DiStefano MT, Hemphill SE, Oza AM, Siegert RK, Grant AR, Hughes MY, Cushman BJ, Azaiez H, Booth KT, Chapin A, et al. ClinGen expert clinical validity curation of 164 hearing loss gene-disease pairs. Genet Med. 2019. https://doi.org/10.1038/s41436-019-0487-0.

    Article  PubMed  Google Scholar 

  6. Phelps IG, Dempsey JC, Grout ME, Isabella CR, Tully HM, Doherty D, Bachmann-Gagescu R. Interpreting the clinical significance of combined variants in multiple recessive disease genes: systematic investigation of Joubert syndrome yields little support for oligogenicity. Genet Med. 2018;20:223–33.

    Article  Google Scholar 

  7. Battelino S, Klancar G, Kovac J, Battelino T, Trebusak Podkrajsek K. TMPRSS3 mutations in autosomal recessive nonsyndromic hearing loss. Eur Arch Otorhinolaryngol. 2016;273:1151–4.

    Article  Google Scholar 

  8. Leone MP, Palumbo P, Ortore R, Castellana S, Palumbo O, Melchionda S, Palladino T, Stallone R, Mazza T, Cocchi R, Carella M. Putative TMPRSS3/GJB2 digenic inheritance of hearing loss detected by targeted resequencing. Mol Cell Probes. 2017;33:24–7.

    Article  CAS  Google Scholar 

  9. Liu W, Bostrom M, Kinnefors A, Rask-Andersen H. Unique expression of connexins in the human cochlea. Hear Res. 2009;250:55–62.

    Article  CAS  Google Scholar 

  10. Liu W, Lowenheim H, Santi PA, Glueckert R, Schrott-Fischer A, Rask-Andersen H. Expression of trans-membrane serine protease 3 (TMPRSS3) in the human organ of Corti. Cell Tissue Res. 2018;372:445–56.

    Article  CAS  Google Scholar 

  11. Lechowicz U, Pollak A, Ozieblo D, Oldak M. Pathogenic p.Cys194Metfs*17 variant argues against TMPRSS3/GJB2 digenic inheritance of hearing loss. Eur Arch Otorhinolaryngol. 2016;273:1327–8.

    Article  Google Scholar 

  12. Lee YJ, Park D, Kim SY, Park WJ. Pathogenic mutations but not polymorphisms in congenital and childhood onset autosomal recessive deafness disrupt the proteolytic activity of TMPRSS3. J Med Genet. 2003;40:629–31.

    Article  CAS  Google Scholar 

  13. Bruzzone R, Veronesi V, Gomes D, Bicego M, Duval N, Marlin S, Petit C, D’Andrea P, White TW. Loss-of-function and residual channel activity of connexin26 mutations associated with non-syndromic deafness. FEBS Lett. 2003;533:79–88.

    Article  CAS  Google Scholar 

  14. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24.

    Article  Google Scholar 

  15. Gazzo AM, Daneels D, Cilia E, Bonduelle M, Abramowicz M, Van Dooren S, Smits G, Lenaerts T. DIDA: a curated and annotated digenic diseases database. Nucleic Acids Res. 2016;44:D900–7.

    Article  CAS  Google Scholar 

  16. Zheng QY, Scarborough JD, Zheng Y, Yu H, Choi D, Gillespie PG. Digenic inheritance of deafness caused by 8 J allele of myosin-VIIA and mutations in other Usher I genes. Hum Mol Genet. 2012;21:2588–98.

    Article  CAS  Google Scholar 

  17. Yan D, Liu XZ. Modifiers of hearing impairment in humans and mice. Curr Genomics. 2010;11:269–78.

    Article  CAS  Google Scholar 

  18. Seo S, Baye LM, Schulz NP, Beck JS, Zhang Q, Slusarski DC, Sheffield VC. BBS6, BBS10, and BBS12 form a complex with CCT/TRiC family chaperonins and mediate BBSome assembly. Proc Natl Acad Sci USA. 2010;107:1488–93.

    Article  CAS  Google Scholar 

  19. Badano JL, Leitch CC, Ansley SJ, May-Simera H, Lawson S, Lewis RA, Beales PL, Dietz HC, Fisher S, Katsanis N. Dissection of epistasis in oligogenic Bardet–Biedl syndrome. Nature. 2006;439:326–30.

    Article  CAS  Google Scholar 

  20. Pollak A, Mueller-Malesinska M, Skorka A, Kostrzewa G, Oldak M, Korniszewski L, Skarzynski H, Ploski R. GJB2 and hearing impairment: promoter defects do not explain the excess of monoallelic mutations. J Med Genet. 2008;45:607–8.

    Article  CAS  Google Scholar 

  21. Vona B, Muller T, Nanda I, Neuner C, Hofrichter MA, Schroder J, Bartsch O, Lassig A, Keilmann A, Schraven S, et al. Targeted next-generation sequencing of deafness genes in hearing-impaired individuals uncovers informative mutations. Genet Med. 2014;16:945–53.

    Article  CAS  Google Scholar 

  22. Lautermann J, ten Cate WJ, Altenhoff P, Grummer R, Traub O, Frank H, Jahnke K, Winterhager E. Expression of the gap-junction connexins 26 and 30 in the rat cochlea. Cell Tissue Res. 1998;294:415–20.

    Article  CAS  Google Scholar 

  23. Guipponi M, Antonarakis SE, Scott HS. TMPRSS3, a type II transmembrane serine protease mutated in non-syndromic autosomal recessive deafness. Front Biosci. 2008;13:1557–67.

    Article  CAS  Google Scholar 

  24. Fasquelle L, Scott HS, Lenoir M, Wang J, Rebillard G, Gaboyard S, Venteo S, Francois F, Mausset-Bonnefont AL, Antonarakis SE, et al. TMPRSS3, a transmembrane serine protease deficient in human DFNB8/10 deafness, is critical for cochlear hair cell survival at the onset of hearing. J Biol Chem. 2011;286:17383–97.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by Polish National Science Centre Grant No. 2012/05/N/NZ5/02629.

Author information

Authors and Affiliations

Authors

Contributions

MO analyzed the data and wrote the manuscript; AP, UL performed genotyping and analyzed the data; MO, UL, AP, DO, HS participated in phenotyping and clinical data collection; MO, DO performed computational analysis and prepared figure and table. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Monika Ołdak.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the bioethics committee at the Institute of Physiology and Pathology of Hearing and performed according to the Declaration of Helsinki. Written informed consent was obtained from each participant.

Consent for publication

A consent for publication has been obtained.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ołdak, M., Lechowicz, U., Pollak, A. et al. Overinterpretation of high throughput sequencing data in medical genetics: first evidence against TMPRSS3/GJB2 digenic inheritance of hearing loss. J Transl Med 17, 269 (2019). https://doi.org/10.1186/s12967-019-2018-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-2018-9

Keywords