The concordance of RNA-sequencing (RNA-seq) with microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed using a range of chemical treatment conditions. ...Here we use a comprehensive study design to generate Illumina RNA-seq and Affymetrix microarray data from the same liver samples of rats exposed in triplicate to varying degrees of perturbation by 27 chemicals representing multiple modes of action (MOAs). The cross-platform concordance in terms of differentially expressed genes (DEGs) or enriched pathways is linearly correlated with treatment effect size (R(2)0.8). Furthermore, the concordance is also affected by transcript abundance and biological complexity of the MOA. RNA-seq outperforms microarray (93% versus 75%) in DEG verification as assessed by quantitative PCR, with the gain mainly due to its improved accuracy for low-abundance transcripts. Nonetheless, classifiers to predict MOAs perform similarly when developed using data from either platform. Therefore, the endpoint studied and its biological complexity, transcript abundance and the genomic application are important factors in transcriptomic research and for clinical and regulatory decision making.
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.
The rat has been used extensively as a model for evaluating chemical toxicities and for understanding drug mechanisms. However, its transcriptome across multiple organs, or developmental stages, has ...not yet been reported. Here we show, as part of the SEQC consortium efforts, a comprehensive rat transcriptomic BodyMap created by performing RNA-Seq on 320 samples from 11 organs of both sexes of juvenile, adolescent, adult and aged Fischer 344 rats. We catalogue the expression profiles of 40,064 genes, 65,167 transcripts, 31,909 alternatively spliced transcript variants and 2,367 non-coding genes/non-coding RNAs (ncRNAs) annotated in AceView. We find that organ-enriched, differentially expressed genes reflect the known organ-specific biological activities. A large number of transcripts show organ-specific, age-dependent or sex-specific differential expression patterns. We create a web-based, open-access rat BodyMap database of expression profiles with crosslinks to other widely used databases, anticipating that it will serve as a primary resource for biomedical research using the rat model.
Stably Expressed Genes (SEGs) whose expression varies within a narrow range may be involved in core cellular processes necessary for basic functions. To identify such genes, we re-analyzed existing ...RNA-Seq gene expression profiles across 11 organs at 4 developmental stages (from immature to old age) in both sexes of F344 rats (n = 4/group; 320 samples). Expression changes (calculated as the maximum expression / minimum expression for each gene) of >19000 genes across organs, ages, and sexes ranged from 2.35 to >109-fold, with a median of 165-fold. The expression of 278 SEGs was found to vary ≤4-fold and these genes were significantly involved in protein catabolism (proteasome and ubiquitination), RNA transport, protein processing, and the spliceosome. Such stability of expression was further validated in human samples where the expression variability of the homologous human SEGs was significantly lower than that of other genes in the human genome. It was also found that the homologous human SEGs were generally less subject to non-synonymous mutation than other genes, as would be expected of stably expressed genes. We also found that knockout of SEG homologs in mouse models was more likely to cause complete preweaning lethality than non-SEG homologs, corroborating the fundamental roles played by SEGs in biological development. Such stably expressed genes and pathways across life-stages suggest that tight control of these processes is important in basic cellular functions and that perturbation by endogenous (e.g., genetics) or exogenous agents (e.g., drugs, environmental factors) may cause serious adverse effects.
The acceptance of microarray technology in regulatory decision-making is being challenged by the existence of various platforms and data analysis methods. A recent report (E. Marshall, Science, 306, ...630-631, 2004), by extensively citing the study of Tan et al. (Nucleic Acids Res., 31, 5676-5684, 2003), portrays a disturbingly negative picture of the cross-platform comparability, and, hence, the reliability of microarray technology.
We reanalyzed Tan's dataset and found that the intra-platform consistency was low, indicating a problem in experimental procedures from which the dataset was generated. Furthermore, by using three gene selection methods (i.e., p-value ranking, fold-change ranking, and Significance Analysis of Microarrays (SAM)) on the same dataset we found that p-value ranking (the method emphasized by Tan et al.) results in much lower cross-platform concordance compared to fold-change ranking or SAM. Therefore, the low cross-platform concordance reported in Tan's study appears to be mainly due to a combination of low intra-platform consistency and a poor choice of data analysis procedures, instead of inherent technical differences among different platforms, as suggested by Tan et al. and Marshall.
Our results illustrate the importance of establishing calibrated RNA samples and reference datasets to objectively assess the performance of different microarray platforms and the proficiency of individual laboratories as well as the merits of various data analysis procedures. Thus, we are progressively coordinating the MAQC project, a community-wide effort for microarray quality control.
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.
Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 ...independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.
We present primary results from the Sequencing Quality Control (SEQC) project, coordinated by the US Food and Drug Administration. Examining Illumina HiSeq, Life Technologies SOLiD and Roche 454 ...platforms at multiple laboratory sites using reference RNA samples with built-in controls, we assess RNA sequencing (RNA-seq) performance for junction discovery and differential expression profiling and compare it to microarray and quantitative PCR (qPCR) data using complementary metrics. At all sequencing depths, we discover unannotated exon-exon junctions, with >80% validated by qPCR. We find that measurements of relative expression are accurate and reproducible across sites and platforms if specific filters are used. In contrast, RNA-seq and microarrays do not provide accurate absolute measurements, and gene-specific biases are observed for all examined platforms, including qPCR. Measurement performance depends on the platform and data analysis pipeline, and variation is large for transcript-level profiling. The complete SEQC data sets, comprising >100 billion reads (10Tb), provide unique resources for evaluating RNA-seq analyses for clinical and regulatory settings.
Serum levels of cardiac troponins serve as biomarkers of myocardial injury. However, troponins are released into the serum only after damage to cardiac tissue has occurred. Here, we report ...development of a mouse model of doxorubicin (DOX)-induced chronic cardiotoxicity to aid in the identification of predictive biomarkers of early events of cardiac tissue injury. Male B6C3F1 mice were administered intravenous DOX at 3mg/kg body weight, or an equivalent volume of saline, once a week for 4, 6, 8, 10, 12, and 14weeks, resulting in cumulative DOX doses of 12, 18, 24, 30, 36, and 42mg/kg, respectively. Mice were sacrificed a week following the last dose. A significant reduction in body weight gain was observed in mice following exposure to a weekly DOX dose for 1week and longer compared to saline-treated controls. DOX treatment also resulted in declines in red blood cell count, hemoglobin level, and hematocrit compared to saline-treated controls after the 2nd weekly dose until the 8th and 9th doses, followed by a modest recovery. All DOX-treated mice had significant elevations in cardiac troponin T concentrations in plasma compared to saline-treated controls, indicating cardiac tissue injury. Also, a dose-related increase in the severity of cardiac lesions was seen in mice exposed to 24mg/kg DOX and higher cumulative doses. Mice treated with cumulative DOX doses of 30mg/kg and higher showed a significant decline in heart rate, suggesting drug-induced cardiac dysfunction. Altogether, these findings demonstrate the development of DOX-induced chronic cardiotoxicity in B6C3F1 mice.
► 24mg/kg was a cumulative cardiotoxic dose of doxorubicin in male B6C3F1 mice. ► Doxorubicin-induced hematological toxicity was in association with splenomegaly. ► Doxorubicin induced severe testicular toxicity in B6C3F1 male mice.