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

Comparative biochemical characterization of pufferfish saxitoxin and tetrodotoxin-binding protein (PSTBP) homologs in the plasma from four Takifugu species: Conservation of heat-stable PSTBP orthologs having three and two tandemly repeated lipocalin domains in genus Takifugu

To study the relationship between domain characteristics of pufferfish saxitoxin and tetrodotoxin binding protein (PSTBP) proteoforms and their thermal stability, a comparative biochemical characterization of PSTBPs from the plasma of four Takifugu species (T. flavipterus, T. pardalis, T. alboplumbeus and T. rubripes) was conducted by Western blot analysis. The heat-tolerance tetrodotoxin (TTX)-binding ability of PSTBP proteoforms in T. rubripes plasma was verified by ultrafiltration and liquid chromatography tandem mass spectrometry (LC-MS/MS). These results suggest that the heat-stable PSTBP proteoforms, composed of three and two tandemly repeated lipocalin domains, are genetically conserved and ubiquitous in the genus Takifugu. This study builds on our knowledge of the structural and functional properties of PSTBP proteoforms, which is vital for understanding how toxins are transmitted and accumulate in organisms and is essential for evaluating the potential risks of toxins in seafood.

A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Improved high-quality reference genome of red drum facilitates the processes of resistance-related gene exploration

Sciaenops ocellatus is among the most important artificially introduced farmed fish across 11 countries and regions. However, the frequent occurrence of extreme weather events and breeding escapes have placed great pressure on local marine biodiversity and ecosystems. We reported the de novo assembly and annotation with a contig N50 of 28.30 Mb using PacBio HiFi sequencing and Hi-C technologies, which resulted in a 283-fold increase in contig N50 length and improvement in continuity and quality in complex repetitive region for S. ocellatus compared to the previous version. In total, 257.36 Mb of repetitive sequences accounted for 35.48% of the genome, and 22,845 protein-coding genes associated with a BUSCO value of 98.32%, were identified by genome annotation. Moreover, 54 hub genes rapidly responding to hypoosmotic stress were identified by WGCNA. The high-quality chromosome-scale S. ocellatus genome and candidate resistance-related gene sets will not only provide a genomic basis for genetic improvement via molecular breeding, but will also lay an important foundation for investigating the molecular regulation of rapid responses to stress.

Affinity Purification and Molecular Characterization of Angiotensin-Converting Enzyme (ACE)-Inhibitory Peptides from Takifugu flavidus

An affinity chromatography filler of CNBr-activated Sepharose 4B-immobilized ACE was used to purify ACE-inhibitory peptides from Takifugu flavidus protein hydrolysate (<1 kDa). Twenty-four peptides with an average local confidence score (ALC) ≥ 80% from bounded components (eluted by 1 M NaCl) were identified by LC-MS/MS. Among them, a novel peptide, TLRFALHGME, with ACE-inhibitory activity (IC50 = 93.5 µmol·L-1) was selected. Molecular docking revealed that TLRFALHGME may interact with the active site of ACE through H-bond, hydrophobic, and electrostatic interactions. The total binding energy (ΔGbinding) of TLRFALHGME was estimated to be -82.7382 kJ·mol-1 by MD simulations, indicating the favorable binding of peptides with ACE. Furthermore, the binding affinity of TLRFALHGME to ACE was determined by surface plasmon resonance (SPR) with a Kd of 80.9 µmol, indicating that there was a direct molecular interaction between them. TLRFALHGME has great potential for the treatment of hypertension.

High-quality wild barley genome assemblies and annotation with Nanopore long reads and Hi-C sequencing data

Wild barley, from "Evolution Canyon (EC)" in Mount Carmel, Israel, are ideal models for cereal chromosome evolution studies. Here, the wild barley EC_S1 is from the south slope with higher daily temperatures and drought, while EC_N1 is from the north slope with a cooler climate and higher relative humidity, which results in a differentiated selection due to contrasting environments. We assembled a 5.03 Gb genome with contig N50 of 3.53 Mb for wild barley EC_S1 and a 5.05 Gb genome with contig N50 of 3.45 Mb for EC_N1 using 145 Gb and 160.0 Gb Illumina sequencing data, 295.6 Gb and 285.35 Gb Nanopore sequencing data and 555.1 Gb and 514.5 Gb Hi-C sequencing data, respectively. BUSCOs and CEGMA evaluation suggested highly complete assemblies. Using full-length transcriptome data, we predicted 39,179 and 38,373 high-confidence genes in EC_S1 and EC_N1, in which 93.6% and 95.2% were functionally annotated, respectively. We annotated repetitive elements and non-coding RNAs. These two wild barley genome assemblies will provide a rich gene pool for domesticated barley.

