RNA-Seq analysis and gene discovery of Andrias davidianus using Illumina short read sequencing

Molecular cloning, expression analysis, and functional characterization of Intelectin-1 from Chinese giant salamanders (Andrias davidianus)

BACKGROUND

Intelectin (ITLN) plays a pivotal role in innate immunity, inflammatory responses, and tumor progression. However, its physiological roles in caudate amphibians, such as Andrias davidianus, remain elusive.

OBJECTIVE

A proto-type of intelectin (AdITLN1) from A. davidianus was characterized, following which an AdITLN1 homology model was generated to analyze its expression profiles and functional characterizations.

METHODS

The intelectin1 gene (AdITLN1) was cloned, and its evolutionary relationship with ITLN1 from other species was explored. Additionally, a homology model was generated using Molecular Operating Environment (MOE) software, and the active binding pocket of AdITLN1 was identified. Infections with Aeromonas hydrophila were conducted to analyze changes in ITLN1 transcript levels in liver tissue. Finally, recombinant A. davidianus ITLN1 protein (rAdITLN1) was synthesized using prokaryotic expression, and its bacterial agglutination activity and impact on macrophage phagocytosis were examined.

RESULTS

Multiple sequence alignment and phylogenetic analysis indicated that AdITLN1 shared the closest evolutionary relationship with amphibians, with its structure being similar to that of human intelectin1 and Xenopus Embryonic Epidermal Lectin. Moreover, AdITLN1 mRNA was expressed in a wide range of tissues and was significantly up-regulated post-A. hydrophila infection. Meanwhile, the AdITLN1 protein was successfully expressed and purified in Escherichia coli (E. coli) BL21 (DE3). The recombinant AdITLN1 (rAdITLN1) displayed strong agglutination activity against different Gram-negative and Gram-positive bacteria. Lastly, The phagocytosis of rAdITLN1-treated E. coli by macrophages was significantly enhanced.

CONCLUSION

The results of the present study demonstrated that AdITLN1 was a multifunctional immune protein with potent immunomodulatory activity. This study offers valuable insights into disease control in giant salamanders and the conservation of natural resources.

Comparative transcriptomes of nine tissues for the Heilongjiang brown frog (Rana amurensis)

The Heilongjiang brown frog (Rana amurensis) is widely used in traditional Chinese medicine. In particular, the oviduct and skin have been developed into various health products. However, limited numbers of complete genomes of amphibian species have been reported, excluding the Heilongjiang brown frog. Here, the transcriptomes of 45 samples from the liver, spleen, heart, ovaries, thigh muscles, skin, oviduct, stomach and intestine of five Heilongjiang brown frog were reassembled and analyzed. A total of 1,085,532 unigenes with an average length of 676.6 bp and N50 of 722 bp were obtained. Comparative transcriptomics of different tissues detected tissue-specific expression. There were 3248 differentially expressed genes (DEGs) in the ovary, and the number of unique DEGs between the ovary and spleen was the largest. The results of DEGs enrichment showed there were many pathways and items related to protein synthesis and metabolism in the oviduct. The DEGs of the skin were enriched with many bacterial defense items, indicating that there were a large number of antimicrobial peptides in the skin. Thus, these were suitable as biological sources for the development and extraction of antimicrobial peptides. Through the assembly of transcriptome sequencing data and functional annotation of the Heilongjiang brown frog genome, this study provides reference materials for further exploring and utilizing functional gene resources of frogs and lays a foundation for medical research and the development of new products.

Transcriptome characterization and SSR discovery in the giant spiny frog Quasipaa spinosa

