Next Article in Journal
Calprotectin, a Promising Serological Biomarker for the Early Diagnosis of Superinfections with Multidrug-Resistant Bacteria in Patients with COVID-19
Previous Article in Journal
Strategies to Maintain Redox Homeostasis in Yeast Cells with Impaired Fermentation-Dependent NADPH Generation
Previous Article in Special Issue
The Importance of Molecular Genetic Testing for Precision Diagnostics, Management, and Genetic Counseling in MODY Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response

by
Igor S. Damarov
1,
Elena E. Korbolina
1,*,
Elena Y. Rykova
1,2 and
Tatiana I. Merkulova
1
1
Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
2
Department of Engineering Problems of Ecology, Novosibirsk State Technical University, 630087 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(17), 9297; https://doi.org/10.3390/ijms25179297
Submission received: 11 July 2024 / Revised: 14 August 2024 / Accepted: 26 August 2024 / Published: 27 August 2024
(This article belongs to the Special Issue Advances in Molecular Research of Diabetes Mellitus)

Abstract

:
The goal of our study was to identify and assess the functionally significant SNPs with potentially important roles in the development of type 2 diabetes mellitus (T2DM) and/or their effect on individual response to antihyperglycemic medication with metformin. We applied a bioinformatics approach to identify the regulatory SNPs (rSNPs) associated with allele-asymmetric binding and expression events in our paired ChIP-seq and RNA-seq data for peripheral blood mononuclear cells (PBMCs) of nine healthy individuals. The rSNP outcomes were analyzed using public data from the GWAS (Genome-Wide Association Studies) and Genotype-Tissue Expression (GTEx). The differentially expressed genes (DEGs) between healthy and T2DM individuals (GSE221521), including metformin responders and non-responders (GSE153315), were searched for in GEO RNA-seq data. The DEGs harboring rSNPs were analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). We identified 14,796 rSNPs in the promoters of 5132 genes of human PBMCs. We found 4280 rSNPs to associate with both phenotypic traits (GWAS) and expression quantitative trait loci (eQTLs) from GTEx. Between T2DM patients and controls, 3810 rSNPs were detected in the promoters of 1284 DEGs. Based on the protein-protein interaction (PPI) network, we identified 31 upregulated hub genes, including the genes involved in inflammation, obesity, and insulin resistance. The top-ranked 10 enriched KEGG pathways for these hubs included insulin, AMPK, and FoxO signaling pathways. Between metformin responders and non-responders, 367 rSNPs were found in the promoters of 131 DEGs. Genes encoding transcription factors and transcription regulators were the most widely represented group and many were shown to be involved in the T2DM pathogenesis. We have formed a list of human rSNPs that add functional interpretation to the T2DM-association signals identified in GWAS. The results suggest candidate causal regulatory variants for T2DM, with strong enrichment in the pathways related to glucose metabolism, inflammation, and the effects of metformin.

1. Introduction

Type 2 diabetes mellitus (T2DM) is the most common form of diabetes and one of the most prevalent metabolic disorders, which affects hundreds of millions of individuals worldwide [1]. The etiology of this multifaceted disease involves an interplay of multiple genetic, epigenetic, developmental, and environmental factors, including age, metabolic syndrome, and lifestyle. T2DM also overlaps in epidemiology, pathogenesis, and genetics as revealed by genome-wide association studies (GWAS) [2] with obesity; notably, there is also a strong immune connection between the two diseases [3,4].
The characteristics of T2DM are the alterations in circulating glucose level regulation, chronic hyperglycemia, and insulin resistance [5,6]. Emerging data suggest that the glucotoxicity caused by chronic hyperglycemia injures many cell types, such as pancreatic cells, leading to loss in functional β-cells [7]; hepatic cells, inducing endoplasmic reticulum stress; hepatocyte cell death; nonalcoholic fatty liver disease [8]; innate immune cells, inducing an inflammatory response [9]; and activation of monocytic lineage cells [10]. In turn, the chronic inflammatory state associated with the onset of T2DM leads to the development of long-term consequences: macrovascular complications, including a range of cardiovascular diseases [11,12] and microvascular complications, including retinopathy, nephropathy, and neuropathy [13]. Furthermore, it may contribute to the chronic inflammation of the central nervous system and neurodegeneration [14], as well as contribute to the association of T2DM with other conditions promoted by inflammation, such as rheumatoid arthritis [15].
The candidate genes [16,17,18] and GWAS [19,20] have linked single nucleotide polymorphisms (SNPs) at more than 250 loci in the human genome to T2DM risk, laying the foundations for functional investigations. Among them, the association of TCF7L2 (transcription factor 7-like 2) with T2DM seems to be the most promising since many TCF7L2 variants that increase the disease risk have been replicated in numerous studies and populations with diverse genetic origins [21]. The gene encodes a downstream effector of the canonical Wnt/β-catenin signaling pathway, which has been associated with several fundamental processes, including adipogenesis [22]. Most of the TCF7L2 variants revealed by GWAS are located in noncoding genome regions suggesting that they exert their effects by modulating expression. The high-performance sequencing of formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) in human islets demonstrated that TCF7L2 at-risk variants were mapped in open chromatin sites. Among them, the heterozygotes of rs7903146 showed that enhancer activity depended on the allele [23]. A therapeutic potential was reported for some TCF7L2 variants, for example, intronic rs290487 influenced the efficacy of repaglinide, an oral hypoglycemic agent, in Chinese T2DM patients [24].
Yet deciphering the functionality of the known risk variants on a genome-wide scale is still challenging. That is why different omics approaches are ever more frequently used to study the functional significance of noncoding variants. In particular, Yan et al. [25] measured the transcription factor–DNA interactions using an ultra-high-throughput multiplex protein–DNA binding assay, termed SNP-SELEX (SNP evaluation by Systematic Evolution of Ligands by Exponential enrichment). The authors examined the in vitro binding of 270 human transcription factors to common sequence variants. The input DNA library contained more than 380,000 oligonucleotides originally designed to represent over 95,000 SNPs. They covered 110 SNPs linked to T2DM susceptibility by GWAS at the start of the project, 6724 SNPs in high LD with them (r2 ≥ 0.8), and 89,162 common SNPs in annotated candidate cis-regulatory sequences located within 500 kb of GWAS-derived SNPs. As a result, approximately 11% of the input SNPs (11,079 SNPs) showed significantly differential binding to at least one transcription factor to predict the potential impact on the involved molecular pathways.
In this study, we used a multi-omics approach to search for the functionally significant SNPs associated with T2DM. First, we obtained the ChIP-seq data for histone modifications, histone H3 lysine K4 trimethylation (H3K4me3), histone H3 lysine 27 acetylation (H3K27ac), and the transcriptome sequencing data (RNA-seq) for peripheral blood mononuclear cells (PBMCs) from nine healthy individuals. After the primary data processing and filtering, we searched for allele-asymmetric binding and expression events in paired ChIP-seq and RNA-seq data for each individual. Thus, we identified 14,796 rSNPs that influenced the expression of 5132 genes, of which about 40% coincided with the corresponding GWAS data. The top GWAS phenotypic associations for these rSNPs included T2DM and a group of closely related characteristics (waist–hip index, waist-to-hip ratio adjusted for BMI, and a body shape index). Then, we used publicly available RNA-seq data from both healthy individuals and the individuals with T2DM (GSE221521), including metformin responders and non-responders (GSE153315), to evaluate the rSNPs potentially involved in T2DM pathogenic mechanisms and drug response.

2. Results

2.1. Algorithm for Searching the rSNPs Associated with Allele-Specific Events in ChIP-seq and RNA-seq Data and Construction of rSNPs Panel

The rSNP search strategy comprised six main stages (Figure 1). At the first stage, PBMCs were isolated from the peripheral blood of nine healthy individuals. At the next stage, the active regulatory genome regions marked with histone modifications H3K4me3 and H3K27ac (ChIP-seq) were sequenced. Concurrently, the transcriptome of these PBMC samples was sequenced (RNA-seq). Allelic asymmetry was separately computed for each heterozygous SNP using the ChIP-seq and RNA-seq data. Next, allele-specific binding (ASB) SNPs located within ±1000 bp from the transcription start sites (TSSs) of the known genes were extracted from ChIP-seq data. At the final stage of the search for rSNPs, the promoter ASB SNPs were intersected with the corresponding genes with determined allele-specific expression (ASE) events. The resulting rSNPs panel was characterized using the GWAS [26,27], GTEx [28,29], and annotation and enrichment analysis of allele-specific transcription factor binding at SNP (ANANASTRA) [30,31] data. Then, the representation of rSNPs in the promoters of the genes with differential expression in the blood cells of T2DM patients versus healthy subjects (GSE221521) and in response to metformin (GSE153315) was assessed followed by the characterization of these genes utilizing the Search Tool for the Retrieval of Interacting Genes (STRING) [32,33], KEGG [34,35], GO [36,37], Reactome [38,39], and literature data.
The total volume of our NGS data is 278.3 GB (ChIP-seq, 142 GB and RNA-seq, 136.3 GB); on average, the information volume per individual is 15.35 GB (15.7 GB, ChIP-seq and 15.1 GB, RNA-seq), which is 177.6 million ChIP-seq paired-end reads and 176 million RNA-seq paired-end reads per individual. The determination of allele-asymmetric events allowed us to find 14,796 rSNPs associated with the expression of 5132 genes (Supplementary Table S1).

2.2. Characterizing Constructed rSNP Panel with GWAS, GTEx, and ANANASTRA

The interception of rSNPs with the genome regions within ±1000 bp from SNPs in the GWAS catalog (as of March 2023) demonstrated the association of 5688 rSNPs (38.4% of all rSNPs of the panel) with different phenotypic traits. Note that 1107 rSNPs (7.5% of all rSNPs of the panel) were directly represented in the GWAS catalog (Supplementary Table S2). The χ2 estimation demonstrated the rSNP panel enrichment in GWAS variants (p-value < 2.2 × 10–16; OR = 1.7; and 95% CI, 1.6–1.8) compared with all heterozygous non-regulatory SNPs.
The resulting rSNP set is enriched in GWAS phenotype associations (Supplementary Table S3) (Fisher’s test with Benjamini–Hochberg. p.adj < 0.1). Preliminarily, we excluded the traits associated with less than 100 GWAS variants from the analysis. In the resulting group of traits ranked according to the number of GWAS-derived rSNPs, most of the traits among the top 20 ones were associated with quantitative blood cell characteristics (Table 1). In addition, this list contains T2DM and the morphometric characteristics directly associated with this disease, namely, waist–hip index, waist-to-hip ratio adjusted for BMI, and a body shape index.
We searched the GTEx Consortium data as of March 2023 for the rSNPs associated with the altered gene expression in different tissues (search for eQTLs) and found that 9871 rSNPs (66.7% of the overall rSNP panel) are associated with eQTLs; notably, 2474 rSNPs were directly represented in the GTEx catalog and 7397 ones situated within ±1000 bp of them (Supplementary Table S4). Of all tissues represented in GTEx, the largest number of rSNPs, 2464, was expectedly detected in the whole blood (the data on individual blood cell fractions are absent in the catalog). The enrichment analysis demonstrated that the rSNP panel was enriched in the GTEx eQTL variants (p-value < 2.2 × 10–16; OR = 2.49; 95% CI, 2.38–2.6) as compared with all heterozygous non-regulatory SNPs identified in our ChIP-seq data.
In total, 4280 rSNPs (28.9% of all detected rSNPs) in 1628 genes were simultaneously associated with an eQTL effect (according to GTEx data) and a certain phenotypic trait (according to GWAS data) (Figure 2). Published data suggest the involvement of many of these genes in inflammatory processes (ADAM17 [40,41,42], AHNAK [43,44], AIF1 [45,46,47], CCDC92 [48], CTBP1 [49], HCP5 [50], MAPKBP1 [51], and MAST3 [52]); obesity (TP53INP1 [53], TKT [54], ADAM17 [40,41], and AHNAK [55]), insulin resistance (TP53INP1 [53], TKT [54], and CCDC92 [48]); insulin secretion (GIPR [56], KCNJ15 [57], and AP3S2 [58], ARAP1 [59]); proliferation of pancreatic β-cells (CDKN1B [60] and CCND2 [61]); mitophagy (TP53INP1 [53]); glycolysis (TKT [54]); and carbohydrate transport (GBA2 [62]). Approximately 14% of these rSNPs (598) were associated with T2DM or tightly correlated with its traits according to GWAS data (Supplementary Table S5).
The intersection of our list of all 14,796 rSNPs with the variants detected in the ANANASTRA allele-specific binding loci of different transcription factors (ChIP-seq data) showed that 30.8% (4560, FDR < 0.05) of our rSNPs were localized to the allele-asymmetric binding sites of different transcription factors (Figure 2). The share of such rSNPs increases to 33.4% (1433, FDR < 0.05) for the list of 4280 rSNPs associated with eQTLs and different GWAS phenotypic traits. The largest number of the 4280 rSNPs was situated in the allele-specific binding sites of ANDR (N = 217), CTCF (N = 215), STAT1 (N = 167), BRD4 (N = 126), ESR1 (N = 121), YY1 (N = 97), GCR (N = 89), SPT5H (N = 76), STAG1 (N = 74), and ZFX (N = 63).

2.3. Analysis of Hub Genes, Modules, and Pathways in Associative Gene Networks of rSNP-Governed Differentially Expressed Genes Related to T2DM

2.3.1. Search for Differentially Expressed Genes Related to T2DM and Harboring rSNPs within Promotors

