Conservation and divergence in the frog immunome: pyrosequencing and de novo assembly of immune tissue transcriptomes

By March 12, 2015 Amphibians, RNAi, Workgroup 2

Author(s): Anna E. Savage, Karen M. Kiemnec-Tyburczy, Amy R. Ellison, Robert C. Fleischer, Kelly R. Zamudio

Publication: Gene

Publication Date: 2014

Background: Frogs are a diverse group of vertebrates forwhich limited genomic resources are available. Natural frog populations face amultitude of threats, including habitat degradation, infectious disease, and environmental change. Characterizing the functional genomics of anuran tissues in general – and the immune system in particular – will enhance our knowledge of genetic and epigenetic responses to environmental threats and inform conservation and recovery efforts.
Results: To increase the number of species with genomic datasets and characterize gene expression in immunerelated tissues, we sequenced the transcriptomes of three tissues from two frogs (Espadarana prosoblepon and Lithobates yavapaiensis) on the Roche 454 GS FLX platform. Our sequencing produced 8881 E. prosoblepon and 5428 L. yavapaiensis annotated gene products after de novo assembly and Gene Ontology classification. Transcripts of the innate and acquired immune system were expressed in all three tissues. Inflammatory response and acquired immunity transcripts were significantly more diverged between E. prosoblepon and L. yavapaiensis compared to innate immunity and immune system development transcripts. Immune-related transcripts did not show an overall elevated rate of functional evolution, with the exception of glycosyl proteases, which include lysozymes, central bacterial and fungal-killing enzymes of the innate immune system.
Conclusions: The three frog transcriptomes provide more than 600 Mbp of new genomic data, and will serve as a valuable framework for future comparative studies of non-model anurans. Additionally, we show that immune gene divergence varies by functional group and that transcriptome studies can be useful in comparing rates of evolutionary change across gene families.

View and Download