Table of Contents Author Guidelines Submit a Manuscript
Journal of Diabetes Research
Volume 2015, Article ID 485132, 6 pages
http://dx.doi.org/10.1155/2015/485132
Research Article

The Genetic Profile from HLA and Non-HLA Loci Allows Identification of Atypical Type 2 Diabetes Patients

1Biodiversity and Genetic Department, Instituto de Investigaciones Biológicas Clemente Estable, 11600 Montevideo, Uruguay
2Teaching Care Unit, Unit of Diabetes “Hospital Pasteur”, ASSE, Ministry of Public Health, 11600 Montevideo, Uruguay
3Diabetology Service of “Centro de Asistencia del Sindicato Médico del Uruguay”, 11800 Montevideo, Uruguay

Received 29 January 2015; Accepted 23 March 2015

Academic Editor: Joseph Fomusi Ndisang

Copyright © 2015 Matias Fabregat et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The complex diagnosis and treatment of diabetes highlight the need for markers to define how to monitor patients correctly during the course of their disease. Different studies demonstrate the existence of patients who cannot be clearly classified. We have previously shown that it is possible to differentiate “atypical diabetic patients” based on genotyping the HLA. In this work we show that the analysis of non-HLA related to type 1 diabetes in the INS-VNTR, SNP rs689, and rs3842753 improves the identification of these patients. We genotyped 913 individuals comprising controls from the general population and “classic” and “atypical” diabetic patients. We compared the distribution of these loci and analyzed linkage disequilibrium. The haplotype was in LD for all the SNPs that were evaluated. Regarding their association with the disease, the haplotype IAC was associated with type 1 (odds 2.60, 1.82–3.72, CI 95%) and “atypical diabetes” (odds 1.50, 1.01–2.23, CI 95%), whereas we did not observe an association with type 2 diabetes. Therefore, our results confirm that atypical diabetes is a different entity of the disease where the patient presents with a genetic background of T1D and a T2D phenotype, findings that are likely to be relevant for patient diagnosis and management in the clinic.

1. Introduction

Uruguay is a little country located in South America with a population of 3.286.314 inhabitants according to the last census [1]. The population has been characterized as of European descent, with a small contribution of African descendants and native amerindian popultion [2]. Diabetes mellitus is a worldwide public health problem and Uruguay is not an exception showing a high prevalence of the disease, 8%, and an equal percent that are likely undiagnosed despite the small size of its population [3]. Diabetes is a multifactorial disease where genetic and environmental factors interact. This complexity could be an important influence in the actual classification and treatment of this pathology. Previously, in an effort to link the current classification with advances at molecular genetics we found that in several studies a significant proportion of diabetic patients did not present a correlation between genotype and phenotype [46]. Type 2 diabetes (T2D) diagnosed patients using international criteria [7] presented Human Leukocyte Antigen (HLA) alleles associated with type 1 diabetes (T1D) and positive antibodies in an intermediate value between both types, thus suggesting the coexistence of both types of diabetes. Several authors have published about this kind of diabetes [814] and these “atypical patients” cannot be grouped into any of the established groups by international guidelines. More importantly, these atypical patients do not have an adequate response when undergoing treatment [15].

Recently Maruthur and colleagues found evidence that in T2D the pharmacogenetic interactions for some antihyperglycemic drugs are consistent with their pharmacokinetics and pharmacodynamics, showing the importance of genetics variants [16]. In line with these observations we hypothesized that in “atypical patients” non-HLA genes, in addition to the HLA susceptibility alleles, could be influencing the lack of genotype-phenotype correlation and the course and development of the disease.

Here we perform an in depth analysis of the presence of other non-HLA loci related to T1D in this atypical diabetic population, particularly focusing on the insulin gene (Insulin-Varaibale number tandem repeat, SNP rs689, and rs3842753). We show that in the Uruguayan population these variants are in complete linkage disequilibrium [17]. This region is located on chromosome 11p15.5, in the promoter region of the insulin gene, and affects the level of transcription of INS and IGF2 genes [10, 18]. Mathematical models estimate that, adding the HLA and VNTR contribution, up to 60% of the total genetic component of T1D is obtained [19]. VNTR polymorphisms located in the insulin promoter, 596 bp upstream of the insulin start codon, consist of a highly polymorphic tandem repeat sequence ACAGGGGTGTGGGG (14 pb) whose variants are grouped into three classes of alleles according to their length [20, 21]. Regarding our population, VNTR is in linkage disequilibrium with rs689 where class I VNTR allele is associated with autoimmunity and inherited linked to the Adenine SNP allele. Furthermore the class III VNTR alleles are associated with autoimmunity protection and inherited linked to the Timina SNP allele. Population studies have shown that class III VNTR alleles have a dominant protective effect and are associated with a 60–70% reduction in the risk of developing T1D [1923]. Even in the presence of autoantibodies and a high-risk HLA genotype, individuals carrying class III VNTR alleles present a significant reduction in the probability of developing the disease [24, 25]. Based on this association, major studies have used 23HphI (rs689) as an associated marker to study the VNTR [21, 26, 27].

