CC BY
55
Haplotype Diversity and Reconstruction of Ancestral Haplotype Associated with the c.35delG Mutation in the GJB2 (Cx26) Gene among the Volgo-Ural Populations of Russia
L. U. Dzhemileva1*, O. L. Posukh2, N. A. Barashkov3, S. A. Fedorova3, F. M. Teryutin3,
V. L. Akhmetova1, I. M. Khidiyatova1, R. I. Khusainova1, S. L. Lobov1, E. K. Khusnutdinova1
institute of Biochemistry and Genetics, Ufa Research Center, Russian Academy of Sciences
institute of Cytology and Genetics, Siberian Branch, Russian Academy of Sciences
3Yakut Research Center of Complex Medical Problems, Siberian Branch, Russian Academy of
Medical Sciences
*E-mail: dzhemilev@anrb.ru
Received 23.03.2011
Copyright © 2011 Park-media, Ltd. 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.
ABSTRACT The mutations in the GJB2 (Cx26) gene make the biggest contribution to hereditary hearing loss. The spectrum and prevalence of the GJB2 gene mutations are specific to populations of different ethnic origins. For several GJB2 mutations, their origin from appropriate ancestral founder chromosome was shown, approximate estimations of “age” obtained, and presumable regions of their origin outlined. This work presents the results of the carrier frequencies’ analysis of the major (for European countries) mutation c.35delG (GJB2 gene) among 2,308 healthy individuals from 18 Eurasian populations of different ethnic origins: Bashkirs, Tatars, Chuvashs, Udmurts, Komi-Permyaks, Mordvins, and Russians (the Volga-Ural region of Russia); Byelorussians, Ukrainians (Eastern Europe); Abkhazians, Avars, Cherkessians, and Ingushes (Caucasus); Kazakhs, Uzbeks, Uighurs (Central Asia); and Yakuts, and Altaians (Siberia). The prevalence of the c.35delG mutation in the studied ethnic groups may act as additional evidence for a prospective role of the founder effect in the origin and distribution of this mutation in various populations worldwide. The haplotype analysis of chromosomes with the c.35delG mutation in patients with nonsyndromic sensorineural hearing loss (N=112) and in population samples (N =358) permitted the reconstruction of an ancestral haplotype with this mutation, established the common origin of the majority of the studied mutant chromosomes, and provided the estimated time of the c.35delG mutation carriers expansion (11,800 years) on the territory of the Volga-Ural region.
KEYWORDS hereditary nonsyndromic sensorineural hearing loss; GJB2 (Cx26) gene; c.35delG mutation; ancestral haplotype; populations of the Volga-Ural region.
INTRODUCTION
Hereditary deafness is a frequent disorder in humans: it is recorded in 1/1,000 newborns. The etiology and pathogenesis of this disease are still to be clarified; however, approximately half of the cases of hereditary deafness are a result of genetic disorders [1].
The hereditary forms of innate hearing loss are characterized by clinical polymorphism and genetic heterogeneity. In the nuclear genome, about 114 loci were mapped and 55 genes were identified, whose mutations to a certain extent cause hearing loss. About 80% of all cases of hereditary nonsyndromic hearing loss fall within the category of autosomal-recessive forms;
15-20% - autosomal-dominant, and about 1% - the
form linked with the X-chromosome and mitochondrial forms of deafness [2].
The most frequent cause of nonsyndromic autosomal recessive hearing loss in humans is the mutations in the GJB2 gene (gap junction P2, subunit P2 of the gap junction protein), localized in the chromosomal region 13q11-q13 and coding connexin 26 (Cx26), which is the transmembrane protein involved in the formation of connexons. Connexons are structures consisting of six protein subunits, which form cellular channels, ensuring a fullblown ion exchange among adjacent cells. This facilitates the maintenance of the homeostasis of endolymph in cochlea tissues. Recently, the fine structure of intercellular channels formed by connexin 26
Fig. 1. Structure of intercellular channels formed by molecules of connexin 26. A, B, C, D, F, E and A', B', C', D', F ', E' - connexin 26 molecules in connexones of neighboring cells; TM 1-4 -transmembrane protein segments of Cx26; NTH - N-terminal helix of protein Cx26. Figure was adapted from [3] by permission from Macmillan Publishers Ltd.
Cytoplasmic diameter 92 Â
Intracellular area 19 Â
Transmembrane area 38 Â
Extracellular area 42 Â
Transmembrane area 38 Â
— 4
Maximum diameter 92 Â
Minimum diameter 51 Â
was reported (Fig. 1) [3]. When there are defects in connexin 26, the functioning of the intercellular channels is irreversibly disrupted in the tissues of the internal ear and endolymph homeostasis is not restored: a factor that is necessary for normal sound perception
[4].
So far, in the GJB2 gene, over 150 pathogenic mutations (mainly recessive), several polymorphisms and sequencing variants whose role in the pathogenesis of hearing loss is still unclear, have been described [2]. The spectrum and frequencies of the GJB2 gene mutations are characterized by significant interpopulation differences. The racial and/or ethnic specificity of the spread of several GJB2 gene mutations is preconditioned in some cases by the founder effect, as well as by, in all likelihood, the geographic and social isolation of some populations. Researchers have managed to show the origin of some GJB2 gene mutations from an ancestor founder chromosome, to obtain approximate estimations of the “age” of the mutations, and to outline suggested regions where they appeared [5-10]. The mutation c.35delG (p.Gly12Valfsx1) is the most frequent in Europe. It first emerged, according to various estimates, 10,000 to 14,000 years ago on the territory of the Middle East or Mediterranean region (probably on the territory of modern Greece) [10-12] and spread in Europe with the migrations of the neolithic population of Homo sapiens [6]. Haplotype analysis of the chromosomes with mutation c.235delC (p.Leu79Cysfsx3) in
populations of Japan, Korea, China, and Mongolia has allowed to put forward a hypothesis about the founder effect regarding the origin and distribution of this mutation in East Asia and to estimate its “age” (~11,500 years) and presumed appearance in the Lake Baikal region, from where it spread, by way of sequential migrations, over Asia [7]. The “age” of the mutation c.71G>A (p.Trp24X), most widespread in India, has also been assessed (~7,880 years) [8]. The ethnic specificity of the mutation c.167delT (p.Leu56Argfsx26) and mutation c.427C > T (p.Arg143Trp) was shown for populations of the Ashkenazi Jews [5] and for some populations of Western Africa (Ghana) [13, 14], respectively.
In Russia, several research groups have studied the hereditary forms of deafness [15-32]. In most cases, they have considered the genetic-epidemiological and clinical-genetic peculiarities of the inherited forms of hearing loss, and a series of papers are devoted to the molecular-genetic analysis of the GJB2 gene or its single mutations [17-21, 25-30]. Some authors have obtained data on the specificity of the range and frequency of separate mutations in the GJB2 gene function in the studied region. For instance, the most frequent mutations in the Siberian populations (Yakuts and Altaians) are IVS1 + 1G>A [27] and c.235delC [25], respectively. The populations of the Volga-Ural region, as generally in the European part of the continent, predominantly exhibit the mutation c.35delG [18-22, 28]. Local differences in the carrier frequency of the
mutation c.35delG are probably related to the genetic history of some populations, and factors of population dynamics, and migration routes of the c.35delG carriers in the world. The available data on the contribution of the GJB2 mutations to the development of a pathology in patients with nonsyndromic sensorineural hearing loss (NSHL), living in the Volga-Ural regions, and the population-based data on the carrier frequency of the most significant recessive mutation c.35delG permitted an adequate assessment of the haplotypic diversity of the chromosomes bearing the mutation c.35delG, the reconstruction of the possible ancestor haplotype linked to that mutations, and the estimated time of its appearance on the territory of the Volga-ural regions, which represents the eastern border of the habitat of the mutation c.35delG.
EXPERIMENT
The material taken for the haplotypic analysis and estimates of the “age” of c.35delG of the gene GJB2 was 56 DNA samples (112 chromosomes) obtained from patients with NSHL, residing in the Volga-ural regions, in which the mutation.35delG was identified in the homozygote state (32 Russians, 10 Tatars, 1 Bashkirs, 4 Ukrainians, 2 Armenians, and 7 individuals of mixed ethnicity). The control group included 179 (358 chromosomes) healthy individuals from three ethno-geograph-ical groups of the Russians (N = 86); Tatars (N = 62); and Bashkirs (N = 31) without this mutation.
To analyze the carrier frequency of c.35delG 2,308 DNA samples were used, which were obtained from healthy individuals-representatives of different populations of the Volga-Ural region, Central Asia, Northern Caucasus, Eastern Europe, and Siberia, belonging to four language families (Table 1).
The blood samples were obtained during expeditions between the years of 2000 and 2010, after receiving the informed written consent of the participants in
the study. The genomic DNA was extracted from peripheral blood using the phenol-chloroform extraction method.
The present scientific-research work was approved by the local committee for biomedical ethics at the IBG UNT of the Russian Academy of Sciences (Ufa).
Screening of the a35delG mutation in the GJB2 gene The screening of the mutation c.35delG in the GJB2 gene was carried out with the allele-specific amplification of fragments of the coding region of the GJB2 gene using the primers listed in Table 2. The results were visualized by vertical electrophoresis in a 10% polyacrylamide gel (PAAG) with further staining of the ethidium bromide solution in the standard concentration and viewing in ultraviolet rays.
Analysis of haplotypes and estimate of the “age” of the mutation c.35delG
For the analysis of the haplotypes and the estimated age of the mutation c.35delG in the GJB2 gene, three high-polymorphic microsatellite CA-markers were used: D13S175, D13S141, and D13S143 [6, 9, 10, 12, 36], flanking the DFNB1 locus, which contains the GJB2 gene. The physical and genetic localization of the markers at chromosome 13 and genetic distances between them, as well as the GJB2 gene, were identified on the basis of the Marshfield genetic linkage map (http:// www.ncbi.nlm.nih.gov/mapview/). The total physical distance of the flanked region was ~ 2 Mb. (Fig. 2). The choice of the markers was preconditioned to fall with the strive to get the possibility of comparable data, as earlier these markers were used for the assessment of the age of mutations in the gene GJB2 in different populations [6, 9, 10, 12, 36].
The STR-markers were genotyped using PCR at the thermocycler (Eppendorf), and appropriate oligonucleotide primers were used (Table 2). Products of the
1.5 Mb.
38.5 Kb.
84 Kb
GJB2
(19659605-19665114 b.p.)
D13S141
(19622459-19622584 b.p.) 1.535 cM (Marshfield)
D13S175
(19746380-19746734 b.p.) 5.94 cM (Marshfield)
D13S143
(21172085-21172213 b.p.) 6.03 cM (Marshfield)
Fig. 2. Localization of the microsatellite markers D13S141, D13S175, and D13S143, flanking the GJB2 gene, on chromosome
13. The distance between the GJB2 gene and the markers is indicated by arrows.
Table 1. Carrier frequencies of the c.35delG mutation in 18 ethnic groups dwelling on the territory of Eurasia
Population Linguistic affiliation (group) Region N Number of heterozygote carriers of the mutation c.35delG / total numbers of individuals tested (carrier frequency)
Eastern Europe
Byelorussians Indo-European / Slavic The Republic of Belarus (disperse sample) 97 6/97 (0.062)
Ukrainians Indo-European / Slavic Kharkov region and Poltava region, Ukraine 90 © 3 ©^ 3 ÇO )
Volga-Ural region
Russians Indo-European / Slavic Yekaterinburg, the Russian Federation (RF) 92 2/92 (0.022)
Bashkirs Altaic / Turkic Baimakskii, Burzyanskii, Abzelilovskii, Kugarchinskii, Salavatskii, and Arkghangelskii districts, the Republic of Bashkortostan, RF 400 1/400 (0.003)
Tatars Altaic / Turkic Al’met’evskii and Yelabuzhskii districts, the Republic of Tatarstan, RF 96 1/96 (0.010)
Chuvashs Altaic / Turkic Morgaushskii district, the Chuvash Republic, RF 100 0/100
Mordvins Uralic / Finno-Ugric Staroshaiginskii district, the Republic of Mordovia, RF 80 5/80 (0.062)
Udmurts Uralic / Finno-ugric Malopurginskii district of the Udmurt Republic and Tatyshlinskii district of the Republic of Bashkortostan, RF 80 3/80 (0.037)
Komi- Permyaks Uralic / Finno-Ugric Kachaevskii district, the Komy-Permyak Autonomous District, RF 80 0/80
Central Asia
Kazakhs Altaic / Turkic Alma-Atinskaya region, Kyzylordinskaya region, and Abaiskii district, Kazakhstan 240 2/240 (0.008)
Uighurs Altaic / Turkic Alma-Atinskaya region, Kazakhstan 116 1/116 (0.009)
Uzbeks Altaic / Turkic The Republic of Uzbekistan (disperse sample) 60 0/60
Caucasus
Abkhazians North-Caucasian / Adygo-Abkhaz Abkhazia and Georgia (disperse sample) 80 3/80 (0.038)
Avars North-Caucasian / Dagestan Gumbetovskii district, the Republic of Dagestan, RF 60 0/60
Cherkessians North-Caucasian / Adygo-Abkhaz The Karachaevo-Cherkess Republic, RF 80 1/80 (0.013)
Ingushes North-Caucasian / Nakh Nazran’ district, the Republic of Ingushetia, RF 80 0/80
Siberia
Altaians Altaic / Turkic The Republic of Altai, RF 230 0/230
Yakuts Altaic / Turkic Megino-Kangalasskii, Amginskii, Churapchinskii, Tattinskii, Verkhnevilyuiskii, Vilyuiskii, Nyurbinskii, and Suntarskii uluses (districts), the Republic of Sakha (Yakutia), RF 247 1/247 (0.004)
Table 2. Sequences of primers used for amplification
Locus Name and nucleotide sequence of primers Detection method Reference
GJB2 (13q11-q12) 35delG F 5’-CTTTTCCAGAGCAAACCGCCC-3’ 35delG R 5’-TGCTGGTGGAGTGTTTGTTCAC-3’ Visualization of PCR fragments in 10% PAAG [15]
D13S141 F- 5’-GTCCTCCCGGCCTAGTCTTA-3’ R-5’-ACCACGGAGCAAAGAACAGA-3’ [33]
D13S143 F-5’-CTC ATG GGC AGT AAC AAC AAAA-3’ R-5’-CTT ATT TCT CTA GGG GCC AGC T-3’ [34]
D13S175 F-5’-TAT TGG ATA CTT GAA TCT GCT G-3’ R-5’-TGC ATC ACC TCA CAT AGG TTA-3’ [35]
PCR-reaction were separated by vertical electrophoresis (glass size 20 x 20 cm, Helicon, Russia) in 10% PAAG and 5% glycerin. The gels were stained with silver ions.
The linkage disequilibrium between alleles of the 13th-chromosome loci was calculated using the following formula:
6 = (Pd- Pn)/(1 - Pn),
where 6 is the measure of linkage disequilibrium, Pd is the frequency of the associated allele among mutant chromosomes, and Pn is the frequency of the same allele among intact chromosomes [37].
The statistical significance of differences in the frequencies of the alleles of the studied markers on 112 chromosomes containing c.35delG and 358 chromosomes without this mutation was assessed using the standard x2 test 2x2 (the MedStat software).
The age of expansion of the founder haplotype bearing the mutation c.35delG in the GJB2 gene was estimated via the genetic clock approach [38], which is based on the definition of the number of generations (q) from the moment of the mutation appearance in the population to the present, proceeding from the ratio of linkage disequilibrium in terms of polymorphous markers linkage with the locus of the disorder. This age was calculated using the following formula:
q = log[1 - Q/(1 - Pn)]/log(1 - 0),
where q is the number of generations since the moment the mutation appeared in the population, Q is the share of mutant chromosomes without the founder haplo-type, Pn is the frequency of the allele of the founder haplotype in the population, and 0 is the recombinant fraction. The value of 0 was computed given the physical distance of the markers from the location of the mutation, stemming from the ratio 1 cM = 1,000,000 b.p.
The value for the allele association was estimated by the coefficient of standard linkage disequilibrium according to [39]:
p - q P + q)(2 - p - q)
where p and q are the frequencies of the alleles or hap-lotypes of the normal (p) and mutant (q) chromosomes.
The haplotypic diversity indicator equivalent to the expected heterozygosity was calculated using the formula:
h = (1 -^ *?) N/(N -1) ,
where x is the frequency of each haplotype in the population, and N is the sample size [40].
RESULTS AND DISCUSSION
In many populations worldwide, the main contributor to the development of nonsyndromic sensorineural hearing loss (NSHL) is the mutations of the GJB2 gene. In most European populations, up to 40-50% of cases of NSHL are preconditionally caused by one of the major recessive mutations of this gene, c.35delG, which is revealed in the homozygous or compound-heterozygous state [41]. In connection with this, many researchers have analyzed the carrier frequency of c.35delG in a variety of populations in the world. A large-scale study which embraced 17 European countries showed that the mean carrier frequency of the c.35delG mutation amounts to 1.96% (1/51), ranging from 2.86% (1/35) in South-European countries to 1.26% (1/79) in Northern Europe [42]. In the Mediterranean region, the highest carrier frequency of c.35delG were observed in Greece (3.5%), in the southern regions of Italy (4.0%), and in France (3.4%) [43]. As a result of the meta-analysis of the c.35delG carrier frequency in over 23,000 individuals from different populations, carried out on the basis of the data published between the years of 1998 and 2008, the average regional frequencies of c.35delG were determined in the European (1.89%), American (1.52%), Asian (0.64%), and African (0.64%) populations, as well as in Oceania (1%). Also, the decreasing gradient of the carrier frequency of c.35delG (from 2.48 to 1.53%) from
south to north in European populations and from west to east (from 1.48 to 0.1%) in Asian populations [44] was confirmed.
Carrier frequency of the mutation c.35delG
An analysis of the carrier frequency of the mutation c.35delG in 18 populations of Russia and the former Soviet states was carried out (Table 1).
High carrier frequencies of the c.35delG were revealed in two Eastern European populations: Ukrainians (3.3%) and Byelorussians (6.2%). In the Turkicspeaking populations of the Volga-Ural region, the carrier frequencies of the mutation c.35delG were 1.0, 0.3, and 0% in Tatars, Bashkirs and Chuvashs, respectively. In the Finno-Ugric populations of the Volga-Ural region, the c.35delG mutation was present with a high carrier frequency of 6.2% in Mordvins, 3.7% in Udmurts, and it was absent in Komi-Permyaks. These frequency fluctuations among the studied populations of the Vol-ga-Ural region are likely due to the specific features of the historic development of these populations in the region, or could be the consequence of the relatively small size of the samples. Earlier, a high carrier frequency of c.35delG (4.4%) was found in Estonians, an apparent exception for Northern European populations, who typically have low frequencies of c.35delG [42]. These data, as well as the data obtained during other studies [15, 20-22, 28, 29], indicate the significant variation in the carrier frequency of c.35delG among indigenous populations of the Volga-Ural region. The carrier frequency of the c.35delG, singled out in Russians (2.2%), was comparable to the results for the Russian population in the Central region of Russia [15, 18, 28]. In the Turkic-speaking populations of central Asia (Kazakhs, Uighurs, and Uzbeks) the mutation c.35delG with a low
•20 0 20 40 60 80 100 120 140 160 180
Fig. 3.The spatial distribution of the carrier frequency of the c.35delG mutation in the GJB2 gene in the populations of Eurasia (performed using the program SURFER 9.0 Golden Software Ink).
carrier frequency was observed in Kazakhs (0.8%), Ui-ghurs (0.9%), and it was absent in Uzbeks. In the Turkic-speaking populations of Siberia (Yakuts and Altaians) a relatively low carrier frequency of the mutation c.35delG (0.4%) was revealed in the Yakut population but was not detected in Altaians. The North Caucasus region was in the past one of the crucial migration corridors in Eurasia. It is characterized by a high diversity of the population and a complex historical development of its resident ethnoses. In the North caucasus populations (Abkhazians, Avars, cherkessians, and Ingushes), the mutation c.35delG was discovered only in Abkhazians (3.8%) and Cherkessians (1.3%).
The spatial distribution of the carrier frequency of the mutation in the Eurasian populations obtained on the basis of our own data and that of relevant published information on the c.35delG available as of 2010 [24] is presented in Fig. 3.
The data obtained significantly added to the picture of the mutation c.35delG distribution in Eurasia: European regions of Russia, just as the European part of the continent on the whole, are characterized by a high frequency of c.35delG, and this mutation is widespread in the polyethnic population of the Volga-Ural region. However, the mechanisms of its spreading and time of appearance in the Volga-Ural region are yet to be studied. To answer these questions, we carried out a hap-lotypic analysis of the chromosomes carrying c.35delG and those without it, using three high-polymorphous microsatellite CA-markers: D13S175, D13S141, and D13S143 [6, 9, 10, 12, 36], which flank the locus DFNB1 containing the GJB2 gene (Fig. 2).
Frequencies of alleles of loci D13S141,
D13S175, and D13S143
Table 3 lists the distribution of frequencies for the alleles of microsatellite markers D13S141, D13S175, and D13S143 in the NSHL patients (c.35delG-mutant chromosomes) and in the control sample (normal chromosomes), including three ethnic groups (Russians, Tatars, and Bashkirs).
In other studies, the area sized 2 Mb and covered by these three markers allowed for the calculation of the approximate number of past generations since the start of the expansion of the proposed founder hap-lotype, including mutations of the GJB2 gene, in the populations of India (mutation p.Trp24X) [8] and Morocco (c.35delG) [9]. A panel consisting of 8 STR-mark-ers (D4S189, D13S1316, D13S141, D13S175, D13S1853, D13S143, D13S1275, and D13S292) and 2 SNP-markers was applied in dating the splicing site IVS1 + 1G>A mutation of the GJB2 gene in the Yakut population [27]. The age of the mutations c.35delG and c.235delC was computed using the 6 SNP-markers [6, 7]. In later stud-
Table 3. Distribution of allele frequencies of microsatellite markers D13S141, D13S175, and D13S143 in patients with NSHL (chromosomes with the c.35delG mutation) and in the control sample (normal chromosomes)
Allele (b.p.) chromosomes with c.35delG mutation (N=112) normal chromosomes (N=358) X2 Р (95% significance level)
number of chromosomes Allele frequency number of chromosomes Allele frequency
D13S141
113 0 0 10 0.027±0.008 0.3131 0.6
123 26 0.232±0.031 211 0.589±0.026 43.458 0.000
125 84 0.750±0.042 112 0.312±0.024 67.058 0.000
127 2 0.017±0.002 25 0.069±0.013 4.2629 0.045
D13S175
101 3 0.026±0.012 21 0.058±0.01 1.793 0.300
103 8 0.071±0.021 88 0.245±0.02 15.875 0.000
105 91 0.812±0.036 157 0.438±0.02 47.866 0.000
107 1 0.008±0.007 30 0.083±0.01 7.763 0.005
109 6 0.053±0.024 38 0.106±0.01 2.783 0.100
111 0 0 7 0.019±0.007 2.222 0.150
113 3 0.026±0.01 17 0.047±0.01 0.897 0.064
D13S143
126 1 0.008±0.007 0 0 3.192 0.04
128 1 0.008±0.007 5 0.013±0.006 0.169 0.65
130 90 0.80±0.048 283 0.79±0.021 0.088 0.81
132 6 0.05±0.026 26 0.07±0.013 0.490 0.54
134 12 0.11±0.016 37 0.10±0.013 0.012 0.89
136 2 0.017±0.021 3 0.008±0.004 0.727 0.55
138 0 0 4 0.011±0.005 3.278 0.05
ies aimed at clarifying the age of the mutation c.35delG in Greece, two STR-markers, D13S175 and D13S141, and six SNP-markers were used [12].
D13S141. The marker D13S141 has seven allelic variants [9, 12]; however, only four of them were revealed in the ethnic groups from the Volga-ural region. The allele 123 (D13S141) frequency is significantly higher (X2 = 43.458; р = 0.000) on chromosomes of individuals from the control group (59%), whereas in NSHL patients, it is only 23%. The allele 125 (D13S141) is observed in the mutant chromosomes with a frequency of 75%, which is significantly higher than in normal chromosomes (31%) (X2 = 67.058; р = 0.000), and it matches the data on the allele 125 (D13S141) dominance on chromosomes with
the mutation a35delG in NSHL patients from Morocco, Greece, Palestine, and Israel [9, 10, 12, 45].
D13S175. The marker D13S175 has eight allelic variants [9, 10, 12], of which seven are present in the ethnic groups of the Volga-ural region. the allele 105 (D13S175) on c.35delG chromosomes is observed at a frequency of 81.2%, which is significantly higher in comparison with the normal chromosomes (43.8%) (X2 = 47.866; p = 0.000), and allele 103 (D13S175) was significantly more frequent in the controls (x2 = 47.866; p = 0.000 and X2 = 15.87; p = 0.000, respectively) (Table. 3). Earlier, it was shown that, on the chromosomes with the mutation a35delG in NSHL individuals from tunisia, Algeria, Morocco, and Greece, allele 105
(D13S175) was also very frequent (from 67 to 100%) [9, 12, 46]. Allele 111 (D13S175) was not observed on the chromosomes of individuals with NSHL but it was present in 2% of the chromosomes of people with no hearing problems.
D13S143. The marker D13S143 has eight allelic variants [9, 12], of which seven are present in the ethnic groups residing in the Volga-Ural region. Allele 130 (D13S143) is the most frequent both on normal chromosomes (79%) and c.35delG chromosomes (80%), and allele 134 (D13S143) is more often detected on chromosomes of c.35delG-mutant individuals (x2 = 9.909; p = 0.005).
An analysis of the distribution of the allele frequencies of three microsatellite loci D13S141, D13S175, and D13S143 on normal and c.35delG chromosomes revealed a pronounced misbalance in linkage between the specific alleles of these markers and mutation c.35delG in the GJB2 gene (Table. 3). The degree of association of the out-gene microsatellite loci under study vividly reflects the standard coefficient of the allele association (ASt) [39]. The greatest degree of linkage with the mutation c.35delG is typical of allele 125 of the marker D13S141 (ASt = -0.438) and allele 105 of the marker D13S175 (ASt = -0.386).
Haplotype analysis and age of c.35delG mutation
Given the data we obtained during the study of the polymorphism of the markers D13S141, D13S175, and D13S143, and the linkage disequilibrium of some of their alleles with the mutation c.35delG in the GJB2 coding region, we suggested that they may be evidence of the presence of a single ancestor haplotype, which carries this mutation. As a result, for three polymorphic
loci, haplotypes of members of each of the 56 families with hereditary deafness and healthy donors were constructed. A precise identification of the haplotype by the alleles D13S175-D13S143-D13S141 is possible for 112 mutant chromosomes with c.35delG and 358 normal chromosomes.
In all of the chromosomes analyzed, 59 different variants of haplotypes were revealed, of which 52 were found on normal chromosomes and 25 on c.35delG-mutant chromosomes (Table 4).
The distribution of the haplotypes on 358 normal chromosomes is characterized by a high value of hap-lotypic diversity (h = 0.943), the frequency of the most wide-spread haplotype 123-105-130 amounts to 17.8%, and 11 other haplotypes are found at frequencies exceeding 2%. For the distribution of haplotype frequencies on 112 chromosomes with the mutation c.35delG, a lower value of the haplotypic diversity (h = 0.645) is typical, the haplotype 125-105-130 is the most frequently found (59%), and the frequency of six haplo-types exceeds 2%. Seven haplotypes rarely occurring on the mutant chromosomes (below 2%) were not detected on normal chromosomes. The graphic mapping of the occurrences of the haplotypes D13S141-D13S175-D13S143 on the normal chromosomes of healthy donors and c.35delG-mutant chromosomes in NSHL patients is illustrated in Fig. 4.
An analysis of the distribution of the haplotypes D13S141-D13S175-D13S143 on normal chromosomes in people of different ethnic origins (Russians, Tatars, and Bashkirs) showed differences in the spectrum of haplotypes in the surveyed ethnic groups and statistically significant differences in the frequencies of the haplotypes (x2 = 57.335; p = 0.000; d.f. = 56). In Russians
Fig. 4. The distribution of frequencies of haplotype D13D141-D13S175-D13S143 on normal chromosomes and chromosomes with a mutation c.35delG in patients with NSHL (non-syndromic sensorineural hearing loss). Along the vertical and horizontal axes the frequency of the haplo-type and the names of the haplotype are indicated, respectively. The arrow shows the haplotype 125105-130.
Haplotype
125-105-130
Chiomosomes with c.35delG mutation ' '
Normal chromosomes '
Table 4. Frequencies of haplotypes D13S141-D13S175-D13S143, identified in chromosomes of patients with the o35delG / o35delG genotype and in normal chromosomes of healthy donors living in the Volga-Ural region
Haplotypes Patients with c.35delG\c.35delG Overall control Russians Tatars Bashkirs
Absolute value Frequency Absolute value Frequency Absolute value Frequency Absolute value Frequency Absolute value Frequency
113-101-130 1 0.009 1 0.002 0 0 1 0.008 0 0
123-101-130 1 0.009 12 0.033 4 0.023 7 0.056 1 0.016
125-101-130 0 0 4 0.011 4 0.023 0 0 0 0
113-101-134 1 0.009 0 0 0 0 0 0 0 0
125-101-136 0 0 1 0.002 1 0.006 0 0 0 0
123-103-128 0 0 1 0.002 0 0 1 0.008 0 0
113-103-130 0 0 3 0.008 2 0.012 1 0.008 0 0
123-103-130 2 0.018 28 0.078 13 0.076 12 0.097 3 0.048
125-103-130 1 0.009 20 0.055 12 0.070 4 0.032 4 0.064
127-103-130 0 0 5 0.013 4 0.023 0 0 1 0.016
113-103-132 0 0 1 0.002 0 0 1 0.008 0 0
123-103-132 0 0 1 0.002 0 0 1 0.008 0 0
125-103-132 0 0 2 0.005 1 0.006 1 0.008 0 0
123-103-134 3 0.027 8 0.022 2 0.012 4 0.032 2 0.032
125-103-134 0 0 3 0.008 3 0.017 0 0 0 0
125-103-136 1 0.009 0 0 0 0 0 0 0 0
123-103-138 0 0 1 0.002 0 0 1 0.008 0 0
125-105-128 1 0.009 0 0 0 0 0 0 0 0
123-105-128 2 0.018 0 0 0 0 0 0 0 0
113-105-130 1 0.009 3 0.008 3 0.017 0 0 0 0
123-105-130 8 0.071 64 0.178 25 0.145 24 0.194 15 0.242
125-105-130 66 0.590 35 0.097 24 0.140 5 0.040 6 0.097
127-105-130 0 0 3 0.008 3 0.017 0 0 0 0
113-105-132 1 0.009 0 0 0 0 0 0 0 0
123-105-132 1 0.009 3 0.008 0 0 2 0.016 1 0.016
125-105-132 4 0.036 1 0.002 1 0.006 0 0 0 0
127-105-132 0 0 2 0.005 2 0.012 0 0 0
123-105-134 3 0.027 9 0.025 3 0.017 4 0.032 2 0.032
125-105-134 5 0.045 4 0.011 3 0.017 1 0.008 0 0
127-105-134 0 0 1 0.002 0 0 1 0.008 0 0
123-105-136 0 0 2 0.005 0 0 2 0.016 0 0
125-105-138 0 0 1 0.002 0 0 1 0.008 0 0
123-105-138 0 0 1 0.002 1 0.006 0 0 0 0
113-107-130 0 0 1 0.002 0 0 1 0.008 0 0
123-107-130 1 0.009 16 0.044 4 0.023 4 0.032 8 0.129
125-107-130 0 0 1 0.002 1 0.006 0 0 0 0
123-107-128 0 0 1 0.002 0 0 0 0 1 0.016
125-107-128 0 0 1 0.002 0 0 1 0.008 0 0
127-107-130 0 0 2 0.005 0 0 2 0.016 0 0
123-107-132 0 0 1 0.002 0 0 1 0.008 0 0
123-107-134 0 0 1 0.002 0 0 1 0.008 0 0
125-109-126 1 0.009 0 0 0 0 0 0 0 0
123-109-130 1 0.009 16 0.044 9 0.052 4 0.032 3 0.048
125-109-130 2 0.018 7 0.019 6 0.035 1 0.008 0 0
127-109-130 0 0 1 0.002 1 0.006 0 0 0 0
123-109-132 0 0 2 0.005 1 0.006 1 0.008 0 0
125-109-132 0 0 1 0.002 1 0.006 0 0 0 0
127-109-134 1 0.009 1 0.002 0 0 0 0 1 0.016
125-109-134 0 0 1 0.002 1 0.006 0 0 0 0
123-109-134 0 0 1 0.002 1 0.006 0 0 0 0
123-109-136 1 0.009 0 0 0 0 0 0 0 0
125-109-138 0 0 1 0.002 1 0.006 0 0 0 0
123-111-130 0 0 3 0.008 0 0 3 0.024 0 0
127-111-130 0 0 2 0.005 1 0.006 0 0 1 0.016
123-111-134 0 0 1 0.002 0 0 1 0.008 0 0
125-111-134 0 0 1 0.002 0 0 0 0 1 0.016
123-113-130 2 0.018 8 0.022 6 0.035 1 0.008 1 0.016
125-113-130 1 0.009 3 0.008 2 0.012 1 0.008 0 0
123-113-134 0 0 1 0.002 0 0 0 0 1 0.016
Haplotypic diversity 0.645 0.943 0.944 0.917 0.946
Total, chromosomes 112 358 172 124 62
Number of haplotype variants 25 52 32 32 16
Table 5. Association and linkage disequilibrium of mark- Table 6. The number of generations which passed since
ers D13S141, D13S175, and D13S143 with the mutation the time of the c.35delG mutation spread in the Volga-
o35delG Ural region
Marker Allele P (95% significance level) x2 А
D13S141 125 <0.001 67.05872 0.636179
D13S175 105 <0.001 47.8665 0.666045
D13S143 130 0.81 0.08801 0.062381
Marker Number of generation since the time of the mutation’s appearance in the population (q) Number of years which passed from the start of expansion Start of expansion
D13S175 (allele 105) 470 11800 9800 B.C.
D13S141 (allele 125) 133 3300 1300 B.C.
Mean value 301 7500 5500 B.C.
(172 analyzed chromosomes), 32 of the 59 haplotypes were revealed in the total sample; in Tatars (124 chromosomes analyzed), 32 of the 59; and in Bashkirs (62 chromosomes), 17 of the 59 haplotypes.
In Russians, the most widespread haplotypes were 123-105-130 (14.5%), 125-105-130 (14.0%), and 123-103130 (7.6%), while the remaining 29 occurred at different frequencies: from 0.6% (12 haplotypes) to 7.0% (haplotype 125-103-130). On normal chromosomes in Tatars, the haplotypes 123-105-130 and 123-103-130 occurred most often at frequencies of 19.4 and 9.7%, respectively, and the frequencies of the remaining 30 varied from 0.8% (19 haplotypes) to 5.6% (haplotype 123-101-130). In Bashkirs, the haplotypes 123-105-130 (24.2%), 123107-130 (12.9%) and 125-105-130 (9.7%) were most often recorded, and their total frequency reached 46.8%. The frequencies of the remaining 14 haplotypes varied from 1.6% (9 haplotypes) to 6.4% (haplotype 125-103130). Besides, in each ethnic group, specific haplotypes (D13S141-D13S175-D13S143), non-occurring in other groups, were observed at low frequencies: in Russians, 14 haplotypes; in Tatars, 15; and in Bashkirs, 4.
An analysis of the haplotypes D13S141-D13S175-D13S143 revealed that there is a higher frequency of haplotype 125-105-130 on chromosomes with the c.35delG mutation compared to the chromosomes of healthy donors (x2 = 64.866, p < 0.001), as well as an ethnic specificity of the spectra and frequencies of occurrence of the haplotypes D13S141-D13S175-D13S143 in three ethnic groups of healthy donors.
Table 5 shows data on the association and linkage disequilibrium of marker D13S175, D13S143, and D13S141 alleles carrying the mutation c.35delG.
The values of the linkage disequilibrium parameter were greatest in alleles 105 and 125 of D13S175 and D13S141 markers, respectively, located proximally relative to the gene GJB2, and the lowest in alleles of the distal D13S143 marker. Stemming from the values of x2 and the linkage disequilibrium parameter Ô, the most
probable founder haplotype (the ancestor haplotype) seems to consist of the alleles 125-105-130 (Fig. 4). The haplotype 125-105-130 was revealed on 59% of all chromosomes carrying c.35delG, which is reliably higher (p < 0.001) than the frequency of this haplotype (9.7%) on normal chromosomes.
To calculate the number of generations since the start of the spread of the c.35delG mutation in the populations of the Volgo-Ural region, we selected two markers: D13S175 and D13S141. The selection criterion for these markers was relatively high values of x2 and measure of linkage disequilibrium Ô; statistically significant differences in the frequency distribution of these marker alleles in chromosomes with and without c.35delG were also considered.
No statistically significant differences in the frequency distribution of alleles of the marker D13S143 between c.35delG and intact chromosomes (Table. 3), and concurrently a minimum value of the linkage disequilibrium parameter, were observed (Table 5). The prevalence of the allele 130 of this STR-mark-er in two groups of chromosomes (with and without the c.35delG mutation) formally explains the absence of statistically significant differences. Nevertheless, given the statistically significant differences obtained for two other STR-markers, which are closer to the c.35delG mutation, the existence of the ancestor haplotype for c.35delG, embracing the area covered by the D13S141-D13S175-D13S143 markers, and its subsequent “tailing” at the expense of recombination and mutation events during the numerous generations seem likely. However, the absence of a statistically significant association of the most frequent allele 130 (D13S143) and the ancestor haplotype provided grounds for excluding D13S143 from markers using which the number of generations was computed (Table 6).
After the beginning of divergence of the c.35delG-mutant ancestor haplotype in the populations of the
Volga-Ural region, from 133 to 470 generations (on average 301 generations) passed. At the estimate of the age (in years) of the ancestor haplotype reconstructed on the territory of the Volga-Ural region, the duration of one generation, as in other studies, was considered equal to 25 years (Table 6).
The time over which the expansion of the c.35delG-mutant chromosomes has been present within the population of the Volga-Ural region ranges from 3,300-
11,800 years (mean is ~ 7,500). However, such estimates of the number of generations (based on the physical distance) often lead to the overestimation of the “age” of the mutation, because it is the probable (not the observed) value for the mutation events that is considered; therefore, with this approach, researchers do not orientate mean values. Instead, they choose the most distant marker, which is linked however to the locus of the disease, i.e. the so-called boundary of the stable haplotype [47]. In the given case, it is the D13S175 marker. If one presumes that the number of generations calculated with the use of this marker is more accurate, then the most likely time of expansion of the founder haplotype with the mutation c.35delG in the populations of the Volga-Ural region is ~ 11,800 years. Such dating of the beginning of c.35delG expansion in the Volga-Ural region matches the results obtained in studies with the use of different DNA markers (SNP-and STR-markers) in other populations of Eurasia (10,000-14,000 years ago) [9-12].
The results of the haplotype analysis, estimates of the age of the c.35delG mutation, and data on the descending gradient to fit carrier frequency from south to north in the European populations allow to suggest that the Middle East and Mediterranean regions (probably, modern Greece) are the most probable centers of c.35delG origin, from where, together with Neolithic migrations of man, it widely propagated throughout Europe [9-12]. An analysis of the haplotypic diversity (using STR-markers) and approximate assessment of the age of c.35delG in the Volga-Ural region favor to a great extent the “traditional” Neolithic hypothesis about the origin and spread of this mutation.
Considering the unified time continuum of the start of c.35delG spread in Eurasia, obtained with the use of various systems of DNA markers (SNP- and STR-markers), an assessment of the world haplotypic diversity of the mutant chromosomes accounting for the integral set of DNA markers is required to get an unambiguous answer to the question of the center of origin of the c.35delG mutation.
CONCLUSIONS
The analysis of the carrier frequency of the mutation c.35delG in the GJB2 gene in Eurasian populations has revealed a tendency towards gradient decrease in the c.35delG frequency from West to East, starting from the populations of Eastern Europe and the Volga-Ural region, with average frequencies of 3.3 and 1.4%, respectively; a low frequency (0.8—0.9%) in Central Asian populations, a minimum frequency (0.4%) in Yakuts residing in Eastern Siberia; and the absence of c.35delG in Altaians (Southern Siberia).
A haplotype analysis of the c.35delG-mutant chromosomes has allowed us to reconstruct the ancestor haplo-type carrying this mutation and to confirm the unified origin of most of the studied mutant chromosomes of NSHL patients living in the Volga-Ural region. The estimated time of expansion of the c.35delG mutation carriers we obtained (11,800 years ago) fits the world values of the “age” of this mutation (10,000-14,000 years).
The body of data collected should help clarify or review existing ideas about the center and time of origination of the c.35delG mutation (GJB2), as well as the factors that define its occurrence worldwide. •
This work was supported by a Grant of the Russian Foundation for Basic Research (№ 09-04-01123-а), Federal Target-Oriented Programs “Scientific and Scientific-Pedagogical Personnel of Innovative Russia in 2009-2013” (№ 16.740.11.0190, 16.740.11.0346) and for 2010-2012 (№ 02.740.11.0701), as well as by the Russian Ministry of Education and Science (Government Contracts P325 and P601) and Government Contract № 16.512.11.2047.
REFERENCES
1. Marazita M.L., Ploughman L.M., Rawlings B., Remington
E., Arnos K.S., Nance W.E. // Am. J. Med. Genet. 1993. V. 46. № 5. P. 486-491.
2. van Camp G., Smith R. // 2009. http://webhost.ua.ac.be/ hhh/
3. Maeda S., Nakagawa S., Suga M., Yamashita E., Oshima A., Fujiyoshi Y., Tsukihara T. // Nature. 2009. V. 458. №7238. P. 597-602.
4. Kikuchi T., Kimura R.S., Paul D.L., Adams J.C. // Anat. Embryol. (Berl.). 1995. V. 191(2). P. 101-118.
5. Morell R.J., Friderici K.H., Wei S., Elfenbein J.L., Friedman T.B., Fisher R.A. // N. Engl. J. Med. 1998. V. 339. P. 1500-1505.
6. van Laer L., Huizing E.H., Verstreken M., van Zuijlen D., Wauters J.G., Bossuyt P. J., van de Heyning P., McGuirt W.T., Smith R.J. // J. Med. Genet. 2001. V. 38. P. 515-518.
7. Yan D., Ke X., Blanton S.H., Ouyang X.M., Pandya A., Du L.L., Nance W.E., Liu X.Z. // Hum. Genet. 2003. V. 114. P. 44-50.
8. RamShankar M., Girirajan S., Dagan O., Ravi Shankar H.M., Jalvi R., Rangasayee R., Avraham K.B., Anand A. // J. Med. Genet. 2003. V. 40. P. 68.
9. Abidi O., Boulouiz R., Nahili H., Ridal M., Noureddine A.M., Tlili A., Rouba H., Masmoudi S., Chafik A., Hassar M., et al. // Genet. Test. Mol. Biomark. 2007. V. 12. №4. P. 569-574.
10. Kokotas H., Grigoriadou M., Villamar M., Giannoulia-Karantana A., del Castillo I., Petersen B.M. // Genet. Test. Mol. Biomark. 2010. V. 14. №2. P. 183-187.
11. Najmabadi H., Cucci R., Sahebjam S., Kouchakian N., Far-hadi M., Kahrizi K., Arzhangi S., Daneshmandan N., Javan K., Smith R.J.H. // Hum. Mutat. 2002. V. 504. P. 135-138.
12. Kokotas H., van Laer L., Grigoriadou M., Ferekidou E., Papadopoulou E., Neou P., Giannoulia-Karantana A., Kan-diloros D., Korres S., Petersen M.B. // Am. J. Med. Genet. 2008. V. 146A. P. 2879-2884.
13. Brobby G., Muller-Myhsok B., Horstmann R. // N. Engl. J. Med. 1998. V. 338 (8). P. 548-550.
14. Hamelmann C., Amedofu G.K., Albrecht K., Muntau B., Gelhaus A., Brobby G.W., Horstmann R.D. // Hum. Mutat. 2001. V. 18. № 1. P. 84-85.
15. Anichkina A., Kulenich T., Zinchenko S., Shagina I., Polyakov A., Ginter E., Evgrafov O., Viktorova T., Khusnit-donova E. // Eur. J. Hum. Genet. 2001. V. 9. P. 151.
16. Nekrasova N.U., Shagina I.A., Petrin A.N., Polyakov A.V.
// Med. genet. 2002. V. 1. № 6. P. 290-294.
17. Zinchenko R.A., Elchinova G.I., Barishnikova N.V., Polyakov A.V., Ginter E.K. // Russian Journal of Genetics. 2007. V.
43. № 9. P. 1246-1254.
18. Zinchenko R.A., Zinchenko S.P., Galkina V.A. // Russian Journal of Genetics. 2003. V. 39. № 9. P. 1275-1284.
19. Osetrova A.A., Sharonova E.I., Rossinskaya T.G., Zinchenko R.A. // Med. genet. 2010. № 9. P. 30-40.
20. Khidyatova I.M, Dzhemileva L.U., Khabibulin R.M., Khusnutdinova E.K. // Mol. biol. 2002. V. 36. № 3. P. 438441.
21. Khusnutdinova E.K., Dzhemileva L.U. // Vestn. Bioteh. and Physico-Chemical Biolog. 2005. № 1. P. 24-31.
22. Dzhemileva L.u., Barashkov N.A., Posukh O.L., Khu-sainova R.I., Akhmetov V.L., Kutuev I.A., Tadinova V.N., Fedorova S.A., Khidyatova I. M., Khusnutdinova E.K. // Med. Genet. 2009. № 8. P. 20-28.
23. Dzhemileva L. u., Posukh O. L., Tazetdinov A. M., Barashkov N. A., Zhuravskii S. A., Ponidelko S. N., Markova T. G., Tadinova V. N., Fedorova S. A., Maksimova N. R., Khusnutdinova E. K. // Russian Journal of Genetics. 2009. № 7. P. 982-991..
24. Dzhemileva L.u., Barashkov N.A., Posukh O.L., Khu-sainova R.I., Akhmetova V.L., Kutuev I.A., Gilyazova I.R., Tadinova V.N., Fedorova S.A., Khidiyatova I.M., et al. // J. Hum. Genet. 2010. V.55. № 11. P. 749-754.
25. Posukh O.L., Pallares-Ruiz N., Tadinova V., Osipova L., Claustres M., Roux A.F. // BMC Med. Genet. 2005. V. 6. № 12. P. 1-7.
26. Barashkov N.A., Dzhemileva L.U., Fedorova S. A., Maksimova N.R., Khusnutdinova E. K. // Vest. otorin. 2008. № 5. P. 23-28.
27. Barashkov N.A., Dzhemileva L.U., Fedorova S. A., Terutin
F.M., Fedorova E.E., Gurinova E.E., Alekseeva S.P., Kononova S.K., Nogovicina A.N., Khusnutdinova E.K. // Med. genet. 2010. V. 9. № 7 (97). P. 22-33.
28. Shokarev R.A., Amelina S.S., Kriventsova N.V. // Med. genet. 2005. V. 4. № 12. P. 556- 567.
29. Shokarev R.A., Amelina S.S., Zinchenko S.N., Elchinova
G.I., Khlebnikova O.O., Blisnez E.A., Tverskaya S.M., Polyakov A.V., Zinchenko R.A. // Med. genet. 2006. V. 5. P. 38-43.
30. Boghkova V.P., Hashaev Z.H., Umankaya T.A. // Biophis. 2010. V. 55. № 3. P.514-525.
31. Tavartkiladze G.A., Polyakov A.V., Markova T.G., Lalay-anz M.R., Blisnez E.A. // Vestn. otorin. 2010. № 3. P.1-18.
32. Tazetdinov A.M., Dzhemileva L.U., Khusnutdinova E.K.
// Russian Journal of Genetics. 2008. V. 44. № 6. P. 725-733.
33. Denoyelle F., Weil D., Maw M. // Hum. Mol. Genet. 1997. V.
12. № 6. P. 2173-2177.
34. Petrukhin K.E., Speer M.C., Cayanis E., Bonaldo M.F., Tantravahi U., Soares M.B., Fischer S.G., Warburton D., Gilliam T.C., Ott J. // Genomics. 1993. V. 15. № 1. P. 76-85.
35. Brown K.A., Janura A., Karbani G., Parrys G., Noble A., Crockford G., Bishop D.T., Newton V.E., Markham A.F., Mueller R.F. // Hum. Mol. Genet. 1996. V. 5. P. 169-173.
36. Rothrock C.R., Murgia A., Sartorato E.L., Leonardi E.,
Wei S., Lebeis S.L., Yu L.E., Elfenbein J.L., Fisher R.A., Friderici K.H. // Hum. Genet. 2003. V. 113. P. 18-23.
37. Bengtsson B.O., Thompson G. // Tissue Antigens. 1981. V.
18. P. 356-363.
38. Risch N., de Leon D., Ozelius L. // Nat. Genet. 1995. V. 9.
№ 2. P. 152-159.
39. Krawczak M., Konecki D.S., Schmidtke I. // Hum. Genet. 1988. V. 80. P. 78-80.
40. Nei M. Molecular Evolutionary Genetics. New York: Columbia University Press, 1987. 275 p.
41. Petersen M., Willems P. // Clin. Genet. 2006. V. 69. P 371-392.
42. Gasparini P., Rabionet R., Barbujani G., Melchionda S., Petersen M., Brondum-Nielsen K., Metspalu A., Oitmaa E., Pisano M., Fortina P., et al. // Eur. J. Hum. Genet. 2000. V. 8. №1. P. 19-23.
43. Lucotte G. // Int. J. Pediatr. Otorhinolaryngol. 2007. V. 71. P. 741-746.
44. Mahdieh N., Rabbani B. // Int. J. Audiol. 2009. V. 48. P. 363-370.
45. Shahin H., Walsh T., Sobe T., Lynch E., •King M.-C., Avraham K.B., Kanaan M. // Hum. Genet. 2002. V. 110. P. 284-289.
46. Belguith H., Hajji S., Salem N., Charfeddine I., Lahmar I., Amor M.B., Ouldim K., Chouery E., Driss N., Drira M., et al. // Clin. Genet. 2005. V. 68. P. 188-189.
47. Slatkin M., Rannala B. // Annu. Rev. Genomics Hum. Genet. 2000. V. 1. P. 225-249.
