BioMed Research International

BioMed Research International / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 6624744 | https://doi.org/10.1155/2021/6624744

Jinying Li, Hongen Xu, Jianfeng Sun, Yongan Tian, Danhua Liu, Yaping Qin, Huanfei Liu, Ruijun Li, Lingling Neng, Xiaohua Deng, Binbin Xue, Changyun Yu, Wenxue Tang, "Missense Variant of Endoplasmic Reticulum Region of WFS1 Gene Causes Autosomal Dominant Hearing Loss without Syndromic Phenotype", BioMed Research International, vol. 2021, Article ID 6624744, 9 pages, 2021. https://doi.org/10.1155/2021/6624744

Missense Variant of Endoplasmic Reticulum Region of WFS1 Gene Causes Autosomal Dominant Hearing Loss without Syndromic Phenotype

Academic Editor: Burak Durmaz
Received24 Nov 2020
Revised01 Feb 2021
Accepted14 Feb 2021
Published04 Mar 2021

Abstract

Objective. Genetic variants in the WFS1 gene can cause Wolfram syndrome (WS) or autosomal dominant nonsyndromic low-frequency hearing loss (HL). This study is aimed at investigating the molecular basis of HL in an affected Chinese family and the genotype-phenotype correlation of WFS1 variants. Methods. The clinical phenotype of the five-generation Chinese family was characterized using audiological examinations and pedigree analysis. Target exome sequencing of 129 known deafness genes and bioinformatics analysis were performed among six patients and four normal subjects to screen suspected pathogenic variants. We built a complete WFS1 protein model to assess the potential effects of the variant on protein structure. Results. A novel heterozygous pathogenic variant NM_006005.3 c.2020G>T (p.Gly674Trp) was identified in the WFS1 gene, located in the C-terminal domain of the wolframin protein. We further showed that HL-related WFS1 missense variants were mainly concentrated in the endoplasmic reticulum (ER) domain. In contrast, WS-related missense variants are randomly distributed throughout the protein. Conclusions. In this family, we identified a novel variant p.Gly674Trp of WFS1 as the primary pathogenic variant causing the low-frequency sensorineural HL, enriching the mutational spectrum of the WFS1 gene.

1. Introduction

