Draft genome of the American Eel (Anguilla rostrata)

High-quality chromosome-level genome assembly of the whitespotted conger (Conger myriaster)

The whitespotted conger (Conger myriaster) is an ecologically and economically significant benthic marine species widely distributed across East Asia's coastal waters. Despite this importance, the genomic resources for this species remain limited, hindering evolutionary and aquaculture research. Here, we present the first high-quality chromosome-level genome assembly of C. myriaster using PacBio CLR, WGS, 10X Genomics and Hi-C data. The resulting 1.09 Gb genome assembly exhibits excellent contiguity, with 97.49% of sequences anchored onto 19 chromosomes. The assembled genome achieved a BUSCO completeness stands at 98.00%, containing 34.80% repetitive sequences and 24,063 predicted protein-coding genes. This foundational genomic resource overcomes a major limitation, providing the essential framework for future investigations into the evolutionary adaptations and for the genetic improvement of C. myriaster in aquaculture.

The First Genome-Wide Survey of Shortbelly Eel (Dysomma anguillare Barnard, 1923) to Provide Genomic Characteristics, Microsatellite Markers and Complete Mitogenome Information

Dysomma anguillare is a demersal eel widespread distributing in tropical waters of the Indo-West Pacific and Atlantic. As an important component of the coastal fishery and marine ecosystem, the lack of genomic information for this species severely restricts the progress of relevant researches. In this study, the abecedarian genome-wide characteristics and phylogenetic relationships analyses were carried out based on next-generation sequencing (NGS) platform. The revised genome size was approximately 1 310 Mb, with the largest scaffold length reaching 23 878 bp through K-mer (K = 17) method. The heterozygosity, repetitive rate and average GC content were about 0.94%, 51.93% and 42.23%, respectively. A total of 1 160 104 microsatellite motifs were identified from the de novo assembled genome of D. anguillare, in which dinucleotide repeats accounted for the largest proportion (592 234, 51.05%), the highest occurrence frequency (14.58%) as well as the largest relative abundance (379.27/Mb). The high-polymorphic and moderate-polymorphic loci composed around 73% of the total single sequence repeats (SSRs), showing a latent capacity for subsequent population genetic structure and genetic diversity appraisal researches. Another byproduct of whole-genome sequencing, the double-stranded and circular mitogenome (16 690 bp) was assembled to investigate the evolutionary relationships of D. anguillare. The phylogenic tree constructed with maximum likelihood (ML) method showed that D. anguillare was closely related to Synaphobranchidae species, and the molecular systematic results further supported classical taxonomy status of D. anguillare.

Panmixia in the American eel extends to its tropical range of distribution: Biological implications and policymaking challenges

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

Genomic microsatellite characteristics analysis of Dysommaanguillare (Anguilliformes, Dysommidae), based on high-throughput sequencing technology

Microsatellite loci were screened from the genomic data of Dysommaanguillare and their composition and distribution were analysed by bioinformatics for the first time. The results showed that 4,060,742 scaffolds with a total length of 1,562 Mb were obtained by high-throughput sequencing and 1,160,104 microsatellite loci were obtained by MISA screening, which were distributed on 770,294 scaffolds. The occurrence frequency and relative abundance were 28.57% and 743/Mb, respectively. Amongst the six complete microsatellite types, dinucleotide repeats accounted for the largest proportion (592,234, 51.05%), the highest occurrence frequency (14.58%) and the largest relative abundance (379.27/Mb). A total of 1488 microsatellite repeats were detected in the genome of D.anguillare, amongst which the hexanucleotide repeat motifs were the most abundant (608), followed by pentanucleotide repeat motifs (574), tetranucleotide repeat motifs (232), trinucleotide repeat motifs (59), dinucleotide repeat motifs (11) and mononucleotide repeat motifs (4). The abundance of microsatellites of the same repeat type decreased with the increase of copy numbers. Amongst the six types of nucleotide repeats, the preponderance of repeated motifs are A (191,390, 43.77%), CA (150,240, 25.37%), AAT (13,168, 14.05%), CACG (2,649, 8.14%), TAATG (119, 19.16%) and CCCTAA (190, 19.16%, 7.65%), respectively. The data of the number, distribution and abundance of different types of microsatellites in the genome of D.anguillare were obtained in this study, which would lay a foundation for the development of high-quality microsatellite markers of D.anguillare in the future.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

