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

Mining expressed sequence tag (EST) microsatellite markers to assess the genetic differentiation of five Hynobius species endemic to Taiwan

Taiwan harbors five endemic species of salamanders (Hynobius spp.) that inhabit distinct alpine regions, contributing to population fragmentation across isolated "sky islands". With an evolutionary history spanning multiple glacial-interglacial cycles, these species represent an exceptional paradigm for exploring biogeography and speciation. However, a lack of suitable genetic markers applicable across species has limited research efforts. Thus, developing cross-amplifying markers is imperative. Expressed sequence-tag simple-sequence repeats (EST-SSRs) that amplify across divergent lineages are ideal for species identification in instances where phenotypic differentiation is challenging. Here, we report a suite of cross-amplifying EST-SSRs from the transcriptomes of the five Hynobius species that exhibit an interspecies transferability rate of 67.67%. To identify individual markers exhibiting cross-species polymorphism and to assess interspecies genetic diversity, we assayed 140 individuals from the five species across 84 sampling sites. A set of EST-SSRs with a high interspecies polymorphic information content (PIC = 0.63) effectively classified these individuals into five distinct clusters, as supported by discriminant analysis of principal components (DAPC), STRUCTURE assignment tests, and Neighbor-joining trees. Moreover, pair-wise FST values > 0.15 indicate notable between-cluster genetic divergence. Our set of 20 polymorphic EST-SSRs is suitable for assessing population structure within and among Hynobius species, as well as for long-term monitoring of their genetic composition.

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.

Re-building limbs, one cell at a time

New techniques for visualizing and interrogating single cells hold the key to unlocking the underlying mechanisms of salamander limb regeneration.

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.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

Transcriptomic Features of Echinococcus granulosus Protoscolex during the Encystation Process

Cystic echinococcosis (CE) is a zoonotic infection caused by Echinococcus granulosus larvae. It seriously affects the development of animal husbandry and endangers human health. Due to a poor understanding of the cystic fluid formation pathway, there is currently a lack of innovative methods for the prevention and treatment of CE. In this study, the protoscoleces (PSCs) in the encystation process were analyzed by high-throughput RNA sequencing. A total of 32,401 transcripts and 14,903 cDNAs revealed numbers of new genes and transcripts, stage-specific genes, and differently expressed genes. Genes encoding proteins involved in signaling pathways, such as putative G-protein coupled receptor, tyrosine kinases, and serine/threonine protein kinase, were predominantly up-regulated during the encystation process. Antioxidant enzymes included cytochrome c oxidase, thioredoxin glutathione, and glutathione peroxidase were a high expression level. Intriguingly, KEGG enrichment suggested that differentially up-regulated genes involved in the vasopressin-regulated water reabsorption metabolic pathway may play important roles in the transport of proteins, carbohydrates, and other substances. These results provide valuable information on the mechanism of cystic fluid production during the encystation process, and provide a basis for further studies on the molecular mechanisms of growth and development of PSCs.

Differential Proteomic Analysis of Chinese Giant Salamander Liver in Response to Fasting

Chinese giant salamander Andrias davidianus has strong tolerance to starvation. Fasting triggers a complex array of adaptive metabolic responses, a process in which the liver plays a central role. Here, a high-throughput proteomic analysis was carried out on liver samples obtained from adult A. davidianus after 3, 7, and 11 months of fasting. As a result, the expression levels of 364 proteins were significantly changed in the fasted liver. Functional analysis demonstrated that the expression levels of key proteins involved in fatty acid oxidation, tricarboxylic acid cycle, gluconeogenesis, ketogenesis, amino acid oxidation, urea cycle, and antioxidant systems were increased in the fasted liver, especially at 7 and 11 months after fasting. In contrast, the expression levels of vital proteins involved in pentose phosphate pathway and protein synthesis were decreased after fasting. We also found that fasting not only activated fatty acid oxidation and ketogenesis-related transcription factors PPARA and PPARGC1A, but also activated gluconeogenesis-related transcription factors FOXO1, HNF4A, and KLF15. This study confirms the central role of lipid and acetyl-CoA metabolism in A. davidianus liver in response to fasting at the protein level and provides insights into the molecular mechanisms underlying the metabolic response of A. davidianus liver to fasting.

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.