Next Article in Journal
The Longitudinal Association of Egg Consumption with Cognitive Function in Older Men and Women: The Rancho Bernardo Study
Previous Article in Journal
The Impact of Obesity on the Health of the Older Population: A Cross-Sectional Study on the Relationship between Health-Related Quality of Life and Body Mass Index across Different Age Groups
Previous Article in Special Issue
Preterm Birth and Its Association with Maternal Diet, and Placental and Neonatal Telomere Length
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

In Vitro Circadian Clock Gene Expression Assessments in Mesenchymal Stem Cells from Human Infants: A Pilot Study

by
Melissa L. Erickson
1,
Devin Dobias
2,
Madeline Rose Keleher
2,
Dana Dabelea
3,
Bryan C. Bergman
4,
Josiane L. Broussard
4,5 and
Kristen E. Boyle
2,3,*
1
Translational Research Institute, AdventHealth, Orlando, FL 32804, USA
2
Department of Pediatrics, University of Colorado Anschutz Medical Campus, Aurora, CO 80045, USA
3
The Lifecourse Epidemiology of Adiposity and Diabetes (LEAD) Center, Aurora, CO 80045, USA
4
Division of Endocrinology, Metabolism and Diabetes, University of Colorado Anschutz Medical Campus, Aurora, CO 80045, USA
5
Department of Health and Exercise Science, Colorado State University, Fort Collins, CO 80011, USA
*
Author to whom correspondence should be addressed.
Nutrients 2024, 16(1), 52; https://doi.org/10.3390/nu16010052
Submission received: 16 November 2023 / Revised: 14 December 2023 / Accepted: 20 December 2023 / Published: 23 December 2023
(This article belongs to the Special Issue Maternal Diet, Epigenetic Mechanisms and Metabolic Programming)

Abstract

:
Background: Exposure to intrauterine obesity can disrupt clock gene rhythmicity in animal models. The aim of this pilot study was to determine if maternal obesity alters rhythmic expression of core clock in mesenchymal stem cells (MSCs) from umbilical cords of human infants born to mothers with obesity (Ob-MSC) vs. normal weight (NW-MSC). Methods: We compared in vitro rhythmic expression patterns of core clock (BMAL1, CLOCK, PER2) and clock-output (NR1D1), components in undifferentiated Ob-MSCs (n = 3) vs. NW-MSCs (n = 3). MSCs were harvested every 2 h, following a dexamethasone shock, for 30 h. Adipogenesis or myogenesis was induced in vitro and markers of adipogenesis and fat storage were assessed, respectively. Results: We detected significant rhythmicity in expression patterns of BMAL1, PER2, and NR1D1 at the group level in Ob- and NW-MSCs (p < 0.05). PER2 oscillatory amplitude was 3-fold higher in Ob-MSCs vs. NW-MSCs (p < 0.006). During adipogenesis, Ob-MSCs had higher PPARγ protein content (p = 0.04) vs. NW-MSC. During myogenesis, Ob-MSCs had higher saturated triacylglycerols (p = 0.04) vs. NW-MSC. Conclusion: Rhythmic expressions of BMAL1, PER2, and NR1D1 are detectable in undifferentiated MSCs. Higher PER2 oscillatory amplitude was paralleled by higher markers of fat storage during differentiation in Ob-MSCs vs. NW-MSCs, and supports that the core clock and cellular metabolism may be linked in infant MSCs.

1. Introduction

