==== Front Vavilovskii Zhurnal Genet Selektsii Vavilovskii Zhurnal Genet Selektsii Vavilov Journal of Genetics and Breeding 2500-0462 2500-3259 The Federal Research Center Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences Russia Website: Vavilov Journal of Genetics and Breeding Original Article A rare splice site mutation in the gene encoding glucokinase/hexokinase 4 in a patient with MODY type 2 Редкий вариант мутации сайта сплайсинга гена, кодирующего глюкокиназу/гексокиназу 4, у пациента с MODY, подтип 2Ivanoshchuk D.E. Shakhtshneider E.V. Research Institute of Internal and Preventive Medicine – Branch of the Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Ovsyannikova A.K. Research Institute of Internal and Preventive Medicine – Branch of the Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Mikhailova S.V. Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Rymar O.D. Research Institute of Internal and Preventive Medicine – Branch of the Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Oblaukhova V.I. Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Yurchenko A.A. Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Voevoda M.I. Research Institute of Internal and Preventive Medicine – Branch of the Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia Correspondence to: D.E. Ivanoshchuk dinara2084@mail.ru 5 2020 24 3 299305 16 10 2019 6 12 2019 13 12 2019 Copyright © AUTHORS, 2020 2020 This work is licensed under a Creative Commons Attribution 4.0 License The article presents a variant of maturity onset diabetes of the young type 2, caused by a rare mutation in the GCK gene. Maturity onset diabetes of the young (MODY) is a hereditary form of diabetes with an autosomal dominant type of inheritance, an onset at a young age, and a primary defect in pancreatic β-cell function. This type of diabetes is different from classical types of diabetes mellitus (DM1 and DM2) in its clinical course, treatment strategies, and prognosis. Clinical manifestations of MODY are heterogeneous and may vary even among members of the same family, i. e., carriers of identical mutations. This phenotypic variation is due to the interaction of mutations with different genetic backgrounds and the influence of environmental factors (e. g., lifestyle). Using next-generation sequencing technology, the c.580–1G>A substitution (IVS5 –1G>A, rs1554335421) located in an acceptor splice site of intron 5 of the GCK gene was found in a proband. The identified variant cosegregated with a pathological phenotype in the examined family members. The GCK gene encodes glucokinase (hexokinase 4), which catalyzes the first step in a large number of glucose metabolic pathways such as glycolysis. Mutations in this gene are the cause of MODY2. The illness is characterized by an insignificant increase in the fasting glucose level, is a well-controlled disease without medication, and has a low prevalence of micro- and macrovascular complications of diabetes. The presented case of MODY2 reveals the clinical significance of a mutation in the splice site of the GCK gene. When nonclassical diabetes mellitus is being diagnosed in young people and pregnant women, genetic testing is needed to verify the diagnosis and to select the optimal treatment method. Key words: human; maturity onset diabetes of the young; MODY2; glucokinase gene; next-generation sequencing; genetic analysis; bioinformatics. В статье рассмотрен вариант развития моногенной формы сахарного диабета (MODY), обуслов- ленный редкой мутацией в гене GCK. Диабет MODY представляет собой сахарный диабет с аутосомно-доми- нантным типом наследования, возникающий в молодом возрасте и проявляющийся в дисфункции β-клеток поджелудочной железы. Этот тип отличается от классических типов сахарного диабета (СД1, СД2) клиниче- ским течением, тактикой лечения и прогнозом для пациента. Клинические проявления MODY гетерогенны и могут различаться даже у представителей одной семьи, носителей одинаковых мутаций. Это обусловлено как сочетанием мутаций в различных генах у индивидуума, так и воздействием внешних факторов. Методом секвенирования нового поколения у пробанда была идентифицирована замена c.580 –1G>A (IVS5 –1G>A, rs1554335421), локализующаяся в акцепторном сайте сплайсинга пятого интрона гена GCK. Обнаруженный вариант сегрегировал с патологическим фенотипом у обследованных членов семьи. Ген GCK кодирует глюко- киназу (гексокиназу 4), которая катализирует первый шаг в различных путях метаболизма глюкозы. Мутации в этом гене ассоциированы с сахарным диабетом взрослого типа у молодых, подтип 2 (MODY2). Заболевание характеризуется незначительным повышением глюкозы натощак, хорошо контролируется медикамента- ми и отличается низкой распространенностью микро- и макрососудистых осложнений. Представленный в исследовании случай MODY2 выявил клиническую значимость мутации в сайте сплайсинга гена GCK. При возникновении у молодых людей и беременных женщин неклассического сахарного диабета проведение генетического тестирования необходимо для подтверждения диагноза и оптимального выбора тактики и способа лечения. human maturity onset diabetes of the young MODY2 glucokinase gene next-generation sequencing genetic analysis bioinformatics. человек диабет взрослого типа у молодых MODY2 ген глюкокиназы секвенирование нового поколения генетический анализ биоинформатика ==== Body Introduction Maturity onset diabetes of the young (MODY) is a hereditary form of diabetes with autosomal dominant inheritance and is characterized by onset at a young age and by the presence of an initial defect in pancreatic β-cell function. This type of diabetes differs from classic types of diabetes mellitus-type 1 (DM1) and type 2 (DM2) in disease progression, in treatment strategies, and prognosis (Anık et al., 2015). Up to 80 % of MODY cases are not detected or are misdiagnosed as DM1 or DM2; therefore, patients with an incorrectly diagnosed type of diabetes are often prescribed inadequate therapy (Shields et al., 2010). On average, MODY is detected in 2–5 % of cases of diabetes (the rest being mostly DM1 and DM2) (Fajans et al., 2001). To reliably diagnose MODY in a patient, molecular genetic analysis should be carried out. To date, 14 types of MODY (MODY1 through MODY14) have been identified, each associated with mutations in a specific gene: HNF4A, GCK, HNF1A, PDX1, HNF1B, NEUROD1, KLF11, CEL, PAX4, INS, BLK, KCNJ11, ABCC8 and APPL1 (Thanabalasingham et al., 2011; Bonnefond et al., 2012; McDonald et al., 2013; Lachance, 2016; Ovsyannikova et al., 2016). Fourteen MODY-associated genes explain 70–85 % of the disease cases and are involved in various stages of glucose metabolism regulation (Thanabalasingham et al., 2011; Bonnefond et al., 2012; Lachance, 2016). According to various researchers, 11 to 30 % of MODY cases are caused by mutation in other genes (Edghill et al., 2010; Bonnefond et al., 2012). These forms of MODY are commonly referred to as MODY-X. Because the vast majority of pathogenic mutations are found in exons and adjacent splicing sites of genes (Stenson et al., 2017), it is reasonable to perform whole-exome sequencing on genomic DNA from individuals with MODY, with subsequent genetic testing of their relatives for the identified mutation. MODY verification allows for successful patient management and ensures healthy pregnancy and provision of genetic counseling to families (Lachance, 2016). Examination of relatives of MODY probands makes it possible to diagnose hyperglycemia in the preclinical phase. In this report, we describe a clinical case of a family with MODY2 associated with a rare splice site mutation in the glucokinase (GCK ) gene identified by the next-generation sequencing technology. Materials and methods The study protocol was approved by the local Ethics Committee of the Institute of Internal and Preventive Medicine (Branch of the Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia, approval # 22.06.2008). Written informed consent to be examined and to participate in the study was obtained from each patient. For individuals younger than 18 years, the informed consent form was signed by a parent or legal guardian. Blood samples were collected from the ulnar vein for biochemical analysis in the morning on an empty stomach. Lipid levels (cholesterol, triglycerides, low-density lipoprotein cholesterol, and high-density lipoprotein cholesterol) and glucose concentration were determined on a KoneLab 300i biochemical analyzer (Thermo Fisher Scientific, Waltham, MA, USA) with Thermo Fisher Scientific reagents. Genomic DNA was isolated from leukocytes of venous blood by phenol-chloroform extraction (Sambrook, Russell, 2006). Quality of the extracted DNA was assessed on a capillary electrophoresis system, Agilent 2100 Bioanalyzer (Agilent Technologies Inc., USA). Sequencing of patients’ DNA was carried out on an Illumina HiSeq1500 instrument (Illumina, San Diego, CA, USA). The enrichment and library preparation were performed with the SureSelectXT Human All Exon V5 + UTRs Kit (Agilent Technologies Inc., USA). Reads were mapped to the reference human genome (GRCh37) by means of the Burrow–Wheeler Alignment tool (BWA v.0.7.12) (Li, Durbin, 2009). Polymerase chain reaction (PCR)-generated duplicates were removed in the PICARD software (https://broadinstitute.github.io/picard/). A search for single-nucleotide variants (SNVs) was conducted using the Genome Analysis Toolkit v.3.3 package by the procedure for local remapping of short insertions/deletions and recalibration of read quality (McKenna et al., 2010). The depth of coverage was 34× to 53×. SNVs with genotype quality scores < 20 and coverage depth < 10× were filtered out and excluded from further analysis. Annotation of the SNVs was performed in the ANNOVAR software (Wang et al., 2010) using the 1000 Genomes Project (The 1000 Genomes Project Consortium…, 2015) and The Genome Aggregation Database (gnomAD) (Karczewski et al., 2019) databases. We selected the spectrum of rare and novel sequence variants in MODY genes (HNF4A, GCK, HNF1A, PDX1, HNF1B, NEUROD1, KLF11, CEL, PAX4, INS, BLK, KCNJ11, ABCC8 and APPL1). Rare variants were selected if their minor allele frequency (MAF) was ≤ 0.5 in the 1000 Genomes Project and gnomAD. Heterozygous substitution c.580–1G>A (IVS5 –1G>A) at an acceptor splice site of intron 5 of the GCK gene was found in the proband and her sister. To predict the possible effect of the SNV on splicing regulation, we employed the SPANR software (Xiong et al., 2015). The substitution was corroborated by Sanger sequencing of the DNA fragment containing exons 5 and 6, intron 5, and parts of introns 4 and 6 using the following forward and reverse primers: 5′-CAGGGAGCCTCAGCAGTCTGGA-3′ and 5′-GCCACGAGGCCTATCTCTCCCC-3′. The oligonucleotides were designed in the Primer-Blast software (https:// www.ncbi.nlm.nih.gov/tools/primer-blast/) and were synthesized by the Biosset company (Russia, Novosibirsk). The sequencing reactions were carried out on an automated ABI 3500 DNA sequencer (Thermo Fisher Scientific, USA) with the BigDye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific, USA). PCR was set up using BioMaster LR HS-PCR (2×) (BioLabMix, Russia), 1 μL of each primer, and 1 μL of DNA, with a total final volume of 25 μL. The thermocycling program consisted of initial denaturation at 94 °C for 3 min and then 35 cycles at 94 °C for 30 s, 66 °C for 30 s, and 72 °C for 50 s. The PCR products were evaluated by electrophoresis in a 5 % polyacrylamide gel after visualization with an ethidium bromide solution. A 100 bp DNA Ladder (BioLabMix) was simultaneously run on the gel as molecular size markers. The amplicons were purified using Agencourt AMPure Xp beads (Beckman Coulter, USA). The sequencing reactions were conducted on an automated ABI 3500 DNA sequencer (Thermo Fisher Scientific, USA) via the BigDye Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific). The sequences were analyzed in the Vector NTI® Advance software (Thermo Fisher Scientific). The hg19 version of the human genome served as a reference sequence for the alignment. Results The white European 44-year-old female proband was under medical observation. When she underwent routine screening in 2012 (at age 40), hyperglycemia at 7.2 mmol/L was revealed. No complaints were registered. During subsequent glycemia control, maximal fasting glucose was 7.2 mmol/L, and the postprandial one was 8.9 mmol/L. C-peptide was 1.83 ng/mL (reference range 0.5–3.2 ng/mL), immunoreactive insulin was 7.9 μU/mL (reference range 2.0–25.0 μU/ mL), and glycated hemoglobin (HbA1c) was 7.1 %. Antibodies to insulin, to pancreatic islet cells, and to glutamic acid decarboxylase were absent. Blood biochemical analysis and determination of thyroid status did not reveal any abnormalities. Ultrasonography of internal organs, echocardiography, and a study of brachiocephalic vessels did not uncover any pathology. The body mass index (BMI) was 20.2 kg/m2. DM2 was diagnosed in the patient, and sitagliptin was prescribed. At the age of 26, the patient spontaneously delivered a healthy girl at 39 weeks of gestation; hyperglycemia was not detected during the pregnancy. The sister of the proband is a white European 35-yearold woman. At age 23, during tests before mastectomy for mastopathy, she got a diagnosis of fasting hyperglycemia (6.3 mmol/L). The patient did not have any complains, and a proper diet was recommended. At the age of 29, during additional examination before cholecystectomy for cholelithiasis, she received a diagnosis of DM2, and vildagliptin was prescribed at the dose of 50 mg twice a day. At age 31, the proband’s sister visited an endocrinologist at the outpatient clinic of the Institute of Internal and Preventive Medicine (Branch of the Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia) with complaints of a failure to get pregnant within a year. On examination, BMI was 20.6 kg/m2, and the objective status was unremarkable. Blood biochemical analysis revealed hypercalcemia (3.25 mmol/L), increased levels of high-density lipoprotein cholesterol (85 mg/dL), hypercholesterolemia (220 mg/ dL), and hyperglycemia (6.8 mmol/L), but other analyzed parameters were within reference ranges. The HbA1c level was 7.1 %. Antibodies to insulin, to pancreatic islet cells, and to glutamic acid decarboxylase were absent. Thyroid-stimulating hormone concentration was 0.759 mU/ mL (reference range 0.4–4.0), whereas the prolactin level was 216 ng/ mL (reference range 1.2–19.5). Echocardiography, Doppler sonography of extracranial parts of cerebral vessels, and abdominal and renal ultrasonographic examination revealed no pathology. Cysts were found in both thyroid lobes during the ultrasonography. Given the existence of the proband’s relatives with impaired glucose metabolism, persistence of normal C-peptide levels, the absence of diabetes-associated autoantibodies, normal BMIs of the proband and her sister, and stable mild hyperglycemia, MODY was assumed. Exons and adjacent splice sites of MODY-associated genes were analyzed by whole-exome sequencing in the proband and her sister. As a result, heterozygous substitution c.580–1G>A (IVS5 –1G>A) at an acceptor splice site of intron 5 of the GCK gene was found in the proband and her sister. The IVS5 (–1G>A) polymorphism of GCK was submitted in ClinVar with an accession number of rs1554335421 (Landrum et al., 2018), but was absent in the 1000 Genomes Project (The 1000 Genomes Project Consortium…, 2015), in gnomAD project databases at the moment of publication. Subsequent genetic analysis by Sanger sequencing of the family members (mother, father, daughter, and nephew of the proband) uncovered segregation of the substitution with DM as an autosomal dominant trait (see the Figure). Fig. 1. Mutation identified in GCK gene. a – family history with inherited diabetes mellitus (DM). Asterisk indicates medically examined family members; b – schematic representation of the mutation in the splicing acceptor site and chromatogram of DNA sequence with mutated allele c.580 –1G>A (IVS5 –1G>A, rs1554335421) of the GCK gene. Our results, literature data, and databases suggest that this splice site mutation is likely pathogenic. After confirmation of GCK-MODY in the proband and her sister, their relatives were screened for carbohydrate metabo lism disorders. The proband’s mother and daughter did not have any abnormalities. The proband’s father showed impaired fasting glucose. No complaints were registered, venous plasma fasting glucose was 6.3 mmol/L, and 2 h after the oral glucose tolerance test, it was 7.5 mmol/L. At present, the man does not take any medication. The same heterozygous substitution rs1554335421 (IVS5 –1G>A) in the proband’s father’s GCK gene was detected by genetic testing. The proband’s sister had her first pregnancy in 2014 (at age 31). In 2015, a boy weighing 3640 g was born by a caesarean section at 39 weeks of gestation. The pregnancy was complicated: premature rupture of membranes, weakness of labor, and fetal hypoxia. After delivery, due to stable glycemic indexes, it was decided that insulin therapy should be discontinued. In January 2018, during treatment with diet, the patient’s HbA1c was 6.4 %. The neonatal period of the proband’s nephew was unremarkable. In 2017, his blood biochemical analysis resulted in a diagnosis of hyperglycemia (6.9 mmol/L). HbA1c was 6.3 %, and the C-peptide level was 0.54 ng/mL. Antibodies to insulin, pancreatic β-cells, and glutamic acid decarboxylase were undetectable. The same heterozygous substitution rs1554335421 (IVS5; –1G>A) in the GCK gene was identified by genetic testing. At present, the child is under medical observation at Almazov Federal Medical Research Centre (Saint Petersburg, Russia); because of GCK-MODY, a balanced diet was recommended. Discussion It is known that in young patients with impaired carbohydrate metabolism, DM1, DM2, or rarer monogenic forms of diabetes may be diagnosed. At the onset of the disease, the proband and her sister had no symptoms characteristic for the common types of diabetes, fasting hyperglycemia was not progressing, and carbohydrate metabolism disorders were detected during routine screening. The presence of DM in the proband’s sister, persistence of normal C-peptide levels, a lack of autoantibodies, and a normal BMI in the proband and her sister pointed to MODY (Chakera et al., 2015). Heterozygous splice site mutation c.580–1G>A (rs1554335421) in intron 5 of their GCK gene was identified by genetic testing. Mutations in this gene are associated with DM2, MODY, and neonatal DM (Plengvidhya et al., 2009; Lachance, 2016). More than 600 variants of the GCK gene associated with MODY have been described, and the list of the mutations is constantly growing. The vast majority of the mutations are missense substitutions, but splice site mutations, deletions, and insertions are reported too (Stenson et al., 2017). The GCK gene is located in chromosomal region 7p15.3- p15.1 and consists of 12 exons that encode a 465-amino-acid protein, glucokinase (Osbak et al., 2009), which is one of four members of the hexokinase family of enzymes. In 1992, GCK was the first gene to be linked to MODY. It plays an important regulatory role in glucose metabolism. Glucokinase catalyzes phosphorylation of glucose to produce glucose-6-phosphate as the first step of glycolysis in pancreatic β-cells (Matschinsky et al., 1993; Iynedjian, 2009). Most individuals with heterozygous GCK mutations show fasting plasma glucose levels between 5.5 and 8.0 mmol/L and a small increase in plasma glucose (< 3 mmol/L in 70 % of the patients) 2 h after the oral glucose test (Stride et al., 2002). This feature also explains asymptomatic fasting hyperglycemia (HbA1c range 5.8–7.6 % (40–60 mmol/mol)) and rare microvascular and macrovascular complications in patients with GCK-MODY (Caetano et al., 2012; Steele et al., 2014). Most patients have an aberrant fasting glucose level or impaired glucose tolerance, and less than 50 % of the affected individuals have diabetes, which is diagnosed during childhood, adolescence, or pregnancy (Caetano et al., 2012). In a study on Italian patients under 18 years of age with incidental hyperglycemia, it was estimated that 15 % of these cases are caused by GCK mutations (Lorini et al., 2009). It was found here that the proband and her sister carry a heterozygous substitution, c.580–1G>A (IVS5 –1G>A, rs1554335421), at an acceptor splice site of GCK intron 5. This allelic variant is of interest because consensus donor (GT dinucleotide) and acceptor (AG dinucleotide) splice sites are highly conserved. Point mutations at these loci can lead to cryptic splice site activation and synthesis of aberrant protein isoforms. In silico analysis of the functional significance of this substitution suggested that the inclusion of exons 5, 6, and 7 in gene transcripts will be reduced in case of the detected variant (see the Table). Table 1. The predicted (by SPANR) effect of c.580 –1G>A (IVS5 –1G>A, rs1554335421) of the GCK gene on splicing Note. PSI – percentage of transcripts with the exon spliced in: * dPSI: a maximal difference in PSI across 16 tissues; ** dPSI_percentile: percentile of mutant dPSI among dPSIs of common SNPs; *** PSI_WT: predicted PSI in the wild type. DM and a family history of gestational diabetes. Experiments on lymphoblastoid cells indicate that rs1554335421 (IVS5 –1G>A) of GCK can activate a cryptic splice site in intron 5 and cause retention of 27 bp of the intron (Toaima et al., 2005). Information about rs1554335421 (IVS5 –1G>A) of GCK is absent in the 1000 Genomes Project, The Exome Aggregation Consortium, and GNOMAD (https://gnomad. broadinstitute.org/) databases; however, taking into account the previous study and the data we obtained here, carriage of the A allele at position –1 of intron 5 is most likely a causative dominant variant of the GCK gene in MODY-affected people. Cryptic splice site activation and formation of several alternative transcripts with intron 7 fragments’ retention were demonstrated in a system of model GCK minigenes with acceptor site mutation IVS7 (–1G>C) (Igudin et al., 2014). Model mice with the homozygous mutation in the splice site of β-cell-specific exon 1 IVS1A (–1G>T) show hyperglycemia, glucosuria, and growth retardation and die within the first week after birth. This phenotype can be explained by exon skipping or intron retention (Inoue et al., 2004). Splicing sites affected by mutations have been described for many pathological phenotypes: neurofibromatosis type 1 (Jang et al., 2016), familial hypercholesterolemia (Shakhtshneider et al., 2017), Wiskott–Aldrich syndrome and chronic colitis (Esmaeilzadeh et al., 2018), hypophosphatemic rickets (Ma et al., 2015), and others. Mutations affecting splicing have been found not only in canonical splicing sites but also in introns and exons and may have a tissue-specific effect, as in familial dysautonomia (Slaugenhaupt et al., 2001; Abramowicz, Gos, 2018). That analysis indicated that the donor splice site mutations were more prevalent than the acceptor splice site variants (ratio 1.5:1.0) (Abramowicz, Gos, 2018). Because the mutations in the GCK gene can cause a mild clinical phenotype, which can vary under the influence of many genetic and lifestyle factors, research on the carriers of these mutations is essential for identifying additional risk factors. GCK-MODY is inherited as an autosomal dominant trait manifested throughout the lifespan as stable, mild fasting hyperglycemia usually reaching 6.7 mmol/L and higher only in middle age (Wedrychowicz et al., 2017). A similar pattern was observed here in the proband and her sister. Nonetheless, the metabolic disturbances in the carriers of GCK mutations are present from birth and can be identified already in the first years of life, almost all of them after puberty (Steele et al., 2014). The proband’s nephew, who is a heterozygous mutation carrier, developed carbohydrate metabolism disorders at two years of age. It has been reported that carriers of GCK gene mutations with a long history of hyperglycemia (48.6 years on average) usually have micro- and macrovascular complications of diabetes and are at a risk of cardiovascular diseases that is identical to that in the general population (Pruhova et al., 2013). Patients with GCK-MODY in childhood and adolescence can be treated only with diet in most cases, and glucoselowering therapy should be considered during pregnancy (Lachance, 2016). The present subtype of MODY2 (in the proband and her sister) currently is treated with diet resulting in sufficient glycemic control. The presence or absence of a GCK mutation in the fetus affects its sensitivity to maternal hyperglycemia (Chakera et al., 2012). If the fetus does not have the mutation, then it will secrete insulin excessively and as a result have a risk of macrosomia (Spyer et al., 2009). In that case, low doses of insulin should be prescribed during pregnancy (Chakera et al., 2014). During her pregnancy, the proband’s sister was given insulin injections in small doses. After delivery, insulin therapy was discontinued. Subsequently, the child was found to carry the same substitution. Conclusion The presented subtype of MODY2 reveals the clinical significance of the mutation in a splice site of the GCK gene. When nonclassical diabetes mellitus is being diagnosed in young people and pregnant women, genetic testing is needed to verify the diagnosis and to select the optimal treatment method. Conflict of interest The authors declare no conflict of interest. References Abramowicz A., Gos M. Splicing mutations in human genetic disorders: examples, detection, and confirmation. J. Appl. Genet. 2018; 59:253-268. DOI 10.1007/s13353-018-0444-7. Anık A., Çatlı G., Abacı A., Böber E. Maturity-onset diabetes of the young (MODY): an update. J. Pediatr. Endocrinol. Metab. 2015;28: 251-263. DOI 10.1515/jpem-2014-0384. Bonnefond A., Philippe J., Durand E., Dechaume A., Huyvaert M., Montagne L., Marre M., Balkau B., Fajardy I., Vambergue A., Vatin V., Delplanque J., Le Guilcher D., De Graeve F., Lecoeur C., Sand O., Vaxillaire M., Froguel P. Whole-exome sequencing and high throughput genotyping identified KCNJ11 as the thirteenth MODY gene. PLoS One. 2012;7(6):e37423. DOI 10.1371/journal. pone.0037423. Caetano L.A., Jorge A.A., Malaquias A.C., Trarbach E.B., Queiroz M.S., Nery M., Teles M.G. Incidental mild hyperglycemia in children: two MODY 2 families identified in Brazilian subjects. Arq. Bras. Endocrinol. Metabol. 2012;56(8):519-24. http://dx.doi.org/10.1590/ S0004-27302012000800010. Chakera A.J., Carleton V.L., Ellard S., Wong J., Yue D.K., Pinner J., Hattersley A.T., Ross G.P. Antenatal diagnosis of fetal genotype determines if maternal hyperglycemia due to a glucokinase mutation requires treatment. Diabetes Care. 2012;35(9):1832-1834. DOI 10.2337/dc12-0151. Chakera A.J., Spyer G., Vincent N., Ellard S., Hattersley A.T., Dunne F.P. The 0.1 % of the population with glucokinase monogenic diabetes can be recognized by clinical characteristics in pregnancy: the Atlantic Diabetes in Pregnancy cohort. Diabetes Care. 2014; 37(5):1230-1236. DOI 10.2337/dc13-2248. Chakera A.J., Steele A.M., Gloyn A.L., Shepherd M.H., Shields B., Ellard S., Hattersley A.T. Recognition and management of individuals with hyperglycemia because of heterozygous glucokinase mutations. Diabetes Care. 2015;38(7):1383-1392. DOI 10.2337/dc14- 2769. Edghill E.L., Minton J.A., Groves C.J., Flanagan S.E., Patch A.M., Rubio- Cabezas O., Shepherd M., Lenzen S., McCarthy M.I., Ellard S., Hattersley A.T. Sequencing of candidate genes selected by beta cell experts in monogenic diabetes of unknown aetiology. J. Pancreas – JOP. 2010;11(1):14-17. http://dx.doi.org/10.6092/1590- 8577/3864. Esmaeilzadeh H., Bordbar M.R., Dastsooz H., Silawi M., Fard M.A.F., Adib A., Kafashan A., Tabatabaei Z., Sadeghipour F., Faghihi M.A. A novel splice site mutation in WAS gene in patient with Wiskott-Aldrich syndrome and chronic colitis: a case report. BMC Med. Genet. 2018;19(1):123. DOI 10.1186/s12881-018-0647-0. Fajans S.S., Bell G.I., Polonsky K.S. Molecular mechanisms and clinical pathophysiology of maturity-onset diabetes of the young. N. Engl. J. Med. 2001;345(13):971-980. DOI 10.1056/NEJMra002168. Igudin E.L., Spirin P.V., Prassolov V.S., Rubtsov P.M., Zubkova N.A., Tyul’pakov A.N., Petryaikina E.E. Functional characterization of two novel splicing mutations of glucokinase gene associated with maturity-onset diabetes of the young type 2 (MODY2). Mol. Biol. 2014;48(2):248-253. DOI 10.1134/S0026893314020071. Inoue M., Sakuraba Y., Motegi H., Kubota N., Toki H., Matsui J., Toyoda Y., Miwa I., Terauchi Y., Kadowaki T., Shigeyama Y., Kasuga M., Adachi T., Fujimoto N., Matsumoto R., Tsuchihashi K., Kagami T., Inoue A., Kaneda H., Ishijima J., Masuya H., Suzuki T., Wakana S., Gondo Y., Minowa O., Shiroishi T., Noda T. A series of maturity onset diabetes of the young, type 2 (MODY2) mouse models generated by a large-scale ENU mutagenesis program. Hum. Mol. Genet. 2004;13(11):1147-1157. DOI 10.1093/hmg/ddh133. Iynedjian P.B. Molecular physiology of mammalian glucokinase. Cell. Mol. Life Sci. 2009;66(1):27-42. DOI 10.1007/s00018-008-8322-9. Jang M.A., Kim Y.E., Kim S.K., Lee M.K., Kim J.W., Ki C.S. Identification and characterization of NF1 splicing mutations in Korean patients with neurofibromatosis type 1. J. Hum. Genet. 2016;61: 705-709. DOI 10.1038/jhg.2016.33. Karczewski K.J., Francioli L.C., Tiao G., Cummings B.B., Alföldi J., Wang Q., Collins R.L., Laricchia K.M., Ganna A., Birnbaum D.P., Gauthier L.D., Brand H., Solomonson M., Watts N.A., Rhodes D., Singer-Berk M., England E.M., Seaby E.G., Kosmicki J.A., Walters R.K., Tashman K., Farjoun Y., Banks E., Poterba T., Wang A., Seed C., Whiffin N., Chong J.X., Samocha K.E., Pierce-Hoffman E., Zappala Z., O’Donnell-Luria A.H., Minikel E.V., Weisburd B., Lek M., Ware J.S., Vittal C., Armean I.M., Bergelson L., Cibulskis K., Connolly K.M., Covarrubias M., Donnelly S., Ferriera S., Gabriel S., Gentry J., Gupta N., Jeandet T., Kaplan D., Llanwarne C., Munshi R., Novod S., Petrillo N., Roazen D., Ruano-Rubio V., Saltzman A., Schleicher M., Soto J., Tibbetts K., Tolonen C., Wade G., Talkowski M.E.; The Genome Aggregation Database Consortium, Neale B.M., Daly M.J., MacArthur D.G. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. bioRxiv. 2019;531210. https://doi.org/10.1101/531210. Lachance C.H. Practical aspects of monogenic diabetes: A clinical point of view. Can. J. Diabetes. 2016;40(5):368-375. DOI 10.1016/j. jcjd.2015.11.004. Landrum M.J., Lee J.M., Benson M., Brown G.R., Chao C., Chitipiralla S., Gu B., Hart J., Hoffman D., Jang W., Karapetyan K., Katz K., Liu C., Maddipatla Z., Malheiro A., McDaniel K., Ovetsky M., Riley G., Zhou G., Holmes J.B., Kattman B.L., Maglott D.R. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018;46(D1):D1062-D1067. DOI 10.1093/nar/ gkx1153. Li H., Durbin R. Fast and accurate short read alignment with Burrows– Wheeler transform. Bioinformatics. 2009;25(14):1754-1760. DOI 10.1093/bioinformatics/btp324. Lorini R., Klersy C., d’Annunzio G., Massa O., Minuto N., Iafusco D., Bellannè-Chantelot C., Frongia A.P., Toni S., Meschi F., Cerutti F., Barbetti F.; Italian Society of Pediatric Endocrinology and Diabetology (ISPED) Study Group. Maturity-onset diabetes of the young in children with incidental hyperglycemia : amulticenter Italian study of 172 families. Diabetes Care. 2009;32(10):1864-1866. DOI 10.2337/dc08-2018. Ma S.L., Vega-Warner V., Gillies C., Sampson M.G., Kher V., Sethi S.K., Otto E.A. Whole exome sequencing reveals novel PHEX splice site mutations in patients with Hypophosphatemic rickets. PLoS One. 2015;10(6):e0130729. DOI 10.1371/journal.pone.0130729. Matschinsky F., Liang Y., Kesavan P., Wang L., Froguel P., Velho G., Cohen D., Permutt M.A., Tanizawa Y., Jetton T.L. Glucokinase as pancreatic cell glucose sensor and diabetes gene. J. Clin. Invest. 1993;92(5):2092-2098. DOI 10.1172/JCI116809. McDonald T.J., Ellard S. Maturity onset diabetes of the young: Identification and diagnosis. Ann. Clin. Biochem. 2013;50:403-415. DOI 10.1177/0004563213483458. McKenna A., Hanna M., Banks E., Sivachenko A., Cibulskis K., Kernytsky A., Garimella K., Altshuler D., Gabriel S., Daly M., DePristo M.A. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297-1303. DOI 10.1101/gr.107524.110. Osbak K.K., Colclough K., Saint-Martin C., Beer N.L., Bellanné-Chantelot C., Ellard S., Gloyn A.L. Update on mutations in glucokinase (GCK), which cause maturity-onset diabetes of the young, permanent neonatal diabetes, and hyperinsulinemic hypoglycemia. Hum. Mutat. 2009;30:1512-1526. DOI 10.1002/humu.21110. Ovsyannikova A.K., Rymar O.D., Shakhtshneider E.V., Klimontov V.V., Koroleva E.A., Myakina N.E., Voevoda M.I. ABCC8-related maturity-onset diabetes of the young (MODY12): clinical features and treatment perspective. Diabetes Ther. 2016;7(3):591-600. DOI 10.1007/s13300-016-0192-9. Plengvidhya N., Boonyasrisawat W., Chongjaroen N., Jungtrakoon P., Sriussadaporn S., Vannaseang S., Banchuin N., Yenchitsomanus P.T. Mutations of maturity-onset diabetes of the young (MODY) genes in Thais with early-onset type 2 diabetes mellitus. Clin. Endocrinol. (Oxf.). 2009;70(6):847-853. DOI 10.1111/j.1365-2265.2008. 03397.x. Pruhova S., Dusatkova P., Kraml P.J., Kulich M., Prochazkova Z., Broz J., Zikmund J., Cinek O., Andel M., Pedersen O., Hansen T., Leb J. Chronic mild hyperglycemia in GCK-MODY patients does not increase carotid intima-media thickness. Int. J. Endocrinol. 2013;2013:718254. DOI 10.1155/2013/718254. Sambrook J., Russell D.W. Purification of nucleic acids by extraction with phenol:chloroform. Cold Spring Harbor Protoc. 2006;2006(1): 4455. DOI 10.1101/pdb.prot4455. Shakhtshneider E.V., Ivanoshchuk D.Е., Makarenkova К.V., Orlov P.S., Timoshchenko О.V., Bazhan S.S., Nikitin Yu.P., Voevoda М.I. Cascade genetic screening in diagnostics of heterozygous familial hypercholesterolemia: clinical case. Russ. J. Cardiol. 2017;6(1460): 178-179. http://dx.doi.org/10.15829/1560-4071-2017-6-178-179. (in Russian) Shields B.M., Hicks S., Shepherd M.H., Colclough K., Hattersley A.T., Ellard S. Maturity-onset diabetes of the young (MODY): how many cases are we missing? Diabetologia. 2010;53(12):2504-2508. DOI 10.1007/s00125-010-1799-4. Slaugenhaupt S.A., Blumenfeld A., Gill S.P., Leyne M., Mull J., Cuajungco M.P., Liebert C.B., Chadwick B., Idelson M., Reznik L., Robbins C., Makalowska I., Brownstein M., Krappmann D., Scheidereit C., Maayan C., Axelrod F.B., Gusella J.F. Tissue-specific expression of a splicing mutation in the IKBKAP gene causes familial dysautonomia. Am. J. Hum. Genet. 2001;68(3):598-605. DOI 10.1086/318810. Spyer G., Macleod K.M., Shepherd M., Ellard S., Hattersley A.T. Pregnancy outcome in patients with raised blood glucose due to a heterozygous glucokinase gene mutation. Diabet. Med. 2009;26(1):14-18. DOI 10.1111/j.1464-5491.2008.02622.x. Steele A.M., Shields B.M., Wensley K.J., Colclough K., Ellard S., Hattersley A.T. Prevalence of vascular complications among patients with glucokinase mutations and prolonged, mild hyperglycemia. JAMA. 2014;311(3):279-286. DOI 10.1001/jama.2013.283980. Stenson P.D., Mort M., Ball E.V., Evans K., Hayden M., Heywood S., Hussain M., Phillips A.D., Cooper D.N. The Human Gene Mutation Database: towards a comprehensive repository of inherited mutation data for medical research, genetic diagnosis and next-generation sequencing studies. Hum. Genet. 2017;136(6):665-677. DOI 10.1007/ s00439-017-1779-6. Stride A., Vaxillaire M., Tuomi T., Barbetti F., Njølstad P.R., Hansen T., Costa A., Conget I., Pedersen O., Søvik O., Lorini R., Groop L., Froguel P., Hattersley A.T. The genetic abnormality in the beta cell determines the response to an oral glucose load. Diabetologia. 2002; 45(3):427-435. DOI 10.1007/s00125-001-0770-9. Thanabalasingham G., Owen K.R. Diagnosis and management of maturity onset diabetes of the young (MODY). BMJ. 2011;343:d6044. DOI 10.1136/bmj.d6044. The 1000 Genomes Project Consortium, Auton A., Brooks L.D., Durbin R.M., Garrison E.P., Kang H.M., Korbel J.O., Marchini J.L., McCarthy S., McVean G.A., Abecasis G.R.A global reference for human genetic variation. Nature. 2015;526(7571):68-74. DOI 10.1038/nature15393. Toaima D., Näke A., Wendenburg J., Praedicow K., Rohayem J., Engel K., Galler A., Gahr M., Lee-Kirsch M.A. Identification of novel GCK and HNF1A/TCF1 mutations and polymorphisms in German families with maturity-onset diabetes of the young (MODY). Hum. Mutat. 2005;25(5):503-504. DOI 10.1002/humu.9334. Wang K., Li M., Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164. DOI 10.1093/nar/gkq603. Wędrychowicz A., Tobór E., Wilk M., Ziółkowska-Ledwith E., Rams A., Wzorek K., Sabal B., Stelmach M., Starzyk J.B. Phenotype Heterogeneity in glucokinase-maturity-onset diabetes of the young (GCK-MODY) patients. J. Clin. Res. Pediatr. Endocrinol. 2017; 9(3):246-252. DOI 10.4274/jcrpe.4461. Xiong H.Y., Alipanahi B., Lee L.J., Bretschneider H., Merico D., Yuen R.K., Hua Y., Gueroussov S., Najafabadi H.S., Hughes T.R., Morris Q., Barash Y., Krainer A.R., Jojic N., Scherer S.W., Blencowe B.J., Frey B.J. The human splicing code reveals new insights into the genetic determinants of disease. Science. 2015;347(6218): 1254806. DOI 10.1126/science.1254806. Collection of the samples and clinical data, clinical examination, high-throughput sequencing, and Sanger sequencing have been supported by the Russian Science Foundation (grant No. 14-15-00496-P). Bioinformatic analyses have been supported as part of the budget topic in state assignment No. АААА-А19-119100990053-4. The English language was corrected and certified by shevchuk-editing.com. The authors thank the patients for participation in this study.