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

Perfluorooctanoic Acid (PFOA) Exposure Compromises Fertility by Affecting Ovarian and Oocyte Development

PFOA, a newly emerging persistent organic pollutant, is widely present in various environmental media. Previous reports have proved that PFOA exposure can accumulate in the ovary and lead to reproductive toxicity in pregnant mice. However, the potential mechanism of PFOA exposure on fertility remains unclear. In this study, we explore how PFOA compromises fertility in the zebrafish. The data show that PFOA (100 mg/L for 15 days) exposure significantly impaired fertilization and hatching capability. Based on tissue sections, we found that PFOA exposure led to ovarian damage and a decrease in the percentage of mature oocytes. Moreover, through in vitro incubation, we determined that PFOA inhibits oocyte development. We also sequenced the transcriptome of the ovary of female zebrafish and a total of 284 overlapping DEGs were obtained. Functional enrichment analysis showed that 284 overlapping DEGs function mainly in complement and coagulation cascades signaling pathways. In addition, we identified genes that may be associated with immunity, such as LOC108191474 and ZGC:173837. We found that exposure to PFOA can cause an inflammatory response that can lead to ovarian damage and delayed oocyte development.

Andrias Davidianus Mucus-Based Bioadhesive with Enhanced Adhesion and Wound Healing Properties

The skin secretion of Andrias davidianus (SSAD) is a novel biological adhesive raw material under development. This material exhibits robust adhesion while maintaining the flexibility of the wound. It also has the potential for large-scale production, making it promising for practical application explore. Hence, in-depth research on methods to fine-tune SSAD properties is of great importance to promote its practical applications. Herein, we aim to enhance the adhesive and healing properties of SSAD by incorporating functional components. To achieve this goal, we selected 3,4-dihydroxy-l-phenylalanine and vaccarin as the functional components and mixed them with SSAD, resulting in a new bioadhesive, namely, a formulation termed "enhanced SSAD" (ESSAD). We found that the ESSAD exhibited superior adhesive properties, and its adhesive strength was improved compared with the SSAD. Moreover, ESSAD demonstrated a remarkable ability to promote wound healing. This study presents an SSAD-based bioadhesive formulation with enhanced properties, affirming the feasibility of developing SSAD-based adhesive materials with excellent performance and providing new evidence for the application of SSAD. This study also aims to show that SSAD can be mixed with other substances, and addition of effective components to SSAD can be studied to further adjust or improve its performance.

Isolation and identification of a megalocytivirus strain (SKIV-TJ) from cultured spotted knifejaw (Oplegnathus punctatus) in China and its pathogenicity analysis

The spotted knifejaw (Oplegnathus punctatus) has recently emerged as a highly economically significant farmed fish in China. However, due to increasing environmental pollution and breeding density, a range of infectious diseases, including the iridovirus pathogen, have begun to spread widely. In this study, we isolated and identified a strain of Megalocytivirus, SKIV-TJ, from cultured spotted knifejaw in Tianjin, China. We observed significant cytopathic effects (CPE) in SKIV-TJ-infected spotted knifejaw brain (SKB) cells, and electron microscopy showed numerous virus particles in the cytoplasm of SKB cells 6 days post-infection. The annotated complete genome of SKIV-TJ (GenBank accession number ON075463) contained 112,489 bp and 132 open reading frames. Based on the multigene association evolutionary tree using 26 iridovirus core genes, SKIV-TJ was found to be most closely related to Rock bream iridovirus (RBIV). Cumulative mortality of spotted knifejaw infected with SKIV-TJ reached 100% by day 9. A transcriptomic analysis were conducted and a total of 5517 differentially expressed genes were identified, including 2757 upregulated genes and 2760 downregulated genes. The upregulated genes were associated with viral infection and immune signaling pathways. Our findings provide a valuable genetic resource and a deeper understanding of the immune response to SKIV infection in spotted knifejaw.

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.

Four Mx Genes Identified in Andrias davidianus and Characterization of Their Response to Chinese Giant Salamander Iridovirus Infection

