Draft sequencing and analysis of the genome of pufferfish Takifugu flavidus

Advances in biotechnology and breeding innovations in China's marine aquaculture

Biotechnology is the key driving force behind the sustainable development of aquaculture, as biological innovation would significantly improve the capabilities of aquatic breeding and achieve independent and controllable seeding sources to ensure food safety. In this article, we have analyzed the current status and existing problems of marine aquaculture in China. Based on these data, we have summarized the recent (especially the last 10 years) biotechnological innovation and breeding progress of marine aquaculture in China, including whole genome sequencing, sex-related marker screening, genomic selection, and genome editing, as well as progress of improved marine fish varieties in China. Finally, the perspectives in this field have been discussed, and three future countermeasures have been proposed.

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Anti-Melanogenic Effects of Takifugu flavidus Muscle Hydrolysate in B16F10 Melanoma Cells and Zebrafish

Abnormal melanogenesis can lead to hyperpigmentation. Tyrosinase (TYR), a key rate-limiting enzyme in melanin production, is an important therapeutic target for these disorders. We investigated the TYR inhibitory activity of hydrolysates extracted from the muscle tissue of Takifugu flavidus (TFMH). We used computer-aided virtual screening to identify a novel peptide that potently inhibited melanin synthesis, simulated its binding mode to TYR, and evaluated functional efficacy in vitro and in vivo. TFMH inhibited the diphenolase activities of mTYR, reducing TYR substrate binding activity and effectively inhibiting melanin synthesis. TFMH indirectly reduced cAMP response element-binding protein phosphorylation in vitro by downregulating melanocortin 1 receptor expression, thereby inhibiting expression of the microphthalmia-associated transcription factor, further decreasing TYR, tyrosinase related protein 1, and dopachrome tautomerase expression and ultimately impeding melanin synthesis. In zebrafish, TFMH significantly reduced black spot formation. TFMH (200 μg/mL) decreased zebrafish TYR activity by 43% and melanin content by 52%. Molecular dynamics simulations over 100 ns revealed that the FGFRSP (T-6) peptide stably binds mushroom TYR via hydrogen bonds and ionic interactions. T-6 (400 μmol/L) reduced melanin content in B16F10 melanoma cells by 71% and TYR activity by 79%. In zebrafish, T-6 (200 μmol/L) inhibited melanin production by 64%. TFMH and T-6 exhibit good potential for the development of natural skin-whitening cosmetic products.

Interrelationships and biogeography of the New World pufferfish genus Sphoeroides (Tetraodontiformes: Tetraodontidae) inferred using ultra-conserved DNA elements

Colonization of the New World by marine taxa has been hypothesized to have occurred through the Tethys Sea or by crossing the East Pacific Barrier. To better understand patterns and timing of diversification, geological events can be coupled with time calibrated phylogenetic hypotheses to infer major drivers of diversification. Phylogenetic relationships among members of Sphoeroides, a genus of four toothed pufferfishes (Tetraodontiformes: Tetraodontidae) which are found nearly exclusively in the New World (eastern Pacific and western Atlantic), were reconstructed using sequences from ultra-conserved DNA elements, nuclear markers with clear homology among many vertebrate taxa. Hypotheses derived from concatenated maximum-likelihood and species tree summary methods support a paraphyletic Sphoeroides, with Colomesus deeply nested within the genus. Analyses also revealed S. pachygaster, a pelagic species with a cosmopolitan distribution, as the sister taxon to the remainder of Sphoeroides and recovered distinct lineages within S. pachygaster, indicating that this cosmopolitan species may represent a species complex. Ancestral range reconstruction may suggest the genus colonized the New World through the eastern Pacific before diversifying in the western Atlantic, though date estimates for these events are uncertain due to the lack of reliable fossil record for the genus.

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

BACKGROUND

Interspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.

RESULTS

We firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F1 and allotetraploid (F22 and F24) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.

CONCLUSIONS

Our results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

Genome Analysis of Lagocephalus sceleratus: Unraveling the Genomic Landscape of a Successful Invader

The Tetraodontidae family encompasses several species which attract scientific interest in terms of their ecology and evolution. The silver-cheeked toadfish (Lagocephalus sceleratus) is a well-known “invasive sprinter” that has invaded and spread, in less than a decade, throughout the Eastern and part of the Western Mediterranean Sea from the Red Sea through the Suez Canal. In this study, we built and analysed the first near-chromosome level genome assembly of L. sceleratus and explored its evolutionary landscape. Through a phylogenomic analysis, we positioned L. sceleratus closer to T. nigroviridis, compared to other members of the family, while gene family evolution analysis revealed that genes associated with the immune response have experienced rapid expansion, providing a genetic basis for studying how L. sceleratus is able to achieve highly successful colonisation. Moreover, we found that voltage-gated sodium channel (NaV 1.4) mutations previously connected to tetrodotoxin resistance in other pufferfishes are not found in L. sceleratus, highlighting the complex evolution of this trait. The high-quality genome assembly built here is expected to set the ground for future studies on the species biology.

