Scroll to:
Research of association of the polymorphic locus rs11212617 ATM gene with the response to therapy with metformin in patients with type 2 diabetes
https://doi.org/10.14341/probl20176319-16
Full Text:
Abstract
Rationale. The genetic aspects influencing the effectiveness of metformin (MF) therapy in patients with type 2 diabetes mellitus (DMT2) have recently been intensively studied.
Objective — to study the association between the rs11212617 polymorphism in the ATM gene and response to metformin therapy in DMT2 patients in the Novosibirsk region and to conduct a metaanalysis of the previously reported data.
Material and methods. 460 DMT2 patients (97 males and 363 females) who received MF, both as a part of monotherapy and in combination with sulfonylurea (SU) drugs, were subjected to cross-sectional examination. Depending on HbA1c level, patients were divided into the following groups: patients who have attained the target HbA1c level after MF therapy (n=209) and those who did not attain the target HbA1c level although receiving the maximum dose of MF (n=251). Alleles and genotypes were determined by real-time PCR using TaqMan probes at the Institute of Chemical Biology and Fundamental Medicine (SB RAS).
Results. Frequency of the rare C allele of the rs11212617 polymorphism in the ATM gene in the examined patients was 0.41 and statistically did not differ between the subgroups who received MF monotherapy and combination therapy. Statistically significant association between the genotype of the rs11212617 polymorphism in the ATM gene and the type of response was revealed neither in the total group of patients (OR=0.94, 95% CI 0.73—1.23; p=0.67) nor in the MF monotherapy (OR=0.94, 95% CI 0.73—1.23; p=0.67) or combination therapy subgroups (OR=1.02, 95% CI 0.72—1.43; p=0.92). However, the metaanalysis results verify that the C allele is associated with attainment of the target HbAc1 level (the total OR=1.27, 95% CI 1.10—1.46; p=0.0008).
Conclusions. The rs11212617 polymorphism in the ATM gene can influence the effectiveness of MF therapy in DMT2 patients.
For citations:
Bondar I.A., Shabelnikova O.Y., Sokolova E.A., Filipenko M.L. Research of association of the polymorphic locus rs11212617 ATM gene with the response to therapy with metformin in patients with type 2 diabetes. Problems of Endocrinology. 2017;63(1):9-16. https://doi.org/10.14341/probl20176319-16
BACKGROUND
Currently, metformin (MF) is the fi rst-line medication in the treatment of type 2 diabetes mellitus (T2DM) [1, 2]. It is known that the eff ectiveness of MF therapy pend on renal function, enzyme systems activity, and drug interactions, leading to alteration of MF pharmacokinetics and pharmacodynamics. MF is used in clinical medicine for more than 50 years. However, the mechanisms underlying its therapeutic eff ects are not fully understood.
After entering the cell, MF is accumulated in the mitochondrial matrix. Mitochondrial respiratory chain complex I is a key target of MF, whose inhibition results is decreased production of ATP and increased levels of AMP and ADP. This change in cellular energy metabolism is detected by AMP-activated protein kinase (AMPK), the main cellular energy sensor. MF enhances phosphorylation processes, which leads to increased cellular level of AMPK followed by with inhibition of mTOR (mammalian target of rapamycin) protein kinase, which affects metabolism and growth factor signaling, diff erentiation, aging, apoptosis, and autophagic processes [3, 4]. Activation of AMPK reduces activity of cell cycle enzymes, gluconeogenesis, and glycogenolysis; stimulates anaerobic glycolysis, fatty acid oxidation in liver, and ketogenesis; inhibits lipogenesis, synthesis of cholesterol and triglycerides, inhibits lipolysis and lipogenesis in adipocytes; stimulates fatty acid oxidation and glucose uptake in skeletal muscles through enhanced biosynthesis of GLUT-1 and GLUT-4 glucose transporters; modulates insulin secretion by beta cells [5]. It is assumed that decreased expression and activity of AMPK is associated with obesity, metabolic syndrome, disturbance of carbohydrate metabolism, cancer, and myocardial ischemia [2, 5]. In recent years, genetic aspects aff ecting the efficacy of MF therapy are being extensively studied. The study of the effect of MF on AMPK, the product of ataxiatelangiectasia (ATM) gene expression, is of greatest interest [6]. According to the literature, the minor С allele at the polymorphic locus rs11212617 in the ATM gene is associated with good response to MF [7].
AIM
The study was aimed at assessing association of the polymorphic locus rs11212617 in the ATM gene with the response to MF therapy in T2DM patients in the Novosibirsk region and conducting metaanalysis of previously published data.
METHODS
Study Design
A single-step cross-sectional survey of 2000 T2DM patients (450 males and 1550 females) was carried out on the basis of the mobile diabetes center (Diamobil) in the Novosibirsk region.
Inclusion criteria
A total of 2000 patients were examined and 460 of them were included in the study (97 males and 363 females), who received MF either as monotherapy or in combination with sulfonylurea (SU) medications. Exclusion criteria were as follows: type 1 diabetes mellitus (T1DM), other types of diabetes, T2DM with insulin therapy, cancer, NYHA III—IV heart failure, hepatic insufficiency (more than 7 points on the Child-Pugh scale), decreased glomerular fi ltration rate (GFR) of less than 45 mL/min × 1.73 m2, treatment with corticosteroids or estrogens, alcoholism, drug addiction, dementia, or severe mental disorders.
Study conditions
The study included only patients, who provided written informed consent.
Duration of the study
From 2012 to 2015.
Description of medical intervention
All 460 patients underwent complete clinical examination; the level of glycated hemoglobin (HbA1c) was assessed using immuno-turbidimetric method on BIORAD D10 analyzer (US). DNA isolation and genotyping of rs11212617 polymorphic locus of the ATM gene was carried out in the laboratory of Pharmacogenomics at the Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences (SB RAS) using TaqMan real-time PCR. PCR was performed in a fi nal volume of 25 μl, containing 65 mM Tris-HCl (pH 8.9), 24 mM ammonium sulfate; 3.5 mM MgCl2; 0,05% Tween-20; 300 nM of each primer; 100 nM of each Taq Man-probe conjugated to FAM or R6G; 200 μM of dNTP, 20—100 ng of DNA, and 1 activity unit of Taq-DNA polymerase. Amplifi cation was carried out using CFX384 amplifi er (Bio-Rad, USA) under the following conditions: initial 3’ denaturation at 96 °C followed by 50 cycles, including denaturation at 96°C during 8 sec, primer annealing and subsequent elongation at Tann = 60 °C for 40 seconds (each step was followed by registration of the fl uorescence signal in the range corresponding to fl uorescence intervals of FAM and R6G fl uorophores). Genotype of the polymorphic locusrs11212617 was determined using the following oligonucleotides: forward primer 5’-GATCTACATATACCAATTACAAAGG-3’, reverse primer 5’GGATAACATATAGTTGGGTCTTG-3’, probe for A allele5’-FAM-CAGAGAATGTCAGAGCGG-BHQ-3’, probe for C allele 5’-R6G-CAGAGACTGTCAGAGCGG-BHQ-3’.
The main result of the study
Association of the polymorphic locus rs11212617 in the ATM gene with the response to MF therapy in T2DM patients in the Novosibirsk region was evaluated based on HbA1c level. Patients who have attained target HbA1c values in accordance with the algorithms of specialized medical care for diabetic patients (2015) were included inthe group with good response to MF. Patients who have not attained target HbA1c values were included in the group with poor response to MF.
Subgroup analysis
Depending on the level of HbA1c, patients were divided into groups: patients with target HbA1c level during MF therapy (the group with good response to MF administered both as monotherapy and in combination with SU medications) and the group of patients who have not attained the target HbA1c level, while receiving the maximum dose of MF.
Ethical review
Research protocol was approved by the ethics committee of the Novosibirsk State Medical University (Protocol No 52 of 19.03.13). All patients signed an informed consent before entering the study.
Statistical analysis
We used standard descriptive and comparative analysis. In the case of normal distribution, the data were represented as a mean value (M) and standard deviation (SD), in the case of asymmetric distribution, the data were represented as a median (Me) with 25 and 75 percentiles (25; 75). Pearson’s χ2 test was used to compare the frequency of alleles between the groups. Depending on the type of data distribution, analysis of variance (ANOVA) or Kruskal—Wallis test was used to evaluate the diff erences between the groups. The critical significance level was set at 0.05. Compliance with Hardy— Weinberg equilibrium was assessed using Fisher’s exact test. Statistical analysis was performed using Genetics statistical package as implemented in R-project software (www.r-project.org) and Statistica 6.0 program. Metaanalysis, evaluation of the ≪archival box≫ effect (publication bias), and construction of ≪forest plots≫ and ≪funnel plots≫ were carried out using meta and rmeta packages as implemented in R-project software (www.r-project.org). PubMed database (http://www.ncbi.nlm.nih.gov/pubmed/) was searched for articles published up to and including May 12, 2016 using various combinations of ≪rs11212617≫, ≪T2DM≫, and ≪diabetes mellitus≫ keywords. Inclusion criteria were as follows: 1) the design of the study ≪case-control≫, 2) evaluation of the association of rs11212617 with attainment of target HbA1c level of <7%, 3) presentation of the results in English, 4) availability of all information required to calculate OR. Metaanalysis included the results of two studies, Zhou [8, 11] and van Leeuwen [9], along with our results. We did not include the fi ndings of Zhou (2014, RS11212617 is associated with metformin treatment response in type 2 diabetes in Shanghai local Chinese population) obtained for Chinese population, since the article is not available and the summary does not contain the required information. A total of 5494 T2DM patients, who received MF, were included in metaanalysis.
RESULTS
Objects (participants) of the study
The average age of the examined patients was 58.8Ѓ} 8.6 years, duration of T2DM — 6.2Ѓ} 5.7 years, HbA1c level — 8.4Ѓ} 1.9%, BMI — 34.8Ѓ} 6.7 kg/m2. The patients were divided into groups depending on HbA1c level: the group with target HbA1c level during MF therapy (the group with good response to MF (n=209), of whom 110 patients received MF monotherapy and 99 patients received a combination of the MF with SU; and the group of patients, who have not attained target HbA1c level, while receiving the maximum dose of MF (the group with poor response to MF (n=251), of whom 58 patients received MF monotherapy and 193 patients received a combination of MF with SU. The analysis of clinical data showed that patients in the group with poor response to MF were younger as compared to the group with good response to MF (57.5Ѓ}8.1 vs. 60.4Ѓ}8.9 years; p<0.05), had longer duration of the disease (7.2Ѓ}6.0 vs. 5.1Ѓ}5.1 years; p<0.05), higher HbA1c levels (9.8Ѓ}1.6% vs. 6.7Ѓ}0.4%; p<0.05), earlier onset of the disease (50.2Ѓ}8.6 vs. 55.4Ѓ}9.6 years; p<0.05) and did not diff er signifi cantly in BMI, degree of insulin resistance, glomerular fi ltration rate (GFR), and transaminase levels (Table 1).
Table 1. Clinical and laboratory characteristics of the examined patients
Key findings
The distribution of genotypes of polymorphic locus rs11212617 in the ATM gene demonstrated no statistical deviations from Hardy-Weinberg equilibrium (p=0.21). In the examined patients, the frequency of the minor C allele at the polymorphic locus rs11212617 in the ATM gene was 0.41; in the subgroup of patients with poor response to MF monotherapy, it reached 0.46, and in the subgroup with good response it was 0.39 (p>0.05); in the overall group of patients with poor response to MF in combination with SU, the frequency of C allele was the same, 0.4 and 0.41, respectively and did not diff er signifi cantly from that in the group with good response (Figure 1). Logistic regression analysis of the relationship between the attainment of target HbA1c level of <7% during MF treatment and genotypes of the polymorphic locus rs11212617 in the ATM gene demonstrated no statistically signifi cant association of rs11212617 genotype with the type of response both in the overall group of patients (OR=0.94, 95% CI 0.73—1.23; p=0.67) and in subgroups of patients, who received MF monotherapy (OR=0.90, 95% CI 0.65—1.25; p=0.54) and combination therapy (OR=1.02, 95% CI 0.72—1.43; p = 0.92) (Table 2).
Figure 1. Frequency of minor C allele at the rs11212617 polymorphic locus in the ATM gene in the groups with different responses to metformin therapy.
Table 2. Logistic regression analysis of the relationship between the attainment of target HbA1с level of <7% during metformin therapy and genotypes of polymorphic locus rs11212617 in the ATM gene
Additional findings
We carried out meta-analysis of the results of our investigation and previously published studies. This metaanalysis included a total of 5494 T2DM patients, who received MF (Table 3).
Table 3. Metaanalysis of the studies of the relationship between the attainment of target HbA1с level of <7% during MF therapy and the genotype of polymorphic locus rs11212617 in the ATM gene
According to the results of meta-analysis, minor C allele is a factor associated with good response to MF therapy (overall OR=1.27, 95% CI 1.10—1.46;p=0.0008, the level of heterogeneity p=0.02, the random eff ects model). We found no bias in the results published by other researchers, as shown by the correlation analysis (Begg’s correlation analysis: z = –0.45; p=0.65) and regression test (Egger’s test: t=–0.470; p=0.66). In addition, we carried out metaanalysis in both of the MF monotherapy and combination therapy subgroups. However, while the results in the MF monotherapy subgroup were similar [overall OR = 1.26, 95% CI 1.12—1.42; p=0.0001, the level of heterogeneity p=0.15, fi xed-eff ects model, Begg’s correlation analysis (z=1.32 and p=0.19), Egger’s test (t=0.79; p=0.48), no impact of the polymorphic locus on the attainment of target HbA1c level was observed in the group with combination therapy: overall OR=1.13, 95% CI 0.98—1.29; p=0.09, the level of heterogeneity p=0.75, fi xed-eff ects model, Begg’s correlation analysis (z=–0.56 and p=0.57), Egger’s test 0.73). The results of metaanalysis and bias analysis of publications are plotted in Figure 2.
Figure 2. Graphical representation of the results of meta-analysis and bias analysis in publications.
A — ≪forest plot≫ for rs11212617 association with attainment of target HbA1c level of <7% in the overall group; B — ≪forest plot≫ for rs11212617 association with attainment of target HbA1c level of <7% in the group with metformin monotherapy, C — ≪forest plot≫ for rs11212617 association with attainment of target HbA1c level of <7% in the group with combination therapy, D — ≪funnel plot≫ for rs11212617 association with attainment of target HbA1c level of <7% in the overall group; E — ≪funnel plot≫ for rs11212617 association with attainment of target HbA1c level of <7% in the group with metformin monotherapy, F — ≪funnel plot≫ for rs11212617 association with attainment of target HbA1c level of <7% in the group with combination therapy.
Adverse events
No adverse events were registered during this study.
DISCUSSION
The fi rst pharmacogenetic genome-wide associative study of MF was carried out in Scotland in 1024 T2DM patients involved in the GoDARTS (Genetics of Diabetes Audit and Research Tayside) study [8]. Effi cacy of MF therapy was assessed based on HbA1c levels during the first 18 months from the start of therapy. The target HbA1c level was <7%, which was considered as an evidence of the eff ectiveness of hypoglycemic therapy. The statistical analysis considered the initial level of HbA1c, its dynamics during MF therapy, and creatinine clearance. Researchers have identifi ed 14 polymorphisms in the region, comprising ATM gene. It was found that the carriers of minor C allele at the polymorphic locus rs11212617 in the ATM gene attained target HbA1c level more frequently (OR= 1.64, 95% CI 1.37—1.99; p=1.9× 10–7). The authors genotyped the polymorphic locus in two independent groups consisting of 1783 Scottish T2DM patients (the other part of the GoDARTS study) and 1113 patients involved in the UKPDS study (UK Prospective Diabetes Study) in order to replicate the original results of the genome-wide study. Signifi cant association with response to MF therapy was found in both populations. The cumulative eff ect in attaining the objectives of hypoglycemic MF therapy was statistically signifi cant at the genome-wide level (OR=1.35, 95% CI 1.22—1.49; p=2.9 × 10–7). Later on, the same group of researchers have determined the genotypes of the polymorphic locus rs11212617 in three more cohorts: DCS, Rotterdam Study, and CARDS [8, 9]. Subsequent metaanalysis in fi ve cohorts under study (the fi rst onewas excluded in order to avoid biased results) included 4443 T2DM patients and it also confi rmed the relationship between the minor C allele and attainment of target HbA1c level (OR= 1.25, 95% CI 1.13—1.38; p=7.8 × 10–6). It is noteworthy that, when dividing the group into MF monotherapy and combination therapy subgroups, association of C allele was preserved only in the monotherapy subgroup. Van Leeuwen et al. hypothesized that some eff ects of MF therapy (particularly, AMPK-mediated ones, may partially be neutralized due to SU medications intake [10]. Possibly, this gene portion controls the cell cycle, and its changes can lead to decrease in phosphorylation and activation of AMPK, and therefore to diff erent hypoglycemic eff ect of MF[10]. Other researchers, who studied the eff ect of the polymorphic locus rs11212617 in the ATM gene on the likelihood of developing T2DM in patients with obesity and impaired glucose tolerance (IGT) involved in the Diabetes Prevention Program, DPP, including representatives of five American ethnic groups (n=988), suggested that the time of IGT conversion to T2DM in patients, who receive MF and have diff erent genotypes of the polymorphic locus rs11212617, may diff er. However, no signifi cant difference was found. Florez et al. suggested that the effect of the polymorphic locus during MF treatment is noticeable only in the case of high HbA1c, which is possible only in patients who were already diagnosed with T2DM. In the patients involved in their study, T2DM developed later and HbA1c levels were low [11]. Recently, Vilvanathan et al. evaluated association of polymorphic locus rs11212617 with the development of T2DM in the residents of South India (118 T2DM patients and 112 healthy controls). However, this polymorphiclocus had no effect on predisposition to T2DM: OR=1.09, 95% CI 0.75—1.60 [12].
CONCLUSION
Our study showed no association of polymorphic locus rs11212617 in the ATM gene with response to MF therapy in T2DM patients in the Novosibirsk region. However, meta-analysis of our own data with previously published reports confi rms the association of C allele with attainment of target HbA1c level (overall OR=1.27, 95% CI 1.10—1.46; p=0.0008). Thus, polymorphic locus rs11212617 in the ATM gene may infl uence the eff ectiveness of MF therapy in T2DM patients.
ADDITIONAL INFORMATION
Financial support: The study was supported by Russian Foundation for Basic Research grant No 13-04-00520.
Conflict of interest. The authors declare no confl ict (duality) of interests during writing this article.
References
1. Inzucchi SE, Bergenstal RM, Buse JB, et al. Management of hyperglycemia in type 2 diabetes, 2015: a patient-centered approach: update to a position statement of the American Diabetes Association and the European Association for the study of diabetes. Diabetes Care. 2015;38(1):140-149. doi: 10.2337/dc14-2441
2. Дедов И.И., Шестакова М.В., Галстян Г.Р., и др. Алгоритмы специализированной медицинской помощи больным сахарным диабетом. Под ред. Дедова И.И., Шестаковой М.В. (7-й выпуск). // Сахарный диабет. — 2015. — Т. 18. — №1S. — C. 1-112. [Dedov II, Shestakova MV, Galstyan GR, et al. Standards of specialized diabetes care. Edited by Dedov II, Shestakova MV. (7th edition). Diabetes mellitus. 2015;18(1S):1-112. (in Russ.)]. doi: 10.14341/dm20151s1-112
3. Johnson AB, Webster JM, Sum CF, et al. The impact of metformin therapy on hepatic glucose production and skeletal muscle glycogen synthase activity in overweight type ii diabetic patients. Metabolism. 1993;42(9):1217-1222. doi: 10.1016/0026-0495(93)90284-u
4. Zhou G, Myers R, Li Y, et al. Role of AMP-activated protein kinase in mechanism of metformin action. J Clin Invest. 2001;108(8):1167-1174. doi: 10.1172/jci13505
5. Kickstein E, Krauss S, Thornhill P, et al. Biguanide metformin acts on tau phosphorylation via mTOR/protein phosphatase 2A (PP2A) signaling. Proc Natl Acad Sci USA. 2010;107(50):21830-21835. doi: 10.1073/pnas.0912793107
6. Hawley SA, Ross FA, Chevtzoff C, et al. Use of cells expressing gamma subunit variants to identify diverse mechanisms of AMPK activation. Cell Metab. 2010;11(6):554-565. doi: 10.1016/j.cmet.2010.04.001
7. Foretz M, Hebrard S, Leclerc J, et al. Metformin inhibits hepatic gluconeogenesis in mice independently of the LKB1/AMPK pathway via a decrease in hepatic energy state. J Clin Invest. 2010;120(7):2355-2369. doi: 10.1172/jcI40671
8. GoDarts, Group UDPS, Wellcome Trust Case Control C, et al. Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes. Nat Genet. 2011;43(2):117-120. doi: 10.1038/ng.735
9. Zhou Y, Guo Y, Ye W, et al. RS11212617 is associated with metformin treatment response in type 2 diabetes in Shanghai local chinese population. Int J Clin Pract. 2014;68(12):1462-1466. doi: 10.1111/ijcp.12534
10. vanleeuwen N, Nijpels G, Becker ML, et al. A gene variant near ATM is significantly associated with metformin treatment response in type 2 diabetes: a replication and metaanalysis of five cohorts. Diabetologia. 2012;55(7):1971-1977. doi: 10.1007/s00125-012-2537-x
11. Florez JC, Jablonski KA, Taylor A, et al. The c allele of ATM rs11212617 does not associate with metformin response in the diabetes prevention program. Diabetes Care. 2012;35(9):1864-1867. doi: 10.2337/dc11-2301
12. Vilvanathan S, Gurusamy U, Mukta V, et al. Allele and genotype frequency of a genetic variant in ataxia telangiectasia mutated gene affecting glycemic response to metformin in South Indian population. Indian J Endocrinol Metab. 2014;18(6):850-854. doi: 10.4103/2230-8210.119944.
About the Authors
Irina A. BondarNovosibirsk State Medical University
Russian Federation
MD, PhD, Professor
Competing Interests:
Автор декларирует отсутствие конфликта (двойственности) интересов при написании данной статьи
Olesia Y. Shabelnikova
Novosibirsk State Medical University; Novosibirsk State Regional Clinical Hospital
Russian Federation
MD, PhD
Competing Interests:
Автор декларирует отсутствие конфликта (двойственности) интересов при написании данной статьи
Ekaterina A. Sokolova
Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences; Novosibirsk State University
Russian Federation
PhD
Competing Interests:
Автор декларирует отсутствие конфликта (двойственности) интересов при написании данной статьи
Maksim L. Filipenko
Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Sciences; Novosibirsk State University
Russian Federation
PhD
Competing Interests:
Автор декларирует отсутствие конфликта (двойственности) интересов при написании данной статьи
Supplementary files
![]() |
1. Рисунок 2 | |
Subject | ||
Type | Исследовательские инструменты | |
Download
(211KB)
|
Indexing metadata ▾ |
|
2. Figure 1. Frequency of minor C allele at the rs11212617 polymorphic locus in the ATM gene in the groups with different responses to metformin therapy. | |
Subject | ||
Type | Исследовательские инструменты | |
View
(21KB)
|
Indexing metadata ▾ |
|
3. Figure 2. Graphical representation of the results of meta-analysis and bias analysis in publications. | |
Subject | ||
Type | Исследовательские инструменты | |
View
(193KB)
|
Indexing metadata ▾ |
|
4. Table 2. Logistic regression analysis of the relationship between the attainment of target HbA1с level of <7% during metformin therapy and genotypes of polymorphic locus rs11212617 in the ATM gene | |
Subject | ||
Type | Исследовательские инструменты | |
View
(54KB)
|
Indexing metadata ▾ |
|
5. Table 3. Metaanalysis of the studies of the relationship between the attainment of target HbA1с level of <7% during MF therapy and the genotype of polymorphic locus rs11212617 in the ATM gene | |
Subject | ||
Type | Исследовательские инструменты | |
View
(215KB)
|
Indexing metadata ▾ |
|
6. Table 1. Clinical and laboratory characteristics of the examined patients | |
Subject | ||
Type | Исследовательские инструменты | |
View
(51KB)
|
Indexing metadata ▾ |
Review
For citations:
Bondar I.A., Shabelnikova O.Y., Sokolova E.A., Filipenko M.L. Research of association of the polymorphic locus rs11212617 ATM gene with the response to therapy with metformin in patients with type 2 diabetes. Problems of Endocrinology. 2017;63(1):9-16. https://doi.org/10.14341/probl20176319-16

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0).