Gábor Firneisz1,2*
Klara Rosta3,4
János Rigó4Ákos Nádasdi1
Jürgen Harreiter5
Alexandra Kautzky-Willer5Anikó Somogyi1- 12nd Department of Internal Medicine, Semmelweis University, Budapest, Hungary
- 2MTA-SE Molecular Medicine Research Group, Hungarian Academy of Sciences - Semmelweis University, Budapest, Hungary
- 3Department of Obstetrics and Gynecology, Medical University of Vienna, Vienna, Austria
- 41st Department of Obstetrics and Gynecology, Semmelweis University, Budapest, Hungary
- 5Gender Medicine Unit, Division of Endocrinology and Metabolism, Department of Internal Medicine III, Medical University of Vienna, Vienna, Austria
Introduction
Genome-wide association studies (GWASs) are reliable tools to identify novel gene variant–trait (disease) associations. Recently high ranked journals published exciting reviews on the potential in the health predictions based on GWAS results (Tam et al., 2019). This indicates a renewed interest in the clinical use of polygenic risk and hazard scores (PRS and PHS) for prevalent diseases.
Despite the potential of PRS/PHS that can be applied to evidence the genetic component of specific clinical and treatment phenotypes is suggested to be substantial (Lambert et al., 2019; Loscalzo, 2019; Meigs, 2019; Srinivasan et al., 2019) we are not aware of any clinical praxis guideline in the diabetes mellitus (DM) field today that would have had been implemented to help clinical decisions based on GWAS findings.
This indicates a slow progression of clinical utilization of the genetic association study results. This slow transition to clinics is particularly contrasting when compared to standard diabetes care which is improved greatly due to the number of newly introduced glucose-lowering drugs and also to the novel follow-up/diagnostic techniques (Davies et al., 2018).
There is an enormous context gap between the everyday routine clinics and the 3,730 published (all traits) GWASs identified a total of 37,730 single nucleotide variants (SNVs) and 52,415 unique SNV–trait associations (till January, 2019) (Tam et al., 2019).
Polygenic Risk and Hazard Scores
The most often proposed clinical use of complex genetic data obtained from GWAS is the development of polygenic risk or hazard scores (PRSs/PHSs) (Sarraju and Knowles, 2019). Genome-wide PRS and PHS scores are summed up of genetic risks associated to many individual independent genetic effects and developed to predict a clinical outcome, e.g., a development of a disease with polygenic background (Khera et al., 2018; Florez, 2019).
The literature outlined that the effect of individual gene variants on complex common disease phenotype is relatively weak and in these cases the PRS might be useful (Lambert et al., 2019; Loscalzo, 2019; Meigs, 2019; Srinivasan et al., 2019). Effect sizes of common genetic variants for T2DM are modest ever since the early GWAS era and are usually ranging from 10 to 35% (Lyssenko and Laakso, 2013). Consistently the OR is typically in between 1.1 and 1.3 and the OR for the identified common gene variants for T2DM remains below 1.7 (Prasad and Groop, 2015).
However, due to that GWASs require high sample sizes the studied prevalent disease eventually may be composed of largely heterogeneous subgroups (Ahlqvist et al., 2018). In addition as the sample sizes are increasing further novel rare and/or small effect variants might be found. The task is further complicated by social, population diversity, sex differences, economic, and ethical challenges. Till now the genetic risk predictions have had no added value to prediction models for T2DM based on simple clinical risk factors, provided that GRS have been applied to the “general” population (Meigs et al., 2008). There is a shift in the understanding with a key element that rather the principal disease traits (i.e., sub-phenotypes or “endophenotypes”) should be put into the focus of the genetic prediction models than the “general population” to gain success in the clinical application of GRS. It is estimated that only 30% of (T2DM) patients might fall within the top 10% of a disease trait (“cluster”) associated GRS (Florez, 2019). Therefore, we may conclude that significant associations and the added value of the genetic prediction models may remain hidden with the “usual approach” when the analysis is only applied to less narrow, more heterogenous, “too general” populations, instead of assessing disease subtypes or key traits referring to specific mechanisms. Others also concluded that T2DM genetic risk score is not sufficient to improve health or health-economic, rather the better identification of metabolic sub-phenotypes with stronger genetic background might improve the opportunities for health translation (Meigs, 2019).
Few authors already argue that it is essential to stratify people based on genomic DNA variation and for some common diseases, including coronary artery disease and atrial fibrillation, PRSs reached to a point when it is appropriate to implement into clinical care (Khera et al., 2018). In addition to PRS studies that are designed to identify high-risk groups recently there is also an intriguing discussion about the potential value for the use in evaluating therapeutic response (Gibson, 2019). Nevertheless, the same reviews are also outlining that when PR scores identify the upper two deciles of the risk distribution it is on the cost of declining sensitivity. PRS sensitivity (defined as the approximate proportion of cases captured by the PRS at a given PRS percentile) usually drops below 50% in such typical cases as estimated from a PRS created from many thousands of variants in a population of hundreds of thousands of individuals, such as the UK Biobank (Gibson, 2019).
GWAS Findings as “Clinically Crude Data”
From the clinician's perspective the genetic findings of GWASs might be interpreted as “crude” data. A refined disease subgroup definition and a more precise phenotyping might be applied to get closer to a more meaningful clinical association (Ahlqvist et al., 2018; Meigs, 2019). We highlight that in case of a few—already identified disease associated—gene variants the genetic effect sizes may be substantially higher provided that we assess more precisely defined clinical traits and treatment modalities. This approach that might potentially be applicable only for a few gene variants, but could help more in the clinical use of genetic information.
Association of MTNR1B rs10830963 Gene Variant and Diabetes Mellitus Traits
Type 2 Diabetes Mellitus (T2DM) Development
The MTNR1B rs10830963 associated effects on non-autoimmune diabetes mellitus (DM) traits are discussed as an example, more specifically the differences in the genetic effect sizes for T2DM development and GDM development/therapy.
The odds ratios (ORs)—referring to the genetic effect sizes—associated to the rs10830963 variant are very moderate for T2DM development (1.09–1.12-fold per G allele) (Lyssenko et al., 2009; Prokopenko et al., 2009).
Gestational Diabetes Mellitus (GDM) Development
In contrast, independent case-control genetic studies report OR values between 1.29 and 2. for GDM development in Middle European (Vejrazkova et al., 2014; Rosta et al., 2017), Mediterranean Hispanic (García de la Torre et al., 2019), Middle Eastern (Alharbi et al., 2019) Asian (Wu et al., 2016) and Russian (Popova et al., 2017) pregnant populations. The MTNR1B rs10830963 associated effects appeared to be more preponderant provided that the mean pre-pregnancy BMI was above 25 kg/m2 (Wu et al., 2016). Furthermore the OR values for GDM development were particularly high (1.96) when BMI was above 30 kg/m2 (Alharbi et al., 2019) even after adjustment to BMI. This indicates substantially higher MTNR1B rs10830963 associated genetic effect sizes when GDM compared to T2DM development and therefore suggests additional pregnancy specific factors.
Association of MTNR1B rs10830963 Gene Variant and Prevention and Therapy of GDM
We found that the genetic effect sizes are further increased (OR = 5.2 for G allele carriers) on the need for antenatal insulin therapy (AIT) initiation in Hungarian GDM patients (Firneisz et al., 2018) above the pre-pregnancy BMI threshold of 29 kg/m2. Furthermore, it was independently reported (Grotenfelt et al., 2016) that the rs10830963 G risk allele carrying significantly decreased the odds with a similar effect size (OR for intervention failure in G allele carriers >5) of responding to an early (from the 13th gestational week) medical nutrition therapy and lifestyle intervention in high-risk (pre-pregnancy BMI ≥ 30 kg/m2 and/or prior GDM) Finnish pregnant women. In addition the MTNR1B rs10830963 variant is reported to augment the effect of lifestyle risk factors and dietary habits (in particular sausage consumption) on GDM development (Popova et al., 2017).
These three reports independently indicate clinical consequences associated to the MTNR1B rs10830963 risk genotype in the prevention strategies and also the treatment of high risk disease subpopulations.
Difference in Genetic Effect Sizes and Potential Explanations
The high effect sizes in precisely identified patient sub-populations might be potentially useful for guiding the clinical practice, although it might be relevant for only around 30% of the entire disease population in T2DM (Florez, 2019). A novel approach of associating clinical utility to the crude GWAS findings and the MTNR1B rs10830963 associated genetic effects to distinct endophenotypes within the diabetes mellitus spectrum might demonstrate consistently similar features (Figure 1).
Figure 1. The association of MTNR1B rs10830963 gene variant with different clinical traits in non-autoimmune diabetes mellitus and the corresponding odds ratios.
One possible pregnancy specific factor to explain this phenomenon might be the robust increase (>3x at the end of pregnancy compared to levels before the 24th gestational week) in the night time maternal serum melatonin (physiological ligand of the MTNR1B encoded MT2 receptor) concentrations (Nakamura et al., 2001). The continuous rise of the night time maternal serum melatonin levels from the 24 to 28th gestational week accords with the time-frame of the routine diagnostic OGTT during pregnancy. The expression of MTNR1B mRNA is generally very low in all tissues, yet a higher expression (HPA RNA-seq normal tissues) is found in the placenta where it is also expressed at protein level (MT2) (Li et al., 2018). One may hypothesize that the exceptionally high maternal ligand concentrations during the 3rd trimester are somehow causally related to the multiplied effect sizes associated to the rs1083096 receptor gene variant in GDM development and treatment.
In addition, the difference in age ranges of the study participants in T2DM GWAS (35–70 years of age) vs. GDM genetic association studies (substantially younger, reproductive age) might also contribute to the attenuated MTNR1B rs10830963 associated genetic effect sizes reported for T2DM compared to GDM development. This could potentially be due to the age-associated decline in melatonin production and melatonin receptor expression (Hill et al., 2011) and the ab ovo disrupted central clock driven circadian melatonin output in metabolic diseases in the elderly (Belancio et al., 2014).
If we systematically narrow the studied population to precisely defined subgroups (T2DM development < GDM development < AIT initiation in GDM with BMI over 29 kg/m2), we experience a step-wise increase in the genetic effect size (Figure 1).
The Potential Significance of MTNR1B rs10830963 in GDM Complications
Provided that the observation is approached reversely, from the point of genetic associations to the complications in populations with high enough proportion of undiagnosed/untreated GDM patients (“natural disease course”) then we should—in theory—arrive back to the MTNR1B rs10830963 variant. Recently, a large GWAS reported maternal genetic effects for the neonatal birth weight that were independent of offspring genetics (Beaumont et al., 2018). Among the identified maternal signals the MTNR1B rs10830963 was characterized by the lowest p-value association with neonatal birthweight and had the highest genetic effect size out of the more than 8 million gene variants assessed in the united EGG and UKBB meta-analysis. It is meaningful that this study population was unadjusted to glycemia during pregnancy/presence of GDM (which was not systematically excluded using an OGTT during pregnancy), despite these maternal conditions are well-known to be associated with offspring birth weight. We should hypothesize that this birth weight increasing genetic effect also occurs when the neonatal outcome is expressed in birth weight percentiles and therefore the maternal MTNR1B rs10830963 might be associated with higher occurrence of large for gestational age (LGA) one of the most known GDM complications.
Discussion
The wider clinical use of genetic data in complex diseases could be determined by
• The risk allele frequency (RAF) in the population,
• The genetic effect size related to a clinically significant outcome
• And the appropriate definition of the sub-population (“endophenotype”) in which the prediction of the therapeutic response is aimed.
The proportion of MTNR1B rs10830963 risk G allele carriers reached/exceeded the 50% in the populations reported on Figure 1 that might indicate a rationale for population based screening programs of high risk GDM pregnancies in these regions. Our response as clinicians should already be triggered in any case when a meaningful outcome, such as antenatal insulin treatment in GDM is assessed in a precisely defined clinical trait and the exposure results in an OR value over 5, even in the case of a single risk gene variant without the use of a challenging PRS/PHS. AIT need in GDM might still be a complex trait from the genetic point, however we think that—unlike T2DM development—this genetic complexity could be characterized by a gene variant with the major effect (MTNR1B rs10830963) in populations where it occurs with high enough MAF (such as European, Asian and Middle-Eastern origins).
Interestingly, in parallel with identifying the upper two deciles in the risk distribution in this case this model employing a single, but high effect gene variant could still keep the pick-up rate (equivalents to the sensitivity) above 50%, in contrast to PRS which—in general—still faces the problem of declining sensitivity (defined as the approximate proportion of cases captured by the PRS at a given percentile) as the cost that the higher precision and OR values could only be gained with increased number of gene variants (Gibson, 2019).
We may conclude that an approach employing even a single gene variant might potentially be more informative in predicting the therapeutic responses provided that this model is used in a narrow, clinically well-defined “endophenotype” in patients already diagnosed with a prevalent disease if the genetic effect size is high enough and the variant is not rare.
However, we suspect that such clinical conditions might not be that frequent and therefore the described approach could be considered as auxiliary to the PRS which appears to be more of clinical use (even if it is currently limited) when the aim is to assess the genetic risk in largely different, broader, more heterogeneous populations with individually weak effect gene variants.
Three independent genetic studies confirmed that the maternal MTNR1B rs10830963 risk variant is consistently associated with major forms of clinical consequences (Florez, 2019), such as
• Early prevention failure in high risk pregnant (Grotenfelt et al., 2016)
• Therapy (AIT need) in high risk GDM (Firneisz et al., 2018) subpopulations and
• A neonatal complication related trait in the general population (Beaumont et al., 2018).
Therefore we may not only hypothesize that the precision medicine approach (employing the MTNR1B rs10830963 gene variant) in predicting the AIT need in a narrow subpopulation of patients already diagnosed with GDM might result in better and earlier identification of patients who would require antenatal insulin therapy (Firneisz et al., 2018), but also—consistently with the findings of the maternal GWAS on neonatal birthweight (Beaumont et al., 2018)—we may potentially expect the reduction of complications (e.g., macrosomia/LGA rate).
Nevertheless, as an ultimate gold standard first step in clinics an RCT should be conducted to confirm that this genetic precision medicine based approach might result in any of the expected clinical benefits at all.
Author Contributions
GF conceived the opinion and wrote the manuscript. KR, JR, JH, AK-W, AS, and ÁN reviewed the manuscript.
Funding
The publication costs were covered by the Semmelweis University, Budapest, Hungary.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Ahlqvist, E., Storm, P., Karajamaki, A., Martinell, M., Dorkhan, M., Carlsson, A., et al. (2018). Novel subgroups of adult-onset diabetes and their association with outcomes: a data-driven cluster analysis of six variables. Lancet Diabetes Endocrinol. 6, 361–369. doi: 10.1016/S2213-8587(18)30051-2
Alharbi, K. K., Al-Sulaiman, A. M., Shedaid, K. M. B., Al-Shangiti, A. M., Marie, M., Al-Sheikh, Y. A., et al. (2019). MTNR1B genetic polymorphisms as risk factors for gestational diabetes mellitus: a case-control study in a single tertiary care center. Ann. Saudi Med. 39, 309–318. doi: 10.5144/0256-4947.2019.309
Beaumont, R. N., Warrington, N. M., Cavadino, A., Tyrrell, J., Nodzenski, M., Horikoshi, M., et al. (2018). Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics. Hum. Mol. Genet. 27, 742–756. doi: 10.1093/hmg/ddx429
Belancio, V. P., Blask, D. E., Deininger, P., Hill, S. M., and Jazwinski, S. M. (2014). The aging clock and circadian control of metabolism and genome stability. Front. Genet. 5:455. doi: 10.3389/fgene.2014.00455
Davies, M. J., D'Alessio, D. A., Fradkin, J., Kernan, W. N., Mathieu, C., Mingrone, G., et al. (2018). Management of hyperglycaemia in type 2 diabetes, 2018. A consensus report by the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD). Diabetologia 61, 2461–2498. doi: 10.1007/s00125-018-4729-5
Firneisz, G., Rosta, K., Al-Aissa, Z., Hadarits, O., Harreiter, J., Nadasdi, A., et al. (2018). The MTNR1B rs10830963 variant in interaction with pre-pregnancy BMI is a pharmacogenetic marker for the initiation of antenatal insulin therapy in gestational diabetes mellitus. Int. J. Mol. Sci. 19, 3734–3747. doi: 10.3390/ijms19123734
Florez, J. C. (2019). “A global overview of precision medicine in type 2 diabetes,” in EASD Virtual Meeting (Barcelona).
García de la Torre, N., Valerio, J., Barabash, A., del Valle, L., Assaf-Balut, C., Bordiu, E., et al. (2019). “A genetic risk score combining four risk genes for diabetes modulate the effect of a mediterranean diet on the incidence of gestational diabetes”, in EASD Virtual Meeting (online presentation) (Barcelona).
Gibson, G. (2019). On the utilization of polygenic risk scores for therapeutic targeting. PLoS Genet. 15:e1008060. doi: 10.1371/journal.pgen.1008060
Grotenfelt, N. E., Wasenius, N. S., Rono, K., Laivuori, H., Stach-Lempinen, B., Orho-Melander, M., et al. (2016). Interaction between rs10830963 polymorphism in MTNR1B and lifestyle intervention on occurrence of gestational diabetes. Diabetologia 59, 1655–1658. doi: 10.1007/s00125-016-3989-1
Hill, S. M., Cheng, C., Yuan, L., Mao, L., Jockers, R., Dauchy, B., et al. (2011). Declining melatonin levels and MT1 receptor expression in aging rats is associated with enhanced mammary tumor growth and decreased sensitivity to melatonin. Breast Cancer Res. Treat 127, 91–98. doi: 10.1007/s10549-010-0958-0
Khera, A. V., Chaffin, M., Aragam, K. G., Haas, M. E., Roselli, C., Choi, S. H., et al. (2018). Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations. Nat. Genet. 50, 1219–1224. doi: 10.1038/s41588-018-0183-z
Lambert, S. A., Abraham, G., and Inouye, M. (2019). Towards clinical utility of polygenic risk scores. Hum. Mol. Genet. 28, R133–R142. doi: 10.1093/hmg/ddz187
Li, C., Zhou, Y., Qiao, B., Xu, L., Li, Y., and Li, C. (2018). Association between a melatonin receptor 1B genetic polymorphism and its protein expression in gestational diabetes mellitus. Reprod. Sci. 26, 1382–1388. doi: 10.1177/1933719118765983
Loscalzo, J. (2019). Network medicine and type 2 diabetes mellitus: insights into disease mechanism and guide to precision medicine. Endocrine 66, 456–459. doi: 10.1007/s12020-019-02042-4
Lyssenko, V., and Laakso, M. (2013). Genetic screening for the risk of type 2 diabetes: worthless or valuable? Diabetes Care 36 (Suppl. 2), S120–126. doi: 10.2337/dcS13-2009
Lyssenko, V., Nagorny, C. L., Erdos, M. R., Wierup, N., Jonsson, A., Spegel, P., et al. (2009). Common variant in MTNR1B associated with increased risk of type 2 diabetes and impaired early insulin secretion. Nat Genet. 41, 82–88. doi: 10.1038/ng.288
Meigs, J. B. (2019). The genetic epidemiology of type 2 diabetes: opportunities for health translation. Curr. Diab. Rep. 19:62. doi: 10.1007/s11892-019-1173-y
Meigs, J. B., Shrader, P., Sullivan, L. M., McAteer, J. B., Fox, C. S., Dupuis, J., et al. (2008). Genotype score in addition to common risk factors for prediction of type 2 diabetes. N. Engl. J. Med. 359, 2208–2219. doi: 10.1056/NEJMoa0804742
Nakamura, Y., Tamura, H., Kashida, S., Takayama, H., Yamagata, Y., Karube, A., et al. (2001). Changes of serum melatonin level and its relationship to feto-placental unit during pregnancy. J. Pineal Res. 30, 29–33. doi: 10.1034/j.1600-079x.2001.300104.x
Popova, P. V., Klyushina, A. A., Vasilyeva, L. B., Tkachuk, A. S., Bolotko, Y. A., Gerasimov, A. S., et al. (2017). Effect of gene-lifestyle interaction on gestational diabetes risk. Oncotarget 8, 112024–112035. doi: 10.18632/oncotarget.22999
Prasad, R. B., and Groop, L. (2015). Genetics of type 2 diabetes-pitfalls and possibilities. Genes 6, 87–123. doi: 10.3390/genes6010087
Prokopenko, I., Langenberg, C., Florez, J. C., Saxena, R., Soranzo, N., Thorleifsson, G., et al. (2009). Variants in MTNR1B influence fasting glucose levels. Nat. Genet. 41, 77–81. doi: 10.1038/ng.290
Rosta, K., Al-Aissa, Z., Hadarits, O., Harreiter, J., Nadasdi, A., Kelemen, F., et al. (2017). Association Study with 77 SNPs confirms the robust role for the rs10830963/G of MTNR1B variant and identifies two novel associations in gestational diabetes mellitus development. PLoS ONE 12:e0169781. doi: 10.1371/journal.pone.0169781
Sarraju, A., and Knowles, J. W. (2019). Genetic testing and risk scores: impact on familial hypercholesterolemia. Front. Cardiovasc. Med. 6:5. doi: 10.3389/fcvm.2019.00005
Srinivasan, S., Jablonski, K. A., Knowler, W. C., Dagogo-Jack, S., Kahn, S. E., Boyko, E. J., et al. (2019). A polygenic lipodystrophy genetic risk score characterizes risk independent of BMI in the diabetes prevention program. J. Endocr. Soc. 3, 1663–1677. doi: 10.1210/js.2019-00069
Tam, V., Patel, N., Turcotte, M., Bosse, Y., Pare, G., and Meyre, D. (2019). Benefits and limitations of genome-wide association studies. Nat. Rev. Genet. doi: 10.1038/s41576-019-0127-1
Vejrazkova, D., Lukasova, P., Vankova, M., Vcelak, J., Bradnova, O., Cirmanova, V., et al. (2014). MTNR1B genetic variability is associated with gestational diabetes in czech women. Int. J. Endocrinol. 2014:508923. doi: 10.1155/2014/508923
Keywords: genotype, GDM, therapy, prediction, precision medicine, BMI, LGA
Citation: Firneisz G, Rosta K, Rigó J, Nádasdi Á, Harreiter J, Kautzky-Willer A and Somogyi A (2020) Identification and Potential Clinical Utility of the MTNR1B rs10830963 Core Gene Variant Associated to Endophenotypes in Gestational Diabetes Mellitus. Front. Genet. 11:332. doi: 10.3389/fgene.2020.00332
Received: 16 August 2019; Accepted: 20 March 2020;
Published: 21 April 2020.
Edited by:
Ahmed Rebai, Centre of Biotechnology of Sfax, TunisiaReviewed by:
Luisella Bocchio Chiavetto, Centro San Giovanni di Dio Fatebenefratelli (IRCCS), ItalyYosr Hamdi, Institut Pasteur de Tunis, Tunisia
Lilia Romdhane, University of Carthage, Tunisia
Copyright © 2020 Firneisz, Rosta, Rigó, Nádasdi, Harreiter, Kautzky-Willer and Somogyi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gábor Firneisz, ZmlybmVpc3ouZ2Fib3ImI3gwMDA0MDttZWQuc2VtbWVsd2Vpcy11bml2Lmh1