New algorithms for accurate and efficient de novo genome assembly from long DNA sequencing reads

Building de novo genome assemblies for complex genomes is possible thanks to long-read DNA sequencing technologies. However, maximizing the quality of assemblies based on long reads is a challenging task that requires the development of specialized data analysis techniques. We present new algorithms for assembling long DNA sequencing reads from haploid and diploid organisms. The assembly algorithm builds an undirected graph with two vertices for each read based on minimizers selected by a hash function derived from the k-mer distribution. Statistics collected during the graph construction are used as features to build layout paths by selecting edges, ranked by a likelihood function. For diploid samples, we integrated a reimplementation of the ReFHap algorithm to perform molecular phasing. We ran the implemented algorithms on PacBio HiFi and Nanopore sequencing data taken from haploid and diploid samples of different species. Our algorithms showed competitive accuracy and computational efficiency, compared with other currently used software. We expect that this new development will be useful for researchers building genome assemblies for different species.

Chromosome-level haplotype-resolved genome assembly for Takifugu ocellatus using PacBio and Hi-C technologies

Takifugu species serve as a model system for evolutionary studies due to their compact genomes and diverse phenotypes. The ocellated puffer (Takifugu ocellatus), characterized by special colouration, is a scarce anadromous species in the genus Takifugu. As an ornamental and tasty fish species, T. ocellatus has moderate economic value. However, the available genomic resources for this pufferfish are still limited. Here, a chromosome-level reference genome, as well as two haploid genomes, was constructed by PacBio HiFi long sequencing and Hi-C technologies. The total length of the reference genome was 375.62 Mb with a contig N50 of 11.55 Mb. The assembled sequences were anchored to 22 chromosomes with an integration efficiency of 93.78%. Furthermore, 28,808 protein-coding genes were predicted. The haplotype-resolved reference genome of T. ocellatus provides a crucial resource for investigating the explosive speciation of the Takifugu genus, such as elucidating evolutionary histories, determining the genetic basis of trait evolution, and supporting future conservation efforts.

Comparative genomics analysis and genome assembly integration with the recombination landscape contribute to Takifugu bimaculatus assembly refinement

Takifugu genus has been brought to the fore in scientific and practical research due to its compact genome, explosive speciation progress and economic value. Here we updated the chromosome-level genome of Takifugu bimaculatus by an ultra-high-density linkage map, a classic and accurate way of chromosome assembly. The map constituted a robust assembly frame, with 92.2% (372.77 Mb) of the draft genome cumulatively placed. With intraspecies and interspecies comparative genomic analysis, we developed a criterion to quantify the differences between assemblies and established a novel way to integrate information from multiple assemblies. The integrated assembly rectified potential mis-assemblies, greatly improving the genome contiguity and correctness. Our results rendered profound information on the genetic recombination of T. bimaculatus and provided new insights into effective genome assembly. The consolidated assembly will be a contributory tool of T. bimaculatus and broadly across the Takifugu by providing a convincing reference for genomic research.

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.

Telomere-to-telomere assembly of a fish Y chromosome reveals the origin of a young sex chromosome pair

BACKGROUND

The origin of sex chromosomes requires the establishment of recombination suppression between the proto-sex chromosomes. In many fish species, the sex chromosome pair is homomorphic with a recent origin, providing species for studying how and why recombination suppression evolved in the initial stages of sex chromosome differentiation, but this requires accurate sequence assembly of the X and Y (or Z and W) chromosomes, which may be difficult if they are recently diverged.

RESULTS

Here we produce a haplotype-resolved genome assembly of zig-zag eel (Mastacembelus armatus), an aquaculture fish, at the chromosomal scale. The diploid assembly is nearly gap-free, and in most chromosomes, we resolve the centromeric and subtelomeric heterochromatic sequences. In particular, the Y chromosome, including its highly repetitive short arm, has zero gaps. Using resequencing data, we identify a ~7 Mb fully sex-linked region (SLR), spanning the sex chromosome centromere and almost entirely embedded in the pericentromeric heterochromatin. The SLRs on the X and Y chromosomes are almost identical in sequence and gene content, but both are repetitive and heterochromatic, consistent with zero or low recombination. We further identify an HMG-domain containing gene HMGN6 in the SLR as a candidate sex-determining gene that is expressed at the onset of testis development.

