Transcriptome characterization and SSR discovery in large-scale loach Paramisgurnus dabryanus (Cobitidae, Cypriniformes)

Impacts of dietary supplementation of bamboo vinegar and charcoal powder on growth performance, intestinal morphology, and gut microflora of large-scale loach Paramisgurnus dabryanus

Although accumulating data demonstrate that dietary supplementation of bamboo vinegar and charcoal powder (BVC) can significantly strengthen the immunity and boost the growth of domestic animals, its application potency still awaits verification in commercial fish species such as the large-scale loach Paramisgurnus dabryanus. The impacts of 90-day dietary supplementation of 1% and 2% BVC on survival and growth performance of the loach, and on the intestinal morphological characteristics and gut microflora were analyzed. Our data show that the large-scale loach supplied with BVC at the experimental doses had significant higher survival rates and better growth performance (indicated by greater weight gain (1.13-1.14 times), higher specific growth rate (1.04 times), and lower feed conversion ratio (0.88-0.89 times)) compared to that of the control (P<0.05). Histological examination revealed significant longer villus (3.22-5.54 times), deeper crypt (1.77-1.87 times), and thicker muscle (1.59-3.17 times) in the intestines of large-scale loach fed with BVC (P<0.05). Furthermore, we found that the gut microflora consisted of significantly fewer proportions of potential pathogenic bacterial species (Aeromonas veronii and Escherichia coli), but significantly greater proportions of beneficial microbes (Lactococus raffinolactis and Faecalibacterium prausnitzii). Therefore, dietary intake of BVC can promote intestinal tract development and optimize gut microflora, by which the survival and growth of large-scale loach may be improved.

Development of Microsatellite Markers Based on Transcriptome Sequencing and Evaluation of Genetic Diversity in Swimming Crab (Portunus trituberculatus)

P. trituberculatus is an economically important mariculture species in China. Evaluating its genetic diversity and population structure can contribute to the exploration of germplasm resources and promote sustainable aquaculture production. In this study, a total of 246,243 SSRs were generated by transcriptome sequencing of P. trituberculatus. Among the examined 254,746 unigenes, 66,331 had more than one SSR. Among the different SSR motif types, dinucleotide repeats (110,758, 44.98%) were the most abundant. In 173 different base repeats, A/T (96.86%), AC/GT (51.46%), and ACC/GGT (26.20%) were dominant in mono-, di-, and trinucleotide, respectively. GO annotations showed 87,079 unigenes in 57 GO terms. Cellular process, cell, and binding were the most abundant terms in biological process, cellular component, and molecular function categories separately. A total of 34,406 annotated unigenes were classified into 26 functional categories according to the functional annotation analysis of KOG, of which “general function prediction only” was the biggest category (6,028 unigenes, 17.52%). KEGG pathway annotations revealed the clustering of 34,715 unigenes into 32 different pathways. Nineteen SSRs were identified as polymorphic and, thus, used to assess the genetic diversity and structure of 240 P. trituberculatus individuals from four populations in the Bohai Sea. Genetic parameter analysis showed a similar level of genetic diversity within wild populations, and the cultured population indicated a reduction in genetic diversity compared with wild populations. The pairwise F_ST values were between 0.001 and 0.04 with an average of 0.0205 (p < 0.05), suggesting a low but significant level of genetic differentiation among the four populations. Structure analysis demonstrated that the four populations were classified into two groups including the cultured group and other populations. The phylogenetic tree and PCA revealed that a vast number of samples were clustered together and that cultivated individuals were distributed more centrally than wild individuals. The findings contribute to the further assessment of germplasm resources and assist to provide valuable SSRs for marker-assisted breeding of P. trituberculatus in the future.

Discovery of alternatively spliced isoforms and long non-coding RNA in full length brain transcriptomes of anadromous Hilsa shad, Tenualosa ilisha (Hamilton, 1822)

BACKGROUND

Full length transcriptomes, achieved through long-read sequencing, along with the isoform analysis can reveal complexities in the gene expression profiles, as well as annotate the transcriptomes of non-model organisms.

