Skip to content


  • Original investigation
  • Open Access

Subclinical inflammation associated with prolonged TIMP-1 upregulation and arterial stiffness after gestational diabetes mellitus: a hospital-based cohort study

Cardiovascular Diabetology201716:49

  • Received: 16 January 2017
  • Accepted: 4 April 2017
  • Published:



Gestational diabetes mellitus (GDM) has significant implications for the future health of the mother. Some clinical studies have suggested subclinical inflammation and vascular dysfunction after GDM. We aimed to study whether concentrations of high-sensitivity C-reactive protein (hsCRP), tissue inhibitor of metalloproteinase-1 (TIMP-1), matrix metalloproteinase-8 (MMP-8) and -9, as well as values of arterial stiffness differ between women with and without a history of GDM a few years after delivery. We also investigated possible effects of obesity on the results.


We studied two cohorts—120 women with a history of GDM and 120 controls—on average 3.7 years after delivery. Serum concentrations of hsCRP were determined by immunonephelometric and immunoturbidimetric methods, MMP-8 by immunofluorometric assay, and MMP-9 and TIMP-1 by enzyme-linked immunosorbent assays. Pulse wave velocity (PWV) was determined using the foot-to-foot velocity method from carotid and femoral waveforms by using a SphygmoCor device. Arterial compliance was measured non-invasively by an HDI/PulseWave™CR-2000 arterial tonometer. All 240 women were also included in subgroup analyses to study the effect of obesity on the results. Multiple linear regression analyses were performed with adjustment for confounding factors.


PWV after pregnancy complicated by GDM was significantly higher than after normal pregnancy, 6.44 ± 0.83 (SD) vs. 6.17 ± 0.74 m/s (p = 0.009). Previous GDM was also one of the significant determinants of PWV in multiple linear regression analyses. On the other hand, compliance indices of both large (p = 0.092) and small (p = 0.681) arteries did not differ between the study cohorts. Serum TIMP-1 levels were significantly increased after previous GDM (p = 0.020). However, no differences were found in the serum levels of MMP-8, MMP-9 or hsCRP. In subgroup analyses, there were significantly higher concentrations of hsCRP (p = 0.015) and higher PWV (p < 0.001) among obese women compared with non-obese ones.


PWV values were significantly higher after GDM compared with normoglycemic pregnancies and were associated with prolonged TIMP-1 upregulation. Cardiovascular risk factors were more common in participants with high BMI than in those with previous GDM.


  • Arterial compliance
  • Gestational diabetes mellitus
  • High-sensitivity C-reactive protein
  • Matrix metalloproteinase-8
  • Matrix metalloproteinase-9
  • Pulse wave velocity
  • Subclinical inflammation
  • Tissue inhibitor of matrix metalloproteinase-1


In developed countries, the prevalence of gestational diabetes mellitus (GDM) has increased rapidly in recent decades, along with increasing rates of obesity [1, 2]. In Finland, GDM complicated 15.9% of pregnancies in 2015 [2]. A diagnosis of GDM has significant implications for the future health of the mother. For instance, GDM has been shown to be associated with postpartum insulin resistance, hypertension, and dyslipidemia [35], placing affected women at risk of metabolic syndrome (MetS), type 2 diabetes mellitus (T2DM) and/or cardiovascular disease (CVD) later in life [58]. Incidence of CVD events, and specifically those of coronary artery disease, is known to be increased in women with previous GDM, even in the absence of T2DM [8]. Clinical studies have also revealed subclinical inflammation and vascular dysfunction after GDM [4].

High-sensitivity C-reactive protein (hsCRP) is a well-known acute-phase protein and a sensitive biomarker of systemic inflammation. Elevated levels of hsCRP are a significant risk factor for atherosclerosis [9]. The group of matrix metalloproteinases (MMPs) comprises over 20 structurally and functionally related but genetically distinct members [10, 11]. Expression and activity are normally low, but increased in many pathophysiological conditions. MMPs can modulate immunological responses, and MMPs can be either defensive or destructive [11]. Both upregulation and down-regulation of MMP-8 and -9 have been associated with several noninfectious as well infectious inflammatory states [1218]. MMP-8 may also regulate blood pressure [19]. MMPs and their inhibitors, tissue inhibitors of MMPs (TIMPs) have been related to atherosclerosis development and progression in humans [2022]. It has been suggested that imbalanced concentrations of MMP family members and TIMPs eventually exert an important role in cardiovascular risk [2125].

Inflammation may be pathogenic, by inducing vascular dysfunction [4, 26]. Arterial stiffness has proven to be an important parameter for the assessment of cardiovascular risk, and it has earlier been associated with endothelial dysfunction [27, 28]. Carotid to femoral pulse wave velocity (PWV) has emerged as the gold standard to assess arterial stiffness [29]. When the arteries are stiff or less distensible, PWV increases [30, 31]. PWV increases proportionally to the number of cardiovascular risk factors present, such as diabetes or MetS [27, 32, 33]. In epidemiological studies, increased PWV has been predictive of cardiovascular events [29].

Recently, the implications of GDM as regards women’s future health have been widely discussed. As the prevalence of GDM has increased over the years, a better understanding of the connections between previous GDM and both subclinical inflammation and vascular dysfunction would be of great benefit. In addition, recently it has been suggested that MMP-8 is associated with insulin receptor degradation, and high serum MMP-8 levels with an increased risk of diabetes mellitus type II [17]. In previous studies serum levels of MMP-8, -9, TIMP-1 and hsCRP have been shown to be biomarkers reflecting low-grade inflammation [11, 23, 24, 34, 35]. In addition, TIMP-1 has been shown to exert MMP-independent actions such as pro-inflammatory and growth-factor-like properties [3638].

With this background our aim was to define whether or not cardiovascular risk, assessed by serum concentrations of hsCRP, MMP-8, MMP-9 and TIMP-1, and values of arterial compliance and PWV are enhanced already a few years after GDM. We also evaluated the effect of obesity on the results.