CONCLUSIONS

Our study supports the idea that preexisting regions of low recombination, such as pericentromeric regions, can give rise to SLR in the absence of structural variations between the proto-sex chromosomes.

Improved contiguity of the threespine stickleback genome using long-read sequencing

While the cost and time for assembling a genome has drastically decreased, it still remains a challenge to assemble a highly contiguous genome. These challenges are rapidly being overcome by the integration of long-read sequencing technologies. Here, we use long-read sequencing to improve the contiguity of the threespine stickleback fish (Gasterosteus aculeatus) genome, a prominent genetic model species. Using Pacific Biosciences sequencing, we assembled a highly contiguous genome of a freshwater fish from Paxton Lake. Using contigs from this genome, we were able to fill over 76.7% of the gaps in the existing reference genome assembly, improving contiguity over fivefold. Our gap filling approach was highly accurate, validated by 10X Genomics long-distance linked-reads. In addition to closing a majority of gaps, we were able to assemble segments of telomeres and centromeres throughout the genome. This highlights the power of using long sequencing reads to assemble highly repetitive and difficult to assemble regions of genomes. This latest genome build has been released through a newly designed community genome browser that aims to consolidate the growing number of genomics datasets available for the threespine stickleback fish.

A streamlined workflow for conversion, peer review, and publication of genomics metadata as omics data papers

Background

Data papers have emerged as a powerful instrument for open data publishing, obtaining credit, and establishing priority for datasets generated in scientific experiments. Academic publishing improves data and metadata quality through peer review and increases the impact of datasets by enhancing their visibility, accessibility, and reusability.

Objective

We aimed to establish a new type of article structure and template for omics studies: the omics data paper. To improve data interoperability and further incentivize researchers to publish well-described datasets, we created a prototype workflow for streamlined import of genomics metadata from the European Nucleotide Archive directly into a data paper manuscript.

Methods

An omics data paper template was designed by defining key article sections that encourage the description of omics datasets and methodologies. A metadata import workflow, based on REpresentational State Transfer services and Xpath, was prototyped to extract information from the European Nucleotide Archive, ArrayExpress, and BioSamples databases.

Findings

The template and workflow for automatic import of standard-compliant metadata into an omics data paper manuscript provide a mechanism for enhancing existing metadata through publishing.

Conclusion

The omics data paper structure and workflow for import of genomics metadata will help to bring genomic and other omics datasets into the spotlight. Promoting enhanced metadata descriptions and enforcing manuscript peer review and data auditing of the underlying datasets brings additional quality to datasets. We hope that streamlined metadata reuse for scholarly publishing encourages authors to create enhanced metadata descriptions in the form of data papers to improve both the quality of their metadata and its findability and accessibility.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

Chromosome-level genome assembly of the humpback puffer, Tetraodon palembangensis

The humpback puffer, Tetraodon palembangensis, is a poisonous freshwater pufferfish species mainly distributed in Southeast Asia (Thailand, Laos, Malaysia and Indonesia). The humpback puffer has many interesting biological features, such as inactivity, tetrodotoxin production and body expansion. Here, we report the first chromosome-level genome assembly of the humpback puffer. The genome size is 362 Mb, with a contig N50 value of ∼1.78 Mb and a scaffold N50 value of ∼15.8 Mb. Based on this genome assembly, ∼61.5 Mb (18.11%) repeat sequences were identified, 19,925 genes were annotated, and the function of 90.01% of these genes could be predicted. Finally, a phylogenetic tree of ten teleost fish species was constructed. This analysis suggests that the humpback puffer and T. nigroviridis share a common ancestor 18.1 million years ago (MYA), and diverged from T. rubripes 45.8 MYA. The humpback puffer genome will be a valuable genomic resource to illustrate possible mechanisms of tetrodotoxin synthesis and tolerance.

Molecular Characterization of the Von Willebrand Factor Type D Domain of Vitellogenin from Takifugu flavidus

The von Willebrand factor type D (VWD) domain in vitellogenin has recently been found to bind tetrodotoxin. The way in which this protein domain associates with tetrodotoxin and participates in transporting tetrodotoxin in vivo remains unclear. A cDNA fragment of the vitellogenin gene containing the VWD domain from pufferfish (Takifugu flavidus) (TfVWD) was cloned. Using in silico structural and docking analyses of the predicted protein, we determined that key amino acids (namely, Val¹¹⁵, ASP¹¹⁶, Val¹¹⁷, and Lys¹²²) in TfVWD mediate its binding to tetrodotoxin, which was supported by in vitro surface plasmon resonance analysis. Moreover, incubating recombinant rTfVWD together with tetrodotoxin attenuated its toxicity in vivo, further supporting protein–toxin binding and indicating associated toxicity-neutralizing effects. Finally, the expression profiling of TfVWD across different tissues and developmental stages indicated that its distribution patterns mirrored those of tetrodotoxin, suggesting that TfVWD may be involved in tetrodotoxin transport in pufferfish. For the first time, this study reveals the amino acids that mediate the binding of TfVWD to tetrodotoxin and provides a basis for further exploration of the molecular mechanisms underlying the enrichment and transfer of tetrodotoxin in pufferfish.

Genome-Scale Phylogenetic and Population Genetic Studies Provide Insight Into Introgression and Adaptive Evolution of Takifugu Species in East Asia

As a typical marine adaptive radiation species, most Takifugu species are widely distributed in East Asian offshore, which have diversified morphological characteristics and different ecological habits. The phylogenetic relationship and population structure of the Takifugu species was complicated because of incomplete lineage sorting, widespread hybridization and introgression. Therefore, to systematically clarify the phylogenetic relationships of Takifugu genus, explore the introgression and natural hybridization between different Takifugu species, and detect the selective signatures in the adaptive evolution of diversified traits, whole-genome resequencing was used in 122 Takifugu samples from 10 species. Phylogenetic analysis showed solid sister-group relationships between Takifugu bimaculatus and Takifugu flavidus, Takifugu oblongus, and Takifugu niphobles, Takifugu rubripes, and Takifugu obscurus, Takifugu xanthoptreus, and Takifugu ocellatus. Further admixture analysis indicated the divergence of T. obscurus population and the bidirectional gene flow between T. bimaculatus and T. flavidus. Using species-specific homozygous genetic variance sites, we detected the asymmetric introgression between T. bimaculatus and T. flavidus at China East sea and southern Taiwan Strait. By genome-scale genetic diversity scanning, we detected two copies of syt1, zar1 and tgfbr1 related to the semilunar reproduction rhythm in T. niphobles, involved in memory formation, embryo maturation and female reproduction. Furthermore, we also found lots of T. niphobles specific mutations in CDS region of circadian rhythm related genes and endocrine hormone genes. For Takifugu species, our research provides reliable genetic resources and results for the phylogeny, introgression, hybridization and adaptive evolution, and could be used as a guide for the formulation of the protection and proliferation release policies.

Sequencing smart: De novo sequencing and assembly approaches for a non-model mammal

Background

Whilst much sequencing effort has focused on key mammalian model organisms such as mouse and human, little is known about the relationship between genome sequencing techniques for non-model mammals and genome assembly quality. This is especially relevant to non-model mammals, where the samples to be sequenced are often degraded and of low quality. A key aspect when planning a genome project is the choice of sequencing data to generate. This decision is driven by several factors, including the biological questions being asked, the quality of DNA available, and the availability of funds. Cutting-edge sequencing technologies now make it possible to achieve highly contiguous, chromosome-level genome assemblies, but rely on high-quality high molecular weight DNA. However, funding is often insufficient for many independent research groups to use these techniques. Here we use a range of different genomic technologies generated from a roadkill European polecat (Mustela putorius) to assess various assembly techniques on this low-quality sample. We evaluated different approaches for de novo assemblies and discuss their value in relation to biological analyses.

Results

Generally, assemblies containing more data types achieved better scores in our ranking system. However, when accounting for misassemblies, this was not always the case for Bionano and low-coverage 10x Genomics (for scaffolding only). We also find that the extra cost associated with combining multiple data types is not necessarily associated with better genome assemblies.

Conclusions

The high degree of variability between each de novo assembly method (assessed from the 7 key metrics) highlights the importance of carefully devising the sequencing strategy to be able to carry out the desired analysis. Adding more data to genome assemblies does not always result in better assemblies, so it is important to understand the nuances of genomic data integration explained here, in order to obtain cost-effective value for money when sequencing genomes.

Chromosomal-level assembly of Takifugu obscurus (Abe, 1949) genome using third-generation DNA sequencing and Hi-C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome-level reference genome for T. obscurus was constructed using long-read Pacific Biosciences (PacBio) Sequel sequencing and a Hi-C-based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single-molecule real-time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome-scale scaffolds using a Hi-C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one-to-one orthologous relationship was observed between the two species, supporting the chromosome-level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu-specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.