Hearing loss (HL) is a very common sensory disorder and can be caused by genetic factors, viral infections, drugs, etc. [1]. More than 50% of prelingual HL cases are thought to be related to genetic factors [2]. The World Health Organization estimates that approximately 466 million people worldwide suffer from HL, including 34 million children (https://www.who.int/zh/news-room/fact-sheets/detail/deafness-and-hearing-loss). There are four genetic patterns of deafness: autosomal dominant (15–20%), autosomal recessive (80%), X-linked (1%), and mitochondrial DNA inheritance (1%). According to the presence of abnormalities in other organs, deafness is divided into syndromic (~30%) and nonsyndromic (~70%) HL. Deafness has extensive genetic heterogeneity [3], and 119 genes causing nonsyndromic HL have been identified (https://hereditaryhearingloss.org/). Different variants in the same gene, such as CDH23 [46], SLC26A4 [7, 8], and WFS1 [9, 10], may lead to either syndromic or nonsyndromic HL.

The WFS1 gene encodes a predicted 890-amino-acid transmembrane protein with a calculated molecular mass of approximately 100 kDa, which maps to chromosome 4q16. WFS1 is predicted to have nine central transmembrane domains, with an extracytoplasmic N terminus and an intracytoplasmic C terminus. In 1998, it was identified as the gene causing Wolfram syndrome, an autosomal recessive disorder, also called DIDMOAD (diabetes insipidus, diabetes mellitus, optic atrophy, and deafness) [11]. The deafness associated with DIDMOAD is characterized by high-frequency sensorineural HL. Some variants in this gene can also cause autosomal dominant deafness 6 (DFNA6), also known as DFNA14 or DFNA38, characterized by low-frequency sensorineural hearing impairment. WFS1 plays an important role in maintaining endoplasmic reticulum (ER) homeostasis, and pathogenic variants in WFS1 can lead to ER stress and cause early cell dysfunction and death [12, 13].

Next-generation sequencing technologies have been developed to detect the pathogenic variants underlying Mendelian disorders [14, 15]. We applied target exome sequencing of know deafness genes to determine the variants causing nonsyndromic postlingual deafness in a large Chinese family. A novel missense variant NM_006005.3 c.2020G>T (p.Gly674Trp) was identified in the WFS1 gene that cosegregated with the HL phenotype. By querying the Human Gene Mutation Database (HGMD), ClinVar, and Deafness Variation Database (http://deafnessvariationdatabase.org/) and reviewing the literature, we found that WFS1 missense variants correlated with HL and WS are different in the location relative to the protein domain.

2. Materials and Methods

2.1. Ethics Statement

The study was approved by the Committee of Medical Ethics of Zhengzhou University. Informed consent was obtained from all participants.

2.2. Clinical Evaluations

This study was conducted in a five-generation Chinese family from a village in Henan Province, China. This family has 65 members, including 11 with postlingual nonsyndromic HL. Among them, two patients had died and three patients refused to participate in this study; the remaining patients (II-2, III-2, III-8, III-10, III-12, and IV-3) and four normal subjects (III-1, III-3, III-11, and IV-4) underwent pure-tone audiometry (PTA), genetic test and clinical examinations, including internal auditory canal magnetic resonance imaging (MRI), computed tomography (CT) of the temporal bone, vision test, and fasting blood glucose level test.

2.3. Target Enrichment Sequencing

Genomic DNA was isolated from patient samples using the GenMagBio Genomic DNA Purification kit (GenMagBio, Changzhou, China) as per the manufacturer’s standard procedures. Genomic DNA from affected and unaffected members was fragmented to an average size of 250 bp, and end repair, adapter ligation, and PCR enrichment were performed following the protocol of the VAHTS™ Universal DNA Library Prep Kit for Illumina V3 (Vazyme Biotech, Nanjing, China). Sequence capture was performed using the Human Deafness Panel oto-DA3 (Otogenetics, Atlanta, GA, USA), containing 129 known HL genes, following the manufacturer’s protocol. The resulting libraries were sequenced on an Illumina HiSeq 4000 sequencer (Illumina, San Diego, CA, USA) with 150 bp paired-ends.

2.4. Bioinformatics Analysis and Variant Interpretation

Bioinformatics analysis and variant interpretation were carried out as described previously [16]. Sequencing adapters and low-quality reads were trimmed from the raw reads using Trimmomatic [17]. Clean reads were aligned to the human reference genome (ver. GRCh37) using the Burrow–Wheeler Aligner (ver. 0.7.17-r1188) [18], followed by duplicate reads marking using sambamba (ver. 0.6.6) [19]. Variant (single nucleotide variants and small indels) and genotype calling were performed using the Genome Analysis Toolkit ver. 4 (GATK4) HaplotypeCaller [20]. Variants were annotated by vcfanno [21] using the 1000 Genomes Project database [22], dbSNP [23], Exome Aggregation Consortium (ExAC) [24], Genome Aggregation Database (gnomAD) [25], ClinVar [24], InterVar [26], and dbNSFP [27]; the latter database compiled variant prediction scores from many prediction algorithms [2830]. All analysis steps described above were performed in the framework of bcbio-nextgen (https://github.com/bcbio/bcbio-nextgen). We filtered out variants with minor allele in any general continental population, in which at least 2,000 alleles were observed in the gnomAD database, except those in the American College of Medical Genetics and Genomics (ACMG) benign stand-alone exception list or linked to diseases in the ClinVar database [31]. Variants were interpreted by an expert panel consisting of an otorhinolaryngologist, bioinformatician, and molecular geneticist following the guidelines of the ACMG and Association for Molecular Pathology for clinical sequence interpretation [32]. The variant nomenclature was based on the WFS1 canonical transcript NM_006005.3.

2.5. Sanger Sequencing

To confirm candidate variants detected by next-generation sequencing, we performed PCR amplification and Sanger sequencing of the DNA from 10 individuals. Forward (5-CCGGTGGTTCACGTCTCTGG-3) and reverse (5-GCCAGCAGCTTAAGGCGAC-3) primers were designed using NCBI Primer-BLAST and synthesized by Shangya Biotechnology (Zhengzhou, China). PCR was conducted using the 2x Taq Master Mix kit (Novoprotein, Shanghai, China). The amplified products were subjected to 2.2% agarose gel electrophoresis, purified using a PCR purification kit (Lifefeng Biotechnology, Shanghai, China), and then sequenced using the SeqStudio Genetic Analyzer (Applied Biosystems/Life Technologies, Carlsbad, CA, USA).

2.6. Protein Structure Prediction of WFS1

The three-dimension (3D) model of the wolframin protein was built by a contact-assisted protein structure prediction method of CONFOLD2 [33], which leverages the information of interresidue contacts and secondary structures as well as energy functions. Based on the wolframin multiple sequence alignment (MSA) generated by HHblits [34], the interresidue contacts were then deduced from those residue pairs coupled evolutionarily by using CCMpred [35]. The secondary structures were predicted by SCRATCH1D [36]. We set for each threshold 5 best models and used the top-1 model selected by the final decision of CONFOLD2.

The reported missense variants at position 674 of the wolframin were used to perform in silico mutagenesis. The residue (Gly) was substituted to other four amino acids (Glu, Val, Arg, and Trp) as reported in the HGMD database, and an ensemble of the conformations (the number of conformations limits to 25) was generated for each mutant by low-mode MD (Molecular Dynamics). The force field used for calculation was OPLS-AA, and the implicit solvent was the reaction field (R-Field) model. All calculations were performed in Molecular Operating Environment (MOE) 2018 package [37]. Ab initio structure prediction of WFS1 was used to perform in silico mutagenesis.

3. Results

3.1. Clinical Features of the Family

The five-generation family has 11 affected individuals, ranging in age from 30 to 83 years (Figure 1). The family had a typical autosomal dominant inheritance pattern. In this family, the age of onset ranged from 20 to 30 years. The proband (III-2) developed moderate sensorineural HL involving low frequencies; the average auditory thresholds on PTA for the right and left ears were 67 and 63 dB HL, respectively (Figure 2(a)). In this family, subject II-2 is an 83-year-old lady, and her high-frequency hearing impairment might be related to presbycusis. The other affected members showed postlingual, bilateral nonsyndromic sensorineural HL involving low frequencies that was mild to profound. The patients had flat-sloping audiograms (Figure 2(a)). MRI and temporal bone CT in the proband showed a normal vestibular aqueduct and internal auditory canal (Figure 2(b)). No patient reported vestibular dysfunction, such as balance disorders, vertigo, or Meniere disease. All family members had normal blood glucose levels and visual function. There was no history of exposure to aminoglycosides or noise. Comprehensive examinations showed no evidence of DIDMOAD syndrome. Table 1 shows the main clinical features of the affected individuals.


Subject no.GenderAgeAge of onsetUsage time for HAAudiogram shapeTinnitusVertigo

II-2Female832220Flat-slopingYesNo
III-2Female60208Flat-slopingYesNo
III-8Female55304Flat-slopingNoNo
III-10Female57257.5Flat-slopingNoNo
III-12Female632011Flat-slopingYesNo
IV-3Male30198Flat-slopingYesNo

HA: hearing aid.
3.2. Target Enrichment Sequencing and Sanger Sequencing

Target exome sequencing of 129 known deafness genes was performed in six affected and four unaffected individuals (Figure 1) to identify potentially pathogenic variants. An average of one billion raw base pairs (1 Gbp) was generated for each sample, with more than 87% of the bases having a Phred quality . The mapping rate to the GRCh37 human genome reference sequence exceeded 99.7%, and the average sequencing depth was 164.98x (range from 143.8x to 184.9x), covering 95.8% of the target regions at least 20x. The only variant in WFS1, NM_006005.3, c.2020G>T, cosegregated with HL in the family. Sanger sequencing of this variant in the family verified that the variant was present in the affected members and absent in the unaffected ones (Figure 3).

3.3. Variant Interpretation of the Variants as Pathogenic

The WFS1 (NM_006005.3.) variant c.2020G>T (p.Gly674Trp) is present in exon 8, located in the C-terminal domain of the encoded wolframin protein. According to the standards and guidelines for interpreting genetic variants and the expert specification for genetic hearing loss proposed by the ACMG and the Association for Molecular Pathology [32, 38], this variant was classified as pathogenic (Table 2). This variant is absent from all population databases, including gnomAD and ExAC (Pathogenic Moderate 2, PM2). The variant cosegregated with the DFNA6/14/38 phenotype in the affected individuals (six affected, Pathogenic Supportive 1 Strong, PP1_Strong). The low-frequency HL in this family was highly specific for WFS1 (Pathogenic Supportive 4, PP4). The REVEL score of this variant was ≥0.75 (Pathogenic Supportive 3, PP3). A known likely pathogenic variant (c.2021G>T, p.Gly674Val) is at the same amino acid residue with c.2020G>T (PM5). The evidence of p.Gly674Val is as follows: this variant was also absent from population databases (ACMG PM2), two probands with the variant (ACMG PS4_supporting) [10, 39], the variant segregated with the hearing impairment in a Dutch family (ACMG PP1) [10], REVEL (ACMG PP3), and patient’s phenotype or family history highly specific for a disease with a single genetic etiology (ACMG PP4).


Gene symbolNucleotide changeAmino acid changeSIFTPolyPhen-2MutationTasterPathogenicityACMG evidence

WFS1c.2020G>Tp.Gly674TrpDamagingDamagingDamagingPathogenicPM2; PM5; PS4_S; PP1_S; PP3; PP4

PP1_S: PP1_Strong; PS4_S: PS4_Supporting.
3.4. Missense Variants Leading to HL and WS in Wolframin Protein

Missense variants in wolframin protein can cause both WS and nonsyndromic HL. We classified these variants of WFS1 according to different phenotypes and drew a distribution map of these variants in wolframin protein. We found that the missense variants that caused WS were distributed in all protein regions, while the missense variants that caused HL were clustered in the ER lumen domain (Figure 4).

To reveal the potential mechanism of different missense variants at the same protein position lead to varying diseases, we assessed the potential effects of variants on protein structure by taking four missense variants (nonsyndromic HL: p.Gly674Trp, p.Gly674Glu, and p.Gly674Val; WS: p.Gly674Arg) as an example. We built a complete WFS1 protein model from scratch as there are no three-dimensional molecular structures for any parts of WFS1 available. The three-dimensional structure revealed that all of the four missense variants cause loss of the hydrogen bond interaction between p.Gly674 and p.Thr663, which was observed in wild type. Notably, compared with WS-associated p.Gly674Arg, all the three HL-related variants increased an extra hydrogen bond between p.668Gln and p.675Pro (Figure 5). Changes in hydrogen bonds may affect the stability of the WFS1 protein, suggesting a potential molecular mechanism resulting in varying diseases.

4. Discussions

We examined a five-generation Han Chinese family with postlingual nonsyndromic autosomal dominant HL using target sequencing of 129 genes associated with HL. Extensive clinical evaluation of the family before molecular genetic analysis confirmed nonsyndromic HL with no associated features that segregated with the HL. A novel WFS1 c.2020G>T (p.Gly674Trp) missense variants was the only variant cosegregating with HL identified in this family. We further showed that compared with missense variants associated with WS, HL-related ones are mainly located in the ER domain.

WFS1 encodes a transmembrane protein (wolframin) that is located primarily in the ER; it was highly expressed in the brain, pancreas, heart, and insulinoma beta-cell lines [40] and differentially expressed in inner ear cells [41]. Wolframin plays a role in protein folding [42] and negatively regulates ER stress [13], maintaining endolymphatic ion homeostasis. Once a variant occurs in WFS1, the negative regulation of a feedback loop in the ER stress signaling network disappears, resulting in the accumulation of misfolded and unfolded proteins, leading to cell death [43, 44]. In the inner ear, mutated WFS1 may alter intracellular Ca2+ homeostasis [42], the recycling of K+ back to the endolymph [45], and other endolymphatic ions. p.Gly674Trp is located in the ER lumenal domain. Glycine is most often replaced by charged amino acids, but the location of the variant rather than the properties of the substituting amino acid has a greater impact on disease severity [46]. We speculate that the variants strongly induce ER stress in inner ear cells and disrupt endolymphatic ion composition and homeostasis, which is the leading cause of deafness.

Different variants in WFS1 have different effects on wolframin. Hofmann et al. first show that missense variants do not affect the stability of WFS1 mRNA, but the wolframinR629W protein level was reduced substantially compared with wild-type wolframin [47]. Different amino acid substitutions at the same site can cause different phenotypes. For example, different substitutions at amino acid 678 of WFS1 resulted in three distinct phenotypes: sensorineural HL [48, 49], Wolfram syndrome [50], and optic atrophy and diabetes [51]. The missense variants c.2032T>C (p.Trp678Arg) and c.2033G>T (p.Trp678Leu) lead only to sensorineural HL, whereas the nonsense variants c.2033G>A (p.Trp678Term) and c.2034G>A (p.Trp678Term) lead to more severe symptoms of Wolfram syndrome or optic atrophy and diabetes. The novel missense variant c.2020G>T (p.Gly674Trp) identified in this study caused only nonsyndromic HL in the family. Interestingly, a change from glycine to arginine at the same amino acid position (c.2020G>A, p.Gly674Arg) led to a severe Wolfram syndrome phenotype [52]. These findings can contribute to a better understanding of different variants that lead to phenotypes of different severities.

In conclusion, we identified a novel variant c.2020G>T (p.Gly674Trp) in WFS1 responsible for autosomal dominant low-frequency sensorineural HL in a Chinese family, extending the variant spectrum of WFS1. Based on the variant location in the ER domain, we speculate that the variant might cause hearing impairment due to ER stress. Due to the instability of inner ear cells in vitro, the mechanism of low-frequency sensorineural HL caused by WFS1 variants is unclear. Further studies should explore which domains of the wolframin protein function in hearing and determine the pathogenic effects of genetic variants on the protein.

Abbreviations

HL:Hearing loss
ER:Endoplasmic reticulum
PTA:Pure-tone audiometry
MRI:Magnetic resonance imaging
CT:Computed tomography
DIDMOAD:Diabetes insipidus, diabetes mellitus, optic atrophy, and deafness
DFNA6/14/38:Autosomal dominant deafness 6/14/38
ACMG:American College of Medical Genetics and Genomics
HGMD:Human Gene Mutation Database
MOE:Molecular Operating Environment
MD:Molecular Dynamics.

Data Availability

All data can be viewed in NODE (http://www.biosino.org/node) by pasting the accession (OEP001069) into the text search box or through the URL: https://www.biosino.org/node/project/detail/OEP001069.

Conflicts of Interest

We declare that we have no conflict of interest in connection with the work submitted.

Authors’ Contributions

Changyun Yu performed genetic counseling. Wenxue Tang designed the study. Jinying Li, Hongen Xu, Changyun Yu, and Wenxue Tang performed the original draft preparation, data analysis, and revision of the manuscript. Huanfei Liu, Ruijun Li, Yongan Tian, and Danhua Liu performed the whole-exome sequencing, Sanger sequencing, and variant interpretation. Jianfeng Sun and Yaping Qin built the protein model. Lingling Neng, Xiaohua Deng, and Binbin Xue are responsible for information collection and phenotypic analysis. All authors have read and approved the final manuscript. Jinying Li and Hongen Xu contributed equally to this work. Jinying Li and Hongen Xu are co-first authors.

Acknowledgments

We thank all of the family members for participating in this study. We thank the Supercomputing Center of Zhengzhou University for providing computational and storage resources. The study was funded by the Chinese National Natural Science Foundation (grant no. 81402232), Key Research and Development and Promotion Project in Henan Province (grant no. 202102310116), Key Scientific Research of Colleges and Universities of Henan Province (grant no. 20A320022), and the Collaborative Innovation Project of Zhengzhou (Zhengzhou University) (Grant no. 18XTZX112004).

References

  1. T. B. Friedman and A. J. Griffith, “Human nonsyndromic sensorineural deafness,” Annual Review of Genomics and Human Genetics, vol. 4, no. 1, pp. 341–402, 2003. View at: Publisher Site | Google Scholar
  2. G. Van Camp, P. J. Willems, and R. J. Smith, “Nonsyndromic hearing impairment: unparalleled heterogeneity,” American Journal of Human Genetics, vol. 60, no. 4, pp. 758–764, 1997. View at: Google Scholar
  3. Z. Brownstein, Y. Bhonker, and K. B. Avraham, “High-throughput sequencing to decipher the genetic heterogeneity of deafness,” Genome Biology, vol. 13, no. 5, p. 245, 2012. View at: Publisher Site | Google Scholar
  4. K. Bommakanti, J. S. Iyer, and K. M. Stankovic, “Cochlear histopathology in human genetic hearing loss: state of the science and future prospects,” Hearing Research, vol. 382, p. 107785, 2019. View at: Publisher Site | Google Scholar
  5. J. M. Bork, L. M. Peters, S. Riazuddin et al., “Usher syndrome 1D and nonsyndromic autosomal recessive deafness DFNB12 are caused by allelic mutations of the novel cadherin-like gene CDH23,” American Journal of Human Genetics, vol. 68, no. 1, pp. 26–37, 2001. View at: Publisher Site | Google Scholar
  6. R. Valero, M. de Castro-Miró, S. Jiménez-Ochoa, J. J. Rodríguez-Ezcurra, G. Marfany, and R. Gonzàlez-Duarte, “Aberrant splicing events associated to CDH23 noncanonical splice site mutations in a proband with atypical usher syndrome 1,” Genes, vol. 10, 2019. View at: Publisher Site | Google Scholar
  7. L. A. Everett, H. Morsli, D. K. Wu, and E. D. Green, “Expression pattern of the mouse ortholog of the Pendred’s syndrome gene (Pds) suggests a key role for pendrin in the inner ear,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 17, pp. 9727–9732, 1999. View at: Publisher Site | Google Scholar
  8. Y. Yuan, Y. You, D. Huang et al., “Comprehensive molecular etiology analysis of nonsyndromic hearing impairment from typical areas in China,” Journal of Translational Medicine, vol. 7, no. 1, p. 79, 2009. View at: Publisher Site | Google Scholar
  9. G. Yu, M. L. Yu, J. F. Wang, C. R. Gao, and Z. J. Chen, “WS1 gene mutation analysis of Wolfram syndrome in a Chinese patient and a systematic review of literatures,” Endocrine, vol. 38, no. 2, pp. 147–152, 2010. View at: Publisher Site | Google Scholar
  10. K. Cryns, M. Pfister, R. J. Pennings et al., “Mutations in the WFS1 gene that cause low-frequency sensorineural hearing loss are small non-inactivating mutations,” Human Genetics, vol. 110, no. 5, pp. 389–394, 2002. View at: Publisher Site | Google Scholar
  11. T. M. Strom, K. Hörtnagel, S. Hofmann et al., “Diabetes insipidus, diabetes mellitus, optic atrophy and deafness (DIDMOAD) caused by mutations in a novel gene (wolframin) coding for a predicted transmembrane protein,” Human Molecular Genetics, vol. 7, no. 13, pp. 2021–2028, 1998. View at: Publisher Site | Google Scholar
  12. T. Hara, J. Mahadevan, K. Kanekura, M. Hara, S. Lu, and F. Urano, “Calcium efflux from the endoplasmic reticulum leads to β-cell death,” Endocrinology, vol. 155, no. 3, pp. 758–768, 2014. View at: Publisher Site | Google Scholar
  13. S. G. Fonseca, S. Ishigaki, C. M. Oslowski et al., “Wolfram syndrome 1 gene negatively regulates ER stress signaling in rodent and human cells,” The Journal of Clinical Investigation, vol. 120, no. 3, pp. 744–755, 2010. View at: Publisher Site | Google Scholar
  14. H. L. Rehm, “Disease-targeted sequencing: a cornerstone in the clinic,” Nature Reviews. Genetics, vol. 14, no. 4, pp. 295–300, 2013. View at: Publisher Site | Google Scholar
  15. S. Zou, X. Mei, W. Yang, R. Zhu, T. Yang, and H. Hu, “Whole-exome sequencing identifies rare pathogenic and candidate variants in sporadic Chinese Han deaf patients,” Clinical Genetics, vol. 97, no. 2, pp. 352–356, 2019. View at: Publisher Site | Google Scholar
  16. Y. Mi, D. Liu, B. Zeng et al., “Early truncation of the N-terminal variable region of EYA4 gene causes dominant hearing loss without cardiac phenotype,” Molecular Genetics & Genomic Medicine, no. article e1569, 2020. View at: Publisher Site | Google Scholar
  17. A. M. Bolger, M. Lohse, and B. Usadel, “Trimmomatic: a flexible trimmer for Illumina sequence data,” Bioinformatics, vol. 30, no. 15, pp. 2114–2120, 2014. View at: Publisher Site | Google Scholar
  18. H. Li, “Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM,” http://arxiv.org/abs/1303.3997. View at: Google Scholar
  19. A. Tarasov, A. J. Vilella, E. Cuppen, I. J. Nijman, and P. Prins, “Sambamba: fast processing of NGS alignment formats,” Bioinformatics, vol. 31, no. 12, pp. 2032–2034, 2015. View at: Publisher Site | Google Scholar
  20. M. A. DePristo, E. Banks, R. Poplin et al., “A framework for variation discovery and genotyping using next-generation DNA sequencing data,” Nature Genetics, vol. 43, no. 5, pp. 491–498, 2011. View at: Publisher Site | Google Scholar
  21. B. S. Pedersen, R. M. Layer, and A. R. Quinlan, “Vcfanno: fast, flexible annotation of genetic variants,” Genome Biology, vol. 17, no. 1, p. 118, 2016. View at: Publisher Site | Google Scholar
  22. The 1000 Genomes Project Consortium, “A global reference for human genetic variation,” Nature, vol. 526, no. 7571, pp. 68–74, 2015. View at: Publisher Site | Google Scholar
  23. S. T. Sherry, M. H. Ward, M. Kholodov et al., “dbSNP: the NCBI database of genetic variation,” Nucleic Acids Research, vol. 29, no. 1, pp. 308–311, 2001. View at: Publisher Site | Google Scholar
  24. M. Lek, E. A. Consortium, K. J. Karczewski et al., “Analysis of protein-coding genetic variation in 60,706 humans,” Nature, vol. 536, no. 7616, pp. 285–291, 2016. View at: Publisher Site | Google Scholar
  25. K. J. Karczewski, L. C. Francioli, G. Tiao et al., “Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes,” BioRxiv, vol. 581, no. 7809, pp. 434–443, 2019. View at: Google Scholar
  26. Q. Li and K. Wang, “InterVar: clinical interpretation of genetic variants by the 2015 ACMG-AMP guidelines,” American Journal of Human Genetics, vol. 100, no. 2, pp. 267–280, 2017. View at: Publisher Site | Google Scholar
  27. X. Liu, C. Wu, C. Li, and E. Boerwinkle, “dbNSFP v3.0: a one-stop database of functional predictions and annotations for human nonsynonymous and splice-site SNVs,” Human Mutation, vol. 37, no. 3, pp. 235–241, 2016. View at: Publisher Site | Google Scholar
  28. P. Kumar, S. Henikoff, and P. C. Ng, “Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm,” Nature Protocols, vol. 4, no. 7, pp. 1073–1081, 2009. View at: Publisher Site | Google Scholar
  29. J. M. Schwarz, D. N. Cooper, M. Schuelke, and D. Seelow, “MutationTaster2: mutation prediction for the deep-sequencing age,” Nature Methods, vol. 11, no. 4, pp. 361-362, 2014. View at: Publisher Site | Google Scholar
  30. N. M. Ioannidis, J. H. Rothstein, V. Pejaver et al., “REVEL: an ensemble method for predicting the pathogenicity of rare missense variants,” American Journal of Human Genetics, vol. 99, no. 4, pp. 877–885, 2016. View at: Publisher Site | Google Scholar
  31. R. Ghosh, S. M. Harrison, H. L. Rehm, S. E. Plon, L. G. Biesecker, and on behalf of ClinGen Sequence Variant Interpretation Working Group, “Updated recommendation for the benign stand-alone ACMG/AMP criterion,” Human Mutation, vol. 39, no. 11, pp. 1525–1530, 2018. View at: Publisher Site | Google Scholar
  32. S. Richards, on behalf of the ACMG Laboratory Quality Assurance Committee, N. Aziz 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,” Genetics in Medicine, vol. 17, no. 5, pp. 405–423, 2015. View at: Publisher Site | Google Scholar
  33. B. Adhikari and J. Cheng, “CONFOLD2: improved contact-driven ab initio protein structure modeling,” BMC Bioinformatics, vol. 19, no. 1, p. 22, 2018. View at: Publisher Site | Google Scholar
  34. M. Remmert, A. Biegert, A. Hauser, and J. Söding, “HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment,” Nature Methods, vol. 9, no. 2, pp. 173–175, 2011. View at: Google Scholar
  35. S. Seemayer, M. Gruber, and J. Söding, “CCMpred--fast and precise prediction of protein residue-residue contacts from correlated mutations,” Bioinformatics, vol. 30, no. 21, pp. 3128–3130, 2014. View at: Publisher Site | Google Scholar
  36. C. N. Magnan and P. Baldi, “SSpro/ACCpro 5: almost perfect prediction of protein secondary structure and relative solvent accessibility using profiles, machine learning and structural similarity,” Bioinformatics, vol. 30, no. 18, pp. 2592–2597, 2014. View at: Publisher Site | Google Scholar
  37. ULC C, Molecular Operating Environment MOE, Computing Group ULC, 1010 sherbrooke St.West,Suite #910,Montreal, QC,Canada,H3A 2R7, 2018.
  38. A. M. Oza, M. T. DiStefano, S. E. Hemphill et al., “Expert specification of the ACMG/AMP variant interpretation guidelines for genetic hearing loss,” Human Mutation, vol. 39, no. 11, pp. 1593–1613, 2018. View at: Publisher Site | Google Scholar
  39. M. Kobayashi, M. Miyagawa, S. Y. Nishio et al., “WFS1 mutation screening in a large series of Japanese hearing loss patients: massively parallel DNA sequencing-based analysis,” PLoS One, vol. 13, no. 3, article e0193359, 2018. View at: Publisher Site | Google Scholar
  40. K. Takeda, H. Inoue, Y. Tanizawa et al., “WFS1 (Wolfram syndrome 1) gene product: predominant subcellular localization to endoplasmic reticulum in cultured cells and neuronal expression in rat brain,” Human Molecular Genetics, vol. 10, no. 5, pp. 477–484, 2001. View at: Publisher Site | Google Scholar
  41. K. Cryns, S. Thys, L. van Laer et al., “The WFS1 gene, responsible for low frequency sensorineural hearing loss and Wolfram syndrome, is expressed in a variety of inner ear cells,” Histochemistry and Cell Biology, vol. 119, no. 3, pp. 247–256, 2003. View at: Publisher Site | Google Scholar
  42. D. Takei, H. Ishihara, S. Yamaguchi et al., “WFS1 protein modulates the free Ca(2+) concentration in the endoplasmic reticulum,” FEBS Letters, vol. 580, no. 24, pp. 5635–5640, 2006. View at: Publisher Site | Google Scholar
  43. I. Kim, W. Xu, and J. C. Reed, “Cell death and endoplasmic reticulum stress: disease relevance and therapeutic opportunities,” Nature Reviews. Drug Discovery, vol. 7, no. 12, pp. 1013–1030, 2008. View at: Publisher Site | Google Scholar
  44. L. Rigoli, P. Bramanti, C. di Bella, and F. de Luca, “Genetic and clinical aspects of Wolfram syndrome 1, a severe neurodegenerative disease,” Pediatric Research, vol. 83, no. 5, pp. 921–929, 2018. View at: Publisher Site | Google Scholar
  45. S. S. Spicer and B. A. Schulte, “Golgi-canalicular reticulum system in ion transporting fibrocytes and outer sulcus epithelium of gerbil cochlea,” The Anatomical Record, vol. 249, no. 1, pp. 117–127, 1997. View at: Publisher Site | Google Scholar
  46. M. Jeanne and D. B. Gould, “Genotype-phenotype correlations in pathology caused by collagen type IV alpha 1 and 2 mutations,” Matrix Biology, vol. 57-58, pp. 29–44, 2017. View at: Publisher Site | Google Scholar
  47. S. Hofmann, C. Philbrook, K. D. Gerbitz, and M. F. Bauer, “Wolfram syndrome: structural and functional analyses of mutant and wild-type wolframin, the WFS1 gene product,” Human Molecular Genetics, vol. 12, no. 16, pp. 2003–2012, 2003. View at: Publisher Site | Google Scholar
  48. C. Z. Seco, M. Wesdorp, I. Feenstra et al., “The diagnostic yield of whole-exome sequencing targeting a gene panel for hearing impairment in the Netherlands,” European Journal of Human Genetics, vol. 25, no. 3, pp. 308–314, 2017. View at: Publisher Site | Google Scholar
  49. H. Fukuoka, Y. Kanda, S. Ohta, and S. I. Usami, “Mutations in the _WFS1_ gene are a frequent cause of autosomal dominant nonsyndromic low-frequency hearing loss in Japanese,” Journal of Human Genetics, vol. 52, no. 6, pp. 510–515, 2007. View at: Publisher Site | Google Scholar
  50. D. Astuti, A. Sabir, P. Fulton et al., “Monogenic diabetes syndromes: locus-specific databases for Alström, Wolfram, and thiamine-responsive megaloblastic anemia,” Human Mutation, vol. 38, no. 7, pp. 764–777, 2017. View at: Publisher Site | Google Scholar
  51. A. Galvez-Ruiz, A. Galindo-Ferreiro, and P. Schatz, “Genetic testing for Wolfram syndrome mutations in a sample of 71 patients with hereditary optic neuropathy and negative genetic test results for OPA1/OPA3/LHON,” Neuroophthalmology, vol. 42, no. 2, pp. 73–82, 2018. View at: Publisher Site | Google Scholar
  52. F. Khanim, J. Kirk, F. Latif, and T. G. Barrett, “WFS1/wolframin mutations, Wolfram syndrome, and associated diseases,” Human Mutation, vol. 17, no. 5, pp. 357–367, 2001. View at: Publisher Site | Google Scholar

Copyright © 2021 Jinying Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views136
Downloads115
Citations

Related articles