The giant spiny frog Quasipaa spinosa (Amphibia: Ranidae) is a large unique frog species found mainly in southern China with a low amount of fat and high protein, and it has become one of the most important aquaculture animal species in China. To better understand its genetic background and screen potential molecular markers for artificial breeding and species conservation, we constructed an expression profile of Q. spinosa with high-throughput RNA sequencing and acquired potential SSR markers. Approximately 81.7 Gb of data and 93,887 unigenes were generated. The transcriptome contains 2085 (80.7 %) complete BUSCOs, suggesting that our assembly methods were effective and accurate.These unigenes were functionally classified using 7 functional databases, yielding 17,482 Pfam-, 12,752 Sting-, 17,526 KEGG-, 24,341 Swiss-Prot-, 28,604 Nr-, 16,287 GO- and 12,752 COG-annotated unigenes. Among several amphibian species, Q. spinosa unigenes had the highest number of hits to Xenopus tropicalis (35.25 %), followed by Xenopus laevis (12.68 %). 1417 unigenes were assigned to the immune system. In addition, a total of 33,019 candidate SSR markers were identified from the constructed library. Further tests with 20 loci and 118 large-scale breeding specimens gathered from four culture farms in China showed that 15 (75 %) loci were polymorphic, with the number of alleles per locus varying from 3 to 9 (mean of 4.3). The PIC values for the SSR markers ranged from 0.19 to 0.82, with an average value of 0.43, indicating moderate polymorphism in Q. spinosa. The transcriptomic profile and SSR repertoire obtained in the present study will facilitate population genetic studies and the selective breeding of amphibian species.

Full-length transcriptome assembly of Andrias davidianus (Amphibia: Caudata) skin via hybrid sequencing

The Chinese giant salamander, Andrias davidianus, is the largest amphibian species in the world; it is thus an economically and ecologically important species. The skin of A. davidianus exhibits complex adaptive structural and functional adaptations to facilitate survival in aquatic and terrestrial ecosystems. Here, we report the first full-length amphibian transcriptome from the dorsal skin of A. davidianus, which was assembled using hybrid sequencing and the PacBio and Illumina platforms. A total of 153,038 transcripts were hybrid assembled (mean length of 2039 bp and N50 of 2172 bp), and 133,794 were annotated in at least one database (nr, Swiss-Prot, KEGG, KOGs, GO, and nt). A total of 58,732, 68,742, and 115,876 transcripts were classified into 24 KOG categories, 1903 GO term categories, and 46 KEGG pathways (level 2), respectively. A total of 207,627 protein-coding regions, 785 transcription factors, 27,237 potential long non-coding RNAs, and 8299 simple sequence repeats were also identified. The hybrid-assembled transcriptome recovered more full-length transcripts, had a higher N50 contig length, and a higher annotation rate of unique genes compared with that assembled in previous studies using next-generation sequencing. The high-quality full-length reference gene set generated in this study will help elucidate the genetic characteristics of A. davidianus skin and aid the identification of functional skin proteins.

Salamanders: The molecular basis of tissue regeneration and its relevance to human disease

Salamanders are recognized for their ability to regenerate a broad range of tissues. They have also have been used for hundreds of years for classical developmental biology studies because of their large accessible embryos. The range of tissues these animals can regenerate is fascinating, from full limbs to parts of the brain or heart, a potential that is missing in humans. Many promising research efforts are working to decipher the molecular blueprints shared across the organisms that naturally have the capacity to regenerate different tissues and organs. Salamanders are an excellent example of a vertebrate that can functionally regenerate a wide range of tissue types. In this review, we outline some of the significant insights that have been made that are aiding in understanding the cellular and molecular mechanisms of tissue regeneration in salamanders and discuss why salamanders are a worthy model in which to study regenerative biology and how this may benefit research fields like regenerative medicine to develop therapies for humans in the future.

Proteomic analysis of eleven tissues in the Chinese giant salamander (Andrias davidianus)

The Chinese giant salamander (Andrias davidianus, CGS) is the largest extant amphibian species in the world. Global quantitative proteome analysis of multiple tissues would indicate tissue-specific physiological processes and clarify the function of each protein from a whole-organism perspective. This study performed proteome analysis of eleven tissues collected from adult CGSs using iTRAQ coupled with LC-MS/MS technology. Based on the predicted protein database from previously obtained CGS transcriptome data, 2153 proteins were identified for subsequent analysis. A weighted gene co-expression network analysis (WGCNA) clustered 2153 proteins into 17 co-expressed modules, which will be useful for predicting the functions of unannotated proteins. The protein levels of molecular complexes with housekeeping functions, such as ribosomes, spliceosomes and mitochondrial respiratory chain complexes, were tightly regulated in different tissues of the CGS, as they are in mammalian tissues. Transcription regulator, pathway and bio-functional analysis of tissue-specific proteins showed that highly expressed proteins largely reflected the physiological functions of specific tissues. Our data, as an initial atlas of protein expression of an amphibian species, will be useful for further molecular biology research on CGS.

