Next Article in Journal
Novel Divergent Members of the Kitrinoviricota Discovered through Metagenomics in the Intestinal Contents of Red-Backed Voles (Clethrionomys gapperi)
Next Article in Special Issue
The ABA/LANCL Hormone/Receptor System in the Control of Glycemia, of Cardiomyocyte Energy Metabolism, and in Neuroprotection: A New Ally in the Treatment of Diabetes Mellitus?
Previous Article in Journal
Genome-Wide Analysis of the Rab Gene Family in Melilotus albus Reveals Their Role in Salt Tolerance
Previous Article in Special Issue
Strategies for Therapeutic Amelioration of Aberrant Plasma Zn2+ Handling in Thrombotic Disease: Targeting Fatty Acid/Serum Albumin-Mediated Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study

1
College of Health and Life Sciences (CHLS), Hamad Bin Khalifa University (HBKU), Qatar Foundation (QF), Doha P.O. Box 34110, Qatar
2
Medical and Population Genomics Lab, Sidra Medicine, Doha P.O. Box 26999, Qatar
3
Qatar Genome Program (QGP), Qatar Foundation Research, Development and Innovation, Qatar Foundation (QF), Doha P.O. Box 5825, Qatar
4
Qatar Metabolic Institute, Hamad Medical Corporation, Doha P.O. Box 3050, Qatar
5
Bioinformatics Core, Weill Cornell Medicine-Qatar, Education City, Doha P.O. Box 24144, Qatar
6
Department of Biophysics and Physiology, Weill Cornell Medicine, New York, NY 10065, USA
7
Centre for Genomic and Experimental Medicine, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh EH4 2XU, UK
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(1), 130; https://doi.org/10.3390/ijms24010130
Submission received: 14 November 2022 / Revised: 9 December 2022 / Accepted: 14 December 2022 / Published: 21 December 2022
(This article belongs to the Special Issue Diabetes: Molecular Mechanisms)

Abstract

:
Maturity-onset diabetes of the young (MODY) is a rare monogenic form of diabetes mellitus. In this study, we estimated the prevalence and genetic spectrum of MODY in the Middle Eastern population of Qatar using whole-genome sequencing (WGS) of 14,364 subjects from the population-based Qatar biobank (QBB) cohort. We focused our investigations on 14 previously identified genes ascribed to the cause of MODY and two potentially novel MODY-causing genes, RFX6 and NKX6-1. Genetic variations within the 16 MODY-related genes were assessed for their pathogenicity to identify disease-causing mutations. Analysis of QBB phenotype data revealed 72 subjects (0.5%) with type 1 diabetes, 2915 subjects (20.3%) with type 2 diabetes and 11,377 (79.2%) without diabetes. We identified 22 mutations in 67 subjects that were previously reported in the Human Genetic Mutation Database (HGMD) as disease-causing (DM) or likely disease causing (DM?) for MODY. We also identified 28 potentially novel MODY-causing mutations, predicted to be among the top 1% most deleterious mutations in the human genome, which showed complete (100%) disease penetrance in 34 subjects. Overall, we estimated that MODY accounts for around 2.2–3.4% of diabetes patients in Qatar. This is the first population-based study to determine the genetic spectrum and estimate the prevalence of MODY in the Middle East. Further research to characterize the newly identified mutations is warranted.

1. Introduction

Maturity-onset diabetes of the young (MODY) is a rare monogenic form of diabetes mellitus, accounting for approximately 1–5% of diabetes cases worldwide [1]. MODY is an autosomal dominant condition, clinically characterized by a multigenerational family history of diabetes and absence of pancreatic β-cell autoantibodies, and disease onset before the age of 25 years in most cases. A prominent feature of MODY includes progressive β-cell dysfunction, while in some cases it could be associated with other organ abnormalities, such as kidney abnormalities and fat malabsorption [2]. Patients with MODY are often misdiagnosed due to similarities in their clinical manifestations to those with other forms of diabetes [3]. To date, 14 genes are ascribed to the cause of MODY, each encoding for a different subtype [2]. Recent studies have identified novel variants in two genes encoding regulatory factor X6 (RFX6) and NK6 Homeobox 1 (NKX6-1) as possible MODY contributors [4,5].
MODY screening and profiling in clinic is challenging due to the absence of a single clinical criterion, the occurrence of de novo mutations and the unavailability of suitable genetic testing [6]. The genetic architecture of diabetes is complex. For instance, genetic variations in several MODY genes have been associated with type 2 diabetes mellitus (T2DM) in genome-wide association studies (GWAS) [7,8]. In addition, homozygous mutations in some MODY genes can cause other types of diabetes, such as permanent neonatal diabetes and antibody-negative type 1 diabetes mellitus (T1DM) caused by mutations in the insulin gene (INS; MODY10) [9,10]. Moreover, homozygous mutations in ABCC8 (MODY12) and KCNJ11 (MODY13) also cause neonatal diabetes [11]. Studies have shown that approximately 80% of MODY patients were misclassified as either T1DM or T2DM, leading to suboptimal therapeutic intervention [3]. Prior studies have, however, revealed certain unique phenotypic features that assist in the identification of suspected MODY patients, such as low apolipoprotein levels in MODY1, and renal cysts and diabetes syndrome in MODY5 [12,13,14]. In view of this, understanding the genetic architecture and physiological anomalies associated with MODY are essential and molecular genetic testing is crucial for identifying patients with MODY.
The prevalence of MODY varies amongst different populations [15]. The most frequent subtype in the United States, Japan and Turkey is GCK-MODY2, accounting for approximately 50% of the identified variants, followed by HNF1A-MODY3 and HNF1B-MODY5 [16,17,18,19,20]. In contrast, HNF1A-MODY3 is the most common MODY subtype among the studied populations in South Africa, South Korea and in Europe, including the United Kingdom, Germany and Poland [1,15,21,22]. MODY prevalence in the Middle East and North Africa (MENA) region remains largely unexplored [23]. Although there have been reports of MODY prevalence in Kuwait, Oman and Tunisia from suspected cases, a definitive population-based prevalence estimate has not been determined.
Pathogenic and likely pathogenic variants have been previously identified in 2.1% of pediatric diabetes patients [24]. Notably, a recent clinic-based study of MODY in children from Qatar has shown that MODY affects 4.56 children per 100,000 in Qatar [25]. However, this study focused primarily on children between the ages of 0 and 18 years and included children from different ethnicities attending a single hospital [25]. In this study, we estimate the prevalence and determine the genetic spectrum of MODY in Qatar using whole-genome sequencing (WGS) data from the population-based cohort of Qatar Biobank (QBB) participants (n = 14,364). We sought to determine the spectrum of MODY-related mutations and identify highly penetrant Qatari-specific mutations. Our findings provide important insights into the genetic architecture of MODY in the Middle East and have the potential for improving the diagnosis and clinical management of MODY patients in Qatar.

2. Results

2.1. Study Overview

