Natural variation in SAR11 marine bacterioplankton genomes inferred from metagenomic data

Background  

One objective of metagenomics is to reconstruct information about specific uncultured organisms from fragmentary environmental DNA sequences. We used the genome of an isolate of the marine alphaproteobacterium SAR11 ('Candidatus Pelagibacter ubique'; strain HTCC1062), obtained from the cold, productive Oregon coast, as a query sequence to study variation in SAR11 metagenome sequence data from the Sargasso Sea, a warm, oligotrophic ocean gyre.     

Screening the Sargasso Sea metagenome for data to investigate genome evolution in Ostreococcus (Prasinophyceae, Chlorophyta)

The Sargasso Sea water shotgun sequencing unveiled an unprecedented glimpse of marine prokaryotic diversity and gene content. The sequence data was gathered from 0.8 microm filtered surface water extracts, and revealed picoeukaryotic (cell size<2 microm) sequences alongside the prokaryotic data. We used the available genome sequence of the picoeukaryote Ostreococcus tauri (Prasinophyceae, Chlorophyta) as a benchmark for the eukaryotic sequence content of the Sargasso Sea metagenome. Sequence data from at least two new Ostreococcus strains were identified and analyzed, and showed a bias towards higher coverage of the AT-rich organellar genomes. The Ostreococcus nuclear sequence data retrieved from the Sargasso metagenome is divided onto 731 scaffolds of average size 3917 bp, and covers 23% of the complete nuclear genome and 14% of the total number of protein coding genes in O. tauri. We used this environmental Ostreococcus sequence data to estimate the level of constraint on intronic and intergenic sequences in this compact genome.     

An analysis of the Sargasso Sea resource and the consequences for database composition

Background  

The environmental sequencing of the Sargasso Sea has introduced a huge new resource of genomic information. Unlike the protein sequences held in the current searchable databases, the Sargasso Sea sequences originate from a single marine environment and have been sequenced from species that are not easily obtainable by laboratory cultivation. The resource also contains very many fragments of whole protein sequences, a side effect of the shotgun sequencing method.     

The microbial selenoproteome of the Sargasso Sea

Background  

Selenocysteine (Sec) is a rare amino acid which occurs in proteins in major domains of life. It is encoded by TGA, which also serves as the signal for termination of translation, precluding identification of selenoprotein genes by available annotation tools. Information on full sets of selenoproteins (selenoproteomes) is essential for understanding the biology of selenium. Herein, we characterized the selenoproteome of the largest microbial sequence dataset, the Sargasso Sea environmental genome project.     

Sequence and Expression Analyses of Cytophaga-Like Hydrolases in a Western Arctic Metagenomic Library and the Sargasso Sea

Sequence analysis of environmental DNA promises to provide new insights into the ecology and biogeochemistry of uncultured marine microbes. In this study we used the Sargasso Sea Whole Genome Sequence (WGS) data set to search for hydrolases used by Cytophaga-like bacteria to degrade biopolymers such as polysaccharides and proteins. Analysis of the Sargasso WGS data for contigs bearing both the 16S rRNA genes of Cytophaga-like bacteria and hydrolase genes revealed a cellulase gene (celM) most similar to the gene found in Cytophaga hutchinsonii. A BLAST search of the entire Sargasso Sea WGS data set indicated that celM was the most abundant cellulase-like gene in the Sargasso Sea. However, the similarity between CelM-like cellulases and peptidases belonging to metalloprotease family M42 led us to question whether CelM is involved in the degradation of polysaccharides or proteins. PCR primers were designed for the celM genes in the Sargasso Sea WGS data set and used to identify celM in a fosmid library constructed with prokaryotic DNA from the western Arctic Ocean. Expression analysis of the Cytophaga-like Arctic CelM, which is 63% identical and 77% similar to CelM in C. hutchinsonii, indicated that there was peptidase activity, whereas cellulase activity was not detected. Our analysis suggests that the celM gene plays a role in the degradation of protein by Cytophaga-like bacteria. The abundance of peptidase genes in the Cytophaga-like fosmid clone provides further evidence for the importance of Cytophaga-like bacteria in the degradation of protein in high-molecular-weight dissolved organic matter.     

Swimming suppresses hepatic vitellogenesis in European female silver eels as shown by expression of the estrogen receptor 1, vitellogenin1 and vitellogenin2 in the liver

