ASSOCIATION OF POLYMORPHISM RS10841855 IN THE GLYCOGEN SYNTHASE 2 GENE WITH HYPERCHOLESTEROLEMIA AND TYPE 2 DIABETES MELLITUS Текст научной статьи по специальности «Фундаментальная медицина»

ГЕНЕТИКА GENETICS

И Hi^l

DOI: 10.18413/2658-6533-2022-8-2-0-1 УДК 575.174.015.3:616.379-008.64

Association of polymorphism rs10841855 in the glycogen synthase 2 gene with hypercholesterolemia and type 2 diabetes mellitus

Iuliia E. Azarova , Yaroslava E. Korvyakova , Alexander O. Glotov , Vladimir E. Ivakin , Alexey V. Polonikov

Kursk State Medical University, 3 Karl Marx St., Kursk, 305041, Russia Corresponding author: Iuliia E. Azarova (azzzzar@yandex.ru)

Abstract

Background: Glycogen synthase 2 (GYS2) catalyzes a key step in glycogenesis in the liver. Loss-of-function mutations in the GYS2 gene are associated with type 0a glycogenosis with characteristic fasting hypoglycemia and postprandial hyperglycemia. The association of the single nucleotide polymorphism rs10841855 (G>T) with the risk of type 2 diabetes mellitus (T2D) was first established in the European population, but data on the effect of this variant on the predisposition to T2D in the inhabitants of the Slavic ethnic group of Central Russia are not available in the literature. The aim of the study: To investigate the association of polymorphism rs10841855 (G>T) GYS2 with the risk of T2D in residents of Central Russia. Materials and methods: The study included 2668 unrelated individuals of Slavic origin, including 1387 patients with T2D and 1281 healthy volunteers. Geno-typing of the rs10841855 polymorphism of the GYS2 gene was performed by the MassArray Ana-lyzer-4 genetic analyzer. SNPStats software was used for statistical analysis of the data. Results: Linear regression established an association of the alternative rs10841855-T allele (OR 1.25; 95% CI 1.09-1.43; P=0.001) and the rs10841855-G/T genotype (OR 1.37; 95% CI 1, 13-1.66; P=0.0064) of GYS2 with an increased risk of T2D. The sex- and BMI-stratified analysis showed that the association of the rs10841855-G/T genotype (OR 1.71; 95% CI 1.32-2.22; P=0.0001) occurred only in females with BMI>25 kg/m2. In addition, carriage of the rs10841855-T allele was associated with a higher level of total cholesterol in the blood plasma of T2D patients (P=0.0068). According to the Roadmap Epigenomic Consortium, the rs10841855-T allele is associated with H3K4 histone methylation in the enhancer and promoter regions of the GYS2 gene in the liver. Conclusion: The association of rs10841855 of the GYS2 gene with an increased risk of T2D in residents of Central Russia was validated for the first time. The association may be explained by the low transcriptional activity of the glycogen synthase 2 gene in the carriers of the alternative allele of this polymorphism. Keywords: type 2 diabetes mellitus; single nucleotide polymorphism; GYS2; glycogen synthase; genetic predisposition

For citation: Azarova IE, Korvyakova YaE, Glotov AO, et al. Association of polymorphism rs10841855 in the glycogen synthase 2 gene with hypercholesterolemia and type 2 diabetes mellitus. Research Results in Biomedicine. 2022;8(2):138-147. DOI: 10.18413/2658-6533-2022-8-2-0-1

Introduction. According to the International Diabetes Federation, there are more than 537 million diabetic patients worldwide aged 20 to 79 years, with 90 percent having type 2 diabetes mellitus (T2D) [1]. The pathogenetic basis of the disease is dysfunction of pancreatic beta cells that secrete insulin, as well as insulin resistance of the liver, skeletal muscle, and adipose tissue, resulting in chronic hyperglyce-mia and disruption of all types of metabolism [2]. Type 2 diabetes complications include diabetic retinopathy, nephropathy, neuropathy, and macrovascular complications, which killed 6.7 million people in 2021 [1].