Amphibians, including Andrias davidianus, are declining worldwide partly due to infectious diseases. The Myxovirus resistance (Mx) gene is a typical interferon (IFN)-stimulated gene (ISG) involved in the antiviral immunity. Therefore, knowledge regarding the antiviral immunity of A. davidianus can be used for improved reproduction in captivity and protection in the wild. In this study, we amplified and characterized four different A. davidianus Mx genes (adMx) and generated temporal mRNA expression profiles in healthy and Chinese giant salamander iridovirus (GSIV) infected A. davidianus by qualitative real-time PCR (qPCR). The four adMx genes ranged in length from 2008 to 2840 bp. The sequences revealed conserved protein domains including the dynamin superfamily signature motif and the tripartite guanosine-5-triphosphate (GTP)-binding motif. Gene and deduced amino acid sequence alignment revealed relatively high sequence identity with the Mx genes and proteins of other vertebrates. In phylogenetic analysis, the adMx genes clustered together, but also clustered closely with those of fish species. The four adMx genes were broadly expressed in healthy A. davidianus, but were differentially expressed in the spleen during the GSIV infection. Our results show that the adMx genes share major structural features with their homologs, suggesting similar functions to those in other species.

Bcl-xL Reduces Chinese Giant Salamander Iridovirus-Induced Mitochondrial Apoptosis by Interacting with Bak and Inhibiting the p53 Pathway

Chinese giant salamander iridovirus (GSIV) infection could lead to mitochondrial apoptosis in this animal, a process that involves B-cell lymphoma-2 (BCL-2) superfamily molecules. The mRNA expression level of Bcl-xL, a crucial antiapoptotic molecule in the BCL-2 family, was reduced in early infection and increased in late infection. However, the molecular mechanism remains unknown. In this study, the function and regulatory mechanisms of Chinese giant salamander (Andrias davidianus) Bcl-xL (AdBcl-xL) during GSIV infection were investigated. Western blotting assays revealed that the level of Bcl-xL protein was downregulated markedly as the infection progressed. Plasmids expressing AdBcl-xL or AdBcl-xL short interfering RNAs were separately constructed and transfected into Chinese giant salamander muscle cells. Confocal microscopy showed that overexpressed AdBcl-xL was translocated to the mitochondria after infection with GSIV. Additionally, flow cytometry analysis demonstrated that apoptotic progress was reduced in both AdBcl-xL-overexpressing cells compared with those in the control, while apoptotic progress was enhanced in cells silenced for AdBcl-xL. A lower number of copies of virus major capsid protein genes and a reduced protein synthesis were confirmed in AdBcl-xL-overexpressing cells. Moreover, AdBcl-xL could bind directly to the proapoptotic molecule AdBak with or without GSIV infection. In addition, the p53 level was inhibited and the mRNA expression levels of crucial regulatory molecules in the p53 pathway were regulated in AdBcl-xL-overexpressing cells during GSIV infection. These results suggest that AdBcl-xL plays negative roles in GSIV-induced mitochondrial apoptosis and virus replication by binding to AdBak and inhibiting p53 activation.

The Immune System and the Antiviral Responses in Chinese Giant Salamander, Andrias davidianus

The Chinese giant salamander, belonging to an ancient amphibian lineage, is the largest amphibian existing in the world, and is also an important animal for artificial cultivation in China. However, some aspects of the innate and adaptive immune system of the Chinese giant salamander are still unknown. The Chinese giant salamander iridovirus (GSIV), a member of the Ranavirus genus (family Iridoviridae), is a prominent pathogen causing high mortality and severe economic losses in Chinese giant salamander aquaculture. As a serious threat to amphibians worldwide, the etiology of ranaviruses has been mainly studied in model organisms, such as the Ambystoma tigrinum and Xenopus. Nevertheless, the immunity to ranavirus in Chinese giant salamander is distinct from other amphibians and less known. We review the unique immune system and antiviral responses of the Chinese giant salamander, in order to establish effective management of virus disease in Chinese giant salamander artificial cultivation.

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.

Bid is involved in apoptosis induced by Chinese giant salamander iridovirus and contributes to the viral replication in an amphibian cell line

