A newly identified COCH mutation affects the vWFA2 domain and leads to a

relatively mild DFNA9 phenotype

Supplemental Digital Content

Supplemental Methods

Annotated variants were filtered, based on a population allele frequency of ≤ 0.05% in the gnomAD exome and genome database V.2.1 (https://gnomad.broadinstitute.org/) and our in-house exome database (containing ~15.000 alleles). Variants in coding and splice site regions (-14/+14 nt) were analyzed. Interpretation of missense variants was performed using the in silico pathogenicity prediction tools CADD-PHRED (≥15) (1), SIFT (≤0.05) (2), PolyPhen-2 (PPH2, ≥0.450) (3) and

MutationTaster (deleterious) (4). Only those called variants were considered that had a quality-by- depth >150 and that were associated with an autosomal dominant inherited type of HI. Shared variants in the MIPs data of the selected subjects were analyzed. To address the possible existence of one or more phenocopies, shared variants between two subjects (IV.1 with III.9, IV.1 with III.13 and III.9 with

III.13) were also analyzed separately. Variants were considered if a pathogenic effect was predicted by at least two different tools. Potential effects on splicing of missense and synonymous variants were evaluated using five algorithms (SpliceSiteFinder, MaxEntScan, NNSPLICE, GeneSplicer and Human

Splicing Finder) via the AlamutVisual software (V.2.10, Interactive Biosoftware). A change of ≥5% in splice site scores predicted by at least two algorithms was considered significant.

Table S1. Hearing loss covered by MIPS panel

ACTG1 GRM7 PTPRQ

ADCY1 GRM8 RDX

BDP1 GRXCR1 RIPOR2

BSND GRXCR2 SERPINB6

CABP2 HGF SIX1

CCDC50 ILDR1 SLC17A8

CDH23 KARS SLC26A4

CEACAM16 KCNQ4 SLC26A5

CIB2 LHFPL5 SMPX

CLDN14 LOXHD1 STRC

CLIC5 LRTOMT SYNE4

COCH MARVELD2 TBC1D24

COL11A2 MIR96 TECTA

COL4A6 MSRB3 TJP2

CRYM MYH14 TMC1

DCDC2 MYH9 TMEM132E

GSDME MYO15A TMIE

DFNB31 MYO3A TMPRSS3

DFNB59 MYO6 TNC

DIABLO MYO7A TPRN

DIAPH1 NAT2 TRIOBP

ELMOD3 OSBPL2 TSPEAR

EPS8 OTOA USH1C

ESPN OTOF USH1G

ESRRB OTOG

EYA4 OTOGL

GIPC3 P2RX2

GJB2 PCDH15

GJB3 PNPT1

GJB6 POU3F4

GPSM2 POU4F3

GRHL2 PRPS1 Table S2. Detailed information on COCH, cochlin and pathogenicity predictions on the identified COCH variant.

Genomic RefSeq RefSeq Genomic cDNA GnomAD_ GnomAD_ CADD_ Mutation positions COCH cochlin SIFT PPH2 PhyloP GERP REVEL symbol change change change E AF (%) G AF (%) PHRED Taster COCH transcript protein

chr14:31,34 NM_004 NP_004 g.31,355,3 Disease COCH 3,741- c.1312C>T p.Arg438Cys 3.98E-06 - 23.9 0.05 0.974 1.388 4.06 0.697 086.2 077.1 53C>T causing 31,359,822

Scores that meet the thresholds for a possibly deleterious effect are indicated in red. Thresholds are: CADD-PHRED (≥15)(1), SIFT (<0.05) (2), PolyPhen-2

(≥0.450)(3) and MutationTaster (disease causing)(4), PhyloP (≥2.7)(5), GERP(≥2) (6) and REVEL(>0.5)(7). Genomic positions COCH according to GRCh37/hg19;

GnomAD_E AF and GnomAD_G AF, allele frequencies in gnomAD total exome database and genome database, respectively; CADD_PHRED, Combined

Annotation Dependent Depletion PHRED score; SIFT, Scale-Invariant Feature Transform; PPH2, PolyPhen-2 score; PhyloP, phylogenetic P-values; GERP,

Genomic Evolutionary Rate Profiling; REVEL, Rare Exome Variant Ensemble Learner -, frequency not available.

1. Kircher M, Witten DM, Jain P, O'Roak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nature genetics. 2014;46(3):310-5. 2. Vaser R, Adusumalli S, Leng SN, Sikic M, Ng PC. SIFT missense predictions for genomes. Nature Protocols. 2015;11:1. 3. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nature methods. 2010;7(4):248-9. 4. Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nature Methods. 2014;11:361. 5. Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A. Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010;20(1):110-21. 6. Cooper GM, Stone EA, Asimenos G, Program NCS, Green ED, Batzoglou S, et al. Distribution and intensity of constraint in mammalian genomic sequence. Genome Res. 2005;15(7):901-13. 7. Ioannidis NM, Rothstein JH, Pejaver V, Middha S, McDonnell SK, Baheti S, et al. REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants. Am J Hum Genet. 2016;99(4):877-85.

Page 1 of 2

Figure S1 Linkage analysis

The blue lines define the customary LOD score range. The green line demonstrates the maximum LOD score measured genome-wide in this experiment. It appears after all computations are completed.

Chromosome 1 2 Chromosome 3

Chromosome 4 Chromosome 5

Chromosome 7 Chromosome 8 Chromosome 9

Chromosome 10 Chromosome 11 Chromosome 12

Chromosome 13 Chromosome 14 Chromosome 15 Page 2 of 2

Chromosome 16 Chromosome 17 Chromosome 18

Chromosome 19 Chromosome 20 Chromosome 21

Chromosome 22 Chromosome X W08-2035

I.1 I.2

84 25 II.1 II.2 II.3 II.4 II.5 II.6 II.7 II.9II.8 47 27 15 54 17 67 ? 45 44 77 35 43 49 51 65 64 70 26 ? III.1 III.2 III.3 III.5III.4 III.7III.6 III.9III.8 III.10 III.11 III.12 III.13 III.14 III.15 III.16 III.17 25 34 5 5 26 2 4 2 1 3 4 21 46 21 22 24 22 51 11 1 1 51 5 1 5 3 1 2 33 11 33 51 51 54 66 63 6 6 63 6 4 6 3 8 3 53 66 53 66 66 66

33 35 25 42 26 43 18 33 31 ? Marker Position (cM) IV.1 IV.2 IV.3 IV.4 IV.5 IV.6 IV.7 IV.8 IV.9 IV.10 D14S1021 23 6 2 5 2 2 4 25 1 3 24 26 24 2 2 31.13 COCH D14S54 31.75 51 1 5 1 5 1 1 51 3 1 52 31 31 1 5 D14S257 31.75 66 3 6 6 6 6 4 63 3 8 63 56 56 2 6

Figure S2. Results of VNTR marker genotyping, affected individual III.13 does not carry the disease allele. The shared haplotype is boxed.

The c.1312C>T COCH variant is located between the markers D14S1021 and D14S54. cM, centimorgan.