Next Article in Journal
Genotyping by RAD Sequencing Analysis Assessed the Genetic Distinctiveness of Experimental Lines and Narrowed down the Genomic Region Responsible for Leaf Shape in Endive (Cichorium endivia L.)
Next Article in Special Issue
Differentially Expressed Genes of Natural Killer Cells Can Distinguish Rheumatoid Arthritis Patients from Healthy Controls
Previous Article in Journal
Identification of Genetic Locus Underlying Easy Dehulling in Rice-Tartary for Easy Postharvest Processing of Tartary Buckwheat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrative Phenotype–Genotype Approach Using Phenotypic Characteristics from the UAE National Diabetes Study Identifies HSD17B12 as a Candidate Gene for Obesity and Type 2 Diabetes

1
Sharjah Institute for Medical Research, University of Sharjah, Sharjah 27273, UAE
2
Department of Clinical Sciences, College of Medicine, University of Sharjah, Sharjah 27272, UAE
3
Department of Clinical Sciences, Division of Islets Cell Pathology, Lund University, SE-205 02 Malmö, Sweden
4
Department of Basic Medical Sciences, College of Medicine, University of Sharjah, Sharjah 27272, UAE
5
Department of Family and Community Medicine and Behavioral Sciences, College of Medicine, University of Sharjah, Sharjah 27272, UAE
6
Baker Heart and Diabetes Institute, Melbourne 3004, Australia
*
Authors to whom correspondence should be addressed.
M.Y.H. and H.A. Authors equally contributed to this work.
Genes 2020, 11(4), 461; https://doi.org/10.3390/genes11040461
Submission received: 4 February 2020 / Revised: 15 April 2020 / Accepted: 21 April 2020 / Published: 23 April 2020
(This article belongs to the Special Issue Translational Bioinformatics: From Prediction to Validation)

Abstract

:
The United Arab Emirates National Diabetes and Lifestyle Study (UAEDIAB) has identified obesity, hypertension, obstructive sleep apnea, and dyslipidemia as common phenotypic characteristics correlated with diabetes mellitus status. As these phenotypes are usually linked with genetic variants, we hypothesized that these phenotypes share single nucleotide polymorphism (SNP)-clusters that can be used to identify causal genes for diabetes. We explored the National Human Genome Research Institute-European Bioinformatics Institute Catalog of Published Genome-Wide Association Studies (NHGRI-EBI GWAS) to list SNPs with documented association with the UAEDIAB-phenotypes as well as diabetes. The shared chromosomal regions affected by SNPs were identified, intersected, and searched for Enriched Ontology Clustering. The potential SNP-clusters were validated using targeted DNA next-generation sequencing (NGS) in two Emirati diabetic patients. RNA sequencing from human pancreatic islets was used to study the expression of identified genes in diabetic and non-diabetic donors. Eight chromosomal regions containing 46 SNPs were identified in at least four out of the five UAEDIAB-phenotypes. A list of 34 genes was shown to be affected by those SNPs. Targeted NGS from two Emirati patients confirmed that the identified genes have similar SNP-clusters. ASAH1, LRP4, FES, and HSD17B12 genes showed the highest SNPs rate among the identified genes. RNA-seq analysis revealed high expression levels of HSD17B12 in human islets and to be upregulated in type 2 diabetes (T2D) donors. Our integrative phenotype-genotype approach is a novel, simple, and powerful tool to identify clinically relevant potential biomarkers in diabetes. HSD17B12 is a novel candidate gene for pancreatic β-cell function.

1. Introduction

Type 2 diabetes (T2D) is a multifactorial disorder characterized by the insufficiency of insulin secretion and/or insulin action [1]. The disease is one of the most critical public health challenges of the 21st century [2]. Rapid socio-economic transition is one of the leading causes of the global diabetes epidemic and rising prevalence rates in developing economies [3]. The cost associated with the disease affects more impoverished world regions as well as high-income countries, thus imposing a substantial global economic burden [4]. Recent reports showed a high prevalence of T2D in the United Arab Emirates (UAE) [5]. Thus, effective interventions are urgently needed to slow the diabetes epidemic and reduce diabetes-related complications. An ultimate goal in diabetes management is the identification of novel biomarkers that enable the detection, prevention, treatment of the disease, and its complications long before overt disease development [6].
Genome-wide association studies (GWAS) have identified more than 143 common genetic variants associated with T2D [7]. However, these variants explain only a small proportion of the heritability of the disease [8]. Thus, more extensive studies are needed to identify T2D loci in different populations, as well as novel approaches to make sense of these generated databases [9].
The United Arab Emirates National Diabetes and Lifestyle Study (UAEDIAB) study is a cross-sectional survey designed to investigate the prevalence of diabetes and associated risk factors in 827 Emiratis and 2724 expatriates living in Dubai, Sharjah, and the Northern Emirates [10,11,12]. UAEDIAB study has identified a list of patient’s related measurements including obesity/BMI, hypertension, obstructive sleep apnea, and dyslipidemia that correlated with diabetes. Interestingly, most of these identified phenotypes have previously been associated with genetic variants [13,14,15,16].
Random genetic variants are unlikely to occur in the same genomic position unless there is a positive selection for that mutation, which might represent a fitness advantage [17]. However, genetic variants in terms of single nucleotide polymorphism (SNP) can cluster in particular chromosomal regions [18]. Those regions might represent an increased vulnerability to mutation as a result of their unique sequence, structural, and functional features [19]. Such SNPs clusters are thought to occur preferentially around some functional genes. Therefore, the identification of clustered SNPs might provide insights into the molecular mechanisms of mutagenesis. This may include vulnerable regions, possible positive selection for disease development, and regulatory actions associated with T2D.
In this study, we describe an integrative phenotype-genotype approach to identify novel potential biomarkers for T2D by utilizing the phenotypic characteristics of the UAEDIAB study. This was addressed by a systems genetics approach integrating SNP-clusters in chromosomal regions between the identified UAEDIAB-phenotypes and published SNPs using the National Human Genome Research Institute-European Bioinformatics Institute Catalog of Published Genome-Wide Association Studies (NHGRI-EBI GWAS) Catalog. We believe that our study might lead to novel testable hypotheses for genes involved in β-cell function and will add new insight into the molecular pathogenesis and biomarkers for T2D.

2. Materials and Methods

2.1. Study Design, Population, and Settings

The details of the study design and sampling of the UAEDIAB study are described elsewhere [10,20]. In brief, the UAEDIAB is a cross-sectional study conducted to estimate both prevalence and risk factors of diabetes among UAE nationals and expatriates who have been living in Sharjah, Dubai, and the Northern Emirates for at least four years. The study was approved by the Ethics Committee of Sharjah University and the Ministry of Health Research Ethics Committee (REC number MOHAP/DXB/SUBC/No.14/2017). All participants had consented to the usage of their collected data and samples.

2.2. Identification of SNP-Clusters Associated with the UAEDIAB-Phenotypes