RNA-seq analysis provides insight into molecular adaptations of Andrias davidianus

The Chinese giant salamander Andrias davidianus is regarded as an ideal model for studying local adaptations, such as longevity, tolerance to starvation, and cutaneous respiration. Transcriptome analysis is useful for studying the large and complex genomes of amphibians. Based on the coding gene set of adult A. davidianus, dozens of A. davidianus-specific genes were identified and three signaling pathway (JAK-STAT, HIF-1, and FoxO) genes were expanded as compared with other amphibians. The results of the pathway analysis of A. davidianus-specific genes indicated that the molecular adaptation of A. davidianus may have required a more rapid evolution of the immune system. Additionally, for the first time, the gene expressions in different parts of the skin tissue were compared. The results of the comparison analysis demonstrated that lateral skin could be more focused on mucus secretion, dorsal skin on immunity and melanogenesis, and abdominal skin on water and salt metabolism. This study provides the first insight into studying longevity and starvation tolerance in A. davidianus, and offers a basis for further investigation of the molecular mechanisms of adaptations in amphibians.

Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

BACKGROUND

De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or "pipelines," on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research.

RESULTS

New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla.

CONCLUSIONS

Given current bioinformatics approaches, there is no single "best" reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

Identification of critical sex-biased genes in Andrias davidianus by de novo transcriptome

The Chinese giant salamander Andrias davidianus is a protected amphibian with high nutritional and economic value. Understanding its sex determination mechanism is important for improving culture techniques and sex control in breeding. However, little information on the characterization of critical genes involved in sex is available. Herein, sequencing of ovary and test produced 40,783,222 and 46,128,902 raw reads, respectively, which were jointly assembled into 80,497 unigenes. Of these, 36,609 unigenes were annotated, of which 8907 were female-biased and 10,385 were male-biased. Several sex-related pathways were observed, including the Wnt signaling pathway. After elevated temperature and estrogen exposure, neomale and neofemale specimens were identified by a female-specific marker for the first time. RT-qPCR analysis showed the expression profile of ten selected sex-biased genes to be exhibited consistently in male and neomale and in female and neofemale, with the exception of the Amh and TfIIIa genes. Results suggested that these genes may play important roles in A. davidianus sex determination and gonad development. This provides a basis for further investigation of the molecular mechanisms of sex determination in amphibians.

Identification and Characterization of MicroRNAs in Skin of Chinese Giant Salamander (Andrias davidianus) by the Deep Sequencing Approach

MicroRNAs (miRNA) play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. However, the miRNAs from skin of Andrias davidianus have not been reported. In this study, a small-RNA cDNA library was constructed and sequenced from skin of A. davidianus. A total of 513 conserved miRNAs belonging to 174 families were identified. The remaining 108 miRNAs we identified were novel and likely to be skin tissue-specific but were expressed at low levels. The presence of randomly selected 15 miRNAs identified and their expression in eight different tissues from A. davidianus were validated by stem-loop qRT-PCR. For better understanding the functions of miRNAs, 129,791 predicated target genes were analyzed by GO and their pathways illustrated by KEGG pathway analyses. The results show that these identified miRNAs from A. davidianus skin are involved in a broad range of physiological functions including metabolism, growth, development, and immune responses. This study exhaustively identifies miRNAs and their target genes, which will ultimately pave the way for understanding their role in skin of A. davidianus and other amphibians. Further studies are necessary to better understand miRNA-mediated gene regulation.

Comparative expression analysis identifies the respiratory transition-related miRNAs and their target genes in tissues of metamorphosing Chinese giant salamander (Andrias davidianus)