METHODS AND RESULT

Full length transcripts of brain transcriptome of Tenualosa ilisha, Hilsa shad, were generated through PacBio single molecule real-time sequencing and were characterized. A total of 8.30 Gb clean reads were generated, with PacBio RSII, which resulted in 57,651 high quality consensus transcripts. After removing redundant reads, a total of 19,220 high-quality non-redundant transcripts and 17,341 full length ORF transcripts were classified to 7522 putative ortholog groups. Genes involved in various neural pathways were identified. In addition, isoform clusters and lncRNAs were discovered, along with Hilsa specific transcripts with coding frames and 29,147 SSRs in 944 transcripts (1141 annotated).

CONCLUSION

The present study provided, for the first time, a comprehensive view of the alternative isoforms of genes and transcriptome complexity in Hilsa shad brain and forms a rich resource for functional studies in brain of this anadromous fish.

Comparative transcriptome profiles of large and small bodied large-scale loaches cultivated in paddy fields

Fish culture in paddy fields is a traditional aquaculture mode, which has a long history in East Asia. Large-scale loach (Paramisgurnus dabryanus) fast growth is suitable for paddy fields aquaculture in China. The objective of this study was to identify differential expression genes (DEGs) in the brain, liver and muscle tissues between large (LG, top 5% of maximum total length) and small (SG, top 5% of minimum total length) groups using RNA-seq. In total, 150 fish were collected each week and 450 fish were collected at twelfth week from three paddy fields for all the experimental. Histological observation found that the muscle fibre diameter of LG loaches was greater than that of SG loaches. Transcriptome results revealed that the high expression genes (HEGs) in LG loaches (fold change ≥ 2, p < 0.05) were mainly concentrated in metabolic pathways, such as "Thyroid hormone signalling pathway", "Citrate cycle (TCA cycle)", "Carbon metabolism", "Fatty acid metabolism", and "Cholesterol metabolism", and the HEGs in SG loaches were enriched in the pathways related to environmental information processing such as "Cell adhesion molecules (CAMs)", "ECM- receptor interaction" and "Rap1 signalling pathway"; cellular processes such as "Tight junction", "Focal adhesion", "Phagosome" and "Adherens junction". Furthermore, IGFs gene family may play an important role in loach growth for their different expression pattern between the two groups. These findings can enhance our understanding about the molecular mechanism of different growth and development levels of loaches in paddy fields.

Characterization of Myostain (MSTN) and Myogenic Differentiation Antigen (MyoD) and the Effect of Dexamethasone on Their Expression in Large-Scale Loach Paramisgurnus dabryanus

Myostatin (MSTN) and myogenic differentiation antigen (MyoD) play an essential role in specification and differentiation of skeletal muscle. However, the role of stress in the regulation of MyoD and MSTN has not been fully revealed and more evidence should be provided. Here, we reported the cloning and expressional analyses of MSTN and MyoD in Large-scale Loach Paramisgurnus dabryanus (hereafter PdMSTN and PdMyoD). Injecting individuals with 0, 60, 600, and 1,200 μg/kg dexamethasone (DXM) for five consecutive days resulted in a dose-dependent change of PdMSTN and PdMyoD expression. The expression of PdMSTN was upregulated with increasing DXM concentrations, while PdMyoD expression was downregulated. The changes in the expression of these genes at different time points for 10 consecutive days were studied after individuals were treated with 600 μg/kg DXM. Compared with the control group, PdMSTN expression decreased and PdMyoD expression increased before 12 h, and both PdMSTN and PdMyoD expression levels increased at 24 h, which was significantly higher than those in control group. At a prolonged treatment of 5-10 d, expression levels of PdMSTN and PdMyoD had significantly reduced. The results indicate that both PdMyoD and PdMSTN are involved in DXM-induced stress in Large-scale Loach.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

A High-Density Genetic Linkage Map and QTL Mapping for Sex and Growth-Related Traits of Large-Scale Loach (Paramisgurnus dabryanus)