The digestive tract as an essential organ for water acquisition in marine teleosts: lessons from euryhaline eels

Adaptation to a hypertonic marine environment is one of the major topics in animal physiology research. Marine teleosts lose water osmotically from the gills and compensate for this loss by drinking surrounding seawater and absorbing water from the intestine. This situation is in contrast to that in mammals, which experience a net osmotic loss of water after drinking seawater. Water absorption in fishes is made possible by (1) removal of monovalent ions (desalinization) by the esophagus, (2) removal of divalent ions as carbonate (Mg/CaCO₃) precipitates promoted by HCO₃^- secretion, and (3) facilitation of NaCl and water absorption from diluted seawater by the intestine using a suite of unique transporters. As a result, 70-85% of ingested seawater is absorbed during its passage through the digestive tract. Thus, the digestive tract is an essential organ for marine teleost survival in the hypertonic seawater environment. The eel is a species that has been frequently used for osmoregulation research in laboratories worldwide. The eel possesses many advantages as an experimental animal for osmoregulation studies, one of which is its outstanding euryhalinity, which enables researchers to examine changes in the structure and function of the digestive tract after direct transfer from freshwater to seawater. In recent years, the molecular mechanisms of ion and water transport across epithelial cells (the transcellular route) and through tight junctions (the paracellular route) have been elucidated for the esophagus and intestine. Thanks to the rapid progress in analytical methods for genome databases on teleosts, including the eel, the molecular identities of transporters, channels, pumps and junctional proteins have been clarified at the isoform level. As 10 y have passed since the previous reviews on this subject, it seems relevant and timely to summarize recent progress in research on the molecular mechanisms of water and ion transport in the digestive tract in eels and to compare the mechanisms with those of other teleosts and mammals from comparative and evolutionary viewpoints. We also propose future directions for this research field to achieve integrative understanding of the role of the digestive tract in adaptation to seawater with regard to pathways/mechanisms including the paracellular route, divalent ion absorption, metabolon formation and cellular trafficking of transporters. Notably, some of these have already attracted practical attention in laboratories.

Phylogeny of teleost connexins reveals highly inconsistent intra- and interspecies use of nomenclature and misassemblies in recent teleost chromosome assemblies

BACKGROUND

Based on an initial collecting of database sequences from the gap junction protein gene family (also called connexin genes) in a few teleosts, the naming of these sequences appeared variable. The reasons could be (i) that the structure in this family is variable across teleosts, or (ii) unfortunate naming. Rather clear rules for the naming of genes in fish and mammals have been outlined by nomenclature committees, including the naming of orthologous and ohnologous genes. We therefore analyzed the connexin gene family in teleosts in more detail. We covered the range of divergence times in teleosts (eel, Atlantic herring, zebrafish, Atlantic cod, three-spined stickleback, Japanese pufferfish and spotted pufferfish; listed from early divergence to late divergence).

RESULTS

The gene family pattern of connexin genes is similar across the analyzed teleosts. However, (i) several nomenclature systems are used, (ii) specific orthologous groups contain genes that are named differently in different species, (iii) several distinct genes have the same name in a species, and (iv) some genes have incorrect names. The latter includes a human connexin pseudogene, claimed as GJA4P, but which in reality is Cx39.2P (a delta subfamily gene often called GJD2like). We point out the ohnologous pairs of genes in teleosts, and we suggest a more consistent nomenclature following the outlined rules from the nomenclature committees. We further show that connexin sequences can indicate some errors in two high-quality chromosome assemblies that became available very recently.

CONCLUSIONS

Minimal consistency exists in the present practice of naming teleost connexin genes. A consistent and unified nomenclature would be an advantage for future automatic annotations and would make various types of subsequent genetic analyses easier. Additionally, roughly 5% of the connexin sequences point out misassemblies in the new high-quality chromosome assemblies from herring and cod.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

