Relationships between Uncoupling Protein Genes UCP1, UCP2 and UCP3 and Irisin Levels in Residents of the Coldest Region of Siberia.

Currently, it is known that irisin can participate in the processes of thermoregulation and browning of adipose tissue, and, therefore, it is possible that it is involved in the microevolutionary mechanisms of adaptation to a cold. The aim of this study is to investigate the relationship between the uncoupling protein genes (UCP1, UCP2, UCP3) and the irisin levels in the residents of the coldest region of Siberia. The sample consisted of 279 Yakut people (185 females, 94 males, average age 19.8 ± 2.03 years). The females plasma irisin concentration was 8.33 ± 2.74 mcg/mL and the males was 7.76 ± 1.86 mcg/mL. Comparative analysis of irisin levels with the genotypes of six studied SNP-markers in females revealed a significant association of irisin with rs1800849-UCP3. The TT genotype of rs1800849 was associated with elevated levels of irisin (p = 0.01). It was also found that this TT genotype in females was associated with reduced weight and height (p = 0.03). We searched for natural selection signals for the T-allele rs1800849-UCP3; as a result of which, it was found that this allele has a significantly high frequency of distribution in northern (45%, CI: 0.42-0.484) compared with southern Asian populations (28%, CI: 0.244-0.316) (p = 0.01). The results obtained indicate the probable involvement of irisin and the UCP3 gene in thermoregulation, and the spread of its allelic variants is probably related to adaptation to a cold climate.

The Role of Nonshivering Thermogenesis Genes on Leptin Levels Regulation in Residents of the Coldest Region of Siberia.

Leptin plays an important role in thermoregulation and is possibly associated with the microevolutionary processes of human adaptation to a cold climate. In this study, based on the Yakut population (n = 281 individuals) living in the coldest region of Siberia (t°minimum -71.2 °C), we analyze the serum leptin levels and data of 14 single nucleotide polymorphisms (SNPs) of 10 genes (UCP1, UCP2, UCP3, FNDC5, PPARGC1A, CIDEA, PTGS2, TRPV1, LEPR, BDNF) that are possibly involved in nonshivering thermogenesis processes. Our results demonstrate that from 14 studied SNPs of 10 genes, 2 SNPs (the TT rs3811787 genotype of the UCP1 gene and the GG rs6265 genotype of the BDNF gene) were associated with the elevated leptin levels in Yakut females (p < 0.05). Furthermore, of these two SNPs, the rs3811787 of the UCP1 gene demonstrated more indications of natural selection for cold climate adaptation. The prevalence gradient of the T-allele (rs3811787) of UCP1 increased from the south to the north across Eurasia, along the shore of the Arctic Ocean. Thereby, our study suggests the potential involvement of the UCP1 gene in the leptin-mediated thermoregulation mechanism, while the distribution of its allelic variants is probably related to human adaptation to a cold climate.

Autosomal recessive cataract (CTRCT18) in the Yakut population isolate of Eastern Siberia: a novel founder variant in the FYCO1 gene.

Congenital autosomal recessive cataract with unknown genetic etiology is one of the most common Mendelian diseases among the Turkic-speaking Yakut population (Eastern Siberia, Russia). To identify the genetic cause of congenital cataract spread in this population, we performed whole-exome sequencing (Illumina NextSeq 500) in one Yakut family with three affected siblings whose parents had preserved vision. We have revealed the novel homozygous c.1621C>T transition leading to premature stop codon p.(Gln541*) in exon 8 of the FYCO1 gene (NM_024513.4). Subsequent screening of c.1621C>T p.(Gln541*) revealed this variant in a homozygous state in 25 out of 29 Yakut families with congenital cataract (86%). Among 424 healthy individuals from seven populations of Eastern Siberia (Russians, Yakuts, Evenks, Evens, Dolgans, Chukchi, and Yukaghirs), the highest carrier frequency of c.1621C>T p.(Gln541*) was found in the Yakut population (7.9%). DNA samples of 25 homozygous for c.1621C>T p.(Gln541*) patients with congenital cataract and 114 unaffected unrelated individuals without this variant were used for a haplotype analysis based on the genotyping of six STR markers (D3S3512, D3S3685, D3S3582, D3S3561, D3S1289, and D3S3698). The structure of the identified haplotypes indicates a common origin for all of the studied mutant chromosomes bearing c.1621C>T p.(Gln541*). The age of the с.1621C>T p.(Gln541*) founder haplotype was estimated to be approximately 260 ± 65 years (10 generations). These findings characterize Eastern Siberia as the region of the world with the most extensive accumulation of the unique variant c.1621C>T p.(Gln541*) in the FYCO1 gene as a result of the founder effect.

Comparison of Predictive In Silico Tools on Missense Variants in GJB2, GJB6, and GJB3 Genes Associated with Autosomal Recessive Deafness 1A (DFNB1A).