Exposure to maternal obesity in utero is associated with obesity and metabolic disease later in life [1,2,3,4,5,6,7,8,9]. While the underlying biological mechanisms of these associations are not fully understood, it is posited that factors in the intrauterine environment have a direct effect on the metabolic function of fetal tissues. Mesenchymal stem cells (MSC) are progenitor cells for fetal mesodermal tissues that develop into adipose and skeletal muscle, which play a dominate role in whole-body metabolic health [10]. We have previously shown that MSCs isolated from umbilical cord tissue of infants born to mothers with obesity have lower total fatty acid oxidation during in vitro myogenesis, as well as greater fat content during adipogenesis and myogenesis, as compared to MSCs from infants born to mothers with normal weight [11]. Furthermore, fat content during differentiation corresponds to hypermethylation of genes involved in lipid metabolism (e.g., PRKAG2, SDHC) [11,12,13]. These findings support the notion that intrauterine exposures contribute to metabolic derangements as well as promote a fat storing phenotype in fetal MSCs, thereby increasing metabolic disease risk in offspring.
One potential mechanism by which the intrauterine exposures may impact metabolic function of fetal tissues is by modifying circadian molecular clocks. It is well accepted that maternal nutrition impacts offspring metabolism in important metabolic tissues such as skeletal muscle and adipose [14,15,16,17]. Virtually all cells in the body have an endogenous, self-sustaining molecular clock; the core of the clock mechanism is a network of transcriptional-translational feedback loops that generate rhythmic patterns of gene expression that oscillate in approximate 24 h cycles (e.g., a circadian rhythmic pattern of gene expression). The positive limb of the core clock consists of two promoters, basic helix-loop-helix ARNT like 1 (BMAL1) and clock circadian regulator (CLOCK), that heterodimerize and bind to E-box enhancer sequences to initiate the expression of the negative limb that consists of repressors period circadian regulator (PER)1/2 and cryptochrome circadian regulator (CRY)1/2. This primary loop is supported by a secondary loop consisting of nuclear receptor subfamily 1 group D member (NR1D)1/2 and retinoic acid-related orphan receptor (RORα) [18]. Although the core clock feedback loop is a self-sustaining, it is also regulated by cellular energy status, which has consequences on systemic metabolic physiology [19,20]. For example, core clock disruption in central tissues leads to obesity [21] and BMAL1 knockout in skeletal muscle leads to glucose intolerance [22]. In fact, Hansen et al. report the rhythmic oscillatory amplitude of NR1D1 in human primary myotubes is related to the metabolic health of the donor [23]. These studies in cultured myotubes also show rhythmic expression of metabolic genes, such as Sirtuin (SIRT)1, a protein deacetylase regulating cellular energy metabolism and redox pathway proteins, which is linked to in vivo health outcomes and physical activity levels [23]. Further mechanistic work in mouse embryo fibroblasts has revealed that the core clock loop is regulated by SIRT1, providing the mechanistic link whereby the core clock is intertwined with cellular energy status [24]. The metabolic enzyme glycogen synthase kinase-β (gene name GSK3B) also directly interacts with the core clock feedback loop in neuronal tissue [25], further supporting mechanistic links between core clock activity and energy metabolism.
In vitro, acute palmitate exposure dampens the rhythmic oscillatory amplitude of core clock gene expression in primary myotubes derived from human volunteers [26], suggesting that nutrient exposures can directly influence circadian rhythms. In vivo, some evidence suggests that maternal nutrition during gestation can impact the core clock machinery in offspring peripheral tissues. For example, maternal dietary restriction in goats reduces BMAL1 expression, a core component of the positive limb, in skeletal muscle of young offspring [27]. Likewise, high-fat feeding of pregnant rats reduced the area under the curve of core circadian gene expression (CLOCK, NR1D1, CRY2) in adult offspring liver [17]. Importantly, these animals also displayed disrupted metabolism and altered rhythmicity of SIRT1 [17] and peroxisome proliferator-activated receptor-α (PPARA), a master regulator of lipid metabolism gene expressions [28], suggesting links between core clock gene disruption and metabolic gene dysregulation in the context of maternal nutrition. However, whether maternal obesity during gestation impacts the core clock or rhythms of metabolic genes in progenitor cells of human fetal tissue is unknown.
Based on these previous observations showing rhythms in adult human myotubes, and disruption in animal models of maternal obesity of core clock and metabolic gene rhythmicity, the primary objective of this study was to assess in vitro rhythmic expression patterns of the core clock network in undifferentiated MSCs derived from human umbilical cords, which, to our knowledge, has not been conducted previously. We assessed the expression patterns of the core clock components BMAL1 and CLOCK of the positive limb, PER2 of the negative limb, as well as NR1D1 of the secondary regulatory loop. We hypothesized that MSCs with intrauterine exposure to obesity, which have greater lipid stores and markers of adipogenesis compared with MSCs exposed to normal weight, would display dampened oscillatory patterns of molecular circadian clock gene expression.

2. Materials and Methods

2.1. Mesenchymal Stem Cells

MSCs are from a subset of mother–infant dyads from the Healthy Start Study, in which maternal metabolic phenotyping data [6] was previously published. MSC outcomes including DNA methylation and lineage-specific metabolic phenotyping for lipid accumulation and fatty acid oxidation with adipogenesis and myogenesis were previously assessed [11]. From this initial cohort, we selected MSCs from three participants in the obesity group with previously defined hypermethylation of genes involved in lipid metabolism and MSCs from three participants in the normal weight group with hypomethylation of genes involved in lipid metabolism [13].

2.2. Participant Measures

Data collected in the Healthy Start Study included demographics, tobacco use, height, and weight. As described, women were characterized at mid-gestation (median of 17 weeks) and fasting blood samples were analyzed for glucose, insulin, triglycerides, and free fatty acids (FFA) at the University of Colorado Hospital Clinical and Translational Research Center Core Laboratory [6]. Inclusion criteria were 16 years of age or older, currently pregnant with a singleton carry (i.e., not pregnant with multiple fetuses), and 23 weeks of gestation. Exclusion criteria were prior diabetes, premature birth, and serious psychiatric illness.
Infant birth weight was obtained through medical records. Body composition was assessed within 24–48 h after birth using whole body air plethysmography (PEA POD, COSMED, Inc., Concord, CA, USA). Mixed cord blood was collected at birth and assessed for glucose, insulin, triglycerides, FFA, and adiponectin [6].

2.3. Mesenchymal Stem Cell Culture

The experimental schematic is presented in Figure 1. MSC culture and isolation procedures have been described previously [12]. In brief, plastic adherent progenitor cells were cultured from fresh umbilical cord explants in low-glucose Dulbecco’s Modified Eagle Medium (LG-DMEM, Corning Life Sciences, Corning, NY, USA) supplemented with MSC-qualified fetal bovine serum (FBS; Gibco, ThermoFisher Scientific, Waltham, MA, USA) and antibiotic/antimycotic. Cells were characterized as >98% for MSC markers CD74, CD105, and CD90 and negative for other cell types (CD34, CD45, CD19). Cells were cryopreserved and passaged and were within passages 5–7 for the experiments in this study.

2.4. Rhythmic Gene Expression Patterns in Mesenchymal Stem Cells

We assessed in vitro rhythmic expression patterns of metabolic genes under transcriptional control of the core clock including GSK3B, PPARA, and SIRT1. Undifferentiated MSCs were shocked for 1 h with (100 nM) dexamethasone (Sigma-Aldrich, St. Louis, MO, USA) followed by a 12 h rest period. Cells were then harvested every 2 h for a total of 30 h. Rhythmic expression patterns of components from the circadian clock network were assessed using quantitative PCR. Specifically, this included components from the core clock, BMAL1, well as CLOCK, and PER2, the secondary regulatory loop of the core clock, NRID1, and metabolic genes under transcriptional control of the circadian clock network, SIRT1, PPARA, and GSK3B, using quantitative PCR.

2.5. RNA Isolation and qPCR

