Skip to main content

Novel susceptibility loci identified in a genome-wide association study of type 2 diabetes complications in population of Latvia

Abstract

Background

Type 2 diabetes complications cause a serious emotional and economical burden to patients and healthcare systems globally. Management of both acute and chronic complications of diabetes, which dramatically impair the quality of patients' life, is still an unsolved issue in diabetes care, suggesting a need for early identification of individuals with high risk for developing diabetes complications.

Methods

We performed a genome-wide association study in 601 type 2 diabetes patients after stratifying them according to the presence or absence of four types of diabetes complications: diabetic neuropathy, diabetic nephropathy, macrovascular complications, and ophthalmic complications.

Results

The analysis revealed ten novel associations showing genome-wide significance, including rs1132787 (GYPA, OR = 2.71; 95% CI = 2.02–3.64) and diabetic neuropathy, rs2477088 (PDE4DIP, OR = 2.50; 95% CI = 1.87–3.34), rs4852954 (NAT8, OR = 2.27; 95% CI = 2.71–3.01), rs6032 (F5, OR = 2.12; 95% CI = 1.63–2.77), rs6935464 (RPS6KA2, OR = 2.25; 95% CI = 6.69–3.01) and macrovascular complications, rs3095447 (CCDC146, OR = 2.18; 95% CI = 1.66–2.87) and ophthalmic complications. By applying the targeted approach of previously reported susceptibility loci we managed to replicate three associations: MAPK14 (rs3761980, rs80028505) and diabetic neuropathy, APOL1 (rs136161) and diabetic nephropathy.

Conclusions

Together these results provide further evidence for the implication of genetic factors in the development of type 2 diabetes complications and highlight several potential key loci, able to modify the risk of developing these conditions. Moreover, the candidate variant approach proves a strong and consistent effect for multiple variants across different populations.

Peer Review reports

Background

The past few decades have shown a marked increase in the number of patients with diabetes rising from 151 million (4.6% of the global population) in 2000 to 463 million (9.3%) in 2019 [1]. The risk of type 2 diabetes (T2DM), the most common type of diabetes, is modified by a strong interaction between environmental and genetic factors [2, 3]. T2DM is a multifactorial disease with a population-specific heritability (26% in the European population) [4]. A number of common variants implicated in the pathogenesis and genetic architecture of T2DM have been identified so far, some of them also capable of modifying the pharmacologic response to antidiabetic drugs [5, 6].

Uncontrolled T2DM can lead to long-term illnesses or chronic health conditions divided into microvascular complications, such as diabetic retinopathy, nephropathy, neuropathy, and macrovascular complications, including stroke, heart disease, and peripheral vascular disease [7], accounting for 53% of direct health costs of diabetes with cardiovascular and renal complications contributing to the greatest financial burden [8, 9]. Moreover, diabetic patients show increased all-cause mortality rates, especially cardiovascular deaths (HR 2.6, 95% CI 1.4–4.7) [10].

Studies show evidence of considerable genetic component predisposing to diabetic complications, explaining even around 50% of the risk of proliferative retinopathy [11]. In the last few decades, genetic research including genome-wide association studies (GWAS), linkage analysis, and candidate gene approach has revealed several susceptibility loci for diabetic retinopathy and nephropathy (VEGF, CAT, FTO, UCP1, and INSR), and also macrovascular complications (ADIPOQ). Nevertheless, they explain only a small proportion of the phenotypic variation observed in T2DM patients [12,13,14,15,16,17], justifying a need for identification of novel genetic risk factors for T2DM complications and improvement of knowledge about molecular mechanisms underlying these comorbid conditions.

Since the high impact of population specificity for the discrimination of genetic variants and their contribution to the phenotype of interest is evidenced by a number of SNPs that failed to replicate in different populations, both discovery and replication studies in populations of different ancestries are needed [18, 19]. An example of population-specific allele frequency is rs61736969, the risk variant of T2DM, which is located in the TBC1D4 (TBC1 Domain Family Member 4) gene. It is highly frequent (minor allele frequency of 17%) in the Greenlandic population, nevertheless, it has not been present so far in European individuals, most probably due to different linkage disequilibrium patterns [20]. Although poor glycemic control is considered to be the driving factor of T2DM complications, early genotype-based identification of individuals at high risk of diabetic complications may promote the prevention or, at least, delay of the disease [21]. In this study, we aimed to discriminate novel susceptibility loci for T2DM complications and replicate the findings of other GWAS in the study cohort from the Genome Database of the Latvian Population (LGDB) [22].

Methods

Study group and phenotype definitions

In total, the study cohort consisted of 601 T2DM patients of European ancestry with and without a medical history of diabetes complications, selected from the participants of LGDB (recruited from June 2007 to November 2016) according to the following inclusion criteria: (1) clinically confirmed diagnosis of T2DM (E11 diagnosis code, ICD-10), (2) information on age at diagnosis, sex, weight, and height available, (3) records of national diabetes registry and Latvian hospital inpatient discharges available. Written broad consent was obtained from every subject during the recruitment in LGDB.

