Novel autosomal recessive non-syndromic hearing impairment locus (DFNB71) maps to chromosome 8p22-21.3.

A novel ARNSHI (autosomal recessive non-syndromic hearing impairment) locus, DFNB71, was localized to 8p22-21.3. To map the locus, a whole-genome scan was carried out using DNA samples from a consanguineous seven-generational Pakistani family with profound prelingual ARNSHI. A maximum multipoint LOD score of 4.2 occurred at marker D8S261. The DFNB71 region of homozygosity and 3-unit support interval is flanked by markers D8S1130 and D8S1786. This region has a genetic distance of 19.1 cM and contains 10.6 Mb of sequence.

Introduction

The incidence of profound hearing impairment (HI) is estimated at 1–2 per 1,000 births and ~60% of the congenital cases can be attributed to genetic factors (1, 2). Of the genetic cases, ~30% are syndromic with the remaining cases display non-syndromic HI (NSHI) where the HI phenotype is usually sensorineural. A large number of genes attributable to NSHI have been mapped and identified. For autosomal recessive (AR) NSHI > 60 genes have been localized and 28 genes have been identified (3). Although this magnitude of progress is remarkable and significant advances have been made, many more NSHI loci remain to be discovered. In this article, a new locus, DFNB71, is shown to segregate in a seven-generational consanguineous Pakistani family with ARNSHI. The hearing impaired family members in this pedigree present with profound prelingual HI. The DFNB71 locus maps to a 10.6 Mb region on chromosome 8p22-21.3.

Materials and Methods

Family History

With the prior approval from the Institutional Review Board of Quaid-I-Azam University, Islamabad, Pakistan and Baylor College of Medicine and Affiliated Hospitals the study was initiated and informed consent was obtained from all family members who participated in the study. Pedigree 4216 was constructed on the basis of information collected from the family members (Figure 1). Hearing impaired and non-hearing impaired family members underwent audiometric testing. Hearing impaired family members have prelingual profound HI which affects all frequencies. Figure 2 displays the audiogram for a hearing impaired family member (individual 34) that displays bilateral, profound HI that affects all frequencies and the audiogram for a family member without hearing impairment (individual 38) whose hearing is within the normal range (Figure 2). The hearing impaired members of pedigree 4216 use sign language as their primary means of communication. The family members with hearing impairment also underwent a physical examination for defects in ear morphology, vision, mental retardation and other clinical features that could indicate that the HI is syndromic. There was no evidence in this kindred that the HI belongs to a syndrome or that there is gross vestibular involvement.

Figure 1

Drawing of pedigree 4216. Black symbols represent individuals with hearing impairment and clear symbols represent unaffected individuals. Haplotypes are shown beneath each genotyped individual. The DFNB71 haplotype is enclosed in a rectangle.

Figure 2

Audiograms for a hearing impaired (panel a) and unimpaired family member (panel b) from pedigree 4216. Circles and crosses represent air conduction for the right and left ear, respectively. Panel (a) displays the audiogram for individual 34 at 24 years of age and panel (b) displays the audiogram from individual 38 at 17 years of age.

Genotyping

Venous blood samples were obtained from 11 family members including 5 individuals with profound prelingual HI. Genomic DNA was extracted from whole blood following a standard protocol (4) A genome scan was carried out on all DNA samples at the Center for Inherited Disease Research (CIDR). A total of 396 fluorescently labeled short tandem repeat (STR) markers were genotyped. These markers are spaced ~10 cM apart and are located on the 22 autosomes and the X and Y chromosomes. For fine mapping, microsatellite markers were PCR amplified according to standard procedure in a total volume of 25 μl with 40 ng genomic DNA, 240 nM dNTP, and 1 unit of Taq DNA Polymerase (Fermentas Life Sciences, Burlington, ON, Canada) in GeneAmp® PCR System 9700 (Applied Biosystems, Applera Corp., Foster City, CA, USA). PCR products were resolved on 8% non-denaturing polyacrylamide gel and genotypes were assigned by visual inspection.

Linkage Analysis

PEDCHECK(5) was used to evaluate the genotype data for Mendelian inconsistency and MERLIN (6) was used to assess the data for occurrence of unlikely genotypes, which are most likely due to genotyping error. Two-point linkage analysis was carried out using the MLINK program of the FASTLINK computer package (7). LOD scores were calculated using a fully penetrant autosomal recessive model with a disease allele frequency of 0.001. Multipoint linkage analysis was performed using SUPERLINK (8, 9) and SIMWALK2 (10, 11) utilizing map distances from the Rutgers combined linkage-physical map of the human genome (12). The allele frequencies for the genome scan markers were estimated from the data by means of both observed and reconstructed genotypes of founders from this pedigree and 35 additional families from Pakistan that underwent a genome scan at CIDR at the same time. In order to fine map the DFNB71 loci 12 additional markers were genotyped in the 8p region of which 8 were informative (D8S640, GATA151F02, D8S26, D8S261, KW218, D8S1786, D8S136 and D8S1989). Two-point linkage analysis was carried for the fine mapping markers and multipoint linkage analysis was performed again using both the genome scan and fine mapping markers. Equal allele frequencies were used for the fine mapping markers, since they were only genotyped in family 4216. A sensitivity analysis was carried out by performing multipoint linkage analysis and varying the allele frequency for the alleles segregating with the HI phenotype from 0.2 to 0.6 for the fine mapping markers, because analyses carried out with potentially incorrect allele frequencies (e.g. equal allele frequencies) can increase type I error (13). Haplotypes in the 8p22-21.3 region were reconstructed using SIMWALK2.