Background

Chinese giant salamander (Andrias davidianus) undergoes a metamorphosis from aquatic larvae to terrestrial adults, with concomitant transfer of respiration from gills to lungs prior to metamorphosis. These two tissues, as well as skin, were sampled to identify the differentially expressed miRNAs.

Results

High-coverage reference transcriptome was generated from combined gill, lung and skin tissues of metamorphosing juveniles, and lung tissue of adults: 86,282 unigenes with total length of approximately 77,275,634 bp and N50 of 1732 bp were obtained. Among these, 13,246 unigenes were assigned to 288 pathways. To determine the possible involvement of miRNAs in the respiratory transition, small RNA libraries were sequenced; 282 miRNAs were identified, 65 among which were known and 217 novel. Based on the hierarchical clustering analysis, the twelve studied samples were classified into three major clusters using differentially expressed miRNAs. We have validated ten differentially expressed miRNAs and some of their related target genes using qPCR. These results largely corroborated the results of transcriptomic and miRNA analyses. Finally, an miRNA-gene-network was constructed. Among them, two miRNAs with target genes related to oxygen sensing were differentially expressed between gill and lung tissues. Three miRNAs were differentially expressed between the lungs of larvae and lungs of adults.

Conclusions

This study provides the first large-scale miRNA expression profile overview during the respiration transition from gills to lungs in Chinese giant salamander. Five differentially expressed miRNAs and their target genes were identified among skin, gill and lung tissues. These results suggest that miRNA profiles in respiratory tissues play an important role in the regulation of respiratory transition.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4662-5) contains supplementary material, which is available to authorized users.

Identification of the first cathelicidin gene from skin of Chinese giant salamanders Andrias davidianus with its potent antimicrobial activity

Cathelicidins, as effector molecules, play important roles against infections and represent a crucial component of the innate immune system in vertebrates. They are widely studied in mammals, but little is known in amphibians. In the present study, we report the identification and characterization of a novel cathelicidin from Chinese giant salamander Andrias davidianus, which is the first study in Caudata amphibian. The cDNA sequence encodes a predicted 148-amino-acid polypeptide, which composed of a 20-residue signal peptide, a 94-residue conserved cathelin domain and a 34-residue mature peptide. From the multiple sequence alignments and phylogenetic analysis, AdCath shared conserved structure with other orthologs and clustered with other amphibian peptides. The tissue expression profiles revealed AdCath was highly expressed in skin. To study the function of AdCath gene, the AdCath precursor protein and mature peptide were recombinantly expressed and chemical synthesized respectively. The rAdCath protein could bind to LPS in a dose-dependent manner. When the concentrations of rAdCath protein and mature peptide were up to 22 μg/mL, they showed significantly cytotoxicity to human 293T cell lines. The rAdCath protein and synthetic peptide could exhibit antibacterial activities detected by the minimum inhibitory concentrations assay. From the SEM assay, the synthetic mature peptide could destroy the membrane of bacteria and cause loss of membrane integrity. Collectively, these findings characterized the first cathelicidin from A. davidianus, and highlighted its potential antimicrobial activities, indicating its important roles in the skin immune response against different bacteria.

Transcriptome analysis of the Chinese giant salamander (Andrias davidianus) using RNA-sequencing

The Chinese giant salamander (Andrias davidianus) is an economically important animal on academic value. However, the genomic information of this species has been less studied. In our study, the transcripts of A. davidianus were obtained by RNA-seq to conduct a transcriptomic analysis. In total 132,912 unigenes were generated with an average length of 690 bp and N50 of 1263 bp by de novo assembly using Trinity software. Using a sequence similarity search against the nine public databases (CDD, KOG, NR, NT, PFAM, Swiss-prot, TrEMBL, GO and KEGG databases), a total of 24,049, 18,406, 36,711, 15,858, 20,500, 27,515, 36,705, 28,879 and 10,958 unigenes were annotated in databases, respectively. Of these, 6323 unigenes were annotated in all database and 39,672 unigenes were annotated in at least one database. Blasted with KEGG pathway, 10,958 unigenes were annotated, and it was divided into 343 categories according to different pathways. In addition, we also identified 29,790 SSRs. This study provided a valuable resource for understanding transcriptomic information of A. davidianus and laid a foundation for further research on functional gene cloning, genomics, genetic diversity analysis and molecular marker exploitation in A. davidianus.

