Article Text

Early pregnancy plasma fatty acid profiles of women later diagnosed with gestational diabetes
  1. Ellen Alma Tryggvadottir1,2,
  2. Ingibjorg Gunnarsdottir1,2,
  3. Bryndis Eva Birgisdottir1,2,
  4. Laufey Hrolfsdottir1,3,
  5. Rikard Landberg4,
  6. Ingibjorg Th Hreidarsdottir5,6,
  7. Hildur Hardardottir6,
  8. Thorhallur Ingi Halldorsson1,2
  1. 1Faculty of Food Science and Nutrition, University of Iceland, Reykjavik, Iceland
  2. 2Unit for Nutrition Research, Landspitali University Hospital, Reykjavík, Iceland
  3. 3Department of Education and Science, Akureyri Hospital, Akureyri, Iceland
  4. 4Biology and Biological Engineering, Chalmers University of Technology, Goteborg, Sweden
  5. 5Department of Obstetrics and Gynecology, Landspítali University Hospital, Reykjavík, Iceland
  6. 6Faculty of Medicine, University of Iceland, Reykjavik, Iceland
  1. Correspondence to Ellen Alma Tryggvadottir; eat2{at}hi.is

Abstract

Introduction Fatty acid (FA) concentrations have previously been associated with gestational diabetes mellitus (GDM). However, few studies on GDM have examined FA profiles in early pregnancy or before diagnosis. This study aimed to compare early pregnancy plasma FA profiles of women with and without GDM diagnoses as well as their reported dietary consumption.

Research design and methods The subjects comprised 853 women from the prospective study: Pregnant Women in Iceland II (PREWICE II), attending their 11–14 weeks ultrasound appointment in 2017–2018. During the visit, blood samples were collected for plasma FA analysis, and dietary habits were assessed using a short food frequency questionnaire. Information on GDM diagnoses was then later extracted from medical records. Differences in FA profile between GDM cases and non-cases were evaluated using the Mann-Whitney U test.

Results GDM was diagnosed in 127 women (14.9%). Concentrations of saturated fatty acids, monounsaturated fatty acids, polyunsaturated fatty acids (PUFA) n-6, PUFA n-3 and total FA were higher in the women who later developed GDM compared with those who did not (p≤0.05). The medians for total FA were 2898 μg/mL for the women with GDM and 2681 μg/mL for those without GDM. Mean adjusted difference for total FA between the groups was 133 μg/mL (95% CI 33 to 233). Similar results were observed in prepregnancy normal-weight women and overweight women/women with obesity. Overall diet quality in early pregnancy appeared to be lower among the women later diagnosed with GDM.

Conclusion We found that plasma FA profiles in early pregnancy were different for women later diagnosed with GDM compared with those who were not, independent of the women’s body mass index.

  • fatty acids
  • diet
  • pregnancy
  • diabetes
  • gestational

Data availability statement

No data are available. Data are not in a repository.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Significance of this study

What is already known about this subject?

  • Fatty acid (FA) profiles have been associated with gestational diabetes; however, most studies were performed during or after gestational diabetes mellitus (GDM) diagnoses and few have additionally investigated dietary intake.

What are the new findings?

  • Plasma FA profiles in early pregnancy were different for women with and without GDM diagnoses later in pregnancy.

  • Differences in plasma FA profiles were independent of the women’s body mass index.

  • Women that were later diagnosed with GDM appeared to have lower overall diet quality in early pregnancy.

How might these results change the focus of research or clinical practice?

  • Further research is needed on early pregnancy FA profiles as a predicter for GDM.

  • Increased dietary quality during pregnancy should be a clinical focus in GDM prevention.

Introduction

Gestational diabetes mellitus (GDM) is one of the most common pregnancy complications and a strong risk factor for later development of type 2 diabetes (T2D).1 Potential risk factors for GDM include prepregnancy overweight and obesity, age and unhealthy dietary habits.2 With respect to dietary fat intake, high maternal intake of total fat,3 saturated fatty acids (SFA) as well as cholesterol4 have previously been associated with an increased risk of GDM,3 5 whereas intake of polyunsaturated fatty acid (PUFA) and alpha-linolenic acid (ALA) has been associated with reduced risk.6 Circulating FA can derive from both dietary intake7 and catabolism or endogenous synthesis, both of which being partly influenced by glucose homeostasis.8 As an example, free fatty acid (FFA) levels have been associated with altered secretion of insulin as well as insulin resistance (IR), which may progress to diabetes.9–12

