# Identification of Missense Mutation (I12T) in the BSND Gene and Bioinformatics Analysis

Authors: Hina Iqbal; Tayyba Sarfaraz; Farida Anjum; Zubair Anwar; Asif Mir
Journal: Journal of Biomedicine and Biotechnology
 (2011)

---

## Abstract

Nonsyndromic hearing loss is a paradigm of genetic heterogeneity with 85 loci and 39 nuclear disease genes reported so far. Mutations of BSND have been shown to cause Bartter syndrome type IV, characterized by significant renal abnormalities and deafness and nonsyndromic nearing loss. We studied a Pakistani consanguineous family. Clinical examinations of affected individuals did not reveal the presence of any associated signs, which are hallmarks of the Bartter syndrome type IV. Linkage analysis identified an area of 18.36 Mb shared by all affected individuals between markers D1S2706 and D1S1596. A maximum two-point LOD score of 2.55 with markers D1S2700 and multipoint LOD score of 3.42 with marker D1S1661 were obtained. BSND mutation, that is, p.I12T, cosegregated in all extant members of our pedigree. BSND mutations can cause nonsyndromic hearing loss, and it is a second report for this mutation. The respected protein, that is, BSND, was first modeled, and then, the identified mutation was further analyzed by using different bioinformatics tools; finally, this protein and its mutant was docked with CLCNKB  and REN, interactions of BSND, respectively.


---

## Body

## 1. Introduction