The overview of the study design is illustrated in Figure 1. This study was based on WGS and phenotypic data of 14,364 participants from the population-based QBB cohort [26]. Analysis of phenotype data showed that 72 participants (0.5%) of the study cohort were classified as T1DM, 2915 (20.3%) as T2DM, while 11,377 (79.2%) were classified as non-diabetes. The clinical characteristics of the study population are listed in Table 1. Female participants constituted 55.8% of the study population, while the age of the overall study cohort ranged from 18 to 89 years. Notably, the prevalence of obesity in Qatar is remarkably high (77.6% of the total study population were overweight or obese). Additionally, 74.1% of non-diabetes participants and 91.0% of diabetes patients in our study were overweight or obese. Only 35 subjects (0.24%) met the clinical criteria for diagnosing MODY based on being diabetic with normal BMI and the age of onset before <25 (only 63% of diabetes patients declared age of disease onset). Therefore, the widely accepted clinical criteria for MODY diagnosis was not accurate when applied to our study cohort, given that most patients were overweight or obese, and genetic testing is critical in MODY diagnosis.

2.2. Identification of Previously Known MODY Causative Variants

We first investigated previously known variants associated with MODY onset. We investigated variants that were previously reported as disease-causing (DM) or likely disease-causing (DM?) mutations for MODY in the human gene mutation database (HGMD) and classified as pathogenic or likely pathogenic in Qiagen Clinical Insight (QCI) and/or ClinVar databases (Table 2). QCI classification criteria are based on the 2015/2016 ACMG/AMP Guidelines. Variants were also annotated by SIFT and PolyPhen-2 databases to indicate the potential functional effect of each mutation. We detected 10 previously identified missense MODY mutations in 24 participants (Table 2). The majority of detected variants were located in HNF1A. However, rs587778397 (p.Arg177Trp) in HNF1A and rs766934515 (p.Val113Gly) in BLK were the most frequently occurring mutations. Additional variants in HNF4A, GCK, KLF11 and ABCC8 were also recorded. Moreover, 10 subjects carried mutations in the HNF1A gene, 5 in BLK and 4 in GCK. Complete penetrance was observed in HNF4A (rs1555817727; p.Arg312Cys) and HNF1A (rs137853238; p.Arg272His).
In addition, we also identified variants of uncertain significance (VUS) that are predicted to rank in the top 1% of most deleterious mutations in the human genome based on CADD Phred scores (>20). We identified 12 missense VUS in 43 subjects (Table 3). Most subjects were heterozygous carriers of variants in HNF1A (n = 26), followed by GCK (n = 5) and BLK (n = 5) genes. HNF1A (rs202039659; p.His514Arg) was the most frequent mutation, recorded in 11 subjects. However, complete penetrance was only recorded for HNF4A (rs371124358; p.Arg310Gln and rs377151067; p.V404I). Combined, we found that previously known MODY mutations are mostly detected in the HNF1A gene, accounting for 53.7% of cases, followed by BLK and GCK mutations accounting for approximately 15% and 13% of cases, respectively. Overall, we identified 22 MODY-causing mutations in 8 genes, accounting for 2.2% of diabetes participants (Table 2 and Table 3).

2.3. Identifying Potentially Novel MODY-Causing Mutations

Next, we investigated potentially novel MODY-causing mutations based on complete disease penetrance (100%) and a CADD Phred score >20. We detected 28 mutations in 14 genes in 34 diabetes subjects (Table 4). None of these variants are reported as MODY-causing mutations in the HGMD or ClinVar databases, while the majority are predicted to be damaging by Polyphen-2 and deleterious by SIFT. In addition, 18 mutations are not reported in dbSNP or the Genome Aggregation Database (gnomAD). All detected novel mutations were missense, while variants in HNF1A (MODY3) and CEL (MODY8) were the most common, detected in eight subjects. Moreover, we also detected variants in rare MODY-related genes, including NEUROD1 (MODY6), KCNJ11 (MODY13) and APPL1 (MODY14), among other MODY subtypes. Notably, we also identified three variants in RFX6 and two variants in NKX6-1 in seven subjects (Figure 2). Importantly, all potentially novel mutations were detected in patients classified as T2DM based on phenotype data.

2.4. Estimating the Prevalence of MODY in Qatar

We detected 22 previously known MODY-causing mutations in 67 subjects, yielding a prevalence of MODY at 2.2% in individuals with diabetes (Figure 1). In addition, we identified 28 potentially novel MODY-causing mutations, accounting for ~56% of the overall reported mutations in 34 subjects and 1.14% of diabetes cases. Of note, combining subjects carrying mutations in both known and potentially novel MODY-causing mutations in our study population yielded an overall prevalence of 3.4% of diabetes patients in Qatar (Figure 1). Overall, previously reported and potentially novel MODY mutations were detected in 101 subjects. All detected mutations were heterozygous. The proportions of subjects carrying mutations in each MODY subtype are illustrated in Figure 3. We found that HNF1A-MODY3 is the most prevalent MODY subtype in the Qatari population, accounting for approximately 40% of cases. In addition, we also compared the clinical characteristics of subjects carrying MODY-related mutations (Supplementary Table S1). The identified MODY patients were predominately overweight and obese (83.1%), while they also showed strong diabetic medical history (95%). Moreover, the gender distribution of MODY cases showed higher prevalence in females (54.5%). In addition, we also compared the clinical parameters and genetic matrices of MODY-related genes between non-diabetes, T1DM and T2DM patients (Supplementary Table S2 and Supplementary Figure S1). Based on clinical parameters, half of the participants who carried mutations in MODY-related genes were classified as T2DM patients and half were non-diabetes based on phenotype data, while one patient was classified as T1DM (Supplementary Table S2). The genetic profile of these subjects also showed that the majority of them carried mutations in HNF1A (Supplementary Figure S1).

3. Discussion