As described above, the formed rSNP panel emerged to be enriched in the variants related to T2DM and the associated traits. To further clarify the putative role of the discovered rSNPs in the mechanisms underlying T2DM, we analyzed the RNA-seq (GSE221521) data deposited with the GEO datasets on the PBMCs of healthy subjects (N = 50), patients diagnosed with this disease (N = 74), and T2DM subjects with diabetic retinopathy (N = 69), representing the largest sample of the open access ones [63]. Note that these data were already searched for DEGs [63]; however, the authors did not give their full list. Correspondingly, we recomputed differential expression setting other significance thresholds (|log2FC| > 0.2 and p.adj. < 0.05).
We did not discover any DEGs when comparing the cohorts of healthy subjects and the T2DM patients without retinopathy. However, the comparison of the healthy cohort and the T2DM individuals with retinopathy allowed us to detect 4612 DEGs and 1284 of them contained 3810 rSNPs in their promoters. Moreover, 2481 rSNPs resided in the promoters of 772 upregulated DEGs (log2FC > 0) and 1329 rSNPs in the promoters of 512 downregulated DEGs (log2FC < 0) (Figure 3). Notably, among these rSNPs, rs893617 (AP3S2), rs1552224 (ARAP1), rs3744347 (CBX1), and rs2066827 (CDKN1B) are associated with T2DM according to GWAS data. The genes in which these rSNPs are located play an important role in the pathogenesis of T2DM [58,59,60]. As we mentioned earlier, AP3S2 and ARAP1 are involved in insulin secretion [58,59]. Increased expression of CDKN1B is associated with structural and functional changes in the kidneys in diabetic nephropathy [64], decreased pancreatic β-cell proliferation and serum insulin levels [60], as well as decreased macrophage proliferation and inflammation in atherosclerosis [65].

2.3.2. Identification of Hub Genes and Analysis of Key Modules Using STRING-Based Protein Interactions, KEGG, and GO Enrichment

The detected DEGs carrying rSNPs in their promoters were searched for hub genes, typically defined as genes involved in the regulation of various biological processes via the interaction with their numerous target genes or proteins. We assumed that the rSNPs in the promoters of hub genes influenced their expression and thereby significantly interfered with the regulation of key processes associated with the development of diabetes and its complications.
Using STRING, we constructed a PPI network for the protein products of upregulated DEGs, comprising 735 nodes and 2243 edges. The computations of topological characteristics for each node in this network allowed us to identify the protein products of 31 hub genes (Supplementary Table S6). Hub genes were mainly represented by the genes involved in transcription regulation, namely, coding for transcription factors (TP53, SREBF1, and RXRA), a corepressor (NCOR2), a chromatin-binding protein (BRD4), histone deacetylase (HDAC4), and RNA polymerase II subunit RPB1 (POLR2A). The genes coding for serine/threonine kinases and membrane receptors formed a separate group. Many genes on this list are involved in different processes associated with T2DM development [66,67,68], such as obesity (SREBF1 and H6PD), inflammation (TGFB1), and insulin resistance (RPTOR and AKT1). Notably, according to the analysis of the enrichment of hub genes in KEGG pathways, the insulin signaling pathway (N of genes = 7) displayed the maximum enrichment (p.adj. = 5.6 × 10–7). The top ten pathways with the highest enrichment values include AMPK (p.adj. = 3.93 × 10–5) and FoxO signaling pathways (p.adj. = 3.45 × 10–6), which are responsible for the regulation of glucose metabolism [69,70] (Figure 4).
Three significant modules were discovered in the PPI network constructed for upregulated genes. The most important of them (MCODE score = 5.514) comprised 36 nodes and 102 edges (Figure 5A) and included the protein products of nine hub genes: TGFB1, RPTOR, HDAC4, BRD4, H6PD, NCOR2, RXRA, HTT, and ITGAM. The analysis of enrichment in KEGG demonstrated that the genes coding for the proteins of this PPI module were involved in the PI3K-Akt signaling pathway (p.adj. = 1.98 × 10−6) and insulin signaling pathway (p.adj. = 1.75 × 10−6), directly associated with T2DM development [71]. Of the GO terms associated with T2DM pathogenesis, we identified those related to the immune response, namely, myeloid leukocyte activation (p.adj. = 0.003), inflammatory response (p.adj = 0.02), immune response (p.adj. = 0.04), mononuclear cell differentiation (p.adj. = 0.005), and leukocyte activation (p.adj. = 0.002) (Supplementary Table S7).
The PPI module second in significance comprises 33 nodes and 73 edges (Figure 5B) and includes the protein products of three hub genes: TP53, POLR2A, and MDM2. Analysis of the enrichment in GO terms showed that the genes in this module, similar to those in the first one, are involved in the innate immune response (p.adj. = 0.02). In addition, the response to peptides is among the significant terms (p.adj. = 0.04). A number of genes associated with this term (TP53, MDM2, IRF5, and TYK2) are directly involved in the insulin response [72,73,74,75]. The KEGG enrichment analysis failed to detect any statistically significant pathways associated with T2DM; however, such pathways were discovered using the Reactome data [38,39]: interferon signaling (p.adj. = 3.2 × 10−6), cytokine signaling in immune system (p.adj. = 1.8 × 10−4), and transcriptional regulation by TP53 (p.adj. = 0.005) (Supplementary Table S7).
The third module of 24 nodes and 50 edges (Figure 5C) comprises the protein products of five hub genes: FASN, MAPK3, SREBF1, INSR, and NOTCH1. KEGG pathway enrichment analysis discovered insulin resistance (p.adj. = 0.003), mTOR signaling (p.adj. = 0.003), T2DM (p.adj. = 0.03), insulin signaling (p.adj. = 0.003), and AMPK signaling (p.adj. = 0.01) pathways (Supplementary Table S7).
The constructed PPI network for the downregulated DEGs with the promoters harboring detected rSNPs comprises 485 nodes and 3961 edges. We distinguish 7 modules and 35 hub genes in this network, their protein products mainly belong to ribosomal subunits and chaperons (Supplementary Table S6). The analysis of enrichment in KEGG/GO/Reactome DEG protein products of each module detected the common terms and pathways related to ribosome biogenesis, transcription processes, rRNA processing, etc. (Supplementary Table S7), poorly interpretable with respect to T2DM pathogenesis. Presumably, the inhibition of protein synthesis is associated with the intake of antihyperglycemic drugs [76,77] and is a result of the therapy. The representation analysis demonstrates that the DEGs of this group are involved in the pathways involved in glucose metabolism, namely, insulin (EIF4E2) and PI3K/Akt signaling (EIF4E2, GNB5, HSP90AB1, PPP2CA, CDKN1B, ATF4, BCL2, and YWHAB) pathways. Moreover, BCL2 and HSP90AB1 were identified as hub genes.
Taking into account the upregulated genes, 15 DEGs were represented in the insulin signaling pathway and 22 DEGs, in the PI3K/Akt signaling pathway. Figure 6A,B show the scheme of the PI3K/Akt signaling pathway with highlighted up- and downregulated DEGs harboring rSNPs in their promoters.

2.3.3. Selecting Important Regulators from PPI Network Using ROC Analysis

Using an additional ROC analysis stage, we distinguished the hub genes the expression level of which allowed the healthy individuals to be distinguished from T2DM subjects. Of 66 hub genes, 36 displayed AUC > 0.7, making them valuable T2DM predictors for the studied set of transcriptome data. The promoters of these hub genes contained 117 rSNPs which are putatively associated with T2DM development and are of great interest for further studies (Table 2).

2.4. Search for the rSNPs Potentially Associated with Individual Response to the Antidiabetic Drug Metformin

Using our rSNP panel, we searched for the rSNPs potentially associated with the response to metformin by analyzing the RNA-seq data (GSE153315, GEO DataSets) on the blood cells of T2DM non-responders (N = 10) and T2DM responders (N = 10) after a 3-month medication [78]. In total, 406 DEGs were found for these cohorts. As it emerged, the promoter regions of 131 DEGs harbored 367 rSNPs of our panel (Figure 7).
The group of genes coding for transcription factors emerged to be the largest among those that carried predicted rSNPs in their promoters, namely, (1) FOXP1 (rs9827299 and rs927809457), coding for the well-known glucose homeostasis regulator [79], which is also a known suppressor of macrophage differentiation into an anti-inflammatory phenotype [80]; (2) POU2F2 (rs3826705), the product of which is an important player in glucose homeostasis being an activator of the AKT/mTOR signaling pathway [81] and is involved in the cellular immune response by regulating B-cell proliferation and differentiation genes [82]; (3) YY-1 (rs113799953 and rs760756349), which is expressed in all cell types being involved in a wide range of regulatory processes, including the regulation of lipid and glucose metabolisms and immune response [83]; (4) KLF6 (rs3812715, rs17135808, rs3829201, rs10795076, and rs38127140), the product of which is associated with the differentiation of monocytes to macrophages [84]; (5) NFIX (rs11555274, rs1354007230, rs1293177918, and rs1235554616), the product of which is a putative myelopoiesis regulator [85]; and (6) GMEB1 (rs1013370834, rs1022902336, and rs895890911), the product of which is necessary for a normal function of islet endothelial cells [86]. We also detected a considerable number of rSNPs in the promoter regions of the genes coding for the transcription factors with yet vague functions: TSHZ2 (rs71354397), ZNF121 (rs7253981), ZNF44 (rs386901 and rs422563), ZNF440 (rs1056484229), ZNF611 (rs1350535255, rs8109324, rs11879101, rs707303, rs10419223, rs1044394923, and rs55769230), ZNF714 (rs7258409, rs182211521, and rs61733856), ZNF793 (rs2291003), and ZNF813 (rs1467168591, rs1157700262, rs2015145, rs4599028, rs112581977, rs1968829, and rs181026484).
The genes of other transcription regulators are functionally related to this group. They include (1) TAF3 (rs78092115 and rs2778473), coding for a basal transcription factor involved in the formation of RNA polymerase II initiation complex and histone modifications; the association of TAF3 locus and a decreased Mg2+ level in T2DM were recently demonstrated [87]; (2) NCOR1 (rs8076864, rs1331834429, and rs178797), coding for a corepressor protein involved in insulin resistance [88]; (3) SMARCA4 (rs2288844, rs17766161, rs1057226903, and rs4804550), coding for BRG1, a catalytic subunit of SWI/SNF chromatin remodeling complex involved in the regulation of inflammatory processes in many tissues; in particular, a deficiency in BRG1 decreases the risk of diabetic atherosclerosis in mice by suppressing an inflammatory response in vessels [89]; and (4) PCBD2 (rs319597), the product of which is involved in the transcription upregulation by interacting with the transcription factor HNF1β, mutations in which are among the key genetic causes maturity-onset diabetes of the young [90].
Other distinguished groups of genes are involved in the transport function (NPHP3, SEC62, SLC30A7, SLC35E4, SLC37A2, SPNS1, and USE1; in total, they contain 13 rSNPs in their promoters); are membrane proteins (SMIM14, SPAG1, TGOLN2, TMEM120B, TMEM245, and VAMP8; in total, 14 rSNPs); are involved in ubiquitination (NOSIP, RNF145, TRIM56, TRIM58, and UBE2G2; in total, 12 rSNPs) and inflammation (TNFSF14, PILRB, AIF1, PEBP1, CDE3; in total, 11 rSNPs); and are long noncoding RNAs and RNA-binding proteins (NSRP1, RBM26, SNHG3, SNHG5, TPT1-AS1, and DDX52; in total, 12 rSNPs) and ribosomal proteins or the enzymes involved in ribosome assembly (NSUN4, RPL14, RPL27A, RPL37A, RPLP2, RPS23, and RPS2P32; in total, 24 rSNPs). Noteworthy, an analysis of the relevant literature demonstrated that several of these genes were involved in the processes underlying T2DM development/prevention and its complications [91,92,93,94,95,96,97,98,99,100,101,102].

3. Discussion

