Skip to main content

Two novel variants in DYRK1B causative of AOMS3: expanding the clinical spectrum



We investigated pathogenic DYRK1B variants causative of abdominal obesity-metabolic syndrome 3 (AOMS3) in a group of patients originally diagnosed with type 2 diabetes. All DYRK1B exons were analyzed in a sample of 509 unrelated adults with type 2 diabetes and 459 controls, all belonging to the DMS1 SIGMA-cohort (ExAC). We performed in silico analysis on missense variants using Variant Effect Predictor software. To evaluate co-segregation, predicted pathogenic variants were genotyped in other family members. We performed molecular dynamics analysis for the co-segregating variants.


After filtering, Mendelian genotypes were confirmed in two probands bearing two novel variants, p.Arg252His and p.Lys68Gln. Both variants co-segregated with the AOMS3 phenotype in classic dominant autosomal inheritance with full penetrance. In silico analysis revealed impairment of the DYRK1B protein function by both variants. For the first time, we describe age-dependent variable expressivity of this entity, with central obesity and insulin resistance apparent in childhood; morbid obesity, severe hypertriglyceridemia, and labile type 2 diabetes appearing before 40 years of age; and hypertension emerging in the fifth decade of life. We also report the two youngest individuals suffering from AOMS3.


Monogenic forms of metabolic diseases could be misdiagnosed and should be suspected in families with several affected members and early-onset metabolic phenotypes that are difficult to control. Early diagnostic strategies and medical interventions, even before symptoms or complications appear, could be useful.


Abdominal obesity-metabolic syndrome 3 (AOMS3 [

OMIM:615812]) is a rare autosomal dominant disorder caused by pathogenic variants in the dual-specificity tyrosine phosphorylation-regulated kinase 1B gene (DYRK1B) located on chromosome 19q13.2 [1]. This monogenic form of metabolic syndrome (MetS) is characterized by abdominal obesity, type 2 diabetes, hypertension, and early-onset coronary artery disease [2]. DYRK1B inhibits the Sonic Hedgehog and WNT pathways, increasing the expression of master adipogenic transcription factors CCAAT/enhancer binding protein (C/EBP-alpha) and peroxisome proliferator-activated receptor gamma (PPAR-gamma) [2]. Moreover, DYRK1B induces the expression of glucose-6-phosphatase, a key enzyme in hepatic gluconeogenesis [3].

The identification of carriers of pathogenic variants in genes such as DYRK1B could be useful for establishing early diagnostic strategies and medical interventions in a reasonable number of affected individuals even before symptoms or complications appear. Until now, only two different mutations (p.Arg102Cys and p.His90Pro) in DYRK1B have been described as being causative of AOMS3 in three Iranian families and five unrelated Caucasian individuals [2]. However, families with rare monogenic forms of metabolic diseases could be misdiagnosed or even overlooked if causative variants are not directly explored, especially in populations with a high prevalence of these entities.

Therefore, we searched for DYRK1B variants in the exome sequencing data derived from 968 unrelated individuals (509 with type 2 diabetes) belonging to the DMS1 SIGMA-cohort (ExAC) [4], focusing on variants classified by Variant Effect Predictor (VEP) as deleterious and damaging to confirm their co-segregation with AOMS3. Here, we describe two novel DYRK1B mutations as causative of AOMS3 in two families previously misdiagnosed with type 2 diabetes.


Identification and co-segregation of pathogenic variants in the DYRK1B gene

