Search

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 sexrelated 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. [ABSTRACT FROM AUTHOR]

In studies that use DNA arrays to assess changes in gene expression, it is preferable to measure the significance of treatment effects on a group of genes from a pathway or functional category such as gene ontology terms (GO terms, http://www.geneontology.org) because this facilitates the interpretation of effects and may markedly increase significance. A modified meta-analysis method to combine p-values was developed to measure the significance of an overall treatment effect on such functionally-defined groups of genes, taking into account the correlation structure among genes. For hypothesis testing that allows gene expression to change in both directions, p-values are calculated under the null distribution generated by a Monte Carlo method. As a test of this procedure, we attempted to distinguish altered pathways in microarray studies performed with Mitochips, oligonucleotide microarrays specific to mitochondrial DNA-encoded transcripts. We found that our analytic method improves the specificity of selection for altered pathways, due to incorporation of the inter-gene correlation structure in each pathway. It is thus a practical method to measure treatment effects on GO groups. In many actual applications, microarray experiments measure treatment effects under complicated design structures and with small sample sizes. For such applications to real data of limited statistical power, and also in computer simulations, we demonstrate that our method gives reasonable test results. [ABSTRACT FROM AUTHOR]

Background: Environmental ozone can rapidly degrade cyanine 5 (Cy5), a fluorescent dye commonly used in microarray gene expression studies. Cyanine 3 (Cy3) is much less affected by atmospheric ozone. Degradation of the Cy5 signal relative to the Cy3 signal in 2-color microarrays will adversely reduce the Cy5/Cy3 ratio resulting in unreliable microarray data. Results: Ozone in central Arkansas typically ranges between ∼22 ppb to ∼46 ppb and can be as high as 60-100 ppb depending upon season, meteorological conditions, and time of day. These levels of ozone are common in many areas of the country during the summer. A carbon filter was installed in the laboratory air handling system to reduce ozone levels in the microarray laboratory. In addition, the airflow was balanced to prevent non-filtered air from entering the laboratory. These modifications reduced the ozone within the microarray laboratory to ∼2-4 ppb. Data presented here document reductions in Cy5 signal on both in-house produced microarrays and commercial microarrays as a result of exposure to unfiltered air. Comparisons of identically hybridized microarrays exposed to either carbon-filtered or unfiltered air demonstrated the protective effect of carbon-filtration on microarray data as indicated by Cy5 and Cy3 intensities. LOWESS normalization of the data was not able to completely overcome the effect of ozoneinduced reduction of Cy5 signal. Experiments were also conducted to examine the effects of high humidity on microarray quality. Modest, but significant, increases in Cy5 and Cy3 signal intensities were observed after 2 or 4 hours at 98-99% humidity compared to 42% humidity. Conclusion: Simple installation of carbon filters in the laboratory air handling system resulted in low and consistent ozone levels. This allowed the accurate determination of gene expression by microarray using Cy5 and Cy3 fluorescent dyes. [ABSTRACT FROM AUTHOR]

In studies that use DNA arrays to assess changes in gene expression, it is preferable to measure the significance of treatment effects on a group of genes from a pathway or functional category such as gene ontology terms (GO terms, http://www.geneontology.org) because this facilitates the interpretation of effects and may markedly increase significance. A modified meta-analysis method to combine p-values was developed to measure the significance of an overall treatment effect on such functionally-defined groups of genes, taking into account the correlation structure among genes. For hypothesis testing that allows gene expression to change in both directions, p-values are calculated under the null distribution generated by a Monte Carlo method. As a test of this procedure, we attempted to distinguish altered pathways in microarray studies performed with Mitochips, oligonucleotide microarrays specific to mitochondrial DNA-encoded transcripts. We found that our analytic method improves the specificity of selection for altered pathways, due to incorporation of the inter-gene correlation structure in each pathway. It is thus a practical method to measure treatment effects on GO groups. In many actual applications, microarray experiments measure treatment effects under complicated design structures and with small sample sizes. For such applications to real data of limited statistical power, and also in computer simulations, we demonstrate that our method gives reasonable test results.