One of the main challenges in contemporary human genetics is the discovery of the genome variations influencing different biomedical traits and the elucidation of the molecular mechanisms allowing these variations to lead to phenotypic differences and complex diseases. So far, over half a million variants (mainly, SNPs) associated with various human traits and diseases have been discovered with the help of GWAS [27]. The overwhelming majority of these variants are situated in the noncoding part of the genome, which suggests their regulatory function [103,104]. Since the GWAS technology cannot give information about the functionality of the found variants, the research into a functional interpretation of the GWAS data at the level of individual SNPs and on a mass scale is increasing [103,105], as well as the design of the GWAS-independent large-scale functional approaches based on omics data analysis [29,103,106].
In this work, we used a multi-omics approach to search for the functional variants able to influence gene expression (rSNPs) in PBMCs, which are the most available human tissue involved in immunity and diabetes [107,108,109,110]. We detected allele-asymmetric events in our ChIP-seq data for H3K4me3 and H3K27ac histone marks and the transcriptome data (RNA-seq) for nine healthy individuals and identified 14,796 rSNPs in the promoter regions of 5132 genes. Note that obtaining similar data for other specific affected organs, such as pancreatic beta cells, liver, muscle, etc., would significantly expand the rSNP panel.
We discovered a considerable enrichment of our rSNP panel in GWAS variants as compared with a set of heterozygous non-regulatory SNPs (p-value = 2.8884 × 10–8), which confirmed a significant role of the discovered rSNPs in the formation of the traits. The intersection of rSNPs with the DNA regions within ±1000 bp of GWAS SNPs demonstrated the association of 5688 rSNPs (38.4% of all detected ones) with different phenotypic traits. Of them, 1107 rSNPs (7.5%) were directly represented in the GWAS catalog. This is a good matching of results for principally different approaches. On the one hand, this demonstrates that a considerable part of the detected rSNPs is actually associated with certain phenotypic traits and, on the other hand, uncovers the molecular sense of a large GWAS-identified SNP group.
We obtained a considerably better fit of results when comparing our 14,796 rSNPs to the GTEx project data on eQTL analysis, which is also a functional approach that links different SNP alleles to gene expression level [103]. It was found that 9871 rSNPs (66.7% of the panel) are associated with eQTLs; in addition, 2474 rSNPs were directly contained in the GTEx catalog and 7397 rSNPs fell within the ±1000 bp region from them. On the contrary, the result of the comparison of the functional approach utilizing the search for allele-asymmetric events in the binding sites of different transcription factors [31] was rather modest with only 39.4% (5844) of our rSNPs localized to the allele-asymmetric binding sites of different transcription factors. In our opinion, the better fit of our results to GTEx data compared with ANANASTRA is explainable with that our rSNP search pipeline is based on the analysis of both allele-specific binding (ASB) and ASE events of the corresponding genes, while ANANASTRA utilizes only ASB events.
The computation of the rSNP enrichment in GWAS variants (Table 1) shows the highest enrichment for the rSNPs associated with blood cell quantitative characteristics, which is expected because we used PBMCs. Among the diseases, the largest number of rSNPs were associated with T2DM, matching the published data on the involvement of mononuclear cells in the pathogenesis of this disease [111,112,113]. In addition, the rSNPs associated with several morphometric characteristics (body mass index, body shape index, and waist-to-hip ratio) were also identified; these rSNPs were mainly localized to the promoters of the genes involved in inflammatory processes [40,43,44,45], insulin resistance [53,54], obesity [40,41], and carbohydrate transport [62], suggesting their involvement on T2DM pathogenesis.
To further clarify the potential role of the rSNPs of our panel in the mechanisms of T2DM pathogenesis and its complications, we searched for rSNPs in the promoter regions of the genes the expression of which changed during T2DM development using the largest RNA-seq dataset of the GEO repository (GSE221521) for the PBMCs of healthy subjects and T2DM patients without and with diabetic retinopathy. Retinopathy is a T2DM complication occurring in approximately half of diabetes patients [114]. We did not discover any DEGs when comparing the healthy donors and the T2DM subjects without retinopathy. Presumably, this is associated with the effects of the therapy received by the patients (unfortunately, the authors of the paper did not indicate which particular therapy was used) and the selection of the controls, which included the individuals with HbA1c = 5.80%, typical to a prediabetic state according to certain data [115]. However, 4612 DEGs were identified when comparing the healthy cohort to the T2DM subject with retinopathy and 1284 of them contained 3810 rSNPs in their promoters. The PPI network for these DEGs demonstrated that the genes of key modules were involved in the insulin signal transduction pathways and immune response. In particular, PIK3R2 and PIK3CG, represented in the first module code for the PI3K (phosphoinositide 3-kinase) subunit, an upstream regulator in the PI3K/AKT pathway. PI3K mediates the activation of the key enzyme in this pathway, Akt kinase (AKT serine/threonine kinase), promoting the glucose intake by cells via mediated activation of glucose transporters (GLUT) [116]; activation of glycogenesis via GSK-3 (glycogen synthase kinase 3) inhibition [117]; and negative regulation of gluconeogenesis by downregulating the expression of glucose-6-phosphatase [71] and phosphoenolpyruvate carboxykinase [118]. Note that Akt1 is a key upregulated hub gene in the PPI network based on the highest values of topological parameters. The increase in Akt1 expression in the PBMCs of the T2DM subjects as compared with healthy individuals was also reported by Manoel-Caetano et al. [119]. TP53 (tumor protein p53), identified as the most significant hub gene, also plays an important role in inflammation and T2DM development. This gene codes for transcription factor p53, which as a rule inhibits inflammatory responses. In addition, p53 is able to directly bind to Glut1 and Glut4, thereby directly inhibiting their transcription and promoting the development of insulin resistance. The upregulation of TP53 expression is associated with a positive regulation of gluconeogenesis (via expression upregulation of glucose-6-phosphatase, phosphoenolpyruvate carboxylase, and so on) and a negative regulation of glycolysis (decreasing the levels of fructose-2,6-bisphosphate, phosphoglycerate mutase, glucose-6-phosphate dehydrogenase, and so on) [72]. Thus, our data suggest that the rSNPs influencing gene expression in the considered modules are most likely associated with T2DM development by interfering with glucose metabolism and immune response. As for the rSNPs in hub gene promoters, they can well be the key regulators of these processes and the alteration in their expression can influence their genes/target proteins as well. We distinguished 117 rSNPs harbored in hub gene promoters with a high potential of T2DM prediction in the analyzed dataset according to ROC analysis data. The rSNPs of this type are of the greatest interest for further studies.
We used the formed rSNP panel to analyze the RNA-seq data for the blood cells of T2DM metformin responders and non-responders, which allowed for new insight into the pharmacogenetic factors influencing the metformin response in this disease. Notably, the research into genetic foundations of individual susceptibility to drugs has so far mainly dealt with the search for variations in genes coding for phase I and II drug-metabolizing enzymes, drug transporters, and some upstream transcription regulators of known pharmacogenes [120,121,122]. However, a rather large pool of heterogeneous data on SNPs in the genes belonging to a drug-metabolizing system [123,124,125] sets the challenge to develop the approaches to a systematized search for such variants.
It is known that metformin, a first-line drug for T2DM, is ineffective for approximately 30% of the patients [126]. Association studies succeeded in detecting several tens of the SNPs associated with the individual response to metformin. Since metformin is not metabolized in the human body [127], these studies mainly focused on the genes coding for organic cation transporters (OCTs) and multi-antimicrobial extrusion (MATE) proteins. As a result, numerous SNPs associated with metformin pharmacokinetics and pharmacodynamics in the protein-coding and noncoding regions of several genes of these groups were discovered. These genes include SLC22A1 and SLC22A2, coding for OKT1 and OKT2; SLC47A1 and SLC47A2, coding for MATE1 and MATE2-K; SLC2A, coding for one of the key glucose transporters; and the gene encoding the SP1 transcription factor, which modulates the expression of metformin transporters [127,128,129,130,131]. GWAS supplemented this list with a number of genes the products of which perform other, first and foremost, regulatory functions. In particular, the variations associated with metformin response in the genes of the AMPK signaling pathway, namely, STK11, PRKAA1, and PRKAA2, were found [132]. The SNPs associated with the metformin response are also discovered near the ATM gene, the product of which is involved in the redox homeostasis in the cell, as well as in the PRPF31 (pre-mRNA processing factor 31) gene; CPA6 gene, coding for carboxypeptidase A6), involved in the regulation of Akt/mTOR signaling pathway; and STAT3 gene, coding for the transcription factor actively involved in the regulation of metabolic and immune processes [133]. However, the use of GWAS in pharmacogenomics research requires a rather large cohort of patients with different responses to a drug. A tremendous number of currently used drugs and frequent cases of complex therapies considerably hinder obtaining mass data with the help of this approach [134,135]. An integrated analysis of the current omics functional data can help in resolving this problem, at least in part.
In this study, we, on the one hand, used our panel of 14,796 rSNPs constructed by detecting allele-specific events in ChIP-seq and RNA-seq data for the PBMCs of nine healthy individuals and, on the other hand, the transcriptome data (RNA-seq, GSE153315) for the blood cells of T2DM non-responders to metformin (after 3-month medication, N = 10) and the T2DM responders (N = 10) [78]. An integrated analysis of these data allowed us to obtain a new, wider list of the SNPs potentially associated with the individual response to metformin. Moreover, it has emerged that a considerable number of such SNPs reside in the promoter regions of the genes coding for transcription factors (FOXP1, POU2F2, YY-1, KLF6, NFIX, GMEB1, TSHZ2, ZNF121, ZNF44; ZNF440, ZNF611, ZNF714, ZNF793, and ZNF813) and other regulatory proteins (TAF3, NCOR1, SMARCA4, and PCBD2). The fact that rSNPs are discovered in these genes is important because of the possible modulation of an extensive reaction to the drug because these genes control different networks comprising numerous functionally linked genes. In particular, the YY-1 transcription factor is the best-studied of the mentioned transcription factors. It is ubiquitously expressed in many cell types and is able to activate or repress transcription of the networks comprising the genes involved in cell survival, replication, differentiation, metabolism, and inflammation [136,137]. The multifunctionality of YY-1 also appears as its involvement in the enhancer–promoter and super-enhancer–promoter interactions [138,139], splicing regulation [140], and transcription initiation by binding to the initiator element (INR, a core promoter) [141]. YY-1 has been shown to act in numerous processes the disturbance of which is associated with the development of diabetes and its complications as well as the regulation of systemic inflammation [83,142,143,144]. This suggests that the rSNPs discovered in the promoter-regulatory region of this gene are also potential contributors to the implementation of metformin therapeutic effects by altering its expression level and, as a consequence, the expression of YY-1 target genes. The assumption on an important role of YY-1 in the metformin response agrees well with the results of our analysis of ANANASTRA data [30,31] that the promoter regions of many genes differentially expressed in response to metformin carry the rSNPs influencing YY-1 binding.

4. Materials and Methods

4.1. Subjects of Investigation

The study involved nine somatically healthy individuals living on the territory of the Russian Federation (Novosibirsk, Russia) and displaying no diseases in their life history by the moment of sampling their biological material. The sample comprised males (N = 3) and females (N = 6) aged 26–42 (mean age, 33 years) (Supplementary Table S8). The study was approved by the ethical committee of the Institute of Therapy and Preventive Medicine, Novosibirsk, Russia. Each participant signed the informed consent.

4.2. Isolation of Peripheral Blood Mononuclear Cells

Blood (20 mL) was sampled from the cubital vein of each subject and diluted with an equal volume of phosphate-buffered saline (PBS). Ficoll (12 mL; ρ = 1.077, Biolot, Saint Petersburg, Russia) was added to two clean Falcon tubes to gently layer the diluted blood (20 mL per each Falcon tube) onto the Ficoll solution with a pipette. The samples were centrifuged at 400× g for 15 min at 4 °C. The layer of mononuclear cells from each Falcon was transferred to a new tube, supplemented with PBS (4 mL), mixed, and centrifuged (200× g/7 min/4 °C). The supernatant was removed and the pellet was supplemented with the buffer for lysing erythrocytes (20 mL; 150 mM NH4Cl and 0.01 mM EDTA) to incubate for 10 min at room temperature and centrifuge (200× g/7 min/4 °C). The supernatant was removed and the sediment was washed twice with PBS (3 mL) and centrifuged (200× g/7 min/4 °C). The purified cells from one Falcon tube were used to isolate RNA and from the other one, in the experiments on chromatin immunoprecipitation.

4.3. mRNA Sequencing

The purified cells obtained from each donor were homogenized in Trizol buffer (1 mL; Ambion, Austin, TX, USA), supplemented with chloroform (200 μL), and centrifuged (12,000× g/15 min/4 °C). The aqueous phase was transferred to a new tube, supplemented with 1/5 volume of chloroform, and centrifuged (12,000× g/5 min/4 °C) to transfer the supernatant into a new tube, precipitate it with one volume of isopropanol and 1 μL linear polyacrylamide (BIORON, Römerberg, Germany), and centrifuge (20,238× g/15 min/4 °C). The supernatant was collected to precipitate RNA with 75% alcohol (1 mL) and centrifuge (20,238× g/5 min/4 °C). The alcohol was removed and the pellet was dried at room temperature for 5–10 min to dissolve in 25 μL of ddH2O. RNA concentration was measured in a NanoDrop (Thermo Fisher Scientific, Waltham, MA, USA) spectrophotometer. The complementary DNA (cDNA) libraries were prepared using the NEBNext® Ultra II Directional RNA Library Prep Kit for Illumina® (New England Biolabs, Ipswich, MA, USA) according to the manufacturer’s protocol.

4.4. Chromatin Immunoprecipitation

The cell sediment was dissolved in PBS (5 mL) and 37% formaldehyde (138.8 μL; Sigma, Saint Louis, MO, USA), homogenized, and mixed in a rotator for 10 min at room temperature. The reaction was stopped with 2 M glycine (342 μL), mixed for 5 min in a rotator at room temperature, and centrifuged (200× g/5 min/4 °C); the supernatant was removed. The pellet was twice washed with PBS (3 mL) and supplemented with 800 μL of lysis buffer (10 mM Tris pH 7.5, 1 mM EDTA, 0.5% SDS, 0.1% sodium deoxycholate, and 1% Triton X-100) and 100X protease inhibitor cocktail (Thermo Fisher Scientific, Waltham, MA, USA) to break the chromatin into the fragments of 150–600 bp in an immersion homogenizer (BANDELIN electronic GmbH & Co. KG, Berlin, Germany). The resulting fragments were centrifuged (3000× g/10 min/4 °C) and the supernatant was added to A- and G-protein beads (20 μL; New England Biolabs, Ipswich, MA, USA) and supplemented with 1.5 V of the lysis buffer without SDS and 100X protease inhibitor cocktail. In parallel, G-protein beads (40 μL), lysis buffer (160 μL), lysis buffer without SDS (240 μL), and the antibodies (8 μL) to histone modification H3K4me3 or H3K27ac (Abcam, United Kingdom) were placed into a separate tube to incubate it together with the tube containing chromatin and A- and G-protein beads at a constant rotation for 2 h at 4 °C. The tubes were placed onto a magnetic base and then the supernatant was removed from the tube with antibodies and the sediment was supplemented with the supernatant after the incubation with A- and G-protein beads. The immunoprecipitation reaction was conducted overnight at +4 °C in a rotator. In the morning, the tube was returned onto magnetic base to remove the supernatant and washed twice with an addition of the lysis buffer (1 mL), 1 mL of high salt buffer (10 mM Tris pH 7.5, 1 mM EDTA, 0.1% SDS, 0.1% sodium deoxycholate, 1% Triton X-100, and 0.5 M NaCl), 1 mL of TE buffer (10 mM Tris pH 7.5, 1 mM EDTA, and 1% Triton X-100), and 1 mL of the TE buffer without Triton X-100. After removal of the supernatant, the tube was supplemented with 100 μL of elution buffer (1% SDS, 10 mM Tris pH 7.5, and 1 mM EDTA) and incubated for 1 h at +65 °C in a thermoshaker (800 rpm). The tube was placed onto a magnetic base and the supernatant was transferred into a new tube (the washing with elution buffer was repeated twice). To prevent chromatin cross-linking, the sample was supplemented with 5M NaCl (to a final concentration of 370 mM) and incubated in a thermostat at +65 °C for 4 h. The sample was supplemented with 2.46 volume of ddH2O, 0.5 M EDTA (1/50 of the sample volume), 1 M Tris-HCl pH 6.5 (1/16 of the sample volume), 10% SDS (1/13 of the sample volume) and 1 μL of proteinase K (800 U/mL; New England Biolabs, Ipswich, MA, USA) to incubate in a thermostat at +52 °C for 3 h. Then, RNase A (20 mg/mL; Thermo Fisher Scientific, Waltham, MA, USA) was added at the last 20 min. The sample was taken off the thermostat and supplemented with 1 volume of phenol–chloroform pH 8.0 (1:1) and centrifuged (1000× g/5 min) to transfer the aqueous phase to a new tube. Then, 1 volume of chloroform was added and centrifuged (1000× g/5 min). The aqueous phase was collected and supplemented with 1/10 V of 3 M NaAc, 1 volume of isopropanol, and 1 μL of linear polyacrylamide to incubate for 5–7 h at −70 °C followed by centrifugation (10,000× g/20 min/4 °C), removal of the supernatant, precipitation with 75% ethanol, and centrifugation (10,000× g/5 min/4 °C). The supernatant was removed and the pellet was dissolved in ddH2O (20 μL). DNA libraries were prepared analogously to the protocol for cDNA libraries.

