Skip to main content

Peripheral thyroid hormone levels and hepatic thyroid hormone deiodinase gene expression in dairy heifers on the day of ovulation and during the early peri-implantation period

Abstract

Background

Before the onset of fetal thyroid hormone production, the transplacental delivery of maternal thyroid hormones is necessary for embryonic and fetal development. Therefore, the adaptation of maternal thyroid hormone metabolism may be important for pregnancy success and embryo survival. The aims of this study were to determine the thyroid hormone levels during the early peri-implantation period until day 18 and on the day of ovulation, to determine whether pregnancy success is dependent on a “normothyroid status” and to determine whether physiological adaptations in maternal thyroid hormone metabolism occur, which may be necessary to provide sufficient amounts of biologically active T3 to support early pregnancy. Therefore, blood samples obtained on the day of ovulation (day 0) and days 14 and 18 of the Holstein–Friesian heifers (n = 10) during the respective pregnant, non-pregnant and negative control cycles were analyzed for thyroid-stimulating-hormone (TSH), thyroxine (T4) and triiodothyronine (T3). Liver biopsies (day 18) from pregnant and respective non-pregnant heifers were analyzed for mRNA expression of the most abundant hepatic thyroid hormone deiodinase (DIO1) by real time qPCR.

Results

Although liver DIO1 mRNA expression did not differ between the pregnant and non-pregnant heifers on day 18, the serum concentrations of TSH and T3 on day 18 were higher in non-pregnant heifers compared to pregnant heifers (P < 0.05). Moreover, T3 decreased between day 0 and 18 in pregnant heifers (P < 0.001).

Conclusions

In conclusion, no associations between thyroid hormone patterns on day 18 and pregnancy success were detected. During the early peri-implantation period, TSH and T3 may be affected by the pregnancy status because both TSH and T3 were lower on day 18 in pregnant heifers compared to non-pregnant dairy heifers. In further studies, the thyroid hormone axis should be evaluated throughout the entire gestation to confirm these data and identify other possible effects of pregnancy on the thyroid hormone axis in cattle.

Background

Thyroid hormones play a key role in fertility, conceptus development, tissue differentiation and fetal growth [15]. In humans and rats, the conceptus starts to produce thyroid hormones during the 2nd trimester of pregnancy [1, 5]. Afterwards, the maternal production of thyroid hormones is less essential for fetal development. However, before the onset of fetal thyroid hormone production, the transplacental delivery of maternal thyroid hormones is necessary [6, 7]. The importance of thyroid hormones for the bovine embryo was substantiated by results derived from in vitro studies that show that supplementation with T3 and T4 was associated with an improved viability of bovine embryos [8]. In addition, the hormonal activity of the thyroid gland has an important role in the metabolism of cells, lipids and carbohydrates as well as in the lactation course [9]. According to these data, on the one hand, pregnancy success might be dependent on a “normothyroid status,” and on the other hand, physiological adaptations in maternal thyroid hormone metabolism during early pregnancy in cattle might be necessary to provide sufficient amounts of biologically active T3 to support early pregnancy, including the growth and viability of the conceptus.

The thyroid hormone axis consists of the hypothalamic thyrotropin releasing hormone (TRH) and thyroid stimulating hormone (TSH), which are released by the pituitary. Under the control of TSH, the thyroid gland mainly produces the inactive pro-hormone T4, which can be activated locally through tissue-specific deiodinases (DIO) and deiodinated into the active hormone T3 [10]. By modulating DIO expression levels in different organs, e.g., the liver, adequate T3 concentrations can be achieved even during iodine deficiency [11]. In the liver, DIO1 is expressed in high amounts [11], and its expression level influences the T3 blood concentration. Furthermore, the hepatic expression level of DIO1 is known to be altered by steroid hormones, such as androgens in rats [12], and was previously shown to have higher expressed during the late pregnancy of pluriparous Holstein–Friesian cows than immediately postpartum [13]. Whether sexual steroid hormone patterns or early pregnancy may have an influence on hepatic DIO1 expression has not been examined in cattle until now. McCann and Reimers [14] reported that cyclic heifers had higher serum concentrations of T3 and T4 during estrus compared to diestrus, whereas Elecko et al. [15] showed that the blood concentrations of T4 decreased under estradiol benzoate infusion. However, the effect of early pregnancy on the maternal thyroid hormone axis and hepatic DIO1 expression has not been previously examined in heifers. Therefore, the aim of this study was to examine whether the maternal thyroid hormone axis differs between early pregnant and respective non-pregnant cycles of healthy Holstein–Friesian heifers.

Methods

The experimental setup was previously published by Meyerholz et al. [16]. All of the interventional procedures were performed according to the German legislation on animal welfare and approved by the Lower Saxony Federal State Office for Consumer Protection 279 and Food Safety under the reference number AZ 33.14-42502-04-12/0744. Thirty Holstein–Friesian heifers aged 14.1 ± 1.6 months with a mean body weight of 337.8 ± 23.7 kg were used. The heifers were clinically healthy, and the health status was monitored throughout the entire experimental period. Briefly, three cycles of each of the Holstein–Friesian heifers (n = 10) were evaluated. Each cycle started on the day of ovulation (day 0) and ended on day 18. During the first “negative control” (NC) cycle, no embryo was transferred, whereas throughout the following cycles, the embryo transfer was performed on d 6 ± 12 h. Pregnancy was confirmed by high progesterone concentrations on day 18 and recovery of a trophoblast tissue by uterine flushing. Conversely, non-pregnant heifers showed estrus symptoms on day 18.

Blood samples and liver biopsies

Blood samples were obtained from the jugular vein on days 0, 14 and 18, and due to the circadian and ultradian rhythmicity of thyroid hormones, the samples were always collected in the morning between 7:00 and 9:00 a.m. The 10 mL serum samples were collected in EDTA tubes (Sarstedt, Nümbrecht, Germany), maintained at room temperature (RT) for approximately 2 h and centrifuged for 10 min (1500×g, RT). After centrifugation, the serum and plasma was stored at −20 °C until further analysis. Liver biopsy specimens (10 mm × 3 mm) were collected on day 18 in sterile Eppendorf cups, which were immediately frozen in liquid nitrogen and stored at −80 °C until mRNA extraction. The relative abundance of hepatic DIO1 in comparison with the housekeeping genes GAPDH and RPS9 (Table 1) was measured by real-time PCR as previously described [17].

Table 1 The quantitative real-time PCR primers for DIO1 and two housekeeping genes in the liver biopsy specimens of pregnant and non-pregnant heifers on day 18

Blood parameters

TSH [ng/ml]

For the ELISA, an antibody targeted against bovine TSH (anti-bovine TSH, 1:10 pre-diluted, AFP-642482Rb) was obtained from the National Hormone and Peptide Program (NHPP, National Hormone and Peptide Program, NIDDK and Dr. Parlow) and was diluted and used at a final dilution of 1:2500. The standard curve ranged from 0.2 to 100 ng/ml bovine TSH (AFP-8755B, obtained from the NHPP, NIDDK and Dr. Parlow) and was dissolved in peptide buffer. The standards, controls and plasma samples (in triplicate) were added to a microtiter plate coated with the antibody, and the plate was then incubated for 24 hat RT. After washing the plate, biotin-labeled TSH [biotin labeled TSH (AFP-8755B, obtained from the NHPP, NIDDK and Dr. Parlow)] was added to all of the wells and incubated for 3 h. Then, a streptavidin–horseradish peroxidase solution (Sigma Aldrich, St. Louis, MO, USA) was added, the substrate (containing tetramethylbenzidine, Sigma Aldrich, St. Louis, MO, USA) was pipetted after washing, and the reaction was stopped after 15 min by adding sulfuric acid (2M; Sigma Aldrich, St. Louis, MO, USA). The optical density was obtained at a wavelength of 450 nm, and the concentrations were calculated using Magellan software with the cubic spline modus (Magellan 3.11, Dortmund, Germany). The intra-assay CV was determined by measuring one bovine sample 20 times, and the result was 15.4. The lowest detection limit of the ELISA was 0.6 ng/ml, which was determined by using the last detectable concentration of a bovine serum sample that was serially diluted.

Thyroxine and triiodothyronine [nmol/l]

The total serum T4 and T3 concentrations were determined using radioimmunoassays (TOTAL T4 RIA KIT, IM1447-IM3286 and TOTAL T3 RIA KIT, IM1699-IM3287, Immunotech, Beckman Coulter, CA, USA). The assays were performed according to the manufacturer instructions, and the samples were analyzed in duplicate. The intra-assay CV was determined by analyzing the bovine serum samples within one test run; the CV was 6.2 % for T4 and 6.3 % for T3, and the inter-assay CV was calculated by analyzing the same bovine serum sample in ten different test runs. The inter-assay CV was 8.6 % for T4 and 7.7 % for T3. The lowest detection limits were 13.0 nmol/l for T4 and 0.3 nmol/l for T3.

Statistical analyses

The Statistical Analysis System software (9.3, SAS Inc., Cary, NC, USA) was used to perform all of the statistical analyses. A linear mixed-effect model for repeated measures was chosen to work out the dependence of measured bovine thyroid hormone concentrations (T4, T3 and TSH) on varying groups (p, np, NC) and time (days 0, 14, 18). To meet the primary assumption of the restricted maximum likelihood method (REML), data was tested for normal distribution (PROC UNIVARIATE NORMAL) for each of the nine levels of group x time interaction. As data was predominantly normally distributed, PROC MIXED and method REML with fixed effects “group”, ‘time” and “group × time” and linear regression effect “age” were performed.

The age effect components in the linear mixed-effect model were estimated with the following intercepts and slopes: TSHest(age) = 8.1496–0.01016*age (P < 0.001), T4est(age) =−2.1811 + 0.08157*age (P = 0.019), T3est(age) = 0.7432 + 0.001949*age P = 0.017). The results throughout the manuscript are presented as mean ± standard error of the mean (SEM). Concerning T4, T3 and TSH concentrations, data was adjusted according to the linear regression component of the model (age). Therefore, differences were formed: “hormone concentration at each data point” minus “estimated concentration at age at sample day” resulting in “∆ hormone concentrations” (Fig. 1). Hepatic mRNA expression was compared using a paired Student’s t test (PROC TTEST PAIRED) and results are presented as mean ± SEM as well. For all procedures, the statistical significance was pre-established at P < 0.05. P values between P > 0.05 and P < 0.10 were considered statistical tendencies. F and P values for the fixed effects are presented below each graph (Fig. 1), significant differences between pregnant (p) negative control (NC), and non-pregnant (np) cycles are indicated by different letters (Fig. 1) whereas differences along the time points within each group are described in the text.

Fig. 1
figure 1

Age adjusted A thyroid-stimulating-hormone (∆TSH), B thyroxine (∆T4) and C triiodothyronine (∆T3) blood concentrations (mean ± SEM) of pregnant (p) heifers and respective non-pregnant (np) and negative control cycles (NC); the significant differences between the groups on certain days are marked with different letters. The results of the mixed–model are shown in the table below each graph

Results

Blood parameters

The thyroxine concentrations were higher in non-pregnant compared to pregnant cycles (P = 0.046) on day 14 but comparable at all of the other tested days. Triiodothyronine concentrations decreased in pregnant heifers between days 0 and 14 (P = 0.027) as well as between days 0 and 18 (P = 0.008). Interestingly, in non-pregnant heifers, T3 was higher on d 18 compared to pregnant heifers (P = 0.013). TSH concentrations were higher in non-pregnant compared to pregnant heifers on day 18 (P = 0.008) and during the NC (P < 0.001) (Fig. 1A–C).

Hepatic mRNA expression

The relative abundance of DIO1 mRNA expression on day 18 was comparable between pregnant (ΔCT 33.32 ± 4.35) and non-pregnant heifers (ΔCT 27.02 ± 7.78).

Discussion

The first aim of the present study was to detect the possible associations between pregnancy success and the thyroid status on the day of ovulation; however, no differences in TSH, T3 or T4 levels were detected on the day of ovulation. Second, the maternal thyroid hormone axis was evaluated during early pregnancy in dairy heifers. TSH concentrations were lower in pregnant heifers compared to non-pregnant heifers on day 18; however, the negative controls showed even lower TSH levels compared to pregnant heifers. Therefore, an effect of early pregnancy on TSH is unlikely. However, it should be noted that daily blood sampling may not be adequate to evaluate the TSH concentration in heifers as thyroid hormones have a circadian and ultradian rhythmicity in bovine plasma [18]. In future studies, the TSH concentration and release should be evaluated during the first month of pregnancy because the chosen time period in this study could have been too short to observe the effect of thyroid axis adaptions during pregnancy in cattle. Conversely, collecting more blood samples per day or 24 h measurements could help identify the possible effects of pregnancy on TSH release.

Moreover, due to the experimental setup, the original data could have reflected the growth and pregnancy effects at the same time, indicating that NC differs from np cycles because the NC cycle was the first cycle to be performed during the experiment and the heifers were younger, whereas during the np and p cycles, the animals were comparable in age, but the maintenance of early pregnancy led to differences (data not shown). To minimize the effect of age on the data as a confounding variable, the regression factor “age” was calculated, hormone concentrations are presented as adjusted values and the effect was also statistically considered.

T4 concentrations tended to be lower in p compared to np on day 14. These results were unexpected, as studies in humans show elevated T4 levels during pregnancy [19, 20], which were possibly due to the thyroid stimulating effects of human chorionic gonadotropin (hCG) [21]. In cattle, a comparable effect cannot be indicated by the results of the present study.

Hepatic DIO1 activity did not differ between p and np, irrespective of the higher T3 serum concentrations in non-pregnant heifers on day 18. Cann and Raimers [14] detected lower T4 and T3 levels during the luteal phase in heifers compared to the follicular phase [14]. In the present study, the T3 concentration decreased from days 0 to 18, which indicated a corresponding result concerning the steroid hormone influence on T3; however, no differences between p and np were detected concerning T4 concentrations in the present study. It would be interesting to determine in future studies whether the decreasing T3 concentration is originating from enhanced inactivation of T3 into reverse T3. Those inactive forms of T3 were not detected by the T3 assay used in the present study, and a separate analysis would be of interest in future studies in cattle. Type 2 deiodinase (DIO2) [22] was shown to be more highly expressed in the placental tissues of rats and humans during early pregnancy and could result in an inactivation of T3 to reverse-T3 or 3.3′-diiodothyronine [22, 23]; therefore, DIO2 could represent a possible explanation for the lower T3 concentration in pregnant heifers compared to non-pregnant heifers. Placental DIO2 measurement was not addressed in the present study but should be determined in future studies. In rats, pregnancy results in a decrease of total T4 and T3 concentrations throughout gestation [24], which might be in accordance with the present study; however, the thyroid hormone concentrations should be followed for a longer time period of pregnancy to substantiate the present preliminary data.

Another possible explanation for the lower T3 concentrations during early pregnancy could be an enhanced renal iodine clearance resulting from increased blood flow to the kidney and elevated glomerular filtration as observed in humans [25]. However, data concerning pregnancy-dependent elevated urinary excretion of iodine remain contradictory [2628] and could be addressed in heifers to clarify the causes for T3 changes detected in the present study.

Conclusions

In conclusion, TSH and T3 were lower on day 18 in pregnant compared to non-pregnant dairy heifers and indicate an effect of pregnancy on the thyroid gland axis. However, the underlying cause of the lower T3 concentrations remains unclear as at least the expression of hepatic DIO1, the most abundant deiodinase also responsible for differences in differenced in blood T3, was comparable. Based on our data, it seems interesting to focus on thyroid hormone metabolism in high producing dairy cows.

References

  1. Cassar-Malek I, Picard B, Kahl S, Hocquette JF. Relationships between thyroid status, tissue oxidative metabolism, and muscle differentiation in bovine fetuses. Domest Anim Endocrinol. 2007;33:91–106.

    Article  CAS  PubMed  Google Scholar 

  2. Brent GA. Tissue-specific actions of thyroid hormone: insights from animal models. Rev Endocr Metab Disord. 2000;1:27–33.

    Article  CAS  PubMed  Google Scholar 

  3. Stavreus Evers A. Paracrine interactions of thyroid hormones and thyroid stimulation hormone in the female reproductive tract have an impact on female fertility. Front Endocrinol (Lausanne). 2012;3:50.

    Google Scholar 

  4. Kota SK, Gayatri K, Jammula S, Meher LK, Kota SK, Krishna SVS, Modi KD. Fetal endocrinology. Indian J. Endocrinol Metab. 2013;17:568–79.

    Google Scholar 

  5. Van Herck SLJ, Geysens S, Delbaere J, Darras VM. Regulators of thyroid hormone availability and action in embryonic chicken brain development. Gen Comp Endocrinol. 2013;190:96–104.

    Article  PubMed  Google Scholar 

  6. Obregon MJ, Mallol J, Pastor R, Morreale de Escobar G, Escobar del Rey F. L-thyroxine and 3,5,3′-triiodo-L-thyronine in rat embryos before onset of fetal thyroid function. Endocrinology. 1984;114:305–7.

    Article  CAS  PubMed  Google Scholar 

  7. Owens JA. Endocrine and substrate control of fetal growth: placental and maternal influences and insulin-like growth factors. Reprod Fertil Dev. 1991;3:501–17.

    Article  CAS  PubMed  Google Scholar 

  8. Ashkar FA, Semple E, Schmidt CH, St John E, Bartlewski PM, King W. Thyroid hormone supplementation improves bovine embryo development in vitro. Hum Reprod. 2010;25:334–44.

    Article  CAS  PubMed  Google Scholar 

  9. Fiore E, Piccione G, Giansella M, Praticó V, Dara S, Morgante M. Serum thyroid hormone evaluation during transition periods in dairy cows. Arch Anim Breed. 2005;58:403–6.

    Article  Google Scholar 

  10. Pezzi C, Accorsi PA, Vigo D, Govoni N, Gaiani R. 5′-Deiodinase activity and circulating thyronines in lactating cows. J Dairy Sci. 2003;86:152–8.

    Article  CAS  PubMed  Google Scholar 

  11. Lavado-Autric R, Calvo RM, de Mena RM, de Escobar GM, Obregon MJ. Deiodinase activities in thyroids and tissues of iodine-deficient female rats. Endocrinology. 2013;154:529–36.

    Article  CAS  PubMed  Google Scholar 

  12. Sosic-Jurjevic B, Filipovic B, Renko K, Ajdzanovic V, Manojlovic-Stojanoski M, Milosevic V, et al. Orchidectomy of middle-aged rats decreases liver deiodinase 1 and pituitary deiodinase 2 activity. J Endocrinol. 2012;215:247–56.

    Article  CAS  PubMed  Google Scholar 

  13. Piechotta M, Kedves K, Araujo MG, Hoeflich A, Metzger F, Heppelmann M, et al. Hepatic mRNA expression of acid labile subunit and deiodinase 1 differs between cows selected for high versus low concentrations of insulin-like growth factor 1 in late pregnancy. J Dairy Sci. 2013;96:3737–49.

    Article  CAS  PubMed  Google Scholar 

  14. McCann JP, Reimers TJ. Effects of obesity on insulin and glucose metabolism in cyclic heifers. J Anim Sci. 1986;62:772–82.

    Article  CAS  PubMed  Google Scholar 

  15. Elecko J, Bekeová E, Choma J, Maracek I. Concentration of progesterone and thyroxine in the luteal phase of the cycle in heifers after administration of estradiol benzoate. Vet Med (Praha). 1982;27:705–14.

    CAS  Google Scholar 

  16. Meyerholz MM, Mense K, Lietzau M, Kassens A, Linden M, Knaack H, et al. Serum IGFBP4 concentration decreased in dairy heifers towards day 18 of pregnancy. J Vet Sci. 2015;16:413–21.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Mense K, Meyerholz M, Gil Araujo M, Lietzau M, Knaack H, Wrenzycki C, et al. The somatotropic axis during the physiological estrus cycle in dairy heifers-eEffect on hepatic expression of GHR and SOCS2. J Dairy Sci. 2015;98:2409–18.

    Article  CAS  PubMed  Google Scholar 

  18. Bitman J, Kahl S, Wood DL, Lefcourt AM. Circadian and ultradian rhythms of plasma thyroid hormone concentrations in lactating dairy cows. Am J Physiol. 1994;266:1797–803.

    Google Scholar 

  19. Krassas GE, Poppe K, Glinoer D. Thyroid function and human reproductive health. Endocr Rev. 2010;31:702–55.

    Article  CAS  PubMed  Google Scholar 

  20. Hegedüs L, Karstrup S, Rasmussen N. Evidence of cyclic alterations of thyroid size during the menstrual cycle in healthy women. Am J Obstet Gynecol. 1986;155:142–5.

    Article  PubMed  Google Scholar 

  21. De Escobar GM, Ares S, Berbel P, Obregón MJ, del Rey FE. The changing role of maternal thyroid hormone in fetal brain development. Semin Perinatol. 2008;32:380–6.

    Article  PubMed  Google Scholar 

  22. Galton VA, Martinez E, Hernandez A, St Germain EA, Bates JM, St Germain DL. The type 2 iodothyronine deiodinase is expressed in the rat uterus and induced during pregnancy. Endocrinology. 2001;142:2123–8.

    CAS  PubMed  Google Scholar 

  23. Chan S, Kachilele S, Hobbs E, Bulmer JN, Boelaert K, McCabe CJ, Driver PM, Bradwell AR, Kester M, Visser TJ, Franklyn JA, Kilby MD. Placental iodothyronine deiodinase expression in normal and growth-restricted human pregnancies. J Clin Endocrinol Metab. 2003;88:4488–95.

    Article  CAS  PubMed  Google Scholar 

  24. Calvo R, Obregon MJ, Ruiz de Oña C, Ferreiro B, Escobar Del Rey E, Morreale de Escobar G. Thyroid hormone economy in pregnant rats near term: a “physiological” animal model of nonthyroidal illness? Endocrinology. 1990;127:10–6.

    Article  CAS  PubMed  Google Scholar 

  25. Moleti M, Trimarchi F, Vermiglio F. Thyroid physiology in pregnancy. Endocr Pract. 2014;20:589–96.

    Article  PubMed  Google Scholar 

  26. Dworkin HJ, Jacquez J, Beierwaltes WH. Relationship of iodine ingestion to iodine excretion in pregnancy. J Clin Endocrinol Metab. 1966;26:1329–42.

    Article  CAS  PubMed  Google Scholar 

  27. Smyth PP. Variation in iodine handling during normal pregnancy. Thyroid. 1999;9:637–42.

    Article  CAS  PubMed  Google Scholar 

  28. Caron P, Hoff M, Bazzi S, Dufor A, Faure G, Ghandour I, et al. Urinary iodine excretion during normal pregnancy in healthy women living in the southwest of France: correlation with maternal thyroid parameters. Thyroid. 1997;7:749–54.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

MMM was involved in the realization of the animal model, acquisition of the data, analysis and interpretation of the data, statistical analyses and preparation of the manuscript. KM participated in the acquisition of the data. ML performed the statistical analysis of the data. MR, OS, HJS and MH were involved in the conception and design of the study. OS also participated in the realization of the animal model. MS was responsible for the study conception and design and in the realization of the animal model, interpretation of the data, and preparation of the manuscript. Moreover, all of the authors critically revised the manuscript and provided approval for this version to be published. The authors agree to account for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge the support of Zoetis and the iPUD project under an EMIDA ERA-Net award to Professor Martin Sheldon, Swansea University, UK. Moreover, we would like to thank the animal keepers of the Clinic for Cattle for taking care of the cows and assisting with the project as well as Martina Baumgarten, Angela Jordan and Hendrike Knaak for technical help.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Marion Schmicke.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Meyerholz, M.M., Mense, K., Linden, M. et al. Peripheral thyroid hormone levels and hepatic thyroid hormone deiodinase gene expression in dairy heifers on the day of ovulation and during the early peri-implantation period. Acta Vet Scand 58, 52 (2015). https://0-doi-org.brum.beds.ac.uk/10.1186/s13028-016-0231-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13028-016-0231-6

Keywords