The etiopathogenesis and clinical manifestations of MODY are broadly analogous to other forms of diabetes, such as youth-onset T2DM, which renders clinical risk assessment and genetic testing critical for accurate diagnosis and therapy selection. However, due to the challenges associated with genetic screening, MODY is underreported in the MENA region. In Europe, MODY accounts for 1–5% of all diabetes mellitus cases [28]. We estimated the prevalence of MODY in Qatar at 2.2% of diabetes patients, while inclusion of potentially novel mutations yielded a higher proportion of MODY, accounting for 3.4% of cases. The MENA region has one of the highest and fastest increasing levels of diabetes prevalence compared to global estimates [29]. We estimated the prevalence of diabetes in the population-based QBB cohort as ~21%, in agreement with previous reports (~17–20%) [30,31,32].
The prevalence of MODY mutations among children with diabetes from different ethnicities was estimated at 4.56 per 100,000 individuals (~1.5% of diabetes cases) according to a recent single clinic-based study conducted in Qatar [25]. Among the 12 mutations detected in children, 9 mutations were also detected in our study cohort. These included one previously reported mutation in HNF1A (rs202039659; p.His514Arg) and a previously known mutation in GCK (rs1375656631; p.Ala259Thr). However, seven additional mutations were detected in our study cohort but they did not meet our analysis criteria, such as HGMD and Clinvar classifications. Overall, the differences in reported variants among MODY cases in children compared to Qatari adults could be attributed to differences in study design, given that our study was a large population-based cohort of Qatari adults (18–89 years), whereas the previous study was a small study of children with different ethnicities attending a single diabetes clinic in Qatar.
We identified 22 previously reported and 28 potentially novel MODY-causing mutations in 101 participants, which accounted for 3.4% of diabetes patients. The estimated prevalence is in accordance with the widely accepted overall prevalence of MODY, accounting for 1–5% of all diabetes cases. Combined, mutations in HNF1A were the most common, detected in 40 cases (39.6%); followed by BLK (10.9%) and GCK (9.9%) genes. HNF1A MODY-causing mutations are the most prevalent in European, North American, and Asian populations; accounting for 20–50% of all MODY cases [3,6,33,34,35]. We also detected several known but rare MODY-causative mutations in different genes including PDX1-MODY4, KLF11-MODY7 and ABBC8-MODY12. However, many subjects carrying previously reported mutations did not present the clinical symptoms of diabetes, which could be attributed to incomplete penetrance of MODY mutations, which varies with differences in age and in genetic architecture between populations [6,36]. MODY is well known for its clinical and genetic heterogeneity with variable disease penetrance. Notably, the associations between BLK mutations and MODY penetrance has been indecisive. Five MODY causing mutations in BLK were identified in three families including p.A71T [37]. However, p.A71T was detected in 52 normoglycemic subjects in a European cohort and showed weak associations with increased T2DM risk [38].
We identified 28 potentially novel MODY-causing mutations in 34 participants (~1.1% of diabetes cases), accounting for ~56% of the overall reported mutations. Twenty-three mutations were in MODY-causing genes, while five mutations were detected in potentially novel MODY-causing genes—three variants in RFX6 and two variants in NKX6-1 in seven subjects. Several studies have associated variants in these genes with MODY progression [4]. Notably, the transcription factor RFX6 plays a vital role in the development of endocrine pancreas and in the differentiation of islet cells, including insulin-producing cells [39]. To date, only a few potential MODY-causing variants of RFX6 have been identified [4,5,40]. Mutations in the RFX6 gene, for instance, c.1954C>T (p.Arg652X), have been found to reduce insulin secretion and insulin-dependent insulinotropic polypeptide responses [5,40]. Moreover, GWAS have associated variants of RFX6 with blood glucose levels [41] and T2DM [42]. Similarly, NKX6.1 has been identified as a critical transcription factor in the early and late development of pancreatic β-cells and in maintaining β-cells’ identity, and it regulates many genes involved in glucose uptake, metabolism and insulin biosynthesis [43]. Moreover, variants in NKX6-1 have also been associated with T2DM [44]. Notably, only two mutations in NKX6-1, P329L and S317L have been described in South India and associated with MODY development [4]. These genes could have a profound effect on MODY onset, but further genetic and functional investigations are warranted.
Despite recent advances in technology, the genetic and molecular basis of MODY onset remains to be fully elucidated. The complexity and broad spectrum of phenotypes associated with MODY, including de novo mutations and genetic background, make it challenging to conduct a clinical manifestation profile for all subtypes. However, identification of population-specific mutations attributed to MODY has the potential for the development of diagnostic or prognostic microarrays as a more feasible option. For instance, multiplex ligation-dependent probe amplification led to the identification of two de novo mutations (deletions) in the HNF1B gene in suspected MODY cases in a Greek cohort [45].
Our findings highlight the genetic anomalies associated with MODY in Qatar based on a substantially large cohort. Importantly, our findings also reiterate the importance of genetic testing for MODY diagnosis given that our cohort predominantly comprised subjects classified as overweight or obese. Moreover, of the 101 suspected MODY subjects who carried mutations in MODY-related genes, 84 patients (83.2%) were overweight or obese, consistent with high obesity prevalence in Qatar. Therefore, the widely accepted clinical criteria for MODY diagnosis is not accurate when applied to our study cohort given that most patients were overweight or obese, and genetic testing is critical in MODY diagnosis.
Our study was, however, limited by the lack of detection for structural variants of MODY-associated genes, although MODY is predominantly caused by missense mutations. Moreover, additional clinicopathological data to investigate associations between specific mutations and disease severity, age of diagnosis and mode of inheritance (inherited or de novo) were not available. In addition, details about autoantibody detection to stratify patients into diabetes subtypes were also unavailable. Overall, the reported variants and mutations may be explored further in functional analysis to confirm their pathogenicity in other populations and in the development of diagnostic and prognostic tools for MODY.

4. Materials and Methods

4.1. Study Participants

Qatar Biobank (QBB) was the primary source of study participants. QBB aims to gather information from the native Qatari population, as well as from long-term residents (≥15 years), with follow-ups scheduled every 5 years. The study cohort analyzed in this study consisted of 14,364 Qatari participants aged between 18 to 89 years. Informed consent was obtained from all participants prior to inclusion in the study. The study was approved by the institutional review boards of QBB (E-2019-QF-QBB-RES-ACC-0179-0104) and Hamad Bin Khalifa University (approval no. 2021-03-078). Participants underwent extensive medical examination, filled in a standardized questionnaire and self-reported past medical histories, diet, physical activity and lifestyle [26]. Physical measurements, such as the body mass index (BMI), and biological samples, including blood, urine and saliva, were also collected. Measurements relevant to diabetes were investigated by the Qatar National Health Service (Hamad Medical Corporation, Doha, Qatar) to assess the participants’ health status by evaluating clinical diagnostic biomarkers, including glucose, hemoglobin A1c (HbA1c) and C-peptide levels.

4.2. Phenotypic Data and Patient Classification

Study participants were identified as diabetes or non-diabetes based on blood biomarkers (C-peptide, HbA1c, glucose), self-reported diabetes and self-reported administration of anti-diabetes medications. Serum C-peptide levels were determined by the sandwich electrochemiluminescence immunoassay using Elecsys C-Peptide kit (Roche, Basel, Switzerland), while HbA1c in blood was measured using turbidimetric inhibition immunoassay (TINIA) utilizing Tina-quant HbA1c Gen. 3 kit (Roche) on hemolyzed blood samples and random glucose levels in serum were measured using the enzymatic reference method with GLUC3 glucose hexokinase kit (Roche), utilizing a COBAS instrument (Roche).
We classified our study cohort based on the phenotypic data provided by QBB. We first identified T1DM as self-declared diabetes subjects receiving insulin treatment exclusively and with serum C-peptide levels < 0.5 ng/mL. Subjects who did not meet the criteria for T1DM were classified as T2DM if they met any of the following criteria: self-declared diabetes, self-reported diabetes medication or with HbA1c > 6.5%. Subjects who did not meet the T1DM or T2DM criteria were classified as non-diabetes subjects.

4.3. Whole-Genome Sequencing

The WGS data of the QBB study participants was provided by the Qatar Genome Project (QGP) [46]. Details of WGS and quality control measures have been described previously [47]. Briefly, genomic DNA from peripheral blood samples was extracted using Qiagen MIDI kit (Qiagen, Hilden, Germany), following the manufacturer’s protocol and using the automated QIASymphony SP instrument (Qiagen). DNA integrity was evaluated using Caliper Labchip GXII (Perkin Elmer, Waltham, MA, USA) Genomic DNA assay and was quantified using the Quant-iT dsDNA Assay (Invitrogen, Waltham, MA, USA). Whole-genome libraries were prepared using Illumina TruSeq DNA Nano kit (Illumina, San Diego, CA, USA). HiSeq X Ten (Illumina) was used to sequence genomic libraries to achieve a minimum average coverage of 30×. Library construction and sequencing was performed at Sidra Clinical Genomics Laboratory Sequencing Facility (Sidra Medicine, Doha, Qatar). The generated files were subjected to quality control using FastQC (v0.11.2) and reads were aligned to the CRCh38 reference genome. Quality control on mapped reads was performed using Picard (v1.117). A combined variant call file (gVCF) was generated for all study subjects, which contained all genetic variations detected in the QBB study participants.

