Sex and age differences in the expression of liver microRNAs during the life span of F344 rats
Biology of Sex Differences volume 8, Article number: 6 (2017)
Physiological factors such as age and sex have been shown to be risk factors for adverse effects in the liver, including liver diseases and drug-induced liver injury. Previously, we have reported age- and sex-related significant differences in hepatic basal gene expression in rats during the life span that may be related to susceptibility to such adverse effects. However, the underlying mechanisms of the gene expression changes were not fully understood. In recent years, increasing evidence for epigenetic mechanisms of gene regulation has fueled interest in the role of microRNAs (miRNAs) in toxicogenomics and biomarker discovery. We therefore proposed that significant age and sex differences exist in baseline liver miRNA expression, and that comprehensive profiling of miRNAs will provide insights into the epigenetic regulation of gene expression in rat liver.
To address this, liver tissues from male and female F344 rats were examined at 2, 5, 6, 8, 15, 21, 52, 78, and 104 weeks of age for the expression of 677 unique miRNAs. Following data processing, predictive pathway analysis was performed on selected miRNAs that exhibited prominent age and/or sex differences in expression.
Of the 314 miRNAs found to be expressed, 214 were differentially expressed; 65 and 212 miRNAs showed significant (false discovery rate (FDR) <5% and ≥1.5-fold change) sex- and age-related differences in expression, respectively. Thirty-eight miRNAs showed 2-week-specific expression, of which 31 miRNAs were found to be encoded within the Dlk1-Dio3 cluster located on chromosome 6. This cluster has been associated with tissue proliferation and differentiation, and liver energy homeostasis in postnatal development. Predictive pathway analysis linked sex-biased miRNA expression with sexually dimorphic molecular functions and toxicological functions that may reflect sex differences in hepatic physiology and disease. The expression of miRNAs (miR-18a, miR-99a, and miR-203, miR-451) was also found to associate with specific sexually dimorphic hepatic histopathology. The expression of miRNAs involved in regulating cell death, cell proliferation, and cell cycle was found to change as the rats matured from adult to old age.
Overall, significant age- and sex-related differences in liver miRNA expression were identified and linked to histopathological findings and predicted functional pathways that may underlie susceptibilities to liver toxicity and disease.
The liver plays a crucial role in metabolizing food and drugs, and thus is of great toxicological importance. Physiological factors such as age and sex have been suggested as risk factors for liver diseases and drug-induced liver injury (DILI). For example, anti-epileptic drugs, such as carbamazepine and valproate, have been associated with age-dependent susceptibilities to adverse effects [1, 2]. Elderly patients who take carbamazepine showed an increased risk for adverse liver reactions [3, 4], whereas younger patients taking valproate show increased risk for hepatotoxicity . In addition, other drugs, including erythromycin, halothane, isoniazid, nitrofurantoin, and flucloxacillin, have recently been associated with age-dependent increases in adverse liver effects .
Moreover, sex-related differences in susceptibility to DILI have been observed. For example, female patients show increased incidence of anesthesia-related DILI compared to males . This sex-bias for anesthesia-related DILI was also demonstrated in a mouse model, which was further linked to an estrogen-dependent pro-inflammatory response involving IL-6 suppression of regulatory T-cells . Clinical and epidemiological studies suggest that females show more sensitivity to adverse events after taking acetaminophen [9, 10], and evidence supporting a general female-bias in DILI has recently been reviewed . However, a precise mechanism for the sex difference is not clearly understood. Several other drugs exhibit sex-biased sensitivities to DILI, but it is not clear whether males or females show more sensitivity to DILI in general [12–14], suggesting the sex bias may be drug-dependent. Establishment of the underlying molecular basis for such age- and sex-related susceptibilities to toxicities or diseases will enable safer use of such drugs.
We have previously reported marked age- and sex-related differences in basal hepatic gene expression in rats . Differential expression was observed for over 3500 genes in a variety of different pathways. It is recognized that messenger RNA (mRNA) expression is under the control of various transcriptional and post-transcriptional mechanisms. miRNAs are short, non-coding RNAs of 19–23 nucleotides in length that regulate gene expression, typically by translational repression of target mRNAs [16–19]. A single miRNA may target a variety of mRNA transcripts and single transcripts may be targeted by multiple miRNAs; thus, their role in mRNA regulatory networks is complex. miRNAs have been shown to influence diverse biological pathways and functions, with a prominent role in various cancers [20, 21]. Furthermore, miRNA biomarkers have shown promising utility in distinguishing cancer types , detecting various tissue injuries [23, 24], and making predictions regarding disease outcomes . Altered miRNA expression has been reported in response to liver diseases and injuries, including hepatitis C infection , non-alcoholic fatty liver disease (NAFLD) , and DILI [28, 29]. Recent studies suggest that miRNAs may play an important role as biomarkers and tools in toxicogenomics [30–32].
A few miRNAs have been shown to exhibit age- [33–35] and sex-dependent  patterns of expression. miR-34a was among the first miRNAs to be characterized for its role in aging and longevity in C. elegans [37, 38]. miR-34a is a downstream target of p53 and its expression correlates positively with age . Similarly, increased liver expression of miR-93, miR-214, and miR-669c has been correlated with age. These miRNAs target glutathione-S-transferases, which are crucial in protecting against DILI-related oxidative stress . Li et al.  have also shown age-dependent increases in the expression of miR-34a and miR-93, and that these miRNAs target Mgst1, Sirt1, and Nrf2, three genes that encode proteins important in the defense against oxidative stress, a common feature in DILI. These studies demonstrate the influence of miRNAs on pathways related to DILI in an age-dependent manner and thus may help develop a fuller understanding of age-related susceptibilities to drug toxicities.
The roles of miRNAs in sexually dimorphic physiology and diseases have recently been reviewed . Sex differences in expression can be regulated through both genetic factors (e.g., X-chromosome encoded miRNAs) and hormonal factors (e.g., sex steroids) . The primary direct regulator of sex differences in gene expression in young adult and adult liver is growth hormone . In a carbon tetrachloride-induced liver fibrosis mouse model, females exhibited protection against fibrosis compared to males, which was associated with miR-29 family members. The miR-29 family members are induced by estrogen and reduce fibrosis by inhibiting expression of collagens . Additionally, male-biased expression of miR-216a via the androgen pathway has been reported in early hepatocarcinogenesis . miR-216a has been shown to target and decrease levels of Tslc1, whose expression is decreased in a number of tumors , thus implicating the sex-biased expression of miR-216a in carcinogenesis. Taken together, there is evidence to suggest potential roles for miRNAs in age- and sex-related susceptibilities to liver diseases and toxicities.
Currently, no comprehensive reports characterizing sex and age differences in miRNA expression in rat liver tissue are available. In the present study, whole genome expression profiles of miRNA in the liver of male and female F344 rats during the life span were measured. Results showed substantial sex and age differences in basal liver miRNA expression levels that may have a role in regulating expression of genes associated with hepatic proliferation and differentiation, fibrosis, and steatosis. These data will facilitate a better understanding of epigenetic regulation of sex- and age-related differences in hepatic gene expression and will help elucidate the molecular basis underlying differential susceptibilities to toxicities at different stages of life and between the sexes.
All tissue samples came from a previously described study  and were maintained at −80 °C prior to utilization in the present project. Briefly, F344 rats were obtained from the National Center for Toxicological Research (NCTR) rat breeding colony. Male and female animals were euthanized at 2, 5, 6, 8, 15, 21, 52, 78, and 104 weeks of age as previously described. Female estrus cycles were not evaluated or synchronized. Animal care and procedures were approved by the NCTR Institutional Animal Care and Use Committee. Animals were housed in polycarbonate cages with hardwood chip bedding, two per cage, at 23 °C with a relative humidity of approximately 50%. Rats were euthanized at the same circadian time of day (between 0900 and 1100) by carbon dioxide asphyxiation. Liver tissue sections from rats aged 52, 78, and 104 weeks of age were collected for histological evaluation (Fig. 7). Tissue sections were placed in 10% neutral buffered formalin followed by standard paraffin embedding, sectioning, and hematoxylin and eosin staining. Liver sections were evaluated for lesions, including basophilic foci, hemorrhage, hyperplasia of the bile duct, hepatocyte hypertrophy, lymphocyte infiltration, chronic active inflammation, necrosis, portal vein thrombosis, and hepatocyte cytoplasmic vacuolization by staff board-certified pathologists. Each lesion was graded for severity using the following scale: 1 = minimal, 2 = mild, 3 = moderate, and 4 = marked. Liver sections were not examined in the younger rats because lesions are only generally found in older rats . Contemporaneous National Toxicology Program studies, using F344 rats from the same source and living in the same environmental conditions at NCTR (including feed) found no basophilic foci, hyperplasia of the bile duct, or hepatocyte cytoplasmic vacuolization in rats of 6–7 weeks of age (0/12 males and 0/12 females) and rats of 17–18 weeks of age (0/24 males and 0/24 females) .
Total RNA including the small RNA fraction were isolated from ~30 mg frozen liver tissue that was ground to a powder in liquid nitrogen and homogenized as previously described  following the manufacturer’s protocol (Qiagen miRNeasy Mini Kit, Qiagen Inc. Valencia, CA, USA). RNA samples were evaluated spectrophotometrically on a Nanodrop-1000 (ThermoFisher Scientific, Waltham, MA, USA) by measuring the optical density at 260 nm. RNA quality of each sample was subsequently evaluated on RNA 6000 LabChips using the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). RNA samples with RNA integrity numbers (RINs) greater than 8.0 were used in miRNA microarray experiments; the average RIN was 8.6.
Comprehensive miRNA expression profiling was performed using single color (Cy3) Agilent Rat 8x15k miRNA (version 16) microarrays and labeling reagents according to the manufacturer’s protocol as previously described . A total of 88 arrays were used (n = 5 rat liver RNA samples for all groups, except 5- and 104-week-old males, where n = 4 rat liver RNA samples). 100 ng of total RNA, including the small RNA fraction, was used as input. Agilent spike-in controls were used to evaluate both labeling and hybridization efficiency, and all were within acceptable ranges. BioRad Micro Bio-Spin 6 columns were used to purify labeled probe (BioRad Life Science Research, Hercules, CA) according to manufacturer’s recommendation. Samples were processed on four separate days (four batches). To reduce possible batch effects, samples from each age-sex group were randomized among the four batches so that samples from each group were processed in three or four of the batches, with the exception of the 8 week old female group which was processed in two batches. Hybridized miRNA microarrays were then scanned on an Agilent DNA Microarray Scanner (Agilent Technologies) according to the manufacturer’s instructions.
Microarray data processing
The procedure for processing the miRNA microarray data in the current liver data set is identical to the procedure used in the previously published kidney miRNA microarray data set . This procedure and details pertinent to the liver data are described below for reference. Microarray images (TIF) were analyzed using Feature Extraction (FE) software (V10.7) from Agilent for preliminary image quality control evaluations which utilized the labeling and hybridization spike-in probes and features. FE quality control metrics for each slide were within the ranges designated as “good” or “excellent”. Raw median intensities were calculated for each replicate probe (n = 20 for each miRNA) in FE and used to generate a single value (Total Gene Signal, TGS) for each miRNA. Individual data files were generated in FE for each array (88 total files). The 88 individual files were combined into a single file containing TGS data for each array. TGS values of 0.1 represent a placeholder indicating miRNA signal below the level of detection for the assay, equivalent to “not expressed”. Values of 0.1 were subsequently replaced by the global minimum intensity value (1.84) across all arrays. A single file containing processed TGS values was used as input for further processing in SAS 9.3 (SAS Institute Inc., Cary, NC). SAS was used to remove control spots and target spot redundancy/replication. The Agilent Rat 8x15k miRNA microarray (v16) contains 15,744 array spots; 2204 control spots were removed, leaving 13,540 spots of target miRNA probes (677 unique miRNAs x 20 replicates per miRNA). Redundant replicates containing identical TGS values were removed for each miRNA, resulting in data for 677 unique miRNAs. The 363 non-expressed miRNAs (miRNAs with original unprocessed TGS values of 0.1 across all arrays) were discarded, leaving 314 unique miRNAs expressed in at least one rat liver sample. The TGS values (raw intensities) of these 314 “expressed” miRNAs were then transformed to log2 values in SAS followed by 75th percentile normalization.
To prevent generating false fold differences from placeholders, transformed 1.84 values were replaced by 0.012 (lowest transformed value among all samples with measured values).
Microarray statistical analysis
Statistical analysis was performed in SAS using a generalized linear model procedure (proc glm). Model (RandommodelANOVA) and estimate statements were used to perform ANOVA and pairwise comparisons using least square method to measure the significance of group differences using log2 intensities as the dependent variable and age or sex as independent variables. Significant age differences in miRNA expression were calculated separately for males and females. Significant sex differences in miRNA expression were calculated at each age. The ANOVA p value and p value from pairwise tests were adjusted to false discovery rate (FDR). The magnitude of group differences was calculated as the difference of average log2 intensities; the anti-log of the resulting value was used as the fold change value between groups. When group absolute fold changes were greater than or equal to 1.5 and FDR for ANOVA as well as pairwise tests were significant (p < 0.05), miRNAs were determined to be differentially expressed by age. Differential expression by sex was determined using absolute fold changes greater than or equal to 1.5 between sexes of the same age and a pairwise t test with FDR < 5%. Normalized intensity values (normalized TGS) were subsequently uploaded into ArrayTrack , the Food and Drug Administration’s relational database for genomic data storage, processing, analysis, and visualization that was created at NCTR . ArrayTrack is available for free .
Principal component analysis (PCA) and k-means cluster analysis were performed in ArrayTrack on log2-transformed, normalized intensity values.
Raw expression data, as well as processed data, were deposited in Gene Expression Omnibus (Accession ID: GSE65169).
Taqman qPCR experiments
Quantitative PCR (qPCR) experiments were carried out using Taqman miRNA assays. 5 ng of total RNA from each sample (n = 5 per group) was used in a 15-μl reverse transcription reaction using the Taqman miRNA Reverse Transcription Kit (Life Technologies, Grand Island, NY) along with miRNA-specific primers according to manufacturer’s protocol. The reverse transcription product (1.33 μl) was subsequently used in a 20-μl Taqman MicroRNA Assay reaction with predesigned miRNA-specific probes. The qPCR reactions were performed in MicroAmp Optical 384-well reaction plates on the ABI 7900HT real-time PCR detection system (ThermoFisher). Taqman assays were performed for miR-154* (Assay ID: 000478), miR-29b (Assay ID: 000413), miR-96 (Assay ID: 000186), and miR-34a (Assay ID: 000425) due to their age- and/or sex-related expression differences. U6 snRNA (Assay ID: 001973) was used as endogenous control in the Taqman assay and ∆Ct was further normalized to the global average of each miRNA using the ∆∆Ct method. The average log2 relative expression per group of qPCR and microarray data are displayed for miR-29b, miR-34a, miR-96, and miR-154* such that the global average expression for qPCR (n = 5) and microarray data (n = 4 or 5) is equal to zero. Sex differences in the qPCR data were tested for statistical significance using a simple t test (p < 0.05). Taqman qPCR and microarray expression profiles for each miRNA were compared using a Pearson correlation coefficient (r) in Microsoft Excel.
miRNA pathway analysis
QIAGEN’s Ingenuity Pathway Analysis software (IPA; QIAGEN, Redwood City, CA) was used to associate differentially expressed miRNAs with their target mRNAs, and pathway analysis. Briefly, lists of differentially expressed miRNAs (using the Agilent probe identifiers) were imported into IPA, and the miRNA target filter function in IPA was used to identify their corresponding mRNA targets. The analysis was restricted to experimentally observed interactions retrieved from TarBase , miRecords , and curated peer-reviewed literature. The output of this analysis was exported as an Excel file, duplicate mRNA targets removed, and the file of mRNA targets uploaded to IPA. These lists of target mRNAs were then examined for over-representation in biological pathways curated in IPA using the Core Analysis option. All settings were default settings except the species selected was rat. Molecular and cellular functions, as well as physiological system development and function under the diseases and functions tab were identified. Tox functions, including only hepatotoxicity, were examined as well. Pathways and functions were considered significant if there were at least three mRNAs in the pathway/function, and Fisher’s exact test p value was less than 10−3. The analysis was limited to rat-specific data except for analysis of the 2-week-specific expression class where the more abundant human data was included.
Hepatic miRNA profiles change during the rat life span
The microarrays allowed examination of the expression levels of 677 unique miRNAs. More than half (363 miRNAs) showed no detectable expression at any age or sex and were removed from further analysis. The remaining 314 miRNAs showed a detectable expression level in at least one sample. An ANOVA with an FDR <5% was used to identify age-dependent differences in the expression of the 314 expressed miRNAs in males and females, independently. Subsequently, post hoc multiple comparisons tests (p < 0.05) with a fold change ≥1.5-fold between age groups, identified 176 and 171 miRNAs that showed differential expression by age in females and males, respectively, for a total of 212 age-affected miRNAs in either sex (Additional file 1). miRNAs showing significant sex differences were identified by comparing males and females at each age using a pairwise t test with FDR < 5% and fold change ≥1.5 between the sexes. The application of these criteria identified 65 unique miRNAs that showed sexually dimorphic expression in at least one age. In combination, a total of 214 miRNAs (Additional file 2) were differentially expressed by age and/or sex with 63 miRNAs showing both age and sex differences. These 214 miRNAs were used in further analyses. The expression of more miRNAs was impacted by age than sex, and most sex-related differentially expressed miRNAs (DEMs; 97%; 63/65) also showed an age effect.
PCA, utilizing the expression profiles of these 214 DEMs, was used to evaluate the relationships among the 88 liver samples (Fig. 1). PC1 accounted for a large proportion of the variability (50.7%). Individual animals generally divided into two groups consisting of 2-week-old animals (Fig. 1; blue) and all remaining age groups, suggesting a distinct miRNA expression pattern at 2 weeks of age. The 5 to 104 week age group subdivided slightly between 5 and 21-week-old animals and 52 to 104-week-old animals. Within the 2-week-old group of animals, a slight but consistent division between males (on the left) and females (on the right) was observed. Aside from this slight sexually dimorphic pattern at 2 weeks of age, the remaining age groups showed little to no evidence of sex differences in the PCA plot.
After establishing the global relationships between groups and replicates, k-means cluster analysis was used to visualize patterns among the individual miRNA profiles using the 214 DEMs (Fig. 2; cluster members are shown in Additional file 2). An initial number of k = 10 clusters was used as an estimate of the number of large-scale expression profiles present within the miRNA data. Three trends in the k-means cluster analysis were observed. First, miRNAs generally clustered according to normalized expression level, with highly expressed miRNAs grouping with other highly expressed miRNAs (e.g., clusters 8–10), and low/medium expressed miRNAs with others showing low/medium expression (e.g., clusters 1, 4, 5, and 6). Secondly, miRNAs generally maintained stable expression patterns across age and sex (e.g., clusters 4–10). Thirdly, clusters 2 and 3 diverged from this stable expression profile with a unique 2-week-specific expression pattern in both sexes, followed by little to no expression at subsequent ages. There were a total of 38 miRNAs in clusters 2 and 3 which exhibited a 2-week-specific expression pattern. miRNAs in cluster 3 exhibited the greatest difference in expression between 2 weeks of age and all subsequent ages. Within cluster 3, male and female animals, on average, showed a 28- and 53-fold difference, respectively, in expression between 2 weeks of age and all subsequent ages. Thus, k-means cluster analysis helped identify specific miRNAs that account for the prominent age differences observed in the PCA.
Sex differences in expression of miRNAs
Sex different expression of miRNAs was evaluated at each age (t test with FDR < 5% and fold change ≥ 1.5) and shown in Fig. 3a (females) and Fig. 3b (males), along with the expression fold-differences. A total of 32 miRNAs showed significant female-biased expression at any age while 33 miRNAs exhibited male-biased expression. There was no overlap between the 32 female-biased miRNAs and the 33 male-biased miRNAs. While most of the sex-biased miRNAs occurred at just a single age (54 of the 65 sex-biased miRNAs), 11 (5 female and 6 male) were present at more than one age. Female-biased expression of miR-29b and miR-152 occurred from 15 to 78 weeks of age while male-biased expression of miR-125b-5p and miR-99a occurred during this same age span. A single miRNA (miR-183) showed sex-biased expression at 104 weeks of age and was expressed at levels 17.2 times higher in males than females. Twenty-nine miRNAs were expressed in a sex-biased manner in 2-week-old animals (12 female-biased, 17 male-biased), with differences of 1.5- to 4.8-fold (average = 2.0-fold difference). There were fewer miRNAs showing sex-biased expression at subsequent ages: 4 miRNAs at 8 weeks (3 female, 1 male), 15 miRNAs at 15 weeks (8 female, 7 male), 20 miRNAs at 21 weeks (11 female and 9 male), 9 miRNAs at 52 weeks (2 female and 7 male), and 9 at 78 weeks (5 female and 4 male). Animals between 5 and 6 weeks of age showed no significant sex-biased miRNA expression. Thus, most sex-biased miRNA expression occurred at a single age and occurred predominantly at 2 (29 miRNAs), 15 (15 miRNAs), and 21 (20 miRNAs) weeks of age. Persistent sex-biased expression of 4 miRNAs (miRNAs 29b, 152, 125b-5p, and 99a) occurred in adults at 15–78 weeks of age, with expression fold-differences between the sexes of 1.6–2.9. The sex-biased expression of 9 miRNAs in both sexes during the life span is shown in Fig. 4, with panels a through d showing female-biased miRNA expression, and panels e through i showing male-biased expression of miRNAs.
Age differences in expression of miRNAs
Examination of the PCA (Fig. 1) identified a large difference in liver miRNA expression profiles between 2 weeks of age and the other ages along principal component (PC) 1. In addition, there was a smaller, and incomplete, separation of samples along PC2, with 5–21 weeks of age samples generally grouping away from the more tightly clustered samples from the older ages (52–104 weeks of age). Also, as described below, hepatic histopathology was seen in the livers of the older rats (52–104 weeks of age) and not in the livers of the 5–21 week old rats. For these reasons, the 9 age groups were divided into 3 age classes, young (2 weeks of age), adult (5 to 21 weeks of age) and old (52 to 104 weeks of age), for identification of hepatic age-related DEMs. Age-related DEMs were identified independently in females and males using a t test with an FDR of <5% and requiring a fold change difference of at least 1.5. When comparing young to adult, 159 male DEMs and 131 female DEMs (Fig. 5a) were identified as age-related DEMs. Among these, 115 were in common, suggesting similar miRNA expression patterns in males and females. Indeed, 88% of the female DEMs were also differentially expressed in males. There were also large numbers of DEMs when comparing young to old rats, 161 in males and 148 in females, with 125 (>77%) in common (Fig. 5b). This indicates that there are large differences in the liver miRNA profiles between young rats and adult and old rats, and much smaller differences between adult and old rats. Comparing the miRNA liver expression profiles between adult and old rats identified only 17 DEMs in males and 14 in females, of which 6 were in common (Fig. 5c). The miRNAs in each comparison are shown in Additional file 3.
As noted above, k-means cluster analysis identified 38 miRNAs (clusters 2 and 3 in Fig. 2) that showed a 2-week-specific expression pattern, i.e., there was high expression in the livers of 2-week-old rats but low/undetectable expression at later ages (expression values are shown in Additional file 4). The expression of one of these miRNAs (miR-376b-3p) through the life span is shown in Fig. 4. In males, the expression of miR-376b-3p at 2 weeks of age was >320-fold higher than the average expression at all other ages, while in females it was >150-fold higher. In males, 35 of the 38 miRNAs were expressed at >10-fold greater levels in 2-week-old animals than at all subsequent ages. In females, 28 of the 38 miRNAs were expressed at >10-fold greater levels in 2-week-old animals than at all subsequent ages. Furthermore, male expression was higher than female expression for all except one (miR-675*) of the two-week-specific miRNAs. On average, the expression of these 38 hepatic miRNAs was 1.9-fold higher in males than in females at 2 weeks of age. The genes for 31 of the 38 miRNAs are located on rat chromosome 6 (see Additional file 2 and Additional file 4).
Verification of microarray data by qPCR
Four DEMs that showed significant sex and/or age differences by microarray analysis (miR-29b, miR-96, miR-154*, and miR-34a) were selected for verification by Taqman qPCR assays. Pearson correlation coefficients for the measurements by the two methods for the four miRNAs ranged from 0.835 to 0.966 (shown in Fig. 6), suggesting good agreement between the two methods of expression analysis. Measurement of miR-29b expression by microarray showed female-biased expression at 15, 21, 52, and 78 weeks of age. Significant (p < 0.05) sex difference was confirmed by qPCR for miR-29b only at 52 weeks of age (female > male, Fig. 6a). Microarray data for miR-96 showed male-biased expression at 21, 52, and 78 weeks of age. Taqman qPCR results confirmed this male-biased miR-96 expression at 78 weeks of age and additionally at 104 weeks of age (Fig. 6b). By microarray measurement, miR-154* showed 156- and 383-fold increased expression at 2 weeks of age compared to the average expression at all subsequent ages in females and males, respectively. This large 2-week-specific expression difference was verified in the qPCR data, which showed 2-week expression differences of 1212-fold and 4169-fold in females and males (Fig. 6c), respectively. The larger fold changes determined by qPCR may be due to the larger dynamic range of the qPCR assay. In addition, the significant male-biased expression at 2 weeks of age was confirmed by qPCR. Significant age-related increases in miR-34a expression (7.1- and 13.2-fold maximum age difference in females and males, respectively; microarray data) were also confirmed by qPCR, which showed significant (p < 0.05) 7.2- and 8.7-fold differences in expression across the life span in females and males, respectively (Fig. 6d). Thus, there was good correlation between the expression levels determined by microarray and qPCR, although expression differences of <2-fold were not always confirmed.
Potential pathways and functions of DEMs
Potential mRNA targets of the DEMs were identified using Ingenuity Pathway Analysis software, restricting the analysis to experimentally validated interactions retrieved from TarBase, miRecords, and curated from the peer-reviewed literature. The predicted mRNA targets were then examined for over-representation in biological pathways curated in Ingenuity Pathway Analysis.
2-week-old-specific hepatic miRNA expression
38 miRNAs were identified in the k-means cluster analysis as being expressed in the livers of pre-weaning 2-week-old rats and expressed at reduced levels or not detectable levels at all subsequent ages (Fig. 2, clusters 2 and 3; Additional file 4). Table 1 shows the 38 miRNAs in this 2-week-specific class along with the mRNA targets and associated pathways/functions where the mRNAs were significantly enriched. A relatively small number of mRNA targets (Capn8, Cntn4, Onecut1, Tagln, and Vim) was identified using rat-specific data from IPA and no significant pathways were identified. Human-specific miRNA-mRNA association data were, therefore, included in the analysis and this expanded the mRNA targets by 5, giving a total of 10 mRNA targets for the 38 miRNAs. Only 6 of the miRNAs had associated mRNAs, and these mRNAs were significantly represented in the molecular functions: cellular development, cellular growth and maintenance, cellular growth and proliferation, cell morphology, and cell-to-cell signaling. Notably, 3 of the mRNAs (BCL6, PRDM1, and XBP1), all targets of miR-127, were associated with the hematological system.
Sex different miRNA expression in the liver of adult and old rats
As shown in Fig. 3, 30 hepatic miRNAs were found to be differentially expressed between the sexes in adults (5–21 weeks of age), and 13 were found to be differentially expressed in old rats (52–104 weeks of age). Table 1 shows the sex different miRNAs in these two age groups, as well as the target mRNAs and the significantly affected pathways/functions. Target mRNA information was available for at least half of the miRNAs for these two groups (15/30 for the adults and 9/13 for old), resulting in 149 mRNA targets in the adults and 81 mRNA targets in the old group. The top five most significant molecular functions affected by the mRNA targets of the sexually dimorphic miRNAs were identical between the two age groups: cell death and survival, cellular growth and proliferation, cellular development, gene expression, and cell cycle. These mRNA targets also appeared to be involved in liver toxicology functions, with more significant effects in the adults (significantly affected functions included liver necrosis, liver proliferation, liver fibrosis, hepatocellular carcinoma, and liver hyperplasia) than old rats (significantly affected functions included only liver necrosis).
Hepatic miRNA expression differences between adult and old rats
As shown in Fig. 5c and Additional file 3, there were 25 miRNAs whose hepatic expression differed between adult (5–21 weeks of age) and old (52–104 weeks of age) rats, with 17 different between adult males and old males, and 14 different between adult females and old females. The expression of 6 of these miRNAs differed between adult and old rats in both sexes; the expression of 11 differed only in male rats, and the expression of 8 differed only in female rats. Table 1 shows these age different miRNAs in males and females, along with the target mRNAs and significantly affected pathways/functions. Target mRNA information was available for about half of the miRNAs (8/17 for the males and 8/14 for the females), resulting in 109 mRNA targets for males and 115 mRNA targets for females. The top five molecular functions were identical for both the males and the females (cell death and survival, cellular growth and proliferation, gene expression, cellular development, and cell cycle), and were the same as found for the sex differences in adults and sex differences in old rats. The liver toxicological function differences between adult and old rats were similar in each sex, with liver necrosis and liver proliferation functions found in both sexes. There were additional significant liver toxicological function differences between adult and old males that were not seen between adult and old females, and these included hepatocellular carcinoma, liver hyperplasia, and liver fibrosis. The similarity in pathway/function differences between adult and old in both sexes may be the result of the large overlap in mRNA targets; 47 mRNA targets were in common between the sexes.
Association of hepatic miRNA expression with liver histopathology
Liver histology sections were examined for standard pathological findings at 52, 78, and 104 weeks of age. Sex bias in the incidence of three types of lesions was observed (Table 2). These included basophilic focus, hyperplasia of the bile duct, and vacuolization in hepatocyte cytoplasm. Examples of the histological presentation of these lesions are shown in Fig. 7. Female-biased incidence of basophilic foci was evident at 78 weeks (3 of 8 females showed the lesion compared to 0 of 8 males) and 104 weeks (12 of 16 females compared to 0 of 14 males). As is typical for histopathological findings of basophilic foci, no severity score was provided . Basophilic foci are considered pre-neoplastic lesions; however, no neoplasms were observed in any female animals. One neoplastic lesion (classified as hepatocellular adenoma), however, was observed in a male animal at 104 weeks of age. This is consistent with the slightly higher incidence rate for hepatocellular adenoma in normal males (2.3%) compared to females (0.44%) at two years of age, according to National Toxicology Program historical control animal data for F344 rats . Male-biased incidence of hyperplasia of the bile duct was observed at 52 weeks of age (8 of 10 males compared to 0 of 10 females), 78 weeks of age (7 of 8 males compared to 2 of 8 females), and 104 weeks of age (10 of 14 males compared to 2 of 16 females). Severity scores among animals showing the lesion ranged from minimal to mild in all cases. Male-specific incidence of hepatocyte cytoplasm vacuolization was observed at 52 weeks of age (5 of 10 males), 78 weeks of age (2 of 8 males), and 104 weeks of age (1 of 14 males) compared to zero incidence in females of the same ages. Average severity scores for this lesion also ranged from minimal to mild, with one case of moderate (78 weeks) and one case of marked (104 weeks). All findings were considered common for F344 rats of these ages.
As can be seen in Table 2, there is a substantial sex difference in the prevalence of basophilic foci in the livers of 78- and 104-week-old rats, with females showing this histopathology while males did not. The histopathology was not observed in males or females of 52 weeks of age, or in males of 78 or 104 weeks of age. Basophilic foci were not found (0/12 males and 0/12 females at 6–7 week of age and 0/24 males and 0/24 females at 17–18 weeks of age) in contemporaneous studies at NCTR using F344 rats from the same colony and living in the same environmental conditions at NCTR . Using the expression level from 5 to 21 weeks of age (adults) as a baseline where no liver histopathology was observed, the fold change in expression of each miRNA was calculated at 52, 78, and 104 weeks of age for each sex separately. miRNAs were excluded if the expression value was less than 0.3 (60th percentile of expression level distribution; approximately 40% of expression values were undetectable) in all age groups. Significantly changed expression was defined as p < 0.05 and fold change ≥ 1.5 between the adult age group and 52, 78, or 104 weeks of age groups for each sex independently. A total of 43 miRNAs met these criteria and are shown in Additional file 5. The expression of two miRNAs was associated with the presence and absence of basophilic foci, miR-18a, and miR-451 (Table 3). The expression of miR-18a was reduced relative to the 5–21-week-old baseline group when basophilic foci were present, i.e., significantly decreased expression was seen in the livers of 78 and 104 week old (1.7-fold decrease) females and not in any of the males or 52-week-old female livers. The expression of miR-451 was significantly increased relative to the 5–21-week-old baseline group when basophilic foci were present, i.e., increased expression was seen in the livers of 78 (1.8-fold increase) and 104-week-old (1.5 fold increase) females and not in any of the males or 52-week-old female livers.
Table 2 also shows a substantial sex difference in the prevalence of bile duct hyperplasia, with >80% of males from 52 to 104 weeks of age having this histopathology. It was not seen in females of 52 weeks of age, and occurred in <20% of females of 78 and 104 weeks of age. This lesion was not observed (0/12 males and 0/12 females at 6–7 week of age and 0/24 males and 0/24 females at 17–18 weeks of age) in contemporaneous studies at NCTR using F344 rats from the same colony and living in the same environmental conditions at NCTR . Using the same analysis described above, the expression of 2 miRNAs was found to associate with the presence and absence of liver bile duct hyperplasia, miR-99a and miR-203 (Table 4). The expression of both miRNAs was significantly increased in the livers of males from 52 to 104 weeks of age relative to the expression in the livers of the 5–21-week-old baseline group (1.6- to 2.5-fold). The expression of the 2 miRNAs was not significantly changed in the livers of females of 52 to 104 weeks of age.
General trends in the expression of liver miRNAs during the rat lifespan
Age and sex differences were evaluated independently to identify liver miRNA expression changes during the rat life span. The 214 DEMs included 212 miRNAs which showed significant age differences and 65 miRNAs exhibiting significant sex differences. PCA (Fig. 1) and k-means cluster analysis (Fig. 2) provided visualization of the variability of individuals within age/sex groups, between groups, and among the individual miRNAs. Three general observations followed from these analyses. First, age plays a much larger role than sex in these hepatic miRNA expression changes. Notable difference was observed between 2-week-old animals and all subsequent age groups in the PCA (Fig. 1), and aside from a small sex difference observed between 2-week-old males and females, no large-scale patterns of sex difference were observed. Second, the notable separation in the PCA between 2-week-old animals and all subsequent age groups was found to be caused by 38 two-week-specific miRNAs (Fig. 2 and Additional file 4). Third, aside from miRNAs showing 2-week-specific expression, liver miRNAs exhibited relatively stable expression across both sex and age, with >70% of the miRNAs having coefficient of variation <50% (Fig. 2). These differential expression analyses and data visualization methods revealed key insights into the data that helped focus pathway analysis of the DEMs.
Sex-biased miRNAs linked to key liver pathways
Sex differences in miRNA expression were evaluated at each age group (Fig. 3). The distribution of sex-biased miRNAs identified in the liver showed a similar pattern as the sex-biased miRNAs observed previously in kidney tissues  from the same animal study. Both kidney and liver showed significant sex-biased miRNA expression at 2 weeks of age and then again at ages of 15 weeks and older, with a gap of no observed sex-biased miRNAs between 5 and 6 weeks of age. Previously reported liver mRNA expression study revealed a similar pattern of sex-biased regulation, where 5- and 6-week-old rats showed lower numbers of sex-biased mRNAs compared to neighboring age groups, including the 8–21-week-old groups . This suggests either that sex-biased gene expression change is not as crucial during the pubertal developmental stage compared to other life stages, or that the more subtle or fine-tuned miRNA regulation may not be involved in the hormonally controlled sexually dimorphic expression changes that do exist during this developmental time frame. Among the 65 sex-biased DEMs, only 6 (miR-92a, miR-221, miR-374, miR-505*, miR-532-3p, and miR-652) are encoded on the X-chromosome. Four of these miRNAs exhibited female-biased expression, and two exhibited male biased expression (miR-92a and miR-221). miR-221 showed male-biased expression at 8, 15, and 21 weeks of age, exhibiting a 14-fold difference in expression (21 weeks) compared to females. This was the second highest fold change difference among all sex-biased miRNAs, with only male biased miR-183 having a greater sex-biased expression (17-fold). As there was no consistency in the sex-biased expression of the miRNAs, it is not clear what role X-chromosome inactivation may play in the control of sex-biased miRNA expression. No miRNAs are encoded on the rodent Y-chromosome, and only two miRNAs are Y-chromosome-encoded in humans .
The miRNAs with the most consistent female-biased expression were miR-29b (Fig. 3) and miR-152 (Fig. 3 and Fig. 4) with both showing an approximate 2-fold higher expression in females than males at 15, 21, 52, and 78 weeks of age. miR-29b is a key repressor of renal and pulmonary fibrosis [56, 57] and has recently been shown to inhibit hepatic fibrogenesis by directly targeting PIK3R1 and AKT3 in hepatic stellate cells . The role of the miR-29 family, including miR-29b, in fibrosis disease has been recently reviewed and additional targets appear to be collagens . miR-152 is dysregulated in rat models of NAFLD  and has been characterized to act as a tumor suppressor . Thus, these sexually dimorphic miRNAs may play a role in differential susceptibility to liver pathology and disease.
miRNAs showing the most male-biased expression included miR-125b-5p, miR-99a*, miR-99a, miR-96, miR-221, and miR-183 (Fig. 3). miR-125b-5p and miR-99a showed the most consistent male-biased expression with sex differences observed at 15, 21, 52, and 78 weeks of age and fold changes ranging from 1.5 to 2.3. miR-99a* exhibited similar significant male-biased expression of 1.7- to 2.3-fold at 15, 52, and 78 weeks of age. miR-125b-5p and miR-99a are suggested to be co-expressed in humans from chromosome 21 and play a role in vincristine resistance of leukemia cells . miR-221 and miR-96 showed relatively larger sex differences over the middle portion of the life span, with miR-221 expressed at 3-to 14.3-fold higher levels in males than females from 8 to 21 weeks of age, and miR-96 expressed at 2.2 to 7.2 higher levels in males than females from 21 to 78 weeks of age. miR-125b-5p and miR-221 play a role in glioma carcinogenesis and have been recently proposed as prognostic markers in assessing glioma tumors . The expression of miR-183 exhibited the largest identified sex difference with 17.2-fold higher expression in males compared to females at 104 weeks of age. miR-183 and miR-96 are co-transcribed from the miR-183-96-182 cluster found on chromosome 4 in rats and chromosome 7 in humans , suggesting that post-transcriptional processing may independently control the expression of these miRNAs. MiR-183 has been shown to function in a pro-oncogenic manner by suppressing apoptosis and promoting proliferation in esophageal cancer cells .
The expression of 30 miRNAs was sexually dimorphic in the liver of adult (5 to 21 weeks of age) rats and there is evidence that they are involved in the regulation of 149 mRNAs (Table 1). The top five molecular functions associated with these mRNAs (cellular growth and proliferation, cell death and survival, cellular development, gene expression, and cell cycle) represent important pathways whose disruption or modification may be expected to contribute to disease susceptibility. In fact, by focusing on toxicological functions using the IPA database, these mRNAs appear to be involved in multiple processes related to liver disease, including necrosis, proliferation, fibrosis, hepatocellular carcinoma, and hyperplasia. Thus, this analysis suggests the sexually dimorphic expression of miRNAs described in this study may be reflective of sex differences in liver physiology and disease susceptibility.
In older rats from 52 to 104 weeks of age, the hepatic expression of 13 miRNAs was found to be sexually dimorphic, with evidence that they are involved in the regulation of 81 mRNAs (Table 1). The top five molecular functions associated with these mRNAs are identical to those described above for the sexually dimorphic mRNA targets in adults. Despite only seven sexually dimorphic miRNAs being in common between the adult and the old rats (miR-125b-5p, miR-152, miR-29b, miR-374, miR-96, miR-99a*, and miR-99a), 91% (74/81) of the mRNA targets in the old rats were also mRNA targets in the adult rats. This suggests that multiple miRNAs target the same mRNAs resulting in a similar impact on molecular functions. The toxicological functions of these mRNAs, however, were more limited in old rats compared to the adults. Only one significant toxicological function was associated with the 81 mRNA targets, liver necrosis. Thus, sexually dimorphic hepatic miRNA expression in old rats is mostly a subset of that in adults. The reduction in sexually dimorphic miRNA expression seen in this study is consistent with the finding of increased similarity of the male and female hepatic mRNA expression pattern at the oldest ages .
Two-week-specific miRNAs linked to Dlk1-Dio3 cluster
A prominent feature of the k-means cluster (Fig. 2 and Additional file 4) is the high expression of 38 miRNAs in clusters 2 and 3 at 2 weeks of age with no to little expression at subsequent ages. A review of the miRNA annotations which included chromosome locations revealed a notable trend. Of the 38 miRNAs showing 2-week-specific expression, 31 of them (81%) were encoded on rat chromosome 6. Closer investigation revealed that these chromosome 6-encoded miRNAs were members of the delta-like 1 homolog-deiodinase, iodothyronine 3 (Dlk1-Dio3) miRNA cluster. The conserved DLK1-DIO3 miRNA cluster (located on 14q32 in humans) is the largest miRNA cluster in the human genome and includes 54 encoded miRNAs . The Dlk1-Dio3 genomic region is an imprinted region containing two differentially methylated regions which control its regulation . A subset of miRNAs within the Dlk1-Dio3 miRNA cluster, called the miR-379/410 cluster (containing 38 miRNAs), has been shown to be essential for neonatal liver energy homeostasis (hepatic glucose and lipid metabolism) in mice . miR-379, miR-410, and 16 other miRNAs from this subcluster were among the 2-week-specific miRNAs in our study and were present in k-means clusters 2 or 3. Their association with neonatal energy homeostasis is consistent with their expression at 2 weeks of age, when animals are transitioning from maternal energy sources to independent feeding. miRNA members of the Dlk1-Dio3 cluster have also been shown to be activated in embryonic stem cells  and positively correlate with stem cell pluripotency and proliferation . In a mouse model of myocardial infarction, 29 miRNAs originating from the Dlk1-Dio3 cluster were activated in cardiomyocytes after induced tissue damage , suggesting their role in tissue remodeling efforts to replace damaged cardiac tissue. Taken together, these findings suggest a role for Dlk1-Dio3 cluster miRNAs in liver proliferation, differentiation, and energy homeostasis at 2 weeks of age, consistent with the potentially affected pathways/functions (e.g., cellular growth and proliferation, cellular development, cellular growth and maintenance) identified from possible mRNA targets of these miRNAs (Table 1).
In addition, the miRNAs expressed predominantly at 2 weeks of age have predicted mRNA targets that impact the hematological system (Table 1). These three target mRNAs (PRDM1, BCL6, and XBP1) are important transcription factors and regulators of B and T cell differentiation, with all three encoded proteins interacting with each other [72, 73]. This is consistent with the mid-gestation to early juvenile expression of hematopoiesis in the liver, which ends by postnatal day 35 . Thus, the high level expression of these miRNAs at 2 weeks of age and the dramatically reduced levels by 5 weeks of age correlates with the hematopoietic potential of the liver.
miRNA hepatic expression differences between adult and old rats
The expression of a small number of hepatic miRNAs significantly differed between adult (5–21 weeks of age) and old (52–104 weeks of age) rats, 17 in males and 14 in females. The expression of six (miR-29a, miR-29b, miR-29c, miR-34a, miR-375, and miR-466b-2*) was altered in both males and females (Table 1). There were 109 potential target mRNAs for the 17 male miRNAs and 115 for the 14 female miRNAs, with 77 (~42%) found in both males and females. This high degree of overlap of the mRNA targets resulted in the prediction of similar molecular function changes as both males and females matured from adults to old age. Three of the five most significantly affected functions predicted to be affected included cellular death and survival, cellular growth and proliferation, and cell cycle. The miR-29 family and miR-34a, whose expression was altered in both older males and older females, are relatively well-studied miRNAs, and their roles in cell proliferation and apoptosis are established [75, 76]. In addition, the miR-29 family members appear to play a significant role in the development of liver fibrosis, possibly through the regulation of collagen expression . Thus, miRNA control of important functions of liver maintenance appears to be altered in aged males and females, and may contribute to the age-associated susceptibility to liver disease, including fibrosis .
miRNA expression associates with histopathology
Tables 3 and 4 show the expression of specific miRNAs associated with the presence or absence of specific liver histopathology. Increased expression of miR-451 and decreased expression of miR-18a were associated with the presence of hepatic basophilic loci, which are considered pre-neoplastic lesions . miR-451 has recently been suggested as a prognostic biomarker of hepatocellular carcinoma (HCC), with decreased expression in HCC tissues correlating with advanced stage, metastasis, and worse disease-free or overall survival . In vitro studies demonstrated inhibition of cell growth, possibly through direct suppression of the NF-κB pathway  suggesting that the increased expression associated with the presence of basophilic foci, observed in the present study, may be a compensatory response to the pre-neoplastic lesion. In addition, increased serum levels of miR-451 were recently found in patients with NAFLD . Altered expression of miR-18a has also been implicated in the development of HCC through modulation of the estrogen receptor (ER) alpha protein  with increased expression promoting the downregulation of the ER alpha. Because the estrogen pathway appears to be protective for HCC , this may be a key factor in explaining the increased susceptibility of males to HCC. The decreased expression of miR-18a found to be associated with basophilic foci in the present study may, again, be a compensatory response to the pre-neoplastic basophilic foci.
Increased expression of both miR-99a and miR-203 was associated with the presence of hepatic bile duct hyperplasia, a common aging lesion . miR-99a appears to be a tumor suppressor in the liver  that possibly acts through inhibition of the RNA-induced silencing complex, Ago2 . Several studies suggest that miR-203 also acts as a tumor suppressor in the liver [86–88] by inhibition of liver cell proliferation. Altered expression of miRNA-203 has also been found to be associated with alcoholic steatohepatitis . Thus, even though the exact role of these miRNAs in the development of hepatic basophilic foci and bile duct hyperplasia is not known, as described above, these miRNAs have been shown in a variety of in vitro, in vivo, and bioinformatics studies to have a fundamental impact on liver cell biology. Additional targeted studies will be needed to evaluate these miRNAs as possible biomarkers for these liver lesions and their possible roles in development of age- and sex-dependent susceptibility to liver pathology.
Association of miRNA and mRNA expression
miRNAs have been shown to regulate protein expression by a diverse set of mechanisms, some of which destabilize mRNA targets . Combining miRNA and mRNA expression levels may allow a fuller understanding of the influence of miRNAs on mRNA expression in the rat liver. mRNA expression in the livers of rats from this same study, measured by microarrays, has been previously published  and are available in the gene expression omnibus database (GSE21335). Preliminary analysis of the expression of mRNA targets listed in Table 1 shows that 226 of the 464 mRNA targets (49%) were expressed at detectable levels. Of these detectable mRNA targets, the expression of an average of 26% of mRNAs in each class was significantly (p < 0.05) changed with an inverse relationship to the miRNA expression change; that is, when the miRNA was significantly upregulated, the target mRNA was significantly downregulated, or vice versa. Thus, the stability of a substantial fraction of the detectable target mRNAs appears to be directly affected by the miRNAs listed in Table 1.
The relationship between miRNAs and their targets, however, is complex and not completely understood . For example, multiple miRNAs can affect a single mRNA, and a single miRNA can affect multiple mRNAs. Since the miRNAs may control protein expression by mRNA degradation and/or inhibition of translation, there may not be a direct inverse relationship between miRNA levels and target mRNA levels. In addition, there may be other factors that affect expression, such as DNA methylation, transcription factors, and histone modifications. Also, the current state of knowledge regarding miRNA-mRNA relationships in the liver is incomplete. The relative lack of sensitivity of the older microarray technology used to measure mRNA levels  may be responsible for the low detectability of the predicted target mRNAs. However, better quantitative and sensitive methods, such as RNA-Seq, may allow more accurate evaluation of the correlation of miRNA expression with target mRNA expression. Additional confirmatory computational and experimental analyses, including protein quantitation, will be required to fully understand the role of miRNAs and other gene regulatory mechanisms in the control of hepatic gene expression.
As a key first step in understanding epigenetic mechanisms that may influence sex- and age-related susceptibilities to adverse liver events, genome-wide characterization of liver miRNA expression profiles in both sexes during the rat life span was carried out. miRNAs exhibited largely stable expression between sexes and across the rat life span with the exception of a conspicuous pattern of 38 miRNAs showing 2-week-specific expression. This pattern at 2 weeks of age was explained in large part by co-localization of 31 miRNAs to the Dlk1-Dio3 cluster on chromosome 6. This cluster is putatively associated with early developmental proliferation and differentiation, and postnatal liver energy homeostasis, which is consistent with the developmental events occurring in immature animals at 2 weeks of age. Pathway analysis of sex-biased miRNAs in adult and old rats identified sexually dimorphic molecular functions and toxicological functions that may reflect sex differences in liver physiology and disease susceptibility. The expression of specific miRNAs (miR-451, miR-18a, miR-99a, and miR-203) was found to associate with sexually dimorphic histopathology findings of basophilic foci and bile duct hyperplasia. As rats mature from adults to old-age, miRNAs involved in cell death, cell proliferation, and cell cycle (miR-29 family and miR-34a) were found to change expression. In addition, miR-29 family members may be involved in susceptibility to fibrosis. Thus, the expression of miRNAs involved in the control of important liver functions was shown to change during the rat life span and between the sexes, and may contribute to the age- and sex-associated susceptibility to liver disease.
Differentially expressed miRNAs
Drug-induced liver injury
Delta-like 1 homolog-deiodinase, iodothyronine 3
False discovery rate
Ingenuity Pathway Analysis
Non-alcoholic fatty liver disease
National Center for Toxicological Research
Principal component analysis
RNA integrity numbers
Total gene signal
Perucca E, Beghi E, Dulac O, Shorvon S, Tomson T. Assessing risk to benefit ratio in antiepileptic drug therapy. Epilepsy Res. 2000;41(2):107–39.
Anderson GD. Children versus adults: pharmacokinetic and adverse-effect differences. Epilepsia. 2002;43 Suppl 3:53–9.
Cates M, Powers R. Concomitant rash and blood dyscrasias in geriatric psychiatry patients treated with carbamazepine. Ann Pharmacother. 1998;32(9):884–7.
Harms SL, Eberly LE, Garrard JM, Hardie NA, Bland PC, Leppik IE. Prevalence of appropriate and problematic antiepileptic combination therapy in older people in the nursing home. J Am Geriatr Soc. 2005;53(6):1023–8.
Johnson TN. The development of drug metabolising enzymes and their influence on the susceptibility to adverse drug reactions in children. Toxicology. 2003;192(1):37–48.
Chalasani N, Bjornsson E. Risk factors for idiosyncratic drug-induced liver injury. Gastroenterology. 2010;138(7):2246–59.
Njoku DB, Greenberg RS, Bourdi M, Borkowf CB, Dake EM, Martin JL, et al. Autoantibodies associated with volatile anesthetic hepatitis found in the sera of a large cohort of pediatric anesthesiologists. Anesth Analg. 2002;94(2):243–9.
Cho J, Kim L, Li Z, Rose NR, Talor MV, Njoku DB. Sex bias in experimental immune-mediated, drug-induced liver injury in BALB/c mice: suggested roles for Tregs, estrogen, and IL-6. PLoS One. 2013;8(4):e61186.
Russo MW, Galanko JA, Shrestha R, Fried MW, Watkins P. Liver transplantation for acute liver failure from drug induced liver injury in the United States. Liver Transpl. 2004;10(8):1018–23.
Ostapowicz G, Fontana RJ, Schiodt FV, Larson A, Davern TJ, Han SH, et al. Results of a prospective study of acute liver failure at 17 tertiary care centers in the United States. Ann Intern Med. 2002;137(12):947–54.
Amacher DE. Female gender as a susceptibility factor for drug-induced liver injury. Hum Exp Toxicol. 2013;33(9):928–39.
Andrade RJ, Lucena MI, Fernandez MC, Pelaez G, Pachkoria K, Garcia-Ruiz E, et al. Drug-induced liver injury: an analysis of 461 incidences submitted to the Spanish registry over a 10-year period. Gastroenterology. 2005;129(2):512–21.
Carey EJ, Vargas HE, Douglas DD, Balan V, Byrne TJ, Harrison ME, et al. Inpatient admissions for drug-induced liver injury: results from a single center. Dig Dis Sci. 2008;53(7):1977–82.
Lucena MI, Andrade RJ, Kaplowitz N, Garcia-Cortes M, Fernandez MC, Romero-Gomez M, et al. Phenotypic characterization of idiosyncratic drug-induced liver injury: the influence of age and sex. Hepatology. 2009;49(6):2001–9.
Kwekel JC, Desai VG, Moland CL, Branham WS, Fuscoe JC. Age and sex dependent changes in liver gene expression during the life cycle of the rat. BMC Genomics. 2010;11:675.
Lee Y, Jeon K, Lee JT, Kim S, Kim VN. MicroRNA maturation: stepwise processing and subcellular localization. EMBO J. 2002;21(17):4663–70.
Macrae IJ, Zhou K, Li F, Repic A, Brooks AN, Cande WZ, et al. Structural basis for double-stranded RNA processing by Dicer. Science. 2006;311(5758):195–8.
Kawamata T, Making TY, RISC. Trends Biochem Sci. 2010;35(7):368–76.
Fabian MR, Sonenberg N, Filipowicz W. Regulation of mRNA translation and stability by microRNAs. Annu Rev Biochem. 2010;79:351–79.
Palmero EI, de Campos SG, Campos M, de Souza NC, Guerreiro ID, Carvalho AL, et al. Mechanisms and role of microRNA deregulation in cancer onset and progression. Genet Mol Biol. 2011;34(3):363–70.
Garofalo M, Croce CM. microRNAs: master regulators as potential therapeutics in cancer. Annu Rev Pharmacol Toxicol. 2011;51:25–43.
Wei J, Xie L, Taron M, Rosell R, Liu B. Epigenetic alterations of tumor marker microRNAs: towards new cancer therapies. Drug News Perspect. 2010;23(10):655–61.
Kerr TA, Korenblat KM, Davidson NO. MicroRNAs and liver disease. Transl Res. 2011;157(4):241–52.
Bala S, Marcos M, Szabo G. Emerging role of microRNAs in liver diseases. World J Gastroenterol. 2009;15(45):5633–40.
Reid G, Kirschner MB, van Zandwijk N. Circulating microRNAs: association with disease and potential use as biomarkers. Crit Rev Oncol Hematol. 2011;80(2):193–208.
Bihrer V, Friedrich-Rust M, Kronenberger B, Forestier N, Haupenthal J, Shi Y, et al. Serum miR-122 as a biomarker of necroinflammation in patients with chronic hepatitis C virus infection. Am J Gastroenterol. 2011;106(9):1663–9.
Cermelli S, Ruggieri A, Marrero JA, Ioannou GN, Beretta L. Circulating microRNAs in patients with chronic hepatitis C and non-alcoholic fatty liver disease. PLoS One. 2011;6(8):e23937.
Fukushima T, Hamada Y, Yamada H, Horii I. Changes of micro-RNA expression in rat liver treated by acetaminophen or carbon tetrachloride--regulating role of micro-RNA for RNA expression. J Toxicol Sci. 2007;32(4):401–9.
Wang K, Zhang S, Marzolf B, Troisch P, Brightman A, Hu Z, et al. Circulating microRNAs, potential biomarkers for drug-induced liver injury. Proc Natl Acad Sci U S A. 2009;106(11):4402–7.
Salminen WF, Yang X, Mendrick DL. Using microRNA as biomarkers of drug-induced liver injury. J Mol Biomark Diagn. 2011;2(5):1–8.
Vrijens K, Bollati V, Nawrot TS. microRNAs as potential signatures of environmental exposure or effect: a systematic review. Environ Health Perspect. 2015;123(5):399–411.
Hornby RJ, Starkey Lewis P, Dear J, Goldring C, Park BK. MicroRNAs as potential circulating biomarkers of drug-induced liver injury: key current and future issues for translation to humans. Expert Rev Clin Pharmacol. 2014;7(3):349–62.
Li N, Muthusamy S, Liang R, Sarojini H, Wang E. Increased expression of miR-34a and miR-93 in rat liver during aging, and their impact on the expression of Mgst1 and Sirt1. Mech Ageing Dev. 2011;132(3):75–85.
Kwekel JC, Vijay V, Desai VG, Moland CL, Fuscoe JC. Age and sex differences in kidney microRNA expression during the life span of F344 rats. Biol Sex Differ. 2015;6(1):1.
Zhang J, Liu Q, Zhang W, Li J, Li Z, Tang Z, et al. Comparative profiling of genes and miRNAs expressed in the newborn, young adult, and aged human epididymides. Acta Biochim Biophys Sin (Shanghai). 2010;42(2):145–53.
Sharma S, Eghbali M. Influence of sex differences on microRNA gene regulation in disease. Biol Sex Differ. 2014;5(1):3.
de Lencastre A, Pincus Z, Zhou K, Kato M, Lee SS, Slack FJ. MicroRNAs both promote and antagonize longevity in C. elegans. Curr Biol. 2010;20(24):2159–68.
Kato M, Chen X, Inukai S, Zhao H, Slack FJ. Age-associated changes in expression of small, noncoding RNAs, including microRNAs, in C. elegans. RNA. 2011;17(10):1804–20.
Smith-Vikos T, Slack FJ. MicroRNAs and their roles in aging. J Cell Sci. 2012;125(Pt 1):7–17.
Maes OC, An J, Sarojini H, Wang E. Murine microRNAs implicated in liver functions and aging process. Mech Ageing Dev. 2008;129(9):534–41.
Morgan CP, Bale TL. Sex differences in microRNA regulation of gene expression: no smoke, just miRs. Biol Sex Differ. 2012;3(1):22.
Waxman DJ, O’Connor C. Growth hormone regulation of sex-dependent liver gene expression. Mol Endocrinol. 2006;20(11):2613–29.
Zhang Y, Wu L, Wang Y, Zhang M, Li L, Zhu D, et al. Protective role of estrogen-induced miRNA-29 expression in carbon tetrachloride-induced mouse liver injury. J Biol Chem. 2012;287(18):14851–62.
Chen PJ, Yeh SH, Liu WH, Lin CC, Huang HC, Chen CL, et al. Androgen pathway stimulates microRNA-216a transcription to suppress the tumor suppressor in lung cancer-1 gene in early hepatocarcinogenesis. Hepatology. 2012;56(2):632–43.
Liang QL, Chen GQ, Li ZY, Wang BR. Function and histopathology of a cell adhesion molecule TSLC1 in cancer. Cancer Invest. 2011;29(2):107–12.
Goodman DG, Maronpot RR, Newberne PM, Popp JA, Squire RA. Proliferation and selected other lesions in the liver of rats. Guides for toxologic pathology. Washington: STP/ARP/AFIP; 1994.
National Toxicology Program. Toxicology and carcinogenesis studies of acrylamide (CASRN 79-06-1) in F344/N rats and B6C3F1 mice (feed and drinking water studies). Natl Toxicol Program Tech Rep Ser. 2012;(575):1-234.
Fang H, Harris SC, Su Z, Chen M, Qian F, Shi L, et al. ArrayTrack: an FDA and public genomic tool. Methods Mol Biol. 2009;563:379–98.
Tong W, Cao X, Harris S, Sun H, Fang H, Fuscoe J, et al. ArrayTrack--supporting toxicogenomic research at the U.S. Food and Drug Administration National Center for Toxicological Research. Environ Health Perspect. 2003;111(15):1819–26.
FDA. ArrayTrack- a microarray database, data analysis and interpretation tool 2016. Available from: http://www.fda.gov/ScienceResearch/BioinformaticsTools/Arraytrack/default.htm. Accessed 10 Nov 2016.
Vlachos IS, Paraskevopoulou MD, Karagkouni D, Georgakilas G, Vergoulis T, Kanellos I, et al. DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. 2015;43(Database issue):D153–9.
Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T. miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 2009;37(Database issue):D105–10.
National Toxicology Program NIEHS, NIH. Nonneoplastic Lesion Atlas, Liver - Focus 2015 [updated June 19, 2014]. Available from: http://ntp.niehs.nih.gov/nnl/hepatobiliary/liver/foci/index.htm. Accessed 10 Nov 2016.
National Toxicology Program NIEHS, NIH. Historical Controls 2015 [updated January 2, 2015]. Available from: http://ntp.niehs.nih.gov/results/dbsearch/historical/index.html. Accessed 10 Nov 2016.
Ghorai A, Ghosh U. miRNA gene counts in chromosomes vary widely in a species and biogenesis of miRNA largely depends on transcription or post-transcriptional processing of coding genes. Front Genet. 2014;5:100.
Qin W, Chung AC, Huang XR, Meng XM, Hui DS, Yu CM, et al. TGF-beta/Smad3 signaling promotes renal fibrosis by inhibiting miR-29. J Am Soc Nephrol. 2011;22(8):1462–74.
Xiao J, Meng XM, Huang XR, Chung AC, Feng YL, Hui DS, et al. miR-29 inhibits bleomycin-induced pulmonary fibrosis in mice. Mol Ther. 2012;20(6):1251–60.
Wang J, Chu ES, Chen HY, Man K, Go MY, Huang XR, et al. microRNA-29b prevents liver fibrosis by attenuating hepatic stellate cell activation and inducing apoptosis through targeting PI3K/AKT pathway. Oncotarget. 2014;6(9):7325–38.
He Y, Huang C, Lin X, Li J. MicroRNA-29 family, a crucial therapeutic target for fibrosis diseases. Biochimie. 2013;95(7):1355–9.
Feng YY, Xu XQ, Ji CB, Shi CM, Guo XR, Fu JF. Aberrant hepatic microRNA expression in nonalcoholic fatty liver disease. Cell Physiol Biochem. 2014;34:1983–97.
Huang H, Hu M, Li P, Lu C, Li M. Mir-152 inhibits cell proliferation and colony formation of CD133 liver cancer stem cells by targeting KIT. Tumour Biol. 2014;36(2):921–8.
Kandi R, Gutti U, Saladi RG, Gutti RK. MiR-125b and miR-99a encoded on chromosome 21 co-regulate vincristine resistance in childhood acute megakaryoblastic leukemia. Hematol Oncol Stem Cell Ther. 2014;8(2):95–7.
Li X, Zheng J, Chen L, Diao H, Liu Y. Predictive and prognostic roles of abnormal expression of tissue miR-125b, miR-221, and miR-222 in glioma. Mol Neurobiol. 2014;53(1):577–83.
Li P, Sheng C, Huang L, Zhang H, Cheng Z, Zhu Q. MiR-183/-96/-182 cluster is up-regulated in most breast cancers and increases cell proliferation and migration. Breast Cancer Res. 2014;16(6):473.
Yang M, Liu R, Li X, Liao J, Pu Y, Pan E, et al. miRNA-183 suppresses apoptosis and promotes proliferation in esophageal cancer by targeting PDCD4. Mol Cells. 2014;37(12):873–80.
Benetatos L, Hatzimichael E, Londin E, Vartholomatos G, Loher P, Rigoutsos I, et al. The microRNAs within the DLK1-DIO3 genomic region: involvement in disease pathogenesis. Cell Mol Life Sci. 2013;70(5):795–814.
Benetatos L, Vartholomatos G, Hatzimichael E. DLK1-DIO3 imprinted cluster in induced pluripotency: landscape in the mist. Cell Mol Life Sci. 2014;71(22):4421–30.
Labialle S, Marty V, Bortolin-Cavaille ML, Hoareau-Osman M, Pradere JP, Valet P, et al. The miR-379/miR-410 cluster at the imprinted Dlk1-Dio3 domain controls neonatal metabolic adaptation. EMBO J. 2014;33(19):2216–30.
Gururajan M, Josson S, Chu GC, Lu CL, Lu YT, Haga CL, et al. miR-154* and miR-379 in the DLK1-DIO3 microRNA mega-cluster regulate epithelial to mesenchymal transition and bone metastasis of prostate cancer. Clin Cancer Res. 2014;20(24):6559–69.
Liu L, Luo GZ, Yang W, Zhao X, Zheng Q, Lv Z, et al. Activation of the imprinted Dlk1-Dio3 region correlates with pluripotency levels of mouse stem cells. J Biol Chem. 2010;285(25):19483–90.
Janssen R, Zuidwijk M, Muller A, Mulders J, Oudejans CB, Simonides WS. Cardiac expression of deiodinase type 3 (Dio3) following myocardial infarction is associated with the induction of a pluripotency microRNA signature from the Dlk1-Dio3 genomic region. Endocrinology. 2013;154(6):1973–8.
Boi M, Zucca E, Inghirami G, Bertoni F. PRDM1/BLIMP1: a tumor suppressor gene in B and T cell lymphomas. Leuk Lymphoma. 2015;56(5):1223–8.
Nutt SL, Fairfax KA, Kallies A. BLIMP1 guides the fate of effector B and T cells. Nat Rev Immunol. 2007;7(12):923–7.
Figueiredo JF, Mattix ME, Papenfuss TL. Hematopoietic System. In: Parker GA, Picut CA, editors. Atlas of Histology of the Juvenile Rat. 1st ed. New York: Elsevier Science; 2016. p. 349-71.
Rokavec M, Li H, Jiang L, Hermeking H. The p53/miR-34 axis in development and disease. J Mol Cell Biol. 2014;6(3):214–30.
Wang Y, Zhang X, Li H, Yu J, Ren X. The role of miRNA-29 family in cancer. Eur J Cell Biol. 2013;92(3):123–8.
Kim IH, Kisseleva T, Brenner DA. Aging and liver disease. Curr Opin Gastroenterol. 2015;31(3):184–91.
Huang JY, Zhang K, Chen DQ, Chen J, Feng B, Song H, et al. MicroRNA-451: epithelial-mesenchymal transition inhibitor and prognostic biomarker of hepatocelluar carcinoma. Oncotarget. 2015;6(21):18613–30.
Li HP, Zeng XC, Zhang B, Long JT, Zhou B, Tan GS, et al. miR-451 inhibits cell proliferation in human hepatocellular carcinoma through direct suppression of IKK-beta. Carcinogenesis. 2013;34(11):2443–51.
Yamada H, Suzuki K, Ichino N, Ando Y, Sawada A, Osakabe K, et al. Associations between circulating microRNAs (miR-21, miR-34a, miR-122 and miR-451) and non-alcoholic fatty liver. Clin Chim Acta. 2013;424:99–103.
Li CL, Yeh KH, Liu WH, Chen CL, Chen DS, Chen PJ, et al. Elevated p53 promotes the processing of miR-18a to decrease estrogen receptor-alpha in female hepatocellular carcinoma. Int J Cancer. 2015;136(4):761–70.
Naugler WE, Sakurai T, Kim S, Maeda S, Kim K, Elsharkawy AM, et al. Gender disparity in liver cancer due to sex differences in MyD88-dependent IL-6 production. Science. 2007;317(5834):121–4.
National Toxicology Program NIEHS, NIH. Nonneoplastic Lesion Atlas, Liver, Bile duct- Hyperplasia 2015. Available from: https://ntp.niehs.nih.gov/nnl/hepatobiliary/liver/bdhyperp/liver-bile_duct_hyperplasia_508.pdf. Accessed 26 Jul 2016.
Li D, Liu X, Lin L, Hou J, Li N, Wang C, et al. MicroRNA-99a inhibits hepatocellular carcinoma growth and correlates with prognosis of patients with hepatocellular carcinoma. J Biol Chem. 2011;286(42):36677–85.
Zhang J, Jin H, Liu H, Lv S, Wang B, Wang R, et al. MiRNA-99a directly regulates AGO2 through translational repression in hepatocellular carcinoma. Oncogenesis. 2014;3:e97.
Furuta M, Kozaki KI, Tanaka S, Arii S, Imoto I, Inazawa J. miR-124 and miR-203 are epigenetically silenced tumor-suppressive microRNAs in hepatocellular carcinoma. Carcinogenesis. 2010;31(5):766–76.
Hailer A, Grunewald TG, Orth M, Reiss C, Kneitz B, Spahn M, et al. Loss of tumor suppressor mir-203 mediates overexpression of LIM and SH3 Protein 1 (LASP1) in high-risk prostate cancer thereby increasing cell proliferation and migration. Oncotarget. 2014;5(12):4144–53.
Wei W, Wanjun L, Hui S, Dongyue C, Xinjun Y, Jisheng Z. miR-203 inhibits proliferation of HCC cells by targeting survivin. Cell Biochem Funct. 2013;31(1):82–5.
Chen YP, Jin X, Xiang Z, Chen SH, Li YM. Circulating MicroRNAs as potential biomarkers for alcoholic steatohepatitis. Liver Int. 2013;33(8):1257–65.
Filipowicz W, Sonenberg N. The long unfinished march towards understanding microRNA-mediated repression. RNA. 2015;21(4):519–24.
This research was supported in part by an appointment to the Research Participation Program at the National Center for Toxicological Research administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and the U.S. Food and Drug Administration (VV).
Funding for this research was from appropriated funds to the U.S. Food and Drug Administration. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the FDA.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the Gene Express Omnibus repository, http://www.ncbi.nlm.nih.gov/geo/, Accession ID: GSE65169, and in the Additional files.
JCK carried out the microarray experiments and helped draft the manuscript. VV performed data processing, normalization, and statistical analysis. VGD participated in the study’s conception and in the design and coordination of the in-life study. CLM participated in the in-life study and RNA isolation. JCF conceived, designed and coordinated the study, analyzed and interpreted the data, and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
212 differentially expressed miRNAs by age for each sex. The age contrast with differential expression is given along with fold-ratio, ANOVA FDR, degrees of freedom (df), ANOVA F statistic, and p values for both males and females. (XLSX 232 kb)
214 miRNAs differentially expressed by age and/or sex. Each miRNA is shown as differentially expressed by age, sex, or both age and sex, along with the k-means cluster number from Fig. 2 and chromosome mapping position. (XLSX 20 kb)
Age-related differences in miRNA expression. Identification of miRNAs in Venn diagrams in Fig. 5. (XLSX 26 kb)
38 miRNAs with 2-week-old-specific expression. Identification of 2-week-old-specific miRNAs, along with expression level in each rat liver sample, k-means cluster number (from Fig. 2), differential expression by age and/or sex, and chromosome mapping position. (XLSX 34 kb)
43 significantly differentially expressed miRNAs between 52, 78, and 104 weeks of age males and females compared to 5–21 weeks of age baseline group. Presence or absence of histopathological findings is given with the relative fold change in expression of the miRNAs. (XLSX 10 kb)
About this article
Cite this article
Kwekel, J.C., Vijay, V., Han, T. et al. Sex and age differences in the expression of liver microRNAs during the life span of F344 rats. Biol Sex Differ 8, 6 (2017). https://doi.org/10.1186/s13293-017-0127-9