Genome-Wide Identification, Sequence Alignment, and Transcription of Five Sex-Related Genes in Largemouth Bass (Micropterus Salmoides).

BACKGROUND: Largemouth bass (Micropterus Salmoides) is an economically important fish species in China. Most research has focused on its growth, disease resistance, and nutrition improvement. However, the sex-determining genes in largemouth bass are still unclear. The transforming growth factor-beta (TGF-ß) gene family, including amh, amhr2 and gsdf, plays an important role in the sex determination and differentiation of various fishes. These genes are potentially involved in sex determination in largemouth bass. METHODS: We performed a systematic analysis of 5 sex-related genes (amh, amhr2, gsdf, cyp19a1, foxl2) in largemouth bass using sequence alignment, collinearity analysis, transcriptome, and quantitative real-time polymerase chain reaction (qRT-PCR). This included a detailed assessment of their sequences, gene structures, evolutionary traits, and gene transcription patterns in various tissues including gonads, and at different developmental stages. RESULTS: Comparative genomics revealed that the 5 sex-related genes were highly conserved in various fish genomes. These genes did not replicate, mutate or lose in largemouth bass. However, some were duplicated (amh, amhr2 and gsdf), mutated (gsdf) or lost (amhr2) in other fishes. Some genes (e.g., gsdf) showed significant differences in genomic sequence between males and females, which may contribute to sex determination and sex differentiation in these fishes. qRT-PCR was applied to quantify transcription profiling of the 5 genes during gonadal development and in the adult largemouth bass. Interestingly, amh, amhr2 and gsdf were predominantly expressed in the testis, while cyp19a1 and foxl2 were mainly transcribed in the ovary. All 5 sex-related genes were differentially expressed in the testes and ovaries from the 56th day post-fertilization (dpf). We therefore speculate that male/female differentiation in the largemouth bass may begin at this critical time-point. Examination of the transcriptome data also allowed us to screen out several more sex-related candidate genes. CONCLUSIONS: Our results provide a valuable genetic resource for investigating the physiological functions of these 5 sex-related genes in sex determination and gonadal differentiation, as well as in the control of gonad stability in adult largemouth bass.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation.

INTRODUCTION: Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation. METHODS: Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation. RESULTS: We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation. CONCLUSION: These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

The first Conus genome assembly reveals a primary genetic central dogma of conopeptides in C. betulinus.

Although there are various Conus species with publicly available transcriptome and proteome data, no genome assembly has been reported yet. Here, using Chinese tubular cone snail (C. betulinus) as a representative, we sequenced and assembled the first Conus genome with original identification of 133 genome-widely distributed conopeptide genes. After integration of our genomics, transcriptomics, and peptidomics data in the same species, we established a primary genetic central dogma of diverse conopeptides, assuming a rough number ratio of ~1:1:1:10s for the total genes: transcripts: proteins: post-translationally modified peptides. This ratio may be special for this worm-hunting Conus species, due to the high diversity of various Conus genomes and the big number ranges of conopeptide genes, transcripts, and peptides in previous reports of diverse Conus species. Only a fraction (45.9%) of the identified conotopeptide genes from our achieved genome assembly are transcribed with transcriptomic evidence, and few genes individually correspond to multiple transcripts possibly due to intraspecies or mutation-based variances. Variable peptide processing at the proteomic level, generating a big diversity of venom conopeptides with alternative cleavage sites, post-translational modifications, and N-/C-terminal truncations, may explain how the 133 genes and ~123 transcripts can generate thousands of conopeptides in the venom of individual C. betulinus. We also predicted many conopeptides with high stereostructural similarities to the putative analgesic ω-MVIIA, addiction therapy AuIB and insecticide ImI, suggesting that our current genome assembly for C. betulinus is a valuable genetic resource for high-throughput prediction and development of potential pharmaceuticals.

The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates.

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.

Effect of aroma-producing yeasts in high-salt liquid-state fermentation soy sauce and the biosynthesis pathways of the dominant esters.

