Unique Mutational Spectrum of the GJB2 Gene and its Pathogenic Contribution to Deafness in Tuvinians (Southern Siberia, Russia): A High Prevalence of Rare Variant c.516G>C (p.Trp172Cys)

Functional Consequences of Pathogenic Variants of the GJB2 Gene (Cx26) Localized in Different Cx26 Domains

One of the most common forms of genetic deafness has been predominantly associated with pathogenic variants in the GJB2 gene, encoding transmembrane protein connexin 26 (Cx26). The Cx26 molecule consists of an N-terminal domain (NT), four transmembrane domains (TM1-TM4), two extracellular loops (EL1 and EL2), a cytoplasmic loop, and a C-terminus (CT). Pathogenic variants in the GJB2 gene, resulting in amino acid substitutions scattered across the Cx26 domains, lead to a variety of clinical outcomes, including the most common non-syndromic autosomal recessive deafness (DFNB1A), autosomal dominant deafness (DFNA3A), as well as syndromic forms combining hearing loss and skin disorders. However, for rare and poorly documented variants, information on the mode of inheritance is often lacking. Numerous in vitro studies have been conducted to elucidate the functional consequences of pathogenic GJB2 variants leading to amino acid substitutions in different domains of Cx26 protein. In this work, we summarized all available data on a mode of inheritance of pathogenic GJB2 variants leading to amino acid substitutions and reviewed published information on their functional effects, with an emphasis on their localization in certain Cx26 domains.

The GJB2 (Cx26) Gene Variants in Patients with Hearing Impairment in the Baikal Lake Region (Russia)