Of the 968 unrelated individuals, 52.6% (n = 509) had type 2 diabetes according to the American Diabetes Association criteria [5]. The remaining 459 individuals were healthy subjects > 45 years old with fasting glucose (FG) levels < 100 mg/dL. We were able to identify 29 variants in DYRK1B (Additional file 1: Table S1), including seven missense variants and four variants predicted to be deleterious or damaging by VEP (Table 1). Two of these latter variants, p.Leu28Pro and p.Asp436Asn, were found in six heterozygote individuals (n = 5 and n = 1, respectively), all of them without type 2 diabetes or clinical characteristics of AOMS3. The remaining two missense variants, p.Arg252His and p.Lys68Gln, were found in two unrelated individuals manifesting symptoms suggestive of AOMS3, such as childhood-onset abdominal obesity, type 2 diabetes, hypertriglyceridemia, and arterial hypertension. Both of these individuals also had a family history of obesity and type 2 diabetes, as well as premature death secondary to myocardial infarction in first and second-degree relatives. We were able to recruit three generations of these two families with 26 traceable members, 9 of whom were DYRK1B pathogenic-variant carriers (5 harbored the p.Arg252His variant and 4 p.Lys68Gln). The pedigree analysis revealed co-segregation of the DYRK1B genotype with the AOMS3 phenotype, showing a characteristic autosomal dominant inheritance pattern with full penetrance (Fig. 1) and age-dependent variable expressivity. None of the non-carrier family members had type 2 diabetes or previous diagnosis of dyslipidemia, and only one had arterial hypertension. The carriers had severe metabolic compromise characteristic of AOMS3 (Table 2), and we observed a significant difference between carriers and non-carriers of pathogenic variants; carriers had higher body mass index (BMI), FG, and triglyceride levels. Furthermore, abdominal obesity started in childhood and progressed to morbid obesity in youth. The current mean BMI was 36.6 Kg/m2, but all patients had a higher BMI (mean 45.5 Kg/m2) before the diagnosis of type 2 diabetes. HOMA-IR and the Matsuda index (normal > 3) were impaired in all carriers, even in a child with normoglycemia (Additional file 1: Table S2). Similarly, the mean age of diagnosis of type 2 diabetes was 34 years old, and most of them had difficult to control disease (mean HbA1c: 87 mmol/mol (10.1%) and mean FG: 191.1 mg/dL). The mean serum triglyceride levels were 315.9 mg/dL, though all adult carriers underwent pharmacological treatment for dyslipidemia, except individual II.14 in family 2, who had serum triglyceride levels of 650 mg/dL.

Table 1 Variant Effect Predictor analysis of missense variants in DYRK1B
Fig. 1

Pedigree of two families with an autosomal dominant inheritance pattern of DYRK1B mutations associated with AOMS3. Family 1 with p.Arg252His (Arg: CGC/His: CAC). Family 2 showing the p.Lys68Gln mutation (Lys: AAG/Gln: CAG). The index cases (II.5 in family 1 and II.12 in family 2) are indicated by arrows. Genotypes are presented under each individual. Family members with clinical signs compatible with AOMS3 are indicated by solid symbols. Slashes indicate that the individual is deceased

Table 2 Clinical characteristics of the 26 individuals genotyped in two families

The clinical record documented that arterial hypertension was the last trait to become apparent, and it was present in 4 of 5 carriers older than 40 years old (the mean age of the diagnosis of arterial hypertension was 47 years). The pulse wave velocity increased in all carriers, even in the absence of arterial hypertension (Additional file 1: Table S2). p.Lys68Gln carriers had a higher BMI, lower insulin secretion, worse diabetes control, and an increased urine albumin/creatinine ratio than p.Arg252His carriers (Table 2 and Additional file 1: Figure S1). There were no differences in HDL, apolipoprotein B, creatinine, and hepatic enzymes levels between carriers and non-carriers. The uric acid level was elevated in two carriers, one with the p.Arg252His variant and another with p.Lys68Gln variant (Additional file 1: Table S2).

Interestingly, we were able to document metabolic disorders in two carrier males in their second decade of life. The youngest was an 11-year-old child who carried the p.Arg252His variant (IV.5 in family 1), had obesity (BMI: 30 Kg/m2, weight > 95th percentile of CDC grow charts), and demonstrated insulin resistance (HOMA-IR: 5.44) but not other metabolic abnormalities. The second carrier was an 18-year-old who carried the p.Lys68Gln variant (III.19 in family 2), was overweight (BMI: 26.4 Kg/m2, weight 90-95th percentile of CDC grow charts), serum FG 102 mg/dL, and hypertriglyceridemia since 16 years old; he was undergoing pharmacological treatment. Notably, both carriers had a strictly controlled diet because of the family history of type 2 diabetes.

Molecular dynamics