In this follow-up study of two cohorts, a total of 120 women with a history of GDM during the index pregnancy were compared with 120 age-matched women with normal glucose metabolism during pregnancy. The time from the index pregnancy to the follow-up study was also matched between the study groups. All participants had delivered on average 3.7 (range 2–6) years earlier at Kanta-Häme Central Hospital, Finland, i.e. after the publication of Finnish Current Guidelines for screening GDM. Our national guidelines were published in 2008 and updated in 2013 without any change in the diagnostic criteria of GDM [39]. The complete inclusion and exclusion criteria, with power analysis, have been described earlier [40]. Briefly, GDM was defined (using the diagnostic criteria of Finnish Current Guidelines) as a pathological value in a 2-h 75-g oral glucose tolerance test (OGTT) during pregnancy: venous plasma glucose ≥5.3 mmol/L when fasting, ≥10.0 mmol/L at 1 h or ≥8.6 mmol/L at 2 h [39]. Our national diagnostic thresholds for GDM are similar to those of the International Association of Diabetes and Pregnancy Study Groups (IADPSG): plasma glucose ≥5.1 mmol/L when fasting, ≥10.0 mmol/L at 1 h or ≥8.5 mmol/L at 2 h [41]. Only singleton pregnancies were included. Women were excluded if they had type 1 or type 2 diabetes before the pregnancy, if they were pregnant at time of the study, if they had suspected or verified malignant or endocrine disease, if there was substance abuse or treatment, or a known clinical history of psychiatric illness. Controls had to have normal OGTT results during pregnancy. If the controls had experienced GDM in an earlier pregnancy, or the weight of the newborn was ≥4.5 kg, they were excluded. The electronic database of the hospital was used to pick up the cases and controls. Both recruitment and examinations were accomplished between August 2011 and July 2014.

We interviewed the participants as regards their lifestyle habits. Lifetime tobacco exposure was estimated as pack-years, and one pack-year was defined as 20 cigarettes smoked every day for 1 year [42]. Further, we interviewed the participants as regards their history of trauma or infectious diseases during the previous month. We measured resting heart rate, brachial blood pressure, weight (kg) and height (cm) of the participants, and calculated body mass index (BMI): weight in kilograms divided by height in meters squared (kg/m2).

The study was conducted in accordance with the ethical principles outlined in the Declaration of Helsinki [43], and the protocol was approved by the Ethics Committee of Kanta-Häme Hospital District (reference number 521/2010; date of approval 21.12.2010). Every participant was given both oral and written information on the study before she signed an informed consent document.

Laboratory methods

Serum samples were collected after at least 12 h of fasting and stored at −80 °C until analyzed. Serum concentrations of hsCRP were analyzed according to validated immunonephelometric (United Medix Laboratories Ltd., Espoo, Finland) and immunoturbidimetric (VITA Healthcare Services Ltd., Vita Laboratory, Helsinki, Finland) methods [44, 45]. Concentrations of MMP-8 were determined by immunofluorometric assay (IFMA) (Medix Biochemica, Espoo, Finland), as previously described [25]. Serum levels of MMP-9 and TIMP-1 were analyzed by enzyme-linked immunosorbent assay (ELISA) using commercial kits (Biotrak ELISA System; Amersham Biosciences, GE Healthcare, Buckinghamshire, UK) and according to the manufacturer’s instructions [18]. Fasting serum levels of total cholesterol (TC) and insulin were analyzed according to validated methods as described in detail earlier [40].

Determination of arterial compliance and pulse wave velocity

Three experienced nurses measured the compliance of large and small arteries after at least 10 min of rest in a semi-sitting position. The recording was carried out after an overnight fast. The participants were asked to refrain from eating, having caffeinated drinks, smoking and taking medication for 12 h, and drinking alcohol for 2 days prior to measurement. Radial artery pulse waves were recorded non-invasively with an arterial tonometer (HDI/PulseWave™CR-2000, Hypertension Diagnostics, Inc., Eagan, Minnesota, USA) and the procedure involves the use of a modified Windkessel pulse-contour method [46]. Blood volume inertia and systemic vascular resistance are used to analyze arterial compliance. The capacitive compliance of large arteries (C1), including the aorta, and the endothelial function of small arteries (C2) were automatically assessed as a mean of the five most similar pulse waves appearing during 30-s of measurement. Three consecutive measurements were performed to obtain mean results for every participant.

Carotid-femoral PWV was measured using the foot-to-foot velocity method from carotid and femoral waveforms by employing a SphygmoCor device (AtCor Medical, Sydney, Australia). Transcutaneous readings were obtained at the right common carotid artery and the right femoral artery with the subjects in a supine position with direct-contact pulse sensors. The time delay (Dt or transit time) of the two waveforms was registered, and the distance (D) between carotid and femoral recording sites was obtained by subtracting the carotid measurement site to sternal notch distance from the sternal notch to the femoral measurement site distance. PWV was calculated as follows: D/Dt (m/s) [29, 30]. Three measurements were performed to obtain average results for every participant. Only measurements that met the automatic quality control cutoff were used in the final analysis. All the PWV measurements were performed by two experienced nurses.

Statistical analysis

The data were analyzed by using IBM® SPSS® Statistics Version 23 software (copyright 2015). Variables were tested for normality by way of Shapiro–Wilk or Kolmogorov–Smirnov tests, as appropriate. Data are presented as mean ± standard deviation (SD) if not mentioned otherwise. Differences in continuous variables between GDM participants and controls were studied by using Student’s t test in cases of normality and the Mann–Whitney U test in cases of skewed distribution of measurements.

All 240 women were also included in subgroup analyses to study the effect of obesity on the results. For these analyses, we divided the whole study group into four subgroups according to obesity and previous GDM. Obesity was classified as BMI ≥30 kg/m2 [47]. The clinical characteristics of these four subgroups were studied by way of one-way ANOVA in cases of normality and by using the Kruskal–Wallis test in cases of non-normality. If the overall p value was significant, individual p values between subgroups were also calculated. Post hoc analyses, with a conservative Bonferroni correction factor, were performed in order to correct for multiple testing. The relationships between different cardiovascular risk factors were tested by Pearson’s or Spearman’s correlation analysis, as appropriate.

Further, we conducted univariate linear regression analyses for hsCRP, MMP-8, TIMP-1, PWV and arterial compliance index values to find possible associations with clinically relevant covariates. Then multivariable linear analyses were carried out to examine whether simple associations were changed after adjustment for potential confounders. Finally, stepwise multiple linear regression analyses were done to find out relevant covariates to final models. The selected covariates in all of these analyses were age, BMI, previous GDM, time after the index pregnancy, pack-years of smoking, heart rate, systolic blood pressure, hsCRP, TC and fasting insulin. F-statistics was used to optimize the sequential variable selection procedure. A two-tailed probability value of <0.05 was considered significant.


The basic clinical characteristics of the study participants are summarized in Table 1. There were no significant differences between the two cohorts in self-reported history of respiratory infection, other infectious disease or trauma during the month before follow-up laboratory examinations.
Table 1

Basic clinical characteristics of women with GDM and controls




p value

Average time since delivery, years

3.7 ± 1.0

3.7 ± 0.9


Age, years

35.8 ± 4.4

35.9 ± 4.6


Primiparous, n (%)

23 (19.2%)

23 (19.2%)


Therapy of GDM during pregnancy

 Insulin, n (%)

24 (20.0%)


 Metformin, n (%)

1 (0.8%)


 Dietary therapy, n (%)

95 (79.2%)


Pack-years of smoking

3.8 ± 6.0

2.4 ± 4.6


During the previous month, history of

 Respiratory infection, n (%)

45 (37.5%)

44 (36.7%)


 Other infectious disease, n (%)

18 (15.0%)

10 (8.3%)


 Trauma, n (%)

9 (7.5%)

5 (4.2%)


BMI, kg/m2

28.3 ± 5.0

27.5 ± 5.4


Systolic BP, mmHg

122.4 ± 12.5

119.0 ± 11.5


Diastolic BP, mmHg

73.5 ± 9.0

71.8 ± 8.7


Heart rate, beats per minute

65.9 ± 9.1

63.8 ± 9.6


TC, mmol/L

4.7 ± 0.9

4.6 ± 0.8


F-Gluc, mmol/L

5.6 ± 0.6

5.3 ± 0.3


F-Insu, mU/L

5.2 ± 3.6

4.6 ± 3.6


Data are presented as mean ± SD if not mentioned otherwise

BMI body mass index, BP blood pressure, F-Gluc fasting glucose, F-Insu fasting insulin, TC total cholesterol

Subclinical inflammation

Serum TIMP-1 levels were significantly increased after previous GDM (Table 2). There was a significant positive association between previous GDM and TIMP-1 levels in both univariate and multivariable linear regression analyses (data not shown). There were no differences in the concentrations of MMP-8 and MMP-9 between the groups (Table 2). In stepwise multiple linear regression analyses, hsCRP, previous GDM and TC were important determinants of MMP-8 levels. Likewise, previous GDM, together with BMI and heart rate associated with TIMP-1 in stepwise multiple linear regression analyses. Nevertheless, the significant determinants explained only 13.8% of MMP-8 and 6.7% of TIMP-1 concentrations (Table 3).
Table 2

Results of primary analyses of GDM and control groups




p value

hsCRP, mg/L

2.50 ± 3.69

2.50 ± 4.19


MMP-8, ng/mL

27.83 ± 1.48

32.78 ± 1.90


MMP-9, ng/mL

384.27 ± 13.15

392.15 ± 12.60


TIMP-1, ng/mL

102.80 ± 29.72

94.58 ± 24.51


MMP-8/TIMP-1, mol ratio

0.13 ± 0.009

0.17 ± 0.015


MMP-9/TIMP-1, mol ratio

1.32 ± 0.078

1.43 ± 0.085


C1, mL/mmHg × 10

15.14 ± 3.51

15.85 ± 3.36


C2, mL/mmHg × 100

8.44 ± 3.08

8.60 ± 3.20


PWV, m/s

6.44 ± 0.83

6.17 ± 0.74


Data are presented as mean ± SD

hsCRP high-sensitivity C reactive protein, C1 large artery compliance index, C2 small artery compliance index, PWV pulse wave velocity, MMP-8 matrix metalloproteinase-8, MMP-9 matrix metalloproteinase-9

Table 3

Results of stepwise multiple linear regression analyses


Covariates included in the model

R2 for model

Global p

Standardized β

p value



















Previous GDM













Previous GDM








Heart rate









Systolic BP




Heart rate








Time after the index pregnancy









Systolic BP




Heart rate












Pack-years of smoking









Systolic BP












Previous GDM




Time after the index pregnancy




Covariates in these analyses included age, BMI, previous GDM, pack-years of smoking, time after the index pregnancy, heart rate, systolic blood pressure, hsCRP, TC and fasting insulin. Final models include significant covariates only. Standardized β provides a measure of the relative strength of an association, independent of the measurement units. Standardized β and p values are shown only when p < 0.05

BMI body mass index, BP blood pressure, F-Insu fasting insulin, GDM gestational diabetes mellitus, hsCRP high-sensitivity C reactive protein, C1 large artery compliance index, C2 small artery compliance index, PWV pulse wave velocity, MMP-8 matrix metalloproteinase-8

We found no difference in the concentrations of hsCRP between GDM cases and controls (Table 2), even when participants affected with infections or traumas were excluded (data not shown). In stepwise multiple linear regression analysis (Table 3), only BMI was a significant determinant of hsCRP levels, but the model explained only 9.6% of hsCRP values. Previous GDM did not influence hsCRP concentrations in our data.

Pulse wave velocity and arterial compliance

PWV values differed significantly between the GDM cases and controls (Table 2). In univariate linear regression analysis, there were significant associations with age (p < 0.001), fasting insulin (p < 0.001), previous GDM (p = 0.009), TC (p < 0.001), heart rate (p < 0.001), systolic blood pressure (p < 0.001) and BMI (p < 0.001). In stepwise multiple linear regression analysis, significant determinants of PWV values were systolic BP, age, insulin levels, previous GDM and time after the index pregnancy. Covariates explained 47.0% of PWV (Table 3). In our two study cohorts, there were no interactions between previous GDM and TIMP1 on PWV (data not shown).

There was a nonsignificant difference in C1 values between the study groups. No difference was revealed in C2 values, either. In univariate linear regression analysis, there was no significant association between C2 and BMI (p = 0.726), but an inverse association between C1 and BMI was significant (p = 0.025). In stepwise multiple linear regression analysis, systolic BP, heart rate, BMI and time after the index pregnancy were significant covariates explaining 52.4% of C1 values. Significant determinants of C2 values were systolic BP, heart rate, BMI, age and pack-years of smoking. These covariates explained 31.7% of C2 values (Table 3).

Effect of obesity in subgroups

Altogether, there were 75 women in the obese group (BMI ≥ 30 kg/m2); 43 GDM and 32 control participants. The non-obese group (BMI < 30 kg/m2; n = 165) consisted of 77 GDM and 88 control participants [55]. In subgroup analyses, participants in obese subgroups had higher serum concentrations of hsCRP than those in non-obese subgroups, as shown in Fig. 1. The concentrations of MMP-8 in the four subgroups were as follows: obese GDM cases, 27.76 ± 1.77 ng/mL, obese controls 37.10 ± 4.16 ng/mL, non-obese GDM cases, 27.88 ± 2.08 ng/mL and non-obese controls, 31.21 ± 2.10 ng/mL. The concentration of MMP-8 was highest among obese controls, but the differences between the four subgroups were not significant (p = 0.090). We also found no differences in the levels of MMP-9 or TIMP-1 between these four subgroups (data not shown). Between the subgroups, there were no differences in the MMP-8/TIMP-1 or MMP-9/TIMP-1 ratio either (data not shown). In the four subgroups, differences in PWV values were significant, but differences in both C1 and C2 values were not (Figs. 2, 3).
Fig. 1
Fig. 1

Serum concentrations of hsCRP in the four subgroups. Median values (minimum, maximum) of hsCRP: among obese GDM women 2.1 (0.0, 12.4) mg/mL, obese control women 2.1 (0.3, 18.5) mg/mL, non-obese GDM women 0.9 (0.0, 32.3) mg/mL, and non-obese control women 0.7 (0.0, 25.7) mg/mL. Values of more than 10 mg/mL were measured by turbidimetric immunoassay. The overall p value is given at the bottom. Individual p values for pairwise comparisons are also presented

Fig. 2
Fig. 2

PWV in the four subgroups. Median values (minimum, maximum) of PWV: among obese GDM women 6.8 (5.6, 9.7) m/s, obese control women 6.6 (4.8, 8.5) m/s, non-obese GDM women 6.3 (4.9, 9.2) m/s, and non-obese control women 6.0 (4.5, 7.9) m/s. The overall p value is given at the bottom. Individual p values for pairwise comparisons are also presented

Fig. 3
Fig. 3

Large (a) and small (b) artery compliance index values in the four subgroups. a Median values (minimum, maximum) of the large-artery compliance index (C1): among obese GDM women 13.3 (9.1, 21.8) mL/mmHg × 10, obese control women 14.7 (10.2, 23.5) mL/mmHg × 10, non-obese GDM women 15.2 (7.2, 25.2) mL/mmHg × 10, and non-obese control women 15.9 (7.5, 25.7) mL/mmHg × 10. The overall p value is given. b Median values (minimum, maximum) of the small-artery compliance index (C2): among obese GDM women 8.8 (2.8, 15.2) mL/mmHg × 100, obese control women 8.6 (2.2, 17.7) mL/mmHg × 100, non-obese GDM women 8.1 (1.8, 17.6) mL/mmHg × 100, and non-obese control women 8.1 (2.4, 16.0) mL/mmHg × 100. The overall p value is given


Our main finding was that PWV was significantly higher after GDM than after normoglycemic pregnancy. This was supported by a nonsignificant difference in the large-artery compliance index, C1, which indicates that the arteries of GDM cases were less distensible than those of the controls. Secondly, subclinical low-grade inflammation and reduced arterial compliance especially affected women with high BMI.

Inflammation has been shown to be a strong predictor of women’s cardiovascular complications [48]. We found that levels of TIMP-1 were significantly upregulated after previous GDM, reflecting low-grade inflammation among this relatively healthy and young study population. No differences were found in circulating levels of MMP-8 or MMP-9 between the two study cohorts. In subgroup analyses, the highest levels of MMP-8 were in obese controls, but this did not reach statistical significance either. A search of MEDLINE (English language; 1989–September 2016; search terms: “MMP-8, MMP-9, TIMP-1” and “GDM”) revealed no publications concerning female populations where levels of MMP-8, MMP-9 or TIMP-1 have been studied in connection with previous GDM.

There is evidence that glucose can modulate the expression, production and activity of MMPs. For example, endothelial cells cultured in hyperglycemic conditions present increased expression and activity of MMP-9 [49]. It is a pity that there were no samples left for MMP analysis taken from the patients during the period when they suffered from gestational diabetes. We might postulate, that during the pregnancy GDM increase concentrations of MMPs and they in turn upregulate TIMP-1. After the delivery, the decreasing concentrations of glucose, MMPs and TIMP-1 take place consecutively. The prolonged upregulation of TIMP-1 found in this study without upregulated MMP levels may also be a result of the fact that upregulated TIMP-1 may suppress MMP-8 and MMP-9 levels. Further, third explanation for prolonged TIMP-1 upregulation found in this work may be that prolonged elevation of TIMP-1 levels may mediate MMP-independent pro-inflammatory or growth-factor-like signaling functions contributing to low-grade inflammation [3638].

Recent studies have reported higher CRP and hsCRP levels in women with a history of GDM than in age-matched normal controls after a 1- or 5-year postpartum period [4, 50, 51]. On the contrary, Ajala et al. found no difference in CRP in women after previous GDM compared to controls 4–10 years postpartum [52]. In our study, when hsCRP was determined on average at 3.7 years after delivery, there was no difference between the age-matched study cohorts. However, low-grade inflammation was evident among obese women, in contrast to non-obese participants in subgroup analyses. The GDM and non-GDM women of our study did not differ in BMI, which can partly explain the similar hsCRP levels between the two study cohorts.

Only a few studies have been published concerning a possible relationship between PWV and previous GDM. Lekva et al. reported an enhanced cardiovascular risk at 5-year follow-up as reflected in elevated PWV after previous GDM diagnosed using the old criteria of the World Health Organization (WHO) (OGTT: 2-h plasma glucose ≥7.8 mmol/L). However, they did not find such an association in PWV when using IADPSG diagnostic criteria (OGTT: fasting plasma glucose 5.1–6.9 mmol/L, 1-h plasma glucose ≥10.0 mmol/L or 2-h plasma glucose 8.5–11.0 mmol/L) [41, 53]. Using diagnostic criteria of GDM similar to those of the IADPSG [39], we observed a significant increase in PWV in women with previous GDM. Previous GDM was also a significant determinant of PWV in multiple linear regression analysis. Our results are in accordance with those of Tam et al., who reported higher PWV in women with a history of GDM followed up at a median of 6 years postpartum [54]. In contrast to these findings, Heitritter et al. detected no difference in PWV at an average of 1 year after previous GDM compared with normoglycemic pregnancy [4]. There were no significant differences in C1 or C2 values between the GDM cases and controls. In a recent study, no difference was found in vascular function measured also by using HDI/PulseWave™CR-2000 in women with a history of GDM when compared to healthy controls 4–10 years postpartum, either [52].

Strengths of our study include the fact that we used standardized measurements of arterial stiffness. Determination of systemic arterial stiffness by using HDI/PulseWave™CR-2000 equipment is widely used, and carotid-femoral PWV is accepted as the most reliable measurement of arterial stiffness [29]. We measured the levels of MMP-8, MMP-9 and TIMP-1 by specific immunoassays previously found to be suitable for diagnosis and monitoring of systemic low-grade inflammation associated with cardiovascular and infectious diseases as well as other inflammatory states [11, 1318, 2325]. Further, we performed a well characterized hospital-based study of two cohorts of women with a similar follow-up time and age. Moreover, there was no significant difference in BMI between the study groups, and all participants had undergone OGTT screening during the index pregnancy. Since low-risk parturients do not routinely undergo OGTTs in Finland [39], this last strength may also turn out to be a weakness, because the most low-risk women had to be excluded from our study [40]. Although the relatively short time from delivery to the follow-up study allowed us to observe early cardiovascular changes, it may be one of our study limitations as well, since major differences between the study groups are probably better observable later in their life. For example, within 7 years postpartum, previous GDM was identified as a risk factor of CVD by Goueslard et al. They studied database of more than 1.5 million deliveries and found that the incidence of myocardial infarction was 0.04% in women with a history of GDM and 0.02% without [7].

In our subgroup analyses, obesity was associated with higher levels of hsCRP and higher values of PWV. We have earlier revealed the effect of obesity being similar with many other markers for cardio-metabolic risks among the four subgroups [40, 55]. Earlier, BMI has been shown to associate inversely with arterial compliance [56]. As presented in Fig. 3, this seemed to be the case also in our study in C1 values. Surprisingly, in multiple regression analyses, BMI seemed to be protective as regards arterial compliance (C1 and C2). BMI was significantly correlated with systolic blood pressure and heart rate (data not shown). Hence, adjusted findings concerning C1 and C2 might have been affected by these relationships irrespective of possible biologic associations. In our opinion, this result may be explained by multiple interactions of C1 and C2 measurements with other confounding variables. This was supported by the findings of univariate analysis and stepwise multiple linear regression analysis without systolic BP and heart rate as covariates, where inverse association between BMI and C1 was found and association between BMI and C2 was vanished (data not shown).

The prevalence of obesity is increasing around the world [57]. Specifically, visceral obesity modifies glucose and lipid metabolism. It is associated with increased risk of arterial stiffness and atherosclerosis both in normal-weight subjects and patients with T2DM [58, 59]. Our results imply that in preventing cardiovascular risk among women after delivery, we need a comprehensive attitude in clinical care instead of concentrating on single factors.


When studied 3.7 years after delivery, PWV values were higher in women with previous GDM, indicating that their arteries are less distensible than those in women with previous normoglycemic pregnancy. Among other findings, this relationship was even more evident in obese subjects. We also found that serum levels of TIMP-1 were significantly upregulated after previous GDM, reflecting low-grade inflammation among this relatively healthy and young study population. Altogether, our results demonstrate that previous GDM may reflect a subclinical inflammatory state and together with obesity may contribute to an early stage of the subclinical atherosclerotic process even in relatively young and healthy women.



large artery compliance index


small artery compliance index


gestational diabetes mellitus


high-sensitivity C-reactive protein


matrix metalloproteinase


oral glucose tolerance test


pulse wave velocity


tissue inhibitor of metalloproteinase


Authors’ contributions

TV-K participated in the design of the study, conducted experiments, performed data analyses and drafted the manuscript. AL participated in the design of the study and contributed to drafting the manuscript. TT carried out the analyses of MMP-8, MMP-9 and TIMP-1. OP, JU and TS contributed to drafting the manuscript. AP designed the study, helped to perform data analyses and drafted the manuscript. All authors read and approved the final manuscript.


We appreciate the professional technical aid of Anna Silén, Taru Stranden, Hanna Kujanen, Ari Virta, Nick Bolton, Kirsti Räsänen and Piia Suursalmi. We sincerely acknowledge the work of the clinical staff of Linnan Klinikka and Kanta-Häme Central Hospital.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and analyzed during the current study are not publicly available as a result of the fact that individual privacy could be compromised, but are available from the corresponding author on reasonable request.

Consent for publication

Every participant was given both oral and written information on the study before she signed an informed consent document.

Ethics approval and consent to participate

The study was conducted in accordance with the ethical principles outlined in the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of Kanta-Häme Hospital District (Reference Number 521/2010; date of approval 21.12.2010).


This study was supported by grants from the Finnish Cultural Foundation, Häme Regional Fund and the Ministry of Health and Social Welfare in Finland via Medical Research Funds of Kanta-Häme Central Hospital, Tampere University Hospital, Helsinki University Hospital EVO and Karolinska Institutet.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

School of Medicine, University of Tampere, Tampere, Finland
Department of Obstetrics and Gynecology, Tampere University Hospital, Box 2000, 33521 Tampere, Finland
Department of Infectious Diseases, Inflammation Center, Helsinki University Hospital, Helsinki, Finland
Clinicum, University of Helsinki, Helsinki, Finland
The Social Insurance Institution of Finland, Benefit Services, Helsinki, Finland
Department of Oral and Maxillofacial Diseases, Helsinki University and University Hospital, Helsinki, Finland
Division of Periodontology, Department of Dental Medicine, Karolinska Institutet, Huddinge, Sweden
Department of Emergency Medicine, Kanta-Häme Central Hospital, Hämeenlinna, Finland


  1. Dabelea D, Snell-Bergeon JK, Hartsfield CL, Bischoff KJ, Hamman RF, McDuffie RS. Kaiser Permanente of Colorado GDM Screening Program: increasing prevalence of gestational diabetes mellitus (GDM) over time and by birth cohort: Kaiser Permanente of Colorado GDM Screening Program. Diabetes Care. 2005;28(3):579–84.View ArticlePubMedGoogle Scholar
  2. Vuori E, Gissler M. Perinatal statistics: parturients, deliveries and newborns 2015. Statistical report 16/2016. Helsinki: National Institute for Health and Welfare; 2016.Google Scholar
  3. Gobl CS, Bozkurt L, Yarragudi R, Prikoszovich T, Tura A, Pacini G, Koppensteiner R, Kautzky-Willer A. Biomarkers of endothelial dysfunction in relation to impaired carbohydrate metabolism following pregnancy with gestational diabetes mellitus. Cardiovasc Diabetol. 2014;13(1):138.View ArticlePubMedPubMed CentralGoogle Scholar
  4. Heitritter SM, Solomon CG, Mitchell GF, Skali-Ounis N, Seely EW. Subclinical inflammation and vascular dysfunction in women with previous gestational diabetes mellitus. J Clin Endocrinol Metab. 2005;90(7):3983–8.View ArticlePubMedGoogle Scholar
  5. Vrachnis N, Augoulea A, Iliodromiti Z, Lambrinoudaki I, Sifakis S, Creatsas G. Previous gestational diabetes mellitus and markers of cardiovascular risk. Int J Endocrinol. 2012;2012:458610.PubMedPubMed CentralGoogle Scholar
  6. Shah BR, Retnakaran R, Booth GL. Increased risk of cardiovascular disease in young women following gestational diabetes mellitus. Diabetes Care. 2008;31(8):1668–9.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Goueslard K, Cottenet J, Mariet AS, Giroud M, Cottin Y, Petit JM, Quantin C. Early cardiovascular events in women with a history of gestational diabetes mellitus. Cardiovasc Diabetol. 2016;15:15.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Retnakaran R, Shah BR. Role of type 2 diabetes in determining retinal, renal, and cardiovascular outcomes in women with previous gestational diabetes mellitus. Diabetes Care. 2017;40(1):101–8.View ArticlePubMedGoogle Scholar
  9. Karadeniz M, Duran M, Akyel A, Yarlioglues M, Ocek AH, Celik IE, Kilic A, Yalcin AA, Ergun G, Murat SN. High sensitive crp level is associated with intermediate and high syntax score in patients with acute coronary syndrome. Int Heart J. 2015;56(4):377–80.View ArticlePubMedGoogle Scholar
  10. Lenglet S, Mach F, Montecucco F. Role of matrix metalloproteinase-8 in atherosclerosis. Mediat Inflamm. 2013;2013:659282.View ArticleGoogle Scholar
  11. Sorsa T, Tjaderhane L, Konttinen YT, Lauhio A, Salo T, Lee HM, Golub LM, Brown DL, Mäntylä P. Matrix metalloproteinases: contribution to pathogenesis, diagnosis and treatment of periodontal inflammation. Ann Med. 2006;38(5):306–21.View ArticlePubMedGoogle Scholar
  12. Lauhio A, Salo T, Ding Y, Konttinen YT, Nordström D, Tschesche H, Lähdevirta J, Golub LM, Sorsa T. In vivo inhibition of human neutrophil collagenase (MMP-8) activity during long-term combination therapy of doxycycline and non-steroidal anti-inflammatory drugs (NSAID) in acute reactive arthritis. Clin Exp Immunol. 1994;98(1):21–8.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Lauhio A, Konttinen YT, Tschesche H, Nordström D, Salo T, Lähdevirta J, Golub LM, Sorsa T. Reduction of matrix metalloproteinase 8-neutrophil collagenase levels during long-term doxycycline treatment of reactive arthritis. Antimicrob Agents Chemother. 1994;38(2):400–2.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Lauhio A, Salo T, Tjaderhane L, Lähdevirta J, Golub LM, Sorsa T. Tetracyclines in treatment of rheumatoid arthritis. Lancet. 1995;346(8975):645–6.View ArticlePubMedGoogle Scholar
  15. Lauhio A, Saikku P, Salo T, Tschesche H, Lähdevirta J, Sorsa T. Combination treatment in Chlamydia-triggered reactive arthritis: comment on the article by Carter et al. Arthritis Rheum. 2011;63(1):305–7 (author reply 307–8).View ArticlePubMedGoogle Scholar
  16. Lauhio A, Hästbacka J, Pettilä V, Tervahartiala T, Karlsson S, Varpula T, Varpula M, Ruokonen E, Sorsa T, Kolho E. Serum MMP-8, -9 and TIMP-1 in sepsis: high serum levels of MMP-8 and TIMP-1 are associated with fatal outcome in a multicentre, prospective cohort study. Hypothetical impact of tetracyclines. Pharmacol Res. 2011;64(6):590–4.View ArticlePubMedGoogle Scholar
  17. Lauhio A, Färkkilä E, Pietiläinen KH, Åström P, Winkelmann A, Tervahartiala T, Pirilä E, Rissanen A, Kaprio J, Sorsa TA, Salo T. Association of MMP-8 with obesity, smoking and insulin resistance. Eur J Clin Investig. 2016;46(9):757–65.View ArticleGoogle Scholar
  18. Rautelin HI, Oksanen AM, Veijola LI, Sipponen PI, Tervahartiala TI, Sorsa TA, Lauhio A. Enhanced systemic matrix metalloproteinase response in Helicobacter pylori gastritis. Ann Med. 2009;41(3):208–15.View ArticlePubMedGoogle Scholar
  19. Cena JJ, Lalu MM, Cho WJ, Chow AK, Bagdan ML, Daniel EE, Castro MM, Schulz R. Inhibition of matrix metalloproteinase activity in vivo protects against vascular hyporeactivity in endotoxemia. Am J Physiol Heart Circ Physiol. 2010;298(1):H45–51.View ArticlePubMedGoogle Scholar
  20. Siasos G, Tousoulis D, Kioufis S, Oikonomou E, Siasou Z, Limperi M, Papavassiliou AG, Stefanadis C. Inflammatory mechanisms in atherosclerosis: the impact of matrix metalloproteinases. Curr Top Med Chem. 2012;12(10):1132–48.View ArticlePubMedGoogle Scholar
  21. Paim LR, Schreiber R, Matos-Souza JR, Silva AA, Campos LF, Azevedo ER, Alonso K, de Rossi G, Etchebehere M, Gorla JI, Cliquet A Jr, Nadruz W Jr. Oxidized low-density lipoprotein, matrix-metalloproteinase-8 and carotid atherosclerosis in spinal cord injured subjects. Atherosclerosis. 2013;231(2):341–5.View ArticlePubMedGoogle Scholar
  22. Goncalves FM, Jacob-Ferreira AL, Gomes VA, Casella-Filho A, Chagas AC, Marcaccini AM, Gerlach RF, Tanus-Santos JE. Increased circulating levels of matrix metalloproteinase (MMP)-8, MMP-9, and pro-inflammatory markers in patients with metabolic syndrome. Clin Chim Acta. 2009;403(1–2):173–7.View ArticlePubMedGoogle Scholar
  23. Pussinen PJ, Sarna S, Puolakkainen M, Ohlin H, Sorsa T, Pesonen E. The balance of serum matrix metalloproteinase-8 and its tissue inhibitor in acute coronary syndrome and its recurrence. Int J Cardiol. 2013;167(2):362–8.View ArticlePubMedGoogle Scholar
  24. Sorsa T, Tervahartiala T, Leppilahti J, Hernandez M, Gamonal J, Tuomainen AM, Lauhio A, Pussinen PJ, Mäntylä P. Collagenase-2 (MMP-8) as a point-of-care biomarker in periodontitis and cardiovascular diseases. Therapeutic response to non-antimicrobial properties of tetracyclines. Pharmacol Res. 2011;63(2):108–13.View ArticlePubMedGoogle Scholar
  25. Tuomainen AM, Nyyssönen K, Laukkanen JA, Tervahartiala T, Tuomainen TP, Salonen JT, Sorsa T, Pussinen PJ. Serum matrix metalloproteinase-8 concentrations are associated with cardiovascular outcome in men. Arterioscler Thromb Vasc Biol. 2007;27(12):2722–8.View ArticlePubMedGoogle Scholar
  26. Meigs JB, Hu FB, Rifai N, Manson JE. Biomarkers of endothelial dysfunction and risk of type 2 diabetes mellitus. JAMA. 2004;291(16):1978–86.View ArticlePubMedGoogle Scholar
  27. Boutouyrie P, Tropeano AI, Asmar R, Gautier I, Benetos A, Lacolley P, Laurent S. Aortic stiffness is an independent predictor of primary coronary events in hypertensive patients: a longitudinal study. Hypertension. 2002;39(1):10–5.View ArticlePubMedGoogle Scholar
  28. Nigam A, Mitchell GF, Lambert J, Tardif J. Relation between conduit vessel stiffness (assessed by tonometry) and endothelial function (assessed by flow-mediated dilatation) in patients with and without coronary heart disease. Am J Cardiol. 2003;92(4):395–9.View ArticlePubMedGoogle Scholar
  29. Laurent S, Cockcroft J, Van Bortel L, Boutouyrie P, Giannattasio C, Hayoz D, Pannier B, Vlachopoulos C, Wilkinson I, Struijker-Boudier H. European network for non-invasive investigation of large arteries: expert consensus document on arterial stiffness: methodological issues and clinical applications. Eur Heart J. 2006;27(21):2588–605.View ArticlePubMedGoogle Scholar
  30. Agabiti-Rosei E, Mancia G, O’Rourke MF, Roman MJ, Safar ME, Smulyan H, Wang J, Wilkinson IB, Williams B, Vlachopoulos C. Central blood pressure measurements and antihypertensive therapy. A consensus document. Hypertension. 2007;50:154–60.View ArticlePubMedGoogle Scholar
  31. Nelson MR, Stepanek J, Cevette M, Covalciuc M, Hurst RT, Tajik J. Noninvasive measurement of central vascular pressures with arterial tonometry: clinical revival of the pulse pressure waveform? Mayo Clin Proc. 2010;85(5):460–72.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Kim YK. Impact of the metabolic syndrome and its components on pulse wave velocity. Korean J Intern Med. 2006;21(2):109–15.View ArticlePubMedPubMed CentralGoogle Scholar
  33. Cruickshank K, Riste L, Anderson SG, Wright JS, Dunn G, Gosling RG. Aortic pulse-wave velocity and its relationship to mortality in diabetes and glucose intolerance: an integrated index of vascular function? Circulation. 2002;106(16):2085–90.View ArticlePubMedGoogle Scholar
  34. Kormi I, Alfakry H, Tervahartiala T, Pussinen PJ, Sinisalo J, Sorsa T. The effect of prolonged systemic doxycycline therapy on serum tissue degrading proteinases in coronary bypass patients: a randomized, double-masked, placebo-controlled clinical trial. Inflamm Res. 2014;63(5):329–34.View ArticlePubMedGoogle Scholar
  35. Payne JB, Golub LM, Stoner JA, Lee HM, Reinhardt RA, Sorsa T, Slepian MJ. The effect of subantimicrobial-dose-doxycycline periodontal therapy on serum biomarkers of systemic inflammation: a randomized, double-masked, placebo-controlled clinical trial. J Am Dent Assoc. 2011;142(3):262–73.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Hayakawa T, Yamashita K, Tanzawa K, Uchijima E, Iwata K. Growth-promoting activity of tissue inhibitor of metalloproteinases-1 (TIMP-1) for a wide range of cells. A possible new growth factor in serum. FEBS Lett. 1992;298(1):29–32.View ArticlePubMedGoogle Scholar
  37. Moore CS, Crocker SJ. An alternate perspective on the roles of TIMPs and MMPs in pathology. Am J Pathol. 2012;180(1):12–6.View ArticlePubMedGoogle Scholar
  38. Stetler-Stevenson WG. Tissue inhibitors of metalloproteinases in cell signaling: metalloproteinase-independent biological activities. Sci Signal. 2008;1(27):re6.View ArticlePubMedPubMed CentralGoogle Scholar
  39. Kaaja R, Alenius H, Kinnunen T, Komulainen J, Peränen N, Rönnemaa T, Saramies J, Soukka H, Teramo K, Vuorela P, Vääräsmäki M. Gestational diabetes (online). Current care guidelines. Working group set up by the Finnish Medical Society Duodecim, the Medical Advisory Board of the Finnish Diabetes Association and the Finnish Gynecological Association; 2013. Accessed 25 June 2013.
  40. Vilmi-Kerälä T, Palomäki O, Vainio M, Uotila J, Palomäki A. The risk of metabolic syndrome after gestational diabetes mellitus—a hospital-based cohort study. Diabetol Metab Syndr. 2015;7:43.View ArticlePubMedPubMed CentralGoogle Scholar
  41. Report of a World Health Organization Consultation. Diagnostic criteria and classification of hyperglycaemia first detected in pregnancy: a World Health Organization guideline. Diabetes Res Clin Pract. 2014;103(3):341–63.View ArticleGoogle Scholar
  42. Saquib N, Stefanick ML, Natarajan L, Pierce JP. Mortality risk in former smokers with breast cancer: pack-years vs. smoking status. Int J Cancer. 2013;133(10):2493–7.View ArticlePubMedPubMed CentralGoogle Scholar
  43. World Medical Association Inc. Declaration of Helsinki. Ethical principles for medical research involving human subjects. J Indian Med Assoc. 2009;107(6):403–5.Google Scholar
  44. Chenillot O, Henny J, Steinmetz J, Herbeth B, Wagner C, Siest G. High sensitivity C-reactive protein: biological variations and reference limits. Clin Chem Lab Med. 2000;38(10):1003–11.View ArticlePubMedGoogle Scholar
  45. Sanchez A, Mirabel JL, Barrenechea E, Eugui J, Puelles A, Castaneda A. Evaluation of an improved immunoturbidimetic assay for serum C-reactive protein on a COBAS INTEGRA 400 Analyzer. Clin Lab. 2002;48(5–6):313–7.PubMedGoogle Scholar
  46. Cohn JN, Finkelstein S, McVeigh G, Morgan D, LeMay L, Robinson J, Mock J. Noninvasive pulse wave analysis for the early detection of vascular disease. Hypertension. 1995;26(3):503–8.View ArticlePubMedGoogle Scholar
  47. Obesity: preventing and managing the global epidemic. Report of a WHO consultation. World Health Organ Tech Rep Ser, vol 894. 2000; p. i–xii, 1–253.Google Scholar
  48. Ridker PM, Rifai N, Rose L, Buring JE, Cook NR. Comparison of C-reactive protein and low-density lipoprotein cholesterol levels in the prediction of first cardiovascular events. N Engl J Med. 2002;347(20):1557–65.View ArticlePubMedGoogle Scholar
  49. Berg G, Miksztowicz V, Schreier L. Metalloproteinases in metabolic syndrome. Clin Chim Acta. 2011;412(19–20):1731–9.View ArticlePubMedGoogle Scholar
  50. Ozuguz U, Isik S, Berker D, Arduc A, Tutuncu Y, Akbaba G, Gokay F, Guler S. Gestational diabetes and subclinical inflammation: evaluation of first year postpartum outcomes. Diabetes Res Clin Pract. 2011;94(3):426–33.View ArticlePubMedGoogle Scholar
  51. Lekva T, Michelsen AE, Bollerslev J, Norwitz ER, Aukrust P, Henriksen T, Ueland T. Low circulating pentraxin 3 levels in pregnancy is associated with gestational diabetes and increased apoB/apoA ratio: a 5-year follow-up study. Cardiovasc Diabetol. 2016;15:23.View ArticlePubMedPubMed CentralGoogle Scholar
  52. Ajala O, Jensen LA, Ryan E, Chik C. Women with a history of gestational diabetes on long-term follow up have normal vascular function despite more dysglycemia, dyslipidemia and adiposity. Diabetes Res Clin Pract. 2015;110(3):309–14.View ArticlePubMedGoogle Scholar
  53. Lekva T, Bollerslev J, Norwitz ER, Aukrust P, Henriksen T. Aortic stiffness and cardiovascular risk in women with previous gestational diabetes mellitus. PLoS ONE. 2015;10(8):e0136892.View ArticlePubMedPubMed CentralGoogle Scholar
  54. Tam WH, Ma RC, Chan JC, Lao TT, Chan MH, Li CY. PP103. Arterial stiffness in women with previous GDM—a follow up of Chinese HAPO study cohort. Pregnancy Hypertens. 2012;2(3):295.PubMedGoogle Scholar
  55. Vilmi-Kerälä T, Palomäki O, Kankkunen P, Juurinen L, Uotila J, Palomäki A. Oxidized LDL, insulin resistance and central blood pressure after gestational diabetes mellitus. Acta Obstet Gynecol Scand. 2016;95(12):1425–32.View ArticlePubMedGoogle Scholar
  56. Acree LS, Montgomery PS, Gardner AW. The influence of obesity on arterial compliance in adult men and women. Vasc Med. 2007;12(3):183–8.View ArticlePubMedGoogle Scholar
  57. Huffman MD, Capewell S, Ning H, Shay CM, Ford ES, Lloyd-Jones DM. Cardiovascular health behavior and health factor changes (1988–2008) and projections to 2020: results from the National Health and Nutrition Examination Surveys. Circulation. 2012;125(21):2595–602.View ArticlePubMedPubMed CentralGoogle Scholar
  58. Kim S, Kung C, Park JS, Lee SP, Kim HK, Ahn CW, Kim KR, Kang S. Normal-weight obesity is associated with increased risk of subclinical atherosclerosis. Cardiovasc Diabetol. 2015;14:58.View ArticlePubMedPubMed CentralGoogle Scholar
  59. Bouchi R, Ohara N, Asakawa M, Nakano Y, Takeuchi T, Murakami M, Sasahara Y, Numasawa M, Minami I, Izumiyama H, Hashimoto K, Yoshimoto T, Ogawa Y. Is visceral adiposity a modifier for the impact of blood pressure on arterial stiffness and albuminuria in patients with type 2 diabetes? Cardiovasc Diabetol. 2016;15:10.View ArticlePubMedPubMed CentralGoogle Scholar


© The Author(s) 2017