Fermentation with excellent aroma-producing yeasts can enhance the flavour of soy sauce. In this work, Millerozyma farinosa CS2.23, Zygosaccharomyces rouxii CS2.42, and Candida parapsilosis CS2.53 were added to the high-salt liquid-state moromi to promote soy sauce fermentation. All three yeasts improved the TE of soy sauce, the highest of which reached 1.03 g/L with added CS2.42. Other quality indexes of soy sauce, including RS, TA, and AN, were not greatly affected. The volatile esters of soy sauce added to the three yeasts increased by 108.85%, 166.71%, and 113.61% compared with the control through GC-MS analysis. Obviously, CS2.42 had an excellent ability to produce esters. Studying the biosynthesis pathway of esters, CS2.42 has the best esterification ability, while CS2.53 has the advantage of alcoholysis ability. The exploration of the biosynthetic pathway of acetate and ethyl esters has laid a foundation for regulating esters in soy sauce fermentation.

Comprehensive transcriptional changes in the liver of Kanglang white minnow (Anabarilius grahami) in response to the infection of parasite Ichthyophthirius multifiliis.

Abstract: The notorious parasite Ichthyophthirius multifiliis (Ich) has been recorded worldwide in fish species and causes white spot disease, posing major threats and resulting in severe losses to international fish production. Extensively effective strategies for treating Ich are not available yet, and genetic mechanisms of hosts in response to the parasite are still largely unknown. In this study, we selected Kanglang white minnow (KWM, Anabarilius grahami) to examine its liver transcriptional changes after Ich infection, as white spot disease is one bottleneck problem in exploring this economically important species. We divided the experimental fishes into three groups (control, early-infected, and late-infected) to examine differentially expressed genes (DEGs). A total of 831 DEGs were identified and classified into 128 significantly enriched GO (Gene Ontology) terms and 71 significantly enriched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Most of these terms or pathways were functionally enriched in immunity, inflammatory response, and apoptosis, such as nucleotide-binding oligomerization domain-like (NOD-like) receptor signaling, tumor necrosis factor (TNF) signaling, interleukin-17 (IL-17) signaling, and apoptosis pathways. We also identified 178 putative antimicrobial peptides (AMPs) and AMP precursors based on our previously reported genome assembly of KWM, and revealed that the expressional patterns varied according to different types. In summary, our work reported the first comprehensive transcriptional changes in KWM in response to the exogenous infection of Ich, which would lay a solid foundation for in-depth studies on disease defense or resistant strains selection in this valuable fish.

Genome and Transcriptome Sequencing of casper and roy Zebrafish Mutants Provides Novel Genetic Clues for Iridophore Loss.

casper has been a widely used transparent mutant of zebrafish. It possesses a combined loss of reflective iridophores and light-absorbing melanophores, which gives rise to its almost transparent trunk throughout larval and adult stages. Nevertheless, genomic causal mutations of this transparent phenotype are poorly defined. To identify the potential genetic basis of this fascinating morphological phenotype, we constructed genome maps by performing genome sequencing of 28 zebrafish individuals including wild-type AB strain, roy orbison (roy), and casper mutants. A total of 4.3 million high-quality and high-confidence homozygous single nucleotide polymorphisms (SNPs) were detected in the present study. We also identified a 6.0-Mb linkage disequilibrium block specifically in both roy and casper that was composed of 39 functional genes, of which the mpv17 gene was potentially involved in the regulation of iridophore formation and maintenance. This is the first report of high-confidence genomic mutations in the mpv17 gene of roy and casper that potentially leads to defective splicing as one major molecular clue for the iridophore loss. Additionally, comparative transcriptomic analyses of skin tissues from the AB, roy and casper groups revealed detailed transcriptional changes of several core genes that may be involved in melanophore and iridophore degeneration. In summary, our updated genome and transcriptome sequencing of the casper and roy mutants provides novel genetic clues for the iridophore loss. These new genomic variation maps will offer a solid genetic basis for expanding the zebrafish mutant database and in-depth investigation into pigmentation of animals.