Cells were harvested in Buffer RLT Plus (Qiagen, Germantown, MD, USA). Total RNA was isolated using RNeasy Plus mini kit (Qiagen, Germantown, MD, USA) and cDNA transcribed from 1 μg total RNA using iScript cDNA Synthesis kit (Bio-Rad Laboratories, Hercules, CA, USA). qPCR was performed using Taqman Gene Expression Assays with Taqman Fast Advanced Mastermix (Thermo Fisher, Waltham, MA, USA) for genes of interest with GUSB and PP1B as reference genes with a no-template control per gene. For MSCs differentiating into adipocytes, adiponectin mRNA expression was measured, with GUSB and PP1B as reference genes. Reactions were run in triplicate on a Viia 7 Real Time PCR System (Thermo Fisher, Waltham, MA, USA). Data were normalized using the geometric mean of the reference genes by ddCt method. Gene IDs and primer information are listed in Supplemental Table S1.

2.6. Mesenchymal Stem Cell Differentiation

MSCs underwent myogenesis or adipogenesis differentiation for 21 days as previously described [11,12]. In brief, myogenesis was induced using LG-DMEM supplemented with FBS, horse serum (Gibco, ThermoFisher, Waltham, MA, USA), dexamethasone, and hydrocortisone (Sigma-Aldrich, St. Louis, MO, USA). Adipogenesis was induced with LG-DMEM supplemented with FBS, dexamethasone, indomethacin, 3-isobutyl-1-methylxanthine, and human insulin (all, Sigma-Aldrich, St. Louis, MO, USA), treated alternately for 3 days each with LG-DMEM supplemented with FBS and insulin only. After 21 days, cells were collected for protein or fat content measures or exposed to experimental conditions for FAO assays, as previously described [11,12].

2.7. Triacylglycerol Content

We measured markers of adipogenesis and fat storage in Ob-MSCs and NW-MSCs undergoing differentiation into adipocytes and myotubes. Cell pellets were collected from 21 d myogenically differentiating cells and triacylglycerol (TAG) lipid species analyzed, as previously described [29]. Briefly, cells were lysed in PBS and then fortified with internal standards (ISs). Lipid extracted and analyzed by the Colorado Nutrition Obesity Research Center lipidomics core. Samples were run on an Sciex 2000 triple quadrupole mass spectrometer (Framingham, MA, USA). Concentration was determined by comparing ratios of unknowns to odd chain or deuterated internal standards and compared to standard curves run with standards of each lipid species.

2.8. Protein Content

Cells were harvested from 21 d adipogenic differentiating cells in lysis buffer (CelLyticTM MT, Sigma-Aldrich, St. Louis, MO, USA) supplemented with protease and phosphatase inhibitor cocktails (Sigma-Aldrich, St. Louis, MO, USA). Total protein was determined by BCA assay. Protein content of peroxisome proliferator-activated receptor (PPAR)γ with β-actin (Cell Signaling Technology, Danvers, MA, USA) as reference control was assessed by Simple Western size-based protein assay (WES, ProteinSimple, Santa Clara, CA, USA) following the manufacturer’s protocol. Results from WES were analyzed using ProteinSimple Compass software, v 6.1.0. Antibodies were optimized in-house for this system and antibody specifics and assay conditions are listed in Supplemental Table S2.

2.9. Statistical Analysis

Participant demographic data are presented as mean ± SD. Student’s unpaired, two-tailed t-tests were used to compare participant characteristics in NW vs. Ob groups. All cellular data are presented as mean ± SE unless otherwise noted. Circadian rhythmicity of gene expression was evaluated using CircWave v 1.4, which determines goodness of fit to a cosine wave. For gene targets in which circadian rhythmicity was significantly detected, circadian amplitude was subsequently determined. Student’s two-tailed t-tests were used to compare group differences in circadian oscillatory amplitude of gene expression between NW-MSCs vs. Ob-MSCs. Student’s two-tailed t-tests were also used to compare group differences in cell-based assays between NW-MSCs vs. Ob-MSCs. We conducted preliminary comparisons of expression patterns of core clock components and metabolic clock-output genes in MSCs derived from infants born to mothers with obesity (Ob-MSC) versus normal weight (NW-MSC) to represent diverse in utero exposures using Pearson correlation analyses to test for associations between outcomes of interest. Statistical analyses were conducted with GraphPad Prism (v 9.5.1).

3. Results

3.1. Participant Characteristics

Data are presented from 12 individual participants, resulting in 6 mother–infant dyads. Data on the methylation status of genes from MSCs from these participants were previously published as part of a larger dataset [11]. By design, mothers with obesity (Ob) had significantly higher BMIs than mothers with normal weight (NW; p = 0.002; Table 1). Additionally, maternal insulin was higher in the Ob as compared to the NW group (p = 0.049). Neonatal adiposity (g of fat and % fat mass) tended to be higher in infants born to mothers with obesity (p = 0.07 for both).

3.2. Infant MSCs Have Detectible Rhythmic Expression Patterns in Components of the Circadian Clock Network

Rhythmic gene expression patterns of components from the core clock (BMAL1, CLOCK, PER2) and secondary regulatory loop (NR1D1) from NW-MSCs and Ob-MSCs are shown in Figure 2A–H. Significant circadian rhythms were detected in BMAL1, PER2, and NR1D1 in both NW-MSCs and Ob-MSCs at the group level (Table 2).
Rhythmic gene expression patterns for metabolic gene under transcriptional control of the circadian clock network (GSK3B, PPARA, and SIRT1) from NW-MSCs and Ob-MSCs are shown in Figure 3A–F. Circadian rhythms were not detected in any metabolic genes examined (Table 2).

3.3. PER2 Oscillatory Amplitude Is Higher in Ob-MSCs as Compared to NW-MSCs