Type 2 diabetes mellitus is a classic multifactorial disease caused by genetic and environmental risk factors such as an excessive fat and refined carbohydrate diet [3, 4], sedentary lifestyle [5, 6], overweight and obesity [7, 8], chronic psycho-emotional stress [7, 8], and others. The genetic basis of the disease is being studied extensively: from 2007 to the present, 197 genome-wide association studies (GWAS) for T2D have been conducted worldwide, yielding a total of over four thousand single nucleotide markers (SNPs) [GWAS Catalog https://www.ebi.ac.uk/gwas/]. The majority of the detected loci are in non-coding regions of the genome, pseudogenes, and intergenic regions, which significantly complicates biological interpretation of the existing data layer and impedes clinical application of the results obtained. One of the polymorphic variants identified by R.A. Scott et al. in a genome-wide study in a European population of 159,208 people [9] became the SNP rs10841855 in the intron of the GYS2 gene, which encodes hepatic glycogen synthase 2. The enzyme catalyzes a critical step in hepatic glycogenesis, and its functional deficiency caused by GYS2 gene loss-of-function mutations is known as type 0a glycogenosis, which is characterized by fasting ketotic hypoglycemia and postprandial hyperglycemia due to the liver's inability to store glycogen [10, 11]. It should be noted that the liver plays a critical role in maintaining

normal glucose homeostasis by catalyzing the conversion of glucose into glycogen under hy-perglycemia, as well as the processes of glucose regeneration from non-carbohydrate precursors (gluconeogenesis) and glycogenolysis under glucose deficiency in the blood. Violations of the liver's glucostatic function can result in T2D: specifically, a decrease in the function of pancreatic beta-cell glucokinase leads to the development of maturity-onset diabetes of the young, MODY2 [12]. There is no data in the literature on the effect of the rs10841855 variant in the glycogen synthase 2 gene on predisposition to T2D in residents of the Slavic ethnic group of Central Russia, necessitating this validation study.

The aim of the study. To analyze the association of the polymorphic variant rs10841855 (G>T) of the GYS2 gene with the risk of T2D in residents of Central Russia.

Materials and methods. The KSMU Regional Ethics Committee approved the study protocol. The study included 1387 T2D patients (527 males and 860 females) with a mean age of 61,2±6,9 years who received treatment at the Endocrinology Department of the City Clinical Emergency Hospital, Kursk, from November 2016 to October 2019. Endo-crinologists confirmed the diagnosis of T2D using World Health Organization criteria [9, 10]. As previously described [13, 14, 15], 1281 healthy donors from the Regional Blood Transfusion Station (492 males and 789 females) with an average age of 60.7±5.6 years were invited to participate in the control group. On an empty stomach, all study participants completed a questionnaire to identify the main risk factors and gave 10 ml of venous blood. The phenol-chloroform method was used to extract genomic DNA from whole blood. The GYS2 gene was genotyped for rs10841855 (G>T) using iPLEX technology on a MassArray Analyzer 4 genomic mass spectrometer (Agena Bioscience). The concentrations of carbohydrate and lipid metabolism parameters (glucose, glycated hemoglobin, triglycerides,

total cholesterol, low- and high-density lipoproteins) were measured using a semi-automatic biochemical analyzer Clima MC-15 (RAL). The bioinformatics analysis was carried out using the online resources STRING (https://string-db.org) for visual assessment of the GYS2 interactome at the level of protein gene products; GTEx Portal (https://gtexpor-tal.org) to study tissue-specific expression of the GYS2 gene in carriers of variant rs10841855 alleles; and Haploreg v4.1 (https: //pubs.broadinstitute.org/mammals/) to study epigenetic regulation of the GYS2 gene. The atSNP resource (http://atsnp.biostat.wisc.edu) was used to investigate transcription factor affinity in the SNP region. The collected data were statistically analyzed with the SNPStats program (https://www.snpstats.net/start.htm). Co-dominant, dominant, recessive, over-dominant, and log-additive genetic models were all tested. As the best model, the one with the lowest numerical value of the Akaike criterion was

chosen. At P<0.05, the association was deemed significant. The Kolmogorov-Smirnov test was used to determine the normality of the distribution of biochemical parameters. Non-normally distributed indicators were labeled as median (Median), first (Q1), and third (Q3) quartiles: Median [Q1; Q3]. In such cases, the Kruskal-Wallis test was used to determine statistical significance. The observed differences between groups were considered statistically significant at P<0,05.

Results and Discussion. The studied SNP conformed to the Hardy-Weinberg equilibrium (P>0.05). According to data from the 1000 Genomes project deposited with Ensembl (https://www.ensembl.org), the frequency of the minor allele T of the studied SNP of the GYS2 gene in the Central Russian population was significantly lower than in the European population. Figure 1 shows a diagram of the rs10841855 genotype frequencies.

GV52rs10841855

I controls T2D patients

70,20%

63,10%

G/G

G/T

T/T

OR (95% CI) 1.00 1.37 (1.13-1.66) 1.14 (0.80-1.63) P=0.0064

Fig. 1. Genotype and allele frequencies for rs10841855 in the GYS2 gene in T2D patients

and controls

The genotype rs10841855-G/T GYS2 was associated with an increased risk of developing T2D (OR 1.37; 95% CI 1.13-1.66; P=0.0064). The minor allele rs10841855-T

(OR 1.25; 95 % CI 1.09-1.43; P=0.001) was significantly more common in the DM2 group (21.6%) than in the control group (18.1%).

Table 1

Analysis of frequencies of the GYS2 genotypes and alleles in T2D patients and controls stratified by gender and body mass index

GYS2 SNP Genotype Controls T2D patients OR (95% CI)* P* Controls T2D patients OR (95% CI)* P*

Males

BMI<25 kg/m2 BMI>25 kg/m2

G/G 59 (62.8) 72 (69.2) 1.00 267 (67.1) 273 (64.5) 1.00

G/T 28 (29.8) 23 (22.1) 0.63 (0.33-1.22) 0.36 113 (28.4) 124 (29.3) 1.03 (0.76-1.41) 0.69

T/T 7 (7.5) 9 (8.7) 1.08 (0.38-3.08) 18 (4.5) 26 (6.2) 1.31 (0.70-2.47)

rs10841855 T 22.3 19.7 0.85 (0.53-1.38) 0.60 18.7 20.8 1.14 (0.89-1.46) 0.29

G>T Females

BMI<25 kg/m2 BMI> kg/m2

G/G 107 (72.3) 50 (62.5) 1.00 466 (72.7) 480 (61.5) 1.00

G/T 28 (18.9) 21 (26.2) 1.68 (0.83-3.39) 0.31 132 (20.6) 258 (33.1) 1.71 (1.32-2.22) 0.0001

T/T 13 (8.8) 9 (11.2) 1.12 (0.42-3.00) 43 (6.7) 42 (5.4) 0.86 (0.54-1.37)

T 18.2 24.4 1.44 (0.91-2.30) 0.12 17.0 21.9 1.37 (1.13-1.66) 0.001

Note:* calculations are performed with adjustments for age (codominant model)

Table 2

Associations of the GYS2 genotypes with the total cholesterol plasma level of T2D patients

GYS2 genotype Total cholesterol concentration, Ме (Q1; Q3), mmol/L P

Entire group of T2D patients

G/G 5.02 (4.26; 6.1) 0.0068

G/T 5.36 (4.3; 6.27)

T/T 5.33 (4.3; 6.05)

Males with T2D

G/G 4.84 (4.05; 5.73) 0.16

G/T 5.2 (4.19; 6.22)

T/T 5.1 (4.09; 5.56)

Females with T2D

G/G 5.10 (4.37; 6.29) 0.11

G/T 5.40 (4.33; 6.29)

T/T 5.56 (4.36; 6.4)

Table 1 shows the results of the sex- and BMI-stratified analysis of the associations of the studied GYS2 polymorphism with the risk of T2D. The association of the genotype rs10841855-G/T (OR 1.71; 95% CI 1.32-2.22; P=0.0001) was found to be typical only for overweight and obese females. The described risk association was not observed in females with normal body weight or in males, regardless of BMI (P>0.05). Furthermore, when the relationship between GYS2 genotypes and biochemical parameters of T2D patients' plasma was examined, it was discovered that carrying the rs10841855-T allele was associated with a higher level of total cholesterol in the blood plasma of T2D patients (P=0.0068, Table 2).

The enzyme glycogen synthase 2 (2.4.1.11) was identified on the metabolic map of carbohydrate metabolism using the KEGG Pathways resource (Fig. 2). The enzyme controls the glycogenesis process, which begins in the liver with the phosphorylation of glucose by glucokinase (2.7.1.2). Phosphoglucomutase converts the reaction product, glucose-6-phos-phate, to glucose-1-phosphate (5.4.2.5). UDP-glucosyl pyrophosphorylase (2.4.1.29) uses the energy of one UTP molecule to convert glucose-1-phosphate to UDP-glucose. The enzyme glycogen synthase 2 (2.4.1.11) was identified on the metabolic map of carbohydrate metabolism using the KEGG Pathways resource (Fig. 2). The enzyme controls the gly-cogenesis process, which begins in the liver with the phosphorylation of glucose by gluco-kinase (2.7.1.2). Phosphoglucomutase converts the reaction product, glucose-6-phos-phate, to glucose-1-phosphate (5.4.2.5). UDP-glucosyl pyrophosphorylase (2.4.1.29) uses the energy of one UTP molecule to convert glucose-1-phosphate to UDP-glucose.

Glycogen synthase is involved in the process at 10-20 glucose residues in the primer chain, catalyzing chain elongation by forming alpha1^4-glycosidic bonds. A branching

enzyme (amylo-(1^4)-(1^6)-transglycosyl-ase) that forms alpha1^6 glycosidic bonds creates branch points in glycogen. Glycogen synthase activity is regulated both allosteri-cally and covalently [17]. Simultaneously, insulin is the primary activator of glycogenesis, promoting the dissociation of hepatic glycogen synthase from the intranuclear complex with the regulatory protein and its release into the cytosol for subsequent glucose phosphoryla-tion [18]. Glucose-6-phosphate, a reaction product, is an important allosteric activator of GYS2 [19]. Protein kinase A (PKA), protein kinase C (PKC), AMP kinase (AMPK), glycogen synthase kinase 3 (GSK3), and other kinases phosphorylate the GYS2 enzyme at seven amino acid residues. Insulin inhibits GSK3, which keeps glycogen synthase 2 active (dephosphorylated) [20]. Furthermore, insulin signaling activates PDE3B via PKB, preventing cAMP-mediated phosphorylation of GYS2 [18].

Functional partners of GYS2 (STRING data) form an interactome network of 10 proteins (Fig. 3), including glycogenin 1 (GYG1), glycogenin 2 (GYG2), glycogen-binding domain-containing protein 1 (STBD1), protein phosphatase 1 regulatory subunit 3C (PPP1R3C), regulatory subunit 3B protein phosphatase 1 (PPP1R3B), serine/threonine protein phosphatase 1 catalytic alpha subunit (PPP1CA), serine/threonine protein phosphatase 1 catalytic beta subunit (PPP1CB), ser-ine/threonine protein phosphatase 1 catalytic gamma subunit (PPP1CC), and laforine (EPM2A). According to the Gene Ontology analysis, all of the listed proteins are involved in glycogen metabolism (FDR=3.55*10-22). GYS2, GYG1, GYG2, EPM2A, PPP1R3C, and STBD1 are six network enzymes involved in glycogen biosynthesis (FDR=5.96*10-10), PPP1CA, PPP1CB, and PPP1CC are three proteins involved in the circadian regulation of gene expression (FDR=1.70*10-3).

Й ТА F. ЦЩ АI iï SUCROSE METABOLISM

У

AïtllIW 5'JK3C Uli

г. ос Loot de" гид itr aaotaboûsmi '

D-faœtosc

_. S until f -5Г

\32ï26\-О*—

j О-СНиэда-бГ

,__ 47. trace Euloi )

H2îJ2Li|-О

S

П+-

1.7 : .L

2 ?.. 4

D-Wienas

tïlJb

: J J .Jtl

. J J l45

Tl

L

___ 3-K5toeucraee

ihOlD: О

2-4.1 9

4

2.41.lOh

tew.

131 U

J 4.1.13

¡JCF'-gflKiSSfc

¿4 I J

► О

3*1

Lj ULnccee <

hfj.12®

I 1 5 I-~► О

- l t

^ïMctrat. Û I

-1 i.i.i VI—

L'-Fnrtort

-O*-

ut-1 3 7-L

l^-P-Ghab - —04—

-lïEAiasf

32121

'H :4"|.2пНэ*"

□аРФдае

fï-P-dipncii , j , j| I -iLÏDctiitn^

IUI Jî h

I ? г] ] ° I

