Increased Plasma FGF21 and Activated FGF21 Signaling in Adipose Tissue and Its Possible Association With Insulin Sensitivity in Specic GDM Subtype

Background: To study the discrepancy of the insulin sensitivity alteration pattern, circulating broblast growth factor (FGF21) levels and FGF21 signaling in visceral white adipose tissue (vWAT) of gestational diabetes mellitus (GDM) subtypes. Methods: 26 GDM women with either a predominant of insulin-secretion defect (GDM-dysfunction, n = 9) or insulin-sensitivity defect (GDM-resistance, n = 17)] and 13 normal glucose tolerance (NGT) women scheduled for caesarean-section at term were studied. Plasma and vWAT samples were collected at delivery. Results: The insulin sensitivity was improved from the 2 nd trimester to delivery in the GDM-resistance group. Elevated circulating FGF21 concentration at delivery, increased FGF receptor 1c and decreased klotho beta gene expressions, enhanced ERK1/2 phosphorylation, and increased GLUT1, IR-B, PPAR-γ gene expressions in vWAT were found in the GDM-resistance group as compared with the NGT group. The circulating FGF21 concentration was negatively correlated with fasting blood glucose (r = -0.574, P < 0.001), and associated with the GDM-resistance group (r = 0.574, P < 0.001) in pregnant women at delivery. However, we observed no insulin sensitivity alteration in GDM-dysfunction and NGT groups during pregnancy. No differences of plasma FGF21 level and FGF21 signaling in vWAT at delivery were found between women in the GDM-dysfunction and the NGT group. Conclusions: Women with GDM heterogeneity exhibited different insulin sensitivity alteration patterns. The improvement of insulin sensitivity may relate to the elevated circulating FGF21 concentration and activated FGF21 signaling in vWAT at delivery in the GDM-resistance group. across the three groups (NGT two GDM subtypes) compared using one-way ANOVA for normally distributed continuous variables, Kruskal–Wallis test for non-normally distributed continuous variables, or Chi-squared test for categorical variables. a When P 0.05, pairwise comparisons between the NGT group and each GDM group made using the Tukey’s test, Dunn’s test, or Chi-squared test, respectively. P values for pairwise comparisons were adjusted using the Bonferroni


Introduction
Gestational diabetes mellitus (GDM) is a common obstetric diseases during pregnancy with spontaneous hyperglycemia during the second and third trimesters without pre-gestational diabetes mellitus (1). GDM increases the risk for a variety of maternal and infant metabolic diseases. Women with GDM can be divided into subtypes based on insulin secretion and sensitivity heterogeneity (2). In clinical studies, we found differences in the risk factors of GDM subtypes (3), and having large-for-gestational-age infants is associated with speci c GDM subtypes (4).
FGF21 is a pleiotropic hormone-like protein regulates glucose and lipid metabolism (5), such as increasing tissue glucose uptake, improving insulin sensitivity, and inhibiting lipolysis (6). FGF21 is mainly produced by liver, and acts on target organs such as liver, adipose tissue and skeletal muscle (5,7). The main peripheral site of FGF21 that promotes glucose uptake is the white adipose tissue (WAT) (8) .
Circulating FGF21 concentration was positively associated with many clinical insulin resistance markers in pregnant women (12). However, the relationships of the circulation FGF21 levels and the GDM heterogeneity, insulin sensitivity alterations during pregnancy in women with GDM remain obscure. We aimed to elucidate the relevance of insulin sensitivity alteration pattern in GDM subtypes to their circulating FGF21 concentration at delivery, and to further investigate their FGF21 signaling in visceral WAT (vWAT) at delivery.