Analysis on the expression and function of a chicken-type and goose-type lysozymes in Chinese giant salamanders Andrias davidianus

Lysozymes as an important immune factor, play vital roles in innate immune response against pathogen infection. In the present study, one c-type and g-type lysozymes were identified from Chinese giant salamander (Andrias davidianus). They shared highly conserved structural features with lysozymes from other species. Spatial expression analysis revealed that AdlysC transcript was most abundant in liver and stomach, and least in muscle and brain. In contrast, the expression level of AdlysG was most abundant in liver and least in muscle and skin. The transcription level of c-type and g-type lysozymes were up-regulated after Aeromonas hydrophila infection in liver and spleen, indicating their participations in the immune response. Moreover, the recombinant AdlysC and AdlysG protein were produced and purified, and were used to investigate the lysozyme activity at different pH and temperatures. The optimal lytic activity was determined at pH 6.0 and at a temperature of 30 °C. Through the minimal inhibitory concentration test, the rAdlysC and rAdlysG exhibited apparent antibacterial activity against both Gram-positive and Gram-negative bacteria with a variable concentration. In conclusion, it is the first report of lysozymes in A. davidianus, and c-type and g-type lysozymes should be involved in the innate immune response of A. davidianus.

A reference gene set construction using RNA-seq of multiple tissues of Chinese giant salamander, Andrias davidianus

Background

Chinese giant salamander (CGS) is the largest extant amphibian species in the world. Owing to its evolutionary position and four peculiar phenomenon of life (longevity, starvation tolerance, regenerative ability, and hatch without sunshine), it is an invaluable model species for research. However, lack of genomic resources leads to fewer study progresses in these fields, due to its huge genome of ∼50 GB making it extremely difficult to be assembled.

Results

We reported the sequenced transcriptome of more than 20 tissues from adult CGS using Illumina Hiseq 2000 technology, and a total of 93 366 no-redundancy transcripts with a mean length of 1326 bp were obtained. We developed for the first time an efficient pipeline to construct a high-quality reference gene set of CGS and obtained 26 135 coding genes. BUSCO and homologous assessment showed that our assembly captured 70.6% of vertebrate universal single-copy orthologs, and this coding gene set had a higher proportion of completeness CDS with comparable quality of the protein sets of Tibetan frog.

Conclusions

These highest quality data will provide a valuable reference gene set to the subsequent research of CGS. In addition, our strategy of de novo transcriptome assembly and protein identification is applicable to similar studies.

Comparative microRNAome analysis of the testis and ovary of the Chinese giant salamander

MicroRNAs (miRNAs) are 18-24 nucleotides non-coding RNAs that regulate gene expression by post-transcriptional suppression of mRNA. The Chinese giant salamander (CGS, Andrias davidianus), which is an endangered species, has become one of the important models of animal evolution; however, no miRNA studies on this species have been conducted. In this study, two small RNA libraries of CGS ovary and testis were constructed using deep sequencing technology. A bioinformatics pipeline was developed to distinguish miRNA sequences from other classes of small RNAs represented in the sequencing data. We found that many miRNAs and other small RNAs such as piRNA and tsRNA were abundant in CGS tissue. A total of 757 and 756 unique miRNAs were annotated as miRNA candidates in the ovary and testis respectively. We identified 145 miRNAs in CGS ovary and 155 miRNAs in CGS testis that were homologous to those in Xenopus laevis ovary and testis respectively. Forty-five miRNAs were more highly expressed in ovary than in testis and 21 miRNAs were more highly expressed in testis than in ovary. The expression profiles of the selected miRNAs (miR-451, miR-10c, miR-101, miR-202, miR-7a and miR-499) had their own different roles in other eight tissues and different development stages of testis and ovary, suggesting that these miRNAs play vital regulatory roles in sexual differentiation, gametogenesis and development in CGS. To our knowledge, this is the first study to reveal miRNA profiles that are related to male and female CGS gonads and provide insights into sex differences in miRNA expression in CGS.