RAD-Seq Reveals Patterns of Additive Polygenic Variation Caused by Spatially-Varying Selection in the American Eel (Anguilla rostrata)

The American Eel (Anguilla rostrata) has an exceptional life cycle characterized by panmictic reproduction at the species scale, random dispersal, and selection in a highly heterogeneous habitat extending from subtropical to subarctic latitudes. The genetic consequences of spatially-varying selection in this species have been investigated for decades, revealing subtle clines in allele frequency at a few loci that contrast with complete panmixia on the vast majority of the genome. Because reproduction homogenizes allele frequencies every generation, sampling size, and genomic coverage are critical to reach sufficient power to detect selected loci in this context. Here, we used a total of 710 individuals from 12 sites and 12,098 high-quality single nucleotide polymorphisms to re-evaluate the extent to which local selection affects the spatial distribution of genetic diversity in this species. We used environmental association methods to identify markers under spatially-varying selection, which indicated that selection affects ∼1.5% of the genome. We then evaluated the extent to which candidate markers collectively vary with environmental factors using additive polygenic scores. We found significant correlations between polygenic scores and latitude, longitude and temperature which are consistent with polygenic selection acting against maladapted genotypes in different habitats occupied by eels throughout their range of distribution. Gene functions associated with outlier markers were significantly enriched for the insulin signaling pathway, indicating that the trade-offs inherent to occupying such a large distribution range involve the regulation of metabolism. Overall, this study highlights the potential of the additive polygenic scores approach in detecting selective effects in a complex environment.

Genome Survey of the Freshwater Mussel Venustaconcha ellipsiformis (Bivalvia: Unionida) Using a Hybrid De Novo Assembly Approach

Freshwater mussels (Bivalvia: Unionida) serve an important role as aquatic ecosystem engineers but are one of the most critically imperilled groups of animals. Here, we used a combination of sequencing strategies to assemble and annotate a draft genome of Venustaconcha ellipsiformis, which will serve as a valuable genomic resource given the ecological value and unique “doubly uniparental inheritance” mode of mitochondrial DNA transmission of freshwater mussels. The genome described here was obtained by combining high-coverage short reads (65× genome coverage of Illumina paired-end and 11× genome coverage of mate-pairs sequences) with low-coverage Pacific Biosciences long reads (0.3× genome coverage). Briefly, the final scaffold assembly accounted for a total size of 1.54 Gb (366,926 scaffolds, N50 = 6.5 kb, with 2.3% of “N” nucleotides), representing 86% of the predicted genome size of 1.80 Gb, while over one third of the genome (37.5%) consisted of repeated elements and >85% of the core eukaryotic genes were recovered. Given the repeated genetic bottlenecks of V. ellipsiformis populations as a result of glaciations events, heterozygosity was also found to be remarkably low (0.6%), in contrast to most other sequenced bivalve species. Finally, we reassembled the full mitochondrial genome and found six polymorphic sites with respect to the previously published reference. This resource opens the way to comparative genomics studies to identify genes related to the unique adaptations of freshwater mussels and their distinctive mitochondrial inheritance mechanism.

Whole-Genome Sequencing of 84 Japanese Eels Reveals Evidence against Panmixia and Support for Sympatric Speciation

The Japanese eel (Anguilla japonica), European eel (Anguilla anguilla), and American eel (Anguilla rostrata) are migratory, catadromous, temperate zone fish sharing several common life cycle features. The population genetics of panmixia in these eel species has already been investigated. Our extensive population genetics analysis was based on 1400 Gb of whole-genome sequence (WGS) data from 84 eels. It demonstrated that a Japanese eel group from the Kuma River differed from other populations of the same species. Even after removing the potential adapted/selected single nucleotide polymorphism (SNP) data, and with very small differences (fixation index [Fst] = 0.01), we obtained results consistently indicating that panmixia does not occur in Japanese eels. The life cycle of the Japanese eel is well-established and the Kuma River is in the center of its habitat. Nevertheless, simple reproductive isolation is not the probable cause of non-panmixia in this species. We propose that the combination of spawning area subdivision, philopatry, and habitat preference/avoidance accounts for the non-panmixia in the Japanese eel population. We named this hypothesis the “reproductive isolation like subset mapping” (RISM) model. This finding may be indicative of the initial stages of sympatric speciation in these eels.

