ORIGINAL RESEARCH article

Front. Endocrinol., 01 July 2022

Sec. Diabetes: Molecular Mechanisms

Volume 13 - 2022 | https://doi.org/10.3389/fendo.2022.898893

Rs864745 in JAZF1, an Islet Function Associated Variant, Correlates With Plasma Lipid Levels in Both Type 1 and Type 2 Diabetes Status, but Not Healthy Subjects

  • Department of Endocrinology and Metabolism, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China

Article metrics

View details

4

Citations

2k

Views

967

Downloads

Abstract

Objective:

This study aims to reveal the association between JAZF1 rs864745 A>G variant and type 2 diabetes (T2D), type 1 diabetes (T1D) risk, and their correlation with clinical features, including islet function, islet autoimmunity, and plasma lipid levels.

Methods:

We included 2505 healthy controls based on oral glucose tolerance test (OGTT), 1736 unrelated T2D, and 1003 unrelated autoantibody-positive T1D individuals. Binary logistic regression was performed to evaluate the relationships between rs864745 in JAZF1 and T2D, T1D, and islet-specific autoantibody status under the additive model, while multiple linear regression was used to assess its effect on glycemic-related quantitative traits and plasma lipid levels.

Results:

We did not find any association between rs864745 in JAZF1 and T2D, T1D, or their subgroups (All P > 0.05). For glycemic traits, we found that the G allele of this variant was significantly associated with higher 120 min insulin level, insulinogenic index (IGI), corrected insulin response (CIR), and acute insulin response (BIGTT-AIR) (P = 0.033, 0.006, 0.009, and 0.016, respectively) in healthy individuals. Similar associations were observed in newly diagnosed T2D but not T1D individuals. Although this variant had no impact on islet autoimmunity (All P > 0.05), significant associations with plasma total cholesterol (TC) and low-density lipoprotein (LDL) level stratified by JAZF1 rs864745 variant were observed in the disease status of T2D (P = 0.002 and 0.003) and T1D (P = 0.024 and 0.009), with significant heterogeneity to healthy individuals.

Conclusions:

The common JAZF1 rs864745 variant contributes to islet function and lipid metabolism, which might be put into genetic risk scores to assess the risk of related clinical features.

Introduction

JAZF Zinc Finger 1 (JAZF1), a multifunctional regulatory factor, was initially identified as an orphan nuclear receptor corepressor and involved in gluconeogenesis and lipid metabolism. The human JAZF1 gene, located on chromosome 7p15.2, is predominantly expressed in insulin-responsive organs such as the liver, fat, skeletal muscle, and pancreas (1, 2). Furthermore, JAZF1 is expressed in both human (3, 4) and mouse β-cells (5, 6). JAZF1 is an important regulator of endoplasmic reticulum (ER) stress and ribosome biogenesis via a feedback action preventing the activation of ER and p53 stress-mediated β-cell apoptosis (7), which plays a pivotal role in the differentiation of β-cells and glucose homeostasis (8). Obesity (9) and diabetes status (10, 11) alter the JAZF1 expression pattern.

Genome-wide association studies (GWAS) have shown that the risk A allele of rs864745, an intron variant in JAZF1, correlates with both type 1 (T1D) (12) in Europeans, and type 2 diabetes (T2D) (13) in various ethnicities, including Europeans, Americans, as well as Asians (11, 14–17). Moreover, this variant also correlated with β-cell function (14, 18).

Studies have demonstrated that T1D and T2D have heterogeneous characteristics in terms of both genetic (13, 19) and phenotypic features (20, 21). Previous studies have not fully elucidated the underlying contribution of this variant on T1D and T2D subtypes and related clinical phenotypes. Therefore, this study aimed to shed light on the relationship between rs864745 in JAZF1 and T1D, T2D subgroups, and also β-cell function and lipid metabolism in both healthy and diabetes disease status.

Materials and Methods

Study Population

This study includes unrelated 1736 T2D, 1003 T1D, and 2505 non-diabetic controls. Their clinical characteristics are shown in Table 1. All the unrelated healthy controls and T2D subjects were recruited from Nanjing, one of the 25 communities in the REACTION study (22), and unrelated T1D participants were recruited from the First Affiliated Hospital of Nanjing Medical University between January 2008 and December 2020. T2D and T1D were diagnosed according to the World Health Organization criteria, and diabetes caused by liver dysfunction, pancreatitis, or other autoimmune diseases was excluded. Furthermore, T1D cases were identified as only diabetic patients with at least one autoantibody-positive [zinc transporter 8 antibody (ZnT8A), glutamate decarboxylase antibody (GADA), or insulinoma-associated-2 antibody (IA-2A)]. Non-diabetic controls were extracted based on oral glucose tolerance test (OGTT), and the inclusion criteria for healthy controls were as follows: HbA1c ≤ 6.0% (4.2 mmol/mol), fasting plasma glucose < 6.1 mmol/L, and two-hour plasma glucose < 7.8 mmol/L. The study population was determined as Chinese Han by questionnaire. All samples were collected with appropriate informed consent from all participants and/or their guardians in a written way. This study was approved by the Ethics Committee from the First Affiliated Hospital of Nanjing Medical University and conducted according to the Declaration of Helsinki II principles.

Table 1

T2DT1DControl
n173610032505
Sex (male/female)750/986491/512786/1719
Age (years)62.0 ± 8.023.7 ± 14.854.0 ± 9.0
BMI (kg/m2)25.20 ± 3.4719.84 ± 3.7823.65 ± 3.00
Lipids medication (%)3.1%–1.0%
Plasma glucose (mmol/l)
Fasting7.66 ± 2.27–5.33 ± 0.36
30 min post OGTT11.26 ± 3.06–8.58 ± 1.48
120 min post OGTT12.79 ± 4.17–6.14 ± 0.98
Serum insulin (mIU/L)
Fasting11.32 (8.29, 16.25)–9.86 (7.38, 13.38)
30 min post OGTT23.62 (14.51, 41.79)–57.13 (37.80, 88.59)
120 min post OGTT43.07 (22.74, 83.70)–40.72 (26.25, 63.37)
HbA1c (%)7.10 ± 1.5011.15 ± 3.065.60 ± 0.40
HDL (mmol/L)1.23 ± 0.301.36 ± 0.611.37 ± 0.34
LDL (mmol/L)2.83 ± 0.812.69 ± 1.032.81 ± 0.74
TC (mmol/L)4.86 ± 1.064.31 ± 1.844.88 ± 0.96
TG (mmol/L)1.82 ± 1.321.52 ± 1.911.37 ± 0.85
Serum antibody positive rate
ZnT8A–0.45–
GADA–0.74–
IA-2A–0.50–

Overall study population and metabolic index characteristics.

Data are expressed as mean ± SD or as median (interquartile range).

T2D, type 2 diabetes; T1D, type 1 diabetes; BMI, body mass index; OGTT, oral glucose tolerance test; HbA1c, glycated hemoglobin; HDL, high-density lipoprotein; LDL, low-density lipoprotein; TC, total cholesterol; TG, triglyceride; ZnT8A, zinc transporter 8 antibody; GADA, glutamate decarboxylase antibody; IA-2A, insulinoma-associated-2 antibody.

Laboratory Measurements for Glycemic and Lipid Traits and Islet-Specific Autoantibodies

All non-diabetic controls and T2D individuals measured plasma glucose and serum insulin levels at fasting, 30, and 120 min after OGTT. Serum C-peptide levels were measured in newly diagnosed T1D patients (disease duration less than 3 months) on fasting, 30, 60, 120, and 180 min standard mixed meal tolerance test (MMTT). Serum insulin levels were measured by an insulin radioimmunoassay kit (BNIBT, China). Serum C-peptide levels and plasma lipid levels, including high-density lipoprotein cholesterol (HDL), low-density lipoprotein cholesterol (LDL), total cholesterol (TC), and triglyceride (TG), were determined by the chemiluminescence method (Roche Diagnostics, Switzerland). The methods for detecting islet-specific autoantibodies were described in our previous studies (19).

Genotyping Assay

Genomic DNA was extracted from peripheral blood lymphocytes using the QIAamp DNA Blood Extraction Kit (QIAGEN, Germany). Genotyping of the JAZF1 rs864745 variant was performed using the Sequenom Massarray (BGI CO.LTD, China). The primer sequences and UEP primers were as follows: Forward, 5-ACGTTGGATGCATTGAACATTTCCTACAACC-3; Reverse, 5-ACGTTGGATGCCATATAAGTGATGCTCAAA-3; UEP, 5-ATGCTCAAATATAATTTGAACTGTTA-3. The genotyping success rate was above 98%, and 100 randomly selected samples were repeated for complete agreement. The genotype distribution in control, T2D, and T1D groups were consistent with Hardy-Weinberg equilibrium (P > 0.05).

Statistical Analyses

We measured and calculated insulin release and sensitivity derived from the OGTT in non-diabetic controls and T2D individuals according to the formulas described previously, including homeostasis model assessment of β-cell function (HOMA-B) and insulin resistance (HOMA-IR), insulinogenic index (IGI), acute insulin response index (BIGTT-AIR), insulin sensitivity index (BIGTT-SI), corrected insulin response (CIR), and Matsuda’s insulin sensitivity index (Matsuda ISI) (23). Serum insulin, C-peptide levels, insulin release and resistance index were log-transformed, and then statistically analyzed. Under an appropriately adjusted additive model, a binary logistic regression analysis was performed for the relationship between rs864745 in JAZF1 and T2D and T1D subgroups, and islet-specific autoantibody status. The effect of this variant on glycemic quantitative traits and serum lipid levels were analyzed by multiple linear regression. The glycemic quantitative traits were adjusted by sex, age, and BMI except that BIGTT-AIR and BIGTT-SI were corrected for age. Serum lipid levels were also adjusted by lipids medication in addition to age, gender, and BMI. All P-values were two-sided, and P < 0.05 was considered significant. The heterogeneity was considered significant with I2 > 75% and P < 0.05. Statistical analyses were performed using SPSS 18.0 and STATA 11.0.

Results

JAZF1 rs864745 A>G Variant Is Not Associated With the Risk of Either T2D or T1D Subgroups