4.5. Sequencing and Aligning to Reference Genome

Paired-end sequencing of cDNA and DNA libraries was performed using the NovaSeq 6000 or MGISEQ-2000RS platforms. The Illumina or DNBSEQ adapters and the reads shorter than 36 bp were discarded from the ChIP-seq and RNA-seq data with Trimmomatic v. 3.2.2 [145]. The resulting data were aligned to reference genome assembly GRCh38/hg38 (Genome Reference Consortium Human Build 38), available at the NCBI FTP site [146]. The forward and reverse reads (145–150 bp) were aligned with Bowtie 2 [147]; the threshold value of mapping quality (MAPQ) was >20. The aligned reads were transferred into bam format.

4.6. Search for Heterozygous Positions

To search for heterozygous positions, all ChIP-seq and RNA-seq files for each individual were pooled with the help of the bcftools merge command [148]. A position was regarded as heterozygous if it met the following conditions: coverage depth, >20 and frequency of minor allele, ≥0.2 of the total coverage of the corresponding position. The heterozygous positions were selected using bcftools mpileup [148].

4.7. Assembling Alternative Genome

In order to avoid the alignment bias towards the reads containing reference alleles (resulting from the mapping algorithms that discriminate the alternative allele), an additional individual alternative genome was assembled for each person; in this genome, the reference allele was replaced with the alternative one for each heterozygous position. Then, the ChIP-seq and RNA-seq file carrying the reads for each individual was aligned to the individual alternative sequence.

4.8. Computing Allelic Asymmetry

Allelic asymmetry for all heterozygous positions in the ChIP-seq and RNA-seq data was computed using the binomial test against the null hypothesis on an equal coverage of heterozygous positions by the reference reads aligned to reference genome and the alternative reads aligned to alternative genome with subsequent correction for multiple comparison according to Benjamini–Hochberg (p.adj. < 0.1). The rs ID numbers were assigned to the computed allele-asymmetric positions based on the dbSNP Build 155 [149,150]. Allele-asymmetric SNPs from ChIP-seq data were considered to be associated with allele-specific binding (ASB) events, and allele-asymmetric SNPs from RNA-seq data were considered to be associated with allele-specific expression (ASE) events.

4.9. Searching for Allele-Asymmetric SNPs in the Promoters of the Genes Displaying Allele-Specific Expression

The list of ASB SNPs for the ChIP-seq data was intersected with gene promoter regions. The region of ±1000 bp from the transcription start site of all transcripts (including alternative transcription start sites) contained in the TxDb.Hsapiens.UCSC.hg38.knownGene R annotation package for TxDb object(s) [151] was regarded as promoter. The list of ASE SNPs based on the RNA-seq data was intersected with the coordinates of the genes contained in the TxDb.Hsapiens.UCSC.hg38.knownGene package. The genes containing at least one ASE SNP were regarded as expressed in an allele-specific manner. The total allele-specific expression was computed for the genes carrying more than one ASE SNP using the MBASED approach [152]; this method is based on the merging of the information about allelic asymmetry for each heterozygous SNP in a particular gene. The final stage consisted of the intersection of the ASB SNPs localized to promoters and the target genes with the computed allele-specific expression. The promoter ASB SNPs, for which the target genes are expressed in an allele-specific manner, were named regulatory SNPs (rSNPs).

4.10. Characterizing rSNPs with the Help of Open-Access Data

The data of the GTEx project on eQTLs [28,153], GWAS Catalog on the association of genome loci and phenotypic traits [26,27], and ANANASTRA web server for annotation of the loci of allele-specific transcription factor binding in ChIP-seq data [30,31] were used to characterize the constructed rSNP panel.

4.11. Analyzing DEGs

DEGs were searched for using the DeSeq2 (version 1.44) R package [154] in the RNA-seq data deposited with the Gene Expression Omnibus database [155,156]. These data were obtained for the whole blood samples of the cohorts of (1) healthy subjects (n = 50) and T2DM cases with diagnosed (n = 69) and not diagnosed (n = 74) diabetic retinopathy (GSE221521) [63] and (2) healthy subjects (n = 10) and T2DM responders (n = 10)/non-responders (n = 10) to 3-month metformin treatment (GSE153315) [78]. According to [78] the patients who had a decrease in glycated hemoglobin (HbA1c) level > 1% or a decrease of >20 mg/dL in fasting blood glucose level from baseline after three months of therapy were considered as responders.
Multiple testing correction according to Benjamini–Hochberg procedure setting the false discovery rate (FDR) of <0.05 was used. Genes were regarded as differentially expressed if |log2FC| > 0.2 and p.adj. < 0.05. Then, DEGs were intersected with the genes harboring rSNPs from the constructed panel in their promoters.

4.12. Search for Hub Genes

The hub genes were searched for applying a network of protein interactions constructed with the help of STRING [32,33], which predicts protein–protein interactions (PPI) based on the data on physical (experimental data) and functional (computational predictions, co-expressions, and previous knowledge in databases) associations. The PPI network was separately constructed for the upregulated and downregulated DEG groups, which harbored rSNPs in their promoters. Then, three topological characteristics—degree, stress centrality, and betweenness centrality—were computed using CytoNCA [157,158], a Cytoscape plugin [159,160]. The lists of the top 20 proteins with the highest values of each computed characteristic were formed and merged. The proteins in the consolidated list were referred to as hub proteins and the corresponding genes, as hub genes.

4.13. Identification of Protein Modules

The significant modules in the analyzed PPI networks were identified using the MCODE tool [161,162]; this tool clusters the densely connected regions in large protein–protein interaction networks. The following criteria were used for selecting the significant modules: max depth = 100, degree cut-off = 2, node score cut-off = 0.2, MCODE scores ≥ 4, and K-score = 2.

4.14. Enrichment Analysis of KEGG/Reactome/GO

Functional enrichment analysis of all significant modules was performed. The terms from the GO Knowledgebase [36,37] together with the pathways from the KEGG [34,35] and Reactome Knowledgebase [38,39] were used as functional terms. The enrichment was analyzed using over-representation analysis (ORA) of the GO terms and KEGG/Reactome pathways in the protein modules with the R clusterProfiler package (version 4.12.5) [163]. After the Benjamini–Hochberg procedure, p.adj. < 0.05 was set as the enrichment threshold value.

4.15. ROC Analysis

ROC (Receiver Operator Characteristic) curves were constructed and AUC (Area Under Curve) values were computed from the normalized expression values for each hub gene calculated using DESeq2. The hub genes with AUC > 0.7 were regarded as good predictors of T2DM in the studied transcriptome dataset. The ROC curves were constructed using the pROC package (version 1.18.5) for R [164].

5. Conclusions

The search for allele-specific events in the experimental data on gene expression profiles (RNA-seq) and active chromatin marks H3K4me3 and H3K27ac (ChIP-seq) in the peripheral blood mononuclear cells of nine healthy donors allowed us to construct the panel comprising 14,848 rSNPs in the promoter regions of 5132 genes. Over 38% of the discovered rSNPs are associated with different GWAS phenotypic traits and 66.7%, with GTEx eQTLs, which uncovers the molecular sense of a large GWAS-identified SNP group and is an independent confirmation for the regulatory function of a considerable group of variants from GTEx catalogs. The rSNPs potentially associated with T2DM development and/or diabetic retinopathy as well as with individual specific features of the response to the antidiabetic drug metformin are identified. The genes carrying these rSNPs in their promoters are enriched in the pathways related to glucose metabolism and inflammation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25179297/s1.

Author Contributions

