In silicopredictive software allows assessing the effect of amino acid substitutions on the structure or function of a protein without conducting functional studies. The accuracy ofin silicopathogenicity prediction tools has not been previously assessed for variants associated with autosomal recessive deafness 1A (DFNB1A). Here, we identifyin silicotools with the most accurate clinical significance predictions for missense variants of theGJB2(Cx26),GJB6(Cx30), andGJB3(Cx31) connexin genes associated with DFNB1A. To evaluate accuracy of selectedin silicotools (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 thein silicotools (accuracy, sensitivity, and specificity) by using the missense variants inGJB2(Cx26),GJB6(Cx30), andGJB3(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 theGJB2(Cx26),GJB6(Cx30), andGJB3(Cx31) connexin genes.

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.

Thyroid hormones play a crucial role in regulating normal development, growth, and metabolic function. However, the controversy surrounding seasonal changes in free triiodothyronine (FT3) levels remains unresolved. Therefore, the aim of this study was to conduct a systematic review and meta-analysis of variations in FT3 levels in relation to seasonal air temperatures in the context of current knowledge about its role in nonshivering thermogenesis. Ten eligible articles with a total of 336,755 participants were included in the meta-analysis. The studies were categorized into two groups based on the air temperature: “Cold winter”, where the winter temperature fell below 0 °C, and “Warm winter”, where the winter temperature was above 0 °C. The analysis revealed that in cold regions, FT3 levels decreased in winter compared to summer (I² = 57%, p < 0.001), whereas in warm regions, FT3 levels increased during winter (I² = 28%, p < 0.001). These findings suggest that seasonal variations in FT3 levels are likely to be influenced by the winter temperature. Considering the important role of the FT3 in the nonshivering thermogenesis process, we assume that this observed pattern is probably related to the differences in use of thyroid hormones in the brown adipose tissue during adaptive thermogenesis, which may depend on intensity of cold exposure.

Currently, adipose tissue is considered an endocrine organ that produces hormone-active substances, including leptin, which can play a key role in thermoregulation processes. Therefore, we performed a meta-analysis to investigate the influence of the climatic environment on leptin levels. A systematic literature search in the databases was carried out on 10 January 2020. Finally, 22 eligible articles were included in the current meta-analysis and a total of 13,320 participants were covered in the final analysis. It was shown that males of the “North” subgroup demonstrated significantly higher levels of leptin (10.02 ng/mL; CI: 7.92–12.13) than males of the “South” subgroup (4.9 ng/mL; CI: 3.71–6.25) (p = 0.0001). On the contrary, in the female group, a similar pattern was not detected (p = 0.91). Apparently, in order to maintain body temperature, higher leptin levels are required. The results of the study indicate that such effects are most pronounced in males and to a smaller extent in females, apparently due to a relatively high initial concentration of leptin in females. The correlation between leptin levels and climatic environment data support the hypothesis of leptin-mediated thermoregulation as an adaptive mechanism to cold climates.

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.

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.