The collection of blood samples and relevant anthropometric data was ensured by LGDB according to their standard procedures [22]. Associated clinical data, including the diagnosis date of T2DM, date and type of diabetes complications, HbA1c measures, and medications used, were obtained from the records of Diabetes registry, Latvian hospital inpatient discharges, outpatient progress notes, and pharmacy recipe records, provided by The Centre for Disease Prevention and Control of Latvia and National Health Service of Latvia (Approval No. 3, Decision No. 7.1–3/3). The data about diabetic complications present for T2DM patients involved in LGDB were applied for accurate stratification of 601 T2DM patients in four phenotype-based groups according to the type of complications experienced: diabetic neuropathy, diabetic nephropathy, ophthalmic complications, and macrovascular complications. The definition of phenotypes and patient stratification in different complication groups was done as follows:

  1. 1.

    Diabetic neuropathies: clinical diagnosis codes (ICD-10) E11.4 and E11.5, records of amputation of the leg/toe, gangrene, shunting and angioplasty, and presence of intermittent claudication or fresh ulcers since the diagnosis of T2DM.

  2. 2.

    Diabetic nephropathies: clinical diagnosis code E11.2 or records of kidney failure, kidney transplantation, renal replacement therapy, microalbuminuria, hemodialysis, peritoneal dialysis performed after the diagnosis of T2DM.

  3. 3.

    Ophthalmic complications: clinical diagnosis code E11.3 or records of photocoagulation, maculopathy, retinopathy, operative therapy, blindness made since the diagnosis of T2DM.

  4. 4.

    Macrovascular complications: clinical diagnosis codes I95, I20, I21, I24, I25, I50, I60, I61, I63, I64, and records of coronary shunting and angioplasty after the diagnosis of T2DM.

ICD-10-based phenotype definitions corresponding to the Latvian guidelines of diabetes management are generally used in clinical practice in Latvia. Subjects with the above-mentioned diagnosis codes or medical events recorded were considered as cases in their corresponding complication groups, while T2DM patients with no evidence of complications of interest during their follow-up period were recognized as controls in the particular group. Subjects experiencing specific diabetes complications before the set of T2DM diagnosis were excluded from the analysis of a particular complication group, explaining the variable total number of individuals among all complication groups tested. The follow-up period which coincides with diabetes duration was considered as time since the set of T2DM diagnosis until the date of diabetes complication recorded for cases or the date of the last entry in the National registry for control subjects. Administration of medications was considered in a group-specific manner, accounting for angiotensin II receptor blockers and angiotensin-converting enzyme inhibitors in the analysis of all complication groups analyzed and additional lipid-modifying agents in the analysis of macrovascular complications. In order to adjust for the inter-individual variability of glycaemic control, a key factor in the development of T2DM complications, the median HbA1c level during the observation period was fitted as a covariate, irrespective of the antidiabetic therapy used. Sex, age at the diagnosis of T2DM, body mass index (BMI), diabetes duration, and use of particular medications were also included among the covariates.

DNA extraction and genotyping

Within the framework of this study DNA samples from 601 T2DM patients were used. DNA was isolated from peripheral blood leukocytes using a phenol–chloroform extraction method according to LGDB standard procedures [22]. DNA samples were genotyped with the Infinium Global Screening Array (Illumina, USA) on the iScan System microarray scanner (Illumina, USA). Illumina Genome Studio v2.0 was used to convert raw data into PLINK format and workflow described in Marees et al. [23]. used for data quality control. SHAPEIT v2.r900 [24] and IMPUTE2 [25] were used for genotype phasing and genotype imputing. Imputed data were filtered using the following parameters: marker correlation (INFO) > 0.8, hard call threshold 0.1, minor allele frequency > 1%, Missingness < 2%.

Statistical analysis

Association analyses corresponding to four different complication groups (macrovascular complications, diabetic neuropathy, diabetic nephropathy, ophthalmic complications) were performed using PLINK v1.9 logistic regression with covariates: median HbA1c, sex, age at the diagnosis, diabetes duration, BMI, medications used. A genome-wide significance threshold of P < 5 × 10−8 was defined.