Large-scale loach (Paramisgurnus dabryanus) is a commercially important species in East Asia; however, the cultured population that exhibited degradation of germplasm resource cannot meet the market needs, and the genome resources for P. dabryanus are still lacking. In this study, the first high-density genetic map of P. dabryanus was constructed using 15,830 SNP markers based on high-throughput sequencing with an improved SLAF-seq strategy. The quantitative trait locus (QTL) mapping for sex, growth, and morphology traits was performed for the first time. The genetic map spanned 4,657.64 cM in length with an average inter-marker distance of 0.30 cM. QTL mapping and association analysis identified eight QTLs of growth traits, nine QTLs of morphology traits, and five QTLs of sex-related traits, respectively. Interestingly, the most significant QTLs for almost all the traits were concentrated on the same linkage group LG11. Seven candidate markers and 12 potentially key genes, which were associated with sex determination and growth, were identified within the overlapped QTL regions on LG11. Further, the first genome survey analysis of P. dabryanus was performed which represents the first step toward fully decoding the P. dabryanus genome. The genome scaffolds were anchored to the high-density linkage map, spanning 960.27 Mb of P. dabryanus reference genome. The collinearity analysis revealed a high level of collinearity between the genetic map and the reference genome of P. dabryanus. Moreover, a certain degree of homology was observed between large-scale loach and zebrafish using comparative genomic analysis. The constructed high-density genetic map was an important basis for QTL fine mapping, genome assembly, and genome comparison. The present study will provide a valuable resource for future marker-assisted breeding, and further genetic and genomic researches in P. dabryanus.

Characterization of transcriptional responses mediated by benzo[a]pyrene stress in a new marine fish model of goby, Mugilogobius chulae

Benzo[a]pyrene (BaP) is one of the most studied targets among polycyclic aromatic hydrocarbons (PAHs). Because of the complexity of the toxicity mechanism in BaP, little is known about the molecular mechanism at the level of transcription of BaP in marine fishes. The primary objective of this study was to investigate the molecular basis of the effects of BaP on marine fish, using Mugilogobius chulae (Smith 1932) as the model. A closed colony of M. chulae was used for the BaP toxicity test. Two fish liver samples per replicate from each group were excised and blended into one sample by pooling an equal amount of liver tissue. Total RNA of all samples was extracted separately. Equal quantities of total RNA from the three replicates of the two groups were pooled for sequencing. The sequencing cDNA libraries were sequenced using Illumina HiSeq 2000 system. Differentially expressed genes were detected with the DEGSeq R package. In total, 52,364,032 and 53,771,748 clean nucleotide reads were obtained in the control and BaP-exposed libraries, respectively, with N50 lengths of 1277 and 1288 bp, respectively. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed a significant enrichment of genes related to detoxification, transportation, and lipid metabolism. We also identified, for the first time, an association between endoplasmic reticulum dysfunction and lipid metabolism resulting from BaP exposure. Using quantitative real-time PCR, some effective molecular biomarkers for monitoring of BaP-polluted seawater were identified. The results demonstrate that BaP enhanced the expression of genes involved in detoxification in M. chulae and inhibited that of genes related to lipid metabolism, possibly by suppressing the expression of numerous ER-related genes involved in fat digestion and absorption.

De novo assembly of Schizothorax waltoni transcriptome to identify immune-related genes and microsatellite markers

Schizothorax waltoni (S. waltoni) is one kind of the subfamily Schizothoracinae and an indigenous economic tetraploid fish to Tibet in China. It is rated as a vulnerable species in the Red List of China's Vertebrates, owing to overexploitation and biological invasion. S. waltoni plays an important role in ecology and local fishery economy, but little information is known about genetic diversity, local adaptation, immune system and so on. Functional gene identification and molecular marker development are the first and essential step for the following biological function and genetics studies. For this purpose, the transcriptome from pooled tissues of three adult S. waltoni was sequenced and analyzed. Using paired-end reads from the Illumina Hiseq4000 platform, 83 103 transcripts with an N50 length of 2337 bp were assembled, which could be further clustered into 66 975 unigenes with an N50 length of 2087 bp. The majority of the unigenes (58 934, 87.99%) were successfully annotated by 7 public databases, and 15 KEGG pathways of immune-related genes were identified for the following functional research. Furthermore, 19 497 putative simple sequence repeats (SSRs) of 1-6 bp unit length were detected from 14 690 unigenes (21.93%) with an average distribution density of 1 : 3.28 kb. We identified 3590 unigenes (5.36%) containing more than one SSR, providing abundant potential polymorphic markers in functional genes. This is the first reported high-throughput transcriptome analysis of S. waltoni, and it would provide valuable genetic resources for the functional genes involved in multiple biological processes, including the immune system, genetic conservation, and molecular marker-assisted breeding of S. waltoni.

