A prospective case-control study on miRNA circulating levels in subjects born small for gestational age (SGA) evaluated from childhood into young adulthood

Abstract

Objective

microRNAs (miRNAs) associated with metabolic risk have never been extensively investigated in SGA subjects. The aim of the current study was to evaluate miRNAs in SGA and AGA subjects and their relationships with the metabolic status and growth.

Design and methods

A prospective longitudinal case-control study was performed in 23 SGA with postnatal catch-up growth and 27 AGA subjects evaluated at the age of 9 and 21 years. Circulating levels of miR-122-5p, miR-16-5p, miR-126-3p, and miR-486-5p were assessed by qPCR.

Results

SGA subjects were shorter both at 9 and at 21 years. No significant differences in insulin like growth factors and metabolic profile were found with the exception of basal glycemia at 9 years. miRNA levels did not differ between SGA and AGA subjects, at 9 and 21 years. miR-16-5p and miR-126-3p levels were higher at 9 than at 21 years. In SGA subjects, miR-122-5p at 9 years was inversely related to adiponectin levels at 21 years and miR-486-5p at 9 years was inversely related to whole-body insulin sensitivity at 9 years and directly related to Hb1Ac at 21 years. Regression analyses showed no predictive value of miRNAs for growth parameters in neither SGA nor AGA subjects.

Conclusions

SGA with postnatal catch-up growth did not show any difference in metabolic risk markers or miRNA circulating levels compared to AGA controls in childhood and young adulthood. miR-122-5p during childhood could identify SGA subjects at higher risk of developing insulin resistance and, eventually, type 2 diabetes in adulthood but further studies are needed to confirm it.

About the Authors:

Elena Inzaghi

Roles Conceptualization, Data curation, Formal analysis, Investigation, Writing – original draft

* E-mail: [email protected]

Affiliation: Dipartimento Pediatrico Universitario Ospedaliero, “Bambino Gesù” Children’s Hospital – Tor Vergata University, Rome, Italy

ORCID logo http://orcid.org/0000-0002-8315-5168

Anna Kistner

Roles Data curation, Formal analysis, Funding acquisition, Investigation, Writing – original draft

Affiliations Department of Molecular Medicine and Surgery, Karolinska Institutet, Stockholm, Sweden, Department of Medical Radiation Physics and Nuclear Medicine, Imaging and Physiology, Karolinska University Hospital, Stockholm, Sweden

Daniela Germani

Roles Formal analysis, Investigation, Methodology

Affiliation: Dipartimento di Medicina dei Sistemi, University of Rome Tor vergata, Rome, Italy

Annalisa Deodati

Roles Formal analysis, Investigation, Writing – review & editing

Affiliation: Dipartimento Pediatrico Universitario Ospedaliero, “Bambino Gesù” Children’s Hospital – Tor Vergata University, Rome, Italy

Mireille Vanpee

Roles Data curation, Funding acquisition, Writing – review & editing

Affiliation: Department of Women’s and Children’s Health, Karolinska Institutet and University Hospital, Stockholm, Sweden

Lena Legnevall

Roles Data curation, Formal analysis, Writing – review & editing

Affiliation: Department of Women’s and Children’s Health, Karolinska Institutet and University Hospital, Stockholm, Sweden

Katarina Berinder

Roles Data curation, Formal analysis, Writing – review & editing

Affiliations Department of Molecular Medicine and Surgery, Karolinska Institutet, Stockholm, Sweden, Patient Area Endocrinology and Nephrology, Karolinska University Hospital, Stockholm, Sweden

Stefano Cianfarani

Roles Conceptualization, Formal analysis, Methodology, Project administration, Supervision, Writing – review & editing

Affiliations Dipartimento Pediatrico Universitario Ospedaliero, “Bambino Gesù” Children’s Hospital – Tor Vergata University, Rome, Italy, Department of Women’s and Children’s Health, Karolinska Institutet and University Hospital, Stockholm, Sweden

Introduction

Early life events are critical for the susceptibility to adult chronic diseases and for adult metabolic outcomes [1] and a suboptimal in utero environment leads to low birth weight and intrauterine growth restriction (IUGR), which have been widely associated with the development of cardiovascular disease, insulin resistance, type 2 diabetes (T2D), hypertension, and hyperlipidemia in adulthood [2, 3, 4, 5, 6, 7, 8]. Small for gestational age (SGA) subjects are a heterogenous population in which the pattern and timing of catch-up growth seem to play a major role in influencing the long-term metabolic risk [9, 10].

The molecular mechanisms underlying the development of IUGR-related metabolic anomalies are unknown, and the involvement of epigenetic mechanisms and the regulation of microRNA (miRNA) expression has been hypothesized [11,12].

MiRNAs are short (~22 nucleotides), non-coding single-stranded RNAs, which post-transcriptionally regulate gene expression by inhibiting translation or promoting degradation of the target messenger mRNA. MiRNAs influence many biological processes, regulate cellular function and are involved in metabolic regulation [13,14]. MiRNA circulating profile differs between physiological and pathological states. As serum levels of miRNAs are stable in clinical samples and resistant to the degradation by ribonucleases, they have been proposed as biomarkers and predictors of diseases, including cardiometabolic disease [15,16, 17, 18]. Furthermore, miRNA expression in placenta has been even associated with prenatal and postnatal growth [19]. Pregnancy-related miRNAs may be involved in pregnancy complications [20]. The intrauterine environment is an established key factor potentially affecting postnatal health and miRNAs may mediate the impact of an adverse uterine milieu on metabolic homeostasis.