A Novel Angiotensin-I-Converting Enzyme (ACE) Inhibitory Peptide from Takifugu flavidus

Alcalase, neutral protease, and pepsin were used to hydrolyze the skin of Takifugu flavidus. The T. flavidus hydrolysates (TFHs) with the maximum degree of hydrolysis (DH) and angiotensin-I-converting enzyme (ACE)-inhibitory activity were selected and then ultra-filtered to obtain fractions with components of different molecular weights (MWs) (<1, 1-3, 3-10, 10-50, and >50 kDa). The components with MWs < 1 kDa showed the strongest ACE-inhibitory activity with a half-maximal inhibitory concentration (IC50) of 0.58 mg/mL. Purification and identification using semi-preparative liquid chromatography, Sephadex G-15 gel chromatography, RP-HPLC, and LC-MS/MS yielded one new potential ACE-inhibitory peptide, PPLLFAAL (non-competitive suppression mode; IC50 of 28 μmmol·L-1). Molecular docking and molecular dynamics simulations indicated that the peptides should bind well to ACE and interact with amino acid residues and the zinc ion at the ACE active site. Furthermore, a short-term assay of antihypertensive activity in spontaneously hypertensive rats (SHRs) revealed that PPLLFAAL could significantly decrease the systolic blood pressure (SBP) and diastolic blood pressure (DBP) of SHRs after intravenous administration. These results suggested that PPLLFAAL may have potential applications in functional foods or pharmaceuticals as an antihypertensive agent.

Chromosome-level genome assembly of the greenfin horse-faced filefish (Thamnaconus septentrionalis) using Oxford Nanopore PromethION sequencing and Hi-C technology

The greenfin horse-faced filefish, Thamnaconus septentrionalis, is a valuable commercial fish species that is widely distributed in the Indo-West Pacific Ocean. This fish has characteristic blue-green fins, rough skin and a spine-like first dorsal fin. Thamnaconus septentrionalis is of conservation concern because its population has declined sharply, and it is an important marine aquaculture fish species in China. Genomic resources for the filefish are lacking, and no reference genome has been released. In this study, the first chromosome-level genome of T. septentrionalis was constructed using nanopore sequencing and Hi-C technology. A total of 50.95 Gb polished nanopore sequences were generated and were assembled into a 474.31-Mb genome, accounting for 96.45% of the estimated genome size of this filefish. The assembled genome contained only 242 contigs, and the achieved contig N50 was 22.46 Mb, a surprisingly high value among all sequenced fish species. Hi-C scaffolding of the genome resulted in 20 pseudochromosomes containing 99.44% of the total assembled sequences. The genome contained 67.35 Mb of repeat sequences, accounting for 14.2% of the assembly. A total of 22,067 protein-coding genes were predicted, 94.82% of which were successfully annotated with putative functions. Furthermore, a phylogenetic tree was constructed using 1,872 single-copy orthologous genes, and 67 unique gene families were identified in the filefish genome. This high-quality assembled genome will be a valuable resource for a range of future genomic, conservation and breeding studies of T. septentrionalis.

Nucleotide composition of transposable elements likely contributes to AT/GC compositional homogeneity of teleost fish genomes

BACKGROUND

Teleost fish genome size has been repeatedly demonstrated to positively correlate with the proportion of transposable elements (TEs). This finding might have far-reaching implications for our understanding of the evolution of nucleotide composition across vertebrates. Genomes of fish and amphibians are GC homogenous, with non-teleost gars being the single exception identified to date, whereas birds and mammals are AT/GC heterogeneous. The exact reason for this phenomenon remains controversial. Since TEs make up significant proportions of genomes and can quickly accumulate across genomes, they can potentially influence the host genome with their own GC content (GC%). However, the GC% of fish TEs has so far been neglected.

RESULTS

The genomic proportion of TEs indeed correlates with genome size, although not as linearly as previously shown with fewer genomes, and GC% negatively correlates with genome size in the 33 fish genome assemblies analysed here (excluding salmonids). GC% of fish TE consensus sequences positively correlates with the corresponding genomic GC% in 29 species tested. Likewise, the GC contents of the entire repetitive vs. non-repetitive genomic fractions correlate positively in 54 fish species in Ensembl. However, among these fish species, there is also a wide variation in GC% between the main groups of TEs. Class II DNA transposons, predominant TEs in fish genomes, are significantly GC-poorer than Class I retrotransposons. The AT/GC heterogeneous gar genome contains fewer Class II TEs, a situation similar to fugu with its extremely compact and also GC-enriched but AT/GC homogenous genome.

CONCLUSION

Our results reveal a previously overlooked correlation between GC% of fish genomes and their TEs. This applies to both TE consensus sequences as well as the entire repetitive genomic fraction. On the other hand, there is a wide variation in GC% across fish TE groups. These results raise the question whether GC% of TEs evolves independently of GC% of the host genome or whether it is driven by TE localization in the host genome. Answering these questions will help to understand how genomic GC% is shaped over time. Long-term accumulation of GC-poor(er) Class II DNA transposons might indeed have influenced AT/GC homogenization of fish genomes and requires further investigation.

Chromosome genome assembly and annotation of the yellowbelly pufferfish with PacBio and Hi-C sequencing data

Pufferfish are ideal models for vertebrate chromosome evolution studies. The yellowbelly pufferfish, Takifugu flavidus, is an important marine fish species in the aquaculture industry and ecology of East Asia. The chromosome assembly of the species could facilitate the study of chromosome evolution and functional gene mapping. To this end, 44, 27 and 50 Gb reads were generated for genome assembly using Illumina, PacBio and Hi-C sequencing technologies, respectively. More than 13 Gb full-length transcripts were sequenced on the PacBio platform. A 366 Mb genome was obtained with the contig of 4.4 Mb and scaffold N50 length of 15.7 Mb. 266 contigs were reliably assembled into 22 chromosomes, representing 95.9% of the total genome. A total of 29,416 protein-coding genes were predicted and 28,071 genes were functionally annotated. More than 97.7% of the BUSCO genes were successfully detected in the genome. The genome resource in this work will be used for the conservation and population genetics of the yellowbelly pufferfish, as well as in vertebrate chromosome evolution studies.

Comparative genome analysis of 52 fish species suggests differential associations of repetitive elements with their living aquatic environments

BACKGROUND

Repetitive elements make up significant proportions of genomes. However, their roles in evolution remain largely unknown. To provide insights into the roles of repetitive elements in fish genomes, we conducted a comparative analysis of repetitive elements of 52 fish species in 22 orders in relation to their living aquatic environments.

RESULTS

The proportions of repetitive elements in various genomes were found to be positively correlated with genome sizes, with a few exceptions. More importantly, there appeared to be specific enrichment between some repetitive element categories with species habitat. Specifically, class II transposons appear to be more abundant in freshwater bony fish than in marine bony fish when phylogenetic relationship is not considered. In contrast, marine bony fish harbor more tandem repeats than freshwater species. In addition, class I transposons appear to be more abundant in primitive species such as cartilaginous fish and lamprey than in bony fish.

CONCLUSIONS

The enriched association of specific categories of repetitive elements with fish habitats suggests the importance of repetitive elements in genome evolution and their potential roles in fish adaptation to their living environments. However, due to the restriction of the limited sequenced species, further analysis needs to be done to alleviate the phylogenetic biases.

Evolution of the duplicated intracellular lipid-binding protein genes of teleost fishes

Increasing organismal complexity during the evolution of life has been attributed to the duplication of genes and entire genomes. More recently, theoretical models have been proposed that postulate the fate of duplicated genes, among them the duplication-degeneration-complementation (DDC) model. In the DDC model, the common fate of a duplicated gene is lost from the genome owing to nonfunctionalization. Duplicated genes are retained in the genome either by subfunctionalization, where the functions of the ancestral gene are sub-divided between the sister duplicate genes, or by neofunctionalization, where one of the duplicate genes acquires a new function. Both processes occur either by loss or gain of regulatory elements in the promoters of duplicated genes. Here, we review the genomic organization, evolution, and transcriptional regulation of the multigene family of intracellular lipid-binding protein (iLBP) genes from teleost fishes. Teleost fishes possess many copies of iLBP genes owing to a whole genome duplication (WGD) early in the teleost fish radiation. Moreover, the retention of duplicated iLBP genes is substantially higher than the retention of all other genes duplicated in the teleost genome. The fatty acid-binding protein genes, a subfamily of the iLBP multigene family in zebrafish, are differentially regulated by peroxisome proliferator-activated receptor (PPAR) isoforms, which may account for the retention of iLBP genes in the zebrafish genome by the process of subfunctionalization of cis-acting regulatory elements in iLBP gene promoters.