As shown in Table 2, we did not find any association between JAZF1 rs864745 A>G variant and T2D or T1D risk in total (P > 0.05). Although obesity is a critical risk factor for T2D, no association with this variant was observed in lean or obese T2D subgroups (P > 0.05), stratified by BMI status according to the Chinese criteria (24). As age at diagnosis is an important marker for classifying different subgroups of T1D (21), we divided the T1D population into two subgroups according to age at diagnosis ≥ 18 and < 18 years old, and the results revealed no significant association in either children/adolescents or adults T1D subgroups (P > 0.05).

Table 2

GroupsGenotype distributionORadj (95% CI)P
AAAGGG
T2D
Total1020619970.972 (0.870, 1.086)0.616
BMI < 24kg/m2372228350.958 (0.808, 1.136)0.620
24 ≤ BMI < 28kg/m2464296460.978 (0.845, 1.130)0.760
BMI ≥ 28kg/m218495160.838 (0.620, 1.133)0.250
T1D
Total589360440.930 (0.739, 1.171)0.538
Diagnosed age ≥ 18358196180.922 (0.730, 1.166)0.498
Diagnosed age < 18228157240.613 (0.262, 1.435)0.259
Control1454896155

The association between the JAZF1 rs864745 A>G variant and T1D, T2D risk.

A/G: risk/non-risk alleles for T1D, T2D (ref 11–18). G allele is the effect allele for Oradj. T2D: Binary logistic regression, dependent variable: whether diabetic or not, independent variables: genotype, sex, age. BMI normal < 24kg/m2, overweight ≥ 24kg/m2, < 28kg/m2, obese ≥ 28kg/m2. T1D: Binary logistic regression, dependent variable: whether diabetic or not, correction: sex, age (non-onset age). A total of 12 T1D patients have unknown diagnosed age. T2D, type 2 diabetes; T1D, type 1 diabetes; BMI, body mass index.

JAZF1 rs864745 A>G Variant Significantly Correlates With OGTT-Related Insulin Release, but Not Insulin Sensitivity Indices in Healthy Individuals

We analyzed OGTT-derived indicators of glycemic-related traits in healthy individuals. As shown in Table 3, this variant did not affect fasting, 30 min, or 120 min plasma glucose levels. However, we found that the G allele of this variant significantly improved islet function associated with 120 min insulin level, IGI, CIR, and BIGTT-AIR (P = 0.033, 0.006, 0.009, and 0.016, respectively). We further stratified them by BMI status, the G allele of this variant was significantly associated with high IGI, BIGTT-AIR (P = 0.01 and 0.037, respectively) in normal-weight (BMI < 24kg/m2) individuals, as shown in Table S1. However, we did not find any association of this variant with insulin release or insulin sensitivity in obese or overweight individuals. Subsequently, we also assessed the impact of this variant on the islet function of newly diagnosed T2D and T1D individuals. This variant was associated with 120 min insulin level in newly diagnosed T2D individuals (p = 0.041), as shown in Table S2. However, this variant did not affect fasting or responsive C-peptide levels in newly diagnosed T1D individuals, as shown in Table S3. We also evaluated the association with islet-specific autoantibody status, but we did not find any relationship to the positive rate of ZnT8A (p = 0.723), GADA (p = 0.101), or IA-2A (p = 0.316), as shown in Table 4.

Table 3

ControlAAAGGGβSEPadj
Plasma glucose (mmol/l)
Fasting5.33 (5.10, 5.58)5.32 (5.05, 5.58)5.34 (4.95, 5.62)-0.0290.0110.144
30 min post OGTT8.57 (7.65, 9.59)8.43 (7.58, 9.46)8.43 (7.49, 9.30)-0.0370.0470.057
120 min post OGTT6.21 (5.49, 6.88)6.16 (5.56, 6.85)6.38 (5.58, 7.00)0.0150.0310.433
Serum insulin (mIU/L)
Fasting9.83 (7.36, 13.22)9.93 (7.45, 13.65)10.27 (7.27, 13.64)0.0210.0070.273
30 min post OGTT56.14 (37.34, 85.83)57.49 (37.44, 92.70)61.30 (43.99, 85.60)0.0340.0090.078
120 min post OGTT39.52 (26.12, 61.27)41.93 (26.28, 63.88)43.81 (27.66, 73.57)0.0410.0100.033*
Islet function
HOMA-B105.72 (80.78, 151.27)110.29 (81.39, 154.66)117.08 (80.98, 159.25)0.0360.0080.072
IGI14.77 (8.64, 24.91)16.09 (9.06, 28.59)16.67 (8.70, 27.75)0.0540.0130.006**
BIGTT-AIR7.58 (7.29, 7.94)7.61 (7.34, 8.04)7.68 (7.37, 8.00)0.0480.0090.016*
CIR149.66 (90.35, 238.67)163.72 (96.68, 274.66)168.27 (94.78, 270.96)0.0510.0110.009**
Insulin resistance
HOMA-IR2.31 (1.71, 3.12)2.35 (1.75, 3.26)2.45 (1.67, 3.35)0.0160.0070.423
Matsuda ISI0.07 (0.04, 0.12)0.07 (0.04, 0.11)0.06 (0.04, 0.11)-0.0330.0100.080
BIGTT-SI1.94 (1.54, 2.24)1.92 (1.50, 2.25)1.88 (1.48, 2.16)-0.0240.0090.224
HbA1c%5.60 (5.40, 5.90)5.70 (5.40, 5.90)5.60 (5.40, 5.90)0.0040.0120.826