To gain more insight into the effect of the variants on protein structure and function, we performed molecular dynamics (MD) simulation trajectory analysis of DYRK1B-252Arg and DYRK1B-252His (Fig. 2). The root mean square deviation (RMSD) values indicate that DYRK1B-252His took a longer period of time to reach equilibration (Fig. 2a). The root mean square fluctuation (RMSF) analysis documented instability and structural changes in the N- and C-terminal regions, possibly due to the lack of hydrogen bonds (Fig. 2b, brown and yellow boxes). In addition, analyses revealed that, in contrast, to the results when Arg was present at position 252, His252 lost the hydrogen bonds with Pro552 and Ser554, resulting in conformational changes in the catalytic site and potentially losing functionality (Fig. 2c–f).

Fig. 2

Molecular dynamic analysis of DYRK1B-252His. a Root mean square deviation (RMSD) of atomic positions, analysis of wild-type protein (black line), and DYR1B-252His (red line). b Root mean square fluctuation (RMSF) of the wild-type protein (black line) and DYR1B-252His (red line). c Wild-type structure. d Structure of DYRK1B-252His. The protein structure colors are the same as previously described. e Region of interest zoomed-in for Arg252. The hydrogen bonds are depicted by the blue dots. f Same region of interest with the His252 variant

Despite the DYRK1B-68Gln mutant showed no-structural changes and relative stability, in silico analysis, revealed that this mutation is located within the nuclear localization sequence (NLS) proposed by Jhaisha et al. (positions 66–86) [6], suggesting that, although its catalytic function could not be affected, it could accumulate in the cytosol.


Over the past few years, there have been many attempts to gain more insight into the genetic factors involved in metabolic diseases [7,8,9]. Multiple genetic variants have been shown to participate in the pathogenesis of each of the traits of metabolic syndrome. The polygenic nature of these conditions implies that the effect of the majority of genetic variants in these disorders is small [7]. However, families in which autosomal dominant inheritance is present have been used to search for rare mutations in genes with a strong contribution [10, 11]. The availability of exome sequencing is leading to the rapid identification of new players in the pathogenesis of metabolic diseases. This is the case for DYRK1B mutations, which cause a rare monogenic form of MetS known as AOMS3 [2]. This syndrome has been described as the presence of abdominal obesity, type 2 diabetes, hypertension, and early-onset coronary artery disease [2]. Pathological DYRK1B variants result in the enhanced expression of transcription factors C/EBPalpha and PPARgamma, leading to increased adipogenesis. In addition, DYRK1B increases glucose-6-phosphatase, which is strongly associated with insulin resistance, explaining the metabolic phenotypes characterizing AOMS3 [2, 3].

The participation of DYRK1B in MetS is poorly studied. Six years after DYRK1B was associated with AOMS3, only a few carriers have been reported [1]. Furthermore, only two missense DYRK1B mutations (p.Arg102Cys and p.His90Pro) have been identified in these individuals [2]. It is possible for rare diseases that mimic symptoms of common diseases to be confused with them. The prevalence of metabolic diseases in Mexico is one of the highest in the world [12, 13], and rare metabolic diseases are often hidden behind them; therefore, families with rare monogenic forms can remain unnoticed.

Next generation sequencing technologies have greatly improved the possibility of identifying rare pathogenic variants involved in monogenic diseases [4, 11]. In this study, after analyzing the sequence of all DYRK1B exons in a sample of 968 adult, including 509 with type 2 diabetes (SIGMA-ExAC) [14], we found 29 variants. SIFT and PolyPhen predicted that four of them (p.Leu28Pro, p.Asp436Asn, p.Arg252His, and p.Lys68Gln) have a deleterious and damaging effect. The p.Leu28Pro variant was described previously as having a protective effect against type 2 diabetes in a phenome-wide association study [15]. In this study, we found five heterozygotes individuals with neither type 2 diabetes nor AOMS3, but we were not able to recruit the familial relatives to confirm its protective effect. The p.Asp436Asn variant was found in a male heterozygote, who was 98 years old and metabolically healthy. However, p.Arg252His and p.Lys68Gln exhibited co-segregation with the AOMS3 phenotype with classic dominant autosomal inheritance and full penetrance. Both variants were absent in the 1000 Genomes database [16], and the gnomAD database reports eight individuals with p.Arg252His variant and five with p.Lys68Gln variant, including the ones we report here. In Latinos, these variants were found with a frequency of 0.00008677 and 0.00005783, respectively [17].