In most previous studies on FA profiles in pregnant cohorts, FA analysis was performed during or after GDM diagnoses,12 whereas few studies have investigated the FA profile during the first trimester, before GDM diagnoses.13–17 Their results indicated that women diagnosed later with GDM had higher total concentrations of SFA,13–15 MUFA,13 15 PUFA n-616 and PUFA n-313 and sometimes lower concentrations of PUFA n-613 15 and PUFA n-3.15–17

Our aim was to investigate whether the total plasma FA profile during the first trimester was associated with increased odds of GDM diagnoses. Because overweight and obesity can increase the risk of IR and GDM and is associated with altered FA profiles,12 we also stratified our analysis according to pre-pregnancy body mass index (BMI). Furthermore, self-reported dietary consumption was also compared between the two groups of women.

Subjects and methods

Subjects

Between October 2017 and March 2018, all pregnant women attending routine screening at gestational weeks 11–14 at the Prenatal Diagnostic Unit at The National University Hospital (Reykjavik, Iceland) were invited to become participants in the study. Of the 1684 women scheduled for an appointment, 1350 were eligible by being able to answer an Icelandic Food Frequency Questionnaire (FFQ) in addition to their gestational age being between 11 and 14 weeks. Participant flow chart is shown in figure 1. A total of 1015 women agreed to participate (75%). Most of the women provided a blood sample during this routine visit, at which time additional tubes were drawn for the purpose of FA analysis. A total of 954 participants provided a non-fasting blood sample. We were able to acquire information on GDM diagnoses from medical data for 853 of these women; hence, the entire analysis consisted of 853 women or 84% of the enrolled study participants. The source population and cohort has been described in detail in previous publications.18 19

Figure 1

Participant flow chart. GDM, gestational diabetes mellitus; PREWICE II, Pregnant Women in Iceland II.

Dietary intake and characteristics

During recruitment, the participants answered a short FFQ in electronic format on dietary intake, which also contained questions on age, education, smoking, parity, prepregnancy weight and height. The FFQ assessed the frequency of consumption of 40 food items and beverages and dietary supplement intake, requiring participants to refer to their intake during the previous 3 months. The women chose between 10 potential frequency responses ranging from ‘less than once a month’ to ‘more than 5 times a day’. The FFQ development has previously been described in detail.18 20–22

In brief, the FFQ was pilot-tested and compared with results from a 4-day weighed food record with acceptable correlation for most food groups/items (Spearman’s correlation >0.3).20 Recent publications from our PREWICE II study have described a dose-dependent association between the consumption of dairy and urine iodine concentration,18 wholegrain consumption and alkylresorcinol concentrations (a biomarker for wholegrains)19 and reported intake of vitamin D supplements and 25-OH-D in plasma.22

Measurement of plasma fatty acids in plasma