The association of the JAZF1 rs864745 A>G variant with glycemic quantitative traits in healthy individuals based on OGTT.

Linear regression, with several correction indices calculated separately, with genotype and correction indices as independent variables and islet function as dependent variable. *P < 0.05; **P < 0.01.

OGTT, oral glucose tolerance test; HOMA-B, homeostasis model assessment of β-cell function; IGI, insulinogenic index; BIGTT-AIR, acute insulin response; CIR, corrected insulin response; HOMA-IR, homeostasis model assessment of insulin resistance; Matsuda ISI, Matsuda’s insulin sensitivity index; BIGTT-SI, insulin sensitivity index; HbA1c, glycated hemoglobin.

Table 4

PosNegORadj (95% CI)Padj
ZnT8A
AA265307
AG1611900.960(0.767-1.202)0.723
GG1724
GADA
AA426161
AG276841.240(0.959-1.603)0.101
GG368
IA-2A
AA280308
AG1861731.122(0.896-1.404)0.316
GG2321

The association between the JAZF1 rs864745 A>G variant and T1D risk stratified by islet autoantibody status.

For each autoantibody, association tests included sex, genotype, age at T1D diagnosis, and disease duration as covariates. Numbers in parentheses indicate the frequency of that genotype. P < 0.05 was considered significant.

ZnT8A, zinc transporter 8 antibody; GADA, glutamate decarboxylase antibody; IA-2A, insulinoma-associated-2 antibody.

Significant Associations With TC and LDL Levels Stratified by the JAZF1 rs864745 Variant Are Observed in T2D and T1D Status, With a Significant Heterogeneity Among Healthy Individuals

We did not find any association between this variant and plasma lipid levels (HDL, LDL, TC, TG) in the healthy individuals based on OGTT (All P > 0.05). However, the G allele of this variant was associated with lower plasma LDL and TC levels in both T2D (P = 0.003 and 0.002) and T1D status (P = 0.009 and 0.024). Compared with healthy controls, their impacts on LDL and TC levels had significant heterogeneity in T2D (P = 0.049 and 0.030) and T1D (P = 0.017 and 0.034) status, as shown in Table 5.

Table 5

GroupsAAAGGGβSEadjusted PPhetI2 (%)
Control
HDL1.37 ± 0.341.39 ± 0.341.36 ± 0.340.0060.0100.728
LDL2.81 ± 0.732.81 ± 0.782.75 ± 0.63-0.0160.0240.431
TC4.88 ± 0.964.89 ± 0.984.81 ± 0.82-0.0110.0310.588
TG1.36 ± 0.851.36 ± 0.821.49 ± 0.950.0210.0270.273
T2D
HDL1.24 ± 0.301.21 ± 0.291.22 ± 0.26-0.0430.0110.0610.09364.60
LDL2.87 ± 0.822.77 ± 0.812.67 ± 0.75-0.0710.0320.003**0.049*75.70
TC4.92 ± 1.084.80 ± 1.054.62 ± 0.96-0.0720.0420.002**0.030*78.70
TG1.84 ± 1.361.85 ± 1.331.50 ± 0.71-0.0300.0520.2040.10961.10
T1D
HDL1.36 ± 0.711.34 ± 0.391.75 ± 0.320.0610.0580.1880.22831.30
LDL2.78 ± 1.082.53 ± 0.942.49 ± 0.91-0.1380.0970.009**0.017*82.50
TC4.46 ± 2.004.09 ± 1.493.72 ± 1.94-0.1030.1430.024*0.034*77.70
TG1.61 ± 2.161.35 ± 1.401.74 ± 1.63-0.0610.6050.1850.16947.00

Correlation of plasma TC and LDL in people with different metabolic states and stratified by JAZF1 rs864745 variant.

Linear regression analysis was used, with lipid index as the dependent variable and age, sex, lipids medication and BMI (type 1 not used) as independent variables. All P values were two-tailed and P < 0.05 was considered as significant. *P < 0.05; **P < 0.01.

T2D, type 2 diabetes; T1D, type 1 diabetes; HDL, high-density lipoprotein; LDL, low-density lipoprotein; TC, total cholesterol; TG, triglyceride.

Discussion

Previous studies have found a shared genetic risk factor underlying T1D and T2D etiology (25, 26), with approximately 60 chromosome regions associated with T1D (27) and over 200 associated with T2D (28), reaching genome-wide significance. Therefore, uncovering co-localizing signals could provide biological insights into overlapping disease mechanisms, and potentially reveal therapeutic targets effective for both diseases. The JAZF1 region was reported as the genome-wide significant region in T1D (12) and T2D (13). However, the results were inconsistent. For instance, several studies revealed an association between the JAZF1 rs864745 variant and T2D, whereas others did not find any association (11, 13–17, 29–32). Our study supported that the JAZF1 rs864745 variant did not associate with T1D, T2D, or their subgroups, at least in the Chinese population.