I.S.D.: conceptualization, data curation, formal analysis, investigation, validation, visualization, writing—original draft preparation, writing—review and editing. E.E.K.: validation, writing—original draft preparation, writing—review and editing, funding acquisition, project administration, supervision. E.Y.R.: funding acquisition, writing—original draft, supervision. T.I.M.: conceptualization, writing—original draft preparation, writing—review and editing, funding acquisition, project administration, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the Russian Science Foundation (23-15-00113).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the ethical committee with the Institute of Therapy and Preventive Medicine, Novosibirsk, Russia (18-29-09041).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The original data presented in the study are openly available in Sequence Read Archive (SRA) at https://www.ncbi.nlm.nih.gov/sra, PRJNA1120824 (accessed on 6 June 2024).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reed, J.; Bain, S.; Kanamarlapudi, V. A Review of Current Trends with Type 2 Diabetes Epidemiology, Aetiology, Pathogenesis, Treatments and Future Perspectives. Diabetes Metab. Syndr. Obes. 2021, 14, 3567–3602. [Google Scholar] [CrossRef] [PubMed]
  2. Ingelsson, E.; McCarthy, M.I. Human Genetics of Obesity and Type 2 Diabetes Mellitus. Circ. Genom. Precis. Med. 2018, 11, e002090. [Google Scholar] [CrossRef] [PubMed]
  3. 2. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes—2019. Diabetes Care 2019, 42, S13–S28. [CrossRef] [PubMed]
  4. Ruze, R.; Liu, T.; Zou, X.; Song, J.; Chen, Y.; Xu, R.; Yin, X.; Xu, Q. Obesity and Type 2 Diabetes Mellitus: Connections in Epidemiology, Pathogenesis, and Treatments. Front. Endocrinol. 2023, 14, 1161521. [Google Scholar] [CrossRef]
  5. 2. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes—2018. Diabetes Care 2018, 41, S13–S27. [CrossRef]
  6. Mandrup-Poulsen, T. Type 2 Diabetes Mellitus. Dermatol. Clin. 2013, 31, 495–506. [Google Scholar] [CrossRef]
  7. Eizirik, D.L.; Pasquali, L.; Cnop, M. Pancreatic β-Cells in Type 1 and Type 2 Diabetes Mellitus: Different Pathways to Failure. Nat. Rev. Endocrinol. 2020, 16, 349–362. [Google Scholar] [CrossRef] [PubMed]
  8. Mota, M.; Banini, B.A.; Cazanave, S.C.; Sanyal, A.J. Molecular Mechanisms of Lipotoxicity and Glucotoxicity in Nonalcoholic Fatty Liver Disease. Metabolism 2016, 65, 1049–1061. [Google Scholar] [CrossRef] [PubMed]
  9. Pezhman, L.; Tahrani, A.; Chimen, M. Dysregulation of Leukocyte Trafficking in Type 2 Diabetes: Mechanisms and Potential Therapeutic Avenues. Front. Cell Dev. Biol. 2021, 9, 624184. [Google Scholar] [CrossRef]
  10. Torres-Castro, I.; Arroyo-Camarena, Ú.D.; Martínez-Reyes, C.P.; Gómez-Arauz, A.Y.; Dueñas-Andrade, Y.; Hernández-Ruiz, J.; Béjar, Y.L.; Zaga-Clavellina, V.; Morales-Montor, J.; Terrazas, L.I.; et al. Human Monocytes and Macrophages Undergo M1-Type Inflammatory Polarization in Response to High Levels of Glucose. Immunol. Lett. 2016, 176, 81–89. [Google Scholar] [CrossRef]
  11. Campbell, N.R.C.; Leiter, L.A.; Larochelle, P.; Tobe, S.; Chockalingam, A.; Ward, R.; Morris, D.; Tsuyuki, R. Hypertension in Diabetes: A Call to Action. Can. J. Cardiol. 2009, 25, 299–302. [Google Scholar] [CrossRef]
  12. Beckman, J.A.; Paneni, F.; Cosentino, F.; Creager, M.A. Diabetes and Vascular Disease: Pathophysiology, Clinical Consequences, and Medical Therapy: Part II. Eur. Heart J. 2013, 34, 2444–2452. [Google Scholar] [CrossRef] [PubMed]
  13. Tomic, D.; Shaw, J.E.; Magliano, D.J. The Burden and Risks of Emerging Complications of Diabetes Mellitus. Nat. Rev. Endocrinol. 2022, 18, 525–539. [Google Scholar] [CrossRef]
  14. Nguyen, T.T.; Ta, Q.T.H.; Nguyen, T.K.O.; Nguyen, T.T.D.; Van Giau, V. Type 3 Diabetes and Its Role Implications in Alzheimer’s Disease. Int. J. Mol. Sci. 2020, 21, 3165. [Google Scholar] [CrossRef] [PubMed]
  15. Rohm, T.V.; Meier, D.T.; Olefsky, J.M.; Donath, M.Y. Inflammation in Obesity, Diabetes, and Related Disorders. Immunity 2022, 55, 31–55. [Google Scholar] [CrossRef] [PubMed]
  16. Ding, Q.; Tan, A.L.M.; Parra, E.J.; Cruz, M.; Sim, X.; Teo, Y.-Y.; Long, J.; Alsafar, H.; Petretto, E.; Tai, E.-S.; et al. Genome-Wide Meta-Analysis Associates GPSM1 with Type 2 Diabetes, a Plausible Gene Involved in Skeletal Muscle Function. J. Hum. Genet. 2020, 65, 411–420. [Google Scholar] [CrossRef]
  17. López Rodríguez, M.; Kaminska, D.; Lappalainen, K.; Pihlajamäki, J.; Kaikkonen, M.U.; Laakso, M. Identification and Characterization of a FOXA2-Regulated Transcriptional Enhancer at a Type 2 Diabetes Intronic Locus That Controls GCKR Expression in Liver Cells. Genome Med. 2017, 9, 63. [Google Scholar] [CrossRef]
  18. Locke, J.M.; Hysenaj, G.; Wood, A.R.; Weedon, M.N.; Harries, L.W. Targeted Allelic Expression Profiling in Human Islets Identifies Cis-Regulatory Effects for Multiple Variants Identified by Type 2 Diabetes Genome-Wide Association Studies. Diabetes 2015, 64, 1484–1491. [Google Scholar] [CrossRef]
  19. Shojima, N.; Yamauchi, T. Progress in Genetics of Type 2 Diabetes and Diabetic Complications. J. Diabetes Investig. 2023, 14, 503–515. [Google Scholar] [CrossRef]
  20. Mahajan, A.; Taliun, D.; Thurner, M.; Robertson, N.R.; Torres, J.M.; Rayner, N.W.; Payne, A.J.; Steinthorsdottir, V.; Scott, R.A.; Grarup, N.; et al. Fine-Mapping Type 2 Diabetes Loci to Single-Variant Resolution Using High-Density Imputation and Islet-Specific Epigenome Maps. Nat. Genet. 2018, 50, 1505–1513. [Google Scholar] [CrossRef]
  21. Del Bosque-Plata, L.; Hernández-Cortés, E.P.; Gragnoli, C. The Broad Pathogenetic Role of TCF7L2 in Human Diseases beyond Type 2 Diabetes. J. Cell Physiol. 2022, 237, 301–312. [Google Scholar] [CrossRef] [PubMed]
  22. Geoghegan, G.; Simcox, J.; Seldin, M.M.; Parnell, T.J.; Stubben, C.; Just, S.; Begaye, L.; Lusis, A.J.; Villanueva, C.J. Targeted Deletion of Tcf7l2 in Adipocytes Promotes Adipocyte Hypertrophy and Impaired Glucose Metabolism. Mol. Metab. 2019, 24, 44–63. [Google Scholar] [CrossRef] [PubMed]
  23. Gaulton, K.J.; Nammo, T.; Pasquali, L.; Simon, J.M.; Giresi, P.G.; Fogarty, M.P.; Panhuis, T.M.; Mieczkowski, P.; Secchi, A.; Bosco, D.; et al. A Map of Open Chromatin in Human Pancreatic Islets. Nat. Genet. 2010, 42, 255–259. [Google Scholar] [CrossRef]
  24. Del Bosque-Plata, L.; Martínez-Martínez, E.; Espinoza-Camacho, M.Á.; Gragnoli, C. The Role of TCF7L2 in Type 2 Diabetes. Diabetes 2021, 70, 1220–1228. [Google Scholar] [CrossRef] [PubMed]
  25. Yan, J.; Qiu, Y.; Ribeiro dos Santos, A.M.; Yin, Y.; Li, Y.E.; Vinckier, N.; Nariai, N.; Benaglio, P.; Raman, A.; Li, X.; et al. Systematic Analysis of Binding of Transcription Factors to Noncoding Variants. Nature 2021, 591, 147–151. [Google Scholar] [CrossRef] [PubMed]
  26. GWAS. Available online: https://www.ebi.ac.uk/gwas/ (accessed on 25 March 2023).
  27. Sollis, E.; Mosaku, A.; Abid, A.; Buniello, A.; Cerezo, M.; Gil, L.; Groza, T.; Güneş, O.; Hall, P.; Hayhurst, J.; et al. The NHGRI-EBI GWAS Catalog: Knowledgebase and Deposition Resource. Nucleic Acids Res. 2023, 51, D977–D985. [Google Scholar] [CrossRef]
  28. GTEx. Available online: https://www.gtexportal.org/home/ (accessed on 25 March 2023).
  29. GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020, 369, 1318–1330. [Google Scholar] [CrossRef]
  30. ANANASTRA: ANnotation and Enrichment ANalysis of Allele-Specific TRAnscription Factor Binding at SNPs. Available online: https://ananastra.autosome.org/ (accessed on 14 January 2024).
  31. Boytsov, A.; Abramov, S.; Aiusheeva, A.Z.; Kasianova, A.M.; Baulin, E.; Kuznetsov, I.A.; Aulchenko, Y.S.; Kolmykov, S.; Yevshin, I.; Kolpakov, F.; et al. ANANASTRA: Annotation and Enrichment Analysis of Allele-Specific Transcription Factor Binding at SNPs. Nucleic Acids Res. 2022, 50, W51–W56. [Google Scholar] [CrossRef]
  32. STRING. Available online: https://string-db.org/ (accessed on 14 January 2024).
  33. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING V11: Protein–Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef]
  34. KEGG: Kyoto Encyclopedia of Genes and Genomes. Available online: https://www.genome.jp/kegg/ (accessed on 14 January 2024).
  35. Kanehisa, M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  36. The Gene Ontology Resource. Available online: https://geneontology.org/ (accessed on 14 January 2024).
  37. Aleksander, S.A.; Balhoff, J.; Carbon, S.; Cherry, J.M.; Drabkin, H.J.; Ebert, D.; Feuermann, M.; Gaudet, P.; Harris, N.L.; Hill, D.P.; et al. The Gene Ontology Knowledgebase in 2023. Genetics 2023, 224, iyad031. [Google Scholar] [CrossRef] [PubMed]
  38. Reactome. Available online: https://reactome.org/ (accessed on 14 January 2024).
  39. Milacic, M.; Beavers, D.; Conley, P.; Gong, C.; Gillespie, M.; Griss, J.; Haw, R.; Jassal, B.; Matthews, L.; May, B.; et al. The Reactome Pathway Knowledgebase 2024. Nucleic Acids Res. 2024, 52, D672–D678. [Google Scholar] [CrossRef] [PubMed]
  40. Matthews, J.; Villescas, S.; Herat, L.; Schlaich, M.; Matthews, V. Implications of ADAM17 Activation for Hyperglycaemia, Obesity and Type 2 Diabetes. Biosci. Rep. 2021, 41, BSR20210029. [Google Scholar] [CrossRef] [PubMed]
  41. Kaneko, H.; Anzai, T.; Horiuchi, K.; Morimoto, K.; Anzai, A.; Nagai, T.; Sugano, Y.; Maekawa, Y.; Itoh, H.; Yoshikawa, T.; et al. Tumor Necrosis Factor-α Converting Enzyme Inactivation Ameliorates High-Fat Diet-Induced Insulin Resistance and Altered Energy Homeostasis. Circ. J. 2011, 75, 2482–2490. [Google Scholar] [CrossRef]
  42. Menghini, R.; Fiorentino, L.; Casagrande, V.; Lauro, R.; Federici, M. The Role of ADAM17 in Metabolic Inflammation. Atherosclerosis 2013, 228, 12–17. [Google Scholar] [CrossRef]
  43. Song, X.; Zhang, L.; Du, X.; Zheng, Y.; Jia, T.; Zhou, T.; Che, D.; Geng, S. Neuroblast Differentiation-Associated Protein Derived Polypeptides: AHNAK(5758-5775) Induces Inflammation by Activating Mast Cells via ST2. Immunol. Investig. 2023, 52, 178–193. [Google Scholar] [CrossRef]
  44. Li, L.; Liu, Y.; Feng, T.; Zhou, W.; Wang, Y.; Li, H. The AHNAK Induces Increased IL-6 Production in CD4+ T Cells and Serves as a Potential Diagnostic Biomarker for Recurrent Pregnancy Loss. Clin. Exp. Immunol. 2022, 209, 291–304. [Google Scholar] [CrossRef]
  45. Parikh, D.; Riascos-Bernal, D.F.; Egaña-Gorroño, L.; Jayakumar, S.; Almonte, V.; Chinnasamy, P.; Sibinga, N.E.S. Allograft Inflammatory Factor-1-like Is Not Essential for Age Dependent Weight Gain or HFD-Induced Obesity and Glucose Insensitivity. Sci. Rep. 2020, 10, 3594. [Google Scholar] [CrossRef]
  46. Piotrowska, K.; Słuczanowska-Głabowska, S.; Kurzawski, M.; Dziedziejko, V.; Kopytko, P.; Paczkowska, E.; Rogińska, D.; Safranow, K.; Machaliński, B.; Pawlik, A. Over-Expression of Allograft Inflammatory Factor-1 (AIF-1) in Patients with Rheumatoid Arthritis. Biomolecules 2020, 10, 1064. [Google Scholar] [CrossRef]
  47. Kelemen, S.E.; Autieri, M. V Expression of Allograft Inflammatory Factor-1 in T Lymphocytes: A Role in T-Lymphocyte Activation and Proliferative Arteriopathies. Am. J. Pathol. 2005, 167, 619–626. [Google Scholar] [CrossRef]
  48. Ren, L.; Du, W.; Song, D.; Lu, H.; Hamblin, M.H.; Wang, C.; Du, C.; Fan, G.-C.; Becker, R.C.; Fan, Y. Genetic Ablation of Diabetes-Associated Gene Ccdc92 Reduces Obesity and Insulin Resistance in Mice. iScience 2023, 26, 105769. [Google Scholar] [CrossRef] [PubMed]
  49. Wang, G.; Wu, B.; Zhang, B.; Wang, K.; Wang, H. LncRNA CTBP1-AS2 Alleviates High Glucose-Induced Oxidative Stress, ECM Accumulation, and Inflammation in Diabetic Nephropathy via MiR-155-5p/FOXO1 Axis. Biochem. Biophys. Res. Commun. 2020, 532, 308–314. [Google Scholar] [CrossRef] [PubMed]
  50. Wang, X.; Liu, Y.; Rong, J.; Wang, K. LncRNA HCP5 Knockdown Inhibits High Glucose-Induced Excessive Proliferation, Fibrosis and Inflammation of Human Glomerular Mesangial Cells by Regulating the MiR-93-5p/HMGA2 Axis. BMC Endocr. Disord. 2021, 21, 134. [Google Scholar] [CrossRef] [PubMed]
  51. Fu, L.; Shi, J.; Hu, K.; Wang, J.; Wang, W.; Ke, X. Mitogen-Activated Protein Kinase Binding Protein 1 (MAPKBP1) Is an Unfavorable Prognostic Biomarker in Cytogenetically Normal Acute Myeloid Leukemia. Oncotarget 2015, 6, 8144–8154. [Google Scholar] [CrossRef]
  52. Xu, Q.; Yin, S.; Yao, Y.; Li, X.; Song, B.; Yang, Y.; Liu, Y.; Chen, R.; Li, J.; Ma, T.; et al. MAST3 Modulates the Inflammatory Response and Proliferation of Fibroblast-like Synoviocytes in Rheumatoid Arthritis. Int. Immunopharmacol. 2019, 77, 105900. [Google Scholar] [CrossRef]
  53. Seillier, M.; Pouyet, L.; N’Guessan, P.; Nollet, M.; Capo, F.; Guillaumond, F.; Peyta, L.; Dumas, J.-F.; Varrault, A.; Bertrand, G.; et al. Defects in Mitophagy Promote Redox-Driven Metabolic Syndrome in the Absence of TP53INP1. EMBO Mol. Med. 2015, 7, 802–818. [Google Scholar] [CrossRef] [PubMed]
  54. Tian, N.; Liu, Q.; Li, Y.; Tong, L.; Lu, Y.; Zhu, Y.; Zhang, P.; Chen, H.; Hu, L.; Meng, J.; et al. Transketolase Deficiency in Adipose Tissues Protects Mice From Diet-Induced Obesity by Promoting Lipolysis. Diabetes 2020, 69, 1355–1367. [Google Scholar] [CrossRef]
  55. Ramdas, M.; Harel, C.; Armoni, M.; Karnieli, E. AHNAK KO Mice Are Protected from Diet-Induced Obesity but Are Glucose Intolerant. Horm. Metab. Res. 2015, 47, 265–272. [Google Scholar] [CrossRef]
  56. Erfanian, S.; Mir, H.; Abdoli, A.; Roustazadeh, A. Association of Gastric Inhibitory Polypeptide Receptor (GIPR) Gene Polymorphism with Type 2 Diabetes Mellitus in Iranian Patients. BMC Med. Genom. 2023, 16, 44. [Google Scholar] [CrossRef]
  57. Okamoto, K.; Iwasaki, N.; Doi, K.; Noiri, E.; Iwamoto, Y.; Uchigata, Y.; Fujita, T.; Tokunaga, K. Inhibition of Glucose-Stimulated Insulin Secretion by KCNJ15, a Newly Identified Susceptibility Gene for Type 2 Diabetes. Diabetes 2012, 61, 1734–1741. [Google Scholar] [CrossRef]
  58. Wood, A.R.; Jonsson, A.; Jackson, A.U.; Wang, N.; van Leewen, N.; Palmer, N.D.; Kobes, S.; Deelen, J.; Boquete-Vilarino, L.; Paananen, J.; et al. A Genome-Wide Association Study of IVGTT-Based Measures of First-Phase Insulin Secretion Refines the Underlying Physiology of Type 2 Diabetes Variants. Diabetes 2017, 66, 2296–2309. [Google Scholar] [CrossRef] [PubMed]
  59. Bailetti, D.; Sentinelli, F.; Prudente, S.; Cimini, F.A.; Barchetta, I.; Totaro, M.; Di Costanzo, A.; Barbonetti, A.; Leonetti, F.; Cavallo, M.G.; et al. Deep Resequencing of 9 Candidate Genes Identifies a Role for ARAP1 and IGF2BP2 in Modulating Insulin Secretion Adjusted for Insulin Resistance in Obese Southern Europeans. Int. J. Mol. Sci. 2022, 23, 1221. [Google Scholar] [CrossRef] [PubMed]
  60. Uchida, T.; Nakamura, T.; Hashimoto, N.; Matsuda, T.; Kotani, K.; Sakaue, H.; Kido, Y.; Hayashi, Y.; Nakayama, K.I.; White, M.F.; et al. Deletion of Cdkn1b Ameliorates Hyperglycemia by Maintaining Compensatory Hyperinsulinemia in Diabetic Mice. Nat. Med. 2005, 11, 175–182. [Google Scholar] [CrossRef] [PubMed]
  61. Li, Y.; Chen, Y.; Liu, Z.; Lin, B.; Deng, X.; Xiao, Q.; Chen, Z.; Ye, H.; Chen, D.; Su, Y.; et al. Downregulation of Kcnq1ot1 Attenuates β-Cell Proliferation and Insulin Secretion via the MiR-15b-5p/Ccnd1 and Ccnd2 Axis. Acta Diabetol. 2022, 59, 885–899. [Google Scholar] [CrossRef]
  62. Körschen, H.G.; Yildiz, Y.; Raju, D.N.; Schonauer, S.; Bönigk, W.; Jansen, V.; Kremmer, E.; Kaupp, U.B.; Wachten, D. The Non-Lysosomal β-Glucosidase GBA2 Is a Non-Integral Membrane-Associated Protein at the Endoplasmic Reticulum (ER) and Golgi. J. Biol. Chem. 2013, 288, 3381–3393. [Google Scholar] [CrossRef]
  63. Xiang, Z.-Y.; Chen, S.-L.; Qin, X.-R.; Lin, S.-L.; Xu, Y.; Lu, L.-N.; Zou, H.-D. Changes and Related Factors of Blood CCN1 Levels in Diabetic Patients. Front. Endocrinol. 2023, 14, 1131993. [Google Scholar] [CrossRef]
  64. Wolf, G.; Schanze, A.; Stahl, R.A.K.; Shankland, S.J.; Amann, K. P27Kip1 Knockout Mice Are Protected from Diabetic Nephropathy: Evidence for P27Kip1 Haplotype Insufficiency. Kidney Int. 2005, 68, 1583–1589. [Google Scholar] [CrossRef]
  65. Yamada, S.; Senokuchi, T.; Matsumura, T.; Morita, Y.; Ishii, N.; Fukuda, K.; Murakami-Nishida, S.; Nishida, S.; Kawasaki, S.; Motoshima, H.; et al. Inhibition of Local Macrophage Growth Ameliorates Focal Inflammation and Suppresses Atherosclerosis. Arter. Thromb. Vasc. Biol. 2018, 38, 994–1006. [Google Scholar] [CrossRef]
  66. Sliwinska, A.; Kasznicki, J.; Kosmalski, M.; Mikolajczyk, M.; Rogalska, A.; Przybylowska, K.; Majsterek, I.; Drzewoski, J. Tumour Protein 53 Is Linked with Type 2 Diabetes Mellitus. Indian J. Med. Res. 2017, 146, 237. [Google Scholar] [CrossRef]
  67. Bengal, E.; Aviram, S.; Hayek, T. P38 MAPK in Glucose Metabolism of Skeletal Muscle: Beneficial or Harmful? Int. J. Mol. Sci. 2020, 21, 6480. [Google Scholar] [CrossRef]
  68. Saltiel, A.R. Insulin Signaling in Health and Disease. J. Clin. Investig. 2021, 131, 1710–1711. [Google Scholar] [CrossRef] [PubMed]
  69. Teaney, N.A.; Cyr, N.E. FoxO1 as a Tissue-Specific Therapeutic Target for Type 2 Diabetes. Front. Endocrinol. 2023, 14, 1286838. [Google Scholar] [CrossRef]
  70. Jeon, S.-M. Regulation and Function of AMPK in Physiology and Diseases. Exp. Mol. Med. 2016, 48, e245. [Google Scholar] [CrossRef] [PubMed]
  71. Huang, X.; Liu, G.; Guo, J.; Su, Z. The PI3K/AKT Pathway in Obesity and Type 2 Diabetes. Int. J. Biol. Sci. 2018, 14, 1483–1496. [Google Scholar] [CrossRef]
  72. Kung, C.-P.; Murphy, M.E. The Role of the P53 Tumor Suppressor in Metabolism and Diabetes. J. Endocrinol. 2016, 231, R61–R75. [Google Scholar] [CrossRef]
  73. Mehta, Z.B.; Fine, N.; Pullen, T.J.; Cane, M.C.; Hu, M.; Chabosseau, P.; Meur, G.; Velayos-Baeza, A.; Monaco, A.P.; Marselli, L.; et al. Changes in the Expression of the Type 2 Diabetes-Associated Gene VPS13C in the β-Cell Are Associated with Glucose Intolerance in Humans and Mice. Am. J. Physiol. Endocrinol. Metab. 2016, 311, E488–E507. [Google Scholar] [CrossRef] [PubMed]
  74. Dalmas, E.; Toubal, A.; Alzaid, F.; Blazek, K.; Eames, H.L.; Lebozec, K.; Pini, M.; Hainault, I.; Montastier, E.; Denis, R.G.P.; et al. Irf5 Deficiency in Macrophages Promotes Beneficial Adipose Tissue Expansion and Insulin Sensitivity during Obesity. Nat. Med. 2015, 21, 610–618. [Google Scholar] [CrossRef]
  75. Mori, H.; Takahashi, H.; Mine, K.; Higashimoto, K.; Inoue, K.; Kojima, M.; Kuroki, S.; Eguchi, T.; Ono, Y.; Inuzuka, S.; et al. TYK2 Promoter Variant Is Associated with Impaired Insulin Secretion and Lower Insulin Resistance in Japanese Type 2 Diabetes Patients. Genes 2021, 12, 400. [Google Scholar] [CrossRef]
  76. Lee, M.-J.; Feliers, D.; Mariappan, M.M.; Sataranatarajan, K.; Mahimainathan, L.; Musi, N.; Foretz, M.; Viollet, B.; Weinberg, J.M.; Choudhury, G.G.; et al. A Role for AMP-Activated Protein Kinase in Diabetes-Induced Renal Hypertrophy. Am. J. Physiol. Ren. Physiol. 2007, 292, F617–F627. [Google Scholar] [CrossRef]
  77. Sikka, A.; Kaur, M.; Agarwal, C.; Deep, G.; Agarwal, R. Metformin Suppresses Growth of Human Head and Neck Squamous Cell Carcinoma via Global Inhibition of Protein Translation. Cell Cycle 2012, 11, 1374–1382. [Google Scholar] [CrossRef]
  78. Vohra, M.; Sharma, A.R.; Mallya, S.; Prabhu, N.B.; Jayaram, P.; Nagri, S.K.; Umakanth, S.; Rai, P.S. Implications of Genetic Variations, Differential Gene Expression, and Allele-Specific Expression on Metformin Response in Drug-Naïve Type 2 Diabetes. J. Endocrinol. Investig. 2022, 46, 1205–1218. [Google Scholar] [CrossRef] [PubMed]
  79. Zou, Y.; Gong, N.; Cui, Y.; Wang, X.; Cui, A.; Chen, Q.; Jiao, T.; Dong, X.; Yang, H.; Zhang, S.; et al. Forkhead Box P1 (FOXP1) Transcription Factor Regulates Hepatic Glucose Homeostasis. J. Biol. Chem. 2015, 290, 30607–30615. [Google Scholar] [CrossRef]
  80. Bartolini, R.; Medina-Ruiz, L.; Hayes, A.J.; Kelly, C.J.; Halawa, H.A.; Graham, G.J. Inflammatory Chemokine Receptors Support Inflammatory Macrophage and Dendritic Cell Maturation. Immunohorizons 2022, 6, 743–759. [Google Scholar] [CrossRef] [PubMed]
  81. Yang, R.; Wang, M.; Zhang, G.; Li, Y.; Wang, L.; Cui, H. POU2F2 Regulates Glycolytic Reprogramming and Glioblastoma Progression via PDPK1-Dependent Activation of PI3K/AKT/MTOR Pathway. Cell Death Dis. 2021, 12, 433. [Google Scholar] [CrossRef] [PubMed]
  82. Rahman, M.H.; Peng, S.; Hu, X.; Chen, C.; Rahman, M.R.; Uddin, S.; Quinn, J.M.W.; Moni, M.A. A Network-Based Bioinformatics Approach to Identify Molecular Biomarkers for Type 2 Diabetes That Are Linked to the Progression of Neurological Diseases. Int. J. Environ. Res. Public. Health 2020, 17, 1035. [Google Scholar] [CrossRef]
  83. Qian, Z.; Zhao, J. Silencing YY1 Alleviates Ox-LDL-Induced Inflammation and Lipid Accumulation in Macrophages through Regulation of PCSK9/ LDLR Signaling. J. Microbiol. Biotechnol. 2022, 32, 1406–1415. [Google Scholar] [CrossRef]
  84. Mola, S.; Beauchamp, C.; Boucher, G.; Lesage, S.; Karaky, M.; Goyette, P.; Foisy, S.; Rioux, J.D. Identifying Transcript-Level Differential Expression in Primary Human Immune Cells. Mol. Immunol. 2023, 153, 181–193. [Google Scholar] [CrossRef]
  85. Jentho, E.; Ruiz-Moreno, C.; Novakovic, B.; Kourtzelis, I.; Megchelenbrink, W.L.; Martins, R.; Chavakis, T.; Soares, M.P.; Kalafati, L.; Guerra, J.; et al. Trained Innate Immunity, Long-Lasting Epigenetic Modulation, and Skewed Myelopoiesis by Heme. Proc. Natl. Acad. Sci USA. 2021, 118, e2102698118. [Google Scholar] [CrossRef]
  86. Xiong, L.; Chen, L.; Wu, L.; He, W.; Chen, D.; Peng, Z.; Li, J.; Zhu, X.; Su, L.; Li, Y.; et al. Lipotoxicity-Induced CircGlis3 Impairs Beta Cell Function and Is Transmitted by Exosomes to Promote Islet Endothelial Cell Dysfunction. Diabetologia 2022, 65, 188–205. [Google Scholar] [CrossRef]
  87. Oost, L.J.; Slieker, R.C.; Blom, M.T.; Hart, L.M.; Hoenderop, J.G.J.; Beulens, J.W.J.; de Baaij, J.H.F. Genome-Wide Association Study of Serum Magnesium in Type 2 Diabetes. Genes. Nutr. 2024, 19, 2. [Google Scholar] [CrossRef]
  88. Paluvai, H.; Shanmukha, K.D.; Tyedmers, J.; Backs, J. Insights into the Function of HDAC3 and NCoR1/NCoR2 Co-Repressor Complex in Metabolic Diseases. Front. Mol. Biosci. 2023, 10, 1190094. [Google Scholar] [CrossRef] [PubMed]
  89. You, S.; Zhang, Y.; Xu, J.; Qian, H.; Wu, S.; Wu, B.; Lu, S.; Sun, Y.; Zhang, N. The Role of BRG1 in Antioxidant and Redox Signaling. Oxid. Med. Cell Longev. 2020, 2020, 6095673. [Google Scholar] [CrossRef] [PubMed]
  90. Urakami, T. Maturity-Onset Diabetes of the Young (MODY): Current Perspectives on Diagnosis and Treatment. Diabetes Metab. Syndr. Obes. 2019, 12, 1047–1056. [Google Scholar] [CrossRef] [PubMed]
  91. Li, W.; Yang, S.; Chen, G.; He, S. MiR-200c-3p Regulates Pyroptosis by Targeting SLC30A7 in Diabetic Retinopathy. Hum. Exp. Toxicol. 2022, 41, 096032712210995. [Google Scholar] [CrossRef] [PubMed]
  92. Eldakhakhny, B.M.; Al Sadoun, H.; Choudhry, H.; Mobashir, M. In-Silico Study of Immune System Associated Genes in Case of Type-2 Diabetes With Insulin Action and Resistance, and/or Obesity. Front. Endocrinol. 2021, 12, 641888. [Google Scholar] [CrossRef]
  93. Henriksen, T.I.; Wigge, L.V.; Nielsen, J.; Pedersen, B.K.; Sandri, M.; Scheele, C. Dysregulated Autophagy in Muscle Precursor Cells from Humans with Type 2 Diabetes. Sci. Rep. 2019, 9, 8169. [Google Scholar] [CrossRef]
  94. Flores-León, M.; Alcaraz, N.; Pérez-Domínguez, M.; Torres-Arciga, K.; Rebollar-Vega, R.; De la Rosa-Velázquez, I.A.; Arriaga-Canon, C.; Herrera, L.A.; Arias, C.; González-Barrios, R. Transcriptional Profiles Reveal Deregulation of Lipid Metabolism and Inflammatory Pathways in Neurons Exposed to Palmitic Acid. Mol. Neurobiol. 2021, 58, 4639–4651. [Google Scholar] [CrossRef]
  95. You, Y.-A.; Park, S.; Kwon, E.; Kim, Y.-A.; Hur, Y.M.; Lee, G.I.; Kim, S.M.; Song, J.M.; Kim, M.S.; Kim, Y.J.; et al. Maternal PM2.5 Exposure Is Associated with Preterm Birth and Gestational Diabetes Mellitus, and Mitochondrial OXPHOS Dysfunction in Cord Blood. Environ. Sci. Pollut. Res. 2024, 31, 10565–10578. [Google Scholar] [CrossRef]
  96. Atkin, A.S.; Moin, A.S.M.; Nandakumar, M.; Al-Qaissi, A.; Sathyapalan, T.; Atkin, S.L.; Butler, A.E. Impact of Severe Hypoglycemia on the Heat Shock and Related Protein Response. Sci. Rep. 2021, 11, 17057. [Google Scholar] [CrossRef]
  97. Agostino, M.; Rooney, J.; Herat, L.; Matthews, J.; Simonds, A.; Northfield, S.E.; Hopper, D.; Schlaich, M.P.; Matthews, V.B. TNFSF14-Derived Molecules as a Novel Treatment for Obesity and Type 2 Diabetes. Int. J. Mol. Sci. 2021, 22, 10647. [Google Scholar] [CrossRef]
  98. Fu, Y.; Wang, X.; Zhang, L.; Ren, Y.; Hao, L. Allograft Inflammatory Factor-1 Enhances Inflammation and Oxidative Stress via the NF-ΚB Pathway in Diabetic Kidney Disease. Biochem. Biophys. Res. Commun. 2022, 614, 63–69. [Google Scholar] [CrossRef] [PubMed]
  99. Cheng, Y.; Zhang, Z.; Gao, P.; Lai, H.; Zhong, W.; Feng, N.; Yang, Y.; Yu, H.; Zhang, Y.; Han, Y.; et al. AAV Induces Hepatic Necroptosis and Carcinoma in Diabetic and Obese Mice Dependent on Pebp1 Pathway. EMBO Mol. Med. 2023, 15, e17230. [Google Scholar] [CrossRef] [PubMed]
  100. Zhou, Y.; Li, Z.-L.; Ding, L.; Zhang, X.-J.; Liu, N.-C.; Liu, S.-S.; Wang, Y.-F.; Ma, R.-X. Long Noncoding RNA SNHG5 Promotes Podocyte Injury via the MicroRNA-26a-5p/TRPC6 Pathway in Diabetic Nephropathy. J. Biol. Chem. 2022, 298, 102605. [Google Scholar] [CrossRef]
  101. Kazerouni, F.; Bayani, A.; Asadi, F.; Saeidi, L.; Parvizi, N.; Mansoori, Z. Type2 Diabetes Mellitus Prediction Using Data Mining Algorithms Based on the Long-Noncoding RNAs Expression: A Comparison of Four Data Mining Approaches. BMC Bioinform. 2020, 21, 372. [Google Scholar] [CrossRef]
  102. Wang, X.; Li, X.; Zong, Y.; Yu, J.; Chen, Y.; Zhao, M.; Wu, D.; Liao, Y.; Jiang, C.; Zhu, H. Identification and Validation of Genes Related to RNA Methylation Modification in Diabetic Retinopathy. Curr. Eye Res. 2023, 48, 1034–1049. [Google Scholar] [CrossRef]
  103. Degtyareva, A.O.; Antontseva, E.V.; Merkulova, T.I. Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases. Int. J. Mol. Sci. 2021, 22, 6454. [Google Scholar] [CrossRef]
  104. Fabo, T.; Khavari, P. Functional Characterization of Human Genomic Variation Linked to Polygenic Diseases. Trends Genet. 2023, 39, 462–490. [Google Scholar] [CrossRef]
  105. Alsheikh, A.J.; Wollenhaupt, S.; King, E.A.; Reeb, J.; Ghosh, S.; Stolzenburg, L.R.; Tamim, S.; Lazar, J.; Davis, J.W.; Jacob, H.J. The Landscape of GWAS Validation; Systematic Review Identifying 309 Validated Non-Coding Variants across 130 Human Diseases. BMC Med. Genom. 2022, 15, 74. [Google Scholar] [CrossRef] [PubMed]
  106. Abramov, S.; Boytsov, A.; Bykova, D.; Penzar, D.D.; Yevshin, I.; Kolmykov, S.K.; Fridman, M.V.; Favorov, A.V.; Vorontsov, I.E.; Baulin, E.; et al. Landscape of Allele-Specific Transcription Factor Binding in the Human Genome. Nat. Commun. 2021, 12, 2751. [Google Scholar] [CrossRef]
  107. Haider, M.J.A.; Albaqsumi, Z.; Al-Mulla, F.; Ahmad, R.; Al-Rashed, F. SOCS3 Regulates Dectin-2-Induced Inflammation in PBMCs of Diabetic Patients. Cells 2022, 11, 2670. [Google Scholar] [CrossRef]
  108. Ragavi, R.; Adole, P.S.; Vinod, K.V.; Pillai, A.A. Altered Expression of a Disintegrin and Metalloproteinase 10 in Peripheral Blood Mononuclear Cells in Type 2 Diabetes Mellitus Patients with the Acute Coronary Syndrome: A Pilot Study. Endocrine 2022, 77, 461–468. [Google Scholar] [CrossRef] [PubMed]
  109. Pernold, C.P.S.; Lagumdzic, E.; Stadler, M.; Dolezal, M.; Jäckel, S.; Schmitt, M.W.; Mair, K.H.; Saalmüller, A. Species Comparison: Human and Minipig PBMC Reactivity under the Influence of Immunomodulating Compounds in Vitro. Front. Immunol. 2024, 14, 1327776. [Google Scholar] [CrossRef] [PubMed]
  110. Alexovič, M.; Uličná, C.; Sabo, J.; Davalieva, K. Human Peripheral Blood Mononuclear Cells as a Valuable Source of Disease-related Biomarkers: Evidence from Comparative Proteomics Studies. Proteom. Clin. Appl. 2024, 18, 2300072. [Google Scholar] [CrossRef]
  111. Al Dubayee, M.S.; Alayed, H.; Almansour, R.; Alqaoud, N.; Alnamlah, R.; Obeid, D.; Alshahrani, A.; Zahra, M.M.; Nasr, A.; Al-Bawab, A.; et al. Differential Expression of Human Peripheral Mononuclear Cells Phenotype Markers in Type 2 Diabetic Patients and Type 2 Diabetic Patients on Metformin. Front. Endocrinol. 2018, 9, 537. [Google Scholar] [CrossRef]
  112. Fakhoury, H.M.A.; Elahi, M.A.; Al Sarheed, S.; Al Dubayee, M.; Alshahrani, A.; Zhra, M.; Almassri, A.; Aljada, A. Gene Expression Profiling of Peripheral Blood Mononuclear Cells in Type 2 Diabetes: An Exploratory Study. Medicina 2022, 58, 1829. [Google Scholar] [CrossRef]
  113. Eftekharian, M.M.; Karimi, J.; Safe, M.; Sadeghian, A.; Borzooei, S.; Siahpoushi, E. Investigation of the Correlation between Some Immune System and Biochemical Indicators in Patients with Type 2 Diabetes. Hum. Antibodies 2016, 24, 25–31. [Google Scholar] [CrossRef]
  114. Cohen, S.R.; Gardner, T.W. Diabetic Retinopathy and Diabetic Macular Edema. Retin. Pharmacother. 2016, 55, 137–146. [Google Scholar]
  115. Wang, G.; Chiou, J.; Zeng, C.; Miller, M.; Matta, I.; Han, J.Y.; Kadakia, N.; Okino, M.-L.; Beebe, E.; Mallick, M.; et al. Integration of Single-Cell Multiomic Measurements across Disease States with Genetics Identifies Mechanisms of Beta Cell Dysfunction in Type 2 Diabetes. bioRxiv 2023. [Google Scholar] [CrossRef]
  116. Cruz-Pineda, W.D.; Parra-Rojas, I.; Rodríguez-Ruíz, H.A.; Illades-Aguiar, B.; Matia-García, I.; Garibay-Cerdenares, O.L. The Regulatory Role of Insulin in Energy Metabolism and Leukocyte Functions. J. Leukoc. Biol. 2021, 111, 197–208. [Google Scholar] [CrossRef]
  117. Manning, B.D.; Cantley, L.C. AKT/PKB Signaling: Navigating Downstream. Cell 2007, 129, 1261–1274. [Google Scholar] [CrossRef]
  118. Yu, S.; Meng, S.; Xiang, M.; Ma, H. Phosphoenolpyruvate Carboxykinase in Cell Metabolism: Roles and Mechanisms beyond Gluconeogenesis. Mol. Metab. 2021, 53, 101257. [Google Scholar] [CrossRef]
  119. Manoel-Caetano, F.S.; Xavier, D.J.; Evangelista, A.F.; Takahashi, P.; Collares, C.V.; Puthier, D.; Foss-Freitas, M.C.; Foss, M.C.; Donadi, E.A.; Passos, G.A.; et al. Gene Expression Profiles Displayed by Peripheral Blood Mononuclear Cells from Patients with Type 2 Diabetes Mellitus Focusing on Biological Processes Implicated on the Pathogenesis of the Disease. Gene 2012, 511, 151–160. [Google Scholar] [CrossRef] [PubMed]
  120. Gaedigk, A.; Whirl-Carrillo, M.; Pratt, V.M.; Miller, N.A.; Klein, T.E. PharmVar and the Landscape of Pharmacogenetic Resources. Clin. Pharmacol. Ther. 2020, 107, 43–46. [Google Scholar] [CrossRef] [PubMed]
  121. Gharani, N.; Calendo, G.; Kusic, D.; Madzo, J.; Scheinfeldt, L. Star Allele Search: A Pharmacogenetic Annotation Database and User-Friendly Search Tool of Publicly Available 1000 Genomes Project Biospecimens. BMC Genom. 2024, 25, 116. [Google Scholar] [CrossRef]
  122. Rykova, E.; Ershov, N.; Damarov, I.; Merkulova, T. SNPs in 3′UTR MiRNA Target Sequences Associated with Individual Drug Susceptibility. Int. J. Mol. Sci. 2022, 23, 13725. [Google Scholar] [CrossRef]
  123. Li, H.; Chen, L.; Tong, X.; Dai, H.; Shi, T.; Cheng, X.; Sun, M.; Chen, K.; Wei, Q.; Wang, M. Functional Genetic Variants of CTNNBIP1 Predict Platinum Treatment Response of Chinese Epithelial Ovarian Cancer Patients. J. Cancer 2020, 11, 6850–6860. [Google Scholar] [CrossRef]
  124. Yang, J.; Li, S.; Lv, H.; Wang, W.; Zhang, J.; Chu, L.; Zhang, Y. CREB1 and BDNF Gene Polymorphisms Are Associated with Early Treatment Response to Escitalopram in Panic Disorder. J. Affect. Disord. 2021, 278, 536–541. [Google Scholar] [CrossRef] [PubMed]
  125. Hsieh, T.-Y.; Lin, J.-F.; Liu, F.-C.; Chen, H.-C.; Lui, S.-W.; Chang, Y.-T. Functional Implications of Rs9373441 with FOXP3+Treg and Tr1 for the Clinical Effectiveness of CsDMARDs in Rheumatoid Arthritis. Clin. Chim. Acta 2023, 551, 117612. [Google Scholar] [CrossRef]
  126. Cook, M.N.; Girman, C.J.; Stein, P.P.; Alexander, C.M. Initial Monotherapy with Either Metformin or Sulphonylureas Often Fails to Achieve or Maintain Current Glycaemic Goals in Patients with Type 2 Diabetes in UK Primary Care. Diabet. Med. 2007, 24, 350–358. [Google Scholar] [CrossRef]
  127. Xiao, D.; Guo, Y.; Li, X.; Yin, J.-Y.; Zheng, W.; Qiu, X.-W.; Xiao, L.; Liu, R.-R.; Wang, S.-Y.; Gong, W.-J.; et al. The Impacts of SLC22A1 Rs594709 and SLC47A1 Rs2289669 Polymorphisms on Metformin Therapeutic Efficacy in Chinese Type 2 Diabetes Patients. Int. J. Endocrinol. 2016, 2016, 4350712. [Google Scholar] [CrossRef]
  128. Goswami, S.; Yee, S.W.; Stocker, S.; Mosley, J.D.; Kubo, M.; Castro, R.; Mefford, J.A.; Wen, C.; Liang, X.; Witte, J.; et al. Genetic Variants in Transcription Factors Are Associated With the Pharmacokinetics and Pharmacodynamics of Metformin. Clin. Pharmacol. Ther. 2014, 96, 370–379. [Google Scholar] [CrossRef] [PubMed]
  129. Zhou, Y.; Ye, W.; Wang, Y.; Jiang, Z.; Meng, X.; Xiao, Q.; Zhao, Q.; Yan, J. Genetic Variants of OCT1 Influence Glycemic Response to Metformin in Han Chinese Patients with Type-2 Diabetes Mellitus in Shanghai. Int. J. Clin. Exp. Pathol. 2015, 8, 9533–9542. [Google Scholar] [PubMed]
  130. Nasykhova, Y.A.; Tonyan, Z.N.; Mikhailova, A.A.; Danilova, M.M.; Glotov, A.S. Pharmacogenetics of Type 2 Diabetes—Progress and Prospects. Int. J. Mol. Sci. 2020, 21, 6842. [Google Scholar] [CrossRef] [PubMed]
  131. Nasykhova, Y.; Barbitoff, Y.; Tonyan, Z.; Danilova, M.; Nevzorov, I.; Komandresova, T.; Mikhailova, A.; Vasilieva, T.; Glavnova, O.; Yarmolinskaya, M.; et al. Genetic and Phenotypic Factors Affecting Glycemic Response to Metformin Therapy in Patients with Type 2 Diabetes Mellitus. Genes 2022, 13, 1310. [Google Scholar] [CrossRef]
  132. Jablonski, K.A.; McAteer, J.B.; de Bakker, P.I.W.; Franks, P.W.; Pollin, T.I.; Hanson, R.L.; Saxena, R.; Fowler, S.; Shuldiner, A.R.; Knowler, W.C.; et al. Common Variants in 40 Genes Assessed for Diabetes Incidence and Response to Metformin and Lifestyle Intervention in the Diabetes Prevention Program. Diabetes 2010, 59, 2672–2681. [Google Scholar] [CrossRef]
  133. Rotroff, D.M.; Yee, S.W.; Zhou, K.; Marvel, S.W.; Shah, H.S.; Jack, J.R.; Havener, T.M.; Hedderson, M.M.; Kubo, M.; Herman, M.A.; et al. Genetic Variants in CPA6 and PRPF31 Are Associated With Variation in Response to Metformin in Individuals With Type 2 Diabetes. Diabetes 2018, 67, 1428–1440. [Google Scholar] [CrossRef]
  134. Oni-Orisan, A.; Haldar, T.; Ranatunga, D.K.; Medina, M.W.; Schaefer, C.; Krauss, R.M.; Iribarren, C.; Risch, N.; Hoffmann, T.J. The Impact of Adjusting for Baseline in Pharmacogenomic Genome-Wide Association Studies of Quantitative Change. NPJ Genom. Med. 2020, 5, 1. [Google Scholar] [CrossRef]
  135. Agúndez, J.A.G.; Ayuso, P.; Quiñones, L.A.; García-Martín, E. Editorial: Editor’s Feature: Negative Findings in Pharmacogenetics and Pharmacogenomics. Front. Pharmacol. 2023, 14, 1267344. [Google Scholar] [CrossRef]
  136. Meliala, I.T.S.; Hosea, R.; Kasim, V.; Wu, S. The Biological Implications of Yin Yang 1 in the Hallmarks of Cancer. Theranostics 2020, 10, 4183–4200. [Google Scholar] [CrossRef]
  137. Verheul, T.C.J.; van Hijfte, L.; Perenthaler, E.; Barakat, T.S. The Why of YY1: Mechanisms of Transcriptional Regulation by Yin Yang 1. Front. Cell Dev. Biol. 2020, 8, 592164. [Google Scholar] [CrossRef]
  138. Xiao, R.; Chen, J.-Y.; Liang, Z.; Luo, D.; Chen, G.; Lu, Z.J.; Chen, Y.; Zhou, B.; Li, H.; Du, X.; et al. Pervasive Chromatin-RNA Binding Protein Interactions Enable RNA-Based Regulation of Transcription. Cell 2019, 178, 107–121.E18. [Google Scholar] [CrossRef]
  139. Wang, J.; Wu, X.; Wei, C.; Huang, X.; Ma, Q.; Huang, X.; Faiola, F.; Guallar, D.; Fidalgo, M.; Huang, T.; et al. YY1 Positively Regulates Transcription by Targeting Promoters and Super-Enhancers through the BAF Complex in Embryonic Stem Cells. Stem Cell Rep. 2018, 10, 1324–1339. [Google Scholar] [CrossRef]
  140. Bianchi, M.; Crinelli, R.; Giacomini, E.; Carloni, E.; Radici, L.; Magnani, M. Yin Yang 1 Intronic Binding Sequences and Splicing Elicit Intron-Mediated Enhancement of Ubiquitin C Gene Expression. PLoS ONE 2013, 8, e65932. [Google Scholar] [CrossRef] [PubMed]
  141. Shi, Y.; Lee, J.-S.; Galvin, K.M. Everything You Have Ever Wanted to Know about Yin Yang 1. Biochim. Biophys. Acta (BBA) Rev. Cancer 1997, 1332, F49–F66. [Google Scholar] [CrossRef]
  142. Liu, D.; Yang, K.Y.; Chan, V.W.; Ye, W.; Chong, C.C.N.; Wang, C.C.; Wang, H.; Zhou, B.; Cheng, K.K.Y.; Lui, K.O. YY1 Regulates Glucose Homeostasis Through Controlling Insulin Transcription in Pancreatic β-Cells. Diabetes 2022, 71, 961–977. [Google Scholar] [CrossRef]
  143. Martins Peçanha, F.L.; Jaafar, R.; Werneck-de-Castro, J.P.; Apostolopolou, C.-C.; Bhushan, A.; Bernal-Mizrachi, E. The Transcription Factor YY1 Is Essential for Normal DNA Repair and Cell Cycle in Human and Mouse β-Cells. Diabetes 2022, 71, 1694–1705. [Google Scholar] [CrossRef]
  144. Lu, J.; Jin, K.; Jiao, J.; Liu, R.; Mou, T.; Chen, B.; Zhang, Z.; Jiang, C.; Zhao, H.; Wang, Z.; et al. YY1 (Yin-Yang 1), a Transcription Factor Regulating Systemic Inflammation, Is Involved in Cognitive Impairment of Depression. Psychiatry Clin. Neurosci. 2023, 77, 149–159. [Google Scholar] [CrossRef]
  145. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  146. Genome Assembly GRCh38. Available online: https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_000001405.26/ (accessed on 25 March 2023).
  147. Langmead, B.; Salzberg, S.L. Fast Gapped-Read Alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef]
  148. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve Years of SAMtools and BCFtools. Gigascience 2021, 10, giab008. [Google Scholar] [CrossRef]
  149. Sherry, S.T.; Ward, M.H.; Kholodov, M.; Baker, J.; Phan, L.; Smigielski, E.M.; Sirotkin, K. DbSNP: The NCBI Database of Genetic Variation. Nucleic Acids Res. 2001, 29, 308–311. [Google Scholar] [CrossRef] [PubMed]
  150. SNPlocs.Hsapiens.dbSNP155.GRCh38. Human SNP Locations and Alleles Extracted from dbSNP Build 155 and Placed on the GRCh38/hg38 Assembly. Available online: https://www.bioconductor.org/packages/release/data/annotation/html/SNPlocs.Hsapiens.dbSNP155.GRCh38.html (accessed on 25 March 2023).
  151. Bioconductor Core Team; Bioconductor Package Maintainer. TxDb.Hsapiens.UCSC.Hg38.KnownGene: Annotation Package for TxDb Object(s); R Package Version 3.4.6; Bioconductor Core Team: Winnipeg, MB, Canada, 2019. [Google Scholar]
  152. Mayba, O.; Gilbert, H.N.; Liu, J.; Haverty, P.M.; Jhunjhunwala, S.; Jiang, Z.; Watanabe, C.; Zhang, Z. MBASED: Allele-Specific Expression Detection in Cancer Tissues and Cell Lines. Genome Biol. 2014, 15, 405. [Google Scholar] [CrossRef] [PubMed]
  153. Zhang, J.; Zhao, H. EQTL Studies: From Bulk Tissues to Single Cells. arXiv 2023, arXiv:2302.11662. [Google Scholar] [CrossRef]
  154. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  155. Gene Expression Omnibus. Available online: https://www.ncbi.nlm.nih.gov/geo/ (accessed on 14 January 2024).
  156. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for Functional Genomics Data Sets—Update. Nucleic Acids Res. 2012, 41, D991–D995. [Google Scholar] [CrossRef]
  157. CytoNCA. Available online: https://apps.cytoscape.org/apps/cytonca (accessed on 14 January 2024).
  158. Tang, Y.; Li, M.; Wang, J.; Pan, Y.; Wu, F.-X. CytoNCA: A Cytoscape Plugin for Centrality Analysis and Evaluation of Protein Interaction Networks. Biosystems 2015, 127, 67–72. [Google Scholar] [CrossRef]
  159. Cytoscape. Available online: https://cytoscape.org/ (accessed on 14 January 2024).
  160. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  161. MCODE. Available online: https://apps.cytoscape.org/apps/mcode (accessed on 25 March 2023).
  162. Bader, G.D.; Hogue, C.W. An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinform. 2003, 4, 2. [Google Scholar] [CrossRef]
  163. Wu, T.; Hu, E.; Xu, S.; Chen, M.; Guo, P.; Dai, Z.; Feng, T.; Zhou, L.; Tang, W.; Zhan, L.; et al. ClusterProfiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innov. 2021, 2, 100141. [Google Scholar] [CrossRef]
  164. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.-C.; Müller, M. PROC: An Open-Source Package for R and S+ to Analyze and Compare ROC Curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef]