Deafness or hearing loss can be due to genetic or environmental causes or a combination of both. The genetic hearing loss is classified as syndromic or nonsyndromic. Among the many disorders classified as syndromic hearing loss, the pathology varies widely, but in nonsyndromic hearing loss, the defect is generally sensorineural. Seventy percent of deafness has genetic causes and is classified as nonsyndromic. Autosomal recessive nonsyndromic deafness is genetically heterogeneous and is the most common form of inherited hearing loss. Autosomal recessive genes are responsible for about 77% of the cases of hereditary nonsyndromic deafness, with over 85 loci and 21 different genes identified to date (http://hereditaryhearingloss.org/). The high degree of genetic heterogeneity of deafness reflects the great diversity of specialized proteins that are required to make sense of sound, and continuing discovery of common and rare mutations associated with deafness in humans has provided many serendipitous points of entry into the biology of hearing.Bartter's syndrome (BS) is characterized by hypokalemic, hypochloremic metabolic alkalosis with normal or low blood pressure despite high plasma renin activity and serum aldosterone. The inheritance pattern is autosomal recessive. Antenatal BS with bilateral sensorineural deafness (BSND) was first described in children born to a consanguineous couple from a Bedouin family of Southern Israel [1]. BS type IV, BSND variant, occurs because of a mutation in the BSND gene on chromosome 1p31 coding for protein “barttin” which forms the β subunit of CICKb and CICKa channels located on the basolateral membrane of TAL and inner ear epithelium [2]. The clinical features of Bartter syndrome type IV include sensorineural deafness and peculiar facies, distinguished by the triangular face, large eyes, and protruding ears (1).Mutation in BSND gene is analyzed by model prediction using SAM_T08 server. Comparison of the protein sequence of BSND with closely related species takes place and find out conserved regions. Other bioinformatics analysis such as protein modeling and proteins docking are also carried out for more bioinformation.


## 2. Materials and Methods

Before the onset of the study, approval from the Institutional Review Board, Islamabad, and informed consent was obtained from all individuals. Clinical examination was performed at hospitals of the respective area. Audiometry was performed on selected individuals to detect the level of hearing loss. The families were visited at their places of residence to generate pedigrees [3] and collect other relevant information. Blood samples from available affected and normal individuals of each family were collected for DNA extraction.


### 2.1. Genotyping and Linkage Studies

High molecular weight DNA was extracted from leukocytes following the standard method as described by Sambrook et al. [4]. Genomic DNA was quantified by spectrophotometer readings at OD260 and diluted to 40 ng/uL for amplification by polymerase chain reaction (PCR).To elucidate the gene defect in the family presented here, an initial search for linkage was carried out by using polymorphic microsatellite markers mapped within autosomal recessive nonsyndromic deafness loci listed on the hereditary hearing loss homepage [5]. Genome-wide screening was conducted with microsatellite markers (Linkage Mapping Set 10, Invitrogen, USA). Two-point linkage analysis was carried out using MLINK of the FASTLINK computer package [6]. Multipoint linkage analysis was performed using ALLEGRO [7]. For the analysis, an autosomal recessive mode of inheritance with complete penetrance and a disease allele frequency of 0.001 were assumed. Equal allele frequencies were used in the analysis.  However, since it is well known that using allele frequencies, which are too low, can lead to false positive results, a sensitivity analysis was performed. Haplotypes were constructed using SIMWALK2 [8, 9].


### 2.2. Mutation Screening

To screen for mutation in the BSND gene, exons and splice junction sites were PCR amplified from genomic DNA using the four primer sets. The purified PCR products were subjected to cycle sequencing using big dye terminator V 3.1 ready reaction mix and sequencing buffer (PE Applied Biosystems, Foster City, Calif, USA). The sequencing products were purified to remove unincorporated nucleotides and primers with Centriflex TM Gel Filtration Cartridge (Edge Biosystems, Gaithersburg, Md, USA). These purified products were resuspended in 10μL of TSR (Template Suppression Reagent) and were placed in 0.5 mL septa tubes to be directly sequenced in an ABI Prism 310 Automated Sequencer (PE Applied Biosystems, Foster City, Calif, USA). Chromatograms from normal and affected individuals were compared with the corresponding control gene sequences from NCBI (National Center for Biotechnology Information) database to identify the aberrant nucleotide base-pair change. (http://www.ncbi.nlm.nih.gov/).


### 2.3. Protein Modeling, Docking, and Phylogenetics Analysis

SAM_T08 [10–14] was used to model the BSND protein and its interacting protein with highest confidence score (i.e., 0.990) CLCNKB. For their further evaluation, RAMPAGE server [15] was considered. ClustalW program was subjected to analyze the phylogenetics relationship of different animals having BSND protein. Conserved amino acids were predicted by using UniportKB tool available on Expasy server.Identified BSND mutation was generated by using Swisspdb server. Protein interactor of BSND was  found out by STRING server, an online database of known and predicted protein interactions which includes direct (physical) and indirect (functional) associations [16]. Protein docking of both the normal and the mutant BSND was carried out using HEX tool with its interactive proteins CLCNKB and REN.


## 2.1. Genotyping and Linkage Studies

High molecular weight DNA was extracted from leukocytes following the standard method as described by Sambrook et al. [4]. Genomic DNA was quantified by spectrophotometer readings at OD260 and diluted to 40 ng/uL for amplification by polymerase chain reaction (PCR).To elucidate the gene defect in the family presented here, an initial search for linkage was carried out by using polymorphic microsatellite markers mapped within autosomal recessive nonsyndromic deafness loci listed on the hereditary hearing loss homepage [5]. Genome-wide screening was conducted with microsatellite markers (Linkage Mapping Set 10, Invitrogen, USA). Two-point linkage analysis was carried out using MLINK of the FASTLINK computer package [6]. Multipoint linkage analysis was performed using ALLEGRO [7]. For the analysis, an autosomal recessive mode of inheritance with complete penetrance and a disease allele frequency of 0.001 were assumed. Equal allele frequencies were used in the analysis.  However, since it is well known that using allele frequencies, which are too low, can lead to false positive results, a sensitivity analysis was performed. Haplotypes were constructed using SIMWALK2 [8, 9].


## 2.2. Mutation Screening

To screen for mutation in the BSND gene, exons and splice junction sites were PCR amplified from genomic DNA using the four primer sets. The purified PCR products were subjected to cycle sequencing using big dye terminator V 3.1 ready reaction mix and sequencing buffer (PE Applied Biosystems, Foster City, Calif, USA). The sequencing products were purified to remove unincorporated nucleotides and primers with Centriflex TM Gel Filtration Cartridge (Edge Biosystems, Gaithersburg, Md, USA). These purified products were resuspended in 10μL of TSR (Template Suppression Reagent) and were placed in 0.5 mL septa tubes to be directly sequenced in an ABI Prism 310 Automated Sequencer (PE Applied Biosystems, Foster City, Calif, USA). Chromatograms from normal and affected individuals were compared with the corresponding control gene sequences from NCBI (National Center for Biotechnology Information) database to identify the aberrant nucleotide base-pair change. (http://www.ncbi.nlm.nih.gov/).


## 2.3. Protein Modeling, Docking, and Phylogenetics Analysis

SAM_T08 [10–14] was used to model the BSND protein and its interacting protein with highest confidence score (i.e., 0.990) CLCNKB. For their further evaluation, RAMPAGE server [15] was considered. ClustalW program was subjected to analyze the phylogenetics relationship of different animals having BSND protein. Conserved amino acids were predicted by using UniportKB tool available on Expasy server.Identified BSND mutation was generated by using Swisspdb server. Protein interactor of BSND was  found out by STRING server, an online database of known and predicted protein interactions which includes direct (physical) and indirect (functional) associations [16]. Protein docking of both the normal and the mutant BSND was carried out using HEX tool with its interactive proteins CLCNKB and REN.


## 3. Results

### 3.1. Clinical Assessments

Ages of affected individuals varied between 12 and 65 years at the time of study. Clinical examination of the hearing impaired individuals did not reveal the presence of signs or symptoms, which are hallmarks of the Bartter syndrome type IV. Urinary and blood biochemistry was tested at the Islamabad Diagnostic Centre, including testing for liver function, renal function, electrolytes, and hematology. Blood hematology and biochemistry markers were within the normal range, except that renin level was at borderline in one of the affected subject (VI: 5, 2.53 mL/h). Detailed clinical and biochemical features of normal and affected members are mentioned in Table1. The height and weight were within normal range for the local population. An audiometric evaluation of selected affected members by measuring the threshold of hearing at 250–8000 Hz for pure-tune air conduction and bone conduction showed severe hearing loss across all frequencies.Table 1
Clinical and biochemical features of normal and affected members of studied family.
FeaturesReference rangesNormal membersAffected members from the family with hearing lossV : 3V : 5V : 1VI : 3VI : 5Age (y)3035451210Other symptoms of BSIVnonoNoNoNoSd·Na (mmol/L)136–148141140135137134S·K (mmol/L)3.6–5.03.94.63.43.53.5S·Cl (mmol/L)104–11410510797102101S·Mg (mg/dl)1.9–2.52.32.22.01.92.1S·HCO3 (mmol/L)17.5–27.530.027.530.030.231.0S·Ca (mg/dl)8.6–10.58.79.59.29.810.0S·Creatinine (mg/dL)0.85–1.350.70.80.60.50.4Pe Renin (ng/ml/hr)0.15–2.330.310.62.132.032.53S·Aldosterone (ng/dL)1–6;f 4–31>1.7f13.010.311.96.3S·Osmolality (mosm/kg)275–295293289281290287Ug Na (mmol/L)30–15017.5154.639.116.913.3U·K (mmol/L)20–673.491.3112.138.434.12U·Ca (mg/dL)6.016.015.311.313.0U·Mg (mg/dL)3.614.43.12.41.12U·Cl (mmol/L)46–16816.3174.346.5619.210.5Uosmolality (mosm/kg)50–14007490317910556NephrocacinosisSAbsentAbsentAbsentAbsentAbsentRenal sonography ruled out the presence of nephrocalcinosis or metabolic alkalosis in these affected individuals. Nor have the affected individuals displayed problems with polydipsia, polyuria, nocturnal enuresis, and hypocalciuria.


### 3.2. Genotyping Results

Analysis of the results obtained from genome search performed on the ten individuals of the family (III-1, III-3, IV-1, V-1, V-2, V-4, V-5, VI-1, VI-3, and VI-5) (Figure1) identified an area of interest on chromosome 1. Two-point analysis generated LOD score of 2.027 at marker D1S193 (73.21 cm) and 2.05 at marker D1S3462 (247.23 cm) on chromosome 1 (Table 2). In order to test linkage to these two regions, additional markers located in the vicinity of D1S3462 and D1S193 were chosen from Marshfield map [17] and genotyped in all the ten family members. The maximum two-point LOD score was increased to 2.55 with marker D1S2700. Multipoint linkage analysis for the family derived a maximum LOD score of 3.42 at marker D1S1661 (Table 3). The three-unit multipoint support interval contained a 10.48 Mb region, which span from markers D1S3721 to D1S2690. Haplotypes analysis delimited the centromeric boundary defined by a recombination between markers D1S1596 and D1S2770 in individual V-1. Recombination event between D1S2706 and D1S3721, defined the telomeric boundary in the same individual (V-1). Therefore, the minimum critical region of 18.36 Mb identified for disease locus and shared by all the affected individuals between the markers D1S2706 and D1S1596. The critical interval (17.27 cm) identified in this Pakistani family overlaps with the critical region to which DFNA2 was mapped on 1p34 in a large Indonesian family with autosomal dominant, progressive and sensorineural hearing loss [18]. The DFNA2 locus maps between markers D1S255 (65.47 cm) and D1S211 (73.81 cm), and thus shares a region of 2.68 cm with DFNB interval identified in our family. Riazuddin et al. [19] identified 1.5 Mb in four kindred segregating nonsyndromic deafness at the chromosome 1.Table 2
Two-point LOD score results between the locus identified in family and chromosome 1p34.2–p32.1 markers.
MarkerdeCode map postion1Marshfield map postion2Physical map postion3LOD score ATθ=0.000.010.050.10.20.3D1S47257.9565.4737,102,828−999.99−0.68−0.050.140.200.13D1S18660.1467.2237,900,316−999.99−.069−0.100.050.100.06D1S432—69.8639,046,164−1.334−0.69−0.18−0.020.030.01D1S159863.5170.4140,011,183−1.544−0.61−0.090.030.040.00D1S270665.3171.1340,741,361−999.99−1.34−0.180150.240.14D1S372165.8772.5941,395,5602.0271.971.751.491.490.56D1S19366.6473.2142,688,3692.0271.971.751.490.990.56D1S271367.8773.8144,182,5691.7261.671.471.240.800.43D1S213471.2975.6647,993,2411.6901.661.521.320.910.52D1S166172.4278.2551,142,5511.2281.191.060.900.600.33D1S265277.1580.7755,178,8521.3191.281.130.960.620.34D1S269079.6383.0756,867,5602.0271.971.751.490.990.56D1S270083.0787.3158,828,9552.5562.482.211.881.230.66D1S277085.0288.4059,442,255−999.99−0.260.260.350.290.17D1S1596—89.4959,670,589−1.55−0.52−0.000.100.080.01D1S273785.4790.5861,038,037−999.99−1.100.070.400.460.29D1S278887.2593.8661,568,857−999.99−2.16−0.90−0.46−0.14−0.03D1S23088.1295.3162,314,306−999.99−2.16−0.90−0.46−0.14−0.03D1S164899.62101.4873,125,563−999.99−2.27−1.00−0.55−0.20−0.06D1S166599.62102.0273,941,402−999.99−2.34−1.04−0.57−0.20−0.061Sex-average Kosambi cM map distances from the deCode genetic map [7].2Sex-average Kosambi cM map distances from the Marshfield genetic map [17].3Sequence-based physical map distance in bases according to the Human Genome Project-Santa Cruz (http://genome.ucsc.edu/, Assembly: May 2004).Table 3
Multipoint LOD score results between linked locus identified in family and chromosome 1p34.2–p32.1 markers.
MarkerPhysical positionLocationLODATA79C10—0.0000−2.0630D1S47237,102,8283.0940nanD1S18637,900,3165.1350nanD1S43239,046,1648.2290−2.4708D1S159840,011,1839.2390−0.0755D1S270640,741,36110.2490−8.2103D1S372141,395,56011.25903.0063D1S19342,688,36912.26903.2113D1S271344,182,56914.31003.3672D1S213447,993,24116.35103.4190D1S166151,142,55118.39203.4201D1S265255,178,85221.48603.3868D1S269056,867,56024.58003.2867D1S270058,828,95528.74902.8398D1S277059,442,25529.75901.3404D1S159659,670,58930.7690nanD1S273761,038,03731.7800−infD1S278861,568,85734.8730nanD1S23062,314,30636.9140−infD1S164873,125,56343.3060−infD1S166573,941,40244.3160nanGATA152F05L85,238,77045.3260nanFigure 1
Pedigree of family with nonsyndromic autosomal recessive hearing loss associated haplotypes. Haplotypes for the most closely linked STRPs are shown below each symbol. The genetic map distances according to the Marshfield genetic map in centimorgans (cM) are shown in parenthesis next to the marker name. Arrows adjacent to the haplotypes indicate key recombination events. The alleles are denoted 1–3 according to their allele size.


### 3.3. Mutation Analysis of Candidate Genes

Through a database search, we identified several genes mapping between the linkage interval in the family (Human Genome Project-Santa Cruz;http://genome.ucsc.edu/, May 2004).  Among these, KCNQ4 (MIM 603537), CLDN19, FOXE3 (MIM 601094), FOXD2 (MIM 602211), TSPAN1, and BSND (MIM 606412) were plausible candidate candidates in the interval of our family. We started sequencing of these and completely sequenced two genes (KCNQ4 and CLDIN19), but we found no mutation in both. Recent studies of Riazuddin et al. [19] showed that mutation in BSND can cause nonsyndromic deafness and is the molecular basis of DFNB73 locus in four Pakistani families segregating nonsyndromic deafness [19]. As the critical interval identified in the family contains potential candidate gene barttin (BSND, MIM 606412) as well; therefore, it was sequenced to search for the mutation. Sequence analysis of exon 1 of the BSND gene in our family also revealed missense mutation involving T to C transition at nucleotide number 35 (35T > C), which resulted in substitution of isoleucine to threonine at amino acid position 12 (I12T) (Figure 2). BSND mutation; that is, p.I12T cosegregated in all extant members of the pedigree. This mutation was present in the heterozygous state in obligate carriers within the family.Representative chromatogram generated by Big Dye Terminator, sequencing of translated exon1 of BSND gene indicating a T to C transition at nucleotide 35 (35T > C) from an affected female individual (VI-1) (a), a heterozygous carrier (VI-1) (b) and a normal male individual (V-2) (c) of the family. The T to C transition at nucleotide position 35 results in isoleucine to threonine change at amino acid 12 (I12T). Arrows indicate the nucleotide change in the sequence.
(a)(b)(c)


### 3.4. Bioinformatics Analysis

Models of BSND CLCNKB and REN proteins were generated by SAM_T08 server and visualized by Rasmol (Figure3(a)–3(c)). Rampage values for BSND were number of residues in favored region 95.6%, the number of residues in allowed region 3.1%, and the number of residues in outlier region 1.4% (Figure 4(a)); however, the CLCNKB and REN values were little varying. Its values were the number of residues in favored region:  94.8%, 97%, the number of residues in allowed region:  4.2%, 3%, and the number of residues in outlier region is 0.9% for CLCNKB and 0% for REN (Figure 4(c)).Protein models of BSND (a) CLCNKB (b), and REN (c) generated by SAM_T08 and visualized by Rasmol (display = ribbon; colors = group).
(a)(b)(c)Rampage results of residues for protein model of human BSND protein (a) CLCNKB (b) and REN (c).
(a)(b)(c)Comparison of conserved amino acid of human BSND with closely related species like mouse, rat, and rabbit shows that isoleucine at position 12 is highly conserved, so mutation at this point can be  significant (Figure5(a)). The evolutionary relationship between human, mouse, rat, and rabbit were evaluated by phylogenetic analysis. The analysis represents that the BSND of mouse and rat is closely related to each other and shows homology with human and rabbit (Figure 5(b)).(a) Comparison of conserved regions of BSND protein between human, mouse, rat, and rabbit by using UniportKB ClustalW. Dark region shows conserved amino acid sequence between species. (b) Phylogenetic analysis of BSND protein. It shows evolutionary relationship between human, mouse, rat, and rabbit species.
(a)(b)Our mutation (I12T) of BSND was generated by Swisspdb server and was analyzed by Viewerlite. This substitution of isoleucine with threonine altered the bonding capacity with other side chains as well as other interactors (Figure6(a)-(b)).Figure 6
(a) Normal and mutated model (Il12T) of BSND are superimposed and visulized in Viewerlite. (b) Focused normal model is yellow highlighted and displayed in stick format, and mutated is displayed in ball and stick format.According to STRING database, the physical and functional interaction of BSND was very much significant with CLCNKB protein. Clinical assessment showed the reduced chlorine level in all the affected individuals. According to further clinical assessments, renin level was at border range in one affected individual, so REN protein was also selected to check the effect of the respected mutation on this protein. However, other interactors are COL7A1, ARL2, CLCN5, KCNJ1, CASR, SLC12A1, SLC12A3, and CLCNKA (Figure7). Both the proteins were docked with HEX v 6.1, taking BSND as receptor and CLCNKB and REN as a ligand. The total energy value calculated for BSND-CLCNKB molecules was −229094.6   and BSND-REN was −1145.7. However, the total energy value for the mutant (I12T) on docking with CLCNKB was maximized to −218782.3   and with REN was −1149.8. This values fluctuation ultimately results in reduced binding affinity with CLCNKB and increased binding affinity with REN and so may be the case with other interacting proteins.Figure 7
The interactions of BSND with other functionally similar proteins retrieved by STRING online database. CLCNKB  showed  highest confidence score, that is, 0.990 with BSND.


## 3.1. Clinical Assessments

Ages of affected individuals varied between 12 and 65 years at the time of study. Clinical examination of the hearing impaired individuals did not reveal the presence of signs or symptoms, which are hallmarks of the Bartter syndrome type IV. Urinary and blood biochemistry was tested at the Islamabad Diagnostic Centre, including testing for liver function, renal function, electrolytes, and hematology. Blood hematology and biochemistry markers were within the normal range, except that renin level was at borderline in one of the affected subject (VI: 5, 2.53 mL/h). Detailed clinical and biochemical features of normal and affected members are mentioned in Table1. The height and weight were within normal range for the local population. An audiometric evaluation of selected affected members by measuring the threshold of hearing at 250–8000 Hz for pure-tune air conduction and bone conduction showed severe hearing loss across all frequencies.Table 1
Clinical and biochemical features of normal and affected members of studied family.
FeaturesReference rangesNormal membersAffected members from the family with hearing lossV : 3V : 5V : 1VI : 3VI : 5Age (y)3035451210Other symptoms of BSIVnonoNoNoNoSd·Na (mmol/L)136–148141140135137134S·K (mmol/L)3.6–5.03.94.63.43.53.5S·Cl (mmol/L)104–11410510797102101S·Mg (mg/dl)1.9–2.52.32.22.01.92.1S·HCO3 (mmol/L)17.5–27.530.027.530.030.231.0S·Ca (mg/dl)8.6–10.58.79.59.29.810.0S·Creatinine (mg/dL)0.85–1.350.70.80.60.50.4Pe Renin (ng/ml/hr)0.15–2.330.310.62.132.032.53S·Aldosterone (ng/dL)1–6;f 4–31>1.7f13.010.311.96.3S·Osmolality (mosm/kg)275–295293289281290287Ug Na (mmol/L)30–15017.5154.639.116.913.3U·K (mmol/L)20–673.491.3112.138.434.12U·Ca (mg/dL)6.016.015.311.313.0U·Mg (mg/dL)3.614.43.12.41.12U·Cl (mmol/L)46–16816.3174.346.5619.210.5Uosmolality (mosm/kg)50–14007490317910556NephrocacinosisSAbsentAbsentAbsentAbsentAbsentRenal sonography ruled out the presence of nephrocalcinosis or metabolic alkalosis in these affected individuals. Nor have the affected individuals displayed problems with polydipsia, polyuria, nocturnal enuresis, and hypocalciuria.


## 3.2. Genotyping Results

Analysis of the results obtained from genome search performed on the ten individuals of the family (III-1, III-3, IV-1, V-1, V-2, V-4, V-5, VI-1, VI-3, and VI-5) (Figure1) identified an area of interest on chromosome 1. Two-point analysis generated LOD score of 2.027 at marker D1S193 (73.21 cm) and 2.05 at marker D1S3462 (247.23 cm) on chromosome 1 (Table 2). In order to test linkage to these two regions, additional markers located in the vicinity of D1S3462 and D1S193 were chosen from Marshfield map [17] and genotyped in all the ten family members. The maximum two-point LOD score was increased to 2.55 with marker D1S2700. Multipoint linkage analysis for the family derived a maximum LOD score of 3.42 at marker D1S1661 (Table 3). The three-unit multipoint support interval contained a 10.48 Mb region, which span from markers D1S3721 to D1S2690. Haplotypes analysis delimited the centromeric boundary defined by a recombination between markers D1S1596 and D1S2770 in individual V-1. Recombination event between D1S2706 and D1S3721, defined the telomeric boundary in the same individual (V-1). Therefore, the minimum critical region of 18.36 Mb identified for disease locus and shared by all the affected individuals between the markers D1S2706 and D1S1596. The critical interval (17.27 cm) identified in this Pakistani family overlaps with the critical region to which DFNA2 was mapped on 1p34 in a large Indonesian family with autosomal dominant, progressive and sensorineural hearing loss [18]. The DFNA2 locus maps between markers D1S255 (65.47 cm) and D1S211 (73.81 cm), and thus shares a region of 2.68 cm with DFNB interval identified in our family. Riazuddin et al. [19] identified 1.5 Mb in four kindred segregating nonsyndromic deafness at the chromosome 1.Table 2
Two-point LOD score results between the locus identified in family and chromosome 1p34.2–p32.1 markers.
MarkerdeCode map postion1Marshfield map postion2Physical map postion3LOD score ATθ=0.000.010.050.10.20.3D1S47257.9565.4737,102,828−999.99−0.68−0.050.140.200.13D1S18660.1467.2237,900,316−999.99−.069−0.100.050.100.06D1S432—69.8639,046,164−1.334−0.69−0.18−0.020.030.01D1S159863.5170.4140,011,183−1.544−0.61−0.090.030.040.00D1S270665.3171.1340,741,361−999.99−1.34−0.180150.240.14D1S372165.8772.5941,395,5602.0271.971.751.491.490.56D1S19366.6473.2142,688,3692.0271.971.751.490.990.56D1S271367.8773.8144,182,5691.7261.671.471.240.800.43D1S213471.2975.6647,993,2411.6901.661.521.320.910.52D1S166172.4278.2551,142,5511.2281.191.060.900.600.33D1S265277.1580.7755,178,8521.3191.281.130.960.620.34D1S269079.6383.0756,867,5602.0271.971.751.490.990.56D1S270083.0787.3158,828,9552.5562.482.211.881.230.66D1S277085.0288.4059,442,255−999.99−0.260.260.350.290.17D1S1596—89.4959,670,589−1.55−0.52−0.000.100.080.01D1S273785.4790.5861,038,037−999.99−1.100.070.400.460.29D1S278887.2593.8661,568,857−999.99−2.16−0.90−0.46−0.14−0.03D1S23088.1295.3162,314,306−999.99−2.16−0.90−0.46−0.14−0.03D1S164899.62101.4873,125,563−999.99−2.27−1.00−0.55−0.20−0.06D1S166599.62102.0273,941,402−999.99−2.34−1.04−0.57−0.20−0.061Sex-average Kosambi cM map distances from the deCode genetic map [7].2Sex-average Kosambi cM map distances from the Marshfield genetic map [17].3Sequence-based physical map distance in bases according to the Human Genome Project-Santa Cruz (http://genome.ucsc.edu/, Assembly: May 2004).Table 3
Multipoint LOD score results between linked locus identified in family and chromosome 1p34.2–p32.1 markers.
MarkerPhysical positionLocationLODATA79C10—0.0000−2.0630D1S47237,102,8283.0940nanD1S18637,900,3165.1350nanD1S43239,046,1648.2290−2.4708D1S159840,011,1839.2390−0.0755D1S270640,741,36110.2490−8.2103D1S372141,395,56011.25903.0063D1S19342,688,36912.26903.2113D1S271344,182,56914.31003.3672D1S213447,993,24116.35103.4190D1S166151,142,55118.39203.4201D1S265255,178,85221.48603.3868D1S269056,867,56024.58003.2867D1S270058,828,95528.74902.8398D1S277059,442,25529.75901.3404D1S159659,670,58930.7690nanD1S273761,038,03731.7800−infD1S278861,568,85734.8730nanD1S23062,314,30636.9140−infD1S164873,125,56343.3060−infD1S166573,941,40244.3160nanGATA152F05L85,238,77045.3260nanFigure 1
Pedigree of family with nonsyndromic autosomal recessive hearing loss associated haplotypes. Haplotypes for the most closely linked STRPs are shown below each symbol. The genetic map distances according to the Marshfield genetic map in centimorgans (cM) are shown in parenthesis next to the marker name. Arrows adjacent to the haplotypes indicate key recombination events. The alleles are denoted 1–3 according to their allele size.


## 3.3. Mutation Analysis of Candidate Genes

Through a database search, we identified several genes mapping between the linkage interval in the family (Human Genome Project-Santa Cruz;http://genome.ucsc.edu/, May 2004).  Among these, KCNQ4 (MIM 603537), CLDN19, FOXE3 (MIM 601094), FOXD2 (MIM 602211), TSPAN1, and BSND (MIM 606412) were plausible candidate candidates in the interval of our family. We started sequencing of these and completely sequenced two genes (KCNQ4 and CLDIN19), but we found no mutation in both. Recent studies of Riazuddin et al. [19] showed that mutation in BSND can cause nonsyndromic deafness and is the molecular basis of DFNB73 locus in four Pakistani families segregating nonsyndromic deafness [19]. As the critical interval identified in the family contains potential candidate gene barttin (BSND, MIM 606412) as well; therefore, it was sequenced to search for the mutation. Sequence analysis of exon 1 of the BSND gene in our family also revealed missense mutation involving T to C transition at nucleotide number 35 (35T > C), which resulted in substitution of isoleucine to threonine at amino acid position 12 (I12T) (Figure 2). BSND mutation; that is, p.I12T cosegregated in all extant members of the pedigree. This mutation was present in the heterozygous state in obligate carriers within the family.Representative chromatogram generated by Big Dye Terminator, sequencing of translated exon1 of BSND gene indicating a T to C transition at nucleotide 35 (35T > C) from an affected female individual (VI-1) (a), a heterozygous carrier (VI-1) (b) and a normal male individual (V-2) (c) of the family. The T to C transition at nucleotide position 35 results in isoleucine to threonine change at amino acid 12 (I12T). Arrows indicate the nucleotide change in the sequence.
(a)(b)(c)


## 3.4. Bioinformatics Analysis

Models of BSND CLCNKB and REN proteins were generated by SAM_T08 server and visualized by Rasmol (Figure3(a)–3(c)). Rampage values for BSND were number of residues in favored region 95.6%, the number of residues in allowed region 3.1%, and the number of residues in outlier region 1.4% (Figure 4(a)); however, the CLCNKB and REN values were little varying. Its values were the number of residues in favored region:  94.8%, 97%, the number of residues in allowed region:  4.2%, 3%, and the number of residues in outlier region is 0.9% for CLCNKB and 0% for REN (Figure 4(c)).Protein models of BSND (a) CLCNKB (b), and REN (c) generated by SAM_T08 and visualized by Rasmol (display = ribbon; colors = group).
(a)(b)(c)Rampage results of residues for protein model of human BSND protein (a) CLCNKB (b) and REN (c).
(a)(b)(c)Comparison of conserved amino acid of human BSND with closely related species like mouse, rat, and rabbit shows that isoleucine at position 12 is highly conserved, so mutation at this point can be  significant (Figure5(a)). The evolutionary relationship between human, mouse, rat, and rabbit were evaluated by phylogenetic analysis. The analysis represents that the BSND of mouse and rat is closely related to each other and shows homology with human and rabbit (Figure 5(b)).(a) Comparison of conserved regions of BSND protein between human, mouse, rat, and rabbit by using UniportKB ClustalW. Dark region shows conserved amino acid sequence between species. (b) Phylogenetic analysis of BSND protein. It shows evolutionary relationship between human, mouse, rat, and rabbit species.
(a)(b)Our mutation (I12T) of BSND was generated by Swisspdb server and was analyzed by Viewerlite. This substitution of isoleucine with threonine altered the bonding capacity with other side chains as well as other interactors (Figure6(a)-(b)).Figure 6
(a) Normal and mutated model (Il12T) of BSND are superimposed and visulized in Viewerlite. (b) Focused normal model is yellow highlighted and displayed in stick format, and mutated is displayed in ball and stick format.According to STRING database, the physical and functional interaction of BSND was very much significant with CLCNKB protein. Clinical assessment showed the reduced chlorine level in all the affected individuals. According to further clinical assessments, renin level was at border range in one affected individual, so REN protein was also selected to check the effect of the respected mutation on this protein. However, other interactors are COL7A1, ARL2, CLCN5, KCNJ1, CASR, SLC12A1, SLC12A3, and CLCNKA (Figure7). Both the proteins were docked with HEX v 6.1, taking BSND as receptor and CLCNKB and REN as a ligand. The total energy value calculated for BSND-CLCNKB molecules was −229094.6   and BSND-REN was −1145.7. However, the total energy value for the mutant (I12T) on docking with CLCNKB was maximized to −218782.3   and with REN was −1149.8. This values fluctuation ultimately results in reduced binding affinity with CLCNKB and increased binding affinity with REN and so may be the case with other interacting proteins.Figure 7
The interactions of BSND with other functionally similar proteins retrieved by STRING online database. CLCNKB  showed  highest confidence score, that is, 0.990 with BSND.


## 4. Discussion

Pathophysiological pathway leading from a specific mutation to a specific phenotype has remained elusive in syndromic as well as nonsyndromic hearing loss families. Individuals with the same mutation can fall along a clinical spectrum ranging from asymptomatic to severely affected and can even have completely different diseases with different mutation in the same gene. Our results also support the findings of Riazuddin et al. [19] that pathogenic mutation of BSND gene can cause the nonsyndromic hearing loss. But there are no subclinical renal metabolic changes in our family as reported by Riazuddin et al. [19]. Renin level is at boarder range only in one affected individual, while renal sonography ruled out the presence of nephrocalcinosis in these affected individuals. Nor have the affected individuals displayed problems with polydipsia, polyuria, nocturnal enuresis, metabolic alkalosis, and hypocalciuria. It might be possible that partial loss of function of barttin induced by T12I allele cause only selectively hearing loss and appear to have less pronounced effects. However, the genotype does not always predict the clinical phenotype, which varies both within and between families carrying the mutation in the same gene, implying the existence of other genetic and/or environmental factors that influence phenotype. Several disease-causing BSND mutations have been identified and functionally analyzed. In all cases, there is a genotype-phenotype relationship in that the level of function of mutant barttin predicts the renal phenotype. Our results support the relationship between missense BSND mutation (I12T) and nonsyndromic hearing loss. It is a secondary report that the BSND as causal gene in nonsyndromic deafness.The protein sequence of BSND is highly conserved in closely related species of human, mouse, rat, and rabbit, and it contains 2 putative transmembrane domain starting from 7 to 26 amino acid and 31 to 53 amino acid [20].Bioinformatics analysis of mutation in BSND protein suggest that isoleucine, a hydrophobic amino acid, is being converted into threonine, a polar amino acid, which is more susceptible to posttranslation modifications, affecting the quaternary structure of protein resulting in mutant protein.Barttin (BSND) is an accessory subunit that modifies protein stability, subcellular distribution, and voltage-dependent gating of ClC-K chloride channels expressed in renal and inner ear epithelia. ClC-K channels are double-barreled channels with two identical protopores that may be opened by individual or common gating processes [21].Potassium (K+) secretion by strial marginal cells and vestibular dark cells require Cl− to recycle in the basolateral membrane via a major Cl− conductance. This Cl− conductance is composed of the Cl− channels CLCNKA/BSND and CLCNKB/BSND. Cl− channels CLCNKA/BSND and CLCNKB/BSND consist of the pore-forming CLCNKA α-subunits and CLCNKB and the BSND β-subunit. Mutations of BSND reduce channel conductivity and surface expression [22].Estévez et al. [2] proved that BSND protein is an essential β subunit for CLCNKA and CLCNKB chloride channels, with which it colocalizes in basolateral membranes of renal tubules and of potassium-secreting epithelia of the inner ear. Disease-causing mutations in either CLCNKB or BSND compromise currents through heteromeric channels [2].The CLCNKA and CLCNKB channels are members of the ClC family, which comprises at least 9 mammalian chloride channels. Each is believed to have 12 transmembrane domains and intracellular N and C termini. The prototype of the family in torpedo is gated by both voltage and chloride [23].Docking results of both normal and mutant interaction with CLCNKB suggested that there is a reduced binding affinity on mutation which may results the reduced down regulating or up regulating interactions; however, the mutant does not effect the interaction with renin (REN). So, we can conclude that the boarder range of renin level in one affected individual is not dependent upon this mutation.Identification of mutation in BSND gene in more families will give valuable insight into the genetic mechanisms underlying nonsyndromic hearing disorder. This work is a fundamental step in the bioinformatics analysis of BSND and provides basics for further analysis.


## Electronic-Database Information

(i)
Marshfield Medical Research Foundation database website:http://research.marshfieldclinic.org/genetics/.(ii)
National Center for Biotechnology Information:http://www.ncbi.nlm.nih.gov/.(iii)
UCSC Genome Browser:http://genome.ucsc.edu/.



---

*Source: 304612-2011-04-28.xml*