MD structural analyses predicted that, when DYRK1B-252His was present, the formation of three hydrogen bonds was impaired, with instability in the N- and C-terminal regions. In contrast, the p.Lys68Gln variant did not produce any significant changes in the protein structure, although the NLS motif could be affected, in accordance with Kosugi et al., who showed that the N-terminal basic pattern “Lys68Arg69” is required for a strong NLS activity [18]. More over, the effect of these variants could be similar to those documenting that DYRK1B-102Cys and DYRK1B-90Pro variants cause changes to the structure and perinuclear aggregation but barely affect the kinase activity [6]. We classified p.Arg252His and p.Lys68Gln as causative of AOMS3, in agreement with the American College of Medical Genetics and Genomics standards and guidelines [19].

Compared to non-carriers of pathogenic variants, we found that carriers had higher BMI, and FG and triglyceride levels. Furthermore, pulse wave velocity was increased in all carriers, even in the absence of arterial hypertension, explaining the high cardiovascular risk found in this condition. Insulin action was decreased in all but one case, but the insulinogenic index was significantly decreased in all carriers, even in normoglycemic individuals, suggesting that the remarkable severity of hyperglycemia found in this condition results from a combination of moderate insulin resistance and a moderate to severe defect in insulin secretion. This could be supported by recent findings showing that DYRK1A and DYRK1B play an important role in human β cells proliferation [20].

Notably, we are describing additional features of the disease and manifestations that develop throughout life. Our findings exhibit age-dependent variance in expressivity in all patients, with some clinical features apparent at a very early age and other manifestations appearing later in life. Central obesity and insulin resistance started during childhood and then progressed rapidly to morbid obesity and labile type 2 diabetes. They also had severe hypertriglyceridemia with onset as a teenager. Similarly, AOMS3 patients developed hypertension in the fifth decade of life, and both families had a history of premature death due to cardiovascular events. Another interesting finding was that p.Lys68Gln carriers had a higher BMI, hypoinsulinemia, worse diabetes control and an increased urine albumin/creatinine ratio compared to the p.Arg252His carriers, suggesting an allelic heterogeneity.

Remarkably, none of the non-carriers had diabetes; only a 49-year-old woman (II.2 in family 2) had glucose serum levels of 109 mg/dl, along with obesity and hypertension. Other non-carrier individuals also had low levels of HDL or obesity, however, none of them met all the AOMS3 clinical features. These manifestations could be a reflection of the high prevalence of metabolic diseases in our population [12, 13].


How can these findings improve public health? AOMS3 should be suspected in individuals with the clinical outcome described above and a family history of affected first degree relatives. Assessment of the first-degree relatives should be performed routinely in subjects with extreme obesity. An autosomal dominant inheritance pattern should be systematically sought. Identifying variants involved in rare disorders could bring to light new genes associated with common diseases and may have implications for screening and targeted therapy. Otherwise, cascade genetic screening increases the possibility of finding pre-symptomatic carriers and lead to the implementation of strategies to prevent or delay the disease course, benefiting a reasonable number of affected families. The large number of cases and remarkably large structure of the Mexican families are unique opportunities to identify new genetic variants in understudied populations. Special care should be taken with the youngest members of the families, in whom cardiovascular prevention should be implemented early in life.

Material and methods

Study participants

We included 968 unrelated adult Mexican Mestizos belonging to the DMS1 SIGMA-cohort [14] who were previously sequenced by Sure-Select Human All Exon v2.0 (Illumina) and included in the ExAC project [4]. A peripheral blood sample was collected after fasting for at least 8 h. The following clinical and biochemical data were obtained for all participants of the DMS1 cohort using the Synchron CX5 Analyzer System (Beckman Coulter Fullerton, CA, USA): FG (mg/dL), HDL (mg/dL), and serum triglycerides (mg/dL). HbA1c levels were measured using the IN2it analyzer (Bio-Rad, Hercules, CA, USA). Blood pressure was measured using a digital blood pressure monitor (HEM-907XL, OMRON). Weight and height were measured using a body composition monitor (HBF-500 INT, OMRON) and electronic stadiometer (ADE Germany). Waist circumference was measured midway between the inferior margin of the ribs and the border of the iliac crest using a flexible clinical measuring tape.