We have recently described a specific profile of circulating miRNAs in obese SGA children characterized by upregulation of miR-92a-3p, miR-122-5p, miR-423-5p, miR-484, miR-486-3p and miR-532-5p, and down-regulation of miR-181b-5p [21].

The aim of this study was to assess the circulating levels of miRNAs, previously associated with metabolic alterations [22,23,24,25,26] in SGA and AGA subjects longitudinally evaluated from childhood into adulthood.

Materials and methods

Study population

This prospective longitudinal study was performed in a cohort of 23 SGA (13 females/10 males) and 27 appropriate for gestational age (AGA) (16 females/11 males) subjects, longitudinally evaluated at the age of 9 and 21 years.

The initial study population included 26 SGA children born at term and 30 AGA children born at term (control group), evaluated and described in a previous study [27]. These patients were evaluated at the Karolinska Institutet and Karolinska University hospital in Stockholm at mean age 9.8 years (26) and 24 SGA and 27 AGA agreed to participate to a follow up study after puberty. The final population of the present study consists of 23 SGA children and 27 AGA subjects, as one sample of a SGA subject was not available.

Being born SGA was defined as a birth weight less than—2 SD (i.e. below the sex-specific 2.5th percentile for gestational age) according to Swedish reference data. Ultrasound performed early during pregnancy was used to define gestational age. Standard deviation scores (SDS) for height and Body Mass Index (BMI) are based on the Swedish reference curves at 9 years and in adulthood [28].

All participating subjects were singletons born in Stockholm at Karolinska University Hospital, without chromosomal anomalies, congenital infections or life-threatening congenital anomalies. The SGA group consisted of children born from 1990 to 1993 as well as the control group which consists of healthy normal birth weight subjects born from 1990 to 1993. At follow-up, the entire study group was investigated between the 26th of April 2013 to the 1st of March 2014.

The study was approved by the Regional Ethical Review Board, Stockholm. In the prepubertal study the families were contacted and offered to participate in the study by a letter. After receiving the confirmatory response, signed by the parents and the child, the families were called by the research nurse to explain the study more in detail, to confirm participation by oral consent and to arrange a time for the study (prepubertal). In adulthood all participants gave their written informed consent.

At 9 years of age patients were admitted to the hospital after a 10- to 12-h overnight fast. Height and weight were measured and for each patient body mass index (BMI) was calculated by using the following formula: BMI = weight (kg)/height (m2) and expressed as SD.

Fasting blood samples were obtained to evaluate blood glucose, insulin, IGF-binding protein-1 (IGFBP-1), and leptin. A standard oral glucose tolerance test (OGTT) was performed in all children by giving 1.75 g ⁄ kg body weight glucose (up to 75 g). Thereafter, blood samples were taken at 30 and 120 min for measurements of blood glucose, insulin and IGFBP-1. Glucose tolerance was categorized as normal based on 2-h plasma glucose concentration during an OGTT according to reference criteria [29]. Values of homoeostasis model assessment—insulin resistance (HOMA-IR) and of a modified whole-body insulin sensitivity index (WBISI) were used as insulin resistance ⁄ sensitivity markers [30, 31]. The HOMA-IR was calculated as fasting insulin (mU⁄L) · fasting glucose (mmol ⁄L) ⁄ 22.5. WBISI was calculated as 10 000 ⁄ square root of [(fasting glucose · fasting insulin) · (mean glucose · mean insulin during OGTT, using the three values at 0, 30 and 120 min)] [31]. As glucose and insulin measurements were not taken at 60, 90,150 and 180 min, these values are missing in the WBISI calculations.

At 21 years, patients were re-evaluated. Fasting blood samples were obtained to analyse hemoglobin A1C, (HbA1C), blood glucose, insulin, IGFBP-1, lipid profile (cholesterol, triglycerides, HDL-cholesterol (HDL-c), LDL-cholesterol (LDL-c), Apolipoprotein-A1, ApolipoproteinB1), Leptin, and Adiponectin. In addition, waist circumference (WC)/height ratio, WC/hip ratio, fat %, trunk fat % were evaluated. Height was assessed by the same stable length meter. WC and hip ratio were estimated by measuring tape. Fat % and trunk fat % were evaluated by impedance (Body composition analyser BC-418MA, Tanita, Tanita Corporation of America, Inc, Arlington Heights, Illinois, USA).

Serum samples have been collected in 23 SGA and 27 AGA subjects at 9 and 21 years and have been properly stored at– 80°C and used to evaluate insulin-like growth factors (IGF-I and IGF-II) and miRNAs.

Biochemistry