! Щ

-Q

32.14

If»

_ 4 GDP-skcoiE _ ^J

1 I 4 |—-1"; J 1 O^-[Тп.зЦ-o"

'>] lufal t ûl

f? T I Plj-

I Î.IBUUI ïuçâïaud 1----~ f-—-1

L tiij^Hniidn cîijtjï 1 O** 2 T Г V---' CPP-gliLftli»-

fTTTwl '-'flltbiDef J " ■ 11 _fötWEidai

GÜfcbKectfF

0*-

: 33. U

sä 2

Ci

'jlyïuiïs:. : Lj.___J

ull70Hw0£CEft:L3 11

AEf-jiiKLre

2.1.7. il

JS.I3I

f 2.4,1.1J

I . I {uîrart lliUat: 12 71 Jni|-О

3 2 I Я

■ 1.; i '¿\

'Îlfilibff.sP

Ul.blî

H? 1>ЕЯй1№ ЬЕ-

г-. D-GîTirwtl

- J4.L É*«-г-

lj

2-1124.'h

D-auftKi-tfî

T

ijjii^ j

и i.i

32.1M|

îHJu l-^n^ P-E>-Gl4eüie-IF

Iteiloaljjps^i Ufluilwi-

111 y.: I

|5.-ij|Jljl

. J : 1Д

(.yira.iF

(Я—

Ï.4SSK

j.I i.G J : Jtfùl

l.'-Ji. àcïtjili T

3,2 1 SI

1

ЭД11

1 1 MljlQJÉ

ois-

CyïlcdMUiuO 13.Ü .1 -Qg I

. ¿Л l.ip]

î—i

541 ii

SbrJh, ■jlvra^tt

3-2.1 J

33.13

иа.Ш

Derliiii

Ш 1-33

3 2 12 3i I ':[I

huJtël

KlltQ5î

fafcMetllnlà^

Ц71.Г

7 41

3J.J i

: .149? ir|

Ij :i ml

»-[ГШ4?|-IH

С

Dtj.But о—I ? j i~Tjn—wb

^ I !fl|

n-№HtBB-|T

121 ; 41

1.171 IEB

2 M J

D-fta hifi-JÎT

I : l-l

ГТТзйТ

J.L-jl'MOl'f-l.rjFi

Fig. 2. Glycogen metabolism (KEGG pathways data https://www.genome.jp/kegg/pathway.html)

Fig. 3. GYS2 protein network

The current study successfully validated the association between the polymorphic locus rs10841855 in the regulatory region of the GYS2 gene and an increased risk of T2D in residents of Central Russia. The link between this locus and the disease was discovered in a genome-wide study on a European sample of 159208 people [9], where the rs10841855-G allele was linked to an increased risk of developing the disease. The alternative allele rs10841855-T was found to be risky in our study. The rs10841855-T allele has also been linked to a predisposition to T2D in the Japanese population [21]. The literature describes associations of GYS2 gene loss-of-function mutations with the development of glycoge-nosis type 0 [10] and type 0a [11], in which glycogen synthesis in the liver is disrupted, resulting in ketosis, fasting hypoglycemia with a decrease in lactate and alanine levels, and long-term persistent postprandial hyperglyce-mia. Increased glycogenesis in the liver of Ankita mice with the T2D phenotype, on the other hand, resulted in decreased hyperglycemia and normalization of liver energy homeostasis [22]. J. Chen et al. discovered a link between GYS2 gene polymorphism and blood plasma

insulin levels [23], and a study in the Korean population established an association between GYS2 and obesity, polycystic ovary syndrome, and gestational diabetes mellitus [24].

The current study demonstrated for the first time a link between carrying the rs10841855-T allele and a higher total cholesterol content in the blood plasma of T2D patients, which is explained by the outflow of glucose-6-phosphate accumulating in the liver into lipogenesis, which includes the synthesis of triacylglycerols and cholesterol synthesis.

The rs10841855-T allele is not associated with changes in gene expression in the liver, according to GTEx Portal experimental data. According to the Roadmap Epigenomic Consortium, the rs10841855-T allele is associated with four epigenetic marks in the liver: H3K4me1, H3K27ac in the enhancer region, and H3K4me3, H3K9ac in the GYS2 gene promoter, indicating that the studied variant has significant epigenetic regulatory potential. Interestingly, A. Sharifi-Zarchi discovered that weakly methylated DNA regions are high in H3K4me3, whereas highly methylated genome regions correspond to the H3K4me1 mark [25].

An analysis of transcription factor (TF) binding ability using the bioinformatics tool atSNP revealed that the minor allele rs10841855-T significantly reduces the affinity of three TFs: FOXO1 (P=3,76x10-3), FOXJ1 (P=4,60x10-3), and S0X10 (P=7,29*10-3) preventing them from binding to DNA in the region of that variant allele. It should be noted that FOXO1 in the liver acts as a gluconeogenesis activator [26] and is the primary target of insulin signaling [27].

Conclusion. Thus, the study discovered an association of the GYS2 gene rs10841855 with hypercholesterolemia and an increased risk of T2D in overweight and obese females, which may be due to epigenetic mechanisms of regulation of glycogen synthase 2 gene activity, as well as changes in the network of transcription factors that modulate glucose metabolism in carriers of the alternative allele rs10841855-T. More research, including phar-macogenetic approaches [28], will be needed to determine whether this polymorphism affects the efficacy of type 2 diabetes mellitus treatment.

Financial support

The study was supported by the Russian Science Foundation (N20-15-00227).

Conflict of interests

The authors have no conflict of interest to declare.

References

1. Sun H, Saeedi P, Karuranga S, et al. IDF Diabetes Atlas: Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Research and Clinical Practice. 2022;183:109119. DOI: https://doi.org/ 10.1016/j .diabres.2021.109119

2. DeFronzo RA. Pathogenesis of Type 2 Diabetes Mellitus. In: Bonora E, DeFronzo R, editors. Diabetes. Epidemiology, Genetics, Pathogenesis, Diagnosis, Prevention, and Treatment. Endocrinology. Springer, Cham; 2018:1-74. DOI: https://doi.org/10.1007/978-3-319-27317-4_8-1

3. Skovs0 S. Modeling type 2 diabetes in rats using high fat diet and streptozotocin. Journal

of Diabetes Investigation. 2014;5(4):349-358. DOI: https://doi.org/10.1111/jdi.12235

4. Malik VS, Hu FB. Sweeteners and risk of obesity and type 2 diabetes: the role of sugar-sweetened beverages. Current Diabetes Reports. 2012;12(2): 195-203. DOI: https://doi.org/10.1007/s11892-012-0259-6

5. Li DD, Yang Y, Gao ZY, et al. Sedentary lifestyle and body composition in type 2 diabetes. Diabetology and Metabolic Syndrome. 2022;14(1): 1-11. DOI: https://doi.org/10.1186/s13098-021-00778-6

6. Bertoglia MP, Gormaz JG, Libuy M, et al. The population impact of obesity, sedentary lifestyle, and tobacco and alcohol consumption on the prevalence of type 2 diabetes: Analysis of a health population survey in Chile, 2010. PLoS ONE. 2017;12(5):e0178092. DOI: https://doi.org/ 10.1371/j ournal.pone.0178092

7. Hackett RA, Steptoe A. Type 2 diabetes mellitus and psychological stress—a modifiable risk factor. Nature Reviews Endocrinology. 2017;13(9):547-560. DOI: https://doi.org/ 10.1038/nrendo.2017.64

8. World Health Organization. Global report on diabetes: executive summary (No. WHO/NMH/NVI/16.3). World Health Organization; 2016.

9. Scott RA, Scott LJ, Mägi R, et al. An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans. Diabetes. 2017;66(11):2888-2902. DOI: https://doi.org/10.2337/db16-1253

10. Orho M, Bosshard NU, Buist NRM, et al. Mutations in the liver glycogen synthase gene in children with hypoglycemia due to glycogen storage disease type 0. Journal of Clinical Investigation. 1998;102:507-515. DOI: https://doi.org/10.1172/JCI2890

11. Iijima H, Ago Y, Fujiki R, et al. Novel GYS2 mutations in a Japanese patient with glyco-gen storage disease type 0a. Molecular Genetics and Metabolism Reports. 2021;26:100702. DOI: https://doi.org/10.1016/jymgmr.2020.100702

12. Dedov II, Zubkova NA, Arbatskaya NYu, et al. MODY2: Clinical and molecular genetic characteristics of 13 cases of the disease. The first description of MODY in Russia. Problems of Endocrinology. 2009;55(3):3-7. Russian. DOI: https://doi.org/10.14341/probl20095533-7

13. Azarova I, Bushueva O, Konoplya A, et al. Glutathione S-transferase genes and the risk of

type 2 diabetes mellitus: Role of sexual dimorphism, gene-gene and gene-smoking interactions in disease susceptibility. Journal of Diabetes. 2018;10(5):398-407. DOI:

https://doi.org/10.1111/1753-0407.12623

14. Azarova IE, Klyosova EYu, Samgina TA, et al. Role of cyba gene polymorphisms in pathogenesis of type 2 diabetes mellitus. Medical Genetics. 2019;18(8):37-48. Russian. DOI: https://doi.org/10.25557/2073-7998.2019.08.37-48

15. Azarova IE, Klyosova EYu, Sakali SYu, et al. Contribution of rs11927381 polymorphism of the IGF2BP2 gene to the pathogenesis of type 2 diabetes. Research Results in Biomedicine. 2020;6(1):9-19. Russian. DOI: https://doi.org/10.18413/2658-6533-2020-6-1-0-2

16.Adeva-Andany MM, Gonzâlez-Lucân M, Donapetry-Garcia C, et al. Glycogen metabolism in humans. BBA Clinical. 2016;5:85-100. DOI: https://doi.org/10.1016/j .bbacli.2016.02.001

17.Zhang H, Ma J, Tang K, et al. Beyond energy storage: roles of glycogen metabolism in health and disease. FEBS Journal. 2021;288(12):3772-3783. DOI: https://doi.org/10.1111/febs.15648

18. Petersen MC, Vatner DF, Shulman GI. Regulation of hepatic glucose metabolism in health and disease. Nature Reviews Endocrinology. 2017;13(10):572-587. DOI: https://doi.org/10.1038/nrendo.2017.80

19.Pederson BA. Mammalian Glycogen Metabolism Enzymology, Regulation, and Animal Models of Dysregulated Glycogen Metabolism. Enzymology of Complex Alpha-Glu-cans. CRC Press; 2021:84-137.

20. Lewis GF, Carpentier AC, Pereira S, et al. Direct and indirect control of hepatic glucose production by insulin. Cell Metabolism. 2021;33(4):709-720. DOI: https://doi.org/10.1016/j.cmet.2021.03.007

21. Kanai M, Akiyama M, Takahashi A, et al. Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nature Genetics. 2018;50(3):390-400. DOI: https://doi.org/10.1038/s41588-018-0047-6

22.Lopez-Soldado I, Guinovart JJ, Duran J. Increasing hepatic glycogen moderates the diabetic phenotype in insulin-deficient Akita mice. Journal of Biological Chemistry. 2021;296:100498. DOI: https://doi.org/10.1016/jjbc.2021.100498

23. Chen J, Spracklen CN, Marenne G, et al. The trans-ancestral genomic architecture of glyce-mic traits. Nature Genetics. 2021;53(6):840-860. DOI: https://doi.org/10.1038/s41588-021-00852-9

24.Hwang JY, Lee EJ, Jin Go M, et al. Genome-wide association study identifies GYS2 as a novel genetic factor for polycystic ovary syndrome through obesity-related condition. Journal of Human Genetics. 2012;57(10):660-664. DOI: https://doi.org/10.1038/jhg.2012.92

25. Sharifi-Zarchi A, Gerovska D, Adachi K, et al. DNA methylation regulates discrimination of enhancers from promoters through a H3K4me1-H3K4me3 seesaw mechanism. BMC Genomics. 2017; 18(1): 1-21. DOI: https://doi.org/10.1186/s12864-017-4353-7

26. Oh KJ, Han HS, Kim MJ, et al. CREB and FoxO 1: two transcription factors for the regulation of hepatic gluconeogenesis. BMB Reports. 2013;46(12):567. DOI: https://doi.org/10.5483/BMBRep.2013.46.12.248

27.Puigserver P, Rhee J, Donovan J, et al. Insulin-regulated hepatic gluconeogenesis through FOXO1-PGC-1a interaction. Nature. 2003;423(6939):550-555. DOI: https://doi.org/ 10.1038/nature01667

28. Chernikova NA, Kamynina LL, Ametov AS, et al. The pharmacogenetics of hypoglycemia and the glycemic variability at the patients ith type 2 diabetes mellitus. Therapeutic archive. 2020;24;92(10):54-62. Russian. DOI: https://doi.org/10.26442/00403660.2020.10.000530

Received 20 January 2022 Revised 12 March 2022 Accepted 16 March 2022

Information about the authors Iuliia E. Azarova, Cand. Sci. (Medicine), Head of the Laboratory of Biochemical Genetics and Metabolomics, Research Institute for Genetic and Molecular Epidemiology, Associate Professor at the Department of Biological Chemistry, Kursk State Medical University, Kursk, Russia, E-mail: azzzzar@yandex.ru, ORCID: https://or-cid.org/0000-0001-8098-8052. Yaroslava E. Korvyakova, Student, Kursk State Medical University, Kursk, Russia, E-mail: yasya_kursk@mail.ru, ORCID: https://or-cid.org/0000-0003-4756-9707. Alexander O. Glotov, Student, Kursk State Medical University, Kursk, Russia, E-mail: aoglotoff@yandex.ru, ORCID: https://or-cid.org/0000-0002-4661-752X. Vladimir E. Ivakin, Cand. Sci. (Medicine), Associate Professor at the Department of Internal Diseases N1, Kursk State Medical University, Kursk,

Russia, E-mail: ivakinve@kursksmu.net, ORCID: https://orcid.org/0000-0003-1057-6292. Alexey V. Polonikov, Doct. Sci. (Medicine), Professor, Director of the Research Institute for Genetic and Molecular Epidemiology, Head of the Laboratory of Statistical Genetics and Bioinformatics, Research Institute for Genetic and

Molecular Epidemiology, Professor at the Department of Biology, Medical Genetics and Ecology, Kursk State Medical University, Kursk, Russia, Email: polonikov@rambler.ru, ORCID: https://orcid.org/0000-0001-6280-247X.