The voluntary blood samples provided at Landspitali University Hospital in weeks 11–14 by consenting participants were processed within 1 hour after collection to separate plasma from red blood cells and buffy coat via centrifugation at 3000 rpm for 10 min. Plasma was aliquoted into cryotubes and stored in a freezer at −80°C until shipped for FA analysis at the Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden. The total plasma fatty acid (FA) composition was determined according to the method described by Stråvik et al,23 which is a modification of the method described by Masood et al.24 In total, 24 FA were quantified. An internal standard solution (100 µL of 0.1 mg C23:0 methyl ester/mL toluene) was added to 50 µL of thawed plasma samples, and 1.8 mL of acetyl chloride-MeOH solution (10% (v/v) fortified with butylated hydroxytoluene (2.78 µg/mL) was added. Samples were incubated at 70°C for 60 min in a water bath with shaking. Single extraction of fatty acid methyl esters (FAMEs) was carried out by adding 1.5 mL of hexane. The extraction solvent was evaporated using a vacuum concentrator (125 mbar, 30°C, 30 min). FAME was dissolved in 200 µL of hexane before injection into the GC-FID for analysis. Two water blanks (Millipore-purified water) and six quality controls (pooled plasma) were prepared and run together with the study samples for each batch of a maximum of 50 study samples.

FAMEs were separated using the GC-FID system (Thermo Scientific Focus GC, FID detector, Pal GC-xt autosampler, AD-100 H2 generator and a MicroClip XT hydrogen gas alert and a ZA 1500 zero air generator) equipped with a Zebron ZB-FAME column (20 m×0.18 µm ID×0.15 µm). The oven program was as follows: initial 80°C with a 1.5 min hold; ramp: 40°C/min to 160°C, 5°C/min to 185°C with a 0 min hold and then 30°C/min to 260°C with a 0 min hold. The instrumental condition was as follows: nitrogen as the carrier gas, constant flow, carrier flow 1.25 mL/min, Inlet temperature 260°C, split flow 12.5 mL/min, split ratio: 15. Detector temperature 260°C. Gas flow: air 450 mL/min; hydrogen 35 mL/min; makeup gas 10 mL/min. The injection volume was 1 µL. The concentration of FAs in samples was quantified against external standard calibrations made from GLC-462 mixed FAMEs (Nu-Check Prep, Elysian, Minnesota, USA) dissolved in toluene. The external standard included 24 FAMEs ranging from C12:0 to C24:1. An equal amount of internal standard (C23:0 methyl ester) was added to the external standards as added to the study samples, and all analyte peaks were normalized with the peak of the internal standard before calibration.

Gestational diabetes mellitus diagnoses

The GDM diagnoses are based on The International Association of the Diabetes and Pregnancy Study Groups (IADPSG).25 During the first routine maternal care visit in Iceland (around 10 weeks of pregnancy), women at risk of developing GDM, based on age (≥40 years), BMI (≥30 kg/m2), ethnicity and history of diabetes or macrosomia (≥4500 g), are invited to provide a fasting blood sample to measure fasting blood glucose. If results are ≥5.1 mmol/L, the women are diagnosed with GDM or in some cases T2D, if fasting blood sugar is >7 mmol/L or hemoglobin A1c is ≥48 mmol/mol. The rate of women diagnosed with GDM during this early selective screening is unknown in Iceland. However, a study in France also using selective screening and the IADPSG criteria reported an early hyperglycemia rate of 2.3% in a cohort of almost 800 000 women, corresponding to 26.9% of the women diagnosed with hyperglycemia overall.26

Later in pregnancy, usually at weeks 24–28, the women considered at risk of GDM undergo a 2-hour oral glucose tolerance test (OGTT). During the OGTT blood sugar levels are measured during fasting, after 1 hour and after the second hour.27 Previously reported rates of GDM diagnoses in Iceland have been in the range of 11.8%–16%.28 29 Information on GDM diagnoses for this study was gathered from medical records at the National University Hospital. The records did not always differentiate between early detected GDM, GDMA1 (controlled with diet) or GDMA2 (medication needed); therefore, the GDM diagnoses cover all categories.

Statistical analysis

Data are presented as means and SD for normally distributed variables or median and 10th–90th percentiles for skewed distributions. A t-test for equality of means was used to compare the normally distributed variables, and Pearson’s χ2 test was used to compare dichotomous variables. The Mann-Whitney U test for two independent samples was used to compare differences for skewed variables.

For classification of FA groups, the concentrations of different types of FA were combined to determine the total for SFA, MUFA, PUFA n-6, PUFA n-3 as well as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) especially. Ratios of individual and groups of FA were calculated by dividing their concentration by the total FA concentration.

The value for six of the FA was below the limit of quantification in a few cases. For four of the FA 0.1%–0.2% of the participants had levels below quantification, but for myristoleate (14:1) the rate was 3% and for docosatetraenoic (22:4) it was 10%. Values below limit of quantification were set to half of the lowest value quantified for an FA in our data.

A linear regression model to compare means was used to evaluate the adjusted difference in FA types, between the women who were diagnosed later with GDM and those who were not.