The GJB2 (Cx26) gene pathogenic variants are associated with autosomal recessive deafness type 1A (DFNB1A, OMIM #220290). Direct sequencing of the GJB2 gene among 165 hearing-impaired individuals living in the Baikal Lake region of Russia identified 14 allelic variants: pathogenic/likely pathogenic-nine variants, benign-three variants, unclassified-one variant, and one novel variant. The contribution of the GJB2 gene variants to the etiology of hearing impairment (HI) in the total sample of patients was 15.8% (26 out of 165) and significantly differed in patients of different ethnicity (5.1% in Buryat patients and 28.9% in Russian patients). In patients with DFNB1A (n = 26), HIs were congenital/early onset (92.3%), symmetric (88.5%), sensorineural (100.0%), and variable in severity (moderate-11.6%, severe-26.9% or profound-61.5%). The reconstruction of the SNP haplotypes with three frequent GJB2 pathogenic variants (c.-23+1G>A, c.35delG or c.235delC), in comparison with previously published data, supports a major role of the founder effect in the expansion of the c.-23+1G>A and c.35delG variants around the world. Comparative analysis of the haplotypes with c.235delC revealed one major haplotype G A C T (97.5%) in Eastern Asians (Chinese, Japanese and Korean patients) and two haplotypes, G A C T (71.4%) and G A C C (28.6%), in Northern Asians (Altaians, Buryats and Mongols). The variable structure of the c.235delC-haplotypes in Northern Asians requires more studies to expand our knowledge about the origin of this pathogenic variant.

Evolutionary origin of pathogenic GJB2 alleles in China

The frequency of the pathogenic allele of the autosomal recessive deafness gene GJB2 varies among different populations in the world, and accumulates to a sufficiently high frequency in certain population. The purpose of this study is to investigate the origin and evolution of GJB2 pathogenic alleles in Chinese deaf patients. Children with non-syndromic hearing loss, and their parents, from 295 families were recruited. Customized capture probes targeted at 943 SNPs related to GJB2 gene were designed for sequencing of genomic DNA in blood samples. Haplotypes carrying pathogenic allele were analyzed through linkage disequilibrium block building, ancestry tracing, and extended haplotype heterozygosity calculation. Two pathogenic GJB2 alleles, c.235delC (18.41%) and c.109G>A (15.57%), were observed in 867 donors. For c.235delC allele, 3 different core haplotypes with one major haplotype (97.32%) were found, and their core SNPs were 100% conserved. For c.109G>A allele, 6 different haplotypes with one major haplotype (93.28%) were found and the major c.109G>A allele evolved from a specific ancestral haplotype. Geographical origins of donors carrying GJB2 c.109G>A and c.235delC core haplotypes centered between Qinghai and Neimenggu. GJB2 c.235delC has long-range linkage disequilibrium. No positive selection signature was found for GJB2 c.235delC or c.109G>A in the studied population. In conclusion, we discovered a single origin of GJB2 c.235delC allele and multiple independent origins of GJB2 c.109G>A allele. Alternative to positive selection or multiple independent recurrent mutation event, population bottleneck effect might account for the observed high population frequency of these pathogenic alleles.

Agent-Based Modeling of Autosomal Recessive Deafness 1A (DFNB1A) Prevalence with Regard to Intensity of Selection Pressure in Isolated Human Population

An increase in the prevalence of autosomal recessive deafness 1A (DFNB1A) in populations of European descent was shown to be promoted by assortative marriages among deaf people. Assortative marriages became possible with the widespread introduction of sign language, resulting in increased genetic fitness of deaf individuals and, thereby, relaxing selection against deafness. However, the effect of this phenomenon was not previously studied in populations with different genetic structures. We developed an agent-based computer model for the analysis of the spread of DFNB1A. Using this model, we tested the impact of different intensities of selection pressure against deafness in an isolated human population over 400 years. Modeling of the "purifying" selection pressure on deafness ("No deaf mating" scenario) resulted in a decrease in the proportion of deaf individuals and the pathogenic allele frequency. Modeling of the "relaxed" selection ("Assortative mating" scenario) resulted in an increase in the proportion of deaf individuals in the first four generations, which then quickly plateaued with a subsequent decline and a decrease in the pathogenic allele frequency. The results of neutral selection pressure modeling ("Random mating" scenario) showed no significant changes in the proportion of deaf individuals or the pathogenic allele frequency after 400 years.

Different Rates of the SLC26A4-Related Hearing Loss in Two Indigenous Peoples of Southern Siberia (Russia)

Hereditary hearing loss (HL) is known to be highly locus/allelic heterogeneous, and the prevalence of different HL forms significantly varies among populations worldwide. Investigation of region-specific landscapes of hereditary HL is important for local healthcare and medical genetic services. Mutations in the SLC26A4 gene leading to nonsyndromic recessive deafness (DFNB4) and Pendred syndrome are common genetic causes of hereditary HL, at least in some Asian populations. We present for the first time the results of a thorough analysis of the SLC26A4 gene by Sanger sequencing in the large cohorts of patients with HL of unknown etiology belonging to two neighboring indigenous Turkic-speaking Siberian peoples (Tuvinians and Altaians). A definite genetic diagnosis based on the presence of biallelic SLC26A4 mutations was established for 28.2% (62/220) of all enrolled Tuvinian patients vs. 4.3% (4/93) of Altaian patients. The rate of the SLC26A4-related HL in Tuvinian patients appeared to be one of the highest among populations worldwide. The SLC26A4 mutational spectrum was characterized by the presence of Asian-specific mutations c.919-2A>G and c.2027T>A (p.Leu676Gln), predominantly found in Tuvinian patients, and c.2168A>G (p.His723Arg), which was only detected in Altaian patients. In addition, a novel pathogenic variant c.1545T>G (p.Phe515Leu) was found with high frequency in Tuvinian patients. Overall, based on the findings of this study and our previous research, we were able to uncover the genetic causes of HL in 50.5% of Tuvinian patients and 34.5% of Altaian patients.

A common founder effect of the splice site variant c.-23 + 1G > A in GJB2 gene causing autosomal recessive deafness 1A (DFNB1A) in Eurasia

Mutations in the GJB2 gene are known to be a major cause of autosomal recessive deafness 1A (OMIM 220290). The most common pathogenic variants of the GJB2 gene have a high ethno-geographic specificity in their distribution, being attributed to a founder effect related to the Neolithic migration routes of Homo sapiens. The c.-23 + 1G > A splice site variant is frequently found among deaf patients of both Caucasian and Asian origins. It is currently unknown whether the spread of this mutation across Eurasia is a result of the founder effect or if it could have multiple local centers of origin. To determine the origin of c.-23 + 1G > A, we reconstructed haplotypes by genotyping SNPs on an Illumina OmniExpress 730 K platform of 23 deaf individuals homozygous for this variant from different populations of Eurasia. The analyses revealed the presence of common regions of homozygosity in different individual genomes in the sample. These data support the hypothesis of the common founder effect in the distribution of the c.-23 + 1G > A variant of the GJB2 gene. Based on the published data on the c.-23 + 1G > A prevalence among 16,177 deaf people and the calculation of the TMRCA of the modified f2-haplotypes carrying this variant, we reconstructed the potential migration routes of the carriers of this mutation around the world. This analysis indicates that the c.-23 + 1G > A variant in the GJB2 gene may have originated approximately 6000 years ago in the territory of the Caucasus or the Middle East then spread throughout Europe, South and Central Asia and other regions of the world.

Genetic etiology of hearing loss in Russia

Prevalence and locus/allelic heterogeneity of the hereditary hearing loss (HL) vary significantly in different human populations. Investigation of the hereditary HL diversity and the evaluation of the factors determining the region-specific landscapes of genetic HL are important for local healthcare and medical genetic services. This review presents the summarized data from the published studies concerning the genetic etiology of HL in different populations of Russia. Multiethnic population of Russia (in total, about 146 million on 2021) includes over 180 different ethnic groups, the number of which varies from millions to just several thousand people. Among them, Russians are the largest group (about 111 million). The contribution of GJB2 gene in the HL etiology in patients of different ethnicities and ethnic-specific prevalence of the GJB2 pathogenic variants were studied in many local populations of Russia. However, the investigation of other "deafness" genes is still limited to a relatively small number of studies on patients with HL of unsolved etiology.

Functional Evaluation of a Rare Variant c.516G>C (p.Trp172Cys) in the GJB2 (Connexin 26) Gene Associated with Nonsyndromic Hearing Loss

Mutations in the GJB2 gene encoding transmembrane protein connexin 26 (Cx26) are the most common cause for hearing loss worldwide. Cx26 plays a crucial role in the ionic and metabolic homeostasis in the inner ear, indispensable for normal hearing process. Different pathogenic mutations in the GJB2 gene can affect all stages of the Cx26 life cycle and result in nonsyndromic autosomal recessive (DFNB1) or dominant (DFNA3) deafness and syndromes associating hearing loss with skin disorders. This study aims to elucidate the functional consequences of a rare GJB2 variant c.516G>C (p.Trp172Cys) found with high frequency in deaf patients from indigenous populations of Southern Siberia (Russia). The substitution c.516G>C leads to the replacement of tryptophan at a conserved amino acid position 172 with cysteine (p.Trp172Cys) in the second extracellular loop of Cx26 protein. We analyzed the subcellular localization of mutant Cx26-p.Trp172Cys protein by immunocytochemistry and the hemichannels permeability by dye loading assay. The GJB2 knockout HeLa cell line has been generated using CRISPR/Cas9 genome editing tool. Subsequently, the HeLa transgenic cell lines stably expressing different GJB2 variants (wild type and mutations associated with hearing loss) were established based on knockout cells and used for comparative functional analysis. The impaired trafficking of mutant Cx26-p.Trp172Cys protein to the plasma membrane and reduced hemichannels permeability support the pathogenic effect of the c.516G>C (p.Trp172Cys) variant and its association with nonsyndromic hearing loss. Our data contribute to a better understanding of the role of mutations in the second extracellular loop of Cx26 protein in pathogenesis of deafness.

A new approach to estimating the prevalence of hereditary hearing loss: An analysis of the distribution of sign language users based on census data in Russia

The absence of comparable epidemiological data challenges the correct estimation of the prevalence of congenital hearing loss (HL) around the world. Sign language (SL) is known as the main type of communication of deaf people. We suggest that the distribution of SL can be interpreted as an indirect indicator of the prevalence of congenital HL. Since a significant part of congenital HL is due to genetic causes, an assessment of the distribution of SL users can reveal regions with an extensive accumulation of hereditary HL. For the first time, we analyzed the data on the distribution of SL users that became available for the total population of Russia by the 2010 census. Seventy-three out of 85 federal regions of Russia were ranked into three groups by the 25th and 75th percentiles of the proportion of SL users: 14 regions-"low proportion"; 48 regions-"average proportion"; and 11 regions-"high proportion". We consider that the observed uneven prevalence of SL users can reflect underlying hereditary forms of congenital HL accumulated in certain populations by specific genetic background and population structure. At least, the data from this study indicate that the highest proportions of SL users detected in some Siberian regions are consistent with the reported accumulation of specific hereditary HL forms in indigenous Yakut, Tuvinian and Altaian populations.

Connexin Genes Variants Associated with Non-Syndromic Hearing Impairment: A Systematic Review of the Global Burden

Mutations in connexins are the most common causes of hearing impairment (HI) in many populations. Our aim was to review the global burden of pathogenic and likely pathogenic (PLP) variants in connexin genes associated with HI. We conducted a systematic review of the literature based on targeted inclusion/exclusion criteria of publications from 1997 to 2020. The databases used were PubMed, Scopus, Africa-Wide Information, and Web of Science. The protocol was registered on PROSPERO, the International Prospective Register of Systematic Reviews, with the registration number “CRD42020169697”. The data extracted were analyzed using Microsoft Excel and SPSS version 25 (IBM, Armonk, New York, United States). A total of 571 independent studies were retrieved and considered for data extraction with the majority of studies (47.8% (n = 289)) done in Asia. Targeted sequencing was found to be the most common technique used in investigating connexin gene mutations. We identified seven connexin genes that were associated with HI, and GJB2 (520/571 publications) was the most studied among the seven. Excluding PLP in GJB2, GJB6, and GJA1 the other connexin gene variants (thus GJB3, GJB4, GJC3, and GJC1 variants) had conflicting association with HI. Biallelic GJB2 PLP variants were the most common and widespread variants associated with non-syndromic hearing impairment (NSHI) in different global populations but absent in most African populations. The most common GJB2 alleles found to be predominant in specific populations include; p.Gly12ValfsTer2 in Europeans, North Africans, Brazilians, and Americans; p.V37I and p.L79Cfs in Asians; p.W24X in Indians; p.L56Rfs in Americans; and the founder mutation p.R143W in Africans from Ghana, or with putative Ghanaian ancestry. The present review suggests that only GJB2 and GJB3 are recognized and validated HI genes. The findings call for an extensive investigation of the other connexin genes in many populations to elucidate their contributions to HI, in order to improve gene-disease pair curations, globally.

High Rates of Three Common GJB2 Mutations c.516G>C, c.-23+1G>A, c.235delC in Deaf Patients from Southern Siberia Are Due to the Founder Effect

The mutations in the GJB2 gene (13q12.11, MIM 121011) encoding transmembrane protein connexin 26 (Cx26) account for a significant portion of hereditary hearing loss worldwide. Earlier we found a high prevalence of recessive GJB2 mutations c.516G>C, c.-23+1G>A, c.235delC in indigenous Turkic-speaking Siberian peoples (Tuvinians and Altaians) from the Tyva Republic and Altai Republic (Southern Siberia, Russia) and proposed the founder effect as a cause for their high rates in these populations. To reconstruct the haplotypes associated with each of these mutations, the genotyping of polymorphic genetic markers both within and flanking the GJB2 gene was performed in 28 unrelated individuals homozygous for c.516G>C (n = 18), c.-23+1G>A (n = 6), or c.235delC (n = 4) as well as in the ethnically matched controls (62 Tuvinians and 55 Altaians) without these mutations. The common haplotypes specific for mutations c.516G>C, c.-23+1G>A, or c.235delC were revealed implying a single origin of each of these mutations. The age of mutations estimated by the DMLE+ v2.3 software and the single marker method is discussed in relation to ethnic history of Tuvinians and Altaians. The data obtained in this study support a crucial role of the founder effect in the high prevalence of GJB2 mutations c.516G>C, c.-23+1G>A, c.235delC in indigenous populations of Southern Siberia.