Participants
Pregnant Chinese women were recruited from the Northwest Women and Children's Hospital (NWCH) for participation in this case-control study. The participants were scheduled for delivery by caesarean-section from DEC. 2019 to MAR. 2020. Clinical information was collected from the 2nd-trimester study visit (24th to 28th weeks of gestation) until delivery. The 75-g, 2 h oral glucose tolerance test (OGTT) were performed at the 24th -28th gestational weeks, and GDM was diagnosed based on the criteria of the International Association of Diabetes and Pregnancy Study Groups (37). The exclusion criteria was that: diabetes before pregnancy; fasting blood glucose (FBG) ≥ 7.0 mmol/L, 2-h glucose ≥ 11.1 mmol/L after OGTT; HbA1c ≥ 6.5% in the rst trimester; eclampsia; young maternal age (less than 18); multiple pregnancy; other pre-existing diseases or gestational complications; and loss of vital data. Finally, from a total of 120 pregnant women, 57 subjects met the criteria and were included in further analysis.

Clinical Measurements and De nitions
Pre-pregnancy body mass index (pre-BMI), maternal age at delivery, and the gestational weight gain (GWG), gestational age at delivery and birth weight were noted. ISI composite index and the Stumvoll I index were adopted to evaluate the insulin sensitivity and insulin secretion, respectively. The information used for the calculation is obtained from the OGTT and simultaneous insulin-release test (13)(14)(15).
Women with normal glucose tolerance were de ned as the control group (NGT, n = 13). We classi ed GDM subtypes by Powe's de nition (2), which is based on the distributions of the ISI composite index and the Stumvoll I index in GDM women in contrast to the NGT group. GDM women with insulin sensitivity defect were de ned if the ISI composite index was under the 25th percentile of the range in the NGT group. GDM women with insulin secretion defect were de ned if the Stumvoll I index was under the 25th percentile of the range in the NGT group. Due to the limited number of participants in this study, we used the indexes of the women in the NGT group of our previous study (3) to calculate the normal range, since participants in both studies came from the same population over the same time periods. Therefore, we divided GDM women into the two subgroups: GDM with a primary defect of insulin-sensitivity (GDMresistance, n = 17), GDM with a primary defect of insulin-secretion (GDM-dysfunction, n = 9). We excluded patients who had all the stated characters, or had the two indexes over the 25th percentile.
The Collection of Blood Sample and Adipose Specimen The fasting blood samples were collected using the EDTA-coated tubes (Sarstedt, Newton, NC). At delivery, maternal blood was collected and centrifuged at 1000g for 15 minutes at 4℃. Then, the plasma was stored at -80℃ for following assays. The vWAT specimen was obtained from the greater omentum (visceral) during caesarean-section, quickly frozen in liquid nitrogen and reserved at -80℃ for following studies.
The Biochemical Parameters of Plasma FBG was detected by the glucose oxidase approach (intra-assay variation factor was 2.1% and interassay variation factor was 2.6%). Plasma lipid pro les including plasma triacylglycerol (TG), total cholesterol (CHO), low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) were detected using enzyme catalyzed approach. Plasma insulin levels were detected using commercial available kits (R-C-01-01, 5-180µU/mL, 3V Bioengineering, China). All lab tests were conducted in the certi ed lab of the NWCH with standard laboratory methods.

Quantitative Plasma FGF21 Measurement
The plasma FGF21 levels at delivery were tested via a purchasable ELISA kit (CSB-E16844 h, Cusabio Biotech, Wuhan, China). All measurements were conducted according to the manufacturers' protocol. The reference range of this detection was 15.6-1000 pg/mL with a susceptibility of 3.9 pg/mL. The coe cients of variability intra-assay and inter-assay were < 8% and < 10%, respectively.

Quantitative Polymerase Chain Reaction (qPCR)
The TRIzol reagent (Invitrogen, CA, USA) was employed to isolate total RNA from vWAT. 1µg total mRNA was reverse-transcribed into cDNA with the RT-PCR Kit (Thermo Scienti c, USA). The entire qPCR was conducted with the iQ5 PCR thermocycler (Bio-Rad, USA). The primer sequences for the tested genes were presented in Supplementary Table 1. The LightCycler protocol below was performed: 95℃ for thirty seconds (pre-cultivation); 40 periods of 95℃ for ve seconds and 60℃ for thirty seconds (amplify); and 81 periods of 55℃ for ten seconds (melting curve). We included negative controls in the entire qPCR operations. The -△△Ct method was used to identify the comparative expressing scores. Each sample was analyzed in duplicate. Cyclophilin was used as the housekeeping gene. The e ciency of each primer was coherent within experiment groups.