Although obesity is a critical risk factor for the heterogeneity of T2D, and patients with different BMI have various natural histories, and β-cell mass (33), no association with this variant was reported in T2D subgroups stratified by BMI level in this study. Meanwhile, age at diagnosis is a critical factor for T1D subgroups (21). Our previous studies also confirmed a genetic correlation among T1D subgroups of diagnosed ages (19, 34). Nevertheless, our current study shows that the JAZF1 rs864745 variant is not associated with the onset age of T1D. Explanations for these discrepancies may include differences in sample size, environment, and populations studied. Because minor genetic variation accumulates over time, ancestral groups that became geographically separated many generations ago may yield different genetic risks. Additionally, the subjects in our study included community residents over 40 years old Chinese subjects, and aged participants are expected to have more cumulative environmental exposure and thus be more likely to be affected by gene-environment interaction.

Further studies have suggested that the T allele of rs864745 was associated with abnormal pancreatic β-cell function (14), which plays a role in the pathogenesis of T2D (17). JAZF1 rs864745 T allele is linked to reduced JAZF1 mRNA expression and decreased insulin release (14). Furthermore, an autosomal genomic scan showed that the JAZF1 rs864745 T variant is associated with increased fasting insulin levels (14). Our results are consistent with those of the prior studies showing that the G allele of this variant mainly improves β-cell function in normal glucose tolerance (NGT) and newly diagnosed T2D subjects, especially in those NGT participants with normal BMI levels. However, no effects on β-cell function or insulin sensitivity were observed in the overweight or obese subgroups. We speculate that the impact of the variant on islet β-cell function might be inferior to metabolic impact in people with excess body weight. These indicate that when translating genetic risk loci into clinical features of the disease state, the effect of a single variant is insufficient to achieve significant changes. As a result, multi-genic risk loci should be loaded into the estimation in future studies.

Besides, previous GWAS studies have revealed the genetic associations with islet autoantibody positivity in T1D subjects (35, 36). We evaluated the association with islet-specific autoantibody status with this variant, but no significant relationship to the positive rate of ZnT8A, GADA, or IA-2A was observed. Additionally, this variant did not affect fasting or responsive C-peptide levels in newly diagnosed T1D patients in our study. These findings indicated that the JAZF1 rs864745 A>G variant might influence islet β-cell function, but not islet autoimmunity. Further validation is warranted considering the sample size and selected population in our study.

