Click here to close
Hello! We notice that you are using Internet Explorer, which is not supported by Xenbase and may cause the site to display incorrectly.
We suggest using a current version of Chrome,
FireFox, or Safari.
PLoS One
2015 Jun 18;106:e0130720. doi: 10.1371/journal.pone.0130720.
Show Gene links
Show Anatomy links
De novo Transcriptome Assemblies of Rana (Lithobates) catesbeiana and Xenopus laevis Tadpole Livers for Comparative Genomics without Reference Genomes.
Birol I
,
Behsaz B
,
Hammond SA
,
Kucuk E
,
Veldhoen N
,
Helbing CC
.
???displayArticle.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 tadpoleliver 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.
???displayArticle.pubmedLink???
26121473
???displayArticle.pmcLink???PMC4488148 ???displayArticle.link???PLoS One
Fig 2. Differential expression of assembled transcripts for (A) R. catesbeiana and (B) X. laevis.Three shades of purple designate p-values of differential expression estimates 0.05, 0.02 and 0.002, lighter colours indicating lower thresholds.
Fig 3. Gene ontology classification of DESeq-selected, TH-responsive R. catesbeiana and X. laevis liver transcripts with UniProtKB AC numbers.The two series in each stacked bar plot correspond to differentially expressed R. catesbeiana (CAT) and X. laevis (LAE) transcripts, with a p-value threshold of 5%.
Fig 4. Pathway analysis for liver transcripts from R. catesbeiana (CAT) and X. laevis (LAE).Top 25 impacted pathways after TH treatment for R. catesbeiana ranked by the highest proportion of overall observed genes. The pathway names are indicated in the center of the figure with the total number of genes known in each IGA pathway indicated. The asterisk indicates those pathways that are found in the top 25 list of X. laevis. The colour coded bar plots illustrate the percentage of the total number of gene transcripts in a pathway that are downregulated (blue), non-responsive (yellow), upregulated (red) or not observed in the experiment (gray) relative to the control condition. Differentially expressed transcripts were determined using a p-value threshold of 5%.
Fig 5. Pathway analysis for liver transcripts from R. catesbeiana (CAT) and X. laevis (LAE).Top 25 impacted pathways after TH treatment for X. laevis ranked by the highest proportion of overall observed genes. The asterisk indicates those pathways that are found in the top 25 list of R. catesbeiana. Plot details are as in the Fig 4 legend.
Fig 1. Gene ontology classification of all reconstructed R. catesbeiana and X. laevis liver transcripts with UniProtKB AC numbers.The two series in each stacked bar plot correspond to all R. catesbeiana (CAT) and X. laevis (LAE) transcripts.
Altschul,
Basic local alignment search tool.
1990, Pubmed
Altschul,
Basic local alignment search tool.
1990,
Pubmed
Anders,
Differential expression analysis for sequence count data.
2010,
Pubmed
Bazinet,
Can RNA-Seq resolve the rapid radiation of advanced moths and butterflies (Hexapoda: Lepidoptera: Apoditrysia)? An exploratory study.
2013,
Pubmed
Birol,
De novo transcriptome assembly with ABySS.
2009,
Pubmed
Birol,
Kleat: cleavage site analysis of transcriptomes.
2015,
Pubmed
Blower,
Combining different mRNA capture methods to analyze the transcriptome: analysis of the Xenopus laevis transcriptome.
2013,
Pubmed
,
Xenbase
Camacho,
BLAST+: architecture and applications.
2009,
Pubmed
Cheng,
Differential transcription of fathead minnow immune-related genes following infection with frog virus 3, an emerging pathogen of ectothermic vertebrates.
2014,
Pubmed
Chu,
BioBloom tools: fast, accurate and memory-efficient host species sequence screening using bloom filters.
2014,
Pubmed
Das,
Multiple thyroid hormone-induced muscle growth and death programs during metamorphosis in Xenopus laevis.
2002,
Pubmed
,
Xenbase
Daszak,
Emerging infectious diseases and amphibian population declines.
1999,
Pubmed
Denver,
Thyroid hormone-dependent gene expression program for Xenopus neural development.
1997,
Pubmed
,
Xenbase
Duarte-Guterman,
Regulation of thyroid hormone-, oestrogen- and androgen-related genes by triiodothyronine in the brain of Silurana tropicalis.
2010,
Pubmed
,
Xenbase
Duarte-Guterman,
Transcript profiles and triiodothyronine regulation of sex steroid- and thyroid hormone-related genes in the gonad-mesonephros complex of Silurana tropicalis.
2011,
Pubmed
,
Xenbase
Duarte-Guterman,
Expression and T3 regulation of thyroid hormone- and sex steroid-related genes during Silurana (Xenopus) tropicalis early development.
2010,
Pubmed
,
Xenbase
Eggert,
Sex determination: the amphibian models.
2004,
Pubmed
Fisher,
Emerging fungal threats to animal, plant and ecosystem health.
2012,
Pubmed
Flicek,
Ensembl 2008.
2008,
Pubmed
Frischkorn,
De novo assembly of Aureococcus anophagefferens transcriptomes reveals diverse responses to the low nutrient and low light conditions present during blooms.
2014,
Pubmed
Genome 10K Community of Scientists,
Genome 10K: a proposal to obtain whole-genome sequence for 10,000 vertebrate species.
2009,
Pubmed
Grabherr,
Full-length transcriptome assembly from RNA-Seq data without a reference genome.
2011,
Pubmed
Grayfer,
Immune evasion strategies of ranaviruses and innate immune responses to these emerging pathogens.
2012,
Pubmed
,
Xenbase
Guttman,
Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals.
2009,
Pubmed
Guttman,
Ab initio reconstruction of cell type-specific transcriptomes in mouse reveals the conserved multi-exonic structure of lincRNAs.
2010,
Pubmed
Guttman,
lincRNAs act in the circuitry controlling pluripotency and differentiation.
2011,
Pubmed
Helbing,
The metamorphosis of amphibian toxicogenomics.
2012,
Pubmed
,
Xenbase
Hellsten,
The genome of the Western clawed frog Xenopus tropicalis.
2010,
Pubmed
,
Xenbase
Ichu,
Metabolomic insights into system-wide coordination of vertebrate metamorphosis.
2014,
Pubmed
James-Zorn,
Xenbase: expansion and updates of the Xenopus model organism database.
2013,
Pubmed
,
Xenbase
Kiemnec-Tyburczy,
Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication.
2012,
Pubmed
,
Xenbase
Langlois,
Fadrozole and finasteride exposures modulate sex steroid- and thyroid hormone-related gene expression in Silurana (Xenopus) tropicalis early larval development.
2010,
Pubmed
,
Xenbase
Langlois,
Expression profiles of reproduction- and thyroid hormone-related transcripts in the brains of chemically-induced intersex frogs.
2011,
Pubmed
,
Xenbase
Langmead,
Fast gapped-read alignment with Bowtie 2.
2012,
Pubmed
Li,
Fast and accurate short read alignment with Burrows-Wheeler transform.
2009,
Pubmed
McMahon,
Amphibians acquire resistance to live and dead fungus overcoming fungal immunosuppression.
2014,
Pubmed
McMenamin,
Thyroid hormone-dependent adult pigment cell lineage and pattern in zebrafish.
2014,
Pubmed
Mora,
How many species are there on Earth and in the ocean?
2011,
Pubmed
Mukhi,
Gene switching at Xenopus laevis metamorphosis.
2010,
Pubmed
,
Xenbase
Nakasugi,
Combining transcriptome assemblies from multiple de novo assemblers in the allo-tetraploid plant Nicotiana benthamiana.
2014,
Pubmed
Nordberg,
The genome portal of the Department of Energy Joint Genome Institute: 2014 updates.
2014,
Pubmed
Pagani,
The Genomes OnLine Database (GOLD) v.4: status of genomic and metagenomic projects and their associated metadata.
2012,
Pubmed
Parra,
Assessing the gene space in draft genomes.
2009,
Pubmed
Pooyaei Mehr,
Transcriptome deep-sequencing and clustering of expressed isoforms from Favia corals.
2013,
Pubmed
Pruitt,
NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins.
2007,
Pubmed
Relyea,
The toxicity of Roundup Original Max to 13 species of larval amphibians.
2009,
Pubmed
Ribeiro,
An insight into the sialome of the frog biting fly, Corethrella appendiculata.
2014,
Pubmed
Robert,
Comparative and developmental study of the immune system in Xenopus.
2009,
Pubmed
,
Xenbase
Robertson,
De novo assembly and analysis of RNA-seq data.
2010,
Pubmed
Rollins-Smith,
Metamorphosis and the amphibian immune system.
1998,
Pubmed
Rosenblum,
Genome-wide transcriptional response of Silurana (Xenopus) tropicalis to infection with the deadly chytrid fungus.
2009,
Pubmed
,
Xenbase
Savage,
Conservation and divergence in the frog immunome: pyrosequencing and de novo assembly of immune tissue transcriptomes.
2014,
Pubmed
Schulz,
Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels.
2012,
Pubmed
Simpson,
Efficient de novo assembly of large genomes using compressed data structures.
2012,
Pubmed
Sookruksawong,
RNA-Seq analysis reveals genes associated with resistance to Taura syndrome virus (TSV) in the Pacific white shrimp Litopenaeus vannamei.
2013,
Pubmed
Sumida,
Sequencing and analysis of the internal transcribed spacers (ITSs) and coding regions in the EcoR I fragment of the ribosomal DNA of the Japanese pond frog Rana nigromaculata.
2004,
Pubmed
,
Xenbase
Trapnell,
Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation.
2010,
Pubmed
Ulitsky,
Conserved function of lincRNAs in vertebrate embryonic development despite rapid sequence evolution.
2011,
Pubmed
Veldhoen,
Enabling comparative gene expression studies of thyroid hormone action through the development of a flexible real-time quantitative PCR assay for use across multiple anuran indicator and sentinel species.
2014,
Pubmed
,
Xenbase
Xu,
RNA CoMPASS: a dual approach for pathogen and host transcriptome analysis of RNA-seq datasets.
2014,
Pubmed