Multivariate binary logistic regression was used to evaluate the OR of GDM across quartiles of FA types as well as relative FA. P for trend was evaluated by using the median value in each quartile and modeling the plasma FA and relative FA variables as continuous in the regression model. Results from the regression model are presented as multivariate-adjusted. A Spearman’s correlation was used to investigate correlations between dietary intake and FA. IBM Statistical Package for Social Sciences for Windows, V.24.0 (Armonk, New York, USA), was used to analyze the data. The level of significance was accepted as p<0.05.

Covariates

The covariates included in the adjusted models were selected a priori based on their potential influence on GDM diagnoses.30–34 In the adjusted linear regression model, covariates were age, prepregnancy BMI (kg/m2), weekly weight gain (kg/week) and maternal smoking during pregnancy.

The covariates included in our multivariate adjusted models were age; prepregnancy BMI (kg/m2); weekly weight gain (kg/week); parity (0, 1 and ≥2); maternal smoking during pregnancy (yes/no) and family history of diabetes (yes/no/unknown). In cases when missing values for covariates were low (<5%), they were imputed using the median or the most probable value, which was the case for prepregnancy BMI. When prepregnancy BMI values were missing (2.6% overall, 7.6% of the women diagnosed with GDM), values were imputed using the median or the most probable value based on BMI at the first maternal care visit. For family history of diabetes (missing=15%), missing values were accounted for by use of the missing category for covariate adjustment.

Medical records provided information on gestational age and maternal weight at maternal care visits. From these data, we were able to calculate total weight gain and weight gain per week by dividing total weight gain with the number of weeks between visits.

Self-reported prepregnancy weight and height from the FFQs was used to calculate prepregnancy BMI (kg/m2). A BMI <18.5 kg/m2 was defined as underweight, 18.5–24.9 kg/m2 as normal weight, 25–29.9 kg/m2 as overweight and ≥30.0 kg/m2 as obese.

Results

Characteristics of the participants are presented in table 1.

Table 1

Characteristics of the women in PREWICE II

The mean age of all participants was 30 years, and 44% were nullipara. In total, 59% had a university-level or higher academic education, and 14% smoked before pregnancy. The prevalence of GDM diagnoses was 14.9% (n=127). The women diagnosed with GDM were more likely to be older and have a higher prepregnancy BMI. They were, therefore, also more likely to be overweight or obese prior to pregnancy (59.8%) compared with women without GDM (44.6%).

The concentrations of total and relative SFA, MUFA, PUFA n-6 and n-3 at gestational weeks 11–14 are provided in table 2, whereas the concentrations of individual FAs for women with and without later GDM diagnoses are presented in online supplemental table 1.

Supplemental material

Table 2

Fatty acid concentrations at gestational weeks 11–14 in women with and without GDM diagnosis later in pregnancy, also stratified by BMI*

The total concentration of FA was significantly higher in women diagnosed with GDM, as were the concentrations of all types of FA, except for long n-3 FAs EPA+DHA. When stratified by prepregnancy BMI (<25 kg/m2 vs ≥25 kg/m2), the same tendency toward a higher concentration of total FA and MUFA in women who later were diagnosed with GDM was observed in both groups. When comparing relative FA concentrations of women later diagnosed with GDM with those who were not, MUFA was significantly higher, and PUFA n-6 was significantly lower for the women who later developed GDM. This difference remained after stratifying by prepregnancy BMI (<25 kg/m2 vs ≥25 kg/m2).

Table 3 shows the OR for developing GDM according to quartiles of FA concentrations (total plasma concentration and ratio of total plasma concentration).

Table 3

The associations between quartiles of fatty acid concentrations and ratios with GDM diagnoses

When examining the total concentration, a significant trend toward increased odds of GDM was noted with higher concentrations of MUFA, PUFA n-6 and total FA. However, when FA was expressed as a ratio of total plasma concentration, increased odds for GDM were seen with a higher ratio of MUFA. The odds of GDM were lower in the highest quartile of PUFA n-6 ratio compared with the lowest quartile, although p for trend did not reach significance. No association was observed between PUFA n-3 and GDM; neither for total concentration nor the ratio of total plasma PUFA n-3.