Identification of novel EST-SSR markers by transcriptome sequencing in ridgetail white prawn Exopalaemon carinicauda

The ridgetail white prawn Exopalaemon carinicauda is one of the major commercial mariculture species in eastern China. However, only limited molecular markers are currently available due to the lack of genome information, which hinders its genetic and breeding studies. In this study, we identified new simple sequence repeats from transcriptome sequences by Illumina Hiseq 2500 platform. A total of 14273 SSR loci were identified from 130,082 assembled transcripts, with 6590 pairs of PCR primers designed. A total of 12,155 sequences containing SSR were predicted; and 2764 (22.74%) transcripts had significant matches to the NCBI non redundant protein (Nr) database. 11,563 transcripts were assigned into gene ontology (GO) categories. A set of 200 primers selected randomly were synthesized, of which 152 (76.0%) were successfully amplified. Further test with 60 pairs of polymorphic SSR primers to evaluate the genetic diversity of 30 wild populations and 43 loci were polymorphic, which had a polymorphic information content between 0.204 and 0.911. The results enriched genomic resources of E. carinicauda and provided powerful information for future conservation and breeding researches.

Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing

The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.

De novo assembly, characterization and annotation for the transcriptome of Sarcocheilichthys sinensis

The Chinese lake gudgeon Sarcocheilichthys sinensis is a small cyprinid fish with great aquaculture potential both for its edible and ornamental values. Nevertheless, available genomic and transcriptomic information for this fish is extremely deficient. In this study, a normalized cDNA library was constructed using 13 mixed tissues of an adult male S. sinensis, and was sequenced by the Illumina HiSeq2500 platform. De novo assembly was performed using 38,911,511 obtained clean reads, and a total of 147,282 unigenes with an average length of 900 bp were finally achieved. 96.2% of these unigenes were annotated in 9 public databases, and 16 segments of growth-related genes were identified for future studies. In addition, 28,493 unigenes were assigned to 61 subcategories of Gene Ontology (GO), and 10,483 unigenes were assigned to 25 categories of Cluster of Orthologous Group (COG). Moreover, 14,943 unigenes were classified into 225 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 30,666 microsatellites were detected from 17,627 unigenes with an average distribution density of 1:2405 bp. This transcriptome data set will be valuable for researches on discovery, expression and evolution on genes of interest. Meanwhile, the identified microsatellites would be useful tools for genetic and genomic studies in S. sinensis.

Identification of Immune-Related Genes and Development of SSR/SNP Markers from the Spleen Transcriptome of Schizothorax prenanti

Schizothorax prenanti (S. prenanti) is mainly distributed in the upstream regions of the Yangtze River and its tributaries in China. This species is indigenous and commercially important. However, in recent years, wild populations and aquacultures have faced the serious challenges of germplasm variation loss and an increased susceptibility to a range of pathogens. Currently, the genetics and immune mechanisms of S. prenanti are unknown, partly due to a lack of genome and transcriptome information. Here, we sought to identify genes related to immune functions and to identify molecular markers to study the function of these genes and for trait mapping. To this end, the transcriptome from spleen tissues of S. prenanti was analyzed and sequenced. Using paired-end reads from the Illumina Hiseq2500 platform, 48,517 transcripts were isolated from the spleen transcriptome. These transcripts could be clustered into 37,785 unigenes with an N50 length of 2,539 bp. The majority of the unigenes (35,653, 94.4%) were successfully annotated using non-redundant nucleotide sequence analysis (nt), and the non-redundant protein (nr), Swiss-Prot, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. KEGG pathway assignment identified more than 500 immune-related genes. Furthermore, 7,545 putative simple sequence repeats (SSRs), 857,535 single nucleotide polymorphisms (SNPs), and 53,481 insertion/deletion (InDels) were detected from the transcriptome. This is the first reported high-throughput transcriptome analysis of S. prenanti, and it provides valuable genetic resources for the investigation of immune mechanisms, conservation of germplasm, and molecular marker-assisted breeding of S. prenanti.