High-throughput identification of heavy metal binding proteins from the byssus of chinese green mussel (Perna viridis) by combination of transcriptome and proteome sequencing.

The Byssus, which is derived from the foot gland of mussels, has been proved to bind heavy metals effectively, but few studies have focused on the molecular mechanisms behind the accumulation of heavy metals by the byssus. In this study, we integrated high-throughput transcriptome and proteome sequencing to construct a comprehensive protein database for the byssus of Chinese green mussel (Perna viridis), aiming at providing novel insights into the molecular mechanisms by which the byssus binds to heavy metals. Illumina transcriptome sequencing generated a total of 55,670,668 reads. After filtration, we obtained 53,047,718 clean reads and subjected them to de novo assembly using Trinity software. Finally, we annotated 73,264 unigenes and predicted a total of 34,298 protein coding sequences. Moreover, byssal samples were analyzed by proteome sequencing, with the translated protein database from the foot transcriptome as the reference for further prediction of byssal proteins. We eventually determined 187 protein sequences in the byssus, of which 181 proteins are reported for the first time. Interestingly, we observed that many of these byssal proteins are rich in histidine or cysteine residues, which may contribute to the byssal accumulation of heavy metals. Finally, we picked one representative protein, Pvfp-5-1, for recombinant protein synthesis and experimental verification of its efficient binding to cadmium (Cd2+) ions.

A comparative transcriptomic study on developmental gonads provides novel insights into sex change in the protandrous black porgy (Acanthopagrus schlegelii).

Protandrous black porgy (Acanthopagrus schlegelii) is a popular and valuable commercial marine fish in China and East Asian countries. Controlling and managing its breeding has been an imperative step towards obtaining a sustainable supply of this fish in aquaculture production systems. Therefore, study on the molecular mechanisms of sex change in black porgy has both scientific and commercial importance. Previously, we identified some candidate genes related to sex determination and differentiation from a high-quality genome assembly of the black porgy. In the present study, transcriptome sequencing of developmental gonads (including testis, ovotestis and ovary) of black porgy was performed to further investigate the sex-change mechanisms. Our results showed that the highly expressed male-related genes (dmrt1, piwi1, piwi2, sox9, sox30 and amh) at the male phase were significantly down-regulated to a substantial degree at the intersexual stage, and the female-related genes (jnk1, vasa, wnt4, figla and foxl2) were distinctly up-regulated when the fish grows into a female adult, suggesting the potential roles of these genes in sex change of the black porgy. These data also support a previous hypothesis that the femaleness will be switched on when the testis is entering the degenerated stage through the diminished dmrt1 expression. Our transcriptome data provide a very useful genomic resource for future studies on sex change and practical aquaculture in the black porgy.

A new species of bandy-bandy (Vermicella: Serpentes: Elapidae) from the Weipa region, Cape York, Australia.

Bandy-bandies (genus Vermicella) are small (50-100cm) black and white burrowing elapids with a highly specialised diet of blindsnakes (Typhlopidae). There are currently 5 recognized species in the genus, all located in Australia, with Vermicella annulata the most encountered species with the largest distribution. Morphological and mitochondrial analyses of specimens collected from the Weipa area, Cape York, Queensland reveal the existence of a new species, which we describe as Vermicella parscauda sp. nov. Mitochondrial DNA analysis (16S and ND4) and external morphological characteristics indicate that the closest relatives of the new species are not V. annulata, which also occurs on Cape York, but rather species from Western Australia and the Northern Territory (V. intermedia and V. multifasciata) which, like V. parscauda, occupy monsoon habitats. Internasal scales are present in V. parscauda sp. nov., similar to V. annulata, but V. intermedia and V. multifasciata do not have nasal scales. V. parscauda sp. nov. has 55-94 black dorsal bands and mottled or black ventral scales terminating approximately 2/3rds of the body into formed black rings, suggesting that hyper-banding is a characteristic of the tropical monsoon snakes (V. intermedia, V. multifasciata and V. parscauda). The confined locality, potential habitat disruption due to mining activities, and scarcity of specimens indicates an urgent conservation concern for this species.