4.4. Bioinformatics Analysis to Identify MODY-Causing Mutations

Genetic variants of the following known MODY-related genes—HNF4A, GCK, HNF1A, PDX1, HNF1B, NEUROD1, KLF11, CEL, PAX4, INS, BLK, ABCC8, KCNJ11, APPL1, RFX6 and NKX6-1—were extracted from the gVCF file, based on genomic coordinates obtained from the University of California Santa Cruz (UCSC) Genome Browser [48]. We included genetic variants with minor allele frequency < 0.001 and mean sequencing depth > 20×. The average sequencing depth of MODY-related variants was 32.1×. We extracted genetic variants based on their location in each of the MODY-causing genes and related the genotype data to the clinical phenotype data for stratification into MODY subtypes. Variant annotations were performed using Human Gene Mutation Databases (HGMD) [49] to identify previously reported disease-causing mutations (DM) or likely disease-causing mutations (DM?). The clinical significance of variants was investigated using the QIAGEN Clinical Insight (Qiagen) and ClinVar database [50]. Potentially novel MODY-related mutations were determined based on a Combined Annotation-Dependent Depletion (CADD) [51] Phred score > 20 and complete disease penetrance (100%), implying that the subject had diabetes and the genetic variant showed a CADD Phred score > 20 (predicated to be among the top 1% most deleterious mutations in the human genome). Penetrance for each variant was calculated by dividing the number of carriers showing the phenotype (T1DM or T2DM as described above) by the overall number of mutation carriers (N). The functional annotations of the identified variants was performed using Sorting Intolerant From Tolerant (SIFT) [52] and Polymorphism Phenotyping v2 (Polyphen-2) [53] tools.

4.5. Statistical Analysis

Statistical analyses were performed using GraphPad Prism 9 software (GraphPad Software, San Diego, CA, USA). Shapiro–Wilk normality test was carried out to determine data distribution. Unpaired t-test was used for data comparisons between groups on normally distributed data, while Mann–Whitney test was used for datasets that did not show normal distribution. Chi-square (χ2) test was performed to determine the differences between binary variables within study cohorts. A p value of < 0.05 was considered statistically significant. Data are presented as mean ± standard deviation (SD).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24010130/s1.

Author Contributions

A.A.E.: data curation, formal analysis, investigation and writing—original draft. S.M.T.: investigation, writing—review and editing and visualization. I.D. and Y.A.-S.: investigation and methodology. S.T., K.S. and A.B.A.-S.: funding acquisition, resources, writing—review and editing. O.M.E.A.: conceptualization, funding acquisition, formal analysis, supervision, investigation and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by a grant from the Qatar National Research Fund (QNRF) awarded to OMEA, KS and ABA-S (NPRP11C-0115-180010).

Institutional Review Board Statement

The study was approved by the institutional Review Boards of Qatar Biobank (QBB), Doha, Qatar (Protocol no. E-2019-QF-QBB-RES-ACC-0179-0104) and Hamad Bin Khalifa University (Approval no. 2021-03-078).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data analyzed in this study are subject to the following licenses/restrictions: the raw whole-genome sequence data from Qatar Biobank are protected and are not available for deposition into public databases due to data privacy laws. Access to QBB/QGP phenotype and whole-genome sequence data can be obtained through an ISO-certified protocol, which involves submitting a project request at https://www.qatarbiobank.org.qa/research/how-apply, subject to approval by the Institutional Review Board of the QBB. Requests to access these datasets should be directed to https://www.qatarbiobank.org.qa/research/how-apply.

Acknowledgments

We are grateful to the Qatar Biobank (QBB) and the Qatar Genome Program Research (QGPR) Consortium for providing us the genomic and phenotypic data. We are also thankful to all participants of this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Skoczek, D.; Dulak, J.; Kachamakova-Trojanowska, N. Maturity Onset Diabetes of the Young-New Approaches for Disease Modelling. Int. J. Mol. Sci. 2021, 22, 7553. [Google Scholar] [CrossRef] [PubMed]
  2. Hu, M.; Cherkaoui, I.; Misra, S.; Rutter, G.A. Functional Genomics in Pancreatic beta Cells: Recent Advances in Gene Deletion and Genome Editing Technologies for Diabetes Research. Front. Endocrinol. 2020, 11, 576632. [Google Scholar] [CrossRef] [PubMed]
  3. 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, 2504–2508. [Google Scholar] [CrossRef] [PubMed]
  4. Mohan, V.; Radha, V.; Nguyen, T.T.; Stawiski, E.W.; Pahuja, K.B.; Goldstein, L.D.; Tom, J.; Anjana, R.M.; Kong-Beltran, M.; Bhangale, T.; et al. Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India. BMC Med. Genet. 2018, 19, 22. [Google Scholar] [CrossRef] [Green Version]
  5. Patel, K.A.; Kettunen, J.; Laakso, M.; Stančáková, A.; Laver, T.W.; Colclough, K.; Johnson, M.B.; Abramowicz, M.; Groop, L.; Miettinen, P.J.; et al. Heterozygous RFX6 protein truncating variants are associated with MODY with reduced penetrance. Nat. Commun. 2017, 8, 888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Peixoto-Barbosa, R.; Reis, A.F.; Giuffrida, F.M.A. Update on clinical screening of maturity-onset diabetes of the young (MODY). Diabetol. Metab. Syndr. 2020, 12, 50. [Google Scholar] [CrossRef] [PubMed]
  7. Bonnefond, A.; Boissel, M.; Bolze, A.; Durand, E.; Toussaint, B.; Vaillant, E.; Gaget, S.; Graeve, F.D.; Dechaume, A.; Allegaert, F.; et al. Pathogenic variants in actionable MODY genes are associated with type 2 diabetes. Nat. Metab. 2020, 2, 1126–1134. [Google Scholar] [CrossRef]
  8. Mahajan, A.; Taliun, D.; Thurner, M.; Robertson, N.R.; Torres, J.M.; Rayner, N.W.; Payne, A.J.; Steinthorsdottir, V.; Scott, R.A.; Grarup, N.; et al. Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat. Genet. 2018, 50, 1505–1513. [Google Scholar] [CrossRef] [Green Version]
  9. Stoy, J.; Edghill, E.L.; Flanagan, S.E.; Ye, H.; Paz, V.P.; Pluzhnikov, A.; Below, J.E.; Hayes, M.G.; Cox, N.J.; Lipkind, G.M.; et al. Insulin gene mutations as a cause of permanent neonatal diabetes. Proc. Natl. Acad. Sci. USA 2007, 104, 15040–15044. [Google Scholar] [CrossRef] [Green Version]
  10. Vedovato, N.; Cliff, E.; Proks, P.; Poovazhagi, V.; Flanagan, S.E.; Ellard, S.; Hattersley, A.T.; Ashcroft, F.M. Neonatal diabetes caused by a homozygous KCNJ11 mutation demonstrates that tiny changes in ATP sensitivity markedly affect diabetes risk. Diabetologia 2016, 59, 1430–1436. [Google Scholar] [CrossRef]
  11. Flanagan, S.E.; Dung, V.C.; Houghton, J.A.L.; De Franco, E.; Ngoc, C.T.B.; Damhuis, A.; Ashcroft, F.M.; Harries, L.W.; Ellard, S. An ABCC8 Nonsense Mutation Causing Neonatal Diabetes through Altered Transcript Expression. J. Clin. Res. Pediatr. Endocrinol. 2017, 9, 260–264. [Google Scholar] [CrossRef] [PubMed]
  12. Mughal, S.A.; Park, R.; Nowak, N.; Gloyn, A.L.; Karpe, F.; Matile, H.; Malecki, M.T.; McCarthy, M.I.; Stoffel, M.; Owen, K.R. Apolipoprotein M can discriminate HNF1A-MODY from Type 1 diabetes. Diabet. Med. 2013, 30, 246–250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ræder, H.; Haldorsen, I.S.; Ersland, L.; Grüner, R.; Taxt, T.; Søvik, O.; Molven, A.; Njølstad, P.R. Pancreatic Lipomatosis Is a Structural Marker in Nondiabetic Children with Mutations in Carboxyl-Ester Lipase. Diabetes 2007, 56, 444–449. [Google Scholar] [CrossRef] [Green Version]
  14. Ricci, P.; Magalhães, P.; Krochmal, M.; Pejchinovski, M.; Daina, E.; Caruso, M.R.; Goea, L.; Belczacka, I.; Remuzzi, G.; Umbhauer, M.; et al. Urinary proteome signature of Renal Cysts and Diabetes syndrome in children. Sci. Rep. 2019, 9, 2225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kim, S.H. Maturity-Onset Diabetes of the Young: What Do Clinicians Need to Know? Diabetes Metab. J. 2015, 39, 468–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Relan, S.; Hessel, F.C.; Timmons, K.D.; Sonabend, R.; Redondo, M.J. Maturity-Onset Diabetes of the Young (MODY) in a Racially Diverse Pediatric Population. Diabetes 2018, 67, 1529-P. [Google Scholar] [CrossRef]
  17. Yorifuji, T.; Fujimaru, R.; Hosokawa, Y.; Tamagawa, N.; Shiozaki, M.; Aizu, K.; Jinno, K.; Maruo, Y.; Nagasaka, H.; Tajima, T.; et al. Comprehensive molecular analysis of Japanese patients with pediatric-onset MODY-type diabetes mellitus. Pediatr. Diabetes 2012, 13, 26–32. [Google Scholar] [CrossRef]
  18. Yorifuji, T.; Higuchi, S.; Kawakita, R.; Hosokawa, Y.; Aoyama, T.; Murakami, A.; Kawae, Y.; Hatake, K.; Nagasaka, H.; Tamagawa, N. Genetic basis of early-onset, maturity-onset diabetes of the young-like diabetes in Japan and features of patients without mutations in the major MODY genes: Dominance of maternal inheritance. Pediatr. Diabetes 2018, 19, 1164–1172. [Google Scholar] [CrossRef]
  19. Carmody, D.; Naylor, R.N.; Bell, C.D.; Berry, S.; Montgomery, J.T.; Tadie, E.C.; Hwang, J.L.; Greeley, S.A.W.; Philipson, L.H. GCK-MODY in the US National Monogenic Diabetes Registry: Frequently misdiagnosed and unnecessarily treated. Acta Diabetol. 2016, 53, 703–708. [Google Scholar] [CrossRef] [Green Version]
  20. Doğan, M.; Eröz, R.; Bolu, S.; Yüce, H.; Gezdirici, A.; Arslanoğlu, İ.; Teralı, K. Study of ten causal genes in Turkish patients with clinically suspected maturity-onset diabetes of the young (MODY) using a targeted next-generation sequencing panel. Mol. Biol. Rep. 2022, 49, 7483–7495. [Google Scholar] [CrossRef]
  21. Hwang, J.S.; Shin, C.H.; Yang, S.W.; Jung, S.Y.; Huh, N. Genetic and clinical characteristics of Korean maturity-onset diabetes of the young (MODY) patients. Diabetes Res. Clin. Pract. 2006, 74, 75–81. [Google Scholar] [CrossRef] [PubMed]
  22. Matsha, T.E.; Raghubeer, S.; Tshivhase, A.M.; Davids, S.F.G.; Hon, G.M.; Bjorkhaug, L.; Erasmus, R.T. Incidence of HNF1A and GCK MODY Variants in a South African Population. Appl. Clin. Genet. 2020, 13, 209–219. [Google Scholar] [CrossRef] [PubMed]
  23. Saraswathi, S.; Al-Khawaga, S.; Elkum, N.; Hussain, K. A Systematic Review of Childhood Diabetes Research in the Middle East Region. Front. Endocrinol. 2019, 10, 805. [Google Scholar] [CrossRef] [PubMed]
  24. Johnson, S.R.; Ellis, J.J.; Leo, P.J.; Anderson, L.K.; Ganti, U.; Harris, J.E.; Curran, J.A.; McInerney-Leo, A.M.; Paramalingam, N.; Song, X.; et al. Comprehensive genetic screening: The prevalence of maturity-onset diabetes of the young gene variants in a population-based childhood diabetes cohort. Pediatr. Diabetes 2019, 20, 57–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Haris, B.; Saraswathi, S.; Al-Khawaga, S.; Hasnah, R.; Saeed, A.; Mundekkadan, S.; Hamed, N.; Afyouni, H.; Abdel-Karim, T.; Mohammed, S.; et al. Epidemiology, genetic landscape and classification of childhood diabetes mellitus in the State of Qatar. J. Diabetes Investig. 2021, 12, 2141–2148. [Google Scholar] [CrossRef]
  26. Al Thani, A.; Fthenou, E.; Paparrodopoulos, S.; Al Marri, A.; Shi, Z.; Qafoud, F.; Afifi, N. Qatar Biobank Cohort Study: Study Design and First Results. Am. J. Epidemiol. 2019, 188, 1420–1433. [Google Scholar] [CrossRef]
  27. Weir, C.B.; Jan, A. BMI Classification Percentile and Cut Off Points; StatPearls: Treasure Island, FL, USA, 2022. [Google Scholar]
  28. Nkonge, K.M.; Nkonge, D.K.; Nkonge, T.N. The epidemiology, molecular pathogenesis, diagnosis, and treatment of maturity-onset diabetes of the young (MODY). Clin. Diabetes Endocrinol. 2020, 6, 20. [Google Scholar] [CrossRef]
  29. Sun, H.; Saeedi, P.; Karuranga, S.; Pinkepank, M.; Ogurtsova, K.; Duncan, B.B.; Stein, C.; Basit, A.; Chan, J.C.N.; Mbanya, J.C.; et al. IDF Diabetes Atlas: Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res. Clin. Pract. 2022, 183, 109119. [Google Scholar] [CrossRef]
  30. Al-Thani, M.H.; Al-Mutawa, K.A.; Alyafei, S.A.; Ijaz, M.A.; Khalifa, S.A.H.; Kokku, S.B.; Mishra, A.C.M.; Poovelil, B.V.; Soussi, M.B.; Toumi, A.A.; et al. Characterizing epidemiology of prediabetes, diabetes, and hypertension in Qataris: A cross-sectional study. PLoS ONE 2021, 16, e0259152. [Google Scholar] [CrossRef]
  31. International Diabetes Federation. IDF Diabetes Atlas, 10th ed.; International Diabetes Federation: Brussels, Belgium, 2021. [Google Scholar]
  32. International Diabetes Federation. IDF Middle East and North Africa Region; International Diabetes Federation: Brussels, Belgium, 2022. [Google Scholar]
  33. Pihoker, C.; Gilliam, L.K.; Ellard, S.; Dabelea, D.; Davis, C.; Dolan, L.M.; Greenbaum, C.J.; Imperatore, G.; Lawrence, J.M.; Marcovina, S.M.; et al. Prevalence, characteristics and clinical diagnosis of maturity onset diabetes of the young due to mutations in HNF1A, HNF4A, and glucokinase: Results from the SEARCH for Diabetes in Youth. J. Clin. Endocrinol. Metab. 2013, 98, 4055–4062. [Google Scholar] [CrossRef]
  34. Shields, B.M.; McDonald, T.J.; Ellard, S.; Campbell, M.J.; Hyde, C.; Hattersley, A.T. The development and validation of a clinical prediction model to determine the probability of MODY in patients with young-onset diabetes. Diabetologia 2012, 55, 1265–1272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Yang, Y.S.; Kwak, S.H.; Park, K.S. Update on Monogenic Diabetes in Korea. Diabetes Metab. J. 2020, 44, 627–639. [Google Scholar] [CrossRef] [PubMed]
  36. Sirugo, G.; Williams, S.M.; Tishkoff, S.A. The Missing Diversity in Human Genetic Studies. Cell 2019, 177, 26–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Borowiec, M.; Liew, C.W.; Thompson, R.; Boonyasrisawat, W.; Hu, J.; Mlynarski, W.M.; El Khattabi, I.; Kim, S.H.; Marselli, L.; Rich, S.S.; et al. Mutations at the BLK locus linked to maturity onset diabetes of the young and beta-cell dysfunction. Proc. Natl. Acad. Sci. USA 2009, 106, 14460–14465. [Google Scholar] [CrossRef] [Green Version]
  38. Bonnefond, A.; Yengo, L.; Philippe, J.; Dechaume, A.; Ezzidi, I.; Vaillant, E.; Gjesing, A.P.; Andersson, E.A.; Czernichow, S.; Hercberg, S.; et al. Reassessment of the putative role of BLK-p.A71T loss-of-function mutation in MODY and type 2 diabetes. Diabetologia 2013, 56, 492–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Lu, J.; Cheng, C.; Cheng, Z.C.; Wu, Q.; Shen, H.; Yuan, M.X.; Zhang, B.; Yang, J.K. The dual role of RFX6 in directing beta cell development and insulin production. J. Mol. Endocrinol. 2021, 66, 129–140. [Google Scholar] [CrossRef] [PubMed]
  40. Imaki, S.; Iizuka, K.; Horikawa, Y.; Yasuda, M.; Kubota, S.; Kato, T.; Liu, Y.; Takao, K.; Mizuno, M.; Hirota, T.; et al. A novel RFX6 heterozygous mutation (p.R652X) in maturity-onset diabetes mellitus: A case report. J. Diabetes Investig. 2021, 12, 1914–1918. [Google Scholar] [CrossRef]
  41. Sinnott-Armstrong, N.; Tanigawa, Y.; Amar, D.; Mars, N.; Benner, C.; Aguirre, M.; Venkataraman, G.R.; Wainberg, M.; Ollila, H.M.; Kiiskinen, T.; et al. Genetics of 35 blood and urine biomarkers in the UK Biobank. Nat. Genet. 2021, 53, 185–194. [Google Scholar] [CrossRef]
  42. Ray, D.; Chatterjee, N. A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between Type 2 Diabetes and Prostate Cancer. PLoS Genet. 2020, 16, e1009218. [Google Scholar] [CrossRef]
  43. Aigha, I.I.; Abdelalim, E.M. NKX6.1 transcription factor: A crucial regulator of pancreatic beta cell development, identity, and proliferation. Stem Cell Res. Ther. 2020, 11, 459. [Google Scholar] [CrossRef]
  44. Vujkovic, M.; Keaton, J.M.; Lynch, J.A.; Miller, D.R.; Zhou, J.; Tcheandjieu, C.; Huffman, J.E.; Assimes, T.L.; Lorenz, K.; Zhu, X.; et al. Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis. Nat. Genet. 2020, 52, 680–691. [Google Scholar] [CrossRef] [PubMed]
  45. Tatsi, E.B.; Kanaka-Gantenbein, C.; Scorilas, A.; Chrousos, G.P.; Sertedaki, A. Next generation sequencing targeted gene panel in Greek MODY patients increases diagnostic accuracy. Pediatr. Diabetes 2020, 21, 28–39. [Google Scholar] [CrossRef] [PubMed]
  46. Mbarek, H.; Devadoss Gandhi, G.; Selvaraj, S.; Al-Muftah, W.; Badji, R.; Al-Sarraj, Y.; Saad, C.; Darwish, D.; Alvi, M.; Fadl, T.; et al. Qatar genome: Insights on genomics from the Middle East. Hum. Mutat. 2022, 43, 499–510. [Google Scholar] [CrossRef] [PubMed]
  47. Thareja, G.; Al-Sarraj, Y.; Belkadi, A.; Almotawa, M.; Ismail, S.; Al-Muftah, W.; Badji, R.; Mbarek, H.; Darwish, D.; Fadl, T.; et al. Whole genome sequencing in the Middle Eastern Qatari population identifies genetic associations with 45 clinically relevant traits. Nat. Commun. 2021, 12, 1250. [Google Scholar] [CrossRef] [PubMed]
  48. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Stenson, P.D.; Mort, M.; Ball, E.V.; Chapman, M.; Evans, K.; Azevedo, L.; Hayden, M.; Heywood, S.; Millar, D.S.; Phillips, A.D.; et al. The Human Gene Mutation Database (HGMD®): Optimizing its use in a clinical diagnostic or research setting. Hum. Genet. 2020, 139, 1197–1207. [Google Scholar] [CrossRef] [PubMed]
  50. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.R.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Jang, W.; et al. ClinVar: Improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018, 46, D1062–D1067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Kircher, M.; Witten, D.M.; Jain, P.; O’Roak, B.J.; Cooper, G.M.; Shendure, J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 2014, 46, 310–315. [Google Scholar] [CrossRef] [Green Version]
  52. Ng, P.C.; Henikoff, S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003, 31, 3812–3814. [Google Scholar] [CrossRef] [Green Version]
  53. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Method. 2010, 7, 248–249. [Google Scholar] [CrossRef]