Genetic parameters and quantitative trait loci analysis associated with body size and timing at metamorphosis into glass eels in captive-bred Japanese eels (Anguilla japonica)

The Japanese eel (Anguilla japonica) is among the most important aquaculture fish species in Eastern Asia. The present study aimed to identify the genetic parameters underlying body size and the timing at metamorphosis from leptocephali to glass eels in captive-bred Japanese eels, with the intent to foster sustainable development. Larvae from a partly factorial cross (14 sires × 11 dams) were reared until the point of metamorphosis into glass eels. In these organisms, we observed moderate heritability and mild genetic correlations among traits related to body size (h² = 0.16–0.33) and timing at metamorphosis (h² = 0.36–0.41). In an F₁ full-sib family, quantitative trait loci (QTL) mapping for these traits identified one significant (genome-wide P < 0.05) and five suggestive QTLs (chromosome-wide P < 0.05). These results suggest that in the Japanese eel, metamorphic traits exhibit a polygenic genetic structure comprising many QTLs with small effects. In addition, we updated the genetic linkage map for the Japanese eel and integrated it with our newly constructed de novo genome assembly. The information and tools generated from this study will contribute to the development of freshwater eel genetics and genomics.

Tachykinin-3 Genes and Peptides Characterized in a Basal Teleost, the European Eel: Evolutionary Perspective and Pituitary Role

In mammals, neurokinin B (NKB) is a short peptide encoded by the gene tac3. It is involved in the brain control of reproduction by stimulating gonadotropin-releasing hormone (GnRH) neurons, mainly via kisspeptin. We investigated tac3 genes and peptides in a basal teleost, the European eel, which shows an atypical blockade of the sexual maturation at a prepubertal stage. Two tac3 paralogous genes (tac3a and tac3b) were identified in the eel genome, each encoding two peptides (NKBa or b and NKB-related peptide NKB-RPa or b). Amino acid sequence of eel NKBa is identical to human NKB, and the three others are novel peptide sequences. The four eel peptides present the characteristic C-terminal tachykinin sequence, as well as a similar alpha helix 3D structure. Tac3 genes were identified in silico in 52 species of vertebrates, and a phylogeny analysis was performed on the predicted TAC3 pre-pro-peptide sequences. A synteny analysis was also done to further assess the evolutionary history of tac3 genes. Duplicated tac3 genes in teleosts likely result from the teleost-specific whole genome duplication (3R). Among teleosts, TAC3b precursor sequences are more divergent than TAC3a, and a loss of tac3b gene would have even occurred in some teleost lineages. NKB-RP peptide, encoded beside NKB by tac3 gene in actinopterygians and basal sarcopterygians, would have been lost in ancestral amniotes. Tissue distribution of eel tac3a and tac3b mRNAs showed major expression of both transcripts in the brain especially in the diencephalon, as analyzed by specific qPCRs. Human NKB has been tested in vitro on primary culture of eel pituitary cells. Human NKB dose-dependently inhibited the expression of lhβ, while having no effect on other glycoprotein hormone subunits (fshβ, tshβ, and gpα) nor on gh. Human NKB also dose-dependently inhibited the expression of GnRH receptor (gnrh-r2). The four eel peptides have been synthesized and also tested in vitro. They all inhibited the expression of both lhβ and of gnrh-r2. This reveals a potential dual inhibitory role of the four peptides encoded by the two tac3 genes in eel reproduction, exerted at the pituitary level on both luteinizing hormone and GnRH receptor.

Rapid de novo assembly of the European eel genome from nanopore sequencing reads

We have sequenced the genome of the endangered European eel using the MinION by Oxford Nanopore, and assembled these data using a novel algorithm specifically designed for large eukaryotic genomes. For this 860 Mbp genome, the entire computational process takes two days on a single CPU. The resulting genome assembly significantly improves on a previous draft based on short reads only, both in terms of contiguity (N50 1.2 Mbp) and structural quality. This combination of affordable nanopore sequencing and light weight assembly promises to make high-quality genomic resources accessible for many non-model plants and animals.