The study was carried out according to the Declaration of Helsinki and was approved by the Research, Ethics, and Biosafety Human Committees of the Instituto Nacional de Medicina Genómica (INMEGEN) in Mexico City. All participants provided written informed consent. They were recruited from August 2017 to December 2018, all of them inhabited the Valley of Mexico.

Identification of pathogenic variants

We identified DYRK1B variants by analyzing the exon sequences in each individual. To find the deleterious and damaging variants, we annotated them using the VEP toolset [21]. Missense variants that were predicted as deleterious or affecting the protein structure were used to perform genotype–phenotype linkages. Individuals identified as carriers of these variants and had clinical manifestations suggesting AOMS3 were re-contacted and invited, along with their family members, to participate in a familial co-segregation study. All individuals who participated in the family segregation study provided written informed consent; in the case of children, the parents provided written consent and the children assented.

Genotyping was performed by Sanger sequencing using specific primers: p.Arg252His forward primer, CCTTTCTTCTCTGGCCAT; p.Arg252His reverse primer, ACCCAAACTACTAGCCGTGC; p.Lys69Gln forward primer, TGCCAGCAGCCTTACAGTT; p.Lys69Gln reverse primer, CCACTGCGCAACGATGTAGTC. The obtained sequences were analyzed by 4Peaks V1.8 [22].

The same clinical, demographic, and biochemical data described for index cases were obtained for each family participant. In addition, seven carriers (four p.Arg252His and three p.Lys69Gln) and three non-carriers were clinically and biochemically re-evaluated at the Unidad de Investigación de Enfermedades Metabólicas, Instituto Nacional de Ciencias Médicas y Nutrición. The assessment included fasting biochemical measurements (i.e., clinical chemistry, liver panel, and lipid profile), estimation of body composition, pulse wave velocity measurements, and albumin/creatinine ratio in a spot urine sample. An oral glucose tolerance test was performed using a 75 g oral glucose charge. Serum glucose and insulin were measured at fasting, 30-, 60-, 90-, 120-, and 180-min after oral glucose intake. The glucose concentration was measured by an automated glucose analyzer (Yellow Springs Instruments Co.). The serum insulin concentration was measured by a chemiluminescent immunoassay (Beckman Coulter Access 2) and HbA1c levels by HPLC (Variant II Turbo, BIORAD). Cholesterol, triglycerides, HDL, apolipoprotein B, uric acid, creatinine, and hepatic enzymes levels were measured using colorimetric assays (Unicel DxC 600 Synchron Clinical System Beckman Coulter). LDL was calculated by the Friedewald equation when the triglyceride concentration was < 250 mg/dL.

Molecular dynamics

The wild-type amino acid sequence of DYRK1B (Q9Y463) from the UniProt database [23] was modeled to obtain the 3D protein structure using the I-TASSER server [24]. The structures of mutated proteins were predicted using the predicted wild-type protein and VMD v1.9.3 software [25]. Next, we carried out an atomistic MD simulation with explicit atom representation for proteins, water, and ions under force-field using the CHARMM package and NAMD v2.3 software [26]. Periodic boundary conditions, particle mesh Ewald, and a non-bonded cut-off of 14 Å and 2 fs time step were used. The isothermal-isobaric conditions were maintained with a Langevin thermostat (310 K) and Langevin piston barostat (1 atm). For each model, the system was subjected to energy minimization for 1000 steps, followed by equilibration for 5 ns, and the simulation continued for 50 ns without restraints. In the simulations analyses, we used the VEGAZZ v3.1.2 [27], Carma [28], and R cran project v3.4 [29] programs. The MD analysis included the RMSD and the RMSF.

Statistical analysis

Clinical data are reported as mean ± standard deviation. Data were analyzed using R cran project v3.4 [29] and a general linear model test for comparing metabolite levels between carriers and non-carriers. Gender- and age-adjusted P < 0.05 was considered significant.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available due this work is part of a larger project but are available from the corresponding author on reasonable request.



Abdominal obesity-metabolic syndrome 3


B Dual-specificity tyrosine phosphorylation-regulated kinase 1B