Molecular characterization of the lgals1 gene in large scale loach Paramisgurnus dabryanus

Galectins constitute a group of lectins with binding specificity for β-galactoside sugars. Galectin-1 is a prototype galectin and the multifunctionality of mammalian galectin-1s is well-known, but only a few of fish galectin-1s have been identified. In this study, we obtained the full-length cDNA and genomic sequence of the galectin-1 gene (designated as Pdlgals1) from large scale loach (Paramisgurnus dabryanus), performed phylogenetic analysis, and characterized the expression pattern and the transcriptional activity of its 5' flanking region. The Pdlgals1 gene contains 4 exons that encode a peptide of 132 amino acids with all the galectin signature motifs. Phylogenetic analysis and sequence alignment indicated that Pdlgals1 is a homologue of human LGALS1. RT-PCR and whole-mount in situ hybridization revealed that Pdlgals1 is mainly expressed in the skin, muscle, intestine and cavum oropharyngeum. Transcriptional activity assays demonstrated that the basal promoter of Pdlgals1 is located in a region from -500bp to its transcriptional start site. Potential binding sites for transcription factors including C/EBP, AP-1, GATA, Oct-1, δEF1, NF-κB, c-Myb, SP-1, AP-2, AML-1α, and AP-4 were identified in the basal promoter, suggesting that these factors are associated with the regulation of Pdlgals1. These results provided clues for further investigation of galectin-1 functions in loaches.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.

Evidence for Adaptation to the Tibetan Plateau Inferred from Tibetan Loach Transcriptomes

Triplophysa fishes are the primary component of the fish fauna on the Tibetan Plateau and are well adapted to the high-altitude environment. Despite the importance of Triplophysa fishes on the plateau, the genetic mechanisms of the adaptations of these fishes to this high-altitude environment remain poorly understood. In this study, we generated the transcriptome sequences for three Triplophysa fishes, that is, Triplophysa siluroides, Triplophysa scleroptera, and Triplophysa dalaica, and used these and the previously available transcriptome and genome sequences from fishes living at low altitudes to identify potential genetic mechanisms for the high-altitude adaptations in Triplophysa fishes. An analysis of 2,269 orthologous genes among cave fish (Astyanax mexicanus), zebrafish (Danio rerio), large-scale loach (Paramisgurnus dabryanus), and Triplophysa fishes revealed that each of the terminal branches of the Triplophysa fishes had a significantly higher ratio of nonsynonymous to synonymous substitutions than that of the branches of the fishes from low altitudes, which provided consistent evidence for genome-wide rapid evolution in the Triplophysa genus. Many of the GO (Gene Ontology) categories associated with energy metabolism and hypoxia response exhibited accelerated evolution in the Triplophysa fishes compared with the large-scale loach. The genes that exhibited signs of positive selection and rapid evolution in the Triplophysa fishes were also significantly enriched in energy metabolism and hypoxia response categories. Our analysis identified widespread Triplophysa-specific nonsynonymous mutations in the fast evolving genes and positively selected genes. Moreover, we detected significant evidence of positive selection in the HIF (hypoxia-inducible factor)-1A and HIF-2B genes in Triplophysa fishes and found that the Triplophysa-specific nonsynonymous mutations in the HIF-1A and HIF-2B genes were associated with functional changes. Overall, our study provides new insights into the adaptations and evolution of fishes in the high-altitude environment of the Tibetan Plateau and complements previous findings on the adaptations of mammals and birds to high altitudes.