- 1Center for Biotechnology, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
- 2Department of Biomedical Engineering, College of Engineering, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
- 3Dubai Health Authority, Dubai Diabetes Center, Dubai, United Arab Emirates
- 4Division of Psychiatry, Faculty of Health and Medical Sciences, The University of Western Australia, Crawley, WA, Australia
- 5School of Medical and Health Sciences, Edith Cowan University, Joondalup, WA, Australia
- 6Fakeeh University Hospital, Dubai, United Arab Emirates
Background: Type 2 diabetes mellitus (T2DM) is a chronic, metabolic disorder in which concomitant insulin resistance and β-cell impairment lead to hyperglycemia, influenced by genetic and environmental factors. T2DM is associated with long-term complications that have contributed to the burden of morbidity and mortality worldwide. The objective of this manuscript is to conduct an Exome-Wide Association Study (EWAS) on T2DM Emirati individuals to improve our understanding on diabetes-related complications to improve early diagnostic methods and treatment strategies.
Methods: This cross-sectional study recruited 310 Emirati participants that were stratified according to their medically diagnosed diabetes-related complications: diabetic retinopathy, diabetic neuropathy, diabetic nephropathy, and cardiovascular complications. The Illumina’s Infinium Exome-24 array was used and 39,840 SNPs remained for analysis after quality control.
Findings: The analysis revealed the associations of various genes with each complication category: 1) diabetic retinopathy was associated to SHANK3 gene in locus 22q13.33 (SNP rs9616915; p=5.18 x10-4), ZSCAN5A gene in locus 19q13.43 (SNP rs7252603; p=7.55 x10-4), and DCP1B gene in locus 12p13.33 (SNPs rs715146, rs1044950, rs113147414, rs34730825; p=7.62 x10-4); 2) diabetic neuropathy was associated to ADH4 gene in locus 4q23 (SNP rs4148883; p=1.23 x10-4), SLC11A1 gene in locus 2q35 (SNP rs17235409; p=1.85 x10-4), and MATN4 gene in locus 20q13.12 (SNP rs2072788; p=2.68 x10-4); 3) diabetic nephropathy was associated to PPP1R3A gene in locus 7q31.1 (SNP rs1799999; p=1.91 x10-4), ZNF136 gene in locus 19p13.2 (SNP rs140861589; p=2.80 x10-4), and HSPA12B gene in locus 20p13 (SNP rs6076550; p=2.86 x10-4); and 4) cardiovascular complications was associated to PCNT gene in locus 21q22.3 (SNPs rs7279204, rs6518289, rs2839227, rs2839223; p=2.18 x10-4,3.04 x10-4,4.51 x10-4,5.22 x10-4 respectively), SEPT14 gene in locus 7p11.2 (SNP rs146350220; p=2.77 x10-4), and WDR73 gene in locus 15q25.2 (SNP rs72750868; p=4.47 x10-4).
Interpretation: We have identified susceptibility loci associated with each category of T2DM-related complications in the Emirati population. Given that only 16% of the markers from the Illumina’s Infinium Exome chip passed quality control assessment, this demonstrates that multiple variants were, either, monomorphic in the Arab population or were not genotyped due to the use of a Euro-centric EWAS array that limits the possibility of including targeted ethnic-specific SNPs. Our results suggest the alarming possibility that lack of representation in reference panels could inhibit discovery of functionally important loci associated to T2DM complications. Further effort must be conducted to improve the representation of diverse populations in genotyping and sequencing studies.
1 Introduction
Type 2 Diabetes Mellitus (T2DM) is a chronic, metabolic condition, characterized by elevated blood glucose levels (1). Although the pathogenesis of T2DM is complex, a number of factors that increase the risk for the disease have been identified, including modifiable risk factors (body mass index (BMI), physical inactivity, diet) and nonmodifiable risk factors (age, ethnicity, comorbid diseases, family history and genetic predisposition) (2). The clinical presentation and disease progression of patients with T2DM are heterogeneous, which may lead to a delay of diagnosis, multiple pathophysiological abnormalities, and varying susceptibility to complications. Complications from T2DM can be classified as microvascular complications, such as retinopathy, neuropathy and nephropathy, or macrovascular complications, including cardiovascular, cerebrovascular, and peripheral vascular disease (3). Although there is a strong inheritance of risk of developing T2DM, less is known about the heritability and genetic component of diabetes complications (4). Further studies must be conducted to elucidate the genetic variants associated to each diabetic complication to improve early diagnostic measures and therapeutic strategies.
Genome wide association studies (GWAS) has played a major role in identifying susceptibility loci associated with these various categories of diabetes-driven complications. More than 300 genetic loci have been associated with T2DM, which explain >19% of the phenotypic variance in risk for T2DM risk (5). Early family and twin studies have suggested a high concordance rate of the diabetic complications, with heritability estimated at 18 to 60% (6–10). GWAS studies have identified susceptible loci for diabetic retinopathy (WDR72, NVL, and CCDC146) (11–13), diabetic neuropathy (XIRP2, and APOL1) (13, 14), diabetic nephropathy (GABRR1, and GYPA) (7, 13), and cardiovascular complications (PDE4DIP, NAT8, F5, LPA, and RPS6KA2) (13, 15, 16). However, a number of the single nucleotide polymorphisms (SNPs) that failed to replicate in multiple populations demonstrate the strong influence of population specificity on genetic variation discrimination and contribution to the phenotype of interest. Therefore, discovery and replication investigations in populations of various ancestries are required to identify population-specific traits (17–19). This variability is the leading cause of clinical translation discrepancies due to the scarcity of genetic research specifically to the Middle East region, with multiple countries reporting a T2DM prevalence >20%, including Kuwait, Egypt and the United Arab Emirates (UAE) (20–22).
With the rising prevalence of diabetes-related complications, there is an urgency of conducting genetic studies to uncover new target pathways, and enhance our ability to use precision medicine for targeted therapeutic measures. By identifying new genotypes in an underrepresented region, in this case the UAE, this will yield to the discovery of novel genetic associations in diabetic-related complications. In this study, we aim to conduct an Exome wide association study (EWAS) to identify susceptibility loci associated with diabetic complication development within the Emirati population.
2 Methods
2.1 Ethics approval
An ethical request was submitted to the Dubai Health authority (DHA) whereby it was accepted under reference number DSREC-07/2020_19 and conducted in accordance with the Declaration of Helsinki. All participants provided written informed consent before taking part in this research. All data was de-identified prior to use.
2.2 Study group and phenotype definitions
This prospective, cross-sectional study recruited a total of 338 T2DM patients from the Dubai Diabetes Center (DDC), during the period between October of 2020 and July of 2021. All the patients were diagnosed in accordance to the American Diabetes Association (ADA) diagnosis criteria of a HbA1c ≥ 6.5 and were receiving treatment for their condition. To limit misclassification and ascertainment bias, the patient recruitment process was randomized for a more accurate representation of diabetes within the region.
The blood samples were collected in a sterile 5ml sample tube supplemented with ethylenediaminetetraacetic acid from the cubital vein. Samples were transported in a sealed biohazard bag using a cool transport container to Khalifa University, Center for Biotechnology, in Abu Dhabi for genotypic and analysis. The questionnaire included details on the demographic information, clinical details including physical measurements and medical status, medications prescribed, and biochemical parameters. In this questionnaire, it was ensured that the following clinical data was recorded: date of T2DM diagnosis, presence or absence of a diabetes-related complications, type of complication, and HbA1c measurements attained from the DHA’s Salama electronic medical record system. The patients with the presence of complications were stratified into four different phenotype-based categories: retinopathy, neuropathy, nephropathy, and cardiovascular complication. The group stratification was defined as follows:
1. Retinopathy complication: records of proliferative or non-proliferative retinopathy, or laser since the diagnosis of T2DM.
2. Neuropathy complication: records of foot ulcers, gangrene, amputation of the toe/foot/leg, pain in calf muscle while walking, shunting and angioplasty on artery in the leg since the diagnosis of T2DM.
3. Nephropathy complication: records of protein or albumin in the urine, albuminuria in the range of 30 – 299 mg/g, estimated Glomerular Filtration Rate (eGFR) <30 since the diagnosis of T2DM.
4. Cardiovascular complication: records of coronary artery bypass grafting or a cerebrovascular accident since the diagnosis of T2DM.
2.3 DNA extraction and genotyping
DNA extraction of 338 T2DM patients was conducted, as per the manufacturer’s instructions using the Qiagen DNA extraction kit. DNA samples were genotyped with the Infinium Exome BeadChip (Illumina, USA) scanned with the iScan System microarray scanner (Illumina, USA). This BeadChip has a total of 244,883 fixed markers. The raw data was uploaded onto GenomeStudio 2.0 and converted into PLINK format. Quality control (QC) was done to check for discordant gender information, missing genotype data (<98%), outlying heterozygosity rate (±3), and related individuals (PI_HAT>0.5). This led to the removal of 28 individuals (1 individual had low genotype quality and 27 individuals were related) for not passing the QC. The SNPs were filtered using the following parameters: low minor allele frequency (<0.01), low genotyping rate (<95%), and deviation from Hardy-Weinberg Equilibrium (p<10-6). The number of variants excluded for each filtering parameter was 202075 variants, 2946 variants, and 22 variants, respectively. A total of 310 individuals and 39,840 SNPs passed QC and remained for analysis.
2.4 Statistical analysis
Association analyses corresponding to the following four complication groups were conducted for descriptive statistics and genetic association (EWAS): retinopathy complications, neuropathy complications, nephropathy complications, and cardiovascular complications. For each category, the cases were those that were assigned to that category and the control group were all the remaining individuals that did not experience that particular complication. Statistical analysis of demographic characteristics and anthropometric measurements was conducted. Pearson χ2 was used to measure the association of categorical variables. Independent sample t-test, presented as mean and standard deviation, or nonparametric Mann-Whitney U-test, presented as median and inter-quartile region, were used to study continuous variables. Statistical analysis was performed in R (version 3.4), SPSS (version 46.0) and PLINK (version 1.9).
For the genetic case-control comparisons, logistic regression, assuming additive allelic effects for genotypes SNPs, were conducted, while adjusting for age, sex, and BMI. Exome-wide association markers surpassed a conservative Bonferroni-corrected significance threshold of discovery p<1.2×10-6 (0.05/39,840), whereas markers that identified associations that reached a suggestive association threshold of p<5×10-4. A quantile-quantile (Q-Q) plot analysis was conducted to check whether the distribution of the inflation p-values deviated from the expected distribution under the null hypothesis of no genetic association and the impact of population stratification was evaluated by calculating the genomic control inflation factor [λ GC]. A Manhattan plot was generated with -log10p-values. Q-Q plots and Manhattan plots were generated using the Locuszoom tool. Regional plots were generated by using LocusZoom.
3 Results
A cohort of 310 T2DM patients of which 153 were men and 157 were women aged 14 to 86 years. The cohort was stratified into cases or controls according to four complication groups that are to be tested: retinopathy complications (n=62), neuropathy complications (n=47), nephropathy complications (n=22), and cardiovascular complications (n=42). This classification was done according to diagnosis by the diabetes specialist after the onset of T2DM.
After assessing the anthropometric data of the study cohort (Table 1), it was seen that T2DM patients with neuropathy (p<0.001) and macrovascular (p<0.001) complications were significantly older than the control group. This indicates that T2DM-related complications are more likely to develop with age, providing us with the confidence to adjust for age during the analysis. The gender and mean BMI were not significantly different between cases and control, across all complications. The median glycated hemoglobin levels were significantly higher in the retinopathy cases (p=0.002) compared to controls. The complication groups retinopathy (p<0.001), neuropathy (p<0.001) and cardiovascular complications (p<0.001) were characterized with a longer diabetes duration as opposed to the nephropathy groups (p=0.058).
After performing QC and filtering, 39,840 SNPs were used for further testing in each category of T2DM complication. The total genotyping rate was > 0.995 across all categories. A quantile-quantile (Q-Q) plot analysis was carried out to check whether the distribution of the inflation p-values deviated from the expected distribution under the null hypothesis of no genetic association and investigate if the overall significance of the genome-wide associations is due to potential impact of population stratification. Supplementary Figure 1 presents the Q-Q plot of each respective complication, demonstrating that the genomic inflation factor was negligible in all data sets where it was 1.0 for all the categories based on the chi-squared statistics, after adjustment to age, BMI and gender. Figure 1 demonstrates the Manhattan plot of each complication, and the top 10 SNPs that contributed to the biological relevance of the respective disease is listed in Table 2.
Figure 1 Manhattan plot for diabetes-related (A) retinopathy complications (n=62), (B) neuropathy complications (n=47), (C) nephropathy complications (n=22), and (D) cardiovascular complications (n=42). The GWAS analyses results are shown on the y-axis as -log10(p-value) and on the x-axis is the chromosomal location, adjusted for age, gender, and BMI. The blue horizontal line illustrates the suggestive genome-wide association threshold (p<5×10-4).
Table 2 Top 10 SNPs that were associated with each diabetes-related complication group in the Emirati population.
3.1 Retinopathy complications
Gene ACVR1C is highly expressed in adipose tissue, and has been associated to extraocular retinoblastoma, hyperkeratosis, T2DM, obesity and anthropometric measurements, such as waist-to-hip ratio and body mass index (23–25). Interestingly, ACVRIC is also associated to lipid profile and glycemic markers (26–30). Similarly, gene ZFHX4 is associated to fasting blood glucose measurement and metabolite levels (31, 32). The association with pulse pressure and blood pressure have been associated to diabetic retinopathy through arterial stiffness and vision impairment, which has been identified in multiple genes, including the ZFHX4 gene (26, 33, 34), the SHANK3 gene (35), and the WNT9B gene (34). The SHANK3 gene, expressed in the brain, has also been associated to fibrinogen levels and platelet count, which has been reported to be risk factors in the development and progression of retinopathy (36–40).
The ZSCAN5A gene is expressed in the brain is associated with monocyte count, which may lead to the release of pro-inflammatory factors that interfere with endothelial cell junction integrity of the blood-retinal barrier, resulting in leucocyte infiltration in the retina (26, 37, 40, 41). The DCP1B gene, expressed in the brain, is associated with waist-to-hip ratio, BMI, and obesity-related traits, all risk factors of T2DM (42–44). This DCP1B gene is also associated with Insulin-like growth factors (IGFs), in which transgenic mice models that elucidated that overexpression of IGF-1 in the retina resulted in variations of eye-related diseases similar to that in diabetic humans, through retinal capillaries basement membrane thickening, venule dilation, intra-retinal microvascular abnormalities, and retinal and vitreous cavity neovascularization (44, 45).
3.2 Neuropathy complications
The GFY gene is mainly expressed in brain tissue, and has been associated to atherosclerosis through narrowing of the peripheral arterial vasculature (46). ADH4 gene, expressed in the liver tissue, is associated with eosinophil count, lipid measurements, Apolipoprotein A1 levels (ApoA-I), fibrinogen levels and factor VII levels (38, 39, 47–49). The association with fibrinogen is an important association, as fibrinogen participates in the coagulation process which may lead to an inflammatory process, inhibiting the growth of nerve axons and is closely related to diabetic neuropathy (50, 51). The LRFN2 gene is expressed in the brain, and has been associated to BMI, T2DM, and obesity-related traits (26, 42, 52, 53).
Interestingly, gene PKHD1 has been associated to intraocular pressure, brain measurement, T2DM, and metabolic markers, all risk factors associated to neuropathy (53–57). SLC11A1 gene is expressed in the bone marrow and lymphoid tissues, and has been associated to iron metabolism (58). Using a murine model, Iron’s effect on T2DM was elucidated demonstrating a positive association to motor nerve conduction velocities via a reduction in pro-inflammatory macrophages and an increase in anti-inflammatory macrophages in nerve sections may induce neuropathy (59). The MATN4 encodes a protein that is involved in filamentous networks in the extracellular matrices, which is essential for axonal health and growth and may lead to nerve fiber loss (60). The PPARA gene has been associated to immune and inflammatory responses, as well as lipid markers, glycolytic markers, T2DM and anthropometric measurements, such as waist-to-hip ratio and body mass index, all relevant risk factors for diabetic neuropathy (40, 61–64).
3.3 Nephropathy complications
The TTN gene in the skeletal muscle and has been associated to cardiac serum proteins and fractal structure of the heart, as well as T2DM and nephron-related variables, such glomerular filtration rate (65–67). While gene PI16, DPY6, FRMD4A and CROCC have not been reported to be associated with nephropathy, they have been identified in T2DM (53, 68) and obesity-related traits (69, 70). PPP1R3A gene is associated with T2DM and plays a crucial role in glycogen synthesis in the tubules of the kidney, leading to diabetic nephropathy.
The ZNF136 gene is highly expressed in the kidneys, and encodes a protein that contains a Krüppel-associated box (KRAB) A-box domain, which has been associated to the development of progressive chronic kidney disease (CKD). The Glis2, a Krüppel-like zinc finger protein, mutant mice had increased cell death and basement membrane thickening in the proximal convoluted tubules, resulting in severe renal atrophy with lymphocytic inflammatory cells infiltration and renal failure (71). The HSPA12B gene is expressed in the kidney and urinary bladder whose pathways are related to cellular senescence and cellular response to heat stress. This gene has been associated with gamma-glutamyl transferase (GGT) levels, a marker of oxidative stress that is linked with diabetes and hypertension, both being risk factors of CKD (72, 73).
3.4 Cardiovascular complications
The PKHD1 gene has been associated to T2DM (53), coronary artery disease (49, 74), cardiac troponin T levels (75), and obesity-related traits (44, 57, 76). While the MAST1 gene has not been associated to cardiovascular complications, it has been reported to be linked glycated hemoglobin levels (77). Importantly, gene GFY has been associated to carotid plaque build, leading to cardiovascular complications (46). The SEPT14 gene is expressed in the brain, heart, bone marrow, and lymphoid tissues, encoding a highly conserved septin family of cytoskeletal proteins that represses the accumulation of reactive oxygen species, resulting in cardiac microvascular endothelial cells apoptosis (78).
Multiple signals within the PCNT gene were identified. The PCNT gene is highly expressed in heart, and is an integral component of the microtubule-organizing proteins, which exert compressive forces on cardiomyocytes that drive the development of cardiac disorders and T2DM (53, 79). Interestingly, PCNT was also associated with cataract, indicating how microvascular and macrovascular complications tend to be strongly interrelated as damages of small vessels can ultimately results in heart disease manifestations in diabetes (80, 81). The RILPL2 is highly expressed in lymphocytic cells and artery, and have been associated to obesity-related traits (43, 70), including BMI and waist-to-hip ratio, as well as peripheral arterial disease (82).
4 Discussion
For the first time, we present the top markers identified from an exome-wide association study for T2DM-related complications conducted in the Emirati population. By identifying the susceptible loci associated to high-risk patients that develop complications form T2DM, this may improve targeted therapeutic interventions and early biomarker diagnosis through a panel of genetic markers. Most of the genes identified have been reported in other GWAS studies of different ethnicities, with a biological relevance to the pathogenesis of each respective complication group. These findings provide valuable insight into the pathogenesis of T2DM driven complications and suggest novel candidate genes for future functional studies.
As per the demographic characteristics, T2DM patients with neuropathy and macrovascular complications were significantly older, with a longer diabetes duration, than the control group. The gender and mean BMI were not significantly different between cases and control, across all complications. Interestingly, the median glycated hemoglobin levels was significantly higher in the retinopathy cases (p=0.002) compared to controls, which has been reported in other studies, possibly due to the formation of thrombus, a pathophysiological basis of early diabetic retinopathy (83).
When investigating sub-phenotypes of T2DM, diabetic retinopathy has been identified to be associated with ACVRIC (rs4664229), ZFHX4 (rs61729527), WNT9B (rs4968281), SHANK3 (rs9616915), ZSCAN5A (rs7252603), and DCP1B (rs715146, rs1044950, rs113147414) gene. These genes have intercrossing pathways and similar genetic variants to fibrinogen levels associated to intra-vessel pressure, low platelet count, leukocyte-retinal endothelial cell adhesion, metabolite levels and glycemic markers, all important factors impacting intra-retinal microvascular abnormalities, retinal capillaries and variations of eye-related diseases (23–30, 37, 40–44). For diabetic neuropathy, gene GFY (rs4802605), ADH4 (rs4148883), LRFN2 (rs61731010), PKHD1 (rs2499486), SLC11A1 (rs17235409), MATN4 (rs2072788), and PPARA (rs4253772) were associated or contributed to the biological relevance to the pathogenesis of the complication. Specifically, these markers have been associated to atherosclerosis, immune and inflammatory responses, AST and ApoA-I levels, iron toxicity, intraocular pressure, and compositional changes in extracellular matrices, which is essential for axonal health and growth, and may lead to nerve fiber loss in neuropathic conditions (26, 40, 42, 52–57, 61–64).
The genes that contributed to the biological relevance of diabetic nephropathy, include gene TTN (rs72646845), PI16 (rs113848006), DPY6 (rs36027551), CROCC (rs41272737), PPP1R3A (rs1799999), ZNF136 (rs140861589), HSPA12B (rs6076550), and FRMD4A (rs1541010). The markers identified to the development of diabetic nephropathy have mainly been expressed in the kidney and urine bladder, and have been associated to nephron-related variables, such glomerular filtration rate, glycogen synthesis in the tubules of the kidney and thickening in the proximal convoluted tubules (65–67, 71–73). Cardiovascular complications in T2DM is associated to PKHD1 (rs62406032), MAST1 (rs1078264), GFY (rs480265), SEPT14 (rs146350220), PCNT (rs6518289, rs2839227, rs2839223) and RILPL2 (rs28434767). Interestingly, these markers have been associated to coronary artery disease, glycated hemoglobin levels, cardiac troponin T levels, and obesity-related traits (44, 49, 57, 74–77).
The major limitation in this study is the sample size with an inadequate statistical power to be able to detect rare variants in the population pool. Moreover, the control group of the study included patients with a short duration of illness which could have contributed to a reduced power to the study. However, it is also important to note that the real period of T2DM is usually assumed to be longer than the clinically defined duration by at least several years due to a delay of diagnosis. Future studies, with a larger cohort, should adjust for duration of diabetes as it may serve as a genetic risk factor. Furthermore, the HbA1c levels were recorded only at one time point, at the time of recruitment, which could have been a limiting factor. Another limiting factor is the exome microarray chip where its incompatibility with the Middle Eastern population was seen in the fact that many variants were excluded after quality control due to the identification of monomorphic markers, homozygosity due to high consanguinity, and the accumulation of deleterious recessive alleles within the gene pool of the population. In fact, approximately 82.5% did not pass the MAF cut-off, demonstrating possible missed identification of pathogenic variants. Genetic variation in population arises from new mutations occurring through generations, in which changes in MAF may occur. This is due to genetic drift or differences in fitness levels conferred by different alleles in the presence of certain environment, including population bottleneck due to high consanguinity or migration (84).
Further studies need to be conducted in a large-scale, multi-ethnic cohort to replicate the findings of this study and substantiate our current knowledge of complications associated to T2DM. Given that only 16% of the markers from the Illumina’s Infinium Exome chip passed quality control assessment, this demonstrates that multiple variants were, either, monomorphic in the Arab population or were not genotyped due to the use of a Euro-centric EWAS array that limits the possibility of including targeted ethnic-specific SNPs. Our results suggest the alarming possibility that lack of representation in reference panels could inhibit discovery of functionally important loci associated to T2DM complications. Enabling global equity in the benefits of genomics will be vital for precision medicine initiatives, including risk prediction, development of therapies and implications for screening and diagnostics. Future work in diverse populations should focus on using unbiased approaches, unbiased marker discover and global genome references. This will be beneficial to better understand reproducibility and heterogeneity of effects among populations, improve the power to identify causal drivers of association signals, as well as important resources for fine-mapping of causal and rare variants.
This study has demonstrated that given that the majority of genetic studies, including the genotyping and sequencing panels, are developed based on the European ancestry, it has essentially deemed inapplicable to other ethnic groups. This foreshadows a near future where those genetic tests that are only valid for European descent be used as the blueprint for clinical applications for genetics, creating a skewed standard for ethnic minorities, such as the Middle East population. The scarcity of baseline genetic data is indicative of health inequalities that may be faced, further highlighting the urgency to ensure the inclusion of non-European descents in the genetic research movement. Hence, a microarray chip that is more inclusive to the Arab population needs to be developed and utilized to ensure that a wider spectrum of variants is included to detect rare SNPs associated within this region of the world. Further effort must be conducted to improve the representation of diverse populations in genotyping and sequencing studies to enable the unprecedented characterization of fine-scale genetic architecture and genetic susceptibilities to diseases, globally. This would allow for eventual delving into pharmacogenomics for the development of therapeutic strategies catered to the patient according to the complications experienced.
Data availability statement
The data presented in the study are deposited in the NCBI Gene expression omnibus database (accession number: GSE226084), and are available upon request from the corresponding author.
Ethics statement
An ethical request was submitted to the Dubai Health authority (DHA) whereby it was accepted under reference number DSREC-07/2020_19, and conducted in accordance with the Declaration of Helsinki. All participants provided written informed consent before taking part in this research. All data was de-identified prior to use. The patients/participants provided their written informed consent to participate in this study.
Author contributions
HA, GT, and AH conceived the project to study diabetes-related complications in the UAE. HA, GT, and AH conceived the central research questions for the EWAS data. AM and MM initiated the first draft of the manuscript. MM conducted the analysis of the manuscripts. HA, AM, MM analyses and constructed the Figures and Tables. AM and DA were responsible for the recruitment of the patients and collecting data for the study. AM carried out the laboratory assays used in the study. HA, MM and AM provided critical review during manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
This study has been funded by Abu Dhabi Education and Knowledge.
Acknowledgments
We thank the participants of the study for their generosity in providing samples to advance our understanding of T2DM-related complications. We would like to thank Dubai Diabetes Center Nursing and Research team for their assistance with patient recruitment, sample collection, processing, and data extraction. We are also grateful to Ms. Khayce Juma, Ms. Suna Nazar, and Ms. Hema Vurivi who assisted in processing the samples in the center for Biotechnology at Khalifa University.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2023.1143067/full#supplementary-material
References
1. DeFronzo RA, Ferrannini E, Groop L, Henry RR, Herman WH, Holst JJ, et al. Type 2 diabetes mellitus. Nat Rev Dis Primers (2015) 1(1):1–22. doi: 10.1038/nrdp.2015.19
2. Deshpande AD, Harris-Hayes M, Schootman M. Epidemiology of diabetes and diabetes-related complications. Phys Ther (2008) 88(11):1254–64. doi: 10.2522/ptj.20080020
3. Harding JL, Pavkov ME, Magliano DJ, Shaw JE, Gregg EW. Global trends in diabetes complications: a review of current evidence. Diabetologia (2019) 62(1):3–16. doi: 10.1007/s00125-018-4711-2
4. Cole JB, Florez JC. Genetics of diabetes mellitus and diabetes complications. Nat Rev Nephrol (2020) 16(7):377–90. doi: 10.1038/s41581-020-0278-5
5. Mansour Aly D, Dwivedi OP, Prasad RB, Käräjämäki A, Hjort R, Thangam M, et al. Genome-wide association analyses highlight etiological differences underlying newly defined subtypes of diabetes. Nat Genet (2021) 53(11):1534–42. doi: 10.1038/s41588-021-00948-2
6. Control D, Group CTR. Clustering of long-term complications in families with diabetes in the diabetes control and complications trial. Diabetes (1997) 46(11):1829. doi: 10.2337/diab.46.11.1829
7. Van Zuydam NR, Ahlqvist E, Sandholm N, Deshmukh H, Rayner NW, Abdalla M, et al. A genome-wide association study of diabetic kidney disease in subjects with type 2 diabetes. Diabetes (2018) 67(7):1414–27. doi: 10.2337/db17-0914
8. Leslie R, Pyke D. Diabetic retinopathy in identical twins. Diabetes (1982) 31(1):19–21. doi: 10.2337/diab.31.1.19
9. Simonson MA, Wills AG, Keller MC, McQueen MB. Recent methods for polygenic analysis of genome-wide data implicate an important effect of common variants on cardiovascular disease risk. BMC Med Genet (2011) 12(1):1–9. doi: 10.1186/1471-2350-12-146
10. Zdravkovic S, Wienke A, Pedersen N, Marenberg M, Yashin A, De Faire U. Heritability of death from coronary heart disease: a 36-year follow-up of 20 966 Swedish twins. J Internal Med (2002) 252(3):247–54. doi: 10.1046/j.1365-2796.2002.01029.x
11. 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 Genomic Med (2019) 4(1):1–7. doi: 10.1038/s41525-019-0094-7
12. Pollack S, Igo RP Jr., Jensen RA, Christiansen M, Li X, Cheng C-Y, et al. Multiethnic genome-wide association study of diabetic retinopathy using liability threshold modeling of duration of diabetes and glycemic control. Diabetes (2019) 68(2):441–56. doi: 10.2337/db18-0567
13. Ustinova M, Peculis R, Rescenko R, Rovite V, Zaharenko L, Elbere I, et al. Novel susceptibility loci identified in a genome-wide association study of type 2 diabetes complications in population of Latvia. BMC Med Genomics (2021) 14(1):1–11. doi: 10.1186/s12920-020-00860-4
14. Tang Y, Lenzini PA, Pop-Busui R, Ray PR, Campbell H, Perkins BA, et al. A genetic locus on chromosome 2q24 predicting peripheral neuropathy risk in type 2 diabetes: results from the ACCORD and BARI 2D studies. Diabetes (2019) 68(8):1649–62. doi: 10.2337/db19-0109
15. Fall T, Gustafsson S, Orho-Melander M, Ingelsson E. Genome-wide association study of coronary artery disease among individuals with diabetes: the UK biobank. Diabetologia (2018) 61(10):2174–9. doi: 10.1007/s00125-018-4686-z
16. Qi L, Qi Q, Prudente S, Mendonca C, Andreozzi F, Di Pietro N, et al. Association between a genetic variant related to glutamic acid metabolism and coronary heart disease in individuals with type 2 diabetes. Jama (2013) 310(8):821–8. doi: 10.1001/jama.2013.276305
17. 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. doi: 10.1038/s41586-019-1310-4
18. 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. doi: 10.1016/j.jaut.2009.02.020
19. Mousa M, Vurivi H, Kannout H, Uddin M, Alkaabi N, Mahboub B, et al. Genome-wide association study of hospitalized COVID-19 patients in the united Arab Emirates. EBioMedicine (2021) 74:103695. doi: 10.1016/j.ebiom.2021.103695
20. Al-Rifai RH, Majeed M, Qambar MA, Ibrahim A, AlYammahi KM, Aziz F. Type 2 diabetes and pre-diabetes mellitus: a systematic review and meta-analysis of prevalence studies in women of childbearing age in the middle East and north Africa, 2000–2018. Syst Rev (2019) 8(1):1–32. doi: 10.1186/s13643-019-1187-1
21. Hamoudi R, Saheb Sharif-Askari N, Saheb Sharif-Askari F, Abusnana S, Aljaibeji H, Taneera J, et al. Prediabetes and diabetes prevalence and risk factors comparison between ethnic groups in the united Arab Emirates. Sci Rep (2019) 9(1):1–7. doi: 10.1038/s41598-019-53505-7
22. Sun H, Saeedi P, Karuranga S, Pinkepank M, Ogurtsova K, Duncan BB, et al. IDF diabetes atlas: Global, regional and country-level diabetes prevalence estimates for. World (2021) 965:1053.7. doi: 10.1016/j.diabres.2021.109119
23. Akbari P, Sosina OA, Bovijn J, Landheer K, Nielsen JB, Kim M, et al. Multiancestry exome sequencing reveals INHBE mutations associated with favorable fat distribution and protection from diabetes. Nat Commun (2022) 13(1):1–17. doi: 10.1038/s41467-022-32398-7
24. Justice AE, Karaderi T, Highland HM, Young KL, Graff M, Lu Y, et al. Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution. Nat Genet (2019) 51(3):452–69. doi: 10.1038/s41588-018-0334-2
25. Wang S-H, Su M-H, Chen C-Y, Lin Y-F, Feng Y-CA, Hsiao P-C, et al. Causality of abdominal obesity on cognition: A trans-ethnic Mendelian randomization study. Int J Obes (2022) 46(8):1487–92. doi: 10.1038/s41366-022-01138-8
26. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet (2021) 53(10):1415–24. doi: 10.1038/s41588-021-00931-x
27. Nielsen JB, Rom O, Surakka I, Graham SE, Zhou W, Roychowdhury T, et al. Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease. Nat Commun (2020) 11(1):1–12. doi: 10.1038/s41467-020-20086-3
28. Mozaffarian D, Kabagambe EK, Johnson CO, Lemaitre RN, Manichaikul A, Sun Q, et al. Genetic loci associated with circulating phospholipid trans fatty acids: a meta-analysis of genome-wide association studies from the CHARGE consortium. Am J Clin Nutr (2015) 101(2):398–406. doi: 10.3945/ajcn.114.094557
29. Lee C-J, Chen T-H, Lim AMW, Chang C-C, Sie J-J, Chen P-L, et al. Phenome-wide analysis of Taiwan biobank reveals novel glycemia-related loci and genetic risks for diabetes. Commun Biol (2022) 5(1):1–9. doi: 10.1038/s42003-022-04168-0
30. Klarin D, Damrauer SM, Cho K, Sun YV, Teslovich TM, Honerlaw J, et al. Genetics of blood lipids among~ 300,000 multi-ethnic participants of the million veteran program. Nat Genet (2018) 50(11):1514–23. doi: 10.1038/s41588-018-0222-9
31. Kong S, Cho YS. Identification of female-specific genetic variants for metabolic syndrome and its component traits to improve the prediction of metabolic syndrome in females. BMC Med Genet (2019) 20(1):1–13. doi: 10.1186/s12881-019-0830-y
32. Rhee EP, Ho JE, Chen M-H, Shen D, Cheng S, Larson MG, et al. A genome-wide association study of the human metabolome in a community-based cohort. Cell Metab (2013) 18(1):130–43. doi: 10.1016/j.cmet.2013.06.013
33. Yamamoto M, Fujihara K, Ishizawa M, Osawa T, Kaneko M, Ishiguro H, et al. Pulse pressure is a stronger predictor than systolic blood pressure for severe eye diseases in diabetes mellitus. J Am Heart Assoc (2019) 8(8):e010627. doi: 10.1161/JAHA.118.010627
34. Hoffmann TJ, Ehret GB, Nandakumar P, Ranatunga D, Schaefer C, Kwok P-Y, et al. Genome-wide association analyses using electronic health records identify new loci influencing blood pressure variation. Nat Genet (2017) 49(1):54–64. doi: 10.1038/ng.3715
35. Fujisawa T, Ikegami H, Yamato E, Kawaguchi Y, Ueda H, Shintani M, et al. Association of plasma fibrinogen level and blood pressure with diabetic retinopathy, and renal complications associated with proliferative diabetic retinopathy, in type 2 diabetes mellitus. Diabetic Med (1999) 16(6):522–6. doi: 10.1046/j.1464-5491.1999.00111.x
36. Keşkek NŞ, Gülcan H, Yılmaz G, Akkoyun İ. Impact of platelet count in retinopathy of prematurity. Turkish J Ophthalmol (2020) 50(6):351. doi: 10.4274/tjo.galenos.2020.01058
37. Vuckovic D, Bao EL, Akbari P, Lareau CA, Mousas A, Jiang T, et al. The polygenic and monogenic basis of blood traits and diseases. Cell (2020) 182(5):1214–31. e11. doi: 10.1016/j.cell.2020.08.008
38. Temprano-Sagrera G, Sitlani CM, Bone WP, Martin-Bornez M, Voight BF, Morrison AC, et al. Multi-phenotype analyses of hemostatic traits with cardiovascular events reveal novel genetic associations. J Thromb Haemostasis (2022) 20(6):1331–49. doi: 10.1111/jth.15698
39. De Vries PS, Chasman DI, Sabater-Lleal M, Chen M-H, Huffman JE, Steri M, et al. A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration. Hum Mol Genet (2016) 25(2):358–70. doi: 10.1093/hmg/ddv454
40. Chen M-H, Raffield LM, Mousas A, Sakaue S, Huffman JE, Moscati A, et al. Trans-ethnic and ancestry-specific blood-cell genetics in 746,667 individuals from 5 global populations. Cell (2020) 182(5):1198–213. e14. doi: 10.1016/j.cell.2020.06.045
41. Wan H, Cai Y, Wang Y, Fang S, Chen C, Chen Y, et al. The unique association between the level of peripheral blood monocytes and the prevalence of diabetic retinopathy: a cross-sectional study. J Trans Med (2020) 18(1):1–9. doi: 10.1186/s12967-020-02422-9
42. Pulit SL, Stoneman C, Morris AP, Wood AR, Glastonbury CA, Tyrrell J, et al. Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum Mol Genet (2019) 28(1):166–74. doi: 10.1093/hmg/ddy327
43. Zhu Z, Guo Y, Shi H, Liu C-L, Panganiban RA, Chung W, et al. Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK biobank. J Allergy Clin Immunol (2020) 145(2):537–49. doi: 10.1016/j.jaci.2019.09.035
44. Comuzzie AG, Cole SA, Laston SL, Voruganti VS, Haack K, Gibbs RA, et al. Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population. PloS One (2012) 7(12):e51954. doi: 10.1371/journal.pone.0051954
45. Ruberte J, Ayuso E, Navarro M, Carretero A, Nacher V, Haurigot V, et al. Increased ocular levels of IGF-1 in transgenic mice lead to diabetes-like eye disease. J Clin Invest (2004) 113(8):1149–57. doi: 10.1172/JCI19478
46. Pott J, Burkhardt R, Beutner F, Horn K, Teren A, Kirsten H, et al. Genome-wide meta-analysis identifies novel loci of plaque burden in carotid artery. Atherosclerosis (2017) 259:32–40. doi: 10.1016/j.atherosclerosis.2017.02.018
47. Yu B, Zheng Y, Alexander D, Morrison AC, Coresh J, Boerwinkle E. Genetic determinants influencing human serum metabolome among African americans. PloS Genet (2014) 10(3):e1004212. doi: 10.1371/journal.pgen.1004212
48. Richardson TG, Leyden GM, Wang Q, Bell JA, Elsworth B, Davey Smith G, et al. Characterising metabolomic signatures of lipid-modifying therapies through drug target mendelian randomisation. PloS Biol (2022) 20(2):e3001547. doi: 10.1371/journal.pbio.3001547
49. Kichaev G, Bhatia G, Loh P-R, Gazal S, Burch K, Freund MK, et al. Leveraging polygenic functional enrichment to improve GWAS power. Am J Hum Genet (2019) 104(1):65–75. doi: 10.1016/j.ajhg.2018.11.008
50. Coban E, Sari R, Ozdogan M, Akçit F. Levels of plasma fibrinogen and d-dimer in patients with impaired fasting glucose. Exp Clin Endocrinol Diabetes (2005) 113(01):35–7. doi: 10.1055/s-2004-830524
51. Zhuang Y, Lin X, Chen X, Wu X, Zhang J. Fibrinogen function indexes are potential biomarkers of diabetic peripheral neuropathy. Diabetol Metab Syndrome (2022) 14(1):1–7. doi: 10.1186/s13098-021-00777-7
52. Mahajan A, Spracklen CN, Zhang W, Ng MC, Petty LE, Kitajima H, et al. Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation. Nat Genet (2022) 54(5):560–72. doi: 10.1038/s41588-022-01058-3
53. Vujkovic M, Keaton JM, Lynch JA, Miller DR, Zhou J, Tcheandjieu C, et al. Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis. Nat Genet (2020) 52(7):680–91. doi: 10.1038/s41588-020-0637-y
54. Naqvi S, Sleyp Y, Hoskens H, Indencleef K, Spence JP, Bruffaerts R, et al. Shared heritability of human face and brain shape. Nat Genet (2021) 53(6):830–9. doi: 10.1038/s41588-021-00827-w
55. Ozel AB, Moroi SE, Reed DM, Nika M, Schmidt CM, Akbari S, et al. Genome-wide association study and meta-analysis of intraocular pressure. Hum Genet (2014) 133(1):41–57. doi: 10.1007/s00439-013-1349-5
56. Choquet H, Thai KK, Yin J, Hoffmann TJ, Kvale MN, Banda Y, et al. A large multi-ethnic genome-wide association study identifies novel genetic loci for intraocular pressure. Nat Commun (2017) 8(1):1–9. doi: 10.1038/s41467-017-01913-6
57. Martin S, Cule M, Basty N, Tyrrell J, Beaumont RN, Wood AR, et al. Genetic evidence for different adiposity phenotypes and their opposing influences on ectopic fat and risk of cardiometabolic disease. Diabetes (2021) 70(8):1843–56. doi: 10.2337/db21-0129
58. Wyllie S, Seu P, Goss JA. The natural resistance-associated macrophage protein 1 Slc11a1 (formerly Nramp1) and iron metabolism in macrophages. Microbes Infect (2002) 4(3):351–9. doi: 10.1016/S1286-4579(02)01548-4
59. Paeschke S, Baum P, Toyka KV, Blüher M, Koj S, Klöting N, et al. The role of iron and nerve inflammation in diabetes mellitus type 2-induced peripheral neuropathy. Neuroscience (2019) 406:496–509. doi: 10.1016/j.neuroscience.2019.03.005
60. Hill R. Extracellular matrix remodelling in human diabetic neuropathy. J Anatomy (2009) 214(2):219–25. doi: 10.1111/j.1469-7580.2008.01026.x
61. Barton AR, Sherman MA, Mukamel RE, Loh P-R. Whole-exome imputation within UK biobank powers rare coding variant association and fine-mapping analyses. Nat Genet (2021) 53(8):1260–9. doi: 10.1038/s41588-021-00892-1
62. Thareja G, Belkadi A, Arnold M, Albagha OM, Graumann J, Schmidt F, et al. Differences and commonalities in the genetic architecture of protein quantitative trait loci in European and Arab populations. Hum Mol Genet (2022) 32(6):907–16. doi: 10.1093/hmg/ddac243
63. Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, Ganna A, et al. Discovery and refinement of loci associated with lipid levels. Nat Genet (2013) 45(11):1274. doi: 10.1038/ng.2797
64. Chen J, Sun M, Adeyemo A, Pirie F, Carstensen T, Pomilla C, et al. Genome-wide association study of type 2 diabetes in Africa. Diabetologia (2019) 62(7):1204–11. doi: 10.1007/s00125-019-4880-7
65. Aung N, Vargas JD, Yang C, Fung K, Sanghvi MM, Piechnik SK, et al. Genome-wide association analysis reveals insights into the genetic architecture of right ventricular structure and function. Nat Genet (2022) 54(6):783–91. doi: 10.1038/s41588-022-01083-2
66. Meyer HV, Dawes TJ, Serrani M, Bai W, Tokarczuk P, Cai J, et al. Genetic and functional insights into the fractal structure of the heart. Nature (2020) 584(7822):589–94. doi: 10.1038/s41586-020-2635-8
67. Winkler TW, Justice AE, Graff M, Barata L, Feitosa MF, Chu S, et al. The influence of age and sex on genetic associations with adult body size and shape: a large-scale genome-wide interaction study. PloS Genet (2015) 11(10):e1005378. doi: 10.1371/journal.pgen.1005378
68. Pipal KV, Mamtani M, Patel AA, Jaiswal SG, Jaisinghani MT, Kulkarni H. Susceptibility loci for type 2 diabetes in the ethnically endogamous Indian Sindhi population: A pooled blood genome-wide association study. Genes (2022) 13(8):1298. doi: 10.3390/genes13081298
69. Tachmazidou I, Süveges D, Min JL, Ritchie GR, Steinberg J, Walter K, et al. Whole-genome sequencing coupled to imputation discovers genetic signals for anthropometric traits. Am J Hum Genet (2017) 100(6):865–84. doi: 10.1016/j.ajhg.2017.04.014
70. Christakoudi S, Evangelou E, Riboli E, Tsilidis KK. GWAS of allometric body-shape indices in UK biobank identifies loci suggesting associations with morphogenesis, organogenesis, adrenal cell renewal and cancer. Sci Rep (2021) 11(1):1–18. doi: 10.1038/s41598-021-89176-6
71. Kim Y-S, Kang HS, Herbert R, Beak JY, Collins JB, Grissom SF, et al. Kruppel-like zinc finger protein Glis2 is essential for the maintenance of normal renal functions. Mol Cell Biol (2008) 28(7):2358–67. doi: 10.1128/MCB.01722-07
72. Pazoki R, Vujkovic M, Elliott J, Evangelou E, Gill D, Ghanbari M, et al. Genetic analysis in European ancestry individuals identifies 517 loci associated with liver enzymes. Nat Commun (2021) 12(1):1–12. doi: 10.1038/s41467-021-22338-2
73. Cain L, Shankar A, Ducatman AM, Steenland K. The relationship between serum uric acid and chronic kidney disease among Appalachian adults. Nephrol Dialysis Transplantation (2010) 25(11):3593–9. doi: 10.1093/ndt/gfq262
74. Van Der Harst P, Verweij N. Identification of 64 novel genetic loci provides an expanded view on the genetic architecture of coronary artery disease. Circ Res (2018) 122(3):433–43. doi: 10.1161/CIRCRESAHA.117.312086
75. Yang Y, Bartz TM, Brown MR, Guo X, Zilhao NR, Trompet S, et al. Identification of functional genetic determinants of cardiac troponin T and I in a multiethnic population and causal associations with atrial fibrillation. Circulation: Genomic Precis Med (2021) 14(6):e003460. doi: 10.1161/CIRCGEN.121.003460
76. Koskeridis F, Evangelou E, Said S, Boyle JJ, Elliott P, Dehghan A, et al. Pleiotropic genetic architecture and novel loci for c-reactive protein levels. Nat Commun (2022) 13(1):1–11. doi: 10.1038/s41467-022-34688-6
77. Sinnott-Armstrong N, Tanigawa Y, Amar D, Mars N, Benner C, Aguirre M, et al. Genetics of 35 blood and urine biomarkers in the UK biobank. Nat Genet (2021) 53(2):185–94. doi: 10.1038/s41588-020-00757-z
78. Peng C, Ma J, Gao X, Tian P, Li W, Zhang L. High glucose induced oxidative stress and apoptosis in cardiac microvascular endothelial cells are regulated by FoxO3a. PloS One (2013) 8(11):e79739. doi: 10.1371/journal.pone.0079739
79. Becker R, Leone M, Engel FB. Microtubule organization in striated muscle cells. Cells (2020) 9(6):1395. doi: 10.3390/cells9061395
80. Hsu C-C, Chuang H-K, Hsiao Y-J, Teng Y-C, Chiang P-H, Wang Y-J, et al. Polygenic risk score improves cataract prediction in East Asian population. Biomedicines (2022) 10(8):1920. doi: 10.3390/biomedicines10081920
81. Chawla A, Chawla R, Jaggi S. Microvasular and macrovascular complications in diabetes mellitus: distinct or continuum? Indian J Endocrinol Metab (2016) 20(4):546. doi: 10.4103/2230-8210.183480
82. Ward-Caviness CK, Neas LM, Blach C, Haynes CS, LaRocque-Abramson K, Grass E, et al. Genetic variants in the bone morphogenic protein gene family modify the association between residential exposure to traffic and peripheral arterial disease. PloS One (2016) 11(4):e0152670. doi: 10.1371/journal.pone.0152670
83. Zhang B, Zhang B, Zhou Z, Guo Y, Wang D. The value of glycosylated hemoglobin in the diagnosis of diabetic retinopathy: a systematic review and meta-analysis. BMC Endocrine Disord (2021) 21(1):1–11. doi: 10.1186/s12902-021-00737-2
Keywords: T2DM, diabetes, EWAS, retinopathy, nephropathy, neuropathy, macrovascular complications, microvascular complications
Citation: Mansour A, Mousa M, Abdelmannan D, Tay G, Hassoun A and Alsafar H (2023) Microvascular and macrovascular complications of type 2 diabetes mellitus: Exome wide association analyses. Front. Endocrinol. 14:1143067. doi: 10.3389/fendo.2023.1143067
Received: 12 January 2023; Accepted: 02 March 2023;
Published: 23 March 2023.
Edited by:
Sen Li, Beijing University of Chinese Medicine, ChinaReviewed by:
Jacklyn Hellwege, Vanderbilt University Medical Center, United StatesKyle Flippo, The University of Iowa, United States
Copyright © 2023 Mansour, Mousa, Abdelmannan, Tay, Hassoun and Alsafar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Habiba Alsafar, aGFiaWJhLmFsc2FmYXJAa3UuYWMuYWU=
†These authors have contributed equally to this work