The fish gill, as one of the mucosal barriers, plays an important role in mucosal immune response. The fish swimbladder functions for regulating buoyancy. The fish swimbladder has long been postulated as a homologous organ of the tetrapod lung, but the molecular evidence is scarce. In order to provide new information that is complementary to gill immune genes, initiate new research directions concerning the genetic basis of the gill immune response and understand the molecular function of swimbladder as well as its relationship with lungs, transcriptome analysis of the fugu Takifugu rubripes gill and swimbladder was carried out by RNA-Seq. Approximately 55,061,524 and 44,736,850 raw sequence reads from gill and swimbladder were generated, respectively. Gene ontology (GO) and KEGG pathway analysis revealed diverse biological functions and processes. Transcriptome comparison between gill and swimbladder resulted in 3,790 differentially expressed genes, of which 1,520 were up-regulated in the swimbladder while 2,270 were down-regulated. In addition, 406 up regulated isoforms and 296 down regulated isoforms were observed in swimbladder in comparison to gill. By the gene enrichment analysis, the three immune-related pathways and 32 immune-related genes in gill were identified. In swimbladder, five pathways including 43 swimbladder-enriched genes were identified. This work should set the foundation for studying immune-related genes for the mucosal immunity and provide genomic resources to study the relatedness of the fish swimbladder and mammalian lung.

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies in many species. High-throughput sequencing of RNA was developed primarily to analyze global gene expression, while it is an efficient way to discover SNPs from the expressed genes. In this study, we conducted transcriptome sequencing of the swimbladder of Takifugu rubripes using Illumina HiSeq2000 platform to identify gene-associated SNPs in the swimbladder. A total of 30,312,181 unique-mapped-reads were obtained from 44,736,850 raw reads. A total of 62,270 putative SNPs were discovered, which were located in 11,306 expressed genes and 2,246 scaffolds. The average minor allele frequency (MAF) of the SNPs was 0.26. GO and KEGG pathway analysis were conducted to analyze the genes containing SNPs. Validation of selected SNPs revealed that 54% of SNPs (26/48) were true SNPs. The results suggest that RNA-Seq is an efficient and cost-effective approach to discover gene-associated SNPs. In this study, a large number of SNPs were identified and these data will be useful resources for population genetic study, evolution analysis, resource assessment, genetic linkage analysis and genome-wide association studies.

Papilla and skin are two important organs of the sea cucumber. Both tissues have ectodermic origin, but they are morphologically and functionally very different. In the present study, we performed comparative transcriptome analysis of the papilla and skin from the sea cucumber (Apostichopus japonicus) in order to identify and characterize gene expression profiles by using RNA-Seq technology. We generated 30.6 and 36.4 million clean reads from the papilla and skin and de novo assembled in 156,501 transcripts. The Gene Ontology (GO) analysis indicated that cell part, metabolic process and catalytic activity were the most abundant GO category in cell component, biological process and molecular funcation, respectively. Comparative transcriptome analysis between the papilla and skin allowed the identification of 1,059 differentially expressed genes, of which 739 genes were expressed at higher levels in papilla, while 320 were expressed at higher levels in skin. In addition, 236 differentially expressed unigenes were not annotated with any database, 160 of which were apparently expressed at higher levels in papilla, 76 were expressed at higher levels in skin. We identified a total of 288 papilla-specific genes, 171 skin-specific genes and 600 co-expressed genes. Also, 40 genes in papilla-specific were not annotated with any database, 2 in skin-specific. Development-related genes were also enriched, such asfibroblast growth factor,transforming growth factor-β,collagen-α2andIntegrin-α2, which may be related to the formation of the papilla and skin in sea cucumber. Further pathway analysis identified ten KEGG pathways that were differently enriched between the papilla and skin. The findings on expression profiles between two key organs of the sea cucumber should be valuable to reveal molecular mechanisms involved in the development of organs that are related but with morphological differences in the sea cucumber.