Table 4 shows a comparison of dietary intake in early pregnancy between the women who were later diagnosed with GDM and those who were not.

Table 4

FFQ reported weekly intake of foods at 11th–14th week of pregnancy*

Women who were later diagnosed with GDM had a significantly lower intake of fatty fish and skimmed milk compared with women who did not develop GDM. They also tended to use a higher proportion of saturated fat but a lower proportion of vegetable oil when cooking (50% vs 29% p=0.06) versus (50% vs 71% p=0.06), respectively. Reported intake of omega-3 supplements also tended to be lower for women who were later diagnosed with GDM, with differences close to significance. Other differences observed in dietary intake between the two groups included a lower frequency intake of whole grains, fruit juice, beans, nuts, seeds and coffee during early pregnancy in women who later were diagnosed with GDM.

There were no clear correlations between dietary intake and FA concentrations except in the case of fatty fish and omega-3 supplements, both of which correlated positively with concentrations of PUFA n-3 and EPA+DHA (online supplemental tables 2 and 3) pregnancy (59.8%) compared with women without GDM (44.6%).

One limitation of our study is that we did not have exact data at the time of diagnoses for our GDM cases. As subjects with known risk factors are more likely to be diagnosed early (online supplemental table 4), we examined the concentrations of total and relative SFA, MUFA, PUFA n-6 and n-3, excluding all women with known risk factors that prompt early GDM screening (BMI ≥30 kg/m2, age ≥40 years and parity ≥1, to exclude all women with previous GDM/macrosomia).

Discussion

We found that the total concentrations of plasma FAs, as well as total MUFA and MUFA ratio measured at 11th–14th week of pregnancy, were significantly higher and PUFA n-6 ratios lower in women who were later diagnosed with GDM, independent of the women’s BMI. The fact that stratifying by BMI did not alter our results is important to note because obesity, IR and FA profiles in GDM are strongly inter-related.12

In most previous studies on FA profiles in pregnancy cohorts, FA analysis was performed during or after GDM diagnoses.12 These studies found that SFA concentrations appeared to be higher in women diagnosed with GDM when compared with a control group and that PUFA n-6 and PUFA n-3 concentrations were both lower in women with GDM.12 A recent meta-analysis similarly reported that women with GDM had higher total concentrations of FFAs in the second and third trimester compared with women without GDM, with concentrations decreasing as pregnancy progressed.35

The results of the few studies that have investigated FA concentrations in early pregnancy, prior to GDM diagnosis, agree with the present study’s findings, showing higher total concentrations of SFA13–15 and MUFA13 15 as well as PUFA n-616 and PUFA n-313 in women who were later diagnosed with GDM compared with those who were not. On the other hand, some of these studies reported lower PUFA n-613 15 and PUFA n-315–17 concentrations in women who later received a GDM diagnosis.

FA profiles have been proposed as a means of predicting later T2D diagnosis in non-pregnant populations, where higher relative concentrations of FFA,36 PUFA n-6,37 MUFA36 37 and SFA36 37 have been associated with an increased risk of impaired glucose tolerance and T2D risk. These studies have reported similar results as observed for circulating FA in the pregnant population in our study. Associations between higher PUFA n-6 relative to total FA and increased insulin sensitivity have also been found in previous studies in non-pregnant populations,36 which is in accordance with the results seen in the present pregnancy cohort.

The FA profiles may be a result of differences in intake or absorption of both carbohydrates and fat.12 Plasma MUFA and SFA concentrations do not only represent dietary intake as FAs can be synthesized endogenously, mainly from carbohydrates, which is then referred to as de novo lipogenesis.38 39 Other factors that could influence FA profile are FA synthesis and incorporation of FA into cell membranes. It remains unclear how the plasma total FA profile is physiologically associated with diabetes. Lipogenesis is stimulated by insulin and suppressed by the hormones glucagon and epinephrine. Some studies have suggested that higher FFA may alter insulin signalling, secretion and glucose production.40 41 It is, therefore, possible that an abnormal increase of insulin in the blood may lead to higher FA concentrations and vice versa.