In silico predictive software allows assessing the effect of amino acid substitutions on the structure or function of a protein without conducting functional studies. The accuracy of in silico pathogenicity prediction tools has not been previously assessed for variants associated with autosomal recessive deafness 1A (DFNB1A). Here, we identify in silico tools with the most accurate clinical significance predictions for missense variants of the GJB2 (Cx26), GJB6 (Cx30), and GJB3 (Cx31) connexin genes associated with DFNB1A. To evaluate accuracy of selected in silico tools (SIFT, FATHMM, MutationAssessor, PolyPhen-2, CONDEL, MutationTaster, MutPred, Align GVGD, and PROVEAN), we tested nine missense variants with previously confirmed clinical significance in a large cohort of deaf patients and control groups from the Sakha Republic (Eastern Siberia, Russia): Сх26: p.Val27Ile, p.Met34Thr, p.Val37Ile, p.Leu90Pro, p.Glu114Gly, p.Thr123Asn, and p.Val153Ile; Cx30: p.Glu101Lys; Cx31: p.Ala194Thr. We compared the performance of the in silico tools (accuracy, sensitivity, and specificity) by using the missense variants in GJB2 (Cx26), GJB6 (Cx30), and GJB3 (Cx31) genes associated with DFNB1A. The correlation coefficient (r) and coefficient of the area under the Receiver Operating Characteristic (ROC) curve as alternative quality indicators of the tested programs were used. The resulting ROC curves demonstrated that the largest coefficient of the area under the curve was provided by three programs: SIFT (AUC = 0.833, p = 0.046), PROVEAN (AUC = 0.833, p = 0.046), and MutationAssessor (AUC = 0.833, p = 0.002). The most accurate predictions were given by two tested programs: SIFT and PROVEAN (Ac = 89%, Se = 67%, Sp = 100%, r = 0.75, AUC = 0.833). The results of this study may be applicable for analysis of novel missense variants of the GJB2 (Cx26), GJB6 (Cx30), and GJB3 (Cx31) connexin genes.

Spectrum and Frequency of the GJB2 Gene Pathogenic Variants in a Large Cohort of Patients with Hearing Impairment Living in a Subarctic Region of Russia (the Sakha Republic).

Pathogenic variants in the GJB2 gene, encoding connexin 26, are known to be a major cause of hearing impairment (HI). More than 300 allelic variants have been identified in the GJB2 gene. Spectrum and allelic frequencies of the GJB2 gene vary significantly among different ethnic groups worldwide. Until now, the spectrum and frequency of the pathogenic variants in exon 1, exon 2 and the flanking intronic regions of the GJB2 gene have not been described thoroughly in the Sakha Republic (Yakutia), which is located in a subarctic region in Russia. The complete sequencing of the non-coding and coding regions of the GJB2 gene was performed in 393 patients with HI (Yakuts-296, Russians-51, mixed and other ethnicities-46) and in 187 normal hearing individuals of Yakut (n = 107) and Russian (n = 80) populations. In the total sample (n = 580), we revealed 12 allelic variants of the GJB2 gene, 8 of which were recessive pathogenic variants. Ten genotypes with biallelic recessive pathogenic variants in the GJB2 gene (in a homozygous or a compound heterozygous state) were found in 192 out of 393 patients (48.85%). We found that the most frequent GJB2 pathogenic variant in the Yakut patients was c.-23+1G>A (51.82%) and that the second most frequent was c.109G>A (2.37%), followed by c.35delG (1.64%). Pathogenic variants с.35delG (22.34%), c.-23+1G>A (5.31%), and c.313_326del14 (2.12%) were found to be the most frequent among the Russian patients. The carrier frequencies of the c.-23+1G>A and с.109G>A pathogenic variants in the Yakut control group were 10.20% and 2.80%, respectively. The carrier frequencies of с.35delG and c.101T>C were identical (2.5%) in the Russian control group. We found that the contribution of the GJB2 gene pathogenic variants in HI in the population of the Sakha Republic (48.85%) was the highest among all of the previously studied regions of Asia. We suggest that extensive accumulation of the c.-23+1G>A pathogenic variant in the indigenous Yakut population (92.20% of all mutant chromosomes in patients) and an extremely high (10.20%) carrier frequency in the control group may indicate a possible selective advantage for the c.-23+1G>A carriers living in subarctic climate.

Age-Related Hearing Impairment (ARHI) associated with GJB2 single mutation IVS1+1G>A in the Yakut population isolate in Eastern Siberia.

Age-Related Hearing Impairment (ARHI) is one of the frequent sensory disorders registered in 50% of individuals over 80 years. ARHI is a multifactorial disorder due to environmental and poor-known genetic components. In this study, we present the data on age-related hearing impairment of 48 heterozygous carriers of mutation IVS1+1G>A (GJB2 gene) and 97 subjects with GJB2 genotype wt/wt in the Republic of Sakha/Yakutia (Eastern Siberia, Russia). This subarctic territory was found as the region with the most extensive accumulation of mutation IVS1+1G>A in the world as a result of founder effect in the unique Yakut population isolate. The GJB2 gene resequencing and detailed audiological analysis in the frequency range 0.25, 0.5, 1.0, 2.0, 4.0, 8.0 kHz were performed in all examined subjects that allowed to investigate genotype-phenotype correlations between the presence of single mutation IVS1+1G>A and hearing of subjects from examined groups. We revealed the linear correlation between increase of average hearing thresholds at speech frequencies (PTA0.5,1.0,2.0,4.0 kHz) and age of individuals with GJB2 genotype IVS1+1G>A/wt (rs = 0.499, p = 0.006860 for males and rs = 0.427, p = 0.000277 for females). Moreover, the average hearing thresholds on high frequency (8.0 kHz) in individuals with genotype IVS1+1G>A/wt (both sexes) were significantly worse than in individuals with genotype wt/wt (p<0.05). Age of hearing loss manifestation in individuals with genotype IVS1+1G>A/wt was estimated to be ∼40 years (rs = 0.504, p = 0.003). These findings demonstrate that the single IVS1+1G>A mutation (GJB2) is associated with age-related hearing impairment (ARHI) of the IVS1+1G>A carriers in the Yakuts.