Transcriptomic evidence of adaptive tolerance to high environmental ammonia in mudskippers.

Mudskippers are typical amphibious fishes and possess various strategies to ameliorate ammonia toxicity during exposure to environmental ammonia. The present study aimed to provide transcriptomic evidence through profiling the gill and liver transcriptomes of Boleophthalmus pectinirostris (BP) and Periophthalmus magnuspinnatus (PM), which were subjected to treatment with high environmental ammonia for up to 72 h. The results of gene function annotation showed that most of the differentially expressed genes were involved in metabolic pathways. After ammonia exposure, the protein and amino acid metabolism related genes in mudskippers were down-regulated, and PM had more down-regulated genes than BP. The expression levels of several representative genes involved in ammonia excretion in the gill were commonly increased. Interestingly, NH4+ transporting and H+ excreting related genes, including Na+/K+(NH4+)/2Cl- cotransporter (nkcc), Na+/K+(NH4+)-ATPase (nka), carbonic anhydrase 2 (ca2), H+-ATPase, Na+/H+ (NH4+)-exchanger (nhe), and carbonic anhydrase 15 (ca15), were up-regulated more significantly in BP than PM; however, the transcription levels of Rhesus glucoprotein b (Rhbg) and Rhesus glucoprotein c1 (Rhcg1), which constitute the NH3 transporting channels, were up-regulated more significantly in PM than BP. Furthermore, the present study provides molecular evidence for how mudskippers adopt partial amino acid catabolism to decrease the production of endogenous ammonia under high environmental ammonia loading.

Assessing the genetic diversity of the critically endangered Chinese sturgeon Acipenser sinensis using mitochondrial markers and genome-wide single-nucleotide polymorphisms from RAD-seq.

As a living fossil, the endangered Chinese sturgeon (Acipenser sinensis) has been considered a national treasure in China. Here, the famous Gezhouba Dam and Three Gorges Dam on the Yangtze River were built in 1988 and 2006, for economic purposes. The natural population of Chinese sturgeon has declined since then, as these dams block its migratory route to the original spawning grounds in the middle reaches of the Yangtze River. In 2013 and 2014, there was an absence of spawning where it typically happened near the Gezhouba Dam. Nevertheless, from April to June in 2015, over 1,000 larvae with different body lengths (10-35 cm) were detected along the Shanghai Yangtze Estuary; but only little is currently known about the population genetic structure of the Chinese sturgeon. Herein, we inferred population genetic parameters from 462 available Chinese sturgeon specimens based on a 421-bp fragment of the mitochondrial DNA (mtDNA) D-loop region and 1,481,620 SNPs (single-nucleotide polymorphisms) generated by restriction site-associated DNA sequencing (RAD-seq). For the D-loop dataset, 15 haplotypes were determined. Randomly picked 23 individuals, representing the 15 D-loop haplotype groups, were subsequently used for further RAD-seq validation. The average nucleotide diversity calculated from the mtDNA and RAD datasets was 0.0086 and 0.000478, respectively. The overall effective female population size was calculated to be 1,255 to 2,607, and the long-term effective population size was estimated to range from 11,950 to 119,500. We observed that the genetic variability and the effective female population size of the current population in the Yangtze River are severely low, which are similar to the data reported over 10 years ago. The deduced relatively small effective population of female fish, limiting the genetic connectivity among Chinese sturgeon, should be considered a serious threat to this endangered species.

The Complete Mitochondrial Genome of Glyptothorax macromaculatus Provides a Well-Resolved Molecular Phylogeny of the Chinese Sisorid Catfishes.

Previous phylogenetic analyses of the Chinese sisorid catfishes have either been poorly resolved or have not included all the 12 sisorid genera. Here, we successfully assembled the first complete mitochondrial genome of the sisorid fish Glyptothorax macromaculatus. Based on this novel mitochondrial genome and previously published mitochondrial genomes in the Sisoridae, we generated maximum likelihood and Bayesian phylogenies. We dated our preferred topology using fossil calibration points. We also tested the protein-coding genes in the mitochondrial genomes of the glyptosternoid fishes for signals of natural selection by comparing the nucleotide substitution rate along the branch ancestral to the glyptosternoid fishes to other branches in our topology. The mitochondrial sequence structure of G. macromaculatus was similar to those known from other vertebrates, with some slight differences. Our sisorid phylogenies were well-resolved and well-supported, with exact congruence between the different phylogenetic methods. This robust phylogeny clarified the relationships among the Chinese sisorid genera and strongly supported the division of the family into three main clades. Interestingly, the glyptosternoid divergence time predicted by our molecular dating analysis coincided with the uplift of the Tibetan Plateau, suggesting that geology may have influenced speciation in the Sisoridae. Among the mitochondrial protein-coding genes, atp8 may have most rapidly evolved, and atp6 may have been subjected to positive selection pressure to adapt to high elevations. In summary, this study provided novel insights into the phylogeny, evolution and high-altitude adaptions of the Chinese sisorid fishes.

A Comparative Genomic and Transcriptomic Survey Provides Novel Insights into N-Acetylserotonin Methyltransferase (ASMT) in Fish.

Melatonin is a multifunctional bioactive molecule that plays comprehensive physiological roles in all living organisms. N-acetylserotonin methyltransferase (ASMT, also known as hydroxyindole O-methyltransferase or HIOMT) is the final enzyme for biosynthesis of melatonin. Here, we performed a comparative genomic and transcriptomic survey to explore the ASMT family in fish. Two ASMT isotypes (ASMT1 and ASMT2) and a new ASMT-like (ASMTL) are all extracted from teleost genomes on the basis of phylogenetic and synteny analyses. We confirmed that C-terminal of the ASMTL proteins (ASMTL-ASMT) is homology to the full length of ASMT1 and ASMT2. Our results also demonstrate that the two ASMT isotypes and their distribution in teleosts seem to be the result of combinations of whole-genome duplication (WGD) and gene loss. Differences were also observed in tissue distribution and relative transcript abundances of ASMT1, ASMT2 and ASMTL through transcriptomic analysis. Protein sequence alignment and 3D structure prediction of ASMTs and ASMTL suggest differential roles for these ASMT genes. In summary, our current work provides novel insights into the ASMT genes in fish by combination of genomic and transcriptomic data.

De novo assembly and comparative transcriptome analysis of the foot from Chinese green mussel (Perna viridis) in response to cadmium stimulation.

The Chinese green mussel, Perna viridis, is a marine bivalve with important economic values as well as biomonitoring roles for aquatic pollution. Byssus, secreted by the foot gland, has been proved to bind heavy metals effectively. In this study, using the RNA sequencing technology, we performed comparative transcriptomic analysis on the mussel feet with or without inducing by cadmium (Cd). Our current work is aiming at providing insights into the molecular mechanisms of byssus binding to heavy metal ions. The transcriptome sequencing generated a total of 26.13-Gb raw data. After a careful assembly of clean data, we obtained a primary set of 105,127 unigenes, in which 32,268 unigenes were annotated. Based on the expression profiles, we identified 9,048 differentially expressed genes (DEGs) between Cd treatment (50 or 100 µg/L) at 48 h and the control, suggesting an extensive transcriptome response of the mussels during the Cd stimulation. Moreover, we observed that the expression levels of 54 byssus protein coding genes increased significantly after the 48-h Cd stimulation. In addition, 16 critical byssus protein coding genes were picked for profiling by quantitative real-time PCR (qRT-PCR). Finally, we reached a primary conclusion that high content of tyrosine (Tyr), cysteine (Cys), histidine (His) residues or the special motif plays an important role in the accumulation of heavy metals in byssus. We also proposed an interesting model for the confirmed byssal Cd accumulation, in which biosynthesis of byssus proteins may play simultaneously critical roles since their transcription levels were significantly elevated.

Comparative Transcriptomic Study of Muscle Provides New Insights into the Growth Superiority of a Novel Grouper Hybrid.

Grouper (Epinephelus spp.) is a group of fish species with great economic importance in Asian countries. A novel hybrid grouper, generated by us and called the Hulong grouper (Hyb), has better growth performance than its parents, E. fuscoguttatus (Efu, ♀) and E. lanceolatus (Ela, ♂). We previously reported that the GH/IGF (growth hormone/insulin-like growth factor) system in the brain and liver contributed to the superior growth of the Hyb. In this study, using transcriptome sequencing (RNA-seq) and quantitative real-time PCR (qRT-PCR), we analyzed RNA expression levels of comprehensive genes in the muscle of the hybrid and its parents. Our data showed that genes involved in glycolysis and calcium signaling in addition to troponins are up-regulated in the Hyb. The results suggested that the activity of the upstream GH/IGF system in the brain and liver, along with the up-regulated glycolytic genes as well as ryanodine receptors (RyRs) and troponins related to the calcium signaling pathway in muscle, led to enhanced growth in the hybrid grouper. Muscle contraction inducing growth could be the major contributor to the growth superiority in our novel hybrid grouper, which may be a common mechanism for hybrid superiority in fishes.

High-throughput identification of novel conotoxins from the Chinese tubular cone snail (Conus betulinus) by multi-transcriptome sequencing.

BACKGROUND: The venom of predatory marine cone snails mainly contains a diverse array of unique bioactive peptides commonly referred to as conopeptides or conotoxins. These peptides have proven to be valuable pharmacological probes and potential drugs because of their high specificity and affinity to important ion channels, receptors and transporters of the nervous system. Most previous studies have focused specifically on the conopeptides from piscivorous and molluscivorous cone snails, but little attention has been devoted to the dominant vermivorous species. RESULTS: The vermivorous Chinese tubular cone snail, Conus betulinus, is the dominant Conus species inhabiting the South China Sea. The transcriptomes of venom ducts and venom bulbs from a variety of specimens of this species were sequenced using both next-generation sequencing and traditional Sanger sequencing technologies, resulting in the identification of a total of 215 distinct conopeptides. Among these, 183 were novel conopeptides, including nine new superfamilies. It appeared that most of the identified conopeptides were synthesized in the venom duct, while a handful of conopeptides were identified only in the venom bulb and at very low levels. CONCLUSIONS: We identified 215 unique putative conopeptide transcripts from the combination of five transcriptomes and one EST sequencing dataset. Variation in conopeptides from different specimens of C. betulinus was observed, which suggested the presence of intraspecific variability in toxin production at the genetic level. These novel conopeptides provide a potentially fertile resource for the development of new pharmaceuticals, and a pathway for the discovery of new conotoxins.

Prediction of Toxin Genes from Chinese Yellow Catfish Based on Transcriptomic and Proteomic Sequencing.

Fish venom remains a virtually untapped resource. There are so few fish toxin sequences for reference, which increases the difficulty to study toxins from venomous fish and to develop efficient and fast methods to dig out toxin genes or proteins. Here, we utilized Chinese yellow catfish (Pelteobagrus fulvidraco) as our research object, since it is a representative species in Siluriformes with its venom glands embedded in the pectoral and dorsal fins. In this study, we set up an in-house toxin database and a novel toxin-discovering protocol to dig out precise toxin genes by combination of transcriptomic and proteomic sequencing. Finally, we obtained 15 putative toxin proteins distributed in five groups, namely Veficolin, Ink toxin, Adamalysin, Za2G and CRISP toxin. It seems that we have developed a novel bioinformatics method, through which we could identify toxin proteins with high confidence. Meanwhile, these toxins can also be useful for comparative studies in other fish and development of potential drugs.