For the targeted analysis candidate variants were selected from GWAS Central [26] (http://www.gwascentral.org/) and GWAS Catalog [27] (https://www.ebi.ac.uk/gwas/home), based on previously reported association with T2DM complications: diabetic neuropathy (4 allelic variants selected), macrovascular complications (43 allelic variants selected), ophthalmic complications (98 allelic variants selected) and diabetic nephropathy (49 allelic variants selected). A complete list of selected SNPs, their positions, and associated traits is provided in Additional file 1. False discovery rate (FDR) according to the Benjamini–Hochberg procedure was calculated to account for multiple testing and the threshold was set at 0.05.

Manhattan plots and Q-Q plots were generated in R v3.5.3 using the qqman package, while the Venn diagram was developed in the online visualization tool Venny 2.1.0. Statistical analysis of anthropometric measures and biochemical data was performed in R v3.5.3. by applying the Wilcoxon rank-sum test and Pearson's chi-squared test with a p-value threshold < 0.05. For identification of the functional role of allelic variants expression quantitative trait locus (eQTL) analysis was conducted by using the open-access Genotype Tissue Expression (GTEx) database [28]. The tissue types for eQTL analysis were carefully selected considering the etiology of the diseases (artery, nerve, heart, skin, blood) [29]. The p-value threshold of 0.05 was used to discriminate significant associations. Variant Effect Predictor and Linkage Disequilibrium (LD) data from 1000 Genome project (Utah Residents (CEPH) with Northern and Western European Ancestry) were employed to explore the functional consequences of each variant and other variants in LD [30]. To evaluate the potential SNP effects on quantitative phenotypes, analyses of variant association with HbA1C and BMI were performed using PLINK v1.9—assoc function for quantitative phenotypes, where the genome-wide significance threshold of P < 5 × 10−8 was used to identify the significant hits.

Results

Genome-wide association analysis

We studied a cohort of 601 T2DM patients of which 241 were men and 360 women aged 22 to 82 years (average age 56.86 ± 10.24 years) stratified as controls or cases in four complication groups tested (diabetic nephropathy n = 601, diabetic neuropathy n = 600, ophthalmic complications n = 601 and macrovascular complications n = 559), based on corresponding diagnosis or medical events experienced after the onset of T2DM. A phenotype-based distribution of subjects among different complication groups is shown in Fig. 1. After inspection of patient anthropometric data, we observed that T2DM patients experiencing ophthalmic complications and diabetic neuropathy were significantly younger compared to the group of patients without particular complications while those with diabetic nephropathy, macrovascular and ophthalmic complications were characterized by longer duration of diabetes. As expected, median glycated hemoglobin levels were higher in cases (patients with specific complications) compared to controls (patients without the same diabetes complication) except for the diabetic nephropathy group (Table 1).

Fig. 1
figure1

Venn diagram showing the distribution and overlap of cases among the analyzed T2DM complication groups

Table 1 Characteristics of the study participants

After quality control and filtering 5 378 539 SNPs were used for further testing in each T2DM complication group. The total genotyping rate was > 0.99 in all of the tested T2DM complication groups. The genomic inflation factor was negligible in all data sets based on median chi-squared statistics: 1.02 for neuropathies, 1.02 for macrovascular complications, 1.01 for ophthalmic complications, and 1.00 for nephropathies. After adjustments for age, sex, BMI, diabetes duration, median HbA1c, and medications used ten susceptibility loci were identified for different T2DM complications (Fig. 2), among them rs1132787 (GYPA) and rs522521 (LOC105371557) showed an association with diabetic neuropathy, rs2477088 (PDE4DIP), rs522521 (LOC105371557), rs4852954 (NAT8), rs6032 (F5), rs6935464 (RPS6KA2), rs7236163 (ZNF519), rs3095447 (CCDC146) were significantly associated with macrovascular complications, and only variant rs3095447 (CCDC146) was related to a greater risk of ophthalmic complications, while no significant hits were found for diabetic nephropathy (Table 2, Additional file 2).

Fig. 2
figure2

Manhattan plots for GWAS of T2DM complications. a Diabetic neuropathies, b macrovascular complications, c ophthalmic complications. X-axis shows chromosomal positions. Y-axis shows –log10 P-values. The red lines indicate a genome-wide significant threshold of P < 5 × 10−8, and the blue lines denote a suggestive significance threshold of P < 5 × 10−5. Association signals that reached genome-wide significance are denoted by reference SNP ID number

Table 2 Susceptibility loci associated with type 2 diabetes mellitus complications

The functional consequences of the identified variants were evaluated by the eQTL analysis, which was conducted in GTEx database by focusing only on those tissue types that may be relevant or even damaged according to the etiology of the disease [29]. Three out of eight novel allelic variants identified in this study appeared to be significantly associated with the expression of multiple nearby genes: rs4852954 (ALMS1, DUSP11, NAT8, ALMS1P1, TPRKB, ALMS1-IT1, RP11-434P11.2), rs7236163 (ZNF519, RP11-411B10.2), and rs3095447 (RP11-467H10.1, FGL2, PMS2P9, GSAP, SPDYE18, FAM185BP, CCDC146, UPK3BP1). See Additional file 3 for the full list of the identified associations.

Targeted analysis

In order to investigate if previously reported associations in other populations are also true in our study cohorts, we performed a targeted analysis in the same subjects and stratification in complication groups (diabetic neuropathy, macrovascular complications, ophthalmic complications, and diabetic nephropathy) (Table 1) by performing an association analysis for 194 candidate variants in total. Information on previously known allele-trait associations reported in both GWAS Catalog and GWAS Central was used for the selection of risk variants associated with at least one group of the tested T2DM complications (see Additional file 1 for a full list of selected candidate SNPs). The genotyping rate in all of the complication groups analyzed was > 0.99. By applying the targeted approach we managed to replicate two significant associations for diabetic neuropathy rs3761980 and rs80028505, both mapping to MAPK14 loci and one significant hit (rs136161, APOL1) for diabetic nephropathy (Table 3) in our study cohort.

Table 3 Candidate variants showing a significant association with type 2 diabetes complications

Discussion

Here we present the results of the genome-wide association study for T2DM complications performed in a population of Latvia for the first time, revealing 10 susceptibility loci for T2DM complications, including diabetic neuropathy, macrovascular and ophthalmic complications. As in other reports aimed to identify the risk factors of T2DM complications [15, 32], the control group of our study consisted of T2DM patients with no evidence of the complication type of interest instead of conventional healthy subjects, since the implementation of healthy controls would rather reveal genetic associations with the diagnosis of T2DM itself, not the T2DM complications.

We found two novel variants (rs1132787 and rs522521) associated with diabetic neuropathy and none of them have been linked to any disease or specific phenotype before. Variant rs1132787 is located within the 3′ UTR of a gene coding for glycophorin A coding (GYPA). Glycophorin A is the major erythrocyte membrane sialoglycoprotein. Although it has not been directly associated with susceptibility to any T2DM-related condition before, studies report a significant upregulation of the GYPA gene in the dorsal root ganglia of a mouse model of T2DM and the metabolic syndrome, and even type 1 diabetes with diabetic polyneuropathy [33, 34]. Moreover, copy number variation overlapping GYPA has been already linked to body mass index, obesity, and obesity-related traits, such as weight, hip circumference, and waist circumference, providing more evidence for the potential contribution of GYPA in the development of diabetic complications [35]. Because of the multiple evidence of the implication of GYPA gene in neuropathies and metabolic traits, and the potential functional consequence of the top variant rs1132787 located within the 3′UTR of the gene, we consider the GYPA gene as the first candidate for future functional validation studies. Nevertheless, credible evidence for this association should be established first by performing a replication study with larger sample size. The other identified risk variant rs522521 is located near the poorly characterized gene LOC105371557 with yet unknown function.

In total, seven variants reached genome-wide significance for the association with macrovascular complications of T2DM. The strongest association was exhibited by an intron variant of the Phosphodiesterase 4D Interacting Protein coding gene (PDE4DIP). Although the rs2477088 variant has not been previously linked to any T2DM manifestations, PDE4DIP, also known as myomegalin or cardiomyopathy-associated protein 2, is a well-known contributor of the microtubule control process [36] and some previous evidence exist indicating on the potential role of the gene in macrovascular diseases and T2DM. The exome sequencing has revealed a rare variant of PDE4DIP, which significantly increases the risk of ischemic stroke [37], moreover, CpG island methylation in leukocytes annotated to PDE4DIP contributes to the epigenetic fingerprint of myocardial infarction [38], and finally, the gene is also significantly downregulated in liver of T2DM patients [39]. Another risk allele (rs4852954) identified in the analysis of macrovascular complications is located near the N-Acetyltransferase 8 coding gene (NAT8) and has been previously associated with systolic blood pressure and renal function in the Estonian population [40] which is genetically close to the Latvian population [41]. Although NAT8 gene has not been linked with T2DM macrovascular complications before, it is considered to be a susceptibility locus for diabetic kidney disease [42]. In addition, we found the coagulation factor 5 coding gene (F5) among the risk loci of macrovascular complications. The variant rs6032 is located only around 7 kb from the Factor V Leiden (rs6025), which has been strongly associated with ischemic stroke and incident venous thrombosis before [43, 44], therefore we may speculate that in rs6032 carriers manifestation of the trait are amplified by the presence of T2DM. Although rs6032 is the only missense variant among the top hits, it is categorized as benign according to SIFT and PolyPhen. Nevertheless, the variant is in the LD with rs4524, the risk variant for venous thromboembolism (OR = 1.14; CI = 1.11–1.16) [45].

Finally, 4 more loci showed genome-wide significance for the association with macrovascular complications, among them rs6935464, located within the RPS6KA2 gene coding for Ribosomal Protein S6 Kinase A2, which is involved in cardiac myocyte stress responses and even considered as a therapeutic target for the prevention of heart failure [46], and rs3095447, an intron variant of Coiled-Coil Domain Containing 146 gene (CCDC146) which is also the only significant hit for ophthalmic complications in our data. The last two variants (rs522521 and rs7236163) identified in the analysis of macrovascular complications are located in the intergenic regions near genes LOC105371557 and ZNF519 respectively, and both have not been linked to any disease before.

We noticed two of the variants (rs522521, rs3095447) appearing among the top hits of multiple complication groups, which seems rational since a number of patients had experienced more than one type of complication, resulting in a notable overlap of patients among four tested phenotype groups (Fig. 1). This finding may be also explained by similar etiologic characteristics between microvascular diabetic complications involving small vessels (neuropathy, nephropathy, ophthalmic complications) and macrovascular complications related to large vessel damage, with chronic hyperglycemia as the main cause of all these comorbidities. Moreover, microvascular and macrovascular complications tend to be strongly interconnected, and the damage of small vessels may contribute to the manifestations of heart disease in diabetes [7], which coincides with our data showing only a small number of T2DM patients corresponding to one complication group only, while 18% (n = 109) of T2DM patients included in this study had experienced multiple manifestations during the observation period, corresponding to all four analyzed complication types.

According to Ensembl Variant Effect Predictor, most of the identified variants are intronic, though rs1132787, which is located in the enhancer site, and rs4852954 laying within the promoter flanking region may disrupt the functions of regulatory elements and therefore modulate the gene expression patterns. Some of the genes affected by the identified eQTLs (Additional file 3) are previously associated with the performance of the vascular system, for instance, expression of Fibrinogen-like protein 2 coding gene (FGL2) in endothelial cells has been previously linked to microthrombosis and cardiac impairment in rats with T2DM [47], and Alstrom Syndrome Protein 1 coding gene (ALMS1) is associated with Alström syndrome and characteristic dilated cardiomyopathy [48]. According to GTEx data, rs3095447 negatively correlates with FGL2 expression in heart and artery, while rs4852954 is linked to lower ALMS1 levels in whole blood and skin tissue, suggesting the possible functional implication of these variants in the pathogenesis of the disease. In order to explore additional effects of the top hits, quantitative phenotype analyses were performed, revealing significant associations between rs6935464, rs6032, rs3095447, rs4852954, rs7236163, rs522521, and median HbA1c levels. These data suggest that the identified genetic loci may serve as markers for both, development of T2DM complications and alterations in HbA1c levels, though the particular analysis should be repeated in a different cohort where the variable impact of different anti-diabetic medications should also be considered.

By performing targeted analysis of candidate variants we managed to replicate three associations with T2DM complications in the population of Latvia. Both variants (rs3761980 and rs80028505) showing an association with T2DM neuropathy are located near the Mitogen-activated protein kinase 14 coding gene (MAPK14) and have been previously linked to increased risk of diabetic foot ulcers in the report of The Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) project [49]. In addition, we found a significant association of rs136161 located in the Apolipoprotein L1 coding gene (APOL1) with T2DM nephropathy, which has been already linked to an advanced diabetic kidney disease across multiple ethnic groups [50]. Although we have observed smaller effect sizes for rs3761980 and rs80028505 comparing to other GWAS and much larger effect for rs136161, the previously reported odds ratio values for all three variants fit within the 95% confidence intervals calculated in our study [49, 50].