Bid is a pro-apoptotic BH3-only member of the Bcl-2 superfamily that functions to link the extrinsic apoptotic pathway and the mitochondrial amplification loop of the intrinsic pathway. In this study, the expression and functions of Chinese giant salamander (Andrias davidianus) Bid (AdBid) were investigated. The AdBid cDNA sequence contains an open reading frame (ORF) of 576 nucleotides, encoding a putative protein of 191 aa. AdBid possesses the conserved BH3 interacting domain and shared 34-52% sequence identities with other amphibian Bid. mRNA expression of AdBid was most abundant in muscle. The expression level of AdBid in Chinese giant salamander muscle, kidney and spleen significantly increased after Chinese giant salamander iridovirus (GSIV) infection. Additionally, a plasmid expressing AdBid was constructed and transfected into the Chinese giant salamander muscle cell line (GSM cells). The morphology and cytopathic effect (CPE) and apoptotic process in AdBid over-expressed GSM cells was significantly enhanced during GSIV infection compared with that in control cells. Moreover, a higher level of the virus major capsid protein (MCP) gene copies and protein synthesis was confirmed in the AdBid over-expressed cells. These results indicated that AdBid played a positive role in GSIV induced apoptosis and the viral replication. This study may contribute to the better understanding on the infection mechanism of iridovirus-induced apoptosis.

Co-Expression Network Analysis of Spleen Transcriptome in Rock Bream (Oplegnathus fasciatus) Naturally Infected with Rock Bream Iridovirus (RBIV)

Rock bream iridovirus (RBIV) is a notorious agent that causes high mortality in aquaculture of rock bream (Oplegnathus fasciatus). Despite severity of this virus, no transcriptomic studies on RBIV-infected rock bream that can provide fundamental information on protective mechanism against the virus have been reported so far. This study aimed to investigate physiological mechanisms between host and RBIV through transcriptomic changes in the spleen based on RNA-seq. Depending on infection intensity and sampling time point, fish were divided into five groups: uninfected healthy fish at week 0 as control (0C), heavy infected fish at week 0 (0H), heavy mixed RBIV and bacterial infected fish at week 0 (0MH), uninfected healthy fish at week 3 (3C), and light infected fish at week 3 (3L). We explored clusters from 35,861 genes with Fragments Per Kilo-base of exon per Million mapped fragments (FPKM) values of 0.01 or more through signed co-expression network analysis using WGCNA package. Nine of 22 modules were highly correlated with viral infection (|gene significance (GS) vs. module membership (MM) |> 0.5, p-value < 0.05). Expression patterns in selected modules were divided into two: heavy infected (0H and 0MH) and control and light-infected groups (0C, 3C, and 3L). In functional analysis, genes in two positive modules (5448 unigenes) were enriched in cell cycle, DNA replication, transcription, and translation, and increased glycolysis activity. Seven negative modules (3517 unigenes) built in this study showed significant decreases in the expression of genes in lymphocyte-mediated immune system, antigen presentation, and platelet activation, whereas there was significant increased expression of endogenous apoptosis-related genes. These changes lead to RBIV proliferation and failure of host defense, and suggests the importance of blood cells such as thrombocytes and B cells in rock bream in RBIV infection. Interestingly, a hub gene, pre-mRNA processing factor 19 (PRPF19) showing high connectivity (kME), and expression of this gene using qRT-PCR was increased in rock bream blood cells shortly after RBIV was added. It might be a potential biomarker for diagnosis and vaccine studies in rock bream against RBIV. This transcriptome approach and our findings provide new insight into the understanding of global rock bream-RBIV interactions including immune and pathogenesis mechanisms.

Chinese Giant Salamander (Andrias davidianus) Iridovirus Infection Leads to Apoptotic Cell Death through Mitochondrial Damage, Caspases Activation, and Expression of Apoptotic-Related Genes

Chinese giant salamander iridovirus (GSIV) is the causative pathogen of Chinese giant salamander (Andrias davidianus) iridovirosis, leading to severe infectious disease and huge economic losses. However, the infection mechanism by GSIV is far from clear. In this study, a Chinese giant salamander muscle (GSM) cell line is used to investigate the mechanism of cell death during GSIV infection. Microscopy observation and DNA ladder analysis revealed that DNA fragmentation happens during GSIV infection. Flow cytometry analysis showed that apoptotic cells in GSIV-infected cells were significantly higher than that in control cells. Caspase 8, 9, and 3 were activated in GSIV-infected cells compared with the uninfected cells. Consistently, mitochondria membrane potential (MMP) was significantly reduced, and cytochrome c was released into cytosol during GSIV infection. p53 expression increased at an early stage of GSIV infection and then slightly decreased late in infection. Furthermore, mRNA expression levels of pro-apoptotic genes participating in the extrinsic and intrinsic pathway were significantly up-regulated during GSIV infection, while those of anti-apoptotic genes were restrained in early infection and then rose in late infection. These results collectively indicate that GSIV induces GSM apoptotic cell death involving mitochondrial damage, caspases activation, p53 expression, and pro-apoptotic molecules up-regulation.

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.

Transcriptome analysis of Aeromonas hydrophila infected hybrid sturgeon (Huso dauricus×Acipenser schrenckii)

The hybrid sturgeon (Huso dauricus × Acipenser schrenckii) is an economically important species in China. With the increasing aquaculture of hybrid sturgeon, the bacterial diseases are a great concern of the industry. In this study, de novo sequencing was used to compare the difference in transcriptome in spleen of the infected and mock infected sturgeon with Aeromonas hydrophila. Among 187,244 unigenes obtained, 87,887 unigenes were annotated and 1,147 unigenes were associated with immune responses genes. Comparative expression analysis indicated that 2,723 differently expressed genes between the infected and mock-infected group were identified, including 1,420 up-regulated and 1,303 down-regulated genes. 283 differently expressed anti-bacterial immune related genes were scrutinized, including 168 up-regulated and 115 down-regulated genes. Ten of the differently expressed genes were further validated by qRT-PCR. In this study, toll like receptors (TLRs) pathway, NF-kappa B pathway, class A scavenger receptor pathway, phagocytosis pathway, mannose receptor pathway and complement pathway were shown to be up-regulated in Aeromonas hydrophila infected hybrid sturgeon. Additionally, 65,040 potential SSRs and 2,133,505 candidate SNPs were identified from the hybrid sturgeon spleen transcriptome. This study could provide an insight of host immune genes associated with bacterial infection in hybrid sturgeon.

Rag1 and rag2 gene expressions identify lymphopoietic tissues in juvenile and adult Chinese giant salamander (Andrias davidianus)

Rag1 and rag2 are two closely linked recombination activating genes required for V(D)J recombination of antigen receptors in immature lymphocytes, whose expression can serve as marker to identify the lymphopoietic tissues. To study the development of lymphopoietic tissues in Chinese giant salamander (Andrias davidianus), the Chinese giant salamander rag1 and rag2 coding sequences were cloned and determined. High transcript levels of rag1 and rag2 were co-detected in the thymus before 14 months of age, whereas levels were lower in spleen, liver and kidney at all stage of development. The spatial expression patterns of rag1 and rag2 were studied in combination with igY and tcrβ gene expression using in situ hybridization. Significant transcript signals for rag1, rag2, tcrβ and igY were detected not only in the thymus and spleen but also the liver and kidney of juvenile and adult Chinese giant salamanders, which suggests that cells of lymphocyte lineage are present in multiple tissues of the Chinese giant salamander. This implies that lymphopoiesis may take place in these tissues. The tissue morphology of thymus suggested that the branched thymic primordium developed into mature organ with the development of thymocyte from juvenile to adult. These results not only confirm that as expected the thymus and spleen are primordial lymphopoietic tissues but also suggest that the liver and kidney provide site of lymphocyte differentiation in Chinese giant salamander.

Characterization and function of the T-box 1 gene in Chinese giant salamander Andrias davidianus

