Insulin Resistance in Pregnancy Is Correlated with Decreased Insulin Receptor Gene Expression in Omental Adipose: Insulin Sensitivity and Adipose Tissue Gene Expression in Normal Pregnancy ()
Received 4 August 2015; accepted 25 February 2016; published 29 February 2016
1. Introduction
Various metabolic pathways are involved in glucose homeostasis and in the development of insulin resistance. The insulin signaling pathway is initiated when insulin binds to its receptor resulting in the tyrosine phosphorylation of insulin receptor substrates (IRS) by the insulin receptor tyrosine kinase (INSR). Tyrosine phosphorylated IRS in turn activates phosphoinositide 3-kinase (PI3K), which further activates Akt. Activated Akt induces glycogen synthesis by inhibiting glycogen synthase kinase 3 (GSK-3), and facilitates mTOR-mediated activation of protein synthesis and cell survival. The PI3k/Akt pathway further stimulates the translocation of GLUT4 vesicles to the plasma membrane and allows the uptake of glucose into muscle and adipose cells [1] -[3] .
The role of adipose tissue in the pathophysiology of insulin resistance is well recognized [1] [4] . Pregnancy is associated with insulin resistance and increases risks of gestational diabetes mellitus (GDM) and type II diabetes mellitus (T2DM) [5] [6] . In pregnancy, the cellular mechanisms for decreased insulin sensitivity are multifactorial, involving both skeletal muscle and adipose tissue and remain to be completely elucidated [1] [4] . Human pregnancy is characterized by adipose tissue accretion in early gestation, followed by insulin resistance and facilitated lipolysis in late pregnancy [7] - [9] . Additionally, the adipose tissue transcriptome demonstrates the recruitment of metabolic and immune molecular networks by 8 - 12 weeks of gestation, preceding any pregnancy- related physiological changes associated with insulin resistance [9] .
Studies investigating insulin-signaling intermediates have demonstrated reduced IRS-1 protein levels, reduced GLUT4 translocation and subsequent glucose uptake in adipose tissue of obese women with GDM, and T2DM in comparison to controls [10] - [14] .
These previous reports studied obese women and those with GDM or T2DM, analyzing mainly subcutaneous adipose tissue. There is data demonstrating that visceral adipose tissue accumulation increases the risk for developing obesity-related diseases such as T2DM, and cardiovascular diseases [15] [16] . Consequently, our study focused on the comparative roles of visceral versus subcutaneous adipose tissue. Furthermore, by choosing to include healthy, non-diabetic, non-obese pregnant women in the current analysis, we are able to examine the function of visceral and subcutaneous adipose tissue in pregnancy without confounding factors that may be associated with disease states.
The objectives of our study were two-fold. The first was to characterize in healthy, non-obese pregnant women differentially expressed genes in visceral and subcutaneous adipose tissue with a focus on insulin signaling, carbohydrate metabolism, and diabetic pathways [17] . The second was to determine associations between systemic insulin sensitivity and gene expression in visceral and subcutaneous adipose tissue.
2. Materials and Methods
2.1. Adipose Tissue Collection
Non-obese healthy pregnant women between the ages of 18 to 45 years old were recruited during prenatal care. Women who were scheduled for an elective cesarean delivery and had a normal glucose challenge test were included. After obtaining informed consent and routine prenatal care, during the cesarean delivery, maternal blood, cord blood, subcutaneous adipose tissue samples (from the abdominal wall at the edge of the surgical incision), and omental adipose tissue samples were all obtained. The adipose tissues were flash frozen in liquid nitrogen and stored in a −80˚C freezer. The Cedars-Sinai Medical Center Institutional Review Board approved the study.
2.2. RNA Extraction and Labeling
Total RNA was extracted from adipose tissues using RNeasy Lipid Tissue kit from Qiagen following the manufacturer’s instructions. Biotin-labeled cRNA synthesis and amplification was performed using Illumina Total Prep RNA Amplification kit (Ambion) following manufacturers protocol.
2.3. Microarray Hybridization, Staining and Scanning
Biotin labeled cRNA (750 ng) was loaded on to Illumina HumanHT-12 V4.0 Expression BeadChip in hybridization buffer with hybridization controls and incubated for 16 hours at 58˚C. After hybridization, BeadChips were washed, blocked with blocking buffer, and stained with streptavidin-Cy3 (Amersham Biosciences). Images of stained BeadChips were captured by Illumina BeadArray Reader.
2.4. Microarray Image Processing and Data Analysis
BeadChip images were loaded into Illumina GenomeStudio for quality determination by evaluation of the present call percentages, signal intensities of hybridization control, biotin controls, and negative controls, with quantile normalization. Signal intensities of each gene/probe were exported into Partek Genomic Suite 6.5 and ANOVA model was used to evaluate the impacts of various variables in the study such as different adipose tissues, patient ethnicity, patient-to-patient variability, etc., on the expression data set. Differentially expressed genes of subcutaneous versus omental adipose tissue were defined as a 1.5-fold difference in either direction plus pair-wise t-test. Benjamini & Hochberg procedure was applied with a false discovery rate (FDR) threshold of p < 0.05 to adjust for multiple testing. DAVID Bioinformatics Resources 6.7 was used to conduct functional classification and pathway analysis (National Institute of Allergy and Infectious Diseases, NIAID, NIH) [18] - [20] .
2.5. Real-Time PCR Validation of Microarray Results
Quantitative RT-PCR was used to verify the differential expression level of several genes. A ratio greater than 1 reflected a greater level of transcript expression within the omental adipose group compared with the subcutaneous adipose group. We selected the genes that had the highest fold change in the analysis (CLDN1, CFB) as well as six genes (IL-6, IL-1B, SLC2A1, HSD11B1, ADIPOQ, IRS1, PLA2G2A) that had well recognized insulin signaling, inflammatory and metabolic function.
RNA was first reverse-transcribed to cDNA using iScript Reverse Transcription Supermix for RT-qPCR (Bio-rad). About 10 ng cDNA were mixed with gene specific primers and Platinum Sybr Green qPCR SuperMix-UDG (Invitrogen). Real-time PCR was carried out in ABI 7000 Sequence Detector (AppliedBiosystems) for 40 cycles. The expression level of human GAPDH gene was used as an internal control, and fold changes of omental adipose group versus subcutaneous group were calculated using ΔCT method.
2.6. Matusda-Insulin Sensitivity Calculation
On 11 subjects, a frequently Sampled Oral Glucose Tolerance Test (FSOGTT) was performed. Patients had to be on a 300 grams carbohydrate per day diet. After obtaining a fasting blood sample, a 75 gram glucose solution was ingested orally. Blood samples were obtained every 15 minutes for 3 hours. Glucose and insulin levels were determined and recorded. Insulin sensitivity was calculated using Matsuda index. Matusda et al. developed an index of whole-body insulin sensitivity using 10,000/(square root of fasting glucose × fasting insulin × mean glucose × mean insulin during OGTT), which is highly correlated with the rate of whole-body glucose disposal during the euglycemic insulin clamp [21] . The homeostasis model assessment of insulin resistance (HOMA-IR) index, the product of basal glucose and insulin levels divided by 22.5 [22] , is regarded as a simple, reliable measure of insulin resistance. HOMA-IR was calculating according to the following formula: fasting insulin (μIU/ml) × fasting glucose (mmol/ml)/22.5 [23] .
2.7. Gene Mapping
DEGs were mapped to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways; a collection of pathway maps representing current knowledge on molecular interaction and reaction networks classified as Global, Metabolism, Genetic Information Processing, Environmental Information Processing, Cellular Processes, Organismal Systems, Human Diseases and Drug Development with each molecular interaction and reaction network including multiple pathways [17] . Utilizing NIH David [19] , Gene Card [24] , KEGG pathway database [17] , and PubMed [25] , extensive review was done of known functions of DEGs in regards to gene tendency towards insulin resistance (IR) or insulin sensitivity (IS). In this review, DEGs whose cellular function would cause decreased plasma glucose level or facilitate insulin action were classified as insulin sensitive and genes whose functions would cause increased plasma glucose levels or block insulin action were classified as insulin resistant.
Pearson’s correlation test was applied for transcript expression levels and insulin sensitivity scores. Separate analyses were performed for subcutaneous and visceral adipose tissues. Gene transcripts determined to be significantly correlated with both HOMA-IR score and Matsuda index were selected for pathway analysis. A cutoff p-value of <0.01 was used to determine significance.
3. Results
3.1. Demographic and Categorical Factors
The study is based on 14 patients who had an elective Cesarean section performed between 2008 and 2010. Patients were screened for any pre-existing diabetic condition; oral glucose tolerance test results were normal and excluded gestational diabetes. The mean age was 32.1 years (SD = 6.21), the mean BMI was 23.3 (SD = 2.88), the mean gestational age at delivery was 39 + 0/7 weeks (SD = 4.37 days), and the mean birth weight was 3240 grams (SD = 384.7 grams) (Table 1).
Categorical factors including fat, BMI, and race were analyzed using ANOVA to gauge each factor’s contribution and for allowing the identification of significant sources of variation. Factors having a ratio greater than 1 were considered significant and not attributable to error. Adipose tissue had almost a five-fold effect with an F ratio of 4.77. BMI had a ratio of 2.02, highlighting the impact of BMI. Anticipating this factor, the study was planned to minimize the effect of obesity by excluding from the study any women with a BMI of 30 (classified as obese). Race had a ratio of 1.8. Although the race factor is higher than can be attributed to error, comparing the interaction of race with adipose tissue on differential gene expression using a mixed model ANOVA showed an effect of only 1.25, suggesting minimal impact of race on the differential gene expression of adipose tissue (Figure 1).
Table 1. Demographic data of the study population.
Gestational Age given in number of weeks + days/7. Birth weight in grams. Age in years. BMI = Pre-pregnancy body mass index; Afr Am = African American.
3.2. Upregulated and Downregulated Genes in KEGG Pathways
In the comparison of omental versus subcutaneous adipose tissue, 920 genes showed significant relative fold change of expression level, of which 332 genes showed decreased expression levels and 588 genes showed increased expression levels.
We focused our analysis on the 17 differentially expressed genes in the insulin signaling (Figure 2), diabetic (Figure 3) and carbohydrate metabolism KEGG pathways that were associated in the literature with insulin sensitivity or resistance.
Figure 1. ANOVA model used to evaluate the impacts of the potential confounding variables: fat, race, BMI, fat and race in the differentially expressed genes in the study population. The horizontal column indicates the Mean F Ratio while the various factors analyzed in the multi factor ANOVA are listed horizontally along the x- axis. Any factor that has a Mean F Ratio more than error (1) is considered “significant”.
Figure 2. Graphic representation of KEGG insulin signaling pathway. Genes marked with stars are significantly differentially expressed in omental versus subcutaneous adipose tissue.
Figure 3. Graphic representation of KEGG Type II diabetes mellitus pathway. Genes marked with stars are significantly differentially expressed in omental versus subcutaneous adipose tissue.
3.3. Genes Differentially Expressed Towards Increased Insulin Sensitivity
2 upregulated genes in omental adipose were associated with increased insulin sensitivity: Protein phosphatase 1, regulatory (inhibitor) subunit 3C (PPP1R3C), and Phosphofructokinase, muscle (PFKM).
3 downregulated genes were associated with increased insulin resistance: Phosphorylase kinase, delta (PHK), Protein kinase C delta (PRKCD), Phosphorylase, glycogen, liver (PYGL) (Table 2).
3.4. Genes Differentially Expressed toward Increased Insulin Sensitivity
2 upregulated genes in omental adipose were associated in the literature with increased insulin resistance: FBP1 and GFPT2.
Table 2. Genes differentially expressed towards insulin sensitivity in omental vs. subcutaneous adipose.
10 downregulated DEG associated with increased insulin sensitivity were identified: GBE1, LIPE, SORBS1, IRS1, ADIPOQ, PFKFB3, ME1, IDH1, ACSS2, PFKFB1 (Table 3).
3.5. Correlation of Insulin Sensitivity and HOMA-IR with Differentially Expressed Genes
In the second stage of our study, we evaluated the correlation between insulin sensitivity (IS), HOMA-IR, with gene expression in omental and subcutaneous tissue.
In subcutaneous tissue, expression of 22 genes correlated significantly with IS. Of these, 15 had a negative correlation with IS and 7 had a positive correlation with IS. Only 5 genes matched to known genes in NIH David database with none mapping to the insulin signaling, diabetic, or carbohydrate metabolism pathways.
In omental adipose, 509 genes significantly correlated with IS. Of these, 325 matched to known genes in NIH David database. Of these, 2 genes were in the insulin signaling or type 2 diabetic KEGG pathways.
Two genes, insulin receptor (INSR) and mitogen-activated protein kinase 3 (MAPK1), strongly positively correlated with IS (Table 4).
3.6. Real-Time PCR Validation of Microarray Results
Table 5 shows the quantitative RT-PCR results for the nine positional candidate genes and confirms the differential expression of all upregulated genes and downregulated genes with the microarray analysis. For example, for the CFB gene the microarray fold change was 10.62 and the RT-PCR fold change was 11.6.
4. Discussion
Pregnancy is associated with a decrease in insulin sensitivity of 50% - 60% by the latter half of pregnancy with a 200% - 250% increase in insulin secretion to maintain euglycemia [26] . Adipose tissue plays a central role in the development of insulin resistance because it is insulin responsive and contributes directly, although quantitatively less than skeletal muscle, to the whole-body glucose disposal [1] [4] [8] [27] [28] . Our results suggest that omental tissue in comparison to subcutaneous adipose tissue may have more metabolic pathways alterations that affect carbohydrate metabolism and insulin function, which may contribute to the decreased insulin sensitivity in pregnancy.
In the insulin-signaling pathway, 4 genes were differentially expressed towards decreased insulin sensitivity in omental vs. subcutaneous adipose. IRS1 is downregulated, and it is a crucial gene in insulin sensitivity because it mediates the various cellular processes of insulin [11] [17] [29] [30] . Also downregulated are sorbin and SH3 domain containing 1 (SORBS1), required for insulin-stimulated glucose transport [31] , and Lipase, hormone-sensi- tive (LIPE), shown to have decreased activity in insulin resistant and obese states [32] [33] . Fructose-1,6-
Table 3. Genes differentially expressed towards insulin resistance in omental vs. subcutaneous adipose.
Table 4. Correlations between insulin sensitivity, HOMA-R and microarray differential gene expression mapping to KEGG insulin signaling, type 2 diabetes mellitus and carbohydrate metabolism pathways in omental tissue of healthy non-obese pregnant women.
Legend: ISI matsuda: Matusda-Insulin sensitivity index; HOMA IR: homeostasis model assessment of insulin resistance; Rho: correlation coefficient.
Table 5. Real-time PCR validation of microarray results.
Ref_Seq = reference sequence.
bisphosphatase 1 (FBP1) is upregulated, and it is associated with decreased insulin sensitivity because it regulates gluconeogenesis and has been shown to contribute to insulin secretory dysfunction [34] .
In the insulin signaling pathway, 3 genes were differentially expressed towards increased insulin sensitivity in omental vs. subcutaneous adipose. Two genes are downregulated: phosphorylase, glycogen, liver (PYGL) which is involved in glycogen degradation [35] , and calmodulin 1, phosphorylase kinase, delta (PHK), which increases glycogenolysis [36] . One gene is upregulated: protein phosphatase 1, regulatory (inhibitor) subunit 3C (PPP1R3C), which is associated with insulin sensitivity because it activates glycogen synthase, reduces glycogen phosphorylase activity and limits glycogen breakdown [37] .
In the Type-2 diabetic pathway, Adiponectin was downregulated in omental adipose compared with subcutaneous adipose. Adiponectin is an endogenous insulin-sensitizing hormone, which has anti-diabetic, anti-athero- genic and anti-inflammatory roles [30] . Decreased serum levels of adiponectin in early pregnancy have been shown to predict increased risks of GDM [4] [8] .
Several genes of carbohydrate metabolism were differentially expressed in omental versus subcutaneous adipose towards insulin resistance. 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 (PFKFB3) [38] [39] was downregulated in omental versus subcutaneous adipose. PFKFB3 is associated with increased insulin sensitivity because it is the rate-limiting enzyme in glycolysis. Glycolysis involves the anerobic breakdown of glucose to two pyruvates. Pyruvate is converted to acetyl CoA, which enters the TCA cycle. The TCA cycle is used aerobically to generate energy through the oxidization of acetate, derived from carbohydrates, fats and proteins, to carbon dioxide [17] [40] . Other downregulated genes associated with increased insulin sensitivity include other TCA and glycolytic pathway genes (Table 1), 1,4-alpha-glucan-branching enzyme 1 (GBE1) which is required for glycogen accumulation [41] [42] , and inositol 1,3,4-triphosphate 5/6 kinase (ITPK1) which inhibits TNF alpha induced cell death [43] . Glutamine-fructose-6-phosphate transaminase (GFPT2), associated with insulin resistance and oxidative stress, was upregulated [44] . Further, phosphofructokinase muscle (PFKM), a major insulin glycolytic gene [45] [46] , was downregulated, while aldehyde dehydrogenase 3 family, member B1 (ALDH3B1), which converted acetaldehyde to acetate for glycolysis [47] [48] , was upregulated in omental vs. subcutaneous adipose.
We identified 5 genes involved in the insulin signaling, diabetic and carbohydrate metabolism pathways in omental adipose that showed significant correlations with IS and HOMA-IR but none in subcutaneous adipose. Insulin sensitivity was strongly positively correlated with insulin receptor (INSR) [49] [50] gene expression in omental adipose (r = 0.84). Insulin receptor binds insulin and is thus a key mediator of the metabolic effects of insulin. MAPK1 [51] [52] , which regulates protein synthesis in response to insulin and is associated with cell survival, is also strongly positively correlated (r = 0.78) with insulin sensitivity.
Our findings support and add to previously published studies on adipose tissue. We show a decrease in expression of IRS-1 in omental adipose as compared with subcutaneous adipose, which has been observed in previous studies [10] - [12] [14] [53] . However, we did not find a difference in GLUT4 [13] [14] or in GS activity and GSK-3 expression as in other studies [14] [30] . Notably, previous reports studied diseased states (obesity, GDM and T2DM), and since we investigated only normal non-obese individuals, we would expect the cellular defects to be milder. We further found decreased expressions of adiponectin, glycolysis and TCA genes in omental adipose.
Furthermore, we showed a strong positive correlation between insulin receptor gene expression in omental adipose and systemic insulin sensitivity, with no correlation seen in subcutaneous tissue. This finding suggests the clinical significance of omental adipose of systemic insulin resistance in our study, and calls for further elucidation of the potential role of omental adipose in the pathophysiology of insulin resistance in pregnancy.
Disclosure
None of the authors have a conflict of interest
NOTES
*Corresponding author.