Characterization of the global transcriptome and microsatellite marker information for spotted halibut Verasper variegatus

De Novo Transcriptomic Characterization Enables Novel Microsatellite Identification and Marker Development in Betta splendens

The wild populations of the commercially valuable ornamental fish species, Betta splendens, and its germplasm resources have long been threatened by habitat degradation and contamination with artificially bred fish. Because of the lack of effective marker resources, population genetics research projects are severely hampered. To generate genetic data for developing polymorphic simple sequence repeat (SSR) markers and identifying functional genes, transcriptomic analysis was performed. Illumina paired-end sequencing yielded 105,505,486 clean reads, which were then de novo assembled into 69,836 unigenes. Of these, 35,751 were annotated in the non-redundant, EuKaryotic Orthologous Group, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology databases. A total of 12,751 SSR loci were identified from the transcripts and 7970 primer pairs were designed. One hundred primer pairs were randomly selected for PCR validation and 53 successfully generated target amplification products. Further validation demonstrated that 36% (n = 19) of the 53 amplified loci were polymorphic. These data could not only enrich the genetic information for the identification of functional genes but also effectively facilitate the development of SSR markers. Such knowledge would accelerate further studies on the genetic variation and evolution, comparative genomics, linkage mapping and molecular breeding in B. splendens.

Genome assembly and annotation at the chromosomal level of first Pleuronectidae: Verasper variegatus provides a basis for phylogenetic study of Pleuronectiformes

High quality genome is of great significance for the mining of biological information resources of species. Up to now, the genomic information of several important economic flatfishes has been well explained. All these fishes are eyes on left side-type, and no high-quality genome of eyes on right side-type species has been reported. In this study, we applied a combined strategy involving stLFR and Hi-C technologies to generate sequencing data for constructing the chromosomal genome of Verasper variegates, which belongs to Pleuronectidae with characteristic of eyes on right side. The size of genome of V. variegatus is 556 Mb. More than 97.2% of BUSCO genes were detected, and N50 lengths of the contigs and scaffolds reached 79.8 Kb and 23.8 Mb, respectively, demonstrating the outstanding completeness and sequence continuity of the genome. A total of 22,199 protein-coding genes were predicted in the assembled genome, and more than 95% of those genes could be functionally annotated. Meanwhile, the genomic collinearity, gene family and phylogenetic analyses of similar species in Pleuronectiformes were also investigated and portrayed for metamorphosis and benthic adaptation. Sex related genes mapping has also been achieved at the chromosome level. This study is the first chromosomal level genome of a Pleuronectidae fish (V. variegatus). The chromosomal genome assembly constructed in this work will not only be valuable for conservation and aquaculture studies of the V. variegatus but will also be of general interest in the phylogenetic and taxonomic studies of Pleuronectiformes.

Transcriptome analysis of the spleen provides insight into the immunoregulation of Mastacembelus armatus under Aeromonas veronii infection

Mastacembelus armatus, also known as the zigzag eel, is an economically important species of freshwater fish that is very popular with consumers as a high-grade table fish in China. Recently, the wild population of this fish has declined gradually due to overfishing and various types of ecological imbalance. Meanwhile, the aquaculture of this spiny eel has flourished in southern China. To understand the immune response of zigzag eel to Aeromonas veronii, we carried out transcriptome sequencing of zigzag eel spleens after artificial bacterial infection. After assembly, 110,328 unigenes were obtained with 44.42% GC content. A total of 27,098 unigenes were successfully annotated by four public protein databases, namely, Nr, UniProt, KEGG and KOG. Differential expression analysis revealed the existence of 1278 significantly differentially expressed unigenes at 24 h post infection, with 767 unigenes upregulated and 511 unigenes downregulated. After GO and KEGG enrichment analyses, many immune-related GO categories and pathways were significantly enriched. The typical significantly enriched pathways included toll-like receptor signaling pathway, cytokine-cytokine receptor interaction and TNF signaling pathway. In addition, 40,027 microsatellites (SSRs) and 52,716 candidate single nucleotide polymorphisms (SNPs) were identified from the infection and control transcriptome libraries. Overall, this transcriptomic analysis provided valuable information for studying the immune response of zigzag eels against bacterial infection.