Studies have reported that dyslipidemia is a well-established risk factor for T1D and T2D (37–39). JAZF1 is a metabolic regulator to improve lipid metabolism and resist hyperglycemia through multiple metabolic signaling pathways in T2D (18). The beneficial effects of JAZF1 on lipid metabolism have been observed in hepatocytes and adipocytes (18, 40–42). In the liver, JAZF1 enhances fatty acid β-oxidation. In adipocytes, JAZF1 inhibits the accumulation of fatty acids and triglycerides. A previous multi-ancestry population-based study has revealed multiple variants associated with serum lipid traits (43). Of these, JAZF1 rs864745 A>G variant correlates with TG and HDL. Since the non-diabetic controls recruited in our study were a subgroup of the total population, this locus was not significantly associated with HDL, LDL, TC, or TG in our normoglycemic individuals. However, we found that the G allele of this variant was associated with lower plasma LDL and TC levels in both T2D and T1D status, suggesting the relationship may be more significant during high glucose pathological states. The present findings might carry significant clinical implications. Moreover, databases from GTEx (https://gtexportal.org/), demonstrated that this variant affected JAZF1 gene expression in multi tissues, including the liver, skeletal muscle, and pancreas. Thus, we speculated that two distinct mechanisms for this variant are involved in protecting individuals from diabetes. One is by improving β-cell function, and the other is modulating lipid metabolism (via decreasing TC and LDL levels) in the hyperglycemia state.

In conclusion, we found that the JAZF1 rs864745 variant is linked to the improvement of β-cell function as well as lipid metabolism. Considering the limitation of the reproducibility with a single variant, more T1D or T2D genetic risk loci should be further studied for clinical subgroups and phenotypic of T1D and T2D patients in diabetes screening and precision therapy.

Funding

This study was supported by grants from the National Natural Science Foundation of China (81670715, 81830023, 82070803 and 81900708), and the Young Scholars Fostering Fund of the First Affiliated Hospital of Nanjing Medical University (PY2021005).

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee from the First Affiliated Hospital of Nanjing Medical University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author contributions

HD performed statistical analyses and interpretation of data and drafted the initial manuscript. HL, YQ, and LJ were responsible for the analyses and interpretation of data. HJ, MS, and HC contributed to laboratory measurements. YC and SZ contributed to the collection and selection of samples. QF and TY gave a critical revision of the manuscript. KX directed the study design and provided a critical revision of the manuscript. All the co-authors gave the final approval of the version.

Acknowledgments

The authors appreciate the help and support from all participants who took part in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.898893/full#supplementary-material

Abbreviations

T1D, type 1 diabetes; T2D, type 2 diabetes; NGT, normal glucose tolerance; TC, total cholesterol; LDL, low-density lipoprotein; HDL, high-density lipoprotein; TG, triglyceride; JAZF1, juxtaposed with another zinc finger gene 1; GWAS, genome-wide association studies; OGTT, oral glucose tolerance test; MMTT, mixed meal tolerance test; BMI, body mass index; HbA1c, glycated hemoglobin; ZnT8A, zinc transporter 8 antibody; GADA, glutamate decarboxylase antibody; IA-2A, insulinoma-associated-2 antibody; HOMA-B, homeostasis model assessment of β-cell function; HOMA-IR, homeostasis model assessment of insulin resistance; IGI, insulinogenic index; BIGTT-AIR, acute insulin response; BIGTT-SI, insulin sensitivity index; CIR, corrected insulin response; Matsuda ISI, Matsuda’s insulin sensitivity index.

References

  • 1

    KoontzJISorengALNucciMKuoFCPauwelsPvan den BergheHet al. Frequent Fusion of the JAZF1 and JJAZ1 Genes in Endometrial Stromal Tumors. Proc Natl Acad Sci U S A (2001) 98:6348–53. doi: 10.1073/pnas.101132598

  • 2

    HoMMYoganathanPChuKYKarunakaranSJohnsonJDCleeSM. Diabetes Genes Identified by Genome-Wide Association Studies are Regulated in Mice by Nutritional Factors in Metabolically Relevant Tissues and by Glucose Concentrations in Islets. BMC Genet (2013) 14:10. doi: 10.1186/1471-2582156-14-10

  • 3

    NicaACOngenHIrmingerJCBoscoDBerneyTAntonarakisSEet al. Cell-Type, Allelic, and Genetic Signatures in the Human Pancreatic Beta Cell Transcriptome. Genome Res (2013) 23:1554–62. doi: 10.1101/gr.150706.112

  • 4

    BlodgettDMNowosielskaAAfikSPechholdSCuraAJKennedyNJet al. Novel Observations From Next-Generation RNA Sequencing of Highly Purified Human Adult and Fetal Islet Cell Subsets. Diabetes (2015) 64:3172–81. doi: 10.2337/db15-0039

  • 5

    PetrenkoVSainiCGiovannoniLGobetCSageDUnserMet al. Pancreatic α- and β-Cellular Clocks Have Distinct Molecular Properties and Impact on Islet Hormone Secretion and Gene Expression. Genes Dev (2017) 31:383–98. doi: 10.1101/gad.290379.116

  • 6

    BennerCMeulenTCacéresETigyiKDonaldsonCJHuisingMO. The Transcriptional Landscape of Mouse Beta Cells Compared to Human Beta Cells Reveals Notable Species Differences in Long Non-Coding RNA and Protein-Coding Gene Expression. BMC Genomics (2014) 15:620. doi: 10.1186/1471-2164-15-620

  • 7

    KobiitaAGodbersenSAraldiEGhoshdastiderUSchmidMWSpinasGet al. The Diabetes Gene JAZF1 Is Essential for the Homeostatic Control of Ribosome Biogenesis and Function in Metabolic Stress. Cell Rep (2020) 32:107846. doi: 10.1016/j.celrep.2020.107846

  • 8

    ParkSJKwonWParkSJeongJKimDJangSet al. Jazf1 Acts as a Regulator of Insulin-Producing β-Cell Differentiation in Induced Pluripotent Stem Cells and Glucose Homeostasis in Mice. FEBS J (2021) 288:4412–27. doi: 10.1111/febs.15751

  • 9

    JangWYBaeKBKimSHYuDHKimHJJiYRet al. Overexpression of Jazf1 Reduces Body Weight Gain and Regulates Lipid Metabolism in High Fat Diet. Biochem Biophys Res Commun (2014) 444:296–301. doi: 10.1016/j.bbrc.2013.12.094

  • 10

    TaneeraJLangSSharmaAFadistaJZhouYAhlqvistEet al. A Systems Genetics Approach Identifies Genes and Pathways for Type 2 Diabetes in Human Islets. Cell Metab (2012) 16:122–34. doi: 10.1016/j.cmet.2012.06.006

  • 11

    ZegginiEScottLJSaxenaRVoightBFMarchiniJLHuTet al. Meta-Analysis of Genome-Wide Association Data and Large-Scale Replication Identifies Additional Susceptibility Loci for Type 2 Diabetes. Nat Genet (2008) 40:638–45. doi: 10.1038/ng.120

  • 12

    RobertsonCCInshawJRJOnengut-GumuscuSChenW-MSanta CruzDFYangHet al. Fine-Mapping, Trans-Ancestral and Genomic Analyses Identify Causal Variants, Cells, Genes and Drug Targets for Type 1 Diabetes. Nat Genet (2021) 53:962–71. doi: 10.1038/s41588-021-00880-5

  • 13

    SpracklenCNHorikoshiMKimYJLinKBraggFMoonSet al. Identification of Type 2 Diabetes Loci in 433,540 East Asian Individuals. Nature (2020) 582:240–5. doi: 10.1038/s41586-020-2263-3

  • 14

    GrarupNAndersenGKrarupNTAlbrechtsenASchmitzOJorgensenTet al. Association Testing of Novel Type 2 Diabetes Risk Alleles in the JAZF1, CDC123/CAMK1D, TSPAN8, THADA, ADAMTS9, and NOTCH2 Loci With Insulin Release, Insulin Sensitivity, and Obesity in a Population-Based Sample of 4,516 Glucose-Tolerant Middle-Aged Danes. Diabetes (2008) 57:2534–40. doi: 10.2337/db08-0436

  • 15

    Gamboa-MelendezMAHuerta-ChagoyaAMoreno-MaciasHVazquez-CardenasPOrdonez-SanchezMLRodriguez-GuillenRet al. Contribution of Common Genetic Variation to the Risk of Type 2 Diabetes in the Mexican Mestizo Population. Diabetes (2012) 61:3314–21. doi: 10.2337/db11-0550

  • 16

    ChangY-CLiuP-HYuY-HKuoS-SChangT-JJiangY-Det al. Validation of Type 2 Diabetes Risk Variants Identified by Genome-Wide Association Studies in Han Chinese Population: A Replication Study and Meta-Analysis. PLoS One (2014) 9:e95045. doi: 10.1371/journal.pone.0095045

  • 17

    ZhouDLiuYZhangDLiuSYuLYangYet al. Variations in/Nearby Genes Coding for JAZF1, TSPAN8/LGR5 and HHEX-IDE and Risk of Type 2 Diabetes in Han Chinese. J Hum Genet (2010) 55:810–5. doi: 10.1038/jhg.2010.117

  • 18

    LiaoZWangYQiXXiaoX. JAZF1, a Relevant Metabolic Regulator in Type 2 Diabetes. Diabetes Metab Res Rev (2019) 35. doi: 10.1002/dmrr.3148

  • 19

    ZhuMXuKChenYGuYZhangMLuoFet al. Identification of Novel T1D Risk Loci and Their Association With Age and Islet Function at Diagnosis in Autoantibody-Positive T1D Individuals: Based on a Two-Stage Genome-Wide Association Study. Diabetes Care (2019) 42:1414–21. doi: 10.2337/dc18-2023

  • 20

    PerryJRBVoightBFYengoLAminNDupuisJGanserMet al. Stratifying Type 2 Diabetes Cases by BMI Identifies Genetic Risk Variants in LAMA1 and Enrichment for Risk Variants in Lean Compared to Obese Cases. PLoS Genet (2012) 8:e1002741. doi: 10.1371/journal.pgen.1002741

  • 21

    LileyJToddJAWallaceC. A Method for Identifying Genetic Heterogeneity Within Phenotypically Defined Disease Subgroups. Nat Genet (2017) 49:310–6. doi: 10.1038/ng.3751

  • 22

    NingGThe REACTION Study Group. Risk Evaluation of Cancers in Chinese Diabetic Individuals: A Longitudinal (REACTION) Study: Letter to the Editor. J Diabetes (2012) 4:172–3. doi: 10.1111/j.1753-0407.2012.00182.x

  • 23

    GuYXiaoLGuWChenSFengYWangJet al. Rs2227982 and Rs2227981 in PDCD1 Gene are Functional SNPs Associated With T1D Risk in East Asian. Acta Diabetol (2018) 55:813–9. doi: 10.1007/s00592-018-1152-9

  • 24

    PanX-FWangLPanA. Epidemiology and Determinants of Obesity in China. Lancet Diabetes Endocrinol (2021) 9:373–92. doi: 10.1016/S2213-8587(21)00045-0

  • 25

    InshawJRJSidoreCCuccaFStefanaMICrouchDJMMcCarthyMIet al. Analysis of Overlapping Genetic Association in Type 1 and Type 2 Diabetes. Diabetologia (2021) 64:1342–7. doi: 10.1007/s00125-021-05428-0

  • 26

    AylwardAChiouJOkinoM-LKadakiaNGaultonKJ. Shared Genetic Risk Contributes to Type 1 and Type 2 Diabetes Etiology. Hum Mol Genet (2018). doi: 10.1093/hmg/ddy314

  • 27

    Onengut-GumuscuSChenW-MBurrenOCooperNJQuinlanARMychaleckyjJCet al. Fine Mapping of Type 1 Diabetes Susceptibility Loci and Evidence for Colocalization of Causal Variants With Lymphoid Gene Enhancers. Nat Genet (2015) 47:381–6. doi: 10.1038/ng.3245

  • 28

    MahajanATaliunDThurnerMRobertsonNRTorresJMRaynerNWet 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–13. doi: 10.1038/s41588-018-0241-6

  • 29

    HaimanCAFesinmeyerMDSpencerKLBuzkovaPVorugantiVSWanPet al. Consistent Directions of Effect for Established Type 2 Diabetes Risk Variants Across Populations: The Population Architecture Using Genomics and Epidemiology (PAGE) Consortium. Diabetes (2012) 61:1642–7. doi: 10.2337/db11-1296

  • 30

    NgMCYSaxenaRLiJPalmerNDDimitrovLXuJet al. Transferability and Fine Mapping of Type 2 Diabetes Loci in African Americans: The Candidate Gene Association Resource Plus Study. Diabetes (2013) 62:965–76. doi: 10.2337/db12-0266

  • 31

    TakeuchiFSerizawaMYamamotoKFujisawaTNakashimaEOhnakaKet al. Confirmation of Multiple Risk Loci and Genetic Impacts by a Genome-Wide Association Study of Type 2 Diabetes in the Japanese Population. Diabetes (2009) 58:1690–9. doi: 10.2337/db08-1494

  • 32

    OmoriSTanakaYHorikoshiMTakahashiAHaraKHiroseHet al. Replication Study for the Association of New Meta-Analysis-Derived Risk Loci With Susceptibility to Type 2 Diabetes in 6,244 Japanese Individuals. Diabetologia (2009) 52:1554–60. doi: 10.1007/s00125-009-1397-5

  • 33

    YagihashiSInabaWMizukamiH. Dynamic Pathology of Islet Endocrine Cells in Type 2 Diabetes: β-Cell Growth, Death, Regeneration and Their Clinical Implications. J Diabetes Invest (2016) 7:155–65. doi: 10.1111/jdi.12424

  • 34

    ChenSFanHFengYZhangYChenYGuYet al. The Association Between Rs1893217, Rs478582 in PTPN2 and T1D Risk With Different Diagnosed Age, and Related Clinical Characteristics in Chinese Han Population. Autoimmunity (2019) 52:95–101. doi: 10.1080/08916934.2019.1608191

  • 35

    PlagnolVHowsonJMMSmythDJWalkerNHaflerJPWallaceCet al. Genome-Wide Association Analysis of Autoantibody Positivity in Type 1 Diabetes Cases. PLoS Genet (2011) 7:e1002216. doi: 10.1371/journal.pgen.1002216

  • 36

    BrorssonCAPociotFthe Type 1 Diabetes Genetics Consortium. Shared Genetic Basis for Type 1 Diabetes, Islet Autoantibodies, and Autoantibodies Associated With Other Immune-Mediated Diseases in Families With Type 1 Diabetes. Diabetes Care (2015) 38:S8–13. doi: 10.2337/dcs15-2003

  • 37

    SongYLiuJZhaoKGaoLZhaoJ. Cholesterol-Induced Toxicity: An Integrated View of the Role of Cholesterol in Multiple Diseases. Cell Metab (2021) 33:1911–25. doi: 10.1016/j.cmet.2021.09.001

  • 38

    VergèsB. Dyslipidemia in Type 1 Diabetes: A Masked Danger. Trends Endocrinol Metab (2020) 31:422–34. doi: 10.1016/j.tem.2020.01.015

  • 39

    JamesDEStöckliJBirnbaumMJ. The Aetiology and Molecular Landscape of Insulin Resistance. Nat Rev Mol Cell Biol (2021) 22:751–71. doi: 10.1038/s41580-021-00390-6

  • 40

    MingGXiaoDGongWLiuHLiuJZhouHet al. JAZF1 can Regulate the Expression of Lipid Metabolic Genes and Inhibit Lipid Accumulation in Adipocytes. Biochem Biophys Res Commun (2014) 445:673–80. doi: 10.1016/j.bbrc.2014.02.088

  • 41

    WeiQZhouBYangGHuWZhangLLiuRet al. JAZF1 Ameliorates Age and Diet-Associated Hepatic Steatosis Through SREBP-1c -Dependent Mechanism. Cell Death Dis (2018) 9:859. doi: 10.1038/s41419-018-0923-0

  • 42

    LiLYangYYangGLuCYangMLiuHet al. The Role of JAZF1 on Lipid Metabolism and Related Genes In Vitro. Metabolism (2011) 60:523–30. doi: 10.1016/j.metabol.2010.04.021

  • 43

    GrahamSEClarkeSLWuKHKanoniSZajacGJMRamdasSet al. The Power of Genetic Diversity in Genome-Wide Association Studies of Lipids. Nature (2021) 600:675–9. doi: 10.1038/s41586-021-04064-3

Summary

Keywords

JAZF1, variant, diabetes, islet function, dyslipidemia

Citation

Dai H, Qian Y, Lv H, Jiang L, Jiang H, Shen M, Chen H, Chen Y, Zheng S, Fu Q, Yang T and Xu K (2022) Rs864745 in JAZF1, an Islet Function Associated Variant, Correlates With Plasma Lipid Levels in Both Type 1 and Type 2 Diabetes Status, but Not Healthy Subjects. Front. Endocrinol. 13:898893. doi: 10.3389/fendo.2022.898893

Received

18 March 2022

Accepted

31 May 2022

Published

01 July 2022

Volume

13 - 2022

Edited by

Katsumi Iizuka, Fujita Health University, Japan

Reviewed by

Minako Imamura, University of the Ryukyus, Japan; Junji Kozawa, Osaka University, Japan

Updates

Copyright

*Correspondence: Tao Yang, ; Kuanfeng Xu,

†These authors have contributed equally to this work and share the first authorship

This article was submitted to Diabetes: Molecular Mechanisms, a section of the journal Frontiers in Endocrinology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics