Search

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.

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.

As one economically important fish in the southeastern Himalayas, the giant devil catfish (Bagarius yarrelli) has been known for its extraordinarily large body size. It can grow up to 2 m, whereas the non-Bagarius sisorids only reach 10–30 cm. Another outstanding characteristic of Bagarius species is the salmonids-like reddish flesh color. Both body size and flesh color are interesting questions in science and also valuable features in aquaculture that worth of deep investigations. Bagarius species therefore are ideal materials for studying body size evolution and color depositions in fish muscles, and also potential organisms for extensive utilization in Asian freshwater aquaculture. In a combination of Illumina and PacBio sequencing technologies, we de novo assembled a 571-Mb genome for the giant devil catfish from a total of 153.4-Gb clean reads. The scaffold and contig N50 values are 3.1 and 1.6 Mb, respectively. This genome assembly was evaluated with 93.4% of Benchmarking Universal Single-Copy Orthologs completeness, 98% of transcripts coverage, and highly homologous with a chromosome-level-based genome of channel catfish (Ictalurus punctatus). We detected that 35.26% of the genome assembly is composed of repetitive elements. Employing homology, de novo, and transcriptome-based annotations, we annotated a total of 19,027 protein-coding genes for further use. In summary, we generated the first high-quality genome assembly of the giant devil catfish, which provides an important genomic resource for its future studies such as the body size and flesh color issues, and also for facilitating the conservation and utilization of this valuable catfish.

Haematococcus pluvialis is a freshwater species of Chlorophyta, family Haematococcaceae. It is well known for its capacity to synthesize high amounts of astaxanthin, which is a strong antioxidant that has been utilized in aquaculture and cosmetics. To improve astaxanthin yield and to establish genetic resources for H. pluvialis, we performed whole-genome sequencing, assembly, and annotation of this green microalga. A total of 83.1 Gb of raw reads were sequenced. After filtering the raw reads, we subsequently generated a draft assembly with a genome size of 669.0 Mb, a scaffold N50 of 288.6 kb, and predicted 18,545 genes. We also established a robust phylogenetic tree from 14 representative algae species. With additional transcriptome data, we revealed some novel potential genes that are involved in the synthesis, accumulation, and regulation of astaxanthin production. In addition, we generated an isoform-level reference transcriptome set of 18,483 transcripts with high confidence. Alternative splicing analysis demonstrated that intron retention is the most frequent mode. In summary, we report the first draft genome of H. pluvialis. These genomic resources along with transcriptomic data provide a solid foundation for the discovery of the genetic basis for theoretical and commercial astaxanthin enrichment.

The uplift of the Tibetan Plateau altered the environmental conditions of the local area substantially. Here, we conducted a comprehensive investigation based on morphometrics, population genomics, and climatic factors to evaluate phenotypic and genome-level variations in a radiation of Chimarrichthys catfish endemic to the Plateau. Discriminant function analysis showed phenotypic differences of Chimarrichthys between rivers with respect to elevation. Genetic structure analysis based on 6606 single nucleotide polymorphisms (SNPs) deduced genetic differences between rivers, and species delimitation indicated that the Chimarrichthys fish complex could be divided into three species. Restriction site-associated DNA tags were mapped to the gene sets of Glyptosternon maculatum, and matches were searched against databases for Gene Ontology annotation. Genomic regions exhibiting marked differences among localities represented a range of biological functions, including growth (gdf11), bone development (bmp8a), cellular response to light stimulus (opn3), regulation of the rhodopsin-mediated signalling pathway (grk1), immune response (rag1 and ung), reproductive process (antxr2), and regulation of intracellular iron levels (ireb2). The tag44126, where gene gdf11 is located, was identified as an outlier exhibiting divergence between rivers with altitude differences, and the SNP is thymine (T) in Dadu and Yalong River (~2700 m), but guanine (G) in Jinsha and Qingyi rivers (~2200 and ~ 684 m), suggesting a possible effect of altitude on its differentiation.

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.

Background Seasonal migration is one of the most spectacular events in nature; however, the molecular mechanisms related to this phenomenon have not been investigated in detail. The Chinese tapertail, or Japanese grenadier anchovy, Coilia nasus, is a valuable migratory fish of high economic importance and special migratory dimorphism (with certain individuals as non-migratory residents). Results In this study, an 870.0-Mb high-quality genome was assembled by the combination of Illumina and Pacific Biosciences sequencing. Approximately 812.1 Mb of scaffolds were linked to 24 chromosomes using a high-density genetic map from a family of 104 full siblings and their parents. In addition, population sequencing of 96 representative individuals from diverse areas along the putative migration path identified 150 candidate genes, which are mainly enriched in 3 Ca2+-related pathways. Based on integrative genomic and transcriptomic analyses, we determined that the 3 Ca2+-related pathways are critical for promotion of migratory adaption. A large number of molecular markers were also identified, which distinguished migratory individuals and non-migratory freshwater residents. Conclusions We assembled a chromosome-level genome for the Chinese tapertail anchovy. The genome provided a valuable genetic resource for understanding of migratory adaption and population genetics and will benefit the aquaculture and management of this economically important fish.

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.

Chinese white dolphin (Sousa chinensis), also known as the Indo-Pacific humpback dolphin, has been classified as &ldquo, Vulnerable&rdquo, on the IUCN Red List of Threatened Species. It is a special cetacean species that lives in tropical and subtropical nearshore waters, with significant differences from other cetaceans. Here, we sequenced and assembled a draft genome of the Chinese white dolphin with a total length of 2.3 Gb and annotation of 18,387 protein-coding genes. Genes from certain expanded families are potentially involved in DNA replication and repairing, suggesting that they may be related to adaptation of this marine mammal to nearshore environments. We also discovered that its historical population had undergone a remarkable bottleneck incident before the Mindel glaciation. In addition, a comparative genomic survey on antihypertensive peptides (AHTPs) among five representative mammals with various residential habitats (such as remarkable differences in exogenous ion concentrations and sea depth) revealed that these small bioactive peptides were highly conserved among these examined mammals, and they had the most abundant hits in collagen subunit proteins, especially for two putative AHTP peptides Gly-Leu-Pro (GLP) and Leu-Gly-Pro (LGP). Our genome assembly will be a valuable resource for further genetic researches on adaptive ecology and conservation biology of cetaceans, and for in-depth investigations into bioactive peptides in aquatic and terrestrial mammals for development of peptide-based drugs to treat various human cardiovascular diseases.

Limbs originated from paired fish fins are an important innovation in Gnathostomata. Many studies have focused on limb development-related genes, of which the T-box transcription factor 4 gene (tbx4) has been considered as one of the most essential factors in the regulation of the hindlimb development. We previously confirmed pelvic fin loss in tbx4-knockout zebrafish. Here, we report a high-quality genome assembly of the Japanese eel (Anguilla japonica), which is an economically important fish without pelvic fins. The assembled genome is 1.13 Gb in size, with a scaffold N50 of 1.03 Mb. In addition, we collected 24 tbx4 sequences from 22 teleost fishes to explore the correlation between tbx4 and pelvic fin evolution. However, we observed complete exon structures of tbx4 in several pelvic-fin-loss species such as Ocean sunfish (Mola mola) and ricefield eel (Monopterus albus). More interestingly, an inversion of a special tbx4 gene cluster (brip1-tbx4-tbx2b- bcas3) occurred twice independently, which coincides with the presence of fin spines. A nonsynonymous mutation (M82L) was identified in the nuclear localization sequence (NLS) of the Japanese eel tbx4. We also examined variation and loss of hindlimb enhancer B (HLEB), which may account for pelvic fin loss in Tetraodontidae and Diodontidae. In summary, we generated a genome assembly of the Japanese eel, which provides a valuable genomic resource to study the evolution of fish tbx4 and helps elucidate the mechanism of pelvic fin loss in teleost fishes. Our comparative genomic studies, revealed for the first time a potential correlation between the tbx4 gene cluster and the evolutionary development of toxic fin spines. Because fin spines in teleosts are usually venoms, this tbx4 gene cluster may facilitate the genetic engineering of toxin-related marine drugs.

Blue tilapia (Oreochromis aureus) has been an economically important fish in Asian countries. It can grow and reproduce in both freshwater and brackish water conditions, whereas it is also considered as a significant invasive species around the world. This species has been widely used as the hybridization parent(s) for tilapia breeding with a major aim to produce novel strains. However, available genomic resources are still limited for this important tilapia species. Here, we for the first time sequenced and assembled a draft genome for a seawater cultured blue tilapia (0.92 Gb), with 97.8% completeness and a scaffold N50 of 1.1 Mb, which suggests a relatively high quality of this genome assembly. We also predicted 23,117 protein-coding genes in the blue tilapia genome. Comparisons of predicted antimicrobial peptides between the blue tilapia and its close relative Nile tilapia proved that these immunological genes are highly similar with a genome-wide scattering distribution. As a valuable genetic resource, our blue tilapia genome assembly will benefit for biomedical researches and practical molecular breeding for high resistance to various diseases, which have been a critical problem in the aquaculture of tilapias.

A Yunnan-Guizhou Plateau fish, the Kanglang white minnow (Anabarilius grahami), is a typical “3E” (Endangered, Endemic, and Economic) species in China. Its distribution is limited to Fuxian Lake, the nation’s second deepest lake, with a significant local economic value but a drastically declining wild population. This species has been evaluated as VU (Vulnerable) in the China Species Red List. As one of the “Four Famous Fish” in Yunnan province, the artificial breeding has been achieved since 2003. It has not only re-established its wild natural populations by reintroduction of the artificial breeding stocks, but also brought a wide and popular utilization of this species to the local fish farms. Anabarilius grahami has become one of the main native aquaculture species in Yunnan province, and the artificial production has been emerging in steady growth each year. To promote the conservation and sustainable utilization of this fish, we initiated its whole genome sequencing project using an Illumina Hiseq2500 platform. The assembled genome size of A. grahami is 1.006 Gb, accounting for 98.63% of the estimated genome size (1.020 Gb), with contig N50 and scaffold N50 values of 26.4 kb and 4.41 Mb, respectively. Approximately about 50.38% of the genome was repetitive. A total of 25,520 protein-coding genes were subsequently predicted. A phylogenetic tree based on 4,580 single-copy genes from A. grahami and 18 other cyprinids revealed three well-supported subclades within the Cyprinidae. This is the first inter-subfamily relationship of cyprinids at genome level, providing a simple yet useful framework for understanding the traditional but popular subfamily classification systems. Interestingly, a further population demography of A. grahami uncovered a historical relationship between this fish and Fuxian Lake, suggesting that range expansion or shrinkage of the habitat has had a remarkable impact on the population size of endemic plateau fishes. Additionally, a total of 33,836 simple sequence repeats (SSR) markers were identified, and 11 loci were evaluated for a preliminary genetic diversity analysis in this study, thus providing another useful genetic resource for studying this “3E” species.

Melatonin receptors (MTNRs) play important roles in regulation of circadian rhythms and seasonal reproduction. However, their origin and evolution in vertebrates have not been investigated. Here, we performed a comprehensive examination by comparative genome mining of MTNRs in vertebrates. We successfully extracted 164 putative encoding sequences for MTNRs (including 57 mtnr1a, 59 mtnr1b and 48 mtnr1c) from 45 high-quality representative genomes. Interestingly, the putative expansions of mtnr1a and mtnr1b in zebrafish were also identified in other Cyprinifomes, but not in other orders of teleost. Using phylogenetic interference, we observed this expansion to be clustered into a primitive position of the Actinopterygii, which may be resulted from teleost-specific genome duplication. The C-terminal extension of MTNR1C, predicted to be proteoglycan 4 (PRG4), originated after the speciation of Monotremata or Marsupialia. Our present genomics survey provides novel insights into the evolution of MTNRs in vertebrates and updates our understanding of these proteins.

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.

Bile salt hydrolase (BSH) is a well-known enzyme that has been commonly characterized in probiotic bacteria, as it has cholesterol-lowering effects. However, its molecular investigations are scarce. Here, we build a local database of BSH sequences from Lactobacillaceae (BSH–SDL), and phylogenetic analysis and homology searches were employed to elucidate their comparability and distinctiveness among species. Evolutionary study demonstrates that BSH sequences in BSH–SDL are divided into five groups, named BSH A, B, C, D and E here, which can be the genetic basis for BSH classification and nomenclature. Sequence analysis suggests the differences between BSH-active and BSH-inactive proteins clearly, especially on site 82. In addition, a total of 551 BSHs from 107 species are identified from 451 genomes of 158 Lactobacillaceae species. Interestingly, those bacteria carrying various copies of BSH A or B can be predicted to be potential cholesterol-lowering probiotics, based on the results of phylogenetic analysis and the subtypes that those previously reported BSH-active probiotics possess. In summary, this study elaborates the molecular basis of BSH in Lactobacillaceae systematically, and provides a novel methodology as well as a consistent standard for the identification of the BSH subtype. We believe that high-throughput screening can be efficiently applied to the selection of promising candidate BSH-active probiotics, which will advance the development of healthcare products in cholesterol metabolism.

Melatonin is a pleiotropic molecule with various important physiological roles in vertebrates. l-aromatic amino acid decarboxylase (AAAD) is the second enzyme for melatonin synthesis. By far, a clear-cut gene function of AAAD in the biosynthesis of melatonin has been unclear in vertebrates. Here, we provide novel insights into the evolution of AAAD based on 77 vertebrate genomes. According to our genome-wide alignments, we extracted a total of 151 aaad nucleotide sequences. A phylogenetic tree was constructed on the basis of these sequences and corresponding protein alignments, indicating that tetrapods and diploid bony fish genomes contained one aaad gene and a new aaad-like gene, which formed a novel AAAD family. However, in tetraploid teleosts, there were two copies of the aaad gene due to whole genome duplication. A subsequent synteny analysis investigated 81 aaad sequences and revealed their collinearity and systematic evolution. Interestingly, we discovered that platypus (Ornithorhynchus anatinus), Atlantic cod (Guadus morhua), Mexican tetra (Astyanax mexicanus), and a Sinocyclocheilus cavefish (S. anshuiensis) have long evolutionary branches in the phylogenetic topology. We also performed pseudogene identification and selection pressure analysis; however, the results revealed a deletion of 37 amino acids in Atlantic cod and premature stop codons in the cave-restricted S. anshuiensis and A. mexicanus, suggesting weakening or disappearing rhythms in these cavefishes. Selective pressure analysis of aaad between platypus and other tetrapods showed that rates of nonsynonymous (Ka) and synonymous (Ks) substitutions were higher when comparing the platypus to other representative tetrapods, indicating that, in this semiaquatic mammal, the aaad gene experienced selection during the process of evolution. In summary, our current work provides novel insights into aaad genes in vertebrates from a genome-wide view.

The family of secretory calcium-binding phosphoproteins (SCPPs) have been considered vital to skeletal tissue mineralization. However, most previous SCPP studies focused on phylogenetically distant animals but not on those closely related species. Here we provide novel insights into the coevolution of SCPP genes and fish scales in 10 species from Otophysi. According to their scale phenotypes, these fishes can be divided into three groups, i.e., scaled, sparsely scaled, and scaleless. We identified homologous SCPP genes in the genomes of these species and revealed an absence of some SCPP members in some genomes, suggesting an uneven evolutionary history of SCPP genes in fishes. In addition, most of these SCPP genes, with the exception of SPP1, individually form one or two gene cluster(s) on each corresponding genome. Furthermore, we constructed phylogenetic trees using maximum likelihood method to estimate their evolution. The phylogenetic topology mostly supports two subclasses in some species, such as Cyprinus carpio, Sinocyclocheilus anshuiensis, S. grahamin, and S. rhinocerous, but not in the other examined fishes. By comparing the gene structures of recently reported candidate genes, SCPP1 and SCPP5, for determining scale phenotypes, we found that the hypothesis is suitable for Astyanax mexicanus, but denied by S. anshuiensis, even though they are both sparsely scaled for cave adaptation. Thus, we conclude that, although different fish species display similar scale phenotypes, the underlying genetic changes however might be diverse. In summary, this paper accelerates the recognition of the SCPP family in teleosts for potential scale evolution.

Angiotensin-converting enzymes, ACE and ACE2, are two main elements in the renin&ndash, angiotensin system, with a crucial role in the regulation of blood pressure in vertebrates. Previous studies paid much attention to their physiological functions in model organisms, whereas the studies on other animals and related evolution have been sparse. Our present study performed a comprehensive genomic investigation on ace and ace2 genes in vertebrates. We successfully extracted the nucleotide sequences of ace and ace2 genes from high-quality genome assemblies of 36 representative vertebrates. After construction of their evolutionary tree, we observed that most of the phylogenetic positions are consistent with the species tree, however, certain differences appear in coelacanths and frogs, which may suggest a very slow evolutionary rate in the initial evolution of ace and ace2 in vertebrates. We further compared evolutionary rates within the entire sequences of ace and ace2, and determined that ace2 evolved slightly faster than ace. Meanwhile, we counted that the exon numbers of ace and ace2 in vertebrates are usually 25 and 18 respectively, while certain species may occur exon fusion or disruption to decrease or increase their exon numbers. Interestingly, we found three homologous regions between ace and ace2, suggesting existence of gene duplication during their evolutionary process. In summary, this report provides novel insights into vertebrate ace and ace2 genes through a series of genomic and molecular comparisons.