Sex and age differences in the expression of liver microRNAs during the life span of F344 rats.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Age and sex differences in kidney microRNA expression during the life span of F344 rats.

BACKGROUND: Growing evidence suggests that epigenetic mechanisms of gene regulation may play a role in susceptibilities to specific toxicities and adverse drug reactions. MiRNAs in particular have been shown to be important regulators in cancer and other diseases and show promise as predictive biomarkers for diagnosis and prognosis. In this study, we characterized the global kidney miRNA expression profile in untreated male and female F344 rats throughout the life span. These findings were correlated with sex-specific susceptibilities to adverse renal events, such as male-biased renal fibrosis and inflammation in old age. METHODS: Kidney miRNA expression was examined in F344 rats at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age in both sexes using Agilent miRNA microarrays. Differential expression was determined using filtering criteria of ≥1.5 fold change and ANOVA or pairwise t-test (FDR <5%) to determine significant age and sex effects, respectively. Pathway analysis software was used to investigate the possible roles of these target genes in age- and sex-specific differences. RESULTS: Three hundred eleven miRNAs were found to be expressed in at least one age and sex. Filtering criteria revealed 174 differentially expressed miRNAs in the kidney; 173 and 34 miRNAs exhibiting age and sex effects, respectively. Principal component analysis revealed age effects predominated over sex effects, with 2-week miRNA expression being much different from other ages. No significant sexually dimorphic miRNA expression was observed from 5 to 8 weeks, while the most differential expression (13 miRNAs) was observed at 21 weeks. Potential target genes of these differentially expressed miRNAs were identified. CONCLUSIONS: The expression of 56% of detected renal miRNAs was found to vary significantly with age and/or sex during the life span of F344 rats. Pathway analysis suggested that 2-week-expressed miRNAs may be related to organ and cellular development and proliferation pathways. Male-biased miRNA expression at older ages correlated with male-biased renal fibrosis and mononuclear cell infiltration. These miRNAs showed high representation in renal inflammation and nephritis pathways, and included miR-214, miR-130b, miR-150, miR-223, miR-142-5p, miR-185, and miR-296*. Analysis of kidney miRNA expression throughout the rat life span will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

Sexual dimorphism in the expression of mitochondria-related genes in rat heart at different ages.

Cardiovascular disease (CVD) is the leading cause of mortality worldwide. Moreover, sex and age are considered major risk factors in the development of CVDs. Mitochondria are vital for normal cardiac function, and regulation of mitochondrial structure and function may impact susceptibility to CVD. To identify potential role of mitochondria in sex-related differences in susceptibility to CVD, we analyzed the basal expression levels of mitochondria-related genes in the hearts of male and female rats. Whole genome expression profiling was performed in the hearts of young (8-week), adult (21-week), and old (78-week) male and female Fischer 344 rats and the expression of 670 unique genes related to various mitochondrial functions was analyzed. A significant (p<0.05) sexual dimorphism in expression levels of 46, 114, and 41 genes was observed in young, adult and old rats, respectively. Gene Ontology analysis revealed the influence of sex on various biological pathways related to cardiac energy metabolism at different ages. The expression of genes involved in fatty acid metabolism was significantly different between the sexes in young and adult rat hearts. Adult male rats also showed higher expression of genes associated with the pyruvate dehydrogenase complex compared to females. In young and adult hearts, sexual dimorphism was not noted in genes encoding oxidative phosphorylation. In old rats, however, a majority of genes involved in oxidative phosphorylation had higher expression in females compared to males. Such basal differences between the sexes in cardiac expression of genes associated with energy metabolism may indicate a likely involvement of mitochondria in susceptibility to CVDs. In addition, female rats showed lower expression levels of apoptotic genes in hearts compared to males at all ages, which may have implications for better preservation of cardiac mass in females than in males.

Metzincins and related genes in experimental renal ageing: towards a unifying fibrosis classifier across species.

BACKGROUND: We have previously described a transcriptomic classifier consisting of metzincins and related genes (MARGS) discriminating kidneys and other organs with or without fibrosis from human biopsies. We now apply our MARGS-based algorithm to a rat model of age-associated interstitial renal fibrosis. METHODS: Untreated Fisher 344 rats (n = 76) were sacrificed between 2 to 104 weeks of age. For gene expression studies, we used single colour (Cy3) Agilent Whole Rat Genome 4 × 44k microarrays; 4-5 animals of each sex were profiled at each of the following ages: 2, 5, 6, 8, 15, 21, 78 and 104 weeks. Intensity data were subjected to variance stabilization (www.Partek.com). Data were analysed with ANOVA and other statistical methods. RESULTS: Sixty MARGS were differentially expressed across age groups. More MARGS were differentially expressed in older males than in older females. Principal component analysis showed gene expression induced segregation of age groups by sex from 6 to 104 weeks of age. The expression level of MMP7 correlated best with fibrosis grade. Severity of fibrosis was determined in 20 animals at 78 and 104 weeks of age. Expression values of 15 of 19 genes of the original classifier present on the Agilent array, in conjunction with linear discriminant analysis, was sufficient to correctly classify these 20 samples into non-fibrosis and fibrosis. Overrepresentation of MMP2 protein and CD44 protein in fibrosis was confirmed by immunofluorescence. CONCLUSIONS: Based on these results and our previous work, the MARGS classifier represents a cross-organ and cross-species classifier of fibrosis irrespective of aetiology. This finding provides evidence for a common pathway leading to fibrosis and will help to design a PCR-based clinical test.

Life cycle analysis of kidney gene expression in male F344 rats.

Age is a predisposing condition for susceptibility to chronic kidney disease and progression as well as acute kidney injury that may arise due to the adverse effects of some drugs. Age-related differences in kidney biology, therefore, are a key concern in understanding drug safety and disease progression. We hypothesize that the underlying suite of genes expressed in the kidney at various life cycle stages will impact susceptibility to adverse drug reactions. Therefore, establishing changes in baseline expression data between these life stages is the first and necessary step in evaluating this hypothesis. Untreated male F344 rats were sacrificed at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age. Kidneys were collected for histology and gene expression analysis. Agilent whole-genome rat microarrays were used to query global expression profiles. An ANOVA (p<0.01) coupled with a fold-change>1.5 in relative mRNA expression, was used to identify 3,724 unique differentially expressed genes (DEGs). Principal component analyses of these DEGs revealed three major divisions in life-cycle renal gene expression. K-means cluster analysis identified several groups of genes that shared age-specific patterns of expression. Pathway analysis of these gene groups revealed age-specific gene networks and functions related to renal function and aging, including extracellular matrix turnover, immune cell response, and renal tubular injury. Large age-related changes in expression were also demonstrated for the genes that code for qualified renal injury biomarkers KIM-1, Clu, and Tff3. These results suggest specific groups of genes that may underlie age-specific susceptibilities to adverse drug reactions and disease. This analysis of the basal gene expression patterns of renal genes throughout the life cycle of the rat will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

o-p'-DDT-mediated uterotrophy and gene expression in immature C57BL/6 mice and Sprague-Dawley rats.

1,1,1-Trichloro-2,2-bis(2-chlorophenyl-4-chlorophenyl)ethane (o,p'-DDT) is an organochlorine pesticide and endocrine disruptor known to activate the estrogen receptor. Comprehensive ligand- and species-comparative dose- and time-dependent studies were conducted to systematically assess the uterine physiological, morphological and gene expression responses elicited by o,p'-DDT and ethynyl estradiol (EE) in immature ovariectomized C57BL/6 mice and Sprague-Dawley rats. Custom cDNA microarrays were used to identify conserved and divergent differential gene expression responses. A total of 1256 genes were differentially expressed by both ligands in both species, 559 of which exhibited similar temporal expression profiles suggesting that o,p'-DDT elicits estrogenic effects at high doses when compared to EE. However, 51 genes exhibited species-specific uterine expression elicited by o,p'-DDT. For example, carbonic anhydrase 2 exhibited species- and ligand-divergent expression as confirmed by quantitative real-time PCR. The identification of comparable temporal phenotypic responses linked to gene expression demonstrates that systematic comparative gene expression assessments are valuable for elucidating conserved and divergent estrogen signaling mechanisms in rodent uterotrophy.

Sex differences in kidney gene expression during the life cycle of F344 rats.

BACKGROUND: The kidney functions in key physiological processes to filter blood and regulate blood pressure via key molecular transporters and ion channels. Sex-specific differences have been observed in renal disease incidence and progression, as well as acute kidney injury in response to certain drugs. Although advances have been made in characterizing the molecular components involved in various kidney functions, the molecular mechanisms responsible for sex differences are not well understood. We hypothesized that the basal expression levels of genes involved in various kidney functions throughout the life cycle will influence sex-specific susceptibilities to adverse renal events. METHODS: Whole genome microarray gene expression analysis was performed on kidney samples collected from untreated male and female Fischer 344 (F344) rats at eight age groups between 2 and 104 weeks of age. RESULTS: A combined filtering approach using statistical (ANOVA or pairwise t test, FDR 0.05) and fold-change criteria (>1.5 relative fold change) was used to identify 7,447 unique differentially expressed genes (DEGs). Principal component analysis (PCA) of the 7,447 DEGs revealed sex-related differences in mRNA expression at early (2 weeks), middle (8, 15, and 21 weeks), and late (104 weeks) ages in the rat life cycle. Functional analysis (Ingenuity Pathway Analysis) of these sex-different genes indicated over-representation of specific pathways and networks including renal tubule injury, drug metabolism, and immune cell and inflammatory responses. The mRNAs that code for the qualified urinary protein kidney biomarkers KIM-1, Clu, Tff3, and Lcn2 were also observed to show sex differences. CONCLUSIONS: These data represent one of the most comprehensive in-life time course studies to be published, assessing sex differences in global gene expression in the F344 rat kidney. PCA and Venn analyses reveal specific periods of sexually dimorphic gene expression which are associated with functional categories (xenobiotic metabolism and immune cell and inflammatory responses) of key relevance to acute kidney injury and chronic kidney disease, which may underlie sex-specific susceptibility. Analysis of the basal gene expression patterns of renal genes throughout the life cycle of the rat will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

Characterization of whole genome amplified (WGA) DNA for use in genotyping assay development.

BACKGROUND: Genotyping assays often require substantial amounts of DNA. To overcome the problem of limiting amounts of available DNA, Whole Genome Amplification (WGA) methods have been developed. The multiple displacement amplification (MDA) method using Φ29 polymerase has become the preferred choice due to its high processivity and low error rate. However, the uniformity and fidelity of the amplification process across the genome has not been extensively characterized. RESULTS: To assess amplification uniformity, we used array-based comparative genomic hybridization (aCGH) to evaluate DNA copy number variations (CNVs) in DNAs amplified by two MDA kits: GenomiPhi and REPLI-g. The Agilent Human CGH array containing nearly one million probes was used in this study together with DNAs from a normal subject and 2 cystic fibrosis (CF) patients. Each DNA sample was amplified 4 independent times and compared to its native unamplified DNA. Komogorov distances and Phi correlations showed a high consistency within each sample group. Less than 2% of the probes showed more than 2-fold CNV introduced by the amplification process. The two amplification kits, REPLI-g and GenomiPhi, generate very similar amplified DNA samples despite the differences between the unamplified and amplified DNA samples. The results from aCGH analysis indicated that there were no obvious CNVs in the CFTR gene region due to WGA when compared to unamplified DNA. This was confirmed by quantitative real-time PCR copy number assays at 10 locations within the CFTR gene. DNA sequencing analysis of a 2-kb region within the CFTR gene showed no mutations introduced by WGA. CONCLUSION: The relatively high uniformity and consistency of the WGA process, coupled with the low replication error rate, suggests that WGA DNA may be suitable for accurate genotyping. Regions of the genome that were consistently under-amplified were found to contain higher than average GC content. Because of the consistent differences between the WGA DNA and the native unamplified DNA, characterization of the genomic region of interest, as described here, will be necessary to ensure the reliability of genotyping results from WGA DNA.

Functional comparison of microarray data across multiple platforms using the method of percentage of overlapping functions.

Functional comparison across microarray platforms is used to assess the comparability or similarity of the biological relevance associated with the gene expression data generated by multiple microarray platforms. Comparisons at the functional level are very important considering that the ultimate purpose of microarray technology is to determine the biological meaning behind the gene expression changes under a specific condition, not just to generate a list of genes. Herein, we present a method named percentage of overlapping functions (POF) and illustrate how it is used to perform the functional comparison of microarray data generated across multiple platforms. This method facilitates the determination of functional differences or similarities in microarray data generated from multiple array platforms across all the functions that are presented on these platforms. This method can also be used to compare the functional differences or similarities between experiments, projects, or laboratories.

Age and sex dependent changes in liver gene expression during the life cycle of the rat.

BACKGROUND: Age- and sex-related susceptibility to adverse drug reactions and disease is a key concern in understanding drug safety and disease progression. We hypothesize that the underlying suite of hepatic genes expressed at various life cycle stages will impact susceptibility to adverse drug reactions. Understanding the basal liver gene expression patterns is a necessary first step in addressing this hypothesis and will inform our assessments of adverse drug reactions as the liver plays a central role in drug metabolism and biotransformation. Untreated male and female F344 rats were sacrificed at 2, 5, 6, 8, 15, 21, 52, 78, and 104 weeks of age. Liver tissues were collected for histology and gene expression analysis. Whole-genome rat microarrays were used to query global expression profiles. RESULTS: An initial list of differentially expressed genes was selected using criteria based upon p-value (p < 0.05) and fold-change (+/- 1.5). Three dimensional principal component analyses revealed differences between males and females beginning at 2 weeks with more divergent profiles beginning at 5 weeks. The greatest sex-differences were observed between 8 and 52 weeks before converging again at 104 weeks. K-means clustering identified groups of genes that displayed age-related patterns of expression. Various adult aging-related clusters represented gene pathways related to xenobiotic metabolism, DNA damage repair, and oxidative stress. CONCLUSIONS: These results suggest an underlying role for genes in specific clusters in potentiating age- and sex-related differences in susceptibility to adverse health effects. Furthermore, such a comprehensive picture of life cycle changes in gene expression deepens our understanding and informs the utility of liver gene expression biomarkers.

Tamoxifen-elicited uterotrophy: cross-species and cross-ligand analysis of the gene expression program.

BACKGROUND: Tamoxifen (TAM) is a well characterized breast cancer drug and selective estrogen receptor modulator (SERM) which also has been associated with a small increase in risk for uterine cancers. TAM's partial agonist activation of estrogen receptor has been characterized for specific gene promoters but not at the genomic level in vivo.Furthermore, reducing uncertainties associated with cross-species extrapolations of pharmaco- and toxicogenomic data remains a formidable challenge. RESULTS: A comparative ligand and species analysis approach was conducted to systematically assess the physiological, morphological and uterine gene expression alterations elicited across time by TAM and ethynylestradiol (EE) in immature ovariectomized Sprague-Dawley rats and C57BL/6 mice. Differential gene expression was evaluated using custom cDNA microarrays, and the data was compared to identify conserved and divergent responses. 902 genes were differentially regulated in all four studies, 398 of which exhibit identical temporal expression patterns. CONCLUSION: Comparative analysis of EE and TAM differentially expressed gene lists suggest TAM regulates no unique uterine genes that are conserved in the rat and mouse. This demonstrates that the partial agonist activities of TAM extend to molecular targets in regulating only a subset of EE-responsive genes. Ligand-conserved, species-divergent expression of carbonic anhydrase 2 was observed in the microarray data and confirmed by real time PCR. The identification of comparable temporal phenotypic responses linked to related gene expression profiles demonstrates that systematic comparative genomic assessments can elucidate important conserved and divergent mechanisms in rodent estrogen signalling during uterine proliferation.

Species-specific regulation of PXR/CAR/ER-target genes in the mouse and rat liver elicited by o, p'-DDT.

BACKGROUND: Dichlorodiphenyltrichloroethane (DDT) is a persistent estrogenic organochlorine pesticide that is a rodent hepatic tumor promoter, with inconclusive carcinogenicity in humans. We have previously reported that o, p'-DDT elicits primarily PXR/CAR-mediated activity, rather than ER-mediated hepatic responses, and suggested that CAR-mediated effects, as opposed to ER-mediated effects, may be more important in tumor promotion in the rat liver. To further characterize species-specific hepatic responses, gene expression analysis, with complementary histopathology and tissue level analyses were investigated in immature, ovariectomized C57BL/6 mice treated with 300 mg/kg o, p'-DDT, and compared to Sprague-Dawley rat data. RESULTS: Rats and mice exhibited negligible histopathology with rapid o, p'-DDT metabolism. Gene expression profiles were also similar, exhibiting PXR/CAR regulation with the characteristic induction of Cyp2b10 and Cyp3a11. However, PXR-specific target genes such as Apoa4 or Insig2 exhibited more pronounced induction compared to CAR-specific genes in the mouse. In addition, mouse Car mRNA levels decreased, possibly contributing to the preferential activation of mouse PXR. ER-regulated genes Cyp17a1 and Cyp7b1 were also induced, suggesting o, p'-DDT also elicits ER-mediated gene expression in the mouse, while ER-mediated effects were negligible in the rat, possibly due to the inhibitory effects of CAR on ER activities. In addition, o, p'-DDT induced Gadd45a, Gadd45b and Cdkn1, suggesting DNA damage may be an additional risk factor. Furthermore, elevated blood DHEA-S levels at 12 h after treatment in the mouse may also contribute to the endocrine-related effects of o, p'-DDT. CONCLUSION: Although DDT is known to cause rodent hepatic tumors, the marked species differences in PXR/CAR structure, expression patterns and ligand preference as well as significant species-specific differences in steroidogenesis, especially CYP17A1 expression and activity, confound the extrapolation of these results to humans. Nevertheless, the identification of potential modes of action as well as species-specific responses may assist in the selection and further development of more appropriate models for assessing the toxicity of DDT to humans and wildlife.

o,p'-DDT elicits PXR/CAR-, not ER-, mediated responses in the immature ovariectomized rat liver.

Technical-grade dichlorodiphenyltrichloroethane (DDT) is an agricultural pesticide and malarial vector control agent that has been designated a potential human hepatocarcinogen. The o,p'-enantiomer exhibits estrogenic activity that has been associated with the carcinogenicity of DDT. The temporal and dose-dependent hepatic estrogenicity of o,p'-DDT was investigated using complementary DNA microarrays in immature ovariectomized Sprague-Dawley rats with complementary histopathology and tissue-level analysis. Animals were gavaged with 300 mg/kg o,p'-DDT either once or once daily for 3 consecutive days. Liver samples were examined 2, 4, 8, 12, 18, or 24 h after a single dose or following three daily doses. For dose-response studies, a single dose of 3, 10, 30, 100, or 300 mg/kg body weight o,p'-DTT was administered for 3 consecutive days. Genes associated with drug metabolism (Cyp2b2 and Cyp3a2), the nuclear receptors constitutive androstane receptor (CAR) and pregnane X receptor (PXR), cell proliferation (Ccnd1, Ccnb1, Ccnb2, and Stmn1), and oxidative stress (Gclm and Hmox1) were significantly induced. Cyp2b2 exhibited dose-dependent regulation and was significantly induced across all time points, while cell proliferation- and oxidative stress-related genes exhibited transient induction. The induction of Cyp2b2 and Cyp3a2 mRNA levels suggest PXR/CAR activation, consistent with expression of genes associated with oxidative stress. Few genes known to be estrogen receptor (ER) regulated were differentially expressed when compared to the hepatic gene expression profile elicited by ethynyl estradiol in immature ovariectomized C57BL/6 mice using the same study design and analysis methods. These data indicate that o,p'-DDT elicits PXR/CAR-, not ER-, mediated gene expression in the rat liver. Based on the species-specific differences in CAR regulation, the extrapolation of rodent DDT hepatocarcinogenicity to humans warrants further investigation.

Dioxin induces an estrogen-like, estrogen receptor-dependent gene expression response in the murine uterus.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a ubiquitous environmental contaminant that elicits a broad range of toxicities in a tissue-, sex-, age-, and species-specific manner, including alterations in estrogen signaling. Many, if not all, of these effects involve changes in gene expression mediated via the activation of the aryl hydrocarbon receptor (AhR), a ligand activated transcription factor. Recent data indicate that TCDD may also elicit AhR-mediated estrogenic activity through interactions with the estrogen receptor (ER). In an effort to further characterize the estrogenic activity of TCDD, a comprehensive time-course analysis of uterine gene expression was conducted using ovariectomized C57BL/6 mice. Comparison of the temporal uterine transcriptional response to TCDD with that of ethynyl estradiol (EE) revealed a large proportion of the TCDD-mediated gene expression changes were also responsive to EE. Furthermore, pretreatment of mice with the pure ER antagonist ICI 182 780 (faslodex) inhibited gene expression responses to both EE and TCDD, providing additional evidence that these transcriptional responses involve the ER.

A cross-species analysis of the rodent uterotrophic program: elucidation of conserved responses and targets of estrogen signaling.

Physiological, morphological, and transcriptional alterations elicited by ethynyl estradiol in the uteri of Sprague-Dawley rats and C57BL/6 mice were assessed using comparable study designs, microarray platforms, and analysis methods to identify conserved estrogen signaling networks. Comparative analysis identified 153 orthologous gene pairs that were positively correlated, suggesting conserved transcriptional targets important in uterine proliferation. Functional annotation for these responses were associated with angiogenesis, water and solute transport, cell cycle control, redox control, DNA replication, protein synthesis and transport, xenobiotic metabolism, cell-cell communication, energetics, and cholesterol and fatty acid regulation. The identification of conserved temporal expression patterns of these orthologs provides experimental support for the transfer of functional annotation from mouse orthologs to 44 previously unannotated rat expressed sequence tags based on their homology and co-expression patterns. The identification of comparable temporal phenotypic responses linked to related gene expression profiles demonstrates the ability of systematic comparative genomic assessments to elucidate important conserved mechanisms in rodent estrogen signaling during uterine proliferation.