Circadian oscillatory amplitude was calculated for all genes whose expression displayed a significant circadian rhythm. There were no differences between groups in oscillatory amplitude for BMAL1 or NR1D1 (p = 0.95 and p = 0.13, respectively, Figure 4A,C). In contrast, the oscillatory amplitude of PER2 was significantly higher in MSC-Ob as compared to MSC-NW (p = 0.006, Figure 4B).

3.4. Ob-MSCs Have Greater Markers of Adipogenesis and Higher TAGs during Myogenesis

During adipogenesis, Ob-MSCs had tended to have higher adiponectin mRNA content (p = 0.09; Figure 5A) and had higher protein levels of PPARγ (p = 0.04, Figure 5B) compared to NW-MSCs. During myogenesis, Ob-MSCs had greater TAG content (p = 0.04; Figure 5C) and tended to have higher saturated TAG content (p = 0.09; Figure 5D) compared to NW-MSC.

3.5. PER2 Amplitude May Be Linked to MSC Phenotype as Well as Maternal and Infant Traits

Associations between the oscillatory amplitude of PER2 in undifferentiated MSCs and various outcomes are presented in Supplementary, including MSC markers of adipogenesis (Supplemental Figure S1), maternal insulin (Supplemental Figure S2), and neonatal adiposity (Supplemental Figure S3). PER2 was selected for exploratory analyses based on the observed group differences in rhythmic amplitude between NW-MSCs and Ob-MSCs. p-values not reported due to limited sample size.

4. Discussion

Circadian expression patterns of the molecular circadian clock network including core clock components BMAL1 and PER2, as well as a component of the secondary regulatory loop, NR1D1, were present in MSCs derived from human umbilical cord tissue. To our knowledge, this is the first study to demonstrate circadian gene expression in human-derived primary stem cells prior to differentiation into tissue-specific lineages such as myotubes or adipocytes. One advantage of the MSC model is that it reflects infant phenotype in conjunction with in utero exposures, prior to environmental or lifestyle factors that occur after birth (e.g., diet, sleep/wake cycles) [10]. Thus, oscillatory rhythms of gene expression in MSCs may represent intrinsic circadian rhythms in infant tissues supporting the use of MSCs as an approach to investigate the impact of circadian regulation during human growth and development.
The oscillatory patterns observed in components of the circadian clock network in MSCs were consistent with previous reports of other in vitro human derived cell-based assays. For example, data from experiments in primary human myotubes show that BMAL1, PER2, and NR1D1 gene expression follows circadian rhythmicity in cells derived from groups of metabolically distinct participants [23]. Likewise, work in human progenitor cells derived from adipose tissue show robust oscillatory expression patterns of BMAL1, PER2, and NR1D1 [30]. Similarly, in a rodent model of maternal obesity, adult offspring of obese dams show disruption to CLOCK, NR1D1, and CRY2, though this was mostly evident when offspring were weaned to a high fat diet [17]. Here, we found that CLOCK expression patterns were not rhythmic in MSCs, which is consistent with previous studies in human myotubes [23,26,31,32]. In addition, gene expression patterns of metabolic genes under transcriptional control of the core clock, including SIRT1, PPARA, and GSK3B, were not rhythmic in MSCs. This finding is also consistent with previous studies in human myotubes derived from sedentary lean, obese, and T2D phenotypes. Adult human myotubes and adipocytes originate from the fetal MSC lineage, and perhaps this explains the observed consistencies in cells from infants and adults.
The circadian amplitude of PER2 expression was higher in Ob-MSCs compared to NW-MSCs. This finding is somewhat inconsistent with previous reports of dampened oscillatory amplitude of NR1D1 expression in primary myotubes from participants with T2D relative to myotubes from participants across the metabolic spectrum of obese, sedentary lean, and exercise trained lean [23]. Similarly, maternal obesity combined with a postnatal high fat diet reduces oscillatory amplitude of clock gene expression in mice offspring [17], contrary to the increased amplitude we observed here. The primary difference between the present study and previous reports of dampened oscillatory amplitude of clock gene expression in response to obesogenic exposures is that previous work was conducted in fully differentiated cells, whereas the present study used undifferentiated umbilical cord-derived MSCs. Our results showing increased PPARγ, as well as increased storage of saturated TAG during myogenesis in Ob-MSCs [11,12,13], supports an elevated fat-storing phenotype during adipogenesis and myogenesis in parallel with higher PER2 oscillatory amplitude, compared to NW-MSCs. Thus, perhaps higher PER2 oscillatory amplitude is related more to the increased fat-storing ability of Ob-MSCs rather than maternal obesity per se. In support, increased PER2 amplitude appears to correlate with markers of adipogenesis, including TAG content during myogenesis, as well as adiponectin and PPARγ during adipogenesis. Given that insulin sensitivity is required for lipid storage, the MSC fat storing phenotype may also represent a more insulin sensitive state; thus, a higher PER2 oscillatory amplitude in Ob-MSCs is not inconsistent with a dampened amplitude in myotubes derived from participants with impaired versus normal glucose tolerance [23]. However, we did not assess insulin sensitivity in MSCs in this study.
PER2 is a core component in the repressor limb of the core clock feedback loop [18]. The higher PER2 oscillatory amplitude we observed in undifferentiated Ob-MSCs vs. NW-MSCs, was accompanied by higher fat storage during experimental cellular differentiation. These observations support the relationship between core clock activity and cellular metabolism. Consistent with our cell-based assays, genetic knockout models of PER2 show reduced triacylglycerol and non-esterified fatty acids, as well as a lean phenotype suggestive that the negative limb of core clock plays a role in whole-body fat storage [33]. In addition, PER2 silencing leads to reduced lipid droplet accumulation and reduced cellular triglyceride content [34]. Future studies are warranted to examine the relationship between core clock, specifically robust circadian expression of PER2, and fat storage and adipose tissue insulin sensitivity. While we did not observe statistical links between PER2 oscillatory amplitude and the other metabolic genes tested (SIRT1, PPARA, GSK3B), this may be potentially due to the small sample size of our preliminary study.
Fetal exposures to excess metabolic substrates during maternal obesity are hypothesized to contribute to increased neonatal adiposity at birth [35], and our preliminary findings align with this hypothesis. For example, we observed that maternal insulin during late gestation may correlate with PER2 amplitude in MSCs, suggesting that elevated insulin in utero may be one factor related to increased PER2 circadian expression in fetal tissues. We also observed that PER2 amplitude may correlate with neonatal adiposity at birth. This supports the notion that increased PER2 circadian expression in fetal tissues is associated with fetal fat accretion, subsequently contributing to increased neonatal adiposity detectable at birth. These early findings set the stage for future hypothesis testing regarding the role of circadian pathways in the context of intergenerational obesity transmission. It is worth noting that in addition to maternal obesity-related exposures potentially impacting cellular rhythms, there may be other factors contributing to intrinsic infant MSC rhythm. These may include physical activity, dietary choices, sleep habits, or medications; however, this line of investigation is largely unexplored.

5. Conclusions

This study presents evidence for robust, cell autonomous, rhythmic oscillations of components of the molecular circadian clock network, which are detectable in vitro in undifferentiated MSCs derived from human infants. We also observed increased PER2 oscillatory amplitude in Ob-MSCs compared to NW-MSCs. MSC phenotype reflects infant tissue; and therefore, this experimental approach serves as a useful tool for investigating circadian pathways that may contribute to metabolic disease risk in children exposed to obesity in utero. Future studies using larger sample sizes of MSCs from human infants are needed to better understand the role of the circadian clock network during fetal growth and development and relations to future obesity risk.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu16010052/s1, Table S1: qPCR TaqMan Gene Expression Assay Information; Table S2: Antibodies for simple Western; Figure S1: Associations between PER2 amplitude and cell-based assays during adipogenesis; Figure S2: Associations between maternal insulin and PER2 amplitude; Figure S3: Associations between PER2 amplitude and neonatal adiposity.

Author Contributions

M.L.E., M.R.K., J.L.B. and K.E.B. conceived this project. M.L.E., D.D. (Devin Dobias), M.R.K., B.C.B., J.L.B. and K.E.B. performed experiments and collected and analyzed the data. M.L.E., J.L.B. and K.E.B. interpreted the results. M.L.E., J.L.B. and K.E.B. drafted the manuscript. D.D. (Dana Dabelea) conceptualized and implemented the parent Healthy Start study. K.E.B. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by NIH R01DK117168 and the American Diabetes Association (1-18-ITCS-016) to KEB, R01DK125653 to JLB and a pilot award to JLB and KEB through the NIH/NCATS Colorado CTSA Grant Number UL1 TR002535. The Healthy Start BabyBUMP Project was supported by the American Heart Association (14PRE18230008), the Colorado Nutrition and Obesity Research Center (NORC; NIH P30DK048520), and by the parent Healthy Start study: (NIH R01DK076648, UH3OD023248, DD) and NIH/NCATS Colorado CTSA (UL1TR001082). MLE is supported in part by NIH K01DK134838.

Institutional Review Board Statement

Umbilical cord tissue samples used in this study were originally collected from participants volunteering in the Healthy Start Study (Clinical Trials.gov, NCT02273297; 2014-10-21), which was approved by the Colorado Multiple Institutional Review Board (COMIRB) in July 2009.

Informed Consent Statement

Written informed consent was obtained from all study participants prior to study inclusion.

Data Availability Statement

The datasets available on request: The raw data supporting the conclusions of this article will be made available by the authors, without personal health or identifying information.

Conflicts 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.

References

  1. Boney, C.M.; Verma, A.; Tucker, R.; Vohr, B.R. Metabolic syndrome in childhood: Association with birth weight, maternal obesity, and gestational diabetes mellitus. Pediatrics 2005, 115, e290–e296. [Google Scholar] [CrossRef] [PubMed]
  2. Ismail-Beigi, F.; Catalano, P.M.; Hanson, R.W. Metabolic programming: Fetal origins of obesity and metabolic syndrome in the adult. Am. J. Physiol. Endocrinol. Metab. 2006, 291, E439–E440. [Google Scholar] [CrossRef] [PubMed]
  3. Oken, E.; Taveras, E.M.; Kleinman, K.P.; Rich-Edwards, J.W.; Gillman, M.W. Gestational weight gain and child adiposity at age 3 years. Am. J. Obstet. Gynecol. 2007, 196, 322.e1–e8. [Google Scholar] [CrossRef] [PubMed]
  4. Tam, C.H.T.; Ma, R.C.W.; Yuen, L.Y.; Ozaki, R.; Li, A.M.; Hou, Y.; Chan, M.H.M.; Ho, C.S.; Yang, X.; Chan, J.C.N.; et al. The impact of maternal gestational weight gain on cardiometabolic risk factors in children. Diabetologia 2018, 61, 2539–2548. [Google Scholar] [CrossRef] [PubMed]
  5. Stotland, N.E.; Cheng, Y.W.; Hopkins, L.M.; Caughey, A.B. Gestational weight gain and adverse neonatal outcome among term infants. Obstet. Gynecol. 2006, 108, 635–643. [Google Scholar] [CrossRef]
  6. Starling, A.P.; Brinton, J.T.; Glueck, D.H.; Shapiro, A.L.; Harrod, C.S.; Lynch, A.M.; Siega-Riz, A.M.; Dabelea, D. Associations of maternal BMI and gestational weight gain with neonatal adiposity in the Healthy Start study. Am. J. Clin. Nutr. 2015, 101, 302–309. [Google Scholar] [CrossRef] [PubMed]
  7. Hillier, T.A.; Pedula, K.L.; Vesco, K.K.; Oshiro, C.E.; Ogasawara, K.K. Impact of Maternal Glucose and Gestational Weight Gain on Child Obesity over the First Decade of Life in Normal Birth Weight Infants. Matern. Child. Health J. 2016, 20, 1559–1568. [Google Scholar] [CrossRef]
  8. Whitaker, R.C. Predicting preschooler obesity at birth: The role of maternal obesity in early pregnancy. Pediatrics 2004, 114, e29–e36. [Google Scholar] [CrossRef]
  9. Catalano, P.M.; Farrell, K.; Thomas, A.; Huston-Presley, L.; Mencin, P.; de Mouzon, S.H.; Amini, S.B. Perinatal risk factors for childhood obesity and metabolic dysregulation. Am. J. Clin. Nutr. 2009, 90, 1303–1313. [Google Scholar] [CrossRef]
  10. Marino, L.; Castaldi, M.A.; Rosamilio, R.; Ragni, E.; Vitolo, R.; Fulgione, C.; Castaldi, S.G.; Serio, B.; Bianco, R.; Guida, M.; et al. Mesenchymal Stem Cells from the Wharton’s Jelly of the Human Umbilical Cord: Biological Properties and Therapeutic Potential. Int. J. Stem Cells 2019, 12, 218–226. [Google Scholar] [CrossRef]
  11. Boyle, K.E.; Patinkin, Z.W.; Shapiro, A.L.B.; Bader, C.; Vanderlinden, L.; Kechris, K.; Janssen, R.C.; Ford, R.J.; Smith, B.K.; Steinberg, G.R.; et al. Maternal obesity alters fatty acid oxidation, AMPK activity, and associated DNA methylation in mesenchymal stem cells from human infants. Mol. Metab. 2017, 6, 1503–1516. [Google Scholar] [CrossRef] [PubMed]
  12. Boyle, K.E.; Patinkin, Z.W.; Shapiro, A.L.; Baker, P.R., 2nd; Dabelea, D.; Friedman, J.E. Mesenchymal Stem Cells From Infants Born to Obese Mothers Exhibit Greater Potential for Adipogenesis: The Healthy Start BabyBUMP Project. Diabetes 2016, 65, 647–659. [Google Scholar] [CrossRef] [PubMed]
  13. Erickson, M.L.; Patinkin, Z.W.; Duensing, A.M.; Dabelea, D.; Redman, L.M.; Boyle, K.E. Maternal metabolic health drives mesenchymal stem cell metabolism and infant fat mass at birth. JCI Insight 2021, 6, e146606. [Google Scholar] [CrossRef] [PubMed]
  14. Du, M.; Yan, X.; Tong, J.F.; Zhao, J.; Zhu, M.J. Maternal obesity, inflammation, and fetal skeletal muscle development. Biol. Reprod. 2010, 82, 4–12. [Google Scholar] [CrossRef] [PubMed]
  15. Sandoval, C.; Wu, G.; Smith, S.B.; Dunlap, K.A.; Satterfield, M.C. Maternal Nutrient Restriction and Skeletal Muscle Development: Consequences for Postnatal Health. Adv. Exp. Med. Biol. 2020, 1265, 153–165. [Google Scholar] [CrossRef] [PubMed]
  16. Alfaradhi, M.Z.; Kusinski, L.C.; Fernandez-Twinn, D.S.; Pantaleão, L.C.; Carr, S.K.; Ferland-McCollough, D.; Yeo, G.S.; Bushell, M.; Ozanne, S.E. Maternal Obesity in Pregnancy Developmentally Programs Adipose Tissue Inflammation in Young, Lean Male Mice Offspring. Endocrinology 2016, 157, 4246–4256. [Google Scholar] [CrossRef] [PubMed]
  17. Borengasser, S.J.; Kang, P.; Faske, J.; Gomez-Acevedo, H.; Blackburn, M.L.; Badger, T.M.; Shankar, K. High fat diet and in utero exposure to maternal obesity disrupts circadian rhythm and leads to metabolic programming of liver in rat offspring. PLoS ONE 2014, 9, e84209. [Google Scholar] [CrossRef] [PubMed]
  18. Ko, C.H.; Takahashi, J.S. Molecular components of the mammalian circadian clock. Hum. Mol. Genet. 2006, 15, R271–R277. [Google Scholar] [CrossRef]
  19. Bass, J.; Lazar, M.A. Circadian time signatures of fitness and disease. Science 2016, 354, 994–999. [Google Scholar] [CrossRef]
  20. Panda, S. Circadian physiology of metabolism. Science 2016, 354, 1008–1015. [Google Scholar] [CrossRef]
  21. Turek, F.W.; Joshu, C.; Kohsaka, A.; Lin, E.; Ivanova, G.; McDearmon, E.; Laposky, A.; Losee-Olson, S.; Easton, A.; Jensen, D.R.; et al. Obesity and metabolic syndrome in circadian Clock mutant mice. Science 2005, 308, 1043–1045. [Google Scholar] [CrossRef] [PubMed]
  22. Harfmann, B.D.; Schroder, E.A.; Kachman, M.T.; Hodge, B.A.; Zhang, X.; Esser, K.A. Muscle-specific loss of Bmal1 leads to disrupted tissue glucose metabolism and systemic glucose homeostasis. Skelet. Muscle 2016, 6, 12. [Google Scholar] [CrossRef] [PubMed]
  23. Hansen, J.; Timmers, S.; Moonen-Kornips, E.; Duez, H.; Staels, B.; Hesselink, M.K.; Schrauwen, P. Synchronized human skeletal myotubes of lean, obese and type 2 diabetic patients maintain circadian oscillation of clock genes. Sci. Rep. 2016, 6, 35047. [Google Scholar] [CrossRef] [PubMed]
  24. Nakahata, Y.; Sahar, S.; Astarita, G.; Kaluzova, M.; Sassone-Corsi, P. Circadian control of the NAD+ salvage pathway by CLOCK-SIRT1. Science 2009, 324, 654–657. [Google Scholar] [CrossRef] [PubMed]
  25. Besing, R.C.; Paul, J.R.; Hablitz, L.M.; Rogers, C.O.; Johnson, R.L.; Young, M.E.; Gamble, K.L. Circadian rhythmicity of active GSK3 isoforms modulates molecular clock gene rhythms in the suprachiasmatic nucleus. J. Biol. Rhythms 2015, 30, 155–160. [Google Scholar] [CrossRef] [PubMed]
  26. Sardon Puig, L.; Pillon, N.J.; Naslund, E.; Krook, A.; Zierath, J.R. Influence of obesity, weight loss, and free fatty acids on skeletal muscle clock gene expression. Am. J. Physiol. Endocrinol. Metab. 2020, 318, E1–E10. [Google Scholar] [CrossRef] [PubMed]
  27. Zhou, X.; Yan, Q.; Yang, H.; Ren, A.; He, Z.; Tan, Z. Maternal intake restriction programs the energy metabolism, clock circadian regulator and mTOR signals in the skeletal muscles of goat offspring probably via the protein kinase A-cAMP-responsive element-binding proteins pathway. Anim. Nutr. 2021, 7, 1303–1314. [Google Scholar] [CrossRef]
  28. Bougarne, N.; Weyers, B.; Desmet, S.J.; Deckers, J.; Ray, D.W.; Staels, B.; De Bosscher, K. Molecular Actions of PPARα in Lipid Metabolism and Inflammation. Endocr. Rev. 2018, 39, 760–802. [Google Scholar] [CrossRef]
  29. Chaves, A.B.; Zheng, D.; Johnson, J.A.; Bergman, B.C.; Patinkin, Z.W.; Zaegel, V.; Biagioni, E.M.; Krassovskaia, P.; Broskey, N.T.; May, L.E.; et al. Infant Mesenchymal Stem Cell Insulin Action Is Associated with Maternal Plasma Free Fatty Acids, Independent of Obesity Status: The Healthy Start Study. Diabetes 2022, 71, 1649–1659. [Google Scholar] [CrossRef]
  30. Wu, X.; Zvonic, S.; Floyd, Z.E.; Kilroy, G.; Goh, B.C.; Hernandez, T.L.; Eckel, R.H.; Mynatt, R.L.; Gimble, J.M. Induction of circadian gene expression in human subcutaneous adipose-derived stem cells. Obesity 2007, 15, 2560–2570. [Google Scholar] [CrossRef]
  31. Bonaconsa, M.; Malpeli, G.; Montaruli, A.; Carandente, F.; Grassi-Zucconi, G.; Bentivoglio, M. Differential modulation of clock gene expression in the suprachiasmatic nucleus, liver and heart of aged mice. Exp. Gerontol. 2014, 55, 70–79. [Google Scholar] [CrossRef] [PubMed]
  32. Bjarnason, G.A.; Jordan, R.C.; Wood, P.A.; Li, Q.; Lincoln, D.W.; Sothern, R.B.; Hrushesky, W.J.; Ben-David, Y. Circadian expression of clock genes in human oral mucosa and skin: Association with specific cell-cycle phases. Am. J. Pathol. 2001, 158, 1793–1801. [Google Scholar] [CrossRef] [PubMed]
  33. Grimaldi, B.; Bellet, M.M.; Katada, S.; Astarita, G.; Hirayama, J.; Amin, R.H.; Granneman, J.G.; Piomelli, D.; Leff, T.; Sassone-Corsi, P. PER2 controls lipid metabolism by direct regulation of PPARγ. Cell Metab. 2010, 12, 509–520. [Google Scholar] [CrossRef] [PubMed]
  34. Jing, Y.; Chen, Y.; Wang, S.; Ouyang, J.; Hu, L.; Yang, Q.; Wang, M.; Zhang, B.; Loor, J.J. Circadian Gene PER2 Silencing Downregulates PPARG and SREBF1 and Suppresses Lipid Synthesis in Bovine Mammary Epithelial Cells. Biology 2021, 10, 1226. [Google Scholar] [CrossRef]
  35. Barbour, L.A.; Hernandez, T.L. Maternal Lipids and Fetal Overgrowth: Making Fat from Fat. Clin. Ther. 2018, 40, 1638–1647. [Google Scholar] [CrossRef]
Figure 1. Experimental schematic. MSCs were isolated from umbilical cords of infants from mothers with normal weight or obesity. Gene expression was assessed in undifferentiated NW-MSCs and Ob-MSCs. Cell-based phenotyping was conducted in NW-MSCs and Ob-MSCs differentiating into myotubes and adipocytes.
Figure 1. Experimental schematic. MSCs were isolated from umbilical cords of infants from mothers with normal weight or obesity. Gene expression was assessed in undifferentiated NW-MSCs and Ob-MSCs. Cell-based phenotyping was conducted in NW-MSCs and Ob-MSCs differentiating into myotubes and adipocytes.
Nutrients 16 00052 g001
Figure 2. MSCs were harvested every 2 h for 30. Gene expression was measured, and the temporal patterns of expression were fit to a cosine wave and double plotted with CircWave for components of the circadian clock network, including components of the core clock (BMAL1, CLOCK, and PER2) and the secondary regulatory loop (NR1D1) in NW-MSCs and Ob-MSC at the individual (Top row; (AD)) and group (Bottom row; (EH)) level.
Figure 2. MSCs were harvested every 2 h for 30. Gene expression was measured, and the temporal patterns of expression were fit to a cosine wave and double plotted with CircWave for components of the circadian clock network, including components of the core clock (BMAL1, CLOCK, and PER2) and the secondary regulatory loop (NR1D1) in NW-MSCs and Ob-MSC at the individual (Top row; (AD)) and group (Bottom row; (EH)) level.
Nutrients 16 00052 g002
Figure 3. MSCs were harvested every 2 h for 30. Gene expression was measured, and the temporal patterns of expression were fit to a cosine wave and double plotted with CircWave for components for metabolic targets under transcriptional control of the circadian clock network (e.g., clock output) including GSK3B, PPARA, and SIRT1 in NW-MSCs and Ob-MSC at the individual (Top row (AC)) and group (Bottom row (DF)) level.
Figure 3. MSCs were harvested every 2 h for 30. Gene expression was measured, and the temporal patterns of expression were fit to a cosine wave and double plotted with CircWave for components for metabolic targets under transcriptional control of the circadian clock network (e.g., clock output) including GSK3B, PPARA, and SIRT1 in NW-MSCs and Ob-MSC at the individual (Top row (AC)) and group (Bottom row (DF)) level.
Nutrients 16 00052 g003
Figure 4. Oscillatory amplitude of circadian gene expression of BMAL1 (A), PER2 (B), and NR1D1 (C) in NW-MSCs and Ob-MSCs. ns indicates not significant (p > 0.05) and ** indicates p < 0.01.
Figure 4. Oscillatory amplitude of circadian gene expression of BMAL1 (A), PER2 (B), and NR1D1 (C) in NW-MSCs and Ob-MSCs. ns indicates not significant (p > 0.05) and ** indicates p < 0.01.
Nutrients 16 00052 g004
Figure 5. (A): Adiponectin mRNA content; and (B): PPARγ protein content during adipogenesis; (C): Total triacylglycerols (TAGs) during myogenesis in NW-MSCs and Ob-MSCs; and (D): Saturated triacylglycerols (TAGs) during myogenesis in NW-MSCs and Ob-MSCs. * indicates p < 0.05.
Figure 5. (A): Adiponectin mRNA content; and (B): PPARγ protein content during adipogenesis; (C): Total triacylglycerols (TAGs) during myogenesis in NW-MSCs and Ob-MSCs; and (D): Saturated triacylglycerols (TAGs) during myogenesis in NW-MSCs and Ob-MSCs. * indicates p < 0.05.
Nutrients 16 00052 g005
Table 1. Characteristics of Mother–Infant Dyads.
Table 1. Characteristics of Mother–Infant Dyads.
Normal Weight
n = 3
Obese
n = 3
p
Maternal Characteristics
  Age (yrs)26.3 ± 6.433.7 ± 9.70.335
  Pre-pregnancy BMI (kg/m2)21.2 ± 1.732.3 ± 2.00.002 *
  Race/ethnicity (self-reported), N 0.519
    Hispanic11
    White, non-Hispanic12
    Black, non-Hispanic00
    All others combined10
  Fasting Plasma Glucose (mg/dL)75.0 ± 3.674.3 ± 8.00.902
  Fasting Plasma Insulin (mU/dL)9.0 ± 2.013.7 ± 2.10.049 *
Infant Characteristics
  Sex (M/F)2/12/1--
  Birthweight (kg)3.2 ± 0.093.4 ± 0.30.244
  Fat Free Mass (g)2816.0 ± 154.32906.1 ± 260.70.632
  Fat Mass (g)191.7 ± 54.0395.1 ± 136.00.074
  Fat Mass (%)6.4 ± 2.011.9 ± 3.40.074
Maternal fasting plasma assessments occurred during late gestation. MSC: Mesenchymal stem cells; Values are presented as average ± SD. * p ≤ 0.05.
Table 2. Circadian Rhythmicity of Gene Targets in Undifferentiated MSCs.
Table 2. Circadian Rhythmicity of Gene Targets in Undifferentiated MSCs.
Gene TargetNormal Weight
n = 3
Obese
n = 3
R2pR2p
CLOCK0.1170.1360.0190.716
BMAL10.3470.001 *0.2220.021 *
PER20.3670.001 *0.5050.000 *
NR1D10.2290.012 *0.4440.000 *
GSK3B0.0590.3450.0860.241
PPARA0.1240.0990.0690.318
SIRT10.0360.5300.1580.064
Circadian rhythmicity was assessed at the group level for NW-MSCs and Ob-MSCs. R2 represents goodness of fit to a cosine wave and p-value statistical significance of curve fit. Statistical tests were determined using CircWave. * p ≤ 0.05.
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

Erickson, M.L.; Dobias, D.; Keleher, M.R.; Dabelea, D.; Bergman, B.C.; Broussard, J.L.; Boyle, K.E. In Vitro Circadian Clock Gene Expression Assessments in Mesenchymal Stem Cells from Human Infants: A Pilot Study. Nutrients 2024, 16, 52. https://doi.org/10.3390/nu16010052

AMA Style

Erickson ML, Dobias D, Keleher MR, Dabelea D, Bergman BC, Broussard JL, Boyle KE. In Vitro Circadian Clock Gene Expression Assessments in Mesenchymal Stem Cells from Human Infants: A Pilot Study. Nutrients. 2024; 16(1):52. https://doi.org/10.3390/nu16010052

Chicago/Turabian Style

Erickson, Melissa L., Devin Dobias, Madeline Rose Keleher, Dana Dabelea, Bryan C. Bergman, Josiane L. Broussard, and Kristen E. Boyle. 2024. "In Vitro Circadian Clock Gene Expression Assessments in Mesenchymal Stem Cells from Human Infants: A Pilot Study" Nutrients 16, no. 1: 52. https://doi.org/10.3390/nu16010052

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