Integrative genomic analysis reveals a conserved role for prolactin signalling in the regulation of adrenal function

Dear Editor, We report here that the pituitary hormone prolactin (PRL) has an important, conserved role in regulating adrenal gland function. The adrenal plays a pivotal role in endocrine homeostasis and stress response,which vary across lifespan andhave different features in the two sexes.[1] Furthermore, most adrenal disorders have a higher prevalence in women.[2] To get deeper insight into the mechanisms of ageand sexdependent adrenal function, we performed an integrative analysis of the mouse adrenal transcriptome by RNA-seq and active enhancer (as defined by H3K27ac ChIP-seq) usage at different ages (from E18.5 to P12 weeks) in both sexes. Multivariate analysis showed that age, but not sex, had a significant effect on adrenal global gene expression profiles (Figure and Supporting information S1; Tables S1, S2). Expression of genes related to lipid metabolism and immunity was progressively enriched with age (Figure 1C and Supporting information Table S3), with the proportion of macrophages increasing and neutrophils/mast cells decreasing at older ages (Figure 1D). The expression of different classes of secreted proteins, ion channels and enzymes showed a lower level of sex-dependent enrichment (Figure 1E and Supporting information Table S4). Age-and sex-dependent long non-coding RNA differentially expressed genes (DEG) followed the same pattern as coding DEG, being mostly specific for each age and sex (Supporting information Figure S2 and Table S5). Remarkably, P12 weeks adrenals of both sexes showed increased percentages of novel gene transcripts compared to previous ages (Supporting information Figure S2E). Sexually dimorphic expression pervades all cell populations of the mouse adrenal gland (Supporting information Figure S3 and Table S6). We identified a few thousand enhancers active in the adrenal gland at each temporal stage and in each sex (Supporting information Table S7). A subset among those belong to the super-enhancer class (SEC), genomic regions