Figure 1. Scheme of the main stages in the search for rSNPs and their further analysis. The solid gray frame shows the stages of the search for rSNPs; the dotted gray frame, further analysis of the rSNP panel; and the cyan italic, data sources.
Figure 1. Scheme of the main stages in the search for rSNPs and their further analysis. The solid gray frame shows the stages of the search for rSNPs; the dotted gray frame, further analysis of the rSNP panel; and the cyan italic, data sources.
Ijms 25 09297 g001
Figure 2. Venn diagram showing the number of the rSNPs localized to allele-specific transcription factor binding sites (ANANASTRA), rSNPs associated phenotypic traits (GWAS data), and eQTLs effects (GTEx data).
Figure 2. Venn diagram showing the number of the rSNPs localized to allele-specific transcription factor binding sites (ANANASTRA), rSNPs associated phenotypic traits (GWAS data), and eQTLs effects (GTEx data).
Ijms 25 09297 g002
Figure 3. Volcano plot of DEGs. The horizontal axis stands for log2 fold change and the vertical axis, for –log10 (adjusted p-value). Statistically significant DEGs harboring rSNPs in their promoters are marked red.
Figure 3. Volcano plot of DEGs. The horizontal axis stands for log2 fold change and the vertical axis, for –log10 (adjusted p-value). Statistically significant DEGs harboring rSNPs in their promoters are marked red.
Ijms 25 09297 g003
Figure 4. Network chart illustrating the link of top ten significant KEGG pathways according to the enrichment analysis with upregulated hub genes.
Figure 4. Network chart illustrating the link of top ten significant KEGG pathways according to the enrichment analysis with upregulated hub genes.
Ijms 25 09297 g004
Figure 5. PPI network in the significant modules analyzed for KEGG and GO enrichment. (A) First module. (B) Second module. (C) Third module. In the PPI network, nodes show proteins and edges, their interaction. Hub proteins are denoted with larger symbols and the KEGG- and GO-annotated proteins, with the corresponding color (see the legend).
Figure 5. PPI network in the significant modules analyzed for KEGG and GO enrichment. (A) First module. (B) Second module. (C) Third module. In the PPI network, nodes show proteins and edges, their interaction. Hub proteins are denoted with larger symbols and the KEGG- and GO-annotated proteins, with the corresponding color (see the legend).
Ijms 25 09297 g005
Figure 6. Graphical visualization of DEG representation in (A) insulin and (B) PI3K/Akt signaling pathways (KEGG data). Colors of nodes show the direction and value (log2FC) of expression alteration.
Figure 6. Graphical visualization of DEG representation in (A) insulin and (B) PI3K/Akt signaling pathways (KEGG data). Colors of nodes show the direction and value (log2FC) of expression alteration.
Ijms 25 09297 g006aIjms 25 09297 g006b
Figure 7. Volcano plot of DEGs. The horizontal axis shows the log2 fold change and the vertical axis, –log10 (adjusted p-value). Significant DEGs with detected rSNPs are colored red.
Figure 7. Volcano plot of DEGs. The horizontal axis shows the log2 fold change and the vertical axis, –log10 (adjusted p-value). Significant DEGs with detected rSNPs are colored red.
Ijms 25 09297 g007
Table 1. Top 20 traits ordered according to the number of GWAS-derived rSNPs.
Table 1. Top 20 traits ordered according to the number of GWAS-derived rSNPs.
TraitTotal Number of SNPs per Trait in GWAS CatalogNumber of rSNPs Directly Included in GWAS Catalogp.adj.Odds Ratio *
Waist-to-hip ratio adjusted for BMI3741361.1 × 10−32.2 (1.5–3.1)
Mean corpuscular volume2227311.1 × 10−53.2 (2.2–4.6)
White blood cell count2456316.1 × 10−52.9 (2–4.2)
Waist–hip index2526281.1 × 10−32.5 (1.7–3.7)
Hip circumference adjusted for BMI3359284.1 × 10−21.9 (1.3–2.8)
Platelet count2609273.9 × 10−32.4 (1.6–3.5)
Mean corpuscular hemoglobin2365254.8 × 10−32.4 (1.6–3.6)
Lymphocyte count1651249.8 × 10−53.3 (2.1–5)
Monocyte count1767242.1 × 10−43.1 (2–4.7)
Red blood cell count2487241.9 × 10−22.2 (1.4–3.3)
Neutrophil count1559231.1 × 10−43.4 (2.1–5.1)
Eosinophil count2100234.8 × 10−32.5 (1.6–3.8)
Type 2 diabetes2830237.1 × 10−21.9 (1.2–2.8)
Red cell distribution width1732205.5 × 10−32.6 (1.6–4.1)
Mean platelet volume1402162.1 × 10−22.6 (1.5–4.3)
A body shape index1514163.6 × 10−22.4 (1.4–3.9)
Monocyte percentage of white cells738151.9 × 10−44.6 (2.6–7.7)
Plateletcrit930151.6 × 10−33.7 (2–6.1)
Mean spheric corpuscular volume821141.6 × 10−33.9 (2.1–6.6)
Appendicular lean mass1569148.8 × 10−22 (1.1–3.4)
* 95% confidence interval is parenthesized.
Table 2. rSNPs in promoters of the hub genes with AUC > 0.7.
Table 2. rSNPs in promoters of the hub genes with AUC > 0.7.
GeneRegulationAUCRs_id
NOTCH1Up0.7212rs951509664, rs3013307, and rs3013306
H6PDUp0.8072rs184437520, rs3752547, rs9435144, and rs11121354
POLR2AUp0.7314rs4796424, rs57985740, rs41555218, rs144575559, and rs9901161
NCOR2Up0.7713rs7960906, rs1006100, rs1199426444, rs191752208, rs1432659465, rs998518300, rs906886068, rs1458070990, rs948418315, rs79830634, rs12426514, rs1316249, rs924583078, rs868110059, and rs1407929149
PXNUp0.7181rs7953949 and rs3890165
FASNUp0.7642rs7209621 and rs62078751
SCARB1Up0.7323rs7305310, rs838884, and rs897715
GAKUp0.734rs140032537, rs141564663, rs1403319282, rs182955420, rs3775124 and rs3733352
CTSDUp0.7238rs2292963, rs144932926, rs2292962, and rs35640004
FZR1Up0.7145rs8100223 and rs8644
SMG1Up0.7926rs142606705, rs12929094, and rs560580650
TP53Up0.7522rs1800899
MAN1B1Up0.7762rs4880199 and rs10870178
RPL3Down0.766rs5757613, rs2072872, rs137626, rs2076125, rs143897309, rs969895370, rs84491, rs137627, rs470081, rs754570306, rs6509, rs137620, and rs12484030
RPS3Down0.7119rs186612441
HSP90AB1Down0.8227rs324131
RPL11Down0.801rs3753270, rs111953674, rs878908315, rs558662093, and rs1361739260
RPS11Down0.7358rs739349
RPL5Down0.7411rs34244251
RPS20Down0.7159rs17814456
RPL13ADown0.7101rs11539123
RACK1Down0.8054rs2287715 and rs111326428
RPL14Down0.7908rs62263890 and rs2276869
RPL12Down0.8285rs2247310 and rs2247322
CCT7Down0.8027rs779122697
RUVBL1Down0.7451rs11719546
TCP1Down0.7292rs62621403
PSMA7Down0.8005rs73307256, rs6089665, and rs3746651
PSMC5Down0.7863rs141975038
ATP5F1ADown0.7092rs3753069, rs41274316, and rs34907121
HNRNPKDown0.742rs1011582290, rs796004, and rs296890
SKP1Down0.7163rs56257643
BCL2Down0.7176rs4987834 and rs1160274961
FYNDown0.7438rs62413757, rs189480003, rs71564104, rs9487736, rs1295051482, rs6568706, rs17072881, and rs1057979
SSRP1Down0.8262rs61888886 and rs61888888
CDC42Down0.7247rs61778042, rs866829764, rs11801382, rs2255282, rs12038474, and rs16826302
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Damarov, I.S.; Korbolina, E.E.; Rykova, E.Y.; Merkulova, T.I. Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response. Int. J. Mol. Sci. 2024, 25, 9297. https://doi.org/10.3390/ijms25179297

AMA Style

Damarov IS, Korbolina EE, Rykova EY, Merkulova TI. Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response. International Journal of Molecular Sciences. 2024; 25(17):9297. https://doi.org/10.3390/ijms25179297

Chicago/Turabian Style

Damarov, Igor S., Elena E. Korbolina, Elena Y. Rykova, and Tatiana I. Merkulova. 2024. "Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response" International Journal of Molecular Sciences 25, no. 17: 9297. https://doi.org/10.3390/ijms25179297

APA Style

Damarov, I. S., Korbolina, E. E., Rykova, E. Y., & Merkulova, T. I. (2024). Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response. International Journal of Molecular Sciences, 25(17), 9297. https://doi.org/10.3390/ijms25179297

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