Transcription factors CCAAT/enhancer binding protein


Fasting glucose


Instituto Nacional de Medicina Genómica


Molecular dynamics


Metabolic syndrome


Nuclear localization signal


Peroxisome proliferator-activated receptor gamma


Variant Effect Predictor


Root mean square deviation


Root mean square fluctuation


  1. 1.

    OMIM entry—# 615812—Abdominal obesity-metabolic syndrome 3; AOMS3 [cited 30 Apr 2019].

  2. 2.

    Keramati AR, Fathzadeh M, Go G-W, Singh R, Choi M, Faramarzi S, et al. A Form of the metabolic syndrome associated with mutations in. N Engl J Med. 2014;15(370):1909–19.

    Article  Google Scholar 

  3. 3.

    von Groote-Bidlingmaier F, Schmoll D, Orth HM, Joost HG, Becker W, Barthel A. DYRK1 is a co-activator of FKHR (FOXO1a)-dependent glucose-6-phosphatase gene expression. Biochem Biophys Res Commun. 2003;300(3):764–9.

    Article  Google Scholar 

  4. 4.

    Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;17(536):285.

    Article  Google Scholar 

  5. 5.

    Association AD. 2. Classification and diagnosis of diabetes: standards of medical care in diabetes—2019. Diabetes Care. 2019;42(Supplement 1):S13-28.

    Article  Google Scholar 

  6. 6.

    Jhaisha SA, Widowati EW, Kii I, Sonamoto R, Knapp S, Papadopoulos C, et al. DYRK1B mutations associated with metabolic syndrome impair the chaperone-dependent maturation of the kinase domain. Sci Rep. 2017;7(1):6420.

    Article  Google Scholar 

  7. 7.

    Zafar U, Khaliq S, Ahmad HU, Manzoor S, Lone KP. Metabolic syndrome: an update on diagnostic criteria, pathogenesis, and genetic links. Horm Athens Greece. 2018;17(3):299–313.

    Article  Google Scholar 

  8. 8.

    Saklayen MG. The global epidemic of the metabolic syndrome. Curr Hypertens Rep. 2018;20(2):12.

    Article  Google Scholar 

  9. 9.

    Sherling DH, Perumareddi P, Hennekens CH. Metabolic syndrome: clinical and policy implications of the new silent killer. J Cardiovasc Pharmacol Ther. 2017;22(4):365–7.

    CAS  Article  Google Scholar 

  10. 10.

    Serra-Juhé C, Martos-Moreno GÁ, BoudePieri F, Flores R, Chowen JA, Pérez-Jurado LA, et al. Heterozygous rare genetic variants in non-syndromic early-onset obesity. Int J Obes. 2020;44(4):830–41.

    Article  Google Scholar 

  11. 11.

    DiStefano JK, Kingsley CB. Identification of Disease Susceptibility Alleles in the Next Generation Sequencing Era. Methods Mol Biol Clifton NJ. 2018;1706:3–16.

    CAS  Article  Google Scholar 

  12. 12.

    Rojas R, Aguilar-Salinas CA, Jiménez-Corona A, Shamah-Levy T, Rauda J, Avila-Burgos L, et al. Metabolic syndrome in Mexican adults: results from the National Health and Nutrition Survey 2006. Salud Pública México. 2010;52(Suppl 1):S11-18.

    Article  Google Scholar 

  13. 13.

    Gutiérrez-Solis AL, Datta Banik S, Méndez-González RM. Prevalence of metabolic syndrome in mexico: a systematic review and meta-analysis. Metab Syndr Relat Disord. 2018;16(8):395–405.

    Article  Google Scholar 

  14. 14.

    SIGMA Type 2 Diabetes Consortium, Williams AL, Jacobs SBR, Moreno-Macías H, Huerta-Chagoya A, Churchhouse C, et al. Sequence variants in SLC16A11 are a common risk factor for type 2 diabetes in Mexico. Nature. 2014;506(7486):97–101.

    Article  Google Scholar 

  15. 15.

    Mirshahi T, Murray MF, Carey DJ. The metabolic syndrome and DYRK1B. N Engl J Med. 2014;371(8):784–5.

    Article  Google Scholar 

  16. 16.

    1000 Genomes|A deep catalog of human genetic variation [cited 19 Nov 2019].

  17. 17.

    DYRK1B | gnomAD v2.1.1 | gnomAD. [cited 1 Apr 2021].

  18. 18.

    Kosugi S, Hasebe M, Matsumura N, Takashima H, Miyamoto-Sato E, Tomita M, et al. Six classes of nuclear localization signals specific to different binding grooves of importin alpha. J Biol Chem. 2009;284(1):478–85.

    CAS  Article  Google Scholar 

  19. 19.

    Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med Off J Am Coll Med Genet. 2015;17(5):405–24.

    Google Scholar 

  20. 20.

    Ackeifi C, Swartz E, Kumar K, Liu H, Chalada S, Karakose E et al. Pharmacologic and genetic approaches define human pancreatic β cell mitogenic targets of DYRK1A inhibitors. American Society for Clinical Investigation; 2020 [cited 1 Apr 2021].

  21. 21.

    McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GRS, Thormann A, et al. The ensembl variant effect predictor. Genome Biol. 2016;17(1):122.

    Article  Google Scholar 

  22. 22.

    4Peaks: For peaks, four peaks. The DNA sequence trace viewer for OS X. [cited 19 Nov 2019].

  23. 23.

    UniProt. [cited 19 Nov 2019].

  24. 24.

    Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015;12(1):7–8.

    CAS  Article  Google Scholar 

  25. 25.

    Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14(1):33–8.

    CAS  Article  Google Scholar 

  26. 26.

    Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, et al. Scalable molecular dynamics with NAMD. J Comput Chem. 2005;26(16):1781–802.

    CAS  Article  Google Scholar 

  27. 27.

    Pedretti A, Villa L, Vistoli G. VEGA: a versatile program to convert, handle and visualize molecular structure on Windows-based PCs. J Mol Graph Model. 2002;21(1):47–9.

    CAS  Article  Google Scholar 

  28. 28.

    Glykos NM. Software news and updates. Carma: a molecular dynamics analysis program. J Comput Chem. 2006;27(14):1765–8.

    CAS  Article  Google Scholar 

  29. 29.

    R: The R Project for Statistical Computing. [cited 19 Nov 2019].