Candidate Gene Sequencing

Primers were designed for the exons of the hairless (HR) gene (MIM 602302) isoform a (NM_005144) and isoform b (NM_018411), the platelet-derived growth factor receptor-like (PDGFL) gene (MIM 604584) and the tumor suppressor candidate 3 (TUSC3) gene (MIM 601385) isoform a (NM_006765) and isoform b (NM_178234) using Primer3 software (14). DNA from one unaffected and two HI pedigree members were diluted to 5 μg/nl, amplified by PCR under standard conditions and purified with ExoSAP-IT ® (USB Corp., Cleveland, OH, USA). Sequencing was performed with the BigDye ® Terminator v3.1 Cycle Sequencing Kit together with an Applied Biosystems 3700 DNA Analyzer (Applera Corp., Foster City). Sequence variants were identified via Sequencher™ Version 4.1.4 software (Gene Codes Corp., Ann Arbor, MI, USA).

Results

Analysis of the genome scan markers revealed a maximum multipoint LOD score of 3.1 which occurred at marker D8S1106. The region of homozygosity for the genome scan markers is flanked by markers D8S1130 and D8S1771 which spans a 23.4 cM region according to the Rutgers combined linkage-physical map of the human genome (12) and contains 13.6 Mb of sequence (15). Analysis of the fine mapping markers produced a maximum two-point LOD score of 3.04 (θ=0.00) at marker D8S640 (Table 1) and the maximum multipoint LOD score (4.2) occurred at marker D8S261. The 3-unit support interval spanned from markers D8S1130 to D8S1786. Haplotypes were constructed to determine the critical recombination events. The 8p22-21.3 haplotype for pedigree 4216 is presented in Figure 1. All HI family members are homozygous for a common haplotype, which is also flanked by markers D8S1106 and D8S1786, and agrees with the 3-unit support interval (Table 1). This region is 19.1 cM according to the Rutgers combined linkage-physical map of the human genome and contains 10.6 Mb of sequence. Sensitivity analysis revealed that varying the allele frequencies caused the LOD score to fluctuate. But, a LOD score of greater than 3.3 was maintained within the region of homozygosity for all allele frequencies tested. Sequencing the exons and promoter regions of three candidate genes, TUSC3 (MIM: 601385), PDGFRL (MIM: 604584) and HR (MIM: 602302) using DNA from two HI pedigree members did not reveal any potentially functional sequence variants.

Table 1

Two-point LOD score results between DFNB71 and chromosome 8 markers.

	Genetic map	Physical map	Two Point LOD score at θ =
Marker	PositionI	PositionII	0	0.01	0.05	0.1	0.2	0.3	0.4
D8S1469	19.38	9,127,124	− ∞	−2.76	−0.87	−0.22	0.15	0.13	0.01
D8S1130	22.31	11,871,501	− ∞	−2.53	−0.65	−0.01	0.36	0.31	0.13
D8S1106	24.76	12,880,323	1.58	1.54	1.36	1.15	0.77	0.43	0.16
D8S640	24.76	12,979,668	3.04	2.96	2.65	2.27	1.54	0.87	0.32
GATA151F02	26.35	13,296,758	1.16	1.12	0.99	0.83	0.55	0.30	0.10
D8S26	27.04	13,750,569	1.69	1.65	1.50	1.30	0.92	0.56	0.23
D8S261	32.78	17,870,469	1.69	1.65	1.50	1.30	0.92	0.56	0.23
D8S560	39.68	21,646,126	1.76	1.72	1.56	1.36	0.97	0.60	0.26
KW218	40.71	22,024,717	2.38	2.32	2.08	1.78	1.21	0.68	0.21
D8S1786	41.41	22,489,213	1.62	1.59	1.46	1.29	0.94	0.59	0.25
D8S136	41.41	22,489,339	1.39	1.37	1.25	1.11	0.81	0.51	0.21
D8S1989	45.23	24,697,056	− ∞	−0.51	0.63	0.91	0.85	0.53	0.16
D8S1771	45.71	25,497,030	−2.46	−0.29	0.31	0.48	0.47	0.31	0.10
D8S1477	55.27	32,186,709	− ∞	−4.14	−2.11	−1.28	−0.56	−0.23	−0.06

Markers in italic flank the region of homozygosity and the 3-unit support interval. Genome scan markers are displayed in bold.

The March 2006 human reference sequence (NCBI Build 36.1, Kent et al. 2002)

The Rutgers Combined Linkage-Physical Map of The Human Genome Build 36 (Matise et al. 2007)

Discussion

There are 54 known genes within the 10.6 Mb region of the DFNB71 locus. Expression information on three of these genes indicated that they may be important in hearing impairment physiology. TUSC3 was originally associated with prostate cancer (16) but is of interest as a candidate gene for HI due to its expression in the inner ear (17–19). The platelet-derived growth factor receptor polypeptides alpha, beta, c and d are expressed in the inner ear (17–19), although none of these genes reside in the DFNB71 interval, a related gene the platelet-derived growth factor receptor like (PDGFRL) gene is within the DFNB71 region. HR is a homolog of the mouse HR which displays expression in the mouse organ of corti (20). Sequencing of the exons and promoter regions of TUSC3, PDGRFL and HR and their isoforms did not reveal any potential functional variants and therefore their involvement in the etiology of HI due to DFNB71 is not supported. Eventual identification the gene responsible for DFNB71 will enhance our knowledge of the molecular mechanisms behind hearing and HI.