Another SNP in linkage disequilibrium with the variants described above is rs3842753, a transversion at position +1140 relative to the initial translation codon of the INS gene. This polymorphism is located in the 3′UTR region of the INS mRNA and the cytosine-containing variant may cause instability of the mRNA [28]. In European and Uruguayan populations, the cytosine-containing allele is associated with T1D and is in complete linkage disequilibrium with the VNTR allele of class I and rs698 [17, 25]. Given the likely coexistence of both types of diabetes in our population and the presence of linkage disequilibrium as described in previous papers, we propose the in depth analysis of the genetic profile of these “atypical patients.” The analysis should include patients that present a phenotype characteristic of T2D with HLA susceptibility alleles associated with T1D compared to (a) the general population and (b) type 1 diabetes and type 2 diabetes patients without a HLA association.

2. Materials and Methods

In this unmatched case-control study, a total of 913 individuals, including diabetes patients (413) and controls (500), were enrolled between 2004 and 2012. Recruitment of patients was done by the outpatient health center of Montevideo city.

2.1. Control Samples

We selected 500 unrelated individuals from the DNA bank of the Department of Biodiversity and Genetics, IIBCE. This collection is a representative sample of Montevideo’s general population and was randomly selected from 15 different medical institutions, public and private, when individuals assisted to their annual routine control.

Due to local legislation at the time (ref. 1081/1996), the MSP (Public Health Ministry) did not authorize collecting information about the clinical characteristics of patients when the study involved analyses of DNA samples. Therefore, the only characteristic recorded for this population is the age (>18 years old).

2.2. Diabetic Patient Samples

413 patients from 3rd attention level at Clinics for Diabetes of Reference Health Centers of Montevideo were analyzed. For the preparation of this study we only considered those patients receiving comprehensive care of their diabetes, following a nutritional plan and presenting a good adherence to physical activity according to their functional ability within the recommendations of the American Diabetes Association (ADA) and “Asociación Latinoamericana de Diabetes” (ALAD) [7, 29]. The diabetes samples were subclassified as type 1 diabetes, type 2 diabetes, and “atypical diabetes.”

2.2.1. Type 1 Diabetes Patients (168)

Those patients were defined according to the ADA criteria [22], with an age of onset < 15 years and a body mass index < 25 kg/mts2.

2.2.2. Type 2 Diabetes Patients, Classical and Atypical (153 + 92)

The population diagnosed with T2D was divided into two groups based on the presence or absence of T1D HLA susceptibility alleles described in the Uruguayan population [15, 30] according to the following inclusion criteria.

“Atypical diabetes” (92 individuals) (a) patients who had good adherence to the treatment; (b) they fulfilled the objectives of education and nutrition plans according to international guidelines; (c) present doubts on classification of diabetic type and/or not good therapeutic response (two consecutive measurements of glycated hemoglobin within three months not reduced in 1.5% [31]); (d) patients with susceptibility HLA alleles for autoimmune disease. We considered DQB1 0201-0302 and DR 3-4 as susceptible ones in the Uruguay, (e) body mass index ≥ 25 kg/mts2 [32].

2.2.3. Type 2 Diabetes (153 Individuals)

Those patients were fulfilling the same requirements a and b of atypical patients but without diagnostic doubts, responded to treatment, and do not present HLA alleles associated with autoimmune disease.

Samples from patients who had other endocrine disorders or tumors were excluded.

All subjects were interviewed by specialist medical doctors and gave a written informed consent to participate in the study. The protocol was approved by the Ethical Committee of the Ministry of Public Health (MSP) and the corresponding ethical committee of each participating institution.

2.3. Molecular Analysis