Download references


The authors would like to thank the staff of the Immunogenomics and Metabolic Disease Laboratory/INMEGEN and the Endocrinology and Metabolism Department/Instituto Nacional de Ciencias Médicas y Nutrición for all their support, particularly Juan Luis Jiménez-Ruiz, María Guadalupe Salas-Martínez, Yolanda Saldaña-Alvarez, Luz Elizabeth Guillen-Pineda, María Del Carmen Moreno-Villatoro, and Adriana Cruz-Lopez. We are thankful to the study volunteers for all their work and support throughout the realization of the study.


This study was partially supported by Instituto Carlos Slim de la Salud, A.C.

Author information




ECM-C and FB-O were involved in all aspects of the study and wrote the first draft of the manuscript. LO contributed to the study design, analysis, interpretation of data, and editing of the manuscript. EM and II-F contributed to the molecular dynamics analysis and interpretation of data. CAA-S, DVG-V, IC-A, AR-R, AM-H, CC–C, SI-A, CZ, and HG-O contributed to the study design, interpretation of data, and sample collection. All authors were involved in writing the paper and gave final approval of the submitted and published versions. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lorena Orozco.

Ethics declarations

Ethics approval and consent to participate

The study was carried out according to the Declaration of Helsinki and was approved by the Research, Ethics, and Biosafety Human Committees of the Instituto Nacional de Medicina Genómica (INMEGEN) in Mexico City.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Figure S1. Glucose tolerance curves. Table S1. Variants on DYRK1B in a sample of 968 Mexican Mestizos. Table S2. Pulse wave velocity and other biochemical studies.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mendoza-Caamal, E.C., Barajas-Olmos, F., Mirzaeicheshmeh, E. et al. Two novel variants in DYRK1B causative of AOMS3: expanding the clinical spectrum. Orphanet J Rare Dis 16, 291 (2021).

Download citation


  • AOMS3
  • Diabetes
  • DYRK1B
  • Metabolic syndrome
  • Monogenic, and Obesity