We characterized the Andrias davidianus T-box 1 (Tbx1) gene. Tbx1 expression was high in testis and low in other examined tissues. Immunohistochemistry detected tbx1 expression in somatic and germ cells 62 days post-hatching (dph), prior to gonad differentiation. At 210 dph, after gonad differentiation, tbx1 was expressed in spermatogonia and testis somatic cells and in granulosa cells in ovary. Tbx1 expression was up-regulated in ovary after high temperature treatment. In the neomale, tbx1 expression showed a similar profile to normal males, and vice-versa for genetic male. Over-expression of tbx1 in females after injection of TBX1 protein down-regulated the female-biased genes cyp19a and foxl2 and up-regulated the male-biased amh gene. When tbx1 was knocked down by tbx1/siRNA, cyp19a and foxl2 expression was up-regulated, and expression of amh, cyp26a, dmrt1, and wt1 was down-regulated. Results suggest that tbx1 influenced sex-related gene expression and participates in regulation of A. davidianus testis development.

MicroRNA repertoire and comparative analysis of Andrias davidianus infected with ranavirus using deep sequencing

Andrias davidianus is a large and economically important amphibian in China. Ranavirus infection causes serious losses in A. davidianus farming industry. MicroRNA mediated host-pathogen interactions are important in antiviral defense. In this study, five small-RNA libraries from ranavirus infected and non-infected A. davidianus spleens were sequenced using high throughput sequencing. The miRNA expression pattern, potential functions, and target genes were investigated. In total, 1356 known and 431 novel miRNAs were discovered. GO and KEGG analysis revealed that certain miRNA target genes are associated with apoptotic, signal pathway, and immune response categories. Analysis identified 82 downregulated and 9 upregulated differentially expressed miRNAs, whose putative target genes are involved in pattern-recognition receptor signaling pathways and immune response. These findings suggested miRNAs play key roles in A. davidianus's response to ranavirus and could provide a reference for further miRNA functional identification, leading to novel approaches to improve A. davidianus ranavirus resistance.

Divergent transcriptomic responses underlying the ranaviruses-amphibian interaction processes on interspecies infection of Chinese giant salamander

Background

Ranaviruses (family Iridoviridae, nucleocytoplasmic large DNA viruses) have been reported as promiscuous pathogens of cold-blooded vertebrates. Rana grylio virus (RGV, a ranavirus), from diseased frog Rana grylio with a genome of 105.79 kb and Andrias davidianus ranavirus (ADRV), from diseased Chinese giant salamander (CGS) with a genome of 106.73 kb, contains 99% homologous genes.

Results

To uncover the differences in virus replication and host responses under interspecies infection, we analyzed transcriptomes of CGS challenged with RGV and ADRV in different time points (1d, 7d) for the first time. A total of 128,533 unigenes were obtained from 820,858,128 clean reads. Transcriptome analysis revealed stronger gene expression of RGV than ADRV at 1 d post infection (dpi), which was supported by infection in vitro. RGV replicated faster and had higher titers than ADRV in cultured CGS cell line. RT-qPCR revealed the RGV genes including the immediate early gene (RGV-89R) had higher expression level than that of ADRV at 1 dpi. It further verified the acute infection of RGV in interspecies infection.

The number of differentially expressed genes and enriched pathways from RGV were lower than that from ADRV, which reflected the variant host responses at transcriptional level. No obvious changes of key components in pathway “Antigen processing and presentation” were detected for RGV at 1 dpi. Contrarily, ADRV infection down-regulated the expression levels of MHC I and CD8. The divergent host immune responses revealed the differences between interspecies and natural infection, which may resulted in different fates of the two viruses. Altogether, these results revealed the differences in transcriptome responses among ranavirus interspecies infection of amphibian and new insights in DNA virus-host interactions in interspecies infection.

Conclusion

The DNA virus (RGV) not only expressed self-genes and replicated quickly after entry into host under interspecies infection, but also avoided the over-activation of host responses. The strategy could gain time for the survival of interspecies pathogen, and may provide opportunity for its adaptive evolution and interspecies transmission.

Electronic supplementary material

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

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.

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.

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.

Differentially expressed immune-related genes in hemocytes of the pearl oyster Pinctada fucata against allograft identified by transcriptome analysis