Characterization of galectin-1 from Chinese giant salamanders Andrias davidianus and its involvements during immune response

Galectins are considered as a multifunctional protein which play essential roles in cell adhesion and apoptosis, inflammation, tumor progression and immune response. In spite of extensive studies of galectin importance in immune system among different animals, few studies have been devoted to their functions in amphibian. In the present study, we characterized one proto type of galectin (named AdGal1) from Chinese giant salamander Andrias davidianus and studied its function in immune response. AdGal1 cDNA possesses an open reading frame of 598 bp, which encodes a putative galectin of 134 amino acids containing one carbohydrate recognition domains (CRDs). The constitutive expression of mRNA transcripts was detected in a wide range of tissues, with the highest expression in kidney. Immune challenges with Aeromonas hydrophila and Chinese giant salamander iridovirus (GSIV), the transcript level of AdGal1 in kidney was significantly upregulated. The mature protein of AdGal1 was successfully expressed and purified in Escherichia coli BL21 (DE3). The recombinant AdGal1 (rAdGal1) could show bind activity to different Gram negative and Gram positive bacteria. It could also strongly agglutinate different kinds of bacteria at different concentrations. Collectively, these data from the present study indicate that AdGal1 is a vital pattern recognition receptor to recognize different microbes in the innate immune system of Andrias davidianus.

Sequencing and de novo transcriptome assembly of the Chinese giant salamander (Andrias davidianus)

Next-generation technologies for determination of genomics and transcriptomics composition have a wide range of applications. Andrias davidianus, has become an endangered amphibian species of salamander endemic in China. However, there is a lack of the molecular information. In this study, we obtained the RNA-Seq data from a pool of A. davidianus tissue including spleen, liver, muscle, kidney, skin, testis, gut and heart using Illumina HiSeq 2500 platform. A total of 15,398,997,600 bp were obtained, corresponding to 102,659,984 raw reads. A total of 102,659,984 reads were filtered after removing low-quality reads and trimming the adapter sequences. The Trinity program was used to de novo assemble 132,912 unigenes with an average length of 690 bp and N50 of 1263 bp. Unigenes were annotated through number of databases. These transcriptomic data of A. davidianus should open the door to molecular evolution studies based on the entire transcriptome or targeted genes of interest to sequence. The raw data in this study can be available in NCBI SRA database with accession number of SRP099564.

Proteinaceous Pheromone Homologs Identified from the Cloacal Gland Transcriptome of a Male Axolotl, Ambystoma mexicanum

Pheromones play an important role in modifying vertebrate behavior, especially during courtship and mating. Courtship behavior in urodele amphibians often includes female exposure to secretions from the cloacal gland, as well as other scent glands. The first vertebrate proteinaceous pheromone discovered, the decapeptide sodefrin, is a female attracting pheromone secreted by the cloacal gland of male Cynops pyrrhogaster. Other proteinaceous pheromones in salamanders have been shown to elicit responses from females towards conspecific males. The presence and levels of expression of proteinaceous pheromones have not been identified in the family Ambystomatidae, which includes several important research models. The objective of this research was therefore to identify putative proteinaceous pheromones from male axolotls, Ambystoma mexicanum, as well as their relative expression levels. The results indicate that axolotls possess two different forms of sodefrin precursor-like factor (alpha and beta), as well as a putative ortholog of plethodontid modulating factor. The beta form of sodefrin precursor-like factor was amongst the most highly expressed transcripts within the cloacal gland. The ortholog of plethodontid modulating factor was expressed at a level equivalent to the beta sodefrin precursor-like factor. The results are from a single male axolotl; therefore, we are unable to assess how representative our results may be. Nevertheless, the presence of these highly expressed proteinaceous pheromones suggests that male axolotls use multiple chemical cues to attract female conspecifics. Behavioral assays would indicate whether the putative protein pheromones elicit courtship activity from female axolotls.

Data set for transcriptome analysis of the Chinese giant salamander (Andrias davidianus )

The Chinese giant salamander (Andrias davidianus) occupies a seat at the phylogenetic and species evolution process, which makes it an invaluable model for genetics; however, the genetic information and gene sequences about the Chinese giant salamander in public databases are scanty. Hence, we aimed to perform transcriptome analysis with the help of high-throughput sequencing. In this data, 61,317,940 raw reads were acquired from Chinese giant salamander mRNA using Illumina paired-end sequencing platform. After de novo assembly, a total of 72,072 unigenes were gained, in which 33,834 (46.95%) and 29,479 (40.91%) transcripts exhibited homology to sequences in the Nr database and Swiss-Prot database, (E-value <10(-5)), respectively. In the obtained unigenes, 18,019 (25%) transcripts were assigned with at least one Gene Ontology term, of which 1218 (6.8%) transcripts were assigned to immune system processes. In addition, a total of 17,572 assembled sequences were assigned into 241 predicted KEGG metabolic pathways. Among these, 2552 (14.5%) transcripts were assigned to the immune system relevant pathway and 5 transcripts were identified as potential antimicrobial peptides (AMPs).

Transcriptomic analysis of the host response to an iridovirus infection in Chinese giant salamander, Andrias davidianus

The emergence of an infectious viral disease caused by the Chinese giant salamander iridovirus (GSIV) has led to substantial economic losses. However, no more molecular information is available for the understanding of the mechanisms associated with virus–host interaction. In this study, de novo sequencing was used to obtain abundant high-quality ESTs and investigate differentially-expressed genes in the spleen of Chinese giant salamanders that were either infected or mock infected with GSIV. Comparative expression analysis indicated that 293 genes were down-regulated and 220 genes were up-regulated. Further enrichment analysis showed that the most enriched pathway is “complement and coagulation cascades”, and significantly enriched diseases include “inherited thrombophilia”, “immune system diseases”, “primary immunodeficiency”, “complement regulatory protein defects”, and “disorders of nucleotide excision repair”. Additionally, 30 678 simple sequence repeats (SSRs) from all spleen samples, 26 355 single nucleotide polymorphisms (SNPs) from the spleens of uninfected animals and 36 070 SNPs from the spleens of infected animals were detected. The large amount of variation was specific for the Chinese giant salamanders that were infected with GSIV. The results reported herein provided significant and new EST information that could contribute greatly in investigations into the molecular functions of immune genes in the Chinese giant salamander.

Transcriptome analysis of the endangered Chinese giant salamander (Andrias davidianus): Immune modulation in response to Aeromonas hydrophila infection

The endangered Chinese giant salamander (Andrias davidianus) is the largest extant amphibian species. Disease outbreaks represent one of the major factors threatening A. davidianus populations in the wild and the viability of artificial breeding programmes. Development of future immune therapies to eliminate infectious disease in A. davidianus is dependent on a thorough understanding of the immune mechanisms elicited by pathogen encounters. To this end we have undertaken, for the first time in amphibians, differential transcriptome analysis of the giant salamander response to Aeromonas hydrophila, one of the most devastating pathogens affecting amphibian populations. Out of 87,204 non-redundant consensus unigenes 19,216 were annotated, 6834 of which were upregulated and 906 down-regulated following bacterial infection. 2058 unigenes were involved with immune system processes, including 287 differentially expressed unigenes indicative of the impact of bacterial infection on several innate and adaptive immune pathways in the giant salamander. Other pathways not directly associated with immune-related activity were differentially expressed, including developmental, structural, molecular and growth processes. Overall, this work provides valuable insights into the underlying immune mechanisms elicited during bacterial infection in amphibians that may aid in the future development of disease control measures in protecting the Chinese giant salamander. With the unique position of amphibians in the transition of tetrapods from aquatic to terrestrial habitats, our study will also be invaluable towards the further understanding of the evolution of tetrapod immunity.