This study has several limitations, though the small sample size is the primary limiting factor in the risk variant discovery, which may result in an insufficient statistical power for the detection of rare variants with small effect sizes. This may also explain why variants with large effect sizes (OR up to 2.71) are markedly represented among the significant results of our study. Although the retrospective distribution of cases and controls in the analysis of T2DM neuropathy and ophthalmic complications was relatively balanced, the lack of significant hits in the analysis of diabetic nephropathy may be explained by the high incidence of the specific type of complication among the study participants (86%), leaving the size of the control group too small for the identification of true associations. Additionally, the suboptimal case-to-control ratio in the analysis of macrovascular complications may affect the statistical power of the study and explain the relatively high number of associations identified in the specific complication group. Due to the limited number of study subjects, the follow-up period was not fixed or set as inclusion criteria, though it was fitted as a covariate in the association analysis. Since the duration of diabetes is one of the strongest risk factors for the development of vascular complications [51], the use of an equal observation period in all cases would reduce the residual variability and improve the quality of this study.

Successful integration of genotyping data with longitudinal phenotypic information produced from several national health registries has provided strong support for 10 loci showing a genome-wide significance for the association with T2DM complications, some of them with already known importance to the comorbid conditions analyzed. We believe that these findings provide deeper insight into the pathogenesis of T2DM complications and suggest novel candidate genes for further functional studies, while our targeted approach highlights several susceptibility loci showing a directionally consistent impact on phenotype in multiple populations.

Conclusions

Using the genome-wide genotyping approach this study identified ten novel associations with T2DM complications, including GYPA (rs1132787) in diabetic neuropathy, PDE4DIP (rs2477088), NAT8 (rs4852954), F5 (rs6032), RPS6KA2 (rs6935464) in macrovascular complications, and CCDC146 (rs3095447) in ophthalmic complications. Meanwhile, the candidate gene analysis demonstrated a strong association for diabetic neuropathy (MAPK14: rs3761980, rs80028505), and diabetic nephropathy (APOL1: rs136161), proving the contribution of these risk loci in the pathogenesis of diabetic complications across various populations.

Availability of data and materials

The majority of the data generated and analysed during this study are included in this published article and its supplementary information files. The raw genotyping data are under restricted access from LGDB, nevertheless, they are available for research purposes from the corresponding author on reasonable request.

Abbreviations

BMI:

Body mass index

CI:

Confidence interval

eQTL:

Expression quantitative trait locus

GWAS:

Genome-wide association study

HbA1c:

Hemoglobin A1C

HR:

Hazard ratio

ICD-10:

10Th revision of the International Statistical Classification of Diseases and Related Health Problems.

LGDB:

Genome Database of the Latvian Population

MAF:

Minor allele frequency

OR:

Odds ratio

SD:

Standard deviation

SNP:

Single nucleotide polymorphism

T2DM:

Type 2 diabetes mellitus

UTR:

Untranslated region

References

  1. 1.

    Saeedi P, Petersohn I, Salpea P, Malanda B, Karuranga S, Unwin N, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9th edition. Diabetes Res Clin Pract. 2019;157:107843.

  2. 2.

    Murea M, Ma L, Freedman BI. Genetic and environmental factors associated with type 2 diabetes and diabetic vascular complications. Rev Diabet Stud. 2012;9(1):6–22.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Ali O. Genetics of type 2 diabetes. World J Diabetes. 2013;4(4):114–23.

    PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Poulsen P, Kyvik KO, Vaag A, Beck-Nielsen H. Heritability of type II (non-insulin-dependent) diabetes mellitus and abnormal glucose tolerance–a population-based twin study. Diabetologia. 1999;42(2):139–45.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Zhao W, Rasheed A, Tikkanen E, Lee JJ, Butterworth AS, Howson JMM, et al. Identification of new susceptibility loci for type 2 diabetes and shared etiological pathways with coronary heart disease. Nat Genet. 2017;49(10):1450–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Prudente S, Di Paola R, Pezzilli S, Garofolo M, Lamacchia O, Filardi T, et al. Pharmacogenetics of oral antidiabetes drugs: evidence for diverse signals at the IRS1 locus. Pharmacogenomics J. 2018;18(3):431–5.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Chawla A, Chawla R, Jaggi S. Microvasular and macrovascular complications in diabetes mellitus: Distinct or continuum? Indian J Endocrinol Metab. 2016;20(4):546–51.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    von Ferber L, Koster I, Hauner H. Medical costs of diabetic complications total costs and excess costs by age and type of treatment results of the German CoDiM Study. Exp Clin Endocrinol Diabetes. 2007;115(2):97–104.

    Article  CAS  Google Scholar 

  9. 9.

    Clarke P, Gray A, Legood R, Briggs A, Holman R. The impact of diabetes-related complications on healthcare costs: results from the United Kingdom Prospective Diabetes Study (UKPDS Study No. 65). Diabet Med. 2003;20(6):442–50.

  10. 10.

    Barr EL, Zimmet PZ, Welborn TA, Jolley D, Magliano DJ, Dunstan DW, et al. Risk of cardiovascular and all-cause mortality in individuals with diabetes mellitus, impaired fasting glucose, and impaired glucose tolerance: the Australian Diabetes, Obesity, and Lifestyle Study (AusDiab). Circulation. 2007;116(2):151–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Hietala K, Forsblom C, Summanen P, Groop PH, FinnDiane Study G. Heritability of proliferative diabetic retinopathy. Diabetes. 2008;57(8):2176–80.

    Article  CAS  Google Scholar 

  12. 12.

    Vionnet N, Tregouet D, Kazeem G, Gut I, Groop PH, Tarnow L, et al. Analysis of 14 candidate genes for diabetic nephropathy on chromosome 3q in European populations: strongest evidence for association with a variant in the promoter region of the adiponectin gene. Diabetes. 2006;55(11):3166–74.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Montesanto A, Bonfigli AR, Crocco P, Garagnani P, De Luca M, Boemi M, et al. Genes associated with Type 2 Diabetes and vascular complications. Aging (Albany NY). 2018;10(2):178–96.

    CAS  Article  Google Scholar 

  14. 14.

    Looker HC, Nelson RG, Chew E, Klein R, Klein BE, Knowler WC, et al. Genome-wide linkage analyses to identify Loci for diabetic retinopathy. Diabetes. 2007;56(4):1160–6.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Liu C, Chen G, Bentley AR, Doumatey A, Zhou J, Adeyemo A, et al. Genome-wide association study for proliferative diabetic retinopathy in Africans. NPJ Genom Med. 2019;4:20.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  16. 16.

    Jardeleza MS, Miller JW. Review of anti-VEGF therapy in proliferative diabetic retinopathy. Semin Ophthalmol. 2009;24(2):87–92.

    PubMed  Article  Google Scholar 

  17. 17.

    Cheung CY, Hui EY, Lee CH, Kwok KH, Gangwani RA, Li KK, et al. Impact of genetic loci identified in genome-wide association studies on diabetic retinopathy in Chinese patients with type 2 diabetes. Invest Ophthalmol Vis Sci. 2016;57(13):5518–24.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Wojcik GL, Graff M, Nishimura KK, Tao R, Haessler J, Gignoux CR, et al. Genetic analyses of diverse populations improves discovery for complex traits. Nature. 2019;570(7762):514–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Kochi Y, Suzuki A, Yamada R, Yamamoto K. Genetics of rheumatoid arthritis: underlying evidence of ethnic differences. J Autoimmun. 2009;32(3–4):158–62.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Moltke I, Grarup N, Jorgensen ME, Bjerregaard P, Treebak JT, Fumagalli M, et al. A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes. Nature. 2014;512(7513):190–3.

    CAS  PubMed  Article  Google Scholar 

  21. 21.

    Maghbooli Z, Pasalar P, Keshtkar A, Farzadfar F, Larijani B. Predictive factors of diabetic complications: a possible link between family history of diabetes and diabetic retinopathy. J Diabetes Metab Disord. 2014;13:55.

    PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Rovite V, Wolff-Sagi Y, Zaharenko L, Nikitina-Zake L, Grens E, Klovins J. Genome database of the latvian population (LGDB): design, goals, and primary results. J Epidemiol. 2018;28(8):353–60.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Marees AT, de Kluiver H, Stringer S, Vorspan F, Curis E, Marie-Claire C, et al. A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int J Methods Psychiatr Res. 2018;27(2):e1608.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Delaneau O, Zagury JF, Marchini J. Improved whole-chromosome phasing for disease and population genetic studies. Nat Methods. 2013;10(1):5–6.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012;44(8):955–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Beck T, Shorter T, Brookes AJ. GWAS Central: a comprehensive resource for the discovery and comparison of genotype and phenotype data from genome-wide association studies. Nucleic Acids Res. 2020;48(D1):D933–40.

    CAS  PubMed  Google Scholar 

  27. 27.

    Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47(D1):D1005–12.

    CAS  Article  Google Scholar 

  28. 28.

    Consortium GT. Erratum: Genetic effects on gene expression across human tissues. Nature. 2018;553(7689):530.

    Article  CAS  Google Scholar 

  29. 29.

    Forbes JM, Cooper ME. Mechanisms of diabetic complications. Physiol Rev. 2013;93(1):137–88.

    CAS  Article  Google Scholar 

  30. 30.

    Yates AD, Achuthan P, Akanni W, Allen J, Allen J, Alvarez-Jarreta J, et al. Ensembl 2020. Nucleic Acids Res. 2020;48(D1):D682–8.

    CAS  PubMed  Google Scholar 

  31. 31.

    Phan L,Jin Y,Zhang H, Qiang W, Shekhtman E, Shao D, Revoe D, Villamarin R, Ivanchenko E, Kimura M, Wang ZY, Hao L, Sharopova N, Bihan M, Sturcke A, Lee M, Popova N, Wu W, Bastiani C, Ward M, Holmes JB, Lyoshin V, Kaur K, Moyer E, Feolo M, Kattman BL. ALFA: Allele Frequency Aggregator: National Center for Biotechnology Information, U.S. National Library of Medicine; 2020 [10 Mar]. Available from: www.ncbi.nlm.nih.gov/snp/docs/gsr/alfa/.

  32. 32.

    Meng W, Shah KP, Pollack S, Toppila I, Hebert HL, McCarthy MI, et al. A genome-wide association study suggests new evidence for an association of the NADPH Oxidase 4 (NOX4) gene with severe diabetic retinopathy in type 2 diabetes. Acta Ophthalmol. 2018;96(7):e811–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Hur J, Dauch JR, Hinder LM, Hayes JM, Backus C, Pennathur S, et al. The metabolic syndrome and microvascular complications in a murine model of type 2 diabetes. Diabetes. 2015;64(9):3294–304.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Cheng C, Kobayashi M, Martinez JA, Ng H, Moser JJ, Wang X, et al. Evidence for epigenetic regulation of gene expression and function in chronic experimental diabetic neuropathy. J Neuropathol Exp Neurol. 2015;74(8):804–17.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Zhao W, Wineinger NE, Tiwari HK, Mosley TH, Broeckel U, Arnett DK, et al. Copy number variations associated with obesity-related traits in African Americans: a joint analysis between GENOA and HyperGEN. Obesity (Silver Spring). 2012;20(12):2431–7.

    CAS  Article  Google Scholar 

  36. 36.

    Roubin R, Acquaviva C, Chevrier V, Sedjai F, Zyss D, Birnbaum D, et al. Myomegalin is necessary for the formation of centrosomal and Golgi-derived microtubules. Biol Open. 2013;2(2):238–50.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Auer PL, Nalls M, Meschia JF, Worrall BB, Longstreth WT Jr, Seshadri S, et al. Rare and coding region genetic variants associated with risk of ischemic stroke: the NHLBI exome sequence project. JAMA Neurol. 2015;72(7):781–8.

    PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Ward-Caviness CK, Agha G, Chen BH, Pfeiffer L, Wilson R, Wolf P, et al. Analysis of repeated leukocyte DNA methylation assessments reveals persistent epigenetic alterations after an incident myocardial infarction. Clin Epigenet. 2018;10(1):161.

    CAS  Article  Google Scholar 

  39. 39.

    Li L, Pan Z, Yang X. Key genes and co-expression network analysis in the livers of type 2 diabetes patients. J Diabetes Investig. 2019;10(4):951–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Juhanson P, Kepp K, Org E, Veldre G, Kelgo P, Rosenberg M, et al. N-acetyltransferase 8, a positional candidate for blood pressure and renal regulation: resequencing, association and in silico study. BMC Med Genet. 2008;9:25.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  41. 41.

    Pliss L, Timsa L, Rootsi S, Tambets K, Pelnena I, Zole E, et al. Y-chromosomal lineages of latvians in the context of the genetic variation of the Eastern-Baltic region. Ann Hum Genet. 2015;79(6):418–30.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Kottgen A, Pattaro C, Boger CA, Fuchsberger C, Olden M, Glazer NL, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet. 2010;42(5):376–84.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. 43.

    Simone B, De Stefano V, Leoncini E, Zacho J, Martinelli I, Emmerich J, et al. Risk of venous thromboembolism associated with single and combined effects of Factor V Leiden, Prothrombin 20210A and Methylenetethraydrofolate reductase C677T: a meta-analysis involving over 11,000 cases and 21,000 controls. Eur J Epidemiol. 2013;28(8):621–47.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Casas JP, Hingorani AD, Bautista LE, Sharma P. Meta-analysis of genetic studies in ischemic stroke: thirty-two genes involving approximately 18,000 cases and 58,000 controls. Arch Neurol. 2004;61(11):1652–61.

    PubMed  Article  Google Scholar 

  45. 45.

    Lindstrom S, Wang L, Smith EN, Gordon W, van Hylckama VA, de Andrade M, et al. Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism. Blood. 2019;134(19):1645–57.

    PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Martinez EC, Passariello CL, Li J, Matheson CJ, Dodge-Kafka K, Reigan P, et al. RSK3: a regulator of pathological cardiac remodeling. IUBMB Life. 2015;67(5):331–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Ding Y, Liu K, Wang Y, Su G, Deng H, Zeng Q, et al. Expression and significance of fgl2 prothrombinase in cardiac microvascular endothelial cells of rats with type 2 diabetes. J Huazhong Univ Sci Technolog Med Sci. 2010;30(5):575–81.

    CAS  PubMed  Article  Google Scholar 

  48. 48.

    Hearn T. ALMS1 and Alstrom syndrome: a recessive form of metabolic, neurosensory and cardiac deficits. J Mol Med (Berl). 2019;97(1):1–17.

    CAS  Article  Google Scholar 

  49. 49.

    Meng W, Veluchamy A, Hebert HL, Campbell A, Colhoun HM, Palmer CNA. A genome-wide association study suggests that MAPK14 is associated with diabetic foot ulcers. Br J Dermatol. 2017;177(6):1664–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    Iyengar SK, Sedor JR, Freedman BI, Kao WH, Kretzler M, Keller BJ, et al. Genome-wide association and trans-ethnic meta-analysis for advanced diabetic kidney disease: family investigation of nephropathy and diabetes (FIND). PLoS Genet. 2015;11(8):e1005352.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  51. 51.

    Zoungas S, Woodward M, Li Q, Cooper ME, Hamet P, Harrap S, et al. Impact of age, age at diagnosis and duration of diabetes on the risk of macrovascular and microvascular complications and death in type 2 diabetes. Diabetologia. 2014;57(12):2465–74.

    PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

We thank the participants for their contribution and acknowledge the Genome Database of the Latvian Population for providing biological material and data.

Funding

The study was supported by European Regional Development Fund (ERDF) On Implementation of Activity 1.1.1.2 “Post-doctoral Research Aid” of the Specific Aid Objective 1.1.1 “To increase the research and innovative capacity of scientific institutions of Latvia and the ability to attract external financing, investing in human resources and infrastructure” of the Operational Programme “Growth and Employment”. Project No 1.1.1.2/VIAA/2/18/287 “Identification of clinical subgroups of Type 2 diabetes mellitus and application of pharmacogenetics in the development of personalized antidiabetic therapy”. The funding body was not involved in the design of the study, collection, analysis or interpretation of data, or writing the manuscript.

Author information

Affiliations

Authors

Contributions

MU, RP, VR, and JK made substantial contributions to the design of the work; RP, RR, LZ, and MU performed data analysis and interpretation; VR, IK, IE, LS, and LZ were involved in the DNA processing and acquisition of genotyping data; VP, JS, IE, LS, IK were involved in the acquisition of registry data, analysed clinical outcomes and established definition of phenotypes for analysis; MU and RP drafted the manuscript; VR, JK provided critical revision of the manuscript. All authors read and approved the final manuscript, and have agreed both to be personally accountable for the author's own contributions and to ensure that questions related to the accuracy or integrity of any part of the work, even ones in which the author was not personally involved, are appropriately investigated, resolved, and the resolution documented in the literature. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Janis Klovins.

Ethics declarations

Ethics approval and consent to participate

The study protocol is approved by the Central Medical Ethics Committee of Latvia (No. 01-29.1/2223). The study was conducted in accordance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) and The Convention for the protection of Human Rights and Dignity of the Human Being with regard to the Application of Biology and Medicine: Convention on Human Rights and Biomedicine. All study participants gave broad written informed consent for the recruitment in LGDB involving the use of their samples for research purposes. The acquisition of clinical data from the records of the national diabetes registry and Latvian hospital inpatient discharges was approved by the Latvian Centre for Disease Prevention and Control (No. 3/10.07.2017.) according to the regulations of Latvia.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1:

 List of all candidate SNPs selected for targeted association analysis.

Additional file 2:

Quantile-quantile plots for GWAS of T2DM complications.

Additional file 3:

List of significant associations obtained in eQTL analysis.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ustinova, M., Peculis, R., Rescenko, R. et al. Novel susceptibility loci identified in a genome-wide association study of type 2 diabetes complications in population of Latvia. BMC Med Genomics 14, 18 (2021). https://doi.org/10.1186/s12920-020-00860-4

Download citation

Keywords

  • Type 2 diabetes mellitus
  • Genome-wide genotyping
  • Diabetic complications