Some differences were found in food consumption in the first trimester between women with and without GDM diagnoses, and it is not clear how this difference might be reflected in the plasma FA profile. It could have been expected to see a difference in EPA and DHA concentrations between the two groups because women who later were diagnosed with GDM had a less frequent intake of both fatty fish and omega-3 supplements. However, because FA can be synthesized endogenously from excess carbohydrates, the overall quality of the diet, including carbohydrate quality and amount consumed, might explain some of the difference observed in FA concentration between the two groups, overall carbohydrate quality being one.19 21 42

As we have previously reported in this cohort, the women who were later diagnosed with GDM had lower quality of carbohydrate intake as they had a lower intake of wholegrains (estimated by biomarkers).19 Intake of soft drinks also tended to be greater for women diagnosed with GDM as well as use of saturated fat for cooking, even though the difference was not statistically significant. This could suggest that the overall diet quality in early pregnancy was lower among the women later diagnosed with GDM.

In summary, a stronger association has been observed between plasma FA and diabetes risk compared with dietary intake estimates.37 It is important to note that different sources of FA measurements can represent varying dietary intake periods, such as adipose tissue (long-term FA intake 1–1.5 years), skeletal muscle cells, erythrocytes (120 days), serum, total plasma, phospholipids (1–2 weeks12), cholesteryl esters and FFAs.37 Erythrocytes have been claimed as a preferable option to evaluate differences in recent FA intake, at least regarding PUFA n-3.39 However, FA from erythrocytes and plasma have been found to correlate.37 In our study, we analyzed total plasma FA, which includes the FAs from cholesteryl esters, phospholipids, triacyclglycerols and FFA thought to represent very recent intake (1–2 weeks), whereas the answers in the FFQ covered intake during the previous 3 months.

It might be considered a limitation that our participants were not fasting. This might have resulted in lower concentrations of FFA because they are reduced in response to higher insulin levels postmeals43 and makes a comparison with other studies challenging. However, comparisons of our results with other studies may be affected by the non-fasting state of our participants and by the difference in methodology used when analysing FA profiles. Another limitation is that we did not have exact information on date of diagnoses, thereby part of our GDM cases may have received their diagnosis at or close to drawing of the blood sample in which FA profiles were quantified. However, we did see a similar trend in the results (online supplemental 4) when we ran our analysis excluding all women with known risk factors that prompt early GDM screening (BMI, age and previous GDM/macrosomia). We therefore believe that this limitation did not majorly affect our main results.

The strength of our study is that we analyzed plasma FA early in pregnancy, in addition to acquiring subjective data on dietary intake. Our study also features prospectively collected data from a large sample size of 853 pregnant women with high participation rate (75%). We report both total concentrations of all FA subgroups as well as their ratios. We adjusted for BMI, and stratified our FA results by BMI, a process that, to our knowledge, has not been previously performed in a pregnancy cohort.

In conclusion, we found that women who were later diagnosed with GDM had a higher concentration of total plasma FA, total MUFA and MUFA ratios as well as lower PUFA n-6 ratios in early pregnancy, independent of the women’s BMI, compared with women who remained free of GDM. The women who were not diagnosed with GDM also tended to have better diet quality in early pregnancy. These results suggest that FA biomarkers in early pregnancy may predict GDM. However, further studies are required to confirm this hypothesis.

Data availability statement

No data are available. Data are not in a repository.

Ethics statements

Ethics approval

The study was approved by the National Bioethics Committee (VSN-17-057-S1) and the Medical Directorate of Landspitali University Hospital (LSH 5-17). Written consents were obtained from the participants.

Acknowledgments

The authors would like to thank the staff at the Prenatal Diagnostic Unit at the National University Hospital in Iceland for their hospitality and positive attitudes, which contributed significantly to the recruitment of participants for this study.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors IG, ITH and HH designed the research. EAT and IG conducted the research. RL, LH, ITH and HH provided essential materials. EAT, TIH and IG analyzed data and performed statistical analysis. EAT, TIH and IG wrote the paper. EAT, TIH and IG had primary responsibility for final content. All authors read, edited and approved the final manuscript.

  • Funding This study was funded by the University of Iceland Research Fund (HI18080069) and the Science Fund of The National University Hospital in Iceland (798381 and 893911).

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.