HbA1c was measured by cation exchange chromatography (MonoS column) with HPLC (Bio-Rad) (Roche Diagnostics, Basel, Switzerland) (normal reference value <5.2%), 9 years. In-house RIA was used for IGFBP-1 with individual serum samples in the same assay. The RIA for IGFBP-1 was performed according to Povoa et al 1984 with intra- and inter-assays CV 3% and 10%, respectively. Blood glucose readings were performed and documented by the nursing staff using HemaCue 201`s method (Hemocue AB, Angelholm, Sweden). In children, commercial kits were used for S-insulin by ELISA (Dako, Glostrup, Denmark). The detection limit was 3 pmol/l. The intra- and interassay CV was 6.7 and 7.5% respectively. In follow-up commercial kits were used for S-insulin by electrochemiluminescence immunoassays (Roche Diagnostics GmbH, Mannheim, Germany). Detectionlimit 0.2–1000 ml/L. Reference value 2–25 ml/L.

Cholesterol and triglycerides were analysed by enzymatic method (dxc800 (2020); Beckman AB, Sychron LX, Beckman Synchron LX, Beckman Coulter Inc., Diamond Diagnostics, Holliston, MA, USA), HDL-c and LDL-c by homogenous method (dxc800 (2020); Beckman AB, Sychron LX, Beckman Synchron LX, Beckman Coulter Inc., Diamond Diagnostics, Holliston, MA, USA). Apolipoprotein-A1 and ApolipoproteinB1 levels were determined by turbidimetry assays (Beckman AB, Sychron LX, Beckman Synchron LX, Beckman Coulter Inc., Diamond Diagnostics, Holliston, MA, USA). For Apo A1, the analysis range was 0.21–3.2 g/L, the reference range for men was 1.10–1.80 g/L, and the reference range for women was 1.10–2.10 g/L. The Apo B analysis range was 0.26–3.2 g/L, and the reference range (<40 years) was 0.50–1.50 g/L.

Leptin was measured by RIA, using HL-81K (Linco Research, Inc., St. Charles, MO, USA). The intra-assay coefficient of variation (CV) for the leptin analysis was 5%. The interassay CV was 4.5%. Adiponectin was analyzed by RIA using HADP-61HK (Linco Research, Inc., St. Charles, MO, USA). The samples are diluted 1:500 before the procedure begins. Analyzing interval Ca 0,8–200 μg/L which corresponds to 0.4–100 mg/L after diluting factorisation, detection limit 1.0 ug/L, reference limit 6,0–25 mg/L. The sensitivity of the adiponectin analysis was 1 mg/l. The intra- and interassay CV were 4 and 8.5% respectively.

Assays

IGF-I and IGF-II levels have been evaluated on serum samples of all patients at both 9 and 21 years. IGFs have been measured as previously described [32], by using specific immunoassay kit. IGF-I and IGF-II values were expressed as SDS relative to a reference cohort.

MicroRNA measurement

Circulating RNA was isolated from a serum volume of 100 μl using a miRNeasy Mini Kit according to the manufacturer’s instruction. (Qiagen, Dusseldorf, Germany). Synthetic cel-mir-39 (Applied biosystems, Thermo Fisher scientific, Waltham, MA USA) was used as a control spike-in to help monitor RNA recovery and reverse transcription efficiency. Total RNA was eluted in 14 μl of nuclease free H2O. A fixed volume of 2 μl of RNA eluate was used as input in each reverse transcription reaction using Taqman ADV MicroRNA cDNA synthesis (Applied biosystems, Thermo fisher scientific, Waltham, MA USA) according to the company’s recommendations instructions. Quantitative real-time polymerase chain reaction (qPCR) for miR-486-5p, miR-122-5p, miR-126-5p and miR-16-3p and for miR-484, miR-191, miR-93 and miR-24 were performed by using the Taqman technology (Applied biosystems). A threshold value of 0.2 was set to ensure comparability across plates and a control sample was evaluated in each plates to verify reproducibility. Ct values greater than 32 were excluded. Results were normalized with the mean normalized method [33], which uses the mean of miR-93, miR-191, miR-24 and miR-484 as normalization factor. MiRNA expression levels were calculated according to the ΔCT method.

Statistics

Results are reported as the mean ± SD. Normal distribution was assessed by Shapiro-Wilk test. Differences between means were assessed using unpaired two-tailed t test and one-way ANOVA. The relationships among parameters were evaluated by Pearson correlation. Multiple regression and stepwise regression analyses were used in the selection of predictors of adult height. Significance was assigned for P < 0.05. All analyses were performed using SPSS version 21.0 for windows (SPSS, Chicago, Illinois).

Results

Anthropometry

SGA subjects, despite having experienced a catch-up growth leading to normal stature, overall remained shorter than AGA subjects, both at 9 (0.08 ± 1.06 vs 0.76 ± 1.2 SDS, p<0.05) and at 21 years (-0.21 ± 0.76 vs 0.65± 1.32 SDS, p<0.05). BMI was not significantly different between SGA and AGA subjects, both at 9 and at 21 years. At 21 years, SGA subjects showed no significant differences in WC/height ratio, WC/hip, total body fat percentage, and trunk fat percentage.

Biochemistry

Biochemical and metabolic profile were overall not significantly different between SGA and AGA, both at 9 and at 21 years (Tables 1 and 2). Baseline glucose at 9 years, though in the normal range in all children, was significantly higher in SGA subjects (4.62 ± 0.48 vs 4.18 ± 0.53 mmol/L, p<0.01). At 21 years glycemia was not significantly different (4.77 ± 0.29 vs 4.62 ± 0.44 mmol/L, p = 0.16). Insulin was similar between the groups both at 9 (6.17 ± 2.86 vs 8.32 ± 11.19mU/L, p = 0.27) and 21 years (9.07 ± 3.72 vs 9.07 ± 6.64 p = 0.3), as well as HOMA-IR at 9 years (1.27 ± 0.61 vs 1.53 ± 2.03 p = 0.12) and 21 years (1.95 ± 0.86 vs 1.81 ± 1.37 p = 0.23). No significant differences were found in leptin, adiponectin, lipid profile, apolipoprotein A1 and apolipoprotein B1 both at 9 and 21 years.

[Figure omitted. See PDF.]

Table 1. Biochemical and metabolic profile in small for gestational age (SGA) and appropriate for gestational age (AGA) subjects at 9 years.

https://doi.org/10.1371/journal.pone.0228075.t001

[Figure omitted. See PDF.]

Table 2. Biochemical and metabolic profile between small for gestational age (SGA) and appropriate for gestational age (AGA) subjects at 21 years.

https://doi.org/10.1371/journal.pone.0228075.t002

No significant differences were found in IGF-I, IGF-II and IGFBP-1 levels both at 9 and 21 years.

MiRNAs

MiR-122-5p, miR-16-5p, miR-486-5p and miR-126-3p serum levels were evaluated in all SGA and AGA subjects. MiRNA expression was not influenced by gender. MiR-122-5p and miR-486-5p expression was not influenced by the age (Fig 1a and 1b), whereas mir-16-5p and mir-126-3p levels were higher at 9 than 21 years (Fig 1c and 1d), both in the whole study population and in SGA and AGA groups.

[Figure omitted. See PDF.]

Fig 1. (a-d). MiRNA expression according to age.

MiR-122-5p and miR-486-5p levels were not influenced by the age (p > 0.05), while miR-16-5p and miR-126-3p levels were higher at 9 years than at 21 years both in SGA and in AGA subjects (p<0.05). SGA: small for gestational age. AGA: appropriate for gestational age.

https://doi.org/10.1371/journal.pone.0228075.g001

By comparing miRNAs expression between SGA and AGA groups, no significant difference was found in miRNA circulating levels both at 9 (Fig 2) and 21 years (Fig 3).

[Figure omitted. See PDF.]

Fig 2. MiRNA expression in SGA and AGA subjects.

Comparative analysis of miRNA expression between SGA and AGA subjects at 9 years. SGA: small for gestational age. AGA: appropriate for gestational age. MiR-122: miR-122-5p. MiR-486: miR-486-5p. MiR-16: miR-16-5p. MiR-126: miR-126-3p.

https://doi.org/10.1371/journal.pone.0228075.g002

[Figure omitted. See PDF.]

Fig 3. MiRNA expression in SGA and AGA subjects.

Comparative analysis of miRNA expression between SGA and AGA subjects at 21 years. SGA: small for gestational age. AGA: appropriate for gestational age. MiR-122: miR-122-5p. MiR-486: miR-486-5p. MiR-16: miR-16-5p. MiR-126: miR-126-3p.

https://doi.org/10.1371/journal.pone.0228075.g003

In SGA subjects, mir-122-5p expression at 9 years was inversely related to adiponectin levels at 21 years (r = -0.48, p<0.05) (Fig 4). MiR-122-5p expression was not statistically significantly related with Adiponectin levels at 21 years (r = -0.479, p = 0.071). Mir-486-5p expression at 9 years was inversely related to WBISI at 9 years (r = -0.52, p<0.05) and directly related to Hb1Ac at 21 years (r = 0.52, p<0.05) (Fig 5a and 5b). These associations were not found in AGA subjects.”

[Figure omitted. See PDF.]

Fig 4. MiR-122 expression and Adiponectin levels.

In SGA subjects, miR-122-5p expression at 9 years resulted inversely related to Adiponectin levels at 21 years. SGA: small for gestational age.

https://doi.org/10.1371/journal.pone.0228075.g004

[Figure omitted. See PDF.]

Fig 5. (a and b). MiR-486 expression and WBISI and HbA1c.

Mir-486-5p expression at 9 years was inversely related to WBISI at 9 years (5a) and directly related to Hb1Ac at 21 years (5b).

https://doi.org/10.1371/journal.pone.0228075.g005

Multiple regression analysis

Stepwise regression analysis showed that no miRNA was predictive of insulin sensitivity (HOMA-IR) at 21 years in both SGA and AGA subjects. The only significant predictor of adiponectin levels at the age of 21 years was mir-122-5p at the age of 9 (β = 0.59, adjusted r2 = 0.30, F = 7.43, p<0.02) in SGA subjects, and weight gain between birth and 9 years (β = 0.53, adjusted r2 = 0.24, F = 6.0, p<0.03) in AGA subjects. miRNAs were not significant predictors of growth in SGA or AGA subjects. In AGA subjects the best predictors of height at the age of 21 years were father’s height (β = 0.47), birth length (β = 0.47) and mother’s height (β = 0.34), with an adjusted r2 = 0.83, F = 26.8, p<0.0001. In SGA children the only significant predictor of height at 21 years was mother’s height (β = 0.60, adjusted r2 = 0.32, F = 8.16, p = 0.01).

Discussion

A suboptimal intrauterine environment leads to low birth weight which has been associated with an increased risk of cardiometabolic disease in adult life [2, 3, 4, 5, 6, 7, 34]. Children born SGA comprise a heterogeneous population with different metabolic and growth outcomes depending on the underlying pathophysiology. Unfortunately, to date, no tools are available to identify subjects deserving a more careful metabolic and growth monitoring and a tailored management.

Although circulating miRNAs have been associated with various metabolic alterations, to our knowledge, this is the first prospective longitudinal case-control study from childhood into early adulthood assessing miRNA profile in SGA subjects. Our results show that SGA subjects who have experienced catch-up growth and are metabolically healthy have circulating levels of specific candidate miRNAs such as miR-122-5p, miR-16-5p, miR-126-3p, and miR-486-5p not different from AGA subjects both in childhood and in early adulthood.

MiR-122 has been widely investigated for its role in metabolism [22]. In adults, circulating miR-122 levels have been associated with the severity of NAFLD, lipid metabolism, obesity, insulin resistance, T2D, and metabolic syndrome [35, 36, 37]. In obese children, serum levels of miR-122 represent an extrahepatic fingerprint of NAFLD [38]. In this study we have shown that the expression of miR-122-5p in childhood is the major predictor of adiponectin levels in early adulthood in SGA subjects. MiR-122 has been correlated with insulin resistance and T2D [39]. As adiponectin levels decrease in case of insulin resistance and are reduced prior to the development of T2D [40], it could be speculated that miR-122-5p levels during childhood could identify SGA subjects at higher risk of developing insulin resistance and, eventually, T2D.

MiR-486-5p is involved in glucose metabolism and insulin signaling [41]. MiR-486-5p increases in obese children [24] and adult patients with T2D and miR-486 has been implicated in accelerating preadipocyte proliferation and myotube glucose intolerance [23]. MiR-486-5p has also been associated with insulin resistance and response to thiazolidenidones [41]. Finally, miR-486 modulates PI3K/Akt signaling by targeting the phosphatase and tensin homolog (PTEN) and the forkhead box O1a (FOXO1A) [42]. Therefore, by modulating these targets, miR-486 may affect insulin sensitivity and blood glucose. Our data support an association between miR-486-5p and glucose homeostasis in SGA subjects as we found an inverse relationship between miR-486-5p expression and index of insulin sensitivity and Hb1Ac levels.

MiR-16 affects insulin signaling and is down-regulated in the skeletal muscle of twins with T2D compared with their non-diabetic identical twins [25]. Placental miR-16 expression predicts the risk of being born SGA [26] and miR-16-5p is down-regulated in IUGR pregnancies [43]. MiR-126 expression is down-regulated in patients with coronary artery disease [44] and miR-126-3p is down-regulated in pregnancies complicated by IUGR [43]. Our results show no significant relationship of miR-16-5p and mir-126-3p with metabolic and growth parameters in both SGA and AGA subjects.

Our data indicate an age dependency of both miR-16-5p and miR-126-3p. This finding is in agreement with the recent observation that some miRNAs are differentially expressed according to age [45] being usually downregulated with increasing age.

During the first years of life most SGA subjects experience a spontaneous catch-up growth and achieve normal height although approximately 10% remain permanently short [46]. All SGA subjects investigated in this study had experienced catch-up growth and achieved a normal height in both childhood and adulthood, though being shorter than AGA subjects.

The metabolic profile of SGA and AGA subjects was essentially similar. The only different parameter, was basal glycemia at 9 years. This is in contrast with previous findings showing that SGA subjects have a worse metabolic profile, especially after experiencing an early catch-up growth [47]. Increased fat mass and abdominal fat have been reported in SGA subjects compared to AGA subjects [48]. Notably, early weight gain rather than birth size influences body composition [49]. Our results show no significant differences in body composition between SGA and AGA subjects. This finding contrasting previous reports may be explained either by the small study cohort or by a different timing of catch-up growth, as we have no information regarding early postnatal growth rate.

IGF system has been associated with metabolism and metabolic parameters. Specifically, IGF-I levels have been related with cardiometabolic risk markers both in children and adults [50], and IGF2 polymorphisms have been related to body weight, insulin sensitivity, fat-free mass, lipid profile, and blood pressure [51]. The results of this study show no differences in IGF system as well as no correlation of IGFs with both anthropometric and metabolic parameters.

We are aware of some limitations of our study. First, in SGA subjects the metabolic risk is deeply influenced by the pattern of early catch-up growth, the faster the early catch-up growth the worse is the impact on cardiometabolic risk [8]. Unfortunately, no information on early growth was available in our subjects. Second, the relatively small sample size may have affected the results. Third, we focused the analysis on candidate mi-RNAs and a complete miRNoma analysis was not performed. Fourth, all SGA subjects investigated in this study had experienced catch-up growth, for this reason our results cannot be applied to SGA subjects without postnatal catch-up growth. Finally, as metabolic alterations may require time to occur, a longer longitudinal study would be more appropriate to test whether miRNAs are predictive of the development of metabolic diseases in adulthood and can be used as early markers of an increased cardiometabolic risk in SGA subjects. The major strengths of this study are the longitudinal prospective design and the original approach applying for the first time miRNA serum profile analysis to SGA subjects.

In conclusion, our data do not support the speculation that SGA subjects have a different/peculiar profile of circulating mi-RNAs related to cardiometabolic risk. Further prospective and more prolonged studies in larger populations of SGA subjects including subjects with both catch-up and no catch-up growth are needed to assess whether miRNAs may represent early biomarkers of long-term increased cardiometabolic risk.

Supporting information

[Figure omitted. See PDF.]

S1 Dataset. Includes Ct values of miRNAs in SGA and AGA subjects.

https://doi.org/10.1371/journal.pone.0228075.s001

(XLS)

Acknowledgments

We are grateful to the late Professor Kerstin Hall for helpful discussion and valuable suggestions.

Citation: Inzaghi E, Kistner A, Germani D, Deodati A, Vanpee M, Legnevall L, et al. (2020) A prospective case-control study on miRNA circulating levels in subjects born small for gestational age (SGA) evaluated from childhood into young adulthood. PLoS ONE 15(1): e0228075. https://doi.org/10.1371/journal.pone.0228075

References

1. Langley-Evans SC, McMullen S. Developmental origins of adult disease. Med Princ Pract. 2010;19(2):87–98. Epub 2010/02/06. pmid:20134170.

2. Barker DJ, Winter PD, Osmond C, Margetts B, Simmonds SJ. Weight in infancy and death from ischaemic heart disease. Lancet. 1989;2(8663):577–80. Epub 1989/09/09. pmid:2570282.

3. Barker DJ, Hales CN, Fall CH, Osmond C, Phipps K, Clark PM. Type 2 (non-insulin-dependent) diabetes mellitus, hypertension and hyperlipidaemia (syndrome X): relation to reduced fetal growth. Diabetologia. 1993;36(1):62–7. Epub 1993/01/01. pmid:8436255.

4. Ravelli AC, van der Meulen JH, Michels RP, Osmond C, Barker DJ, Hales CN, et al. Glucose tolerance in adults after prenatal exposure to famine. Lancet. 1998;351(9097):173–7. Epub 1998/02/05. pmid:9449872.

5. Barker DJ, Martyn CN, Osmond C, Hales CN, Fall CH. Growth in utero and serum cholesterol concentrations in adult life. Bmj. 1993;307(6918):1524–7. Epub 1993/12/11. pmid:8274921.

6. Ramadhani MK, Grobbee DE, Bots ML, Castro Cabezas M, Vos LE, Oren A, et al. Lower birth weight predicts metabolic syndrome in young adults: the Atherosclerosis Risk in Young Adults (ARYA)-study. Atherosclerosis. 2006;184(1):21–7. Epub 2005/12/06. pmid:16326169.

7. van der Steen M, Hokken-Koelega AC. Growth and Metabolism in Children Born Small for Gestational Age. Endocrinol Metab Clin North Am. 2016;45(2):283–94. Epub 2016/06/01. pmid:27241965.

8. Mericq V, Martinez-Aguayo A, Uauy R, Iñiguez G, Van der Steen M, Hokken-Koelega A. Long-term metabolic risk among children born premature or small for gestational age. Nat Rev Endocrinol. 2017;13(1):50–62. Epub 2016/08/19. pmid:27539244.

9. Cianfarani S, Germani D, Branca F. Low birthweight and adult insulin resistance: the "catch-up growth" hypothesis. Arch Dis Child Fetal Neonatal Ed. 1999;81(1):F71–3. Epub 1999/06/22. pmid:10375369.

10. Leunissen RW, Kerkhof GF, Stijnen T, Hokken-Koelega A. Timing and tempo of first-year rapid growth in relation to cardiovascular and metabolic risk profile in early adulthood. JAMA: the journal of the American Medical Association. 2009;301(21):2234–42. pmid:19491185

11. Mansego ML, Garcia-Lacarte M, Milagro FI, Marti A, Martinez JA. DNA methylation of miRNA coding sequences putatively associated with childhood obesity. Pediatr Obes. 2017;12(1):19–27. Epub 2016/01/20. pmid:26780939.

12. Sookoian S, Gianotti TF, Burgueno AL, Pirola CJ. Fetal metabolic programming and epigenetic modifications: a systems biology approach. Pediatr Res. 2013;73(4 Pt 2):531–42. Epub 2013/01/15. pmid:23314294.

13. Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15(8):509–24. Epub 2014/07/17. pmid:25027649.

14. Fischer-Posovszky P, Roos J, Kotnik P, Battelino T, Inzaghi E, Nobili V, et al. Functional Significance and Predictive Value of MicroRNAs in Pediatric Obesity: Tiny Molecules with Huge Impact? Horm Res Paediatr. 2016;86(1):3–10. Epub 2016/05/11. pmid:27161162.

15. Seyhan AA. microRNAs with different functions and roles in disease development and as potential biomarkers of diabetes: progress and challenges. Mol Biosyst. 2015;11(5):1217–34. Epub 2015/03/15. pmid:25765998.

16. Kitsiou-Tzeli S, Tzetis M. Maternal epigenetics and fetal and neonatal growth. Curr Opin Endocrinol Diabetes Obes. 2017;24(1):43–6. Epub 2016/11/30. pmid:27898587.

17. Romaine SP, Tomaszewski M, Condorelli G, Samani NJ. MicroRNAs in cardiovascular disease: an introduction for clinicians. Heart. 2015;101(12):921–8. Epub 2015/03/31. pmid:25814653.

18. Parrizas M, Novials A. Circulating microRNAs as biomarkers for metabolic disease. Best Pract Res Clin Endocrinol Metab. 2016;30(5):591–601. Epub 2016/12/08. pmid:27923453.

19. Carreras-Badosa G, Bonmati A, Ortega FJ, Mercader JM, Guindo-Martinez M, Torrents D, et al. Dysregulation of Placental miRNA in Maternal Obesity Is Associated With Pre- and Postnatal Growth. J Clin Endocrinol Metab. 2017;102(7):2584–94. Epub 2017/04/04. pmid:28368446.

20. Cai M, Kolluru GK, Ahmed A. Small Molecule, Big Prospects: MicroRNA in Pregnancy and Its Complications. J Pregnancy. 2017;2017:6972732. Epub 2017/07/18. pmid:28713594.

21. Marzano F, Faienza MF, Caratozzolo MF, Brunetti G, Chiara M, Horner DS, et al. Pilot study on circulating miRNA signature in children with obesity born small for gestational age and appropriate for gestational age. Pediatr Obes. 2018;13(12):803–11. Epub 2018/08/31. pmid:30160046.

22. Rottiers V, Naar AM. MicroRNAs in metabolism and metabolic disorders. Nat Rev Mol Cell Biol. 2012;13(4):239–50. Epub 2012/03/23. pmid:22436747.

23. Cui X, You L, Zhu L, Wang X, Zhou Y, Li Y, et al. Change in circulating microRNA profile of obese children indicates future risk of adult diabetes. Metabolism. 2018;78:95–105. Epub 2017/10/03. pmid:28966078.

24. Prats-Puig A, Ortega FJ, Mercader JM, Moreno-Navarrete JM, Moreno M, Bonet N, et al. Changes in circulating microRNAs are associated with childhood obesity. J Clin Endocrinol Metab. 2013;98(10):E1655–60. Epub 2013/08/10. pmid:23928666.

25. Bork-Jensen J, Scheele C, Christophersen DV, Nilsson E, Friedrichsen M, Fernandez-Twinn DS, et al. Glucose tolerance is associated with differential expression of microRNAs in skeletal muscle: results from studies of twins with and without type 2 diabetes. Diabetologia. 2015;58(2):363–73. Epub 2014/11/19. pmid:25403480.

26. Maccani MA, Padbury JF, Marsit CJ. miR-16 and miR-21 expression in the placenta is associated with fetal growth. PLoS One. 2011;6(6):e21210. Epub 2011/06/24. pmid:21698265.

27. Kistner A, Rakow A, Legnevall L, Marchini G, Brismar K, Hall K, et al. Differences in insulin resistance markers between children born small for gestational age or born preterm appropriate for gestational age. Acta Paediatr. 2012;101(12):1217–24. Epub 2012/08/29. pmid:22924816.

28. Wikland KA, Luo ZC, Niklasson A, Karlberg J. Swedish population-based longitudinal reference values from birth to 18 years of age for height, weight and head circumference. Acta Paediatr. 2002;91(7):739–54. Epub 2002/08/31. pmid:12200898.

29. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2010;33 Suppl 1:S62–9. Epub 2010/01/29. pmid:20042775.

30. Yeckel CW, Weiss R, Dziura J, Taksali SE, Dufour S, Burgert TS, et al. Validation of insulin sensitivity indices from oral glucose tolerance test parameters in obese children and adolescents. J Clin Endocrinol Metab. 2004;89(3):1096–101. Epub 2004/03/06. pmid:15001593.

31. Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care. 1999;22(9):1462–70. Epub 1999/09/10. pmid:10480510.

32. Cianfarani S, Inzaghi E, Alisi A, Germani D, Puglianiello A, Nobili V. Insulin-like growth factor-I and -II levels are associated with the progression of nonalcoholic fatty liver disease in obese children. J Pediatr. 2014;165(1):92–8. Epub 2014/03/13. pmid:24607243.

33. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):Research0034. Epub 2002/08/20. pmid:12184808.

34. Mericq V, Martinez-Aguayo A, Uauy R, Iniguez G, Van der Steen M, Hokken-Koelega A. Long-term metabolic risk among children born premature or small for gestational age. Nat Rev Endocrinol. 2017;13(1):50–62. Epub 2016/11/04. pmid:27539244.

35. Wang R, Hong J, Cao Y, Shi J, Gu W, Ning G, et al. Elevated circulating microRNA-122 is associated with obesity and insulin resistance in young adults. Eur J Endocrinol. 2015;172(3):291–300. Epub 2014/12/18. pmid:25515554.

36. Miyaaki H, Ichikawa T, Kamo Y, Taura N, Honda T, Shibata H, et al. Significance of serum and hepatic microRNA-122 levels in patients with non-alcoholic fatty liver disease. Liver Int. 2014;34(7):e302–7. Epub 2013/12/10. pmid:24313922.

37. Yamada H, Suzuki K, Ichino N, Ando Y, Sawada A, Osakabe K, et al. Associations between circulating microRNAs (miR-21, miR-34a, miR-122 and miR-451) and non-alcoholic fatty liver. Clin Chim Acta. 2013;424:99–103. Epub 2013/06/04. pmid:23727030.

38. Brandt S, Roos J, Inzaghi E, Kotnik P, Kovac J, Battelino T, et al. Circulating levels of miR-122 and nonalcoholic fatty liver disease in pre-pubertal obese children. Pediatr Obes. 2018;13(3):175–82. Epub 2017/12/23. pmid:29271122.

39. Castano C, Kalko S, Novials A, Parrizas M. Obesity-associated exosomal miRNAs modulate glucose and lipid metabolism in mice. Proc Natl Acad Sci U S A. 2018;115(48):12158–63. Epub 2018/11/16. pmid:30429322.

40. Lihn AS, Pedersen SB, Richelsen B. Adiponectin: action, regulation and association to insulin sensitivity. Obes Rev. 2005;6(1):13–21. Epub 2005/01/19. pmid:15655035.

41. Flowers E, Aouizerat BE, Abbasi F, Lamendola C, Grove KM, Fukuoka Y, et al. Circulating microRNA-320a and microRNA-486 predict thiazolidinedione response: Moving towards precision health for diabetes prevention. Metabolism. 2015;64(9):1051–9. Epub 2015/06/03. pmid:26031505.

42. Small EM, O’Rourke JR, Moresi V, Sutherland LB, McAnally J, Gerard RD, et al. Regulation of PI3-kinase/Akt signaling by muscle-enriched microRNA-486. Proc Natl Acad Sci U S A. 2010;107(9):4218–23. Epub 2010/02/10. pmid:20142475.

43. Hromadnikova I, Kotlabova K, Hympanova L, Krofta L. Cardiovascular and Cerebrovascular Disease Associated microRNAs Are Dysregulated in Placental Tissues Affected with Gestational Hypertension, Preeclampsia and Intrauterine Growth Restriction. PLoS One. 2015;10(9):e0138383. Epub 2015/09/24. pmid:26394310.

44. Wang X, Lian Y, Wen X, Guo J, Wang Z, Jiang S, et al. Expression of miR-126 and its potential function in coronary artery disease. Afr Health Sci. 2017;17(2):474–80. Epub 2017/10/25. pmid:29062343.

45. Huan T, Chen G, Liu C, Bhattacharya A, Rong J, Chen BH, et al. Age-associated microRNA expression in human peripheral blood is associated with all-cause mortality and age-related traits. Aging Cell. 2018;17(1). Epub 2017/10/19. pmid:29044988.

46. Clayton PE, Cianfarani S, Czernichow P, Johannsson G, Rapaport R, Rogol A. Management of the child born small for gestational age through to adulthood: a consensus statement of the International Societies of Pediatric Endocrinology and the Growth Hormone Research Society. J Clin Endocrinol Metab. 2007;92(3):804–10. Epub 2007/01/04. pmid:17200164.

47. Finken MJJ, van der Steen M, Smeets CCJ, Walenkamp MJE, de Bruin C, Hokken-Koelega ACS, et al. Children Born Small for Gestational Age: Differential Diagnosis, Molecular Genetic Evaluation, and Implications. Endocr Rev. 2018;39(6):851–94. Epub 2018/07/10. pmid:29982551.

48. Meas T, Deghmoun S, Armoogum P, Alberti C, Levy-Marchal C. Consequences of being born small for gestational age on body composition: an 8-year follow-up study. J Clin Endocrinol Metab. 2008;93(10):3804–9. Epub 2008/07/17. pmid:18628518.

49. Leunissen RW, Stijnen T, Hokken-Koelega AC. Influence of birth size on body composition in early adulthood: the programming factors for growth and metabolism (PROGRAM)-study. Clin Endocrinol (Oxf). 2009;70(2):245–51. Epub 2008/07/12. pmid:18616715.

50. Inzaghi E, Baldini Ferroli B, Fintini D, Grossi A, Nobili V, Cianfarani S. Insulin-Like Growth Factors and Metabolic Syndrome in Obese Children. Horm Res Paediatr. 2017;87(6):400–4. Epub 2017/06/02. pmid:28571015.

51. Cianfarani S. Insulin-like growth factor-II: new roles for an old actor. Front Endocrinol (Lausanne). 2012;3:118. Epub 2012/10/13. pmid:23060858.