Figure 1. Study Design. This study was based on whole-genome sequencing (WGS) data of the Qatar Biobank (QBB) participants. The analysis cohort comprised n = 14,364 subjects, classified into diabetes mellitus subtypes (T1DM; n = 72 and T2DM; n = 2915) and non-diabetes subjects (n = 11,377) based on phenotypic data. The detected mutations were classified into previously known or potentially novel MODY-causing mutations based on pathogenicity annotations in the human gene mutation database (HGMD), QIAGEN Clinical Insight and ClinVar databases or Combined Annotation-Dependent Depletion (CADD) scores, as well as on disease (diabetes mellitus) penetrance.
Figure 1. Study Design. This study was based on whole-genome sequencing (WGS) data of the Qatar Biobank (QBB) participants. The analysis cohort comprised n = 14,364 subjects, classified into diabetes mellitus subtypes (T1DM; n = 72 and T2DM; n = 2915) and non-diabetes subjects (n = 11,377) based on phenotypic data. The detected mutations were classified into previously known or potentially novel MODY-causing mutations based on pathogenicity annotations in the human gene mutation database (HGMD), QIAGEN Clinical Insight and ClinVar databases or Combined Annotation-Dependent Depletion (CADD) scores, as well as on disease (diabetes mellitus) penetrance.
Ijms 24 00130 g001
Figure 2. Identified mutations in RFX6 and NKX6-1. The locations of potentially novel MODY-causing mutations in (A) RFX6 and (B) NKX6-1 genes. The number of participants carrying each mutation is presented by a dot above each gene mutation.
Figure 2. Identified mutations in RFX6 and NKX6-1. The locations of potentially novel MODY-causing mutations in (A) RFX6 and (B) NKX6-1 genes. The number of participants carrying each mutation is presented by a dot above each gene mutation.
Ijms 24 00130 g002
Figure 3. Proportion of subjects carrying MODY-causing mutations. Donut charts represent the proportions of subjects carrying: (A) previously reported pathogenic and likely pathogenic MODY-causing mutations (n = 28 subjects); (B) previously reported disease-causing (DM) or likely disease-causing (DM?) variants in the Human Gene Mutation Database, predicted to be among the top 1% most deleterious mutations in the human genome (n = 43 subjects); (C) novel fully penetrant, potentially MODY-causing mutations detected and predicted to be among the top 1% most deleterious mutations in the human genome (n = 34 subjects). (D) Donut chart illustrates the overall proportions of subjects (n = 101) carrying the previously known and potentially novel MODY-causing mutations. All genes are color coded, while percentages of subjects carrying each mutation subtype are also shown.
Figure 3. Proportion of subjects carrying MODY-causing mutations. Donut charts represent the proportions of subjects carrying: (A) previously reported pathogenic and likely pathogenic MODY-causing mutations (n = 28 subjects); (B) previously reported disease-causing (DM) or likely disease-causing (DM?) variants in the Human Gene Mutation Database, predicted to be among the top 1% most deleterious mutations in the human genome (n = 43 subjects); (C) novel fully penetrant, potentially MODY-causing mutations detected and predicted to be among the top 1% most deleterious mutations in the human genome (n = 34 subjects). (D) Donut chart illustrates the overall proportions of subjects (n = 101) carrying the previously known and potentially novel MODY-causing mutations. All genes are color coded, while percentages of subjects carrying each mutation subtype are also shown.
Ijms 24 00130 g003
Table 1. Clinical characteristics of study cohort.
Table 1. Clinical characteristics of study cohort.
Type 1 DiabetesType 2 DiabetesNon-Diabetes
Number of subjects 72 (0.5%)2915 (20.3%)11,377 (79.2%)
Age (mean ± SD)38.43 ± 14.7 *51.76 ± 11.8 *35.05 ± 11.7
Male34 (47.2%)1216 (41.7%)5104 (44.9%)
Female38 (52.8%)1699 (58.3%)6273 (55.1%)
BMI (kg/m²) **28.6 ± 5.132.3 ± 6.0 *28.9 ± 6.0
Underweight-3 (0.1%)262 (2.3%)
Normal weight18 (25.0%)236 (8.1%) *2672 (23.5%)
Overweight28 (38.9%)889 (30.5%) *4047 (35.6%)
Obese25 (34.7%)1777 (61.0%) *4385 (38.5%)
N/A1 (1.4%)10 (0.3%)11 (0.1%)
HbA1c (%)8.8 ± 1.5 *7.4 ± 1.8 *5.3 ± 0.4
C-peptide (ng/mL)0.32 ± 0.14 *2.83 ± 1.67 *2.02 ± 1.4
Family History of Diabetes
Father33 (45.8%) *1357 (46.5%) *4790 (42.1%)
Mother33 (45.8%) *1810 (62.1%) *5113 (44.9%)
Both18 (25.0%) *966 (33.1) *2504 (22.0%)
Treatment
Insulin72 (100.0%)645 (22.1%)-
Tablets-2011 (69.0%)-
Study cohort comprised of 14,364 subjects. * Statistically significant (p < 0.001) compared to non-diabetes controls. BMI: body mass index. ** Subjects with BMI below 18.5 kg/m2 were classified as underweight, between 18.5 to 24.9 kg/m2 as normal weight, between 25 to 29.9 kg/m2 as overweight and greater than or equal to 30 kg/m2 as obese [27]. N/A, not available.
Table 2. Identification of previously reported pathogenic and likely pathogenic MODY-causing mutations.
Table 2. Identification of previously reported pathogenic and likely pathogenic MODY-causing mutations.
Gene *MODYIDRef/AltProtein ChangeNMAFPenetrance (%)Diagnosis (n)CADDHGMDQCIClinVarSIFTPolyPhen-2
HNF4A1rs1555817727C/Tp.Arg312cys10.000035100T2DM30DMPLPDP.D
GCK2rs1375656631C/Tp.Ala259Thr30.00010466.6N.D (1)
T2DM (2)
24.4DMPPDBenign
GCK2rs104894005C/Gp.Glu279Gln10.0000350N.D22.2DMLPPTBenign
HNF1A3rs137853238G/Ap.Arg272his10.000035100T2DM31DMPPDP.D
HNF1A3rs587778397C/Tp.Arg177Trp50.00017420N.D (4)
T2DM (1)
25.7DMLPVUSTP.D
HNF1A3rs371717826C/Gp.Pro379arg30.0001040N.D (3)26.7DMPN/ADP.D
HNF1A3rs754729248C/Gp.Pro379Ala10.0000350N.D (1)25DMVUSPDP.D
KLF117rs121912645G/Tp.Ala347Ser10.0000350N.D (1)< 10DMPPTBenign
BLK11rs766934515T/Gp.Val113Gly50.0001740N.D (5)29.9DM?LPN/ADP.D
ABCC812rs761862121T/Gp.Lys889Thr30.0001040N.D (3)27.5DMPVUSDP.D
* Identified variants are reported as disease-causing (DM) or likely disease-causing (DM?) in the Human Gene Mutation Database (HGMD) and classified as pathogenic (P) or likely pathogenic (LP) in Qiagen Clinical Insight (QCI) and/or ClinVar databases. Penetrance is shown as percentage of mutation carriers with diabetes. Diagnosis is based on phenotype data. Ref, reference allele; Alt, alternative allele; N, number of heterozygous subjects; MAF, minor allele frequency; T2DM, type 2 diabetes; N.D, nondiabetes; CADD, Combined Annotation-Dependent Depletion Phred score; QCI, QIAGEN Clinical Insight classification based on the 2015/2016 ACMG/AMP Guidelines; VUS, variant of uncertain significance; D, deleterious; T, tolerated; P.D, probably damaging; N/A, not available.
Table 3. Identification of previously reported disease-causing or likely disease-causing variants in the Human Gene Mutation Database, predicted to be among the top 1% most deleterious mutations in the human genome.
Table 3. Identification of previously reported disease-causing or likely disease-causing variants in the Human Gene Mutation Database, predicted to be among the top 1% most deleterious mutations in the human genome.
GeneMODYIDRef/AltProtein ChangeNMAFPenetrance (%) *Diagnosis (n)CADDHGMDQCIClinVarSIFTPolyPhen-2
HNF4A1rs769007443G/Ap.Gly64Arg10.0000350N.D25.9DMVUSN/ADP.D
HNF4A1rs371124358G/Ap.Arg310Gln10.000035100T2DM23.5DMVUSVUSTBenign
HNF4A1rs377151067G/Ap.V404I10.000035100T2DM22.8DMVUSN/ATBenign
GCK2rs193922285C/Ap.Met462Ile50.00017440N.D (3)
T2DM (2)
22.9DMVUSVUSTBenign
HNF1A3rs201095611G/Ap.Ala161Thr60.00020916.67N.D (5)
T1DM (1)
27.3DMVUSVUSDP.D
HNF1A3rs867576513T/Cp.Met283Thr60.00020933.3N.D (4)
T2DM (2)
25.6DMVUSN/ADP.D
HNF1A3rs368683806G/Ap.Gly415Arg10.0000350N.D31DMVUSVUSDP.D
HNF1A3rs577078110C/Ap.His505Asn20.0000700N.D (2)26.3DMVUSN/ADP.D
HNF1A3rs202039659A/Gp.His514Arg110.00038327.2N.D (8)
T2DM (3)
24.8DMVUSVUSDP.D
PDX14rs753249965G/Ap.Gly55Asp30.0001040N.D (3)23.7DMVUSN/ADBenign
HNF1B5rs113042313C/Tp.Gly370Ser10.0000350N.D21DMVUSN/ATBenign
BLK11rs368427116C/Tp.Thr270Met50.0001740N.D (5)25.4DMVUSN/ATP.D
* Penetrance is shown as percentage of mutation carriers with diabetes. Diagnosis is based on phenotype data. Ref, reference allele; Alt, alternative allele; N, number of heterozygous subjects; MAF, minor allele frequency; T2DM, type 2 diabetes; N.D, nondiabetes; CADD, Combined Annotation-Dependent Depletion Phred score; HGMD, Human Gene Mutation Database; DM, disease-causing mutation; QCI, QIAGEN Clinical Insight classification based on the 2015/2016 ACMG/AMP Guidelines; VUS, variant of uncertain significance; N/A, not available; D, deleterious; T, tolerated; P.D, probably damaging.
Table 4. Discovery of potentially novel fully penetrant MODY-causing mutations in the Qatar Biobank (QBB) cohort.
Table 4. Discovery of potentially novel fully penetrant MODY-causing mutations in the Qatar Biobank (QBB) cohort.
GeneMODYChr:PosIDRef/
Alt
Protein ChangeN *MAFCADDgnomAD
MAF
SIFTPolyPhen-2
HNF4A120:44418486-T/Cp.Met237Thr10.00003526.2-DP.D
HNF4A120:44424203rs776656815C/Gp.Leu360Val10.00003522.80.0000239DBenign
GCK27:44145572-A/Tp.Met393Lys10.00003526.2-DBenign
HNF1A312:120993541-G/Ap.Gly183Glu10.00003526.1-DP.D
HNF1A312:120979082-A/Tp.Glu105Val20.00007025.2-DBenign
HNF1A312:120999538-C/Tp.Ser591Phe10.00003524.1-DP.D
PDX1413:27924690rs1357043267G/Cp.Glu281Gln30.00010423.10.000119DBenign
HNF1B517:37744622-G/Ap.Thr88Ile10.00003525.4-DBenign
HNF1B517:37731607rs755951130T/Cp.Asn345Asp10.00003522.60.0000081TBenign
NEUROD162:181678436-G/Ap.Thr142Ile10.00003528.2-DP.D
KLF1172:10046337-A/Gp.Asp77Gly10.00003527.9-DP.D
KLF1172:10048367-C/Tp.Pro344Ser10.00003523.2-DP.D
KLF1172:10046240-A/Tp.Met45Leu10.00003520.1-TBenign
CEL89:133070564rs766487195G/Ap.Gly467Arg10.00003526.10.0000080DP.D
CEL89:133064413rs748643667G/Ap.Val29Met10.00003525.30.0000040DP.D
CEL89:133066888rs773198000C/Gp.Ser243Arg10.00003522.90.0000242DP.D
CEL89:133069087-A/Gp.Thr375Ala10.00003522.3-DBenign
BLK118:11554800rs775313404A/Gp.Tyr177Cys10.00003529.40.0000119DP.D
ABCC81211:17404576rs769818698C/Tp.Val1165Met20.00007023.10.0000955DBenign
ABCC81211:17427908-C/Tp.Gly692Glu10.00003522.7-TBenign
KCNJ111311:17387440-G/Tp.Gln131Lys10.00003524.2-TP.D
KCNJ111311:17387761rs867211548C/Tp.Val24Ile10.00003520.7-TBenign
APPL1143:57269587-A/Tp.Asn677Ile10.00003521.4-TBenign
RFX6-6:116922098rs1485759457G/Ap.Val462Met10.00003526.1-DP.D
RFX6-6:116927531rs762356403C/Tp.Ser797Leu10.00003522.90.0000081DBenign
RFX6-6:116924786-A/Cp.Asn558Thr10.00003522.4-TBenign
NKX6-1-4:84498074-G/Ap.Ser52Phe20.00007024.2-DBenign
NKX6-1-4:84498145rs369821275C/Gp.Met28Ile20.00007023.50.0000525DBenign
* All identified variants showed complete disease penetrance (100%) and mutation carriers were classified as type 2 diabetes based on phenotype data. Penetrance is shown as a percentage of mutation carriers with diabetes. Chr: Pos, Chromosome: Position of identified variant based on genome build CRCh38; Ref, reference allele; Alt, alternative allele; N, number of heterozygous subjects; MAF, minor allele frequency; CADD, Combined Annotation-Dependent Depletion Phred score; gnomAD, Genome Aggregation Database; D, deleterious; T, tolerated; P.D, probably damaging.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Elashi, A.A.; Toor, S.M.; Diboun, I.; Al-Sarraj, Y.; Taheri, S.; Suhre, K.; Abou-Samra, A.B.; Albagha, O.M.E. The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study. Int. J. Mol. Sci. 2023, 24, 130. https://doi.org/10.3390/ijms24010130

AMA Style

Elashi AA, Toor SM, Diboun I, Al-Sarraj Y, Taheri S, Suhre K, Abou-Samra AB, Albagha OME. The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study. International Journal of Molecular Sciences. 2023; 24(1):130. https://doi.org/10.3390/ijms24010130

Chicago/Turabian Style

Elashi, Asma A., Salman M. Toor, Ilhame Diboun, Yasser Al-Sarraj, Shahrad Taheri, Karsten Suhre, Abdul Badi Abou-Samra, and Omar M. E. Albagha. 2023. "The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study" International Journal of Molecular Sciences 24, no. 1: 130. https://doi.org/10.3390/ijms24010130

APA Style

Elashi, A. A., Toor, S. M., Diboun, I., Al-Sarraj, Y., Taheri, S., Suhre, K., Abou-Samra, A. B., & Albagha, O. M. E. (2023). The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study. International Journal of Molecular Sciences, 24(1), 130. https://doi.org/10.3390/ijms24010130

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop