De novo Transcriptome Assemblies of Rana (Lithobates) catesbeiana and Xenopus laevis Tadpole Livers for Comparative Genomics without Reference Genomes

PLoS One. 2015 Jun 29;10(6):e0130720. doi: 10.1371/journal.pone.0130720. eCollection 2015.

Abstract

In this work we studied the liver transcriptomes of two frog species, the American bullfrog (Rana (Lithobates) catesbeiana) and the African clawed frog (Xenopus laevis). We used high throughput RNA sequencing (RNA-seq) data to assemble and annotate these transcriptomes, and compared how their baseline expression profiles change when tadpoles of the two species are exposed to thyroid hormone. We generated more than 1.5 billion RNA-seq reads in total for the two species under two conditions as treatment/control pairs. We de novo assembled these reads using Trans-ABySS to reconstruct reference transcriptomes, obtaining over 350,000 and 130,000 putative transcripts for R. catesbeiana and X. laevis, respectively. Using available genomics resources for X. laevis, we annotated over 97% of our X. laevis transcriptome contigs, demonstrating the utility and efficacy of our methodology. Leveraging this validated analysis pipeline, we also annotated the assembled R. catesbeiana transcriptome. We used the expression profiles of the annotated genes of the two species to examine the similarities and differences between the tadpole liver transcriptomes. We also compared the gene ontology terms of expressed genes to measure how the animals react to a challenge by thyroid hormone. Our study reports three main conclusions. First, de novo assembly of RNA-seq data is a powerful method for annotating and establishing transcriptomes of non-model organisms. Second, the liver transcriptomes of the two frog species, R. catesbeiana and X. laevis, show many common features, and the distribution of their gene ontology profiles are statistically indistinguishable. Third, although they broadly respond the same way to the presence of thyroid hormone in their environment, their receptor/signal transduction pathways display marked differences.

Publication types

  • Comparative Study
  • Research Support, Non-U.S. Gov't

MeSH terms

  • Animals
  • Gene Expression Profiling
  • Gene Expression Regulation, Developmental
  • Gene Ontology
  • Genome*
  • Genomics*
  • High-Throughput Nucleotide Sequencing
  • Larva / genetics
  • Liver / metabolism*
  • Molecular Sequence Annotation
  • RNA, Messenger / genetics
  • RNA, Messenger / metabolism
  • Rana catesbeiana / genetics*
  • Reference Standards
  • Signal Transduction / genetics
  • Transcriptome / genetics*
  • Xenopus laevis / genetics*

Substances

  • RNA, Messenger

Grants and funding

This work was funded by the following: Genome British Columbia (IB and CCH), Mitacs Accelerate Internship (SAH), and Natural Sciences and Engineering Research Council of Canada (IB and CCH). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.