Background  

When European silver eels (Anguilla anguilla) venture into the Atlantic Ocean for their 6,000 km semelparous spawning run to the Sargasso Sea, they are still in a prepubertal stage. Further sexual development appears to be blocked by dopaminergic inhibition of hypothalamus and pituitary activity. Recently, we found that swimming for several weeks in freshwater stimulated the incorporation of fat droplets in the oocytes. So, it was hypothesized that long term swimming in seawater would release the inhibition further and would also stimulate the production of vitellogenin by the liver.     

Comparison of prokaryotic diversity at offshore oceanic locations reveals a different microbiota in the Mediterranean Sea

The bacterial and archaeal assemblages at two offshore sites located in polar (Greenland Sea; depth: 50 and 2000 m) and Mediterranean (Ionian Sea; depth 50 and 3000 m) waters were studied by PCR amplification and sequencing of the last 450-500 bp of the 16S rRNA gene. A total of 1621 sequences, together with alignable 16S rRNA gene fragments from the Sargasso Sea metagenome database, were analysed to ascertain variations associated with geographical location and depth. The Ionian 50 m sample appeared to be the most diverse and also had remarkable differences in terms of the prokaryotic groups retrieved; surprisingly, however, many similarities were found at the level of large-scale diversity between the Sargasso database fragments and the Greenland 50 m sample. Most sequences with more than 97% sequence similarity, a value often taken as indicative of species delimitation, were only found at a single location/depth; nevertheless, a few examples of cosmopolitan sequences were found in all samples. Depth was also an important factor and, although both deep-water samples had overall similarities, there were important differences that could be due to the warmer waters at depth of the Mediterranean Sea.     

Metagenomic analysis reveals the contribution of anaerobic methanotroph-1b in the oxidation of methane at the Ulleung Basin,East Sea of Korea

We have previously identified a sulfate methane transition zone (SMTZ) within the methane hydrate-bearing sediment in the Ulleung Basin, East Sea of Korea, and the presence of ANME-1b group in the sediment has been shown by phylogenetic analysis of a 16S rRNA gene. Herein, we describe taxonomic and functional profiling in the SMTZ sample by metagenomic analysis, comparing with that of surface sediment. Metagenomic sequences of 115 Mbp and 252 Mbp were obtained from SMTZ and surface sediments, respectively. The taxonomic profiling using BLASTX against the SEED within MG-RAST showed the prevalence of methanogens (19.1%), such as Methanosarcinales (12.0%) and Methanomicrobiales (4.1%) predominated within the SMTZ metagenome. A number of 185,200 SMTZ reads (38.9%) and 438,484 surface reads (62.5%) were assigned to functional categories, and methanogenesis-related reads were statistically significantly overrepresented in the SMTZ metagenome. However, the mapping analysis of metagenome reads to the reference genomes, most of the sequences of the SMTZ metagenome were mapped to ANME-1 draft genomes, rather than those of methanogens. Furthermore, the two copies of the methyl-coenzyme M reductase gene (mcrA) segments of the SMTZ metagenome were clustered with ANME-1b in the phylogenetic cluster. These results indicate that ANME-1b reads were miss-annotated to methanogens due to limitation of database. Many of key genes necessary for reverse methanogenesis were present in the SMTZ metagenome, except for N⁵,N¹⁰-methenyl-H₄MPT reductase (mer) and CoB-CoM heterodisulfide reductase subunits D and E (hdrDE). These data suggest that the ANME-1b represents the primary player the anaerobic methane oxidation in the SMTZ, of the methane hydrate-bearing sediment at the Ulleung Basin, East Sea of Korea.     

Distribution analysis of hydrogenases in surface waters of marine and freshwater environments

Background

Surface waters of aquatic environments have been shown to both evolve and consume hydrogen and the ocean is estimated to be the principal natural source. In some marine habitats, H₂ evolution and uptake are clearly due to biological activity, while contributions of abiotic sources must be considered in others. Until now the only known biological process involved in H₂ metabolism in marine environments is nitrogen fixation.

Principal Findings

We analyzed marine and freshwater environments for the presence and distribution of genes of all known hydrogenases, the enzymes involved in biological hydrogen turnover. The total genomes and the available marine metagenome datasets were searched for hydrogenase sequences. Furthermore, we isolated DNA from samples from the North Atlantic, Mediterranean Sea, North Sea, Baltic Sea, and two fresh water lakes and amplified and sequenced part of the gene encoding the bidirectional NAD(P)-linked hydrogenase. In 21% of all marine heterotrophic bacterial genomes from surface waters, one or several hydrogenase genes were found, with the membrane-bound H₂ uptake hydrogenase being the most widespread. A clear bias of hydrogenases to environments with terrestrial influence was found. This is exemplified by the cyanobacterial bidirectional NAD(P)-linked hydrogenase that was found in freshwater and coastal areas but not in the open ocean.

Significance

This study shows that hydrogenases are surprisingly abundant in marine environments. Due to its ecological distribution the primary function of the bidirectional NAD(P)-linked hydrogenase seems to be fermentative hydrogen evolution. Moreover, our data suggests that marine surface waters could be an interesting source of oxygen-resistant uptake hydrogenases. The respective genes occur in coastal as well as open ocean habitats and we presume that they are used as additional energy scavenging devices in otherwise nutrient limited environments. The membrane-bound H₂-evolving hydrogenases might be useful as marker for bacteria living inside of marine snow particles.     

The Distribution and Ecology of Ariosoma Balearicum (Congridae) Leptocephali in the Western North Atlantic

A total of 4589 leptocephali of the congrid eel, Ariosoma balearicum, were examined from 17 cruises in the western North Atlantic Ocean. Myomere counts made on 915 of these indicated there were two ranges of number of myomeres that appear to be associated with separate spawning populations. Those with the higher range (high count: 128–137) were consistently 70–100mm in length in the Sargasso Sea from February to April and 20–80mm in length in the northern Sargasso Sea and Gulf Stream from September to October. Those with the lower range (low count: 120–130) were rare in the northern and eastern Sargasso Sea where they had consistently greater lengths than high count leptocephali and were most abundant in the Florida Current and Providence Channel. The geographic distributions of size and myomere ranges in relation to hydrography provide strong support for the hypothesis that high count eels found along the South Atlantic Bight (SAB) migrate across the Florida Current to spawn in the northwest Sargasso Sea. This migratory pattern is similar to those of Anguilla rostrata and Conger oceanicus, which use the southern Sargasso Sea for development as larvae. However, the distribution of high count leptocephali suggests that they use the entire Sargasso Sea gyre as a development area as larvae before crossing the Florida Current and recruiting to the SAB. The low count eels inhabiting the Bahamas appear to spawn near the banks and their abundance in the Providence Channel and southwest Sargasso Sea suggests most are retained close to the Bahamas. These two distinct styles of spawning, distribution and recruitment of larvae are hypothesized to be related to the different hydrographic regimes of the two juvenile habitats and the resulting constraints on growth and recruitment of larvae. Vertebral and myomere counts reported from other areas suggest there are distinct populations in other regions of the North Atlantic Ocean.     

Sequence and expression analyses of Cytophaga-like hydrolases in a Western arctic metagenomic library and the Sargasso Sea

Sequence analysis of environmental DNA promises to provide new insights into the ecology and biogeochemistry of uncultured marine microbes. In this study we used the Sargasso Sea Whole Genome Sequence (WGS) data set to search for hydrolases used by Cytophaga-like bacteria to degrade biopolymers such as polysaccharides and proteins. Analysis of the Sargasso WGS data for contigs bearing both the 16S rRNA genes of Cytophaga-like bacteria and hydrolase genes revealed a cellulase gene (celM) most similar to the gene found in Cytophaga hutchinsonii. A BLAST search of the entire Sargasso Sea WGS data set indicated that celM was the most abundant cellulase-like gene in the Sargasso Sea. However, the similarity between CelM-like cellulases and peptidases belonging to metalloprotease family M42 led us to question whether CelM is involved in the degradation of polysaccharides or proteins. PCR primers were designed for the celM genes in the Sargasso Sea WGS data set and used to identify celM in a fosmid library constructed with prokaryotic DNA from the western Arctic Ocean. Expression analysis of the Cytophaga-like Arctic CelM, which is 63% identical and 77% similar to CelM in C. hutchinsonii, indicated that there was peptidase activity, whereas cellulase activity was not detected. Our analysis suggests that the celM gene plays a role in the degradation of protein by Cytophaga-like bacteria. The abundance of peptidase genes in the Cytophaga-like fosmid clone provides further evidence for the importance of Cytophaga-like bacteria in the degradation of protein in high-molecular-weight dissolved organic matter.     

Patterns of ecological specialization among microbial populations in the Red Sea and diverse oligotrophic marine environments

Large swaths of the nutrient‐poor surface ocean are dominated numerically by cyanobacteria (Prochlorococcus), cyanobacterial viruses (cyanophage), and alphaproteobacteria (SAR11). How these groups thrive in the diverse physicochemical environments of different oceanic regions remains poorly understood. Comparative metagenomics can reveal adaptive responses linked to ecosystem‐specific selective pressures. The Red Sea is well‐suited for studying adaptation of pelagic‐microbes, with salinities, temperatures, and light levels at the extreme end for the surface ocean, and low nutrient concentrations, yet no metagenomic studies have been done there. The Red Sea (high salinity, high light, low N and P) compares favorably with the Mediterranean Sea (high salinity, low P), Sargasso Sea (low P), and North Pacific Subtropical Gyre (high light, low N). We quantified the relative abundance of genetic functions among Prochlorococcus, cyanophage, and SAR11 from these four regions. Gene frequencies indicate selection for phosphorus acquisition (Mediterranean/Sargasso), DNA repair and high‐light responses (Red Sea/Pacific Prochlorococcus), and osmolyte C1 oxidation (Red Sea/Mediterranean SAR11). The unexpected connection between salinity‐dependent osmolyte production and SAR11 C1 metabolism represents a potentially major coevolutionary adaptation and biogeochemical flux. Among Prochlorococcus and cyanophage, genes enriched in specific environments had ecotype distributions similar to nonenriched genes, suggesting that inter‐ecotype gene transfer is not a major source of environment‐specific adaptation. Clustering of metagenomes using gene frequencies shows similarities in populations (Red Sea with Pacific, Mediterranean with Sargasso) that belie their geographic distances. Taken together, the genetic functions enriched in specific environments indicate competitive strategies for maintaining carrying capacity in the face of physical stressors and low nutrient availability.     

CLAME: a new alignment-based binning algorithm allows the genomic description of a novel Xanthomonadaceae from the Colombian Andes

Background

Hot spring bacteria have unique biological adaptations to survive the extreme conditions of these environments; these bacteria produce thermostable enzymes that can be used in biotechnological and industrial applications. However, sequencing these bacteria is complex, since it is not possible to culture them. As an alternative, genome shotgun sequencing of whole microbial communities can be used. The problem is that the classification of sequences within a metagenomic dataset is very challenging particularly when they include unknown microorganisms since they lack genomic reference. We failed to recover a bacterium genome from a hot spring metagenome using the available software tools, so we develop a new tool that allowed us to recover most of this genome.

Results

We present a proteobacteria draft genome reconstructed from a Colombian’s Andes hot spring metagenome. The genome seems to be from a new lineage within the family Rhodanobacteraceae of the class Gammaproteobacteria, closely related to the genus Dokdonella. We were able to generate this genome thanks to CLAME. CLAME, from Spanish “CLAsificador MEtagenomico”, is a tool to group reads in bins. We show that most reads from each bin belong to a single chromosome. CLAME is very effective recovering most of the reads belonging to the predominant species within a metagenome.

Conclusions

We developed a tool that can be used to extract genomes (or parts of them) from a complex metagenome.

     

DNA barcoding from the hooked squids (Cephalopods: Onychoteuthidae) of the Sargasso Sea

Squids of the family Onychoteuthidae are ecologically important in pelagic food webs and have been reported from every ocean except the Arctic. Although they are abundant and caught frequently as bycatch in fisheries, the biogeography of many species remains poorly understood. Species identification within the Atlantic Ocean is usually restricted to two species: Onychoteuthis banksii and Onykia carriboea. Here, we report the occurrence of four species of the family Onychoteuthidae (Onychoteuthis cf. banksii, Onykia carriboea, Walvisteuthis jeremiahi, and Onychoteuthis sp. AL 2) from the Sargasso Sea in the western Atlantic, identified using DNA barcoding (cytochrome c oxidase subunit I) and morphology. Our results have doubled the known onychoteuthid biodiversity in the Sargasso Sea, which has implications for the ecology of this oceanic region.     

MetCap: a bioinformatics probe design pipeline for large-scale targeted metagenomics

BackgroundMassive sequencing of genes from different environments has evolved metagenomics as central to enhancing the understanding of the wide diversity of micro-organisms and their roles in driving ecological processes. Reduced cost and high throughput sequencing has made large-scale projects achievable to a wider group of researchers, though complete metagenome sequencing is still a daunting task in terms of sequencing as well as the downstream bioinformatics analyses. Alternative approaches such as targeted amplicon sequencing requires custom PCR primer generation, and is not scalable to thousands of genes or gene families.ResultsIn this study, we are presenting a web-based tool called MetCap that circumvents the limitations of amplicon sequencing of multiple genes by designing probes that are suitable for large-scale targeted metagenomics sequencing studies. MetCap provides a novel approach to target thousands of genes and genomic regions that could be used in targeted metagenomics studies. Automatic analysis of user-defined sequences is performed, and probes specifically designed for metagenome studies are generated. To illustrate the advantage of a targeted metagenome approach, we have generated more than 300,000 probes that match more than 400,000 publicly available sequences related to carbon degradation, and used these probes for target sequencing in a soil metagenome study. The results show high enrichment of target genes and a successful capturing of the majority of gene families. MetCap is freely available to users from: http://soilecology.biol.lu.se/metcap/.ConclusionMetCap is facilitating probe-based target enrichment as an easy and efficient alternative tool compared to complex primer-based enrichment for large-scale investigations of metagenomes. Our results have shown efficient large-scale target enrichment through MetCap-designed probes for a soil metagenome. The web service is suitable for any targeted metagenomics project that aims to study several genes simultaneously. The novel bioinformatics approach taken by the web service will enable researchers in microbial ecology to tap into the vast diversity of microbial communities using targeted metagenomics as a cost-effective alternative to whole metagenome sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0501-8) contains supplementary material, which is available to authorized users.     

MetaSim: a sequencing simulator for genomics and metagenomics

Background

The new research field of metagenomics is providing exciting insights into various, previously unclassified ecological systems. Next-generation sequencing technologies are producing a rapid increase of environmental data in public databases. There is great need for specialized software solutions and statistical methods for dealing with complex metagenome data sets.

Methodology/Principal Findings

To facilitate the development and improvement of metagenomic tools and the planning of metagenomic projects, we introduce a sequencing simulator called MetaSim. Our software can be used to generate collections of synthetic reads that reflect the diverse taxonomical composition of typical metagenome data sets. Based on a database of given genomes, the program allows the user to design a metagenome by specifying the number of genomes present at different levels of the NCBI taxonomy, and then to collect reads from the metagenome using a simulation of a number of different sequencing technologies. A population sampler optionally produces evolved sequences based on source genomes and a given evolutionary tree.

Conclusions/Significance

MetaSim allows the user to simulate individual read datasets that can be used as standardized test scenarios for planning sequencing projects or for benchmarking metagenomic software.     

The mysterious ecosystem at the ocean’s surface

Life on the ocean’s surface connects worlds. From shallow waters to the deep sea, the open ocean to rivers and lakes, numerous terrestrial and marine species depend on the surface ecosystem and the organisms found therein. Organisms that live freely at the surface, termed “neuston,” include keystone organisms like the golden seaweed Sargassum that makes up the Sargasso Sea, floating barnacles, snails, nudibranchs, and cnidarians. Many ecologically and economically important fish species live as or rely upon neuston. Species at the surface are not distributed uniformly; the ocean’s surface harbors unique neustonic communities and ecoregions found at only certain latitudes and only in specific ocean basins. But the surface is also on the front line of climate change and pollution. Despite the diversity and importance of the ocean’s surface in connecting disparate habitats, and the risks it faces, we know very little about neustonic life. This Essay will introduce you to the neuston, their connections to diverse habitats, the threats they face, and new opportunities for research and discovery at the air-sea interface.

The mysterious ’neuston’ ecosystem at the ocean’s surface includes keystone organisms like the golden seaweed Sargassum that makes up the Sargasso Sea, floating barnacles, snails, nudibranchs, and cnidarians; this Essay explores threats to its wellbeing and the importance of further research.     

NeSSM: A Next-Generation Sequencing Simulator for Metagenomics

Background

Metagenomics can reveal the vast majority of microbes that have been missed by traditional cultivation-based methods. Due to its extremely wide range of application areas, fast metagenome sequencing simulation systems with high fidelity are in great demand to facilitate the development and comparison of metagenomics analysis tools.

Results

We present here a customizable metagenome simulation system: NeSSM (Next-generation Sequencing Simulator for Metagenomics). Combining complete genomes currently available, a community composition table, and sequencing parameters, it can simulate metagenome sequencing better than existing systems. Sequencing error models based on the explicit distribution of errors at each base and sequencing coverage bias are incorporated in the simulation. In order to improve the fidelity of simulation, tools are provided by NeSSM to estimate the sequencing error models, sequencing coverage bias and the community composition directly from existing metagenome sequencing data. Currently, NeSSM supports single-end and pair-end sequencing for both 454 and Illumina platforms. In addition, a GPU (graphics processing units) version of NeSSM is also developed to accelerate the simulation. By comparing the simulated sequencing data from NeSSM with experimental metagenome sequencing data, we have demonstrated that NeSSM performs better in many aspects than existing popular metagenome simulators, such as MetaSim, GemSIM and Grinder. The GPU version of NeSSM is more than one-order of magnitude faster than MetaSim.

Conclusions

NeSSM is a fast simulation system for high-throughput metagenome sequencing. It can be helpful to develop tools and evaluate strategies for metagenomics analysis and it’s freely available for academic users at http://cbb.sjtu.edu.cn/~ccwei/pub/software/NeSSM.php.     

Spawning ofConger oceanicus andConger triporiceps (Congridae) in the Sargasso Sea and subsequent distribution of leptocephali

Synopsis Distribution of leptocephali ofConger in the Western North Atlantic Ocean was studied using specimens from our collections, specimens from other collections, and various existing collection records. The presence of leptocephali ofConger oceanicus andConger triporiceps < 30 mm long over deep water in the southwestern Sargasso Sea in autumn and winter implies a protracted spawning period there. The subtropical convergence zone, meandering east-west across the Sargasso Sea, is probably the northern limit of spawning of both species. Spawning may also occur close to the Bahamas and Antilles.C. triporiceps may spawn also in the Caribbean Sea judging by the capture of small leptocephali in the western Caribbean and of the more southerly continental distribution of its juveniles. The claim of Johannes Schmidt in 1931 that the EuropeanC. conger spawns across the North Atlantic into the western Sargasso Sea is probably incorrect, because leptocephali ofConger are rare in the eastern Sargasso Sea and becauseC. triporiceps, with myomere numbers overlapping those ofC. conger, was recently described in the western North Atlantic. With increasing size, leptocephali ofC. oceanicus and a portion ofC. triporiceps spread westward and northward in the Florida Current and Gulf Stream, but larger leptocephali especially ofC. triporiceps are found also in the Caribbean and Gulf of Mexico. Spawning ofC. oceanicus in the Sargasso Sea indicates that adults cross the Florida Current-Gulf Stream, and successful leptocephali cross the current in the opposite direction to colonize juvenile habitat on the continental shelf, a migratory pattern similar to that of the American eelAnguilla rostrata (Anguillidae).     

Qualitative assessment of the diet of European eel larvae in the Sargasso Sea resolved by DNA barcoding

European eels (Anguilla anguilla) undertake spawning migrations of more than 5000 km from continental Europe and North Africa to frontal zones in the Sargasso Sea. Subsequently, the larval offspring are advected by large-scale eastward ocean currents towards continental waters. However, the Sargasso Sea is oligotrophic, with generally low plankton biomass, and the feeding biology of eel larvae has so far remained a mystery, hampering understanding of this peculiar life history. DNA barcoding of gut contents of 61 genetically identified A. anguilla larvae caught in the Sargasso Sea showed that even the smallest larvae feed on a striking variety of plankton organisms, and that gelatinous zooplankton is of fundamental dietary importance. Hence, the specific plankton composition seems essential for eel larval feeding and growth, suggesting a linkage between eel survival and regional plankton productivity. These novel insights into the prey of Atlantic eels may furthermore facilitate eel larval rearing in aquaculture, which ultimately may replace the unsustainable use of wild-caught glass eels.