The pearl oyster Pinctada fucata is commonly cultured for marine pearls in China. To culture pearls, a mantle piece from a donor pearl oyster is grafted with a nucleus into a receptor. This transplanted mantle piece may be rejected by the immune system of the recipient oyster, thus reducing the success of transplantation. However, there have been limited studies about the oyster's immune defense against allograft. In this study, hemocyte transcriptome analysis was performed to detect the immune responses to allograft in P. fucata at 0 h and 48 h after a transplant. The sequencing reaction produced 92.5 million reads that were mapped against the reference genome sequences of P. fucata. The Gene Ontology (GO) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to identify all immune-related differentially expressed genes (DEGs). Compared with patterns at 0 h, a total of 798 DEGs were identified, including 410 up-regulated and 388 down-regulated genes at 48 h. The expression levels of interleukin receptor and toll-like receptor in hemocytes were increased significantly 48 h post-transplant, indicating that the oyster immune response was induced. Finally, altered levels of 18 randomly selected immune-related DEGs were confirmed by quantitative real-time PCR (qRT-PCR). Our results provide the basis for further analysis of the immune rejection of allotransplantation.

Transcriptomic analysis of the head kidney of Topmouth culter (Culter alburnus) infected with Flavobacterium columnare with an emphasis on phagosome pathway

Flavobacterium columnare (FC) has caused worldwide fish columnaris disease with high mortality and great economic losses in cultured fish, including Topmouth culter (Culter alburnus). However, the knowledge about the host factors involved in FC infection is little known. In this study, the transcriptomic profiles of the head kidney from Topmouth culter with or without FC infection were obtained using HiSeq™ 2500 (Illumina). Totally 79,641 unigenes with high quality were obtained. Among them, 4037 differently expressed genes, including 1217 up-regulated and 2820 down-regulated genes, were identified and enriched using databases of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The differently expressed genes were mainly associated with pathways such as immune response, carbohydrate metabolism, amino acid metabolism, and lipid metabolism. Since phagocytosis is a central mechanism of innate immune response by host cells to defense against infectious agents, genes related to the phagosome pathway were scrutinized and 9 differently expressed phagosome-related genes were identified including 3 up-regulated and 6 down-regulated genes. Five of them were further validated by quantitative real-time polymerase chain reaction (qRT-PCR). This transcriptomic analysis of host genes in response to FC infection provides data towards understanding the infection mechanisms and will shed a new light on the prevention of columnaris.

De novo annotation of the immune-enriched transcriptome provides insights into immune system genes of Chinese sturgeon (Acipenser sinensis)

Chinese sturgeon (Acipenser sinensis), one of the oldest extant actinopterygian fishes with very high evolutionary, economical and conservation interest, is considered to be one of the critically endangered aquatic animals in China. Up to date, the immune system of this species remains largely undetermined with little sequence information publicly available. Herein, the first comprehensive transcriptome of immune tissues for Chinese sturgeon was characterized using Illumina deep sequencing. Over 67 million high-quality reads were generated and de novo assembled into the final set of 91,739 unique sequences. The annotation pipeline revealed that 25,871 unigenes were successfully annotated in the public databases, of which only 2002 had significant match to the existing sequences for the genus Acipenser. Overall 22,827 unigenes were categorized into 52 GO terms, 12,742 were classified into 26 KOG categories, and 4968 were assigned to 339 KEGG pathways. A more detailed annotation search showed the presence of a notable representation of immune-related genes, which suggests that this non-teleost actinopterygian fish harbors the same intermediates as in the well known immune pathways from mammals and teleosts, such as pattern recognition receptor (PRR) signaling pathway, JAK-STAT signaling pathway, complement and coagulation pathway, T-cell receptor (TCR) and B-cell receptor (BCR) signaling pathways. Additional genetic marker discovery led to the retrieval of 20,056 simple sequence repeats (SSRs) and 327,140 single nucleotide polymorphisms (SNPs). This immune-enriched transcriptome of Chinese sturgeon represents a rich resource that adds to the currently nascent field of chondrostean fish immunogenetics and furthers the conservation and management of this valuable fish.