High molecular weight DNA extraction was performed from peripheral blood by standard phenol chloroform protocol.

The HLA typing was performed by reverse ASO technique (Innogenetics Ltd., Belgium, UE). INS 5VNTR, rs3842753, and HLA were processed in a previous work [4, 5, 15, 17, 31].

The SNP rs689 (−23 HPh1) for controls and T2D were processed in a previous work [16, 33]: Atypical diabetes patients were genotyped by PCR-RFLP with the following primers: forward-5′AGCAGGTCTGTTCCAAGG-3′ and reverse-5′CTTGGGTGTGTAGAAGAAGC-3′ which amplifies a 360 bp fragment. The identification of the genotypes was performed by digesting the DNA fragments with BsmAI. Sequencing was used for confirmation of genotypes (Macrogen, Ltda, Korea).

2.4. Statistical Analysis

(a) Calculations of power sample size were done with the Epi Info 3.4.3 database and the statistics software for public health professionals and the Quanto statistical package (http://biostats.usc.edu/software) considering the prevalence of diabetes in Uruguay (8%) [3]. We assumed 95% of confidential interval and power 80% in an unmatched case-control design.

(b) The statistical analysis of a polymorphism was done with the web tool for SNP analysis SNPStats (http://bioinfo.iconcologia.net/snpstats/start.htm).

(c) Selection of inheritance model: best inheritance mode was selected according to SNPStat program. The statistic follows a chi-square distribution with degrees of freedom equal to the number of additional parameters in the more complex model and Akaike information criterion (AIC) and Bayesian information criterion (BIC). The statistical test values were calculated via an exact test. The reference category used by program was the homozygous form.

(d) Hardy-Weinberg equilibrium (HWE) for allelic and genotype frequencies was tested by chi-square test.

(e) Association between one polymorphism and disease: we make the contingency table and then apply a chi-square test and the estimation of the OR (odds ratio) for each genotype with respect to the reference genotype (Epi Info and SNPStat).

(f) Haplotype analysis: D statistic (under the assumption of no association) and correlation coefficient between alleles and the observed frequency were done. Linkage disequilibrium and haplotype were calculated with SNPStat program. Analysis of multiple SNPs and haplotype and analysis of association between haplotypes and disease were done.

3. Results

Nine hundred and thirteen DNAs were analyzed in this study. Five hundred correspond to samples obtained from the general population of Uruguay and four hundred and thirteen correspond to diabetes patients. First we analyzed the clinical characteristics of diabetes patients and performed the comparison between the three groups defined previously (Table 1). We found a significant statistical difference when we compared T1D with any other subgroup in all variables considered. However, the only difference between T2D and “atypical diabetes” was in the level of HDL. These results are in accordance with previous reports [14]. Our genetic analysis showed that all SNPs analyzed in all groups were in Hardy-Weinberg equilibrium (HWE). The best heredity model for each SNP analyzed was the log-additive model, according to the AIC and BIC criteria of the SNPstat program. Allelic frequencies for all samples and association with the disease are shown in Tables 2(a), 2(b), and 2(c).

Table 1: Clinical characteristics of diabetes patients.
Table 2: Allelic frequencies of the SNPs: (a) type 1 diabetes, (b) type 2 diabetes, and (c) atypical diabetes.

According to the D′ value and the correlation coefficient between alleles and the observed frequency, the evaluated haplotype was in linkage disequilibrium with all SNPs analyzed. The most frequent haplotype was IAC in controls and diabetes patients (Table 3). Haplotype with frequencies minor 1% was considered rare. Regarding the association with the disease, the high protective effect (IIITA) was present in T1D. The contingency table for haplotypes associated with T1D versus no presence of this haplotype revealed an association in T1D and “atypical diabetes” but resulted in no significant association in T2D (Table 4).

Table 3: Haplotype frequencies for all groups.
Table 4: Haplotype association with the disease.

4. Discussion

Our study allowed us to analyze the genetic profile of “atypical diabetes” comparing the phenotype of type 2 diabetes with HLA susceptibility alleles in the general population, T1D, and T2D without HLA associated. Based on our results we propose that this marker can aid in the diagnosis and management of patients who present difficulties in control and follow-up. In addition, our analysis of the INS VNTR locus, as the second marker of importance for T1D susceptibility, further supports the notion that patients with difficulties in all approach areas (diagnosis, treatment, and evolution) have a strong genetic basis.

The last ADA guidelines classified diabetes in four clinical categories but recognized that some patients cannot be clearly classified as type 1 or 2 diabetes. Importantly, such difficulties in classification may occur at any age and genetic studies could improve the timely diagnosis in these kinds of patients, without having to wait for the disease to progress.

Mathematical models estimate that together the HLA and INS VNTR contribution may account for near 60% of the total genetic basis for developing T1D [19]. Therefore these markers, typically associated with autoimmune diabetes, when present in T2D patients, can modify the expected development of the pathology.

It has been shown that the second susceptibility locus for autoimmune disease, after HLA, corresponds to a mini satellite in the insulin gene (INS-VNTR) and in Caucasians populations this VNTR is in linkage disequilibrium with −23 HphI SNP. Previous studies estimate that 10% of the genetic susceptibility to T1D corresponds to this mini satellite [20].

In previous works, we showed that there are patients in our population who cannot be classified as either type 1 or type 2 diabetes according to the international guidelines. The clinical presentation, evolution, and difficulty in reaching expected therapeutic goals make them “atypical patients” in relation with HLA alleles [14]. Now we have gained insight into understanding the genetic basis of “atypical diabetes” by showing an association with genes different from HLA. It is noteworthy that, in these patients, there are no homozygous individuals for allele III of 5′VNTR. The protective allele III in the promoter region of INS is associated with a greater functional reserve of β cells in response to a direct hyperglycemic stimulus [34]. Although we did not analyze the functional reserve of beta cells, it is striking to find a high frequency of susceptibility allele I in “atypical patients” and lower presence of III alleles in their genotypes. This observation could suggest an earlier beta cell failure in these patients and therefore could be one possible explanation to understand the cellular basis of their bad treatment response, a possibility that should be analyzed in future investigations.

In conclusion, our results confirm that atypical diabetes is a different entity of the disease where patients combine a genetic background of T1D with a T2D phenotype, thus highlighting the importance of genetics as a new tool in clinical practice.

Conflict of Interests

There are no financial competing interests.

Authors’ Contribution

Matias Fabregat conducted experimental procedures and validations, data acquisition, compilation, analysis, and interpretation and also was involved in retrieving conceptual information and draft paper preparation; Mariana Fernandez contributed to draft paper preparation and performed critical revision to structure the content intellectually. Gerardo Javiel provided substantial contributions to the concept and design of the work. He was in charge of the admission protocol and monitoring of patients. He was involved in retrieving medical conceptual information and draft paper preparation, performed a critical revision to structure the content intellectually, and gave approval for the final version to be published. Graciela Vitarella was in charge of the admission protocol and monitoring of patients. She was involved in retrieving medical conceptual information and draft paper preparation, performed a critical revision to structure the content intellectually, and gave approval for the final version to be published. Adriana Mimbacas provided substantial contributions to the concept and design of the work and controlled experimental procedures, validations of the results, and statistical analysis. She was involved in retrieving conceptual information and draft paper preparation, performed critical revision to structure the content intellectually, and gave approval for the final version to be published.

Acknowledgments

The authors thank all the patients that participated in the study and the health personnel who contributed to the development of this study. They thank Dr. Jose Luis Badano for correcting the English paper.

References

  1. INE (Instituto Nacional de Estadística), 2011, http://www.ine.gub.uy/censos2011/index.html.
  2. M. Sans, G. Figueiro, and P. C. Hidalgo, “A new mitochondrial C1 lineage from the prehistory of Uruguay: population genocide, ethnocide, and continuity,” Human Biology, vol. 84, no. 3, pp. 287–305, 2012. View at Publisher · View at Google Scholar · View at Scopus
  3. R. Ferrero and T. Garcia, “Survey on the prevalence of diabetes in Uruguay,” Archives of Internal Medicine, vol. 27, pp. 1–12, 2005. View at Google Scholar
  4. A. Mimbacas, L. García, P. Zorrilla et al., “Genotype and phenotype correlations in diabetic patients in Uruguay,” Genetics and Molecular Research, vol. 8, no. 4, pp. 1352–1358, 2009. View at Publisher · View at Google Scholar · View at Scopus
  5. A. Mimbacas and G. Javiel, “Diabetes type 1 and type 2: what is behind a classification?” in Type 1 Diabetes. Complications, Pathogenesis, and Alternative Treatments, C.-P. Liu, Ed., pp. 287–304, InTech, 2011. View at Google Scholar
  6. A. Mimbacas, G. Vitarella, J. Souto et al., “The phenotype masks the genotype: a possible new expression of diabetes,” Journal of Pediatric Genetics, vol. 1, no. 2, pp. 131–134, 2012. View at Publisher · View at Google Scholar
  7. American Diabetes Association. Position Statement, “Standards of medical care in diabetes,” Diabetes Care, vol. 35, supplement 1, pp. S11–S63, 2012. View at Publisher · View at Google Scholar
  8. I. M. Libman and D. J. Becker, “Coexistence of type 1 and type 2 diabetes mellitus: ‘double’ diabetes?” Pediatric Diabetes, vol. 4, no. 2, pp. 110–113, 2003. View at Publisher · View at Google Scholar · View at Scopus
  9. H. Kolb and T. Mandrup-Poulsen, “An immune origin of type 2 diabetes?” Diabetologia, vol. 48, no. 6, pp. 1038–1050, 2005. View at Publisher · View at Google Scholar · View at Scopus
  10. T. Tuomi, “Type 1 and type 2 diabetes: what do they have in common?” Diabetes, vol. 54, supplement 2, pp. S40–S45, 2005. View at Publisher · View at Google Scholar · View at Scopus
  11. P. Pozzilli and R. Buzzetti, “A new expression of diabetes: double diabetes,” Trends in Endocrinology & Metabolism, vol. 18, no. 2, pp. 52–57, 2007. View at Publisher · View at Google Scholar · View at Scopus
  12. P. Pozzilli, C. Guglielmi, E. Pronina, and E. Petraikina, “Double or hybrid diabetes associated with an increase in type 1 and type 2 diabetes in children and youths,” Pediatric Diabetes, vol. 8, no. 9, pp. 88–95, 2007. View at Publisher · View at Google Scholar · View at Scopus
  13. M. A. Stone, J. Camosso-Stefinovic, J. Wilkinson, S. de Lusignan, A. T. Hattersley, and K. Khunti, “Incorrect and incomplete coding and classification of diabetes: a systematic review,” Diabetic Medicine, vol. 27, no. 5, pp. 491–497, 2010. View at Publisher · View at Google Scholar · View at Scopus
  14. S. E. Inzucchi, R. M. Bergenstal, J. B. Buse et al., “Management of hyperglycemia in type2 diabetes: a patient-centered approach,” Diabetes Care, vol. 35, pp. 1364–1379, 2012. View at Google Scholar
  15. M. Fernández, M. Fabregat, G. Javiel, and A. Mimbacas, “HLA alleles may serve as a tool to discriminate atypical type 2 diabetic patients,” World Journal of Diabetes, vol. 5, no. 5, pp. 711–716, 2014. View at Publisher · View at Google Scholar
  16. N. M. Maruthur, M. O. Gribble, W. L. Bennett et al., “The pharmacogenetics of type 2 Diabetes: a systematic review,” Diabetes Care, vol. 37, no. 3, pp. 876–886, 2014. View at Publisher · View at Google Scholar · View at Scopus
  17. C. Gascue, Análisis de polimorfismos de riesgo a la diabetes tipo 1 en una población trihíbrida, Uruguay Montevideo, Uruguay [M.S. thesis], LAP LAMBERT Academic Publishing GmbH & Company KG, 2007, http://www.bod.de/shop/buch.html?&bookId=587987.
  18. M. I. McCarthy and P. Froguel, “Genetic approaches to the molecular understanding of type 2 diabetes,” American Journal of Physiology - Endocrinology and Metabolism, vol. 283, no. 2, pp. E217–E225, 2002. View at Publisher · View at Google Scholar · View at Scopus
  19. S. T. Bennett and J. A. Todd, “Human type 1 diabetes and the insulin gene: principles of mapping polygenes,” Annual Review of Genetics, vol. 30, pp. 343–370, 1996. View at Publisher · View at Google Scholar · View at Scopus
  20. G. I. Bell, M. J. Selby, and W. J. Rutter, “The highly polymorphic region near the human insulin gene is composed of simple tandemly repeating sequences,” Nature, vol. 295, no. 5844, pp. 31–35, 1982. View at Publisher · View at Google Scholar · View at Scopus
  21. D. Cimponeriu, P. Apostol, I. Radu et al., “The insulin polymorphism -23Hph increases the risk for type 1 diabetes mellitus in the Romanian population,” Genetics and Molecular Biology, vol. 33, no. 4, pp. 610–614, 2010. View at Publisher · View at Google Scholar
  22. S. T. Bennett, A. M. Lucassen, S. C. L. Gough et al., “Susceptibility to human type 1 diabetes at IDDM2 is determined by tandem repeat variation at the insulin gene minisatellite locus,” Nature Genetics, vol. 9, no. 3, pp. 284–292, 1995. View at Publisher · View at Google Scholar · View at Scopus
  23. S. T. Bennett, A. J. Wilson, F. Cucca et al., “IDDM2-VNTR-encoded susceptibility to type 1 diabetes: dominant protection and parental transmission of alleles of the insulin gene-linked minisatellite locus,” Journal of Autoimmunity, vol. 9, no. 3, pp. 415–421, 1996. View at Publisher · View at Google Scholar · View at Scopus
  24. A. Pugliese, Z. L. Awdeh, C. A. Alper, R. A. Jackson, and G. S. Eisenbarth, “The INS B allele (1,428 Fok I) is associated with reduced IDDM incidence in autoantibody-positive first degree relatives despite the presence of autoantibodies and high-risk HLA alleles,” Transplantation Proceedings, vol. 27, no. 6, p. 3392, 1995. View at Google Scholar · View at Scopus
  25. A. M. Lucassen, G. R. Screaton, C. Julier, T. J. Elliott, M. Lathrop, and J. I. Bell, “Regulation of insulin gene expression by the IDDM associated, insulin locus haplotype,” Human Molecular Genetics, vol. 4, no. 4, pp. 501–506, 1995. View at Publisher · View at Google Scholar · View at Scopus
  26. J. D. H. Stead and A. J. Jeffreys, “Allele diversity and germline mutation at the insulin minisatellite,” Human Molecular Genetics, vol. 9, no. 5, pp. 713–723, 2000. View at Publisher · View at Google Scholar · View at Scopus
  27. P. Cejkova, P. Novota, M. Cerna et al., “HLA DRB1, DQB1 and insulin promoter VNTR polymorphisms: interactions and the association with adult-onset diabetes mellitus in Czech patients,” International Journal of Immunogenetics, vol. 35, no. 2, pp. 133–140, 2008. View at Publisher · View at Google Scholar · View at Scopus
  28. C. Julier, A. Lucassen, P. Villedieu et al., “Multiple DNA variant association analysis: application to the insulin gene region in type I diabetes,” The American Journal of Human Genetics, vol. 55, no. 6, pp. 1247–1254, 1994. View at Google Scholar · View at Scopus
  29. Asociación Latinoamericana de Diabetes (ALAD), 2012 Guias ALAD de diagnóstico, control y tratamiento de la diabetes tipo 2, http://alad-latinoamerica.org/.
  30. A. Mimbacas, F. Perez-Bravo, P. C. Hidalgo et al., “Association between diabetes type 1 and DQB1 alleles in a case-control study conducted in Montevideo, Uruguay,” Genetics and Molecular Research, vol. 2, pp. 29–35, 2003. View at Google Scholar
  31. D. M. Nathan, J. B. Buse, M. B. Davidson et al., “Medical management of hyperglycemia in type 2 diabetes: a consensus algorithm for the initiation and adjustment of therapy,” Diabetes Care, vol. 32, no. 1, pp. 193–203, 2009. View at Publisher · View at Google Scholar · View at Scopus
  32. World Health Organization, “BMI classification,” http://www.who.int/gho/ncd/risk_factors/bmi_text/en/.
  33. P. Zorrilla, Polimorfismos asociados a Diabetes Mellitus Tipo 1 en el gen PTPN22 en la población de Montevideo, Editorial Academica Espanola, 2011, https://www.eae-publishing.com/catalog/search.
  34. W. Fendler, K. Wyka, A. Cieślik-Heinrich, E. Polakowska, J. Jastrzebska, and W. Młynarski, “The 5′VNTR proinsulin gene (INS) polymorphism and the functional reserve of β cells in the intravenous glucose tolerance test,” Pediatric Endocrinology, Diabetes and Metabolism, vol. 17, no. 1, pp. 5–9, 2011. View at Google Scholar · View at Scopus