Western Blotting
The vWAT was homogenized with RIPA buffer (Beyotime, China) with the protease inhibitor and phosphatase inhibitor (Roche, Germany). Equal amount of protein was loaded in the 10% TGX stain-free gels (Bio-Rad) and was then transferred to PVDF membranes (Millipore). Images captured of stain-free gels were used to determine the total protein amount. The membranes were blocked with 5% non-fat dry milk and then were incubated with the primary antibodies at 4℃ overnight. The Cell Signaling Technology offered antibodies below: phosphorylated and total protein kinase B [Akt (4723/4550)], phosphorylated and total extracellular signal-regulated kinase 1/2 [Erk1/2 (4695/4370)]. Then the membranes were treated with secondary antibodies at room temperature for 1h and processed for enhanced chemiluminescence detection. The ChemiDoc Touch Imaging System (Bio-Rad) was used to visualize the total protein quantity and the targeted protein signals. Image Lab software (Bio-Rad) was used to perform the densitometric assay of the total protein and the targeted protein signals in all lanes. The ratio of phosphorylated protein value to total protein value was employed to express the changes in protein activation.

Statistical Analysis
Statistical analyses were done using SPSS 22.0 (SPSS Inc., USA). Data were displayed as means (SD or SEM) or median (IQR). One-way ANOVA was used for normal distribution data, the Kruskal-Wallis test was used for non-normal distribution continuous variables, and the Chi-squared test (or Fisher's exact possibility test) was used for class variables to compare the differences across the three groups (the NGT group and the two GDM subgroups). Tukey's test, Dunn's test and Chi-squared test were carried out to conduct the pairwise comparisons between the NGT group and the experimental groups when the P-value from either of the above tests was < 0.05. The Bonferroni correction was used to modify the P-value for Chi-squared test paired contrasts. Table 1 shows that women in the GDM-resistance group had higher pre-BMI (P = 0.006, Table 1) as compared with the NGT group. We observed no statistically signi cant differences in maternal age, GWG and infant birth weight among the three groups during pregnancy. Table 1 Clinical characteristics of women in the second and third trimester. Data are presented as n (%) for categorical variables, median (IQR, interquartile range) or mean (SD, standard deviation) for continuous variables. At the second trimester, compared with the women in the NGT group, women of the two GDM subtypes showed higher circulating levels of blood glucose during the OGTT test and larger glycemic area under the curve (AUC) (all P < 0.05). The insulin levels at all the time points and AUC (for insulin) exhibited a statistically signi cant increase (P < 0.001) in the GDM-resistance group when compared with the NGT group. By contrast, insulin level at the 2nd hour and AUC (insulin) were decreased in the GDM-dysfunction group (P < 0.001). The levels of insulin secretion indicators (HOMA2-β, Stumvoll I index) were higher in the GDM-resistance group but lower in the GDM-dysfunction group when compared with the NGT group (all P < 0.05). However, the insulin sensitivity indicators [HOMA2-S (P < 0.05) and ISI composite index (P < 0.05)] were decreased in the GDM-resistance group when compared with the NGT group. No statistically signi cant differences were found in the insulin sensitivity indicators between the women in the GDMdysfunction group and the NGT group.

Clinical characteristics of women at the second and third trimester
Before delivery, compared with women in the NGT group, women in the GDM-resistance group had larger AUC (for glucose) (P = 0.035). However, the blood glucose levels showed no statistically signi cant differences at all of the time points during the OGTT test between the GDM-resistance group and the NGT group (all P > 0.05). By contrast, women in the GDM-dysfunction group showed higher blood glucose levels during the OGTT test and larger AUC (for glucose) when compared with the NGT group (all P < 0.05). The overall tendency of other OGTT related indicators (insulin levels during OGTT test, AUC for insulin, insulin secretion and sensitivity indicators) in women of all the groups were similar from the 2nd trimester to delivery. Women in the GDM-resistance group exhibited lower level of plasma LDL-C (P = 0.033) as compared with the NGT group. Whereas, we found no statistically signi cant differences of plasma HDL-C, CHO and TG levels between women in the GDM-resistance group and the NGT group. Women in the GDM-dysfunction group had comparable lipid pro le with the NGT group.
The circulating levels of plasma FGF21 in women of the GDM-dysfunction group were similar with the NGT group. However, women in the GDM-resistance group had elevated plasma FGF21 levels when compared with the NGT group (P < 0.001).

The comparison of ISI Composite Index and Stumvoll I Index between the second trimester and before delivery
To investigate the alteration of insulin sensitivity and insulin secretion during pregnancy, we compared the ISI Composite Index and Stumvoll I Index between the second trimester and before delivery. The ISI composite index before delivery was improved when compared to that at the 2nd trimester in the GDMresistance group (P = 0.015, Fig. 1A). The ISI composite index at these two time points had no statistically signi cant differences in the GDM-dysfunction group and the NGT group. Meanwhile, we did not observe any statistically signi cant differences of the Stumoll index between the 2nd trimester and before delivery in each group (Fig. 1B).

Linear correlation analysis between the clinical indicators and circulating FGF21 at delivery
In order to nd whether plasma FGF21 level is related with the clinical indicators, we performed linear correlation analyses. The plasma FGF21 levels were positively correlated with pre-BMI (r = 0.361, P = 0.028, Table 2) and negatively related to FBG (r = -0.718, P < 0.001, Table 2) and LDL-C (r = -0.438, P = 0.007, Table 2) at delivery. However, the plasma FGF21 levels were found negatively correlated with FBG (r = -0.574, P < 0.001, Table 3) and associated with the GDM-resistance group (r = 0.574, P < 0.001, Table 3) after adjusted by each other in the multiple linear regression. Plasma FGF21 concentrations had no linear relationships with GWG, CHO, TG, HDL-C.  Table 3 Multiple linear regression analysis between the clinical indicators and circulating FGF21 at delivery. a Relationships between FGF21 and indicators correlated with FGF21 in linear correlation were conducted using multiple linear regression analysis. The NGT group is the reference in the classi cation variable. n = 39.

Relative gene and protein expression of FGF21 receptors and signaling pathways in vWAT
In order to further study the FGF21 signaling in target organs, we detected the relative gene and protein expression of FGF21 receptors and downstream signaling pathways in vWAT. Compared with the NGT group, gene expression of broblast growth factor receptor 1c (FGFR1c) was increased (P < 0.05, Fig. 2A), while gene expression of β-Klotho (KLB) was decreased (P < 0.05) in the GDM-resistance group (Fig. 2B).
Furthermore, we found no statistically signi cant differences in the phosphorylation level of Akt (Ser 473) among the three groups (Fig. 2C). However, the phosphorylation level of Erk1/2 was increased in the GDM-resistance group (P < 0.05, Fig. 2D) when compared with the NGT group. The Akt and Erk1/2 signaling was not signi cantly altered in the vWAT of GDM-dysfunction group as compared with the NGT group. 6. Linear correlation analysis between circulating FGF21 at delivery and relative expression of genes involved in glucose uptake, insulin sensitivity and lipolysis

Relative expression of genes involved in glucose uptake, insulin sensitivity and lipolysis in vWAT
In Table 4, circulating FGF21 levels were positively correlated with the relative gene expression of GLUT1 (r = 0.383, P = 0.021), IR-B (r = 0.355, P = 0.042) and PPAR-γ (r = 0.402, P = 0.015). Plasma FGF21 concentrations at delivery had no linear relationships with the relative gene expressions of GLUT4, ADIPOQ, ATGL and PLIN-1. Table 4 Linear correlation analysis between circulating FGF21 at delivery and relative expression of genes involved in glucose uptake, insulin sensitivity and lipolysis in vWAT. a Relationships between FGF21 and relative expression of genes involved in glucose uptake, insulin sensitivity and lipolysis using linear correlation analysis. n = 39.

Discussion
We found that women of the GDM-resistance group showed improved insulin sensitivity before delivery compared with that at the 2nd trimester. At the same time, increased plasma FGF21 concentrations and activated FGF21 signaling in the vWAT were found in the GDM-resistance group at delivery. Interestingly, these manifestations were not found in the women of the GDM-dysfunction group.
According to the classi cation of GDM subtypes in this study at the 2nd trimester, decreased insulin sensitivity and increased insulin secretion were clinical characters of women in the GDM-resistance group, while insu cient insulin secretion and normal insulin sensitivity were the main manifestations of women in the GDM-dysfunction group. Compared with women in the NGT group, women in the GDMresistance group showed higher pre-BMI, while women in the GDM-dysfunction group had comparable pre-BMI. Pregnant women with higher pre-BMI may be more prone to other metabolic syndrome during pregnancy.
Insulin resistance and obesity (16) affect the secretion of FGF21. Consistently, women in the GDMresistance group but not the GDM-dysfunction group manifested increased plasma FGF21 concentrations compared with the NGT group in this study. A clinical study found human BMI were positively correlated with their circulating FGF21 concentrations (17), and the increased FGF21 concentrations compensate for the insulin resistance induced by obesity and other factors (5). These ndings also explain the comparable circulating FGF21 levels between women in the GDM-dysfunction group and the NGT group, since they had similar pre-BMI and insulin sensitivity indexes.
Insulin resistance often accompanies with impaired FGF21 signal transduction (also referred as FGF21 resistance) in obese T2DM patients (18). However, the GDM-resistance group showed improved insulin sensitivity and elevated circulating FGF21 concentrations as compared with the NGT group at delivery. As we known, the increased circulating FGF21 levels were positively correlated with metabolic syndrome in obese population (19), because the physiological increased dose of circulating FGF21 helps to maintain insulin sensitivity in speci c tissues during the early stages of these diseases (20). In animal studies, increased FGF21 expression in liver and adipose tissue was found in db/db mice (19). Besides, regular exercise helps to maintain metabolic homeostasis of the GDM-resistance women. Exercise increases the sensitivity of FGF21 in adipose tissue, then improves insulin sensitivity by sending humoral signals to coordinate multi-organs (21). Moreover, unlike T2DM patients, the participants of the GDM-resistance group in our study did not have severe metabolic disorders. At the beginning of the diagnosis of GDM, obstetricians often provide exercise instruction to control their weight gain. These may explain the inconsistency of the insulin resistance and FGF21signal transduction between GDM-resistance women and T2DM patients.
Plasma FBG levels and the GDM-resistance subtype were independently correlated with plasma FGF21 concentrations. Rikke and colleagues (22) demonstrated that the physiological range of insulin increased serum FGF21 level through dose-dependent way during the euglycaemic hyperinsulinaemic clamp test.
Studies found that FGF21 was negatively correlated with FBG after adjusted by age, sex, BMI and other confounding indexes (23), and was positively correlated with adiposity and fasting plasma insulin levels in healthy subjects after adjusted by BMI and age (17). These ndings suggest that the elevated plasma FGF21 concentrations may associate with plasma insulin level, and thus play a role in improving plasma FBG.
We observed no FGF21 resistance in vWAT of the GDM-resistance group at delivery. FGFR1c and KLB are known to be the receptors of FGF21 (24). Interestingly, we found markedly increased FGFR1c expression, but reduced KLB expression in vWAT of the GDM-resistance group. Rikke and colleagues also found overweight/obesity led to decreased KLB but increased FGFR1c expression in WAT (22). FGF21 signaling through KLB in WAT may be primarily related with obesity, as decreased KLB expression was observed in WAT of obese mice (25), non-human primates fed with high-fat diet (26), and obese subjects with different levels of abnormal glucose homeostasis (27). However, FGFR1c is the predominant FGFR involved in FGF21 signaling (28). The increased FGFR1c expression could activate FGF21 signaling in a KLB-independent manner, and could compensate for the reduced KLB expression (29). Meanwhile, the GDM-resistance group showed increased Erk1/2 phosphorylation, and increased GLUT1, IR-B and PPAR-γ mRNA expression. The activation of Erk1/2 signaling in WAT could increase GLUT1 mRNA expression (5,30), stimulate PPAR-γ transcriptional activity, promote insulin-independent glucose uptake, improve insulin sensitivity and inhibit lipolysis (31,32). Meanwhile, we found FGF21 concentration was positively correlated with the relative expression of GLUT1, IR-B and PPAR-γ by the linear regression analysis. These data suggest that the FGF21 signaling pathway was activated in the vWAT of the GDM-resistance group.
We could not detect the mRNA expression of FGF21 in the vWAT, and failed to found signi cant differences in the mRNA expression of lipolysis genes and ADIPOQ in the vWAT of the GDM-resistance group. Sara and colleagues (33) also found that basal FGF21 mRNA expression was hardly detected in adipose tissue of young men, while the expression of FGF21 was signi cantly increased under supraphysiological insulin level during hyperinsulinemic euglycemic clamp. FGF21 inhibits lipolysis by reducing expression of lipid droplet-associated phosphoprotein, but not affects expression of lipolysis regulatory genes (34). FGF21 increases plasma adiponectin levels by enhancing both its gene transcription and protein secretion in adipocytes (35), however, many other factors could negatively mediate the expression of ADIPOQ, such as pro-in ammatory cytokines, endoplasmic reticulum stress, and oxidative stress(36).
The limitations of our research included that it is a single-center study with a limited sample size. Furthermore, we failed to collect plasma and vWAT samples of pregnant women at the 2nd trimester, so that we could not delineate the dynamic changes of circulating FGF21 concentrations and FGF21 signaling in the vWAT during pregnancy. Additionally, we lack the data of hyperinsulinemic and hyperglycemic clamp test during pregnancy, which prevents us from clarifying the interaction between circulating FGF21 concentrations and blood glucose/insulin levels.

Conclusions
We found women with GDM heterogeneity exhibited different insulin sensitivity alteration patterns from the 2nd trimester to delivery. The improved insulin sensitivity in the women of the GDM-resistance group may be associated with the increased FGF21 level and activated FGF21 signaling in the vWAT (Fig. 4). Our results point out a new direction for understanding the function of FGF21 in GDM women. Based on the etiology and pathogenesis of GDM heterogeneity, the dynamic balance of plasma FGF21 may help to understand the insulin sensitivity alteration patterns during pregnancy. Abbreviations GDM, gestational diabetes mellitus; NGT, normal glucose tolerance; GDM-dys, gestational diabetes mellitus dysfunction group; GDM-res, gestational diabetes mellitus resistance group; pre-BMI: prepregnancy body mass index; GWG, gestational weight gain; OGTT, oral glucose tolerance test; FBG, fasting blood glucose; AUC, area under the curve; HOMA, homeostasis model assessment; TG, triglyceride; CHO, total cholesterol; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; FGF21, broblast growth factor 21; GWG, gestational weight gain; GLUT1, glucose transporter-1; GLUT4, glucose transporter-4; IR-B, insulin receptor subunit beta; PPAR-γ, peroxisome proliferator activated receptor gamma; ADIPOQ, adiponectin, C1Q and collagen domain containing; ATGL, adipose triglyceride lipase; PLIN1, perilipin 1; FGFR1, broblast growth factor receptor 1; KLB, klotho beta.

Declarations 1. Ethics approval and consent to participate
We performed this study as per the Declaration of Helsinki. The Ethics Committee of the First A liated Hospital of Xi'an Jiaotong University approved and supervised this study (XJTU1AF2019LSL-007), and the NWCH accepted and documented this research.

Consent for publication
All subjects signed a written informed consent before performing any procedures.

Availability of data and materials
Data that support the ndings of this study are available from the corresponding authors doctorweiweixjtu@126.com upon reasonable request.