We utilized the NHGRI-EBI GWAS Catalog (https://www.ebi.ac.uk/gwas/, v. 1.0) of published genome-wide association studies [21]. The file downloaded was (gwas_catalog_v1.0-associations_e96_r2019-07-12). The entire catalog was downloaded and prioritized for SNPs associated with the 5 given phenotypes; (1) diabetes (type 1 diabetes (T1D) and T2D), (2) obesity (obesity and body mass index), (3) dyslipidemia (triglyceride levels, triglycerides, lipid traits, high-density lipoprotein (HDL) cholesterol-triglycerides (HDLC-TG), low-density lipoprotein (LDL) cholesterol levels, HDL cholesterol levels, and total cholesterol levels), (4) obstructive sleep apnea (snoring, obstructive sleep apnea, obstructive sleep apnea trait (apnea–hypopnea index), (average oxygen saturation during sleep), or (average respiratory event duration), and (5) hypertension (systolic blood pressure ≥ 140/diastolic blood pressure ≥ 90. SNPs with a significant association with each of the given phenotypes associated with diabetes (adjusted p-value < 0.05) were identified along with their corresponding chromosomal regions.

2.3. SNP to SNP Functional and Pathway Enrichment

To investigate whether the identified genes share common pathways that may link the given phenotypes, we used eXploring Genomic Relations (XGR) for enhanced interpretation web tools (http://galahad.well.ox.ac.uk:3040/). XGR provides enhanced interpretation of GWAS via comprehensively utilizing ontology and network information [22]. The tool uses the Experimental Factor Ontology (EFO) option to provide a systematic description of many experimental variables available in EBI databases and the NHGRI-EBI Catalog [23]. Enrichment analysis is based on the hypergeometric/binomial distribution or Fisher′s exact test which tests the statistical significance of the observed number of SNPs overlapped between an input group of SNPs and SNPs annotated by EFO. As validation of the results, Enriched Ontology Clustering for the identified genes and SNPs associated with the 5 identified UAEDIAB-phenotypes study were generated using the Metascape (a web-based tool used for comprehensive gene list annotation and analysis resource) and Reactome option in SNPnexus [24]. Socialiser for SNPs in XGR was used to assess the degree of relatedness between any two SNPs in terms of annotation profiles and network were performed by the tool.
Similarity functions serve to conduct similarity analysis calculating semantic similarity—a type of comparison to assess the degree of relatedness between two entities (e.g., genes or SNPs) based on their annotation profiles (by ontology terms) [25]. To do so, information content (IC) of a term is first defined to measure how informative a term is to being used for annotating genes: −log10 (frequency of genes annotated to this term). The similarity between the two terms is then measured based on IC, usually at the most informative common ancestor (MICA). Finally, the similarity between two entities (e.g., genes) was derived from pairwise term similarity using best-matching based methods: average, maximum, and complete. To identify phenotypically important SNPs, we used SNPnexus (http://www.snp-nexus.org) to assess the potential significance of identified SNPs on the major transcriptome, proteome, regulatory, and structural variation models [26].

2.4. Targeted Next-Generation Sequencing

As proof of concept to validate the identified SNP-clusters and genes related to diabetes in our cohort, two Emirati diabetic patients (1 male and 1 female; age 60 ± 1) were selected as they have an exceptionally high level of glycated hemoglobin (HbA1c > 10.5%), fasting blood glucose values > 13 mmol/L, BMI > 40, weight 120 ± 3 kg, with dyslipidemia profile, hypertension, and obstructive sleep apnea/snoring. Extracted blood DNA was subjected to targeted DNA next-generation sequencing using the S5 semi-conductor based DNA sequencer as described previously [26]. Briefly, we used the targeted coding-exome sequencing approach by designing multiplex primers that map across key genes along the coding region of the genome. We multiplexed the two patients on one Ion Chip. A pooled barcoded amplicon-tagged library generated using Fluidigm Access Array (Fluidigm Europe B.V, Amsterdam, The Netherlands) was diluted and subjected to emulsion polymerase chain reaction (PCR) with Ion SphereTM particles with Ion Template OT2 200 Kit (Ion OneTouchTM system) following the manufacturer’s instructions (Thermo Fisher, Waltham, MA, USA). The pooled samples were sequenced using 540 Ion Chip with Ion 200 Sequencing kit on the S5 sequencer following manufacturer instructions (Thermo Fisher). The total mapped reads obtained were around 40 million per patient, with 109 mean depth and 94% uniformity coverage. 92.4% of the amplicons were aligned to the reference genome (HG19 build). The AQ20 quality mean length was 138 bp and the mean raw accuracy was 98.6%. The mean coverage across the amplicon was around ×130. The bioinformatics analysis was carried out by first aligning the data using the BWA alignment algorithm, followed by sequence filtering using SAMtools. Mutations file (VCF) was generated using vcftools. The mutations were visualized using Integrated Genome Viewer.

2.5. RNA-Seq from Human Islets

RNA-seq data were extracted from human pancreatic islets obtained from 75 cadaver donors (nondiabetic HbA1c < 6%, n = 63 and T2D/hyperglycemic HbA1c ≥ 6.3%, n = 12) as previously described [27,28]. RNA-seq was performed using Illumina’s TruSeq RNA Sample Preparation Kit as described previously [27,28]. Data normalization was processed using a trimmed mean of M-values and presented as Fragments/Kilobase of Exon Per Million Fragments Mapped (FPKM) or transformed into log2 counts per million using the voom-function (edgeR/limma R-packages). RNA-seq are deposited in a MIAME database (GSE50398).

2.6. In Silico Validation

To link the genomic SNPs to the functional and dynamic mRNA transcriptomics of the identified genes, we searched publicly available transcriptomics database “GEO Omnibus” of pancreatic islets isolated from nondiabetic and compared to diabetic samples. Then, we extracted blood transcriptomic datasets from T2D, pre-diabetic, and healthy controls. Details of the used datasets are listed in Supplementary Table S2.

2.7. Statistical Methods

GraphPad Prism version 7.00 for Windows (GraphPad Software, La Jolla, CA, USA) was used for statistical analysis. First, the D′Agostino–Pearson normality test was used to determine whether to perform parametric or nonparametric tests. One-way ANOVA test was performed to determine whether there are any statistically significant differences between the mean values of the controls and different groups for the gene expression and protein levels. For nonparametric tests, the Kruskal–Wallis test and Dunn’s correction for multiple testing were used. The same software was used to examine the correlations between the different parameters using linear regression. A Student’s t-test was used to look for the difference between two groups under a given experiment or treatment. p < 0.05 was considered significant.

3. Results

3.1. NHGRI-EBI GWAS Catalog Analysis Identifies Eight SNP-Clusters in Chromosome Regions with Frequent SNPs Associated with Most of the UAEDIAB-Phenotypes

SNPs related to the different UAEDIAB-phenotypes and diabetes were explored using the NHGRI-EBI GWAS Catalog. Our analysis showed a total of 708 chromosomal regions associated with most SNPs that have been associated with UAEDIAB-phenotypes. Of them, 211 chromosomal regions were associated with diabetes, 270 SNPs with BMI and obesity, 155 regions with lipid profile, and 48 regions with hypertension, while 24 regions were associated with obstructive sleep apnea (Figure 1).
Interestingly, the intersection of chromosomal regions associated with UAEDIAB-phenotypes with diabetes identified eight common chromosomal regions between the phenotypes. Of them, four chromosomal regions (8p22, 1q32.3, 12q24.13, and 7p15.2) were overlapped between BMI/obesity, dyslipidemia, hypertension, and T1D/T2D, three regions (11p11.2, 6q21, and 17q12) between BMI/obesity, dyslipidemia, sleep apnea, and T1D/T2D and only, one common region (15q26.1) was found between BMI/obesity, hypertension, sleep apnea, and T1D/T2D (Figure 2 and Table 1). Analysis of those eight shared chromosomal regions revealed 46 potential SNPs (Table 2). Data from the NHGRI-EBI GWAS Catalog revealed that 34 genes were reported to be affected by those SNPs (Table 2).

3.2. SNPs Enrichment Analysis Showed Significant Ontology Similarity and Metabolic Pathways Enrichment

XGR web tool was used to explore the correlation of the 46 identified SNPs as a group to each other or to UAEDIAB-phenotypes. This analysis further confirms the association between the identified 46 SNPs and the given phenotypes (FDR < 0.05). Additionally, our analysis showed a significant association between the 46 identified SNPs and other phenotypes (Table 3). Of them, 15 SNPs were associated with metabolic disease, 12 SNPs with diabetes mellitus, 13 SNPs with body mass index, and 10 SNPs with T2D.
Next, we harnessed the Socialiser tool to assess the degree of relatedness or network between any two SNPs in the meaning of annotation profiles. Our results showed that 25 out of the 46 SNPs had a similarity score of more than 0.5 in scale range from 0 to 1 with at least five other SNPs (Figure 3 and Table 4). These results indicate a possible co-occurrence or functional relationship of related SNPs. Moreover, clustering of the reported 34 genes in the shared chromosomal regions according to pathway enrichment revealed a significant enrichment of four major pathways, namely developmental biology, signal transduction, and metabolism of lipids and steroids (Supplementary Table S1).

3.3. Targeted Next-Generation Sequencing (NGS) Confirms the Presence of Similar SNP-Clusters in Emirati Diabetic Patients

To confirm that our pipeline was able to identify similar SNP-clusters in the Emirati population, we performed NGS on two diabetic patients (as described in the methods section). As shown in Table 5 and Supplementary Figure S1, NGS analysis confirmed the presence of a high rate of SNPs (defined as more than one SNP) in the same genes identified earlier.
Out of the 34 identified genes by our pipeline, four genes (ASAH1, LRP4, FES, and HSD17B12) showed to have at least four SNPs from each patient (four was the median of the SNPs number per patients after excluding those who have one SNP only).

3.4. RNA-Seq Analysis in Human Pancreatic Islets Showed that HSD17B12 is Novel Candidate Gene for Pancreatic β Cell Function

To gain more insights about the role of the four genes (ASAH1, LRP4, FES, and HSD17B12) in pancreatic islets, we used our published RNA-seq data to investigate their expression in human islets [27,28]. RNA-seq mean expression analysis from nondiabetic islets showed that ASAH1 and HSD17B12 are highly expressed at a high level in human islets as compared FES, LRP4, or to the ion channel gene KCNJ11, a functional marker for pancreatic β-cell function. (Figure 4A). Differential expression analysis exhibited that expression of HSD17B12 is significantly reduced (p = 0.03) in diabetic islets (n = 12) when compared with nondiabetic islets (n = 63, Figure 4E). Expression of ASAH1, FES, and LRP4 was not affected by diabetes status (Figure 4B–D). Collectively, based on the expression level and differential expression in human islets, our results suggest that HSD17B12 is a candidate gene for pancreatic β-cell function.

3.5. Annotations of SNPs in HSD17B12

Targeted NGS performed in two local Emirati diabetic patients identified a total of four intronic SNPs in the HSD17B12 gene (Supplementary Table S2). Interestingly, the two patients showed one common genotype of rs4573668 (homozygote C/C instead of the reference G). To identify if this SNP is phenotypically important, we used the SNPnexus tool to explore the location and functional consequences of rs4573668. Interestingly, we found that rs4573668 to be located in the 5’ untranslated region of the HSD17B12 transcript which could affect its expression, while it does not alter the amino acid sequence of the mature protein (Supplementary Figure S2). This might indicate that rs4573668 plays a role in controlling the activity and function of the HSD17B12 in terms of transcription or epigenetic modifications.

4. Discussion

In this study, we employed an integrative phenotype-genotype systems approach to investigate whether UAEDIAB-phenotypes shared chromosomal regions containing SNP-clusters, which can be exploited as candidate biomarkers for pancreatic β-cell function and risk of diabetes. Our analysis identified eight shared chromosomal regions in four UAE-DIAB-phenotypes (Figure 2 and Table 1). The identified regions contain 46 SNPs with reported effect on 34 genes (Table 2), hence, putatively involved in the pathogenesis processes leading to the disease. In vivo validation on Emirati diabetic patients, showed that ASAH1, LRP4, FES, and HSD17B12 have the highest number of SNPs among the 34 genes.
An important part of our study is providing more knowledge for the association of the five identified candidate genes with diabetes in human pancreatic islets from donors with and without diabetes (Figure 4). ASAH1 gene is a member of the ceramidases family which degrade ceramide to free fatty acids and sphingosine (SPH) [29]. It has been shown that the inhibitory effects of accumulated saturated fatty acids on insulin signaling were prevented by ASAH1 overexpression [30]. Other studies have revealed that endurance exercise was associated with a reduction in ceramide levels in the skeletal muscles of obese individuals, hence leading to the improvement of insulin sensitivity [31]. Although the expression of ASAH1 in human islets was high, we were not able to observe any differential expression of the gene in diabetic islets (Figure 4). This suggests that ASAH1 is probably not involved in β-cell function but rather insulin signaling as evident by other reports [30,31].
LRP4 is a single-pass transmembrane protein belonging to the LDL receptor-related protein (LRP) family. Members of this family were found to be involved in various processes, including development and physiology [32]. LRP4 loss of function has been associated with developmental anomalies associated with Cenani–Lenz syndrome (CLS) disease, which includes limb malformation and renal agenesis [33]. Furthermore, adipocytes in LRP4 deficient mice were characterized by an improvement of glucose and insulin tolerance, lipid homeostasis, less adipocyte hypertrophy as well as reduction of serum fatty acids, thus confirming the role of LRP4 in regulating glucose metabolism [34]. FES gene encodes for a non-receptor protein with tyrosine-specific activity. A recent investigation of the genetic determinants of cardiovascular diseases has found FES to be associated with hypertension and BMI [35]. The findings that expression of LPR4 and FES is not affected by diabetes status and low expressed in human islets indicate that those genes might not be potential important players in β-cell function (Figure 4).
HSD17B12 (Hydroxysteroid 17-β dehydrogenase) is an essential molecule in the elongation of very-long-chain fatty acids (VLCFAs) and the production of arachidonic acid through the conversion of 17-keto and 17-hydroxysteroids [36]. The gene has been recently linked to BMI and obesity susceptibility [37]. We showed that HSD17B12 expression was very high in human islets and significantly decreased in diabetic islets compared to nondiabetic (Figure 4). The findings may suggest a role of HSD17B12 gene in pancreatic β-cell function. In support of our RNA-seq data, it was reported that tissues involved in lipid metabolisms such as liver, kidney, muscle, and adipose tissues have a high expression level of HSD17B12 in humans and mice [36]. Interestingly, pathway enrichment analyses of the 34 genes pointed out that HSD17B12 is involved in signal transduction and metabolism of lipids and steroids, which fits well with studies of phenotypes.
A previous study showed that individuals carrying the T2D risk allele T for the intronic SNP rs11037579 had lower expression of HSD17B12 in adipose tissue of insulin-resistant subjects [38]. Additionally, rs1061810 was documented to be associated with T2D indicating a role for HSD17B12 in diabetes [39]. Our results indicated that genotyping of rs4573668 (C/C) is common in the local Emirati patients. Since the rs4573668 variant might lead to a change in the protein amino acid sequence, therefore, it could play a significant role in controlling the activity and function of the HSD17B12 in terms of transcription, translation, and epigenetic modifications. Moreover, the findings that HSD17B12 expression is reduced in diabetic islets, highlights the potential of utilization of this gene as an early biomarker in blood samples for pre-diabetic patients.
In summary, our study proposes HSD17B12 as a causal gene for T2D and pancreatic β-cell function. More future functional studies still needed to validate the finding.

5. Conclusions

In conclusion, our combined approach supports the existence of common chromosomal regions and SNP-clusters among obesity/BMI, hypertension, obstructive sleep apnea, dyslipidemia, and diabetes which might be involved in the pathogenesis of these clinically related phenotypes.
We were able to replicate these SNP-clusters on the local Emirati population different from the ethnic groups available in GWAS. HSD17B12 was identified as a candidate gene for β-cell function. Our data are the implementation of genomics-based approaches to chronic disease detection using bioinformatics in silico approach, and, therefore, further functional validation is still needed to elucidate the role of HSD17B12 and rs4573668 in the pathogenesis of T2D. In the end, we believe that such knowledge could increase our understanding of diabetes and facilitate the development of drugs.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/4/461/s1, Figure S1: Screenshots of examples SNPs in the genes identified generated by targeted DNA next-generation sequencing in Emirati diabetic patients, Figure S2: List of the location and functional consequences of rs4573668 as predicted by SNPnexus tool, Table S1: Enriched Ontology Clusters of the identified 34 genes associated with the UAEDIA-phenotypes, Table S2: List of the 4 identified SNPs in HSD17B12 gene in two Emirati diabetic patients, with their location and significance in gene expression.

Author Contributions

Conceptualization, M.Y.H., N.S. and J.T.; Methodology, M.Y.H., H.A., A.S. and J.T.; Software, M.Y.H., I.Y.H., R.A.H.; Writing—Original Draft Preparation, M.Y.H., H.A. and J.T.; Writing—Review & Editing, M.Y.H., N.M.E., A.K.M., N.S. and J.T.; Supervision, N.S., R.A.H. and J.T.; Project Administration, A.K.M., H.A.; Funding Acquisition, N.S. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

J.T. is supported by grants from the University of Sharjah (1701090119-P and 1701090121-P) and the Al-Jalila Foundation (AJF201723). R.A.H. is funded by Al-Jalila Foundation (AJF201741) and the Sharjah Research Academy (Grant code: MED001) and University of Sharjah (Grant code: 1901090254). N.S. is funded by the Al-Jalila Foundation (AJF201726), Diabetes and Metabolic Syndrome Research Group (University of Sharjah, Code: 150306), National Diabetes and Lifestyle Survey (MOH Code: 120301).

Acknowledgments

We would like to thank Alagappan Kumarappan in sample processing and Thenmozhi Venkatachalam for running the targeted NGS.

Conflicts of Interest

The authors declare no conflicts of interest in this work.

References

  1. Tripathi, B.K.; Srivastava, A.K. Diabetes mellitus: Complications and therapeutics. Med. Sci. Monit. 2006, 12, Ra130–Ra147. [Google Scholar] [PubMed]
  2. Zimmet, P.; Alberti, K.G.; Magliano, D.J.; Bennett, P.H. Diabetes mellitus statistics on prevalence and mortality: Facts and fallacies. Nat. Rev. Endocrinol. 2016, 12, 616–622. [Google Scholar] [CrossRef] [PubMed]
  3. Pradeepa, R.; Mohan, V. Prevalence of type 2 diabetes and its complications in India and economic costs to the nation. Eur. J. Clin. Nutr. 2017, 71, 816–824. [Google Scholar] [CrossRef] [PubMed]
  4. Bommer, C.; Heesemann, E.; Sagalova, V.; Manne-Goehler, J.; Atun, R.; Barnighausen, T.; Vollmer, S. The global economic burden of diabetes in adults aged 20–79 years: A cost-of-illness study. Lancet Diabetes Endocrinol. 2017, 5, 423–430. [Google Scholar] [CrossRef]
  5. Meo, S.A.; Usmani, A.M.; Qalbani, E. Prevalence of type 2 diabetes in the Arab world: Impact of GDP and energy consumption. Eur. Rev. Med. Pharmacol. Sci. 2017, 21, 1303–1312. [Google Scholar]
  6. Lyons, T.J.; Basu, A. Biomarkers in diabetes: Hemoglobin A1c, vascular and tissue markers. Transl. Res. J. Lab. Clin. Med. 2012, 159, 303–312. [Google Scholar] [CrossRef] [Green Version]
  7. Xue, A.; Wu, Y.; Zhu, Z.; Zhang, F.; Kemper, K.E.; Zheng, Z.; Yengo, L.; Lloyd-Jones, L.R.; Sidorenko, J.; Wu, Y. Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes. Nat. Commun. 2018, 9, 1–14. [Google Scholar] [CrossRef] [Green Version]
  8. Billings, L.K.; Florez, J.C. The genetics of type 2 diabetes: What have we learned from GWAS? Ann. N. Y. Acad. Sci. 2010, 1212, 59–77. [Google Scholar] [CrossRef]
  9. Flannick, J.; Florez, J.C. Type 2 diabetes: Genetic data sharing to advance complex disease research. Nat. Rev. Genet. 2016, 17, 535–549. [Google Scholar] [CrossRef]
  10. Sulaiman, N.; Albadawi, S.; Abusnana, S.; Fikri, M.; Madani, A.; Mairghani, M.; Alawadi, F.; Zimmet, P.; Shaw, J. Novel approach to systematic random sampling in population surveys: Lessons from the United Arab Emirates National Diabetes Study (UAEDIAB). J. Diabetes 2015, 7, 642–648. [Google Scholar] [CrossRef] [Green Version]
  11. Sulaiman, N.; Albadawi, S.; Abusnana, S.; Mairghani, M.; Hussein, A.; Al Awadi, F.; Madani, A.; Zimmet, P.; Shaw, J. High prevalence of diabetes among migrants in the United Arab Emirates using a cross-sectional survey. Sci. Rep. 2018, 8, 6862. [Google Scholar] [CrossRef] [PubMed]
  12. Sulaiman, N.; Elbadawi, S.; Hussein, A.; Abusnana, S.; Madani, A.; Mairghani, M.; Alawadi, F.; Sulaiman, A.; Zimmet, P.; Huse, O.; et al. Prevalence of overweight and obesity in United Arab Emirates Expatriates: The UAE National Diabetes and Lifestyle Study. Diabetol. Metab. Syndr. 2017, 9, 88. [Google Scholar] [CrossRef] [Green Version]
  13. Beaney, K.E.; Ward, C.E.; Bappa, D.A.S.; McGale, N.; Davies, A.K.; Hirani, S.P.; Li, K.; Howard, P.; Vance, D.R.; Crockard, M.A.; et al. A 19-SNP coronary heart disease gene score profile in subjects with type 2 diabetes: The coronary heart disease risk in type 2 diabetes (CoRDia study) study baseline characteristics. Cardiovasc. Diabetol. 2016, 15, 141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cheng, M.; Mei, B.; Zhou, Q.; Zhang, M.; Huang, H.; Han, L.; Huang, Q. Computational analyses of obesity associated loci generated by genome-wide association studies. PLoS ONE 2018, 13, e0199987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Li, W.; Zhu, L.; Huang, H.; He, Y.; Lv, J.; Li, W.; Chen, L.; He, W. Identification of susceptible genes for complex chronic diseases based on disease risk functional SNPs and interaction networks. J. Biomed. Inf. 2017, 74, 137–144. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, H.; Wang, W.; Zhang, C.; Xu, C.; Duan, H.; Tian, X.; Zhang, D. Heritability and Genome-Wide Association Study of Plasma Cholesterol in Chinese Adult Twins. Front. Endocrinol. 2018, 9. [Google Scholar] [CrossRef] [PubMed]
  17. Desai, M.M.; Fisher, D.S. Beneficial mutation–selection balance and the effect of linkage on positive selection. Genetics 2007, 176, 1759–1798. [Google Scholar] [CrossRef] [Green Version]
  18. Calarco, L.; Barratt, J.; Ellis, J. Genome wide identification of mutational hotspots in the apicomplexan parasite Neospora caninum and the implications for virulence. Genome Biol. Evol. 2018, 10, 2417–2431. [Google Scholar] [CrossRef]
  19. Amos, W. Even small SNP clusters are non-randomly distributed: Is this evidence of mutational non-independence? Proc. R. Soc. B Biol. Sci. 2010, 277, 1443–1449. [Google Scholar] [CrossRef] [Green Version]
  20. Hamoudi, R.; Saheb Sharif-Askari, N.; Saheb Sharif-Askari, F.; Abusnana, S.; Aljaibeji, H.; Taneera, J.; Sulaiman, N. Prediabetes and diabetes prevalence and risk factors comparison between ethnic groups in the United Arab Emirates. Sci. Rep. 2019, 9, 17437. [Google Scholar] [CrossRef]
  21. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; McMahon, A.; Morales, J.; Mountjoy, E.; Sollis, E.; et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Fang, H.; Knezevic, B.; Burnham, K.L.; Knight, J.C. XGR software for enhanced interpretation of genomic summary data, illustrated by application to immunological traits. Genome Med. 2016, 8, 129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Malone, J.; Holloway, E.; Adamusiak, T.; Kapushesky, M.; Zheng, J.; Kolesnikov, N.; Zhukova, A.; Brazma, A.; Parkinson, H. Modeling sample variables with an Experimental Factor Ontology. Bioinformatics 2010, 26, 1112–1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 2019, 10, 1523. [Google Scholar] [CrossRef] [PubMed]
  25. Pesquita, C.; Faria, D.; Falcao, A.O.; Lord, P.; Couto, F.M. Semantic similarity in biomedical ontologies. PLoS Comput. Biol. 2009, 5, e1000443. [Google Scholar] [CrossRef] [PubMed]
  26. Presneau, N.; Baumhoer, D.; Behjati, S.; Pillay, N.; Tarpey, P.; Campbell, P.J.; Jundt, G.; Hamoudi, R.; Wedge, D.C.; Loo, P.V.; et al. Diagnostic value of H3F3A mutations in giant cell tumour of bone compared to osteoclast-rich mimics. J. Pathol. Clin. Res. 2015, 1, 113–123. [Google Scholar] [CrossRef]
  27. Fadista, J.; Vikman, P.; Laakso, E.O.; Mollet, I.G.; Esguerra, J.L.; Taneera, J.; Storm, P.; Osmark, P.; Ladenvall, C.; Prasad, R.B.; et al. Global genomic and transcriptomic analysis of human pancreatic islets reveals novel genes influencing glucose metabolism. Proc. Natl. Acad. Sci. USA 2014, 111, 13924–13929. [Google Scholar] [CrossRef] [Green Version]
  28. Taneera, J.; Lang, S.; Sharma, A.; Fadista, J.; Zhou, Y.; Ahlqvist, E.; Jonsson, A.; Lyssenko, V.; Vikman, P.; Hansson, O.; et al. A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets. Cell Metab. 2012, 16, 122–134. [Google Scholar] [CrossRef] [Green Version]
  29. Lucki, N.C.; Bandyopadhyay, S.; Wang, E.; Merrill, A.H.; Sewer, M.B. Acid ceramidase (ASAH1) is a global regulator of steroidogenic capacity and adrenocortical gene expression. Mol. Endocrinol. 2012, 26, 228–243. [Google Scholar] [CrossRef] [Green Version]
  30. Chavez, J.A.; Holland, W.L.; Bar, J.; Sandhoff, K.; Summers, S.A. Acid ceramidase overexpression prevents the inhibitory effects of saturated fatty acids on insulin signaling. J. Biol. Chem. 2005, 280, 20148–20153. [Google Scholar] [CrossRef] [Green Version]
  31. Bruce, C.R.; Thrush, A.B.; Mertz, V.A.; Bezaire, V.; Chabowski, A.; Heigenhauser, G.J.; Dyck, D.J. Endurance training in obese humans improves glucose tolerance and mitochondrial fatty acid oxidation and alters muscle lipid content. Am. J. Physiol. Endocrinol. Metab. 2006, 291, E99–E107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Weatherbee, S.D.; Anderson, K.V.; Niswander, L.A. LDL-receptor-related protein 4 is crucial for formation of the neuromuscular junction. Development 2006, 133, 4993–5000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Khan, T.N.; Klar, J.; Ali, Z.; Khan, F.; Baig, S.M.; Dahl, N. Cenani-Lenz syndrome restricted to limb and kidney anomalies associated with a novel LRP4 missense mutation. Eur. J. Med. Genet. 2013, 56, 371–374. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, S.P.; Da, H.; Li, Z.; Kushwaha, P.; Beil, C.; Mei, L.; Xiong, W.C.; Wolfgang, M.J.; Clemens, T.L.; Riddle, R.C. Lrp4 expression by adipocytes and osteoblasts differentially impacts sclerostin’s endocrine effects on body composition and glucose metabolism. J. Biol. Chem. 2019, 294, 6899–6911. [Google Scholar] [CrossRef]
  35. Rankinen, T.; Sarzynski, M.A.; Ghosh, S.; Bouchard, C. Are there genetic paths common to obesity, cardiovascular disease outcomes, and cardiovascular risk factors? Circ. Res. 2015, 116, 909–922. [Google Scholar] [CrossRef] [Green Version]
  36. Rantakari, P.; Lagerbohm, H.; Kaimainen, M.; Suomela, J.P.; Strauss, L.; Sainio, K.; Pakarinen, P.; Poutanen, M. Hydroxysteroid (17β) dehydrogenase 12 is essential for mouse organogenesis and embryonic survival. Endocrinology 2010, 151, 1893–1901. [Google Scholar] [CrossRef] [Green Version]
  37. Locke, A.E.; Kahali, B.; Berndt, S.I.; Justice, A.E.; Pers, T.H.; Day, F.R.; Powell, C.; Vedantam, S.; Buchkovich, M.L.; Yang, J.; et al. Genetic studies of body mass index yield new insights for obesity biology. Nature 2015, 518, 197–206. [Google Scholar] [CrossRef] [Green Version]
  38. Das, S.K.; Sharma, N.K. Expression quantitative trait analyses to identify causal genetic variants for type 2 diabetes susceptibility. World J. Diabetes 2014, 5, 97–114. [Google Scholar] [CrossRef]
  39. Scott, R.A.; Scott, L.J.; Magi, R.; Marullo, L.; Gaulton, K.J.; Kaakinen, M.; Pervjakova, N.; Pers, T.H.; Johnson, A.D.; Eicher, J.D.; et al. An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans. Diabetes 2017, 66, 2888–2902. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart to identify hot spots of chromosomal regions associated with single nucleotide polymorphisms (SNPs) related to the phenotypes identified using the National Human Genome Research Institute-European Bioinformatics Institute Catalog of Published Genome-Wide Association Studies (NHGRI-EBI GWAS) Catalog. T1D: type 1 diabetes; T2D: type 2 diabetes; BMI: body mass index; GWAS: Genome-Wide Association Studies; HDL: high-density lipoprotein; HDLC-TG: cholesterol-triglycerides; LDL: low-density lipoprotein; DM: diabetes mellitus.
Figure 1. Flowchart to identify hot spots of chromosomal regions associated with single nucleotide polymorphisms (SNPs) related to the phenotypes identified using the National Human Genome Research Institute-European Bioinformatics Institute Catalog of Published Genome-Wide Association Studies (NHGRI-EBI GWAS) Catalog. T1D: type 1 diabetes; T2D: type 2 diabetes; BMI: body mass index; GWAS: Genome-Wide Association Studies; HDL: high-density lipoprotein; HDLC-TG: cholesterol-triglycerides; LDL: low-density lipoprotein; DM: diabetes mellitus.
Genes 11 00461 g001
Figure 2. Numbers of chromosomal regions sharing SNPs correlated with the five UAEDIAB-phenotypes (BMI/obesity, dyslipidemia, hypertension, sleep apnea, and T1D/T2D). The regions were extracted from the NHGRI-EBI GWAS Catalog. The figure was generated using the InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams.
Figure 2. Numbers of chromosomal regions sharing SNPs correlated with the five UAEDIAB-phenotypes (BMI/obesity, dyslipidemia, hypertension, sleep apnea, and T1D/T2D). The regions were extracted from the NHGRI-EBI GWAS Catalog. The figure was generated using the InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams.
Genes 11 00461 g002
Figure 3. SNP-based Similarity Analysis of the identified SNPs using the “Socialiser for SNPs” option in the eXploring Genomic Relations (XGR) web tool. Seven SNPs showed similarity to at least 10 other SNPs with a similarity score of more than 0.5 in scale range from 0 to 1.
Figure 3. SNP-based Similarity Analysis of the identified SNPs using the “Socialiser for SNPs” option in the eXploring Genomic Relations (XGR) web tool. Seven SNPs showed similarity to at least 10 other SNPs with a similarity score of more than 0.5 in scale range from 0 to 1.
Genes 11 00461 g003
Figure 4. Gene expression in human pancreatic islets. (A) RNA-seq mean expression of FES, LRP4, ASAH1, and HSD17B12 in non-diabetic human pancreatic islets (n = 63). KCNJ11 gene was used as functional marker for pancreatic islets; (BE) Whisker boxes of differential expression analysis of ASAH1, (B) LRP4, (C) FEX, and (D) HSD17B12; (E) in pancreatic islets of diabetics (n = 12) and non-diabetics (n = 63). HSD17B12 was significantly downregulated in T2D islets (p = 0.03).
Figure 4. Gene expression in human pancreatic islets. (A) RNA-seq mean expression of FES, LRP4, ASAH1, and HSD17B12 in non-diabetic human pancreatic islets (n = 63). KCNJ11 gene was used as functional marker for pancreatic islets; (BE) Whisker boxes of differential expression analysis of ASAH1, (B) LRP4, (C) FEX, and (D) HSD17B12; (E) in pancreatic islets of diabetics (n = 12) and non-diabetics (n = 63). HSD17B12 was significantly downregulated in T2D islets (p = 0.03).
Genes 11 00461 g004
Table 1. List of overlapped chromosomal regions that correlated with UAEDIAB-phenotypes and diabetes.
Table 1. List of overlapped chromosomal regions that correlated with UAEDIAB-phenotypes and diabetes.
Number of Shared PhenotypesShared PhenotypesNumber of Shared RegionsShared Regions/Bands
4 Phenotypes(BMI/obesity) and (dyslipidemia) and (hypertension) and (T1D/T2D)48p22,1q32.3,12q24.13,7p15.2
(BMI/obesity) and (dyslipidemia) and [sleep apnea] and (T1D/T2D)311p11.2,6q21,17q12
(BMI/obesity) and (hypertension) and (sleep apnea) and (T1D/T2D)115q26.1
3 Phenotypes(BMI/obesity) and (dyslipidemia) and (T1D/T2D)3416q12.2,19q13.32,16p11.2,5q13.3,11p15.4,
10p13,16q23.2,18q21.32,18q11.2,6p22.3,
2p23.3,19q13.11,10p15.1,22q12.3,15q15.1,
1p31.3,19p13.2,2q36.3,11p15.1,12q24.12,
12q24.11,11q13.1,2q24.3,1p32.3,1q21.3,
10q25.2,6p21.32,6p21.1,6q23.3,17p13.2,
10q21.3,12q24.31,3p25.2,3q21.1
(dyslipidemia) and (hypertension) and (T1D and T2D)510q23.33,1q41,1q43,6p21.33,8q24.12
(BMI/obesity) and (hypertension) and (T1D/T2D)216p12.3,1p13.2
(BMI/obesity) and (sleep apnea) and (T1D/T2D)31q32.1,1q42.2,11q13.4
2 Phenotypes(dyslipidemia) and (T1D/T2D)176q13,6q27,20q13.12,1q42.13,5q11.2,4p16.3,
6q24.1,12p13.31,9q34.2,8q24.3,17p13.1,
9p24.2,22q12.2,1p34.3,7q32.2,1p22.1,2q33.2
(BMI/obesity) and (T1D/T2D)476q23.1,10q22.3,8q21.13,2p25.3,4p12,
12q13.12,14q31.1,3q27.2,2p16.1,9p21.1,4q28.2,
17q21.32,18q12.3,1p21.3,9q22.31,15q14,
9q22.2,5q33.2,14q11.2,3p14.1,12q12,8q22.3,
21q22.3,2p23.2,17p11.2,2q21.3,9q31.3,
10q26.13,17p13.3,11p15.5,11q13.3,8q24.21,
7q36.3,1p12,3q23,6p21.2,7p12.1,9p24.1,
14q24.1,18p11.21,13q31.1,5q21.1,2p21,12p12.1,
7p14.3,8q22.1,9p21.3
(hypertension) and (T1D/T2D)65q31.1,18p11.31,7q22.1,21q22.11,8p11.21,
9q21.32
(sleep apnea) and (T1D/T2D) 12p24.3
1 Phenotype(T1D/T2D)898q24.11,11q14.3,4q35.1,6q12,17q21.33,
3q26.2,17q11.2,4q22.2,22q13.33,7p21.2,5q22.2,
3q13.31,11p12,15q22.2,2q23.3,11q24.3,
3q26.33,1q32.2,10q26.3,19q13.2,4q32.3,
14q32.2,2q24.2,7q32.1,9q34.3,3p24.3,Xq28,
2q33.1,3q27.3,13q21.31,6q25.1,13q14.13,
13q21.33,13q22.1,16p13.12,20q11.21,1q21.2,
5q14.2,20p12.2,2q14.3,14q23.1,15q24.3,
12q13.2,1p22.3,6q15,16p13.13,18q22.2,2q11.2,
5p13.2,15q25.1,4q27,4p15.2,6q22.32,10q23.31,
17q21.1,17q21.2,20p13,13q22.2,7p15.1,
10q24.2,10q26.11,13q12.12,8q24.22,11p15.4,
11p15.5,3p23,12q21.2,1p22.2,12q21.1,2q12.1,10q22.1,
3q12.3,9p23,9q21.31,12q14.3,
4p16.1,4q31.3,6p24.3,13q12.13,12p11.22,
2p23.1,7p14.1,8q13.2,12p11.21,2p16.2,
10q26.12,16q24.1,3p14.3,20q13.31
Table 2. List of the 46 SNPs located in the identified eight shared chromosomal regions and their reported affected genes. Data extracted from the NHGRI-EBI GWAS Catalog.
Table 2. List of the 46 SNPs located in the identified eight shared chromosomal regions and their reported affected genes. Data extracted from the NHGRI-EBI GWAS Catalog.
SNP IDMapped GeneChro. RegionLocationReported Gene
rs3817334MTCH211p11.211:47629441MTCH2
rs7124681CELF111p11.211:47508395CUGBP1
rs11066280HECTD412q24.1312:112379979HECTD4
rs4430796HNF1B17q1217:37738049HNF1B
rs2176598HSD17B1211p11.211:43842728 HSD17B12
rs17696736NAA2512q24.1312:112049014Not Reported
rs2028299AP3S215q26.115:89831025AP3S2
rs9400239FOXO36q216:108656460FOXO3
rs17126232AC124242.38p228:18120141ASAH1
rs6990042SGCZ8p228:14316465 SGCZ
rs10742752AC103855.311p11.211:45416824SYT13
rs17630235TRAFD1, HECTD412q24.1312:112153882Not reported
rs1439620AC013394.1, LINC0157815q26.115:92886416LOC100507217
rs12150665GGNBP217q1217:36558947GGNBP2
rs3800229FOXO36q21 6:108675760FOXO3
rs35424364CCDC162P6q216:109322403C6ORF183, CCDC162P
rs1495741PSD3, NAT28p228:18415371 NAT2
rs10838738MTCH211p11.211:47641497MTCH2
rs326214MADD11p11.211:47276809LRP4
rs74472562TSPAN1811p11.211:44741205RP11-45A12.2, TSPAN18
rs1061810HSD17B12, AC087521.2, AC087521.411p11.211:43856384HSD17B12
rs936674AC091078.115q26.115:93360368RP11-266O8.1
rs148024591AC091078.115q26.115:93371222RP11-266O8.1
rs2521501FES15q26.115:90894158FURIN, FES
rs8042680PRC1, PRC1-AS115q26.115:90978107PRC1
rs12899811VPS33B15q26.115:91000846PRC1
rs79548680RCCD115q26.115:90962549PRC1
rs1877031STARD317q1217:39657827STARD3
rs4796285AC243830.1, LHX1-DT17q1217:36824731Not reported
rs10908278HNF1B17q1217:37739961HNF1B, TCF2
rs1704198AC096639.11q32.31:213737151PROX1
rs340839PROX11q32.31:213988477PROX1
rs7526425AC105275.1, SLC30A1, RD31q32.31:211527316SLC30A1
rs2075423PROX1-AS11q32.31:213981376PROX1
rs884366CCDC162P6q216:109252892LOC100996634
rs149358103RPS27AP11, LINC025416q216:11358684SOCS5P5, MARCKS
rs10261878AC010719.1, AC018706.17p15.27:25910925NFE2L3, MIR148A
rs4719841NFE2L3, MIR148A7p15.27:25957916MIR148A
rs4722551NFE2L3, MIR148A7p15.27:25952206MIR148A
rs6969780HOXA3, HOXA-AS27p15.27:27119517HOXA3
rs10279895HNRNPA1P73, RPL35P47p15.27:27288591EVX1, HOXA
rs7804356SKAP27p15.27:26852046Not reported
rs4921914PSD3, NAT28p228:18414928NAT2
rs1961456NAT28p228:18398199NAT2
rs115706913SGCZ8p228:14224308Not Reported
rs2946504TRMT9B8p228:12954071KIAA1456
Table 3. Enrichment analysis of the 46 SNPs found in eight shared chromosomal regions in at least four of the UAEDIAB-phenotypes.
Table 3. Enrichment analysis of the 46 SNPs found in eight shared chromosomal regions in at least four of the UAEDIAB-phenotypes.
Term NameNumber of SNPs OverlappedSNP IDsZ-Score p-ValueFDR
1Metabolic Disease15rs1061810, rs10908278, rs11066280, rs12899811, rs149358103, rs1704198, rs17126232, rs17696736, rs2028299, rs2075423, rs2946504, rs4430796, rs7804356, rs79548680, rs804268013.28.20 × 10−152.60 × 10−13
2Diabetes Mellitus12rs1061810, rs10908278, rs12899811,rs149358103, rs17696736, rs2028299, rs2075423, rs2946504, rs4430796, rs7804356, rs79548680, rs804268012.21.50 × 10−122.30 × 10−11
3Body Mass Index13rs10261878, rs10742752, rs10838738, rs12150665, rs1439620, rs17630235, rs2176598, rs3800229, rs3817334, rs6990042, rs7124681, rs936674, rs940023910.96.40 × 10−124.80 × 10−11
4Sleep Apnea4rs148024591, rs35424364, rs4796285, rs74472562 29.77.80 × 10−124.80 × 10−11
5Sleep Apnea4rs148024591, rs35424364, rs4796285, rs7447256229.77.80 × 10−124.80 × 10−11
6Type II Diabetes Mellitus10rs1061810, rs10908278, rs12899811, rs2946504, rs149358103, rs8042680, rs2028299, rs2075423, rs4430796, rs7954868011.92.30 × 10−111.20 × 10−10
7Sleep Apnea Measurement4rs148024591, rs35424364, rs4796285, rs7447256220.53.90 × 10−101.70 × 10−9
8Sleep Disorder4rs148024591, rs35424364, rs4796285, rs7447256215.56.50 × 10−92.50 × 10−8
9Hypertension5rs10279895, rs11066280, rs115706913, rs2521501, rs696978012.11.40 × 10−84.80 × 10−8
10Triglyceride Measurement6rs11066280, rs1495741, rs340839, rs4719841, rs4722551, rs49219147.95.70 × 10−71.8 × 10−6
11Lipid Measurement9rs11066280, rs1495741, rs1877031, rs326214, rs340839, rs4719841, rs4722551, rs4921914, rs8843666.211.9 × 10−65.4 × 10−6
12Lipoprotein Measurement7rs11066280, rs1495741, rs1877031, rs1961456, rs326214, rs4722551, rs8843665.451.9 × 10−50.00005
13Diastolic Blood Pressure5rs10279895,rs11066280, rs17696736, rs2521501, rs69697805.54.1 × 10−59.7 × 10−5
14Mean Arterial Pressure3rs17696736, rs2521501, rs69697806.25.3 × 10−50.00012
15Physical Activity Measurement3rs3800229, rs3817334, rs71246815.310.000150.00028
16High-Density Lipoprotein Cholesterol Measurement4rs11066280, rs1877031, rs326214, rs8843664.920.000150.00028
17Drinking Behavior3rs11066280, rs17696736, rs25215015.270.000160.00028
18Obesity2rs1704198, rs171262325.40.000250.00043
19Systolic Blood Pressure4rs11066280, rs17696736, rs2521501, rs69697804.210.000460.00074
20Parental Longevity2rs17630235, rs176967364.370.000730.0011
21Total Cholesterol Measurement3rs1495741, rs1961456, rs47225514.010.000810.0012
22Type I Diabetes Mellitus2rs17696736, rs78043564.260.000830.0012
23Blood Pressure5rs10279895, rs11066280, rs17696736, rs2521501, rs69697803.680.000890.0012
24Behavior5rs11066280, rs17696736, rs2521501, rs3817334, rs94002393.610.0010.0013
25Alcohol Drinking2rs17696736, rs25215013.540.00190.0023
26Longevity2rs17630235, rs176967363.530.0020.0023
27Blood Metabolite Measurement2rs1495741, rs49219143.530.0020.0023
28Vital Signs5rs10279895, rs11066280, rs17696736, rs2521501, rs69697803.10.00260.0029
29Smoking Behavior2rs3817334, rs94002391.710.0260.028
Table 4. List of identified SNPs with related genes that showed similarity to each other using the Socialiser tool.
Table 4. List of identified SNPs with related genes that showed similarity to each other using the Socialiser tool.
SNPGene Related to the SNPNumber of SNPs with Similarity Score > 0.5 in a Scale of 0–1
rs11066280ALDH223
rs17696736C12orf3015
rs17126232ASAH113
rs1704198PROX113
rs9400239FOXO310
rs3800229FOXO39
rs1961456NAT29
rs1495741NAT29
rs4722551MIR148A8
rs884366LOC1009966347
rs326214LRP47
rs1877031STARD37
rs4719841MIR148A6
rs340839PROX16
rs2028299AP3S5
rs10908278HNF1B, TCF25
rs1061810HSD17B125
rs2946504KIAA14565
rs4921914NAT25
rs12899811PRC15
rs79548680PRC15
rs8042680PRC15
rs2075423PROX15
rs2075423PROX15
rs149358103SOCS5P5, MARCKS5
Table 5. List of SNPs identified by next-generation sequencing (NGS) from two diabetic patients in genes located in the shared chromosomal region.
Table 5. List of SNPs identified by next-generation sequencing (NGS) from two diabetic patients in genes located in the shared chromosomal region.
Identified GenesNumber of SNPs per Sample
Patient 1Patient 2
ASAH12016
LRP487
FES96
HSD17B1274
HNF1B44
HECTD426
SH2B344
NAT242
SGCZ31
FURIN31
GGNBP251
TSPAN1822
ALDH211
STARD311
MTCH212
NDUFS311
PRC100
PTPN1110
CELF130
NAA2520

Share and Cite

MDPI and ACS Style

Hachim, M.Y.; Aljaibeji, H.; Hamoudi, R.A.; Hachim, I.Y.; Elemam, N.M.; Mohammed, A.K.; Salehi, A.; Taneera, J.; Sulaiman, N. An Integrative Phenotype–Genotype Approach Using Phenotypic Characteristics from the UAE National Diabetes Study Identifies HSD17B12 as a Candidate Gene for Obesity and Type 2 Diabetes. Genes 2020, 11, 461. https://doi.org/10.3390/genes11040461

AMA Style

Hachim MY, Aljaibeji H, Hamoudi RA, Hachim IY, Elemam NM, Mohammed AK, Salehi A, Taneera J, Sulaiman N. An Integrative Phenotype–Genotype Approach Using Phenotypic Characteristics from the UAE National Diabetes Study Identifies HSD17B12 as a Candidate Gene for Obesity and Type 2 Diabetes. Genes. 2020; 11(4):461. https://doi.org/10.3390/genes11040461

Chicago/Turabian Style

Hachim, Mahmood Y., Hayat Aljaibeji, Rifat A. Hamoudi, Ibrahim Y. Hachim, Noha M. Elemam, Abdul Khader Mohammed, Albert Salehi, Jalal Taneera, and Nabil Sulaiman. 2020. "An Integrative Phenotype–Genotype Approach Using Phenotypic Characteristics from the UAE National Diabetes Study Identifies HSD17B12 as a Candidate Gene for Obesity and Type 2 Diabetes" Genes 11, no. 4: 461. https://doi.org/10.3390/genes11040461

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop