Skip to main content

Chronic cortisol exposure in early development leads to neuroendocrine dysregulation in adulthood



Chronic early life stress can affect development of the neuroendocrine stress system, leading to its persistent dysregulation and consequently increased disease risk in adulthood. One contributing factor is thought to be epigenetic programming in response to chronic cortisol exposure during early development. We have previously shown that zebrafish embryos treated chronically with cortisol develop into adults with constitutively elevated whole-body cortisol and aberrant immune gene expression. Here we further characterize that phenotype by assessing persistent effects of the treatment on cortisol tissue distribution and dynamics, chromatin accessibility, and activities of glucocorticoid-responsive regulatory genes klf9 and fkbp5. To that end cortisol levels in different tissues of fed and fasted adults were measured using ELISA, open chromatin in adult blood cells was mapped using ATAC-seq, and gene activity in adult blood and brain cells was measured using qRT-PCR.


Adults derived from cortisol-treated embryos have elevated whole-body cortisol with aberrantly regulated tissue distribution and dynamics that correlate with differential activity of klf9 and fkbp5 in blood and brain.


Epidemiological studies have shown that early life stress (ELS) can increase later life risk of developing a variety of physical and mental health problems, many of which are linked to chronic inflammation [1,2,3]. This long-term ‘developmental programming’ has been hypothesized to in part be an epigenetic effect of chronic exposure to elevated cortisol levels [4, 5]. Consistent with this, chronic glucocorticoid (GC) exposure, such as occurs with chronic stress, genetic diseases such as Cushing’s syndrome, or with extended GC therapy, promotes development of metabolic and inflammatory disease [6]. Developmental programming in response to ELS can be adaptive to the extent that it tunes the responsiveness of the stress system to the stressfulness of the environment, but can be maladaptive when there is a mismatch between the adult environment and that which was encountered during earlier development [1, 7]. The cost of such programming can be an increased allostatic load that promotes unhealthy aging [8].

Zebrafish are emerging as an excellent model system for experimental studies of stress-induced developmental programming [9,10,11,12,13,14,15,16,17,18]. We recently showed that zebrafish embryos treated chronically with 1 μM cortisol develop into adults that maintain elevated whole body cortisol levels and aberrantly express pro-inflammatory genes, with higher basal expression levels in peripheral tissues but a blunted response to tailfin injury or lipopolysaccharide injection [13], indicating that the treatment has long-term effects on the neuroendocrine stress axis and GC-regulated gene expression. Those results are extended by the experiments reported here.

Main text


Adults derived from cortisol-treated embryos display aberrant cortisol tissue distribution and dynamics

Zebrafish embryos were treated with cortisol or vehicle as described previously [13], then raised to adulthood (5 + months), at which time cortisol levels in different tissues were measured in fish that had recently been fed or that had been fasted for 24 h (Fig. 1a, b). As we observed previously [13], the fish in the cortisol treated group had higher whole-body cortisol levels than their untreated siblings, which in fed fish correlated with higher levels in kidney (which when dissected out likely contains the cortisol-producing interrenal gland) and blood (Fig. 1b). However, in fish that had been fasted for 24 h, cortisol levels were lower in blood of the cortisol-treated group compared to controls, but higher in the brain, skin and gut (Fig. 1b). Plotting of the ratios of cortisol levels in fed and fasted fish reveals that while 24 h of fasting did not affect brain cortisol levels in control fish, it caused brain cortisol levels to quadruple in the treated fish (Fig. 1b, right panel). A three-way ANOVA (Additional file 1) showed significant variance in cortisol levels across different tissues is (p = 0.005), as well as a significant interaction effect for tissue type and fed state (p = 0.035). Post-hoc two-way ANOVA indicated that fed state significantly affected cortisol levels in kidney, blood, brain, and muscle (Additional file 1).

Fig. 1
figure 1

Embryos treated chronically with cortisol develop into adults with aberrant cortisol levels, tissue distribution, and dynamics. a Schematic of experimental design; b Cortisol levels in different tissues of fed and fasted (24 h) fish (left panel) and the inferred change in those levels induced by fasting (right panel). Each measurement was taken from pooled tissues of 6 fish, with equal numbers of males and females. Error bars are the standard deviations of technical replicates. c Cortisol levels in blood draws from single fish at night and the following morning, with the amount of time since the last feeding indicated. Thin lines between data points indicate repeat measurements from the same fish. Averages ± standard errors of the mean (SEM) are also shown. d Fasting cortisol levels in blood of individual fish (data points) following 8 h exposure to 1 μM dexamethasone (DEX) or vehicle. The averages ± SEM for each group are also indicated. The fact that the levels are somewhat lower than those shown in (b) may reflect circadian fluctuation (these samples were taken in late afternoon whereas those shown in panel (b) were taken mid-morning, ~ 8 h earlier). e Relative expression of pomca in brain tissue in the late morning, as well as later the same day following 8 h exposure to DEX or vehicle. The bars indicate averages ± the standard deviation of 3 qPCR readings (technical replicates). For each measurement brain tissue was pooled from 5–6 fish, mixed males and females in equivalent proportions within each comparison group (controls 1 male and 4–5 females; cortisol-treated 2 males and 3 females)

In multiple experiments, fasted adults derived from cortisol-treated embryos consistently displayed lower blood cortisol on average, with a compressed dynamic range compared to their control siblings (Fig. 1b–d). Stress-induced spikes in brain cortisol (as appears to occur in the treated fish in response to fasting, Fig. 1b) would be expected to downregulate the hypothalamus–pituitary–adrenal/interrenal axis via glucocorticoid receptor (GR)-mediated negative feedback. Consistent with this, in fasted fish derived from cortisol-treated embryos, blood cortisol levels were similar to those of their untreated siblings that had been exposed to Dexamethasone (Dex) for 8 h and were not further reduced by the Dex treatment (Fig. 1d). Moreover, brain expression of the ACTH-encoding gene pomca, a target of GR-mediated repression, was significantly lower in fasted fish derived from cortisol-treated embryos and is less affected by Dex (Fig. 1e). Altogether these data confirm our earlier finding that whole-body cortisol levels are elevated in adult fish derived from cortisol-treated embryos, and show further that cortisol is aberrantly regulated in those adults, with a compressed dynamic range in blood that correlates with an expanded dynamic range in brain compared to controls. One possible explanation for this is that the cortisol treatment may have persistent effects on the system of cortisol buffering and transport in plasma, for example by affecting the synthesis or metabolism of plasma cortisol-binding proteins that regulate steroid delivery to tissues [19].

Adults derived from cortisol-treated embryos differentially express the glucocorticoid-responsive regulatory genes klf9 and fkbp5 in blood and brain

To determine if developmental exposure to cortisol has persistent effects on gene activity at the level of chromatin we performed ATAC-seq [20, 21] on blood cells isolated from three independent replicates comparing ~ 1-year old adults derived from cortisol-treated embryos and their vehicle-treated (control) siblings. Each ATAC-seq peak was scored according to its size in cortisol-treated fish relative to controls, with higher scores indicating proportionally more sequence reads and hence chromatin with a more open configuration in the cortisol treated group. Of ~ 20,000 scored peaks, that with the 3rd highest score encompassed the promoter of klf9 (Fig. 2a; Table 1), a known GR target gene that functions as a feedforward regulator of GR signaling [22,23,24,25,26]. Several additional genes with high peak scores are also known GR targets in mammals, including fkbp5 (Fig. 2a; Table 1), a clinically important feedback regulator of the GR [23, 27,28,29,30]. Three of the genes in the top 35 (chac1, klf9, and fkbp5) were found in our previous study to be highly upregulated in larvae chronically exposed to cortisol (Table 1). Interestingly, in a HOMER motif enrichment analysis [31] of the 251 peaks with scores  ≥ 100, four of the top ten scoring motifs were consensus binding sequences for krüppel-like factors, including Klf9 (Additional file 2: Figure S1).

Fig. 2
figure 2

The regulatory genes klf9 and fkbp5 have higher on average activity in blood cells of adults derived from cortisol-treated embryos. a ATAC-seq peaks associated with klf9 and fkbp5, from three biological replicates of each treatment (r1-r3). For each replicate sample blood was pooled from 6 individuals mixed sex with equal representation of males and females. b Relative expression of klf9 and fkbp5 in 17 blood samples of adults from different experimental cohorts of cortisol-treated embryos compared to their control siblings sampled at the same time. The averages ± SEM are also shown. For each experimental sample blood was pooled from 6 individuals of each group (control and treated), of mixed sex with equal representation of males and females. c ChIP-qPCR of H3K4me3 levels in the promoter regions of klf9 and fkbp5 from a single sample of blood pooled from 6 individuals (mixed sex, equal representation). The averages ± SD of three replicate qPCR measurements (technical replicates) are shown. d Relative expression of klf9 and fkbp5 in blood cells immediately after being drawn, and then after 6 or 24 h of ex vivo culture in the absence of cortisol. The plots represent the grand means ± the SEM of three biological replicates, each done on blood samples pooled from 6 individuals of mixed sex

Table 1 Top 35 ATAC-seq peaks in blood cells of 1-year old fish derived from cortisol-treated embryos

The ATAC-seq results suggest that adults derived from cortisol-treated embryos persistently upregulate the activities of both klf9 and fkbp5, possibly in response to the persistently elevated cortisol levels. To further test this, quantitative reverse transcription and polymerase chain reaction (qRT-PCR) of both genes was carried out on blood cells of adults derived from cortisol-treated embryos and their vehicle-treated siblings. In 14 out of 17 independent blood cell samples, klf9 and/or fkbp5 expression was elevated in the cortisol-treated fish, both genes being overexpressed ~ twofold on average (Fig. 2b). The fact that fkbp5 and klf9 expression was elevated in most, but not all the samples from treated fish might relate to our observation that blood cortisol levels of the treated fish were lower than those of their untreated siblings after fasting (Fig. 1b–d), suggesting that the three samples that display lower klf9 and fkbp5 expression (8, 16, and 17) may have been recently stressed. Chromatin immunoprecipitation detected higher levels of H3K4 trimethylation in the promoter regions of both genes in blood cells of the treated fish (Fig. 2c), suggesting that the on-average elevated expression in those fish is transcriptional. This is likely a direct effect of elevated blood cortisol, as, when blood cells were cultured ex vivo the expression of both genes decreased over time (Fig. 2d).

In an initial experiment that surveyed the expression of 95 genes in adult brains, fkpb5 and klf9 stood out as the two that were most highly overexpressed in the cortisol-treated group compared to their control siblings (Additional file 3: Figure S2a). However, this result was not reproduced in measurements from additional experimental cohorts (Additional file 3: Figure S2b–f), suggesting that the fish used for the initial survey may have recently experienced a stressor that produced a spike in brain cortisol levels of the treated group (e.g., as observed with fasting; Fig. 1b). One strikingly consistent finding from the gene expression measurements made from both blood cells (Fig. 2b) and brain (Additional file 3: Figure S2) was that expression of fkbp5 and klf9 was highly correlated. This could simply reflect that they are both glucocorticoid-responsive genes, but it could also be indicative of additional regulatory interactions between the two genes. Finally, although most of the gene expression measurements were made from pools of mixed sexes containing equal numbers of matched-sized males and females, in an experiment where the sexes were segregated we found that brain expression levels for both genes were higher in males, and that the differential expression in treated and control groups was greater in males than in females (Additional file 3: Figure S2f), indicating that both genes are also regulated by sex.


These results extend our earlier findings showing that in zebrafish, chronic glucocorticoid exposure during early development leads to persistent dysregulation of the neuroendocrine stress system and associated gene expression in adulthood. The data indicate that the neuroendocrine dysregulation occurs a multiple levels, including cortisol production (higher in treated fish), delivery of cortisol from circulation to tissues (lower blood retention in the treated fish after fasting, with concomitantly higher levels in brain and peripheral tissues), and activity of key genes that regulate GC signaling dynamics in receiving cells (klf9 and fkbp5). These results underscore the relevance of zebrafish as a model organism for investigating the mechanisms underlying developmental programming and neuroendocrine dysregulation associated with ELS and chronic glucocorticoid exposure in humans.


Zebrafish procedures

Wild-type zebrafish of the AB strain maintained in the MDI Biological Laboratory zebrafish facility as previously described [13] were used for all of the reported experiments. Cortisol treatment was carried out as previously described [13]. For the dexamethasone treatments adult zebrafish were removed from the re-circulating system and transferred into 1 L holding tanks containing either 1 µM dexamethasone or 2.5 ppm vehicle (DMSO) in system water. For tissue collection fish were euthanized in tricaine as previously described [13].

Ex vivo blood cell culture

Blood was collected by a nonlethal venous puncture under the operculum. 1µL samples of blood from 6 fish of mixed sex per condition were pooled in heparin (5 I.U./mL) in PBS, washed, and plated at a density of 1 million cells/mL in supplemented L-15 media (l-glutamine, 10% FBS, 25 mM HEPES, 2.2 g/L sodium bicarbonate, Primocin). Cells were cultured at 28° C in 5% CO2. Blood cell viability was assessed by cell count with a hemocytometer. Zebrafish blood cells were determined to be stable in culture for at least 5 days.

Cortisol measurements

Cortisol was measured by ELISA as previously described [13], using dissected tissues or blood samples. In some cases, blood was collected by nonlethal venous puncture under the operculum, allowing repeated measurements (e.g. as in Fig. 1c).


ATAC-seq libraries from blood cells (~ 200,000 cells per sample) from three biological replicates of adults derived from cortisol-treated embryos and their control (vehicle-treated) siblings (10, 11, and 12 months old) were constructed and sequenced and the sequence reads processed as described in the Additional file 4: Methods [20, 21].

RNA isolation and quantitative reverse transcription—polymerase chain reaction qRT-PCR

Total RNA was isolated from blood cells or dissected brain tissue using Trizol and ethanol precipitation. qRT-PCR was performed and the data quantified by the ddCt method as previously described [13].

Chromatin immunoprecipitation (ChIP)

ChIP was performed on formaldehyde-crosslinked chromatin from blood cells, essentially using the procedure described by Lindeman et al. [32] an anti-H3K4me3 antibody (Abcam Cat. No. Ab8580) (detailed ChIP protocol provided on request). Input and immunoprecipitated DNA samples were analyzed by qPCR with Sybr Green Fastmix. The primers were targeted to promoter regions of klf9 (CACATCGACACCGCCTTCATCTG; ATCTTGAACCTCGCCGCTGATTG) and fkbp5 (ACACAACGGGTTACGGGTC; AAATTAAGGGCAGGCCTTGGA). A standard curve was made using serial dilutions of control input DNA and regression analysis was used to calculate percent recovery (immunoprecipitated DNA/input DNA) × 100%.

Statistical analysis

Cohort sizes were determined by power analysis in R with the “pwr” package, using a Cohen’s d effect size calculated from previous experiments of cortisol measurements in adult zebrafish treated as larvae with 1 uM cortisol. This analysis indicated that a minimum n of 5 is needed to detect statistically significant (p < 0.05) differences between treatment and controls. Some experiments had multiple biological replicates, whereas others were done only once with technical replicates, as noted in the Figure legends. Tissue cortisol measurements in Fig. 1b were analyzed with three-way ANOVA to assess the interactions of tissue type, fed/fasted state, and larval cortisol treatment on adult tissue cortisol level. Post-hoc two-way ANOVA was then performed for each tissue type. Grand means in Fig. 2b were assessed with a two-tailed t test.


  • The measurements shown in Fig. 1 were done only once.

  • When interpreting the single time point measurements reported here it is important to consider that cortisol levels are dynamic, with circadian rhythms as well as ultradian fluctuations with an order of magnitude dynamic range.

Availability of data and materials

The ATAC-seq dataset generated during this study is available in the NCBI Gene Expression Omnibus (GEO) repository, under accession number GSE137987.



Early life stress




Glucocorticoid receptor


Assay for transposase accessible chromatin–sequencing


Quantitative reverse transcription and polymerase chain reaction


Enzyme linked immunosorbant assay


  1. Harris A, Seckl J. Glucocorticoids, prenatal stress and the programming of disease. Horm Behav. 2011;59(3):279–89.

    CAS  PubMed  Google Scholar 

  2. Khulan B, Drake AJ. Glucocorticoids as mediators of developmental programming effects. Best Pract Res Clin Endocrinol Metab. 2012;26(5):689–700.

    CAS  PubMed  Google Scholar 

  3. Chen T, Liu HX, Yan HY, Wu DM, Ping J. Developmental origins of inflammatory and immune diseases. Mol Hum Reprod. 2016;22(8):858–65.

    PubMed  PubMed Central  Google Scholar 

  4. Moisiadis VG, Matthews SG. Glucocorticoids and fetal programming part 1: outcomes. Nat Rev Endocrinol. 2014;10(7):391–402.

    CAS  PubMed  Google Scholar 

  5. Moisiadis VG, Matthews SG. Glucocorticoids and fetal programming part 2: mechanisms. Nat Rev Endocrinol. 2014;10(7):403–11.

    CAS  PubMed  Google Scholar 

  6. Oster H, Challet E, Ott V, Arvat E, de Kloet ER, Dijk DJ, Lightman S, Vgontzas A, Van Cauter E. The functional and clinical significance of the 24-h rhythm of circulating glucocorticoids. Endocr Rev. 2016.

    Article  PubMed Central  Google Scholar 

  7. Spencer KA. Developmental stress and social phenotypes: integrating neuroendocrine, behavioural and evolutionary perspectives. Philos Trans R Soc Lond B Biol Sci. 2017;372(1727):20160242.

    PubMed  PubMed Central  Google Scholar 

  8. Danese A, McEwen BS. Adverse childhood experiences, allostasis, allostatic load, and age-related disease. Physiol Behav. 2012;106(1):29–39.

    CAS  PubMed  Google Scholar 

  9. Cachat J, Stewart A, Grossman L, Gaikwad S, Kadri F, Chung KM, Wu N, Wong K, Roy S, Suciu C, et al. Measuring behavioral and endocrine responses to novelty stress in adult zebrafish. Nat Protoc. 2010;5(11):1786–99.

    CAS  PubMed  Google Scholar 

  10. Krug RG 2nd, Poshusta TL, Skuster KJ, Berg MR, Gardner SL, Clark KJ. A transgenic zebrafish model for monitoring glucocorticoid receptor activity. Genes Brain Behav. 2014;13(5):478–87.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Griffiths BB, Schoonheim PJ, Ziv L, Voelker L, Baier H, Gahtan E. A zebrafish model of glucocorticoid resistance shows serotonergic modulation of the stress response. Front Behav Neurosci. 2012;6:68.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Nesan D, Vijayan MM. The transcriptomics of glucocorticoid receptor signaling in developing zebrafish. PLoS ONE. 2013;8(11):e80726.

    PubMed  PubMed Central  Google Scholar 

  13. Hartig EI, Zhu S, King BL, Coffman JA. Cortisol-treated zebrafish embryos develop into pro-inflammatory adults with aberrant immune gene regulation. Biol Open. 2016;5(8):1134–41.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Weger BD, Weger M, Nusser M, Brenner-Weiss G, Dickmeis T. A chemical screening system for glucocorticoid stress hormone signaling in an intact vertebrate. ACS Chem Biol. 2012;7(7):1178–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Benato F, Colletti E, Skobo T, Moro E, Colombo L, Argenton F, Dalla Valle L. A living biosensor model to dynamically trace glucocorticoid transcriptional activity during development and adult life in zebrafish. Mol Cell Endocrinol. 2014;392(1–2):60–72.

    CAS  PubMed  Google Scholar 

  16. De Marco RJ, Groneberg AH, Yeh CM, Castillo Ramirez LA, Ryu S. Optogenetic elevation of endogenous glucocorticoid level in larval zebrafish. Front Neural Circuits. 2013;7:82.

    PubMed  PubMed Central  Google Scholar 

  17. Castillo-Ramirez LA, Ryu S, De Marco RJ. Active behaviour during early development shapes glucocorticoid reactivity. Sci Rep. 2019;9(1):12796.

    PubMed  PubMed Central  Google Scholar 

  18. Gutierrez-Triana JA, Herget U, Castillo-Ramirez LA, Lutz M, Yeh CM, De Marco RJ, Ryu S. Manipulation of interrenal cell function in developing zebrafish using genetically targeted ablation and an optogenetic tool. Endocrinology. 2015;156(9):3394–401.

    CAS  PubMed  Google Scholar 

  19. Hammond GL. Plasma steroid-binding proteins: primary gatekeepers of steroid hormone action. J Endocrinol. 2016;230(1):R13–25.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Buenrostro JD, Wu B, Chang HY, Greenleaf WJ. ATAC-seq: a method for assaying chromatin accessibility genome-wide. Curr Protoc Mol Biol. 2015;109:21–9.

    PubMed  PubMed Central  Google Scholar 

  21. Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10(12):1213–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Cui A, Fan H, Zhang Y, Zhang Y, Niu D, Liu S, Liu Q, Ma W, Shen Z, Shen L, et al. Dexamethasone-induced Kruppel-like factor 9 expression promotes hepatic gluconeogenesis and hyperglycemia. J Clin Invest. 2019;130:2266–78.

    Google Scholar 

  23. Juszczak GR, Stankiewicz AM. Glucocorticoids, genes and brain function. Prog Neuropsychopharmacol Biol Psychiatry. 2018;82:136–68.

    CAS  PubMed  Google Scholar 

  24. Shewade LH, Schneider KA, Brown AC, Buchholz DR. In-vivo regulation of Kruppel-like factor 9 by corticosteroids and their receptors across tissues in tadpoles of Xenopus tropicalis. Gen Comp Endocrinol. 2017;248:79–86.

    CAS  PubMed  Google Scholar 

  25. Chinenov Y, Coppo M, Gupte R, Sacta MA, Rogatsky I. Glucocorticoid receptor coordinates transcription factor-dominated regulatory network in macrophages. BMC Genomics. 2014;15:656.

    PubMed  PubMed Central  Google Scholar 

  26. Sporl F, Korge S, Jurchott K, Wunderskirchner M, Schellenberg K, Heins S, Specht A, Stoll C, Klemz R, Maier B, et al. Kruppel-like factor 9 is a circadian transcription factor in human epidermis that controls proliferation of keratinocytes. Proc Natl Acad Sci U S A. 2012;109(27):10903–8.

    PubMed  PubMed Central  Google Scholar 

  27. Binder EB. The role of FKBP5, a co-chaperone of the glucocorticoid receptor in the pathogenesis and therapy of affective and anxiety disorders. Psychoneuroendocrinology. 2009;34(Suppl 1):S186–95.

    CAS  PubMed  Google Scholar 

  28. Klengel T, Mehta D, Anacker C, Rex-Haffner M, Pruessner JC, Pariante CM, Pace TW, Mercer KB, Mayberg HS, Bradley B, et al. Allele-specific FKBP5 DNA demethylation mediates gene-childhood trauma interactions. Nat Neurosci. 2012;16(1):33–41.

    PubMed  PubMed Central  Google Scholar 

  29. Zannas AS, Wiechmann T, Gassen NC, Binder EB. Gene-stress-epigenetic regulation of FKBP5: clinical and translational implications. Neuropsychopharmacology. 2016;41(1):261–74.

    CAS  PubMed  Google Scholar 

  30. Touma C, Gassen NC, Herrmann L, Cheung-Flynn J, Bull DR, Ionescu IA, Heinzmann JM, Knapman A, Siebertz A, Depping AM, et al. FK506 binding protein 5 shapes stress responsiveness: modulation of neuroendocrine reactivity and coping behavior. Biol Psychiatry. 2011;70(10):928–36.

    CAS  PubMed  Google Scholar 

  31. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38(4):576–89.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Lindeman LC, Vogt-Kielland LT, Alestrom P, Collas P. Fish’n ChIPs: chromatin immunoprecipitation in the zebrafish embryo. Methods Mol Biol. 2009;567:75–86.

    CAS  PubMed  Google Scholar 

Download references


Not applicable.


Research reported in this publication was supported by Institutional Development Awards (IDeA) from the National Institute of General Medical Sciences of the National Institutes of Health under grant numbers P20-GM104318 and P20-GM103423. The funding sources played no role in the design of the study, in collection, analysis, and interpretation of data, or in writing the manuscript.

Author information

Authors and Affiliations



EIH designed the fed versus fasting cortisol experiment, carried out all the zebrafish manipulations, tissue dissections, blood draws, cortisol measurements, qRT-PCR and ChIP analyses from blood, and statistical analyses.SZ constructed the ATAC-seq libraries, developed the ChIP protocol that was used in the lab and carried out the qRT-PCR measurements from brain tissue. BLK carried out the ATAC-seq data analyses. JAC conceived and designed the overall project, drafted the figures, and wrote the manuscript. All authors read and approved the final manuscript prior to submission.

Corresponding author

Correspondence to James A. Coffman.

Ethics declarations

Ethics approval and consent to participate

All animal procedures were approved by the Institutional Animal Use and Care Committee of the MDI Biological Laboratory. This research did not involve human subjects.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1:

ANOVA of tissue cortisol levels shown in Fig. 1b.

Additional file 2: Figure S1.

Top ten transcription factor binding motifs identified by HOMER motif enrichment analysis of sequences from the 251 peaks that scored ≥ 100.

Additional file 3: Figure S2.

Survey of gene expression by 96-well qPCR array and relative expression levels of klf9 and fkbp5 in brain tissue dissected from different experimental cohorts. Each measurement was made from pooled brain RNA of 6 individuals of mixed sex (equal representation), except where indicated (panel f). The qPCR measurement shown in panel (b) is from the same RNA samples that were used to generate the data depicted in panel (a).

Additional file 4: Methods.

ATAC-seq library construction protocol, sequencing and read processing.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hartig, E.I., Zhu, S., King, B.L. et al. Chronic cortisol exposure in early development leads to neuroendocrine dysregulation in adulthood. BMC Res Notes 13, 366 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: