The round goby genome provides insights into mechanisms that may facilitate biological invasions

Reference genome for the endangered, genetically subdivided, northern tidewater goby, Eucyclogobius newberryi

The federally endangered sister species, Eucyclogobius newberryi (northern tidewater goby, NTG) and E. kristinae (southern tidewater goby) comprise the California endemic genus Eucyclogobius, which historically occurred in all coastal California counties. Isolated lagoons that only intermittently connect to the sea are their primary habitat. Reproduction occurs during lagoon closure, minimizing marine dispersal and generating the most genetically subdivided vertebrate genus on the California coast. We present a new genome assembly for E. newberryi using HiFi long reads and Hi-C chromatin-proximity sequencing. The 980 Mb E. newberryi reference genome has an N50 of 34 Mb with 22 well-described scaffolds comprising 88% of the genome and a complete BUSCO (Benchmarking Universal Single-Copy Orthologs) score of 96.7%. This genome will facilitate studies addressing selection, drift, and metapopulation genetics in subdivided populations, as well as the persistence of the critically endangered E. kristinae, where reintroduction will be an essential element of conservation actions for recovery. It also provides tools critical to the recovery of the genetically distinct management units in the NTG, as well as broader ecological and evolutionary studies of gobies, the most speciose family of fishes in the world.

Genomic investigations of successful invasions: the picture emerging from recent studies

Invasion biology aims to identify traits and mechanisms that contribute to successful invasions, while also providing general insights into the mechanisms underlying population expansion and adaptation to rapid climate and habitat changes. Certain phenotypic attributes have been linked to successful invasions, and the role of genetics has been critical in understanding adaptation of invasive species. Nevertheless, a comprehensive summary evaluating the most common evolutionary mechanisms associated with successful invasions across species and environments is still lacking. Here we present a systematic review of studies since 2015 that have applied genomic tools to investigate mechanisms of successful invasions across different organisms. We examine demographic patterns such as changes in genomic diversity at the population level, the presence of genetic bottlenecks and gene flow in the invasive range. We review mechanisms of adaptation such as selection from standing genetic variation and de novo mutations, hybridisation and introgression, all of which can have an impact on invasion success. This comprehensive review of recent articles on the genomic diversity of invasive species led to the creation of a searchable database to provide researchers with an accessible resource. Analysis of this database allowed quantitative assessment of demographic and adaptive mechanisms acting in invasive species. A predominant role of admixture in increasing levels of genetic diversity enabling molecular adaptation in novel habitats is the most important finding of our study. The "genetic paradox" of invasive species was not validated in genomic data across species and ecosystems. Even though the presence of genetic drift and bottlenecks is commonly reported upon invasion, a large reduction in genomic diversity is rarely observed. Any decrease in genetic diversity is often relatively mild and almost always restored via gene flow between different invasive populations. The fact that loci under selection are frequently detected suggests that adaptation to novel habitats on a molecular level is not hindered. The above findings are confirmed herein for the first time in a semi-quantitative manner by molecular data. We also point to gaps and potential improvements in the design of studies of mechanisms driving rapid molecular adaptation in invasive populations. These include the scarcity of comprehensive studies that include sampling from multiple native and invasive populations, identification of invasion sources, longitudinal population sampling, and the integration of fitness measures into genomic analyses. We also note that the potential of whole genome studies is often not exploited fully in predicting invasive potential. Comparative genomic studies identifying genome features promoting invasions are underrepresented despite their potential for use as a tool in invasive species control.

Does swimming at the bottom serve as a hydraulic advantage for benthic fish Neogobius melanostomus Pallas (1814) in flowing water?

Benthic fish, such as the round goby (Neogobius melanostomus Pallas, 1814) tend to swim near the bottom, especially at increased water velocities. To test whether these fish have a hydraulic advantage from swimming near the bottom and how the substrate affects the forces experienced, we measured the hydraulic forces experienced by preserved fish in a flow channel. The fish were tested 5.0 mm above the bottom at smooth and rough surface, and in the water column (10.0 cm elevation) above smooth and rough surface at 0.95 m/s water velocity. No significant effect among the mean hydraulic forces was observed between both fish positions, whereas the mean hydraulic forces in the water column were significantly higher (P<0.05) above the rough surface (mean 0.077 N±0.025 s.d.) than above the smooth surface (mean 0.068 N±0.021 s.d.). A convolutional neural network (CNN) predicted the column smooth treatment was the most characteristic force data time series (mean F1=0.88±0.03 s.d.). We conclude that the body posture and body movements of the fish are more relevant for the hydraulic forces experienced by the fish than the vertical position in the water column. Further factors explaining the affinity to swimming near the bottom are discussed.

Round gobies (Neogobius melanostomus) in the River Rhine: Population genetic support for invasion via two different routes

The round goby, Neogobius melanostomus, is a successful invasive fish species. Originating from the Caspian and Black Sea, it is now distributed widely within European fresh- and brackish waters. The River Rhine was colonized in 2008 only a few years after the opening of the Rhine-Main-Danube canal and only four years after N. melanostomus was first reported in the upper Danube River. Its invasion history of the River Rhine was unclear because the species was first detected close to the Rhine river delta which would suggest a route of invasion other than via the Rhine-Main-Danube canal. To investigate the colonization history of N. melanostomus in the Rhine, we combined abundance estimates with molecular analysis. Abundance estimates found N. melanostomus to be dominant in the Lower Rhine (> 90% of all catches). Molecular analysis was done on 286 individuals from four different sites. Using the mitochondrial control region (d-loop), we found three different haplotypes in both Rhine sites. None of the potential invasive source populations in the rivers Danube and Trave exhibited all three haplotypes. The molecular data therefore supported a scenario of two different colonization directions. Our results show that the invasion history of the River Rhine is complex and warrants further investigation.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and noncoding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole-genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared with other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared with a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

Copy number variation and population-specific immune genes in the model vertebrate zebrafish

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Effects of CRISPR/Cas9 targeting of the myo-inositol biosynthesis pathway on hyper-osmotic tolerance of tilapia cells

Myo-inositol is an important compatible osmolyte in vertebrates. This osmolyte is produced by the myo-inositol biosynthesis (MIB) pathway composed of myo-inositol phosphate synthase and inositol monophosphatase. These enzymes are among the highest upregulated proteins in tissues and cell cultures from teleost fish exposed to hyperosmotic conditions indicating high importance of this pathway for tolerating this type of stress. CRISPR/Cas9 gene editing of tilapia cells produced knockout lines of MIB enzymes and control genes. Metabolic activity decreased significantly for MIB KO lines in hyperosmotic media. Trends of faster growth of the MIB knockout lines in isosmotic media and faster decline of MIB knockout lines in hyperosmotic media were also observed. These results indicate a decline in metabolic fitness but only moderate effects on cell survival when tilapia cells with disrupted MIB genes are exposed to hyperosmolality. Therefore MIB genes are required for full osmotolerance of tilapia cells.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and non-coding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared to other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared to a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

Genomic Characteristics of Okamejei kenojei and the Implications to Its Evolutionary Biology Study

Okamejei kenojei is an important economic species widely distributed in shallow coastal waters of the western North Pacific. In this study, the whole-genome survey analysis of O. kenojei was conducted to reveal its genomic characteristics. The genome size was estimated to be 2027.44 Mb, the repeat sequence content was 44.90%, and the heterozygous ratio was 1.04%. The mitochondrial genome excavated from the sequencing data was 16,974 bp, and it can form the closed circular molecule. The phylogenetic tree based on 13 protein-coding gene sequences supported the validity of Okamejei and assisted the conclusion that Raja porosa was the junior synonym of O. kenojei. Plenty of potential microsatellite loci were identified, and the distribution frequency was estimated to be approximately 236.3 SSRs per Mb. Among all motif types of microsatellites, the dinucleotide repeats were dominant (82.59%), followed by the trinucleotide repeats (8.05%), tetranucleotide repeats (5.80%), pentanucleotide repeats (2.83%), and hexanucleotide repeats (0.72%). The results of the present study could not only provide useful information for understanding the genome structure and functional characteristics of O. kenojei, but also lay the foundation for the subsequent mapping of the whole genome.

Genome size estimation and its associations with body length, chromosome number and evolution in teleost fishes

Precise estimation of genome size (GS) is vital for various genomic studies, such as deciding genome sequencing depth, genome assembly, biodiversity documentation, evolution, genetic disorders studies, duplication events etc. Animal Genome Size Database provides GS of over 2050 fish species, which ranges from 0.35 pg in pufferfish (Tetraodon nigroviridis) to 132.83 pg in marbled lungfish (Protopterus aethiopicus). The GS of majority of the fishes inhabiting waters of Indian subcontinent are still missing. In present study, we estimated GS of 51 freshwater teleost (31 commercially important, 7 vulnerable and 13 ornamental species) that ranged from 0.58 pg in banded gourami (Trichogaster fasciata) to 1.92 pg in scribbled goby (Awaous grammepomus). Substantial variation in GS was observed within the same fish orders (0.64-1.45 pg in cypriniformes, 0.70-1.41 pg in siluriformes and 0.58-1.92 pg in perciformes). We examined the relationship between the GS, chromosome number and body length across all the fishes. Body length was found to be associated with GS, whereas no relationship was noticed between the GS and the chromosome number. The analysis using ancestral information revealed haploid chromosome number 25, 27 and 24 for the most recent common ancestor of cypriniformes, siluriformes and perciformes, respectively. The study led to generation of new records on GS of 43 fish species and revalidated records for 8 species. The finding is valuable resource for further research in the areas of fish genomics, molecular ecology and evolutionary conservation genetics.

Ten years of marine evolutionary biology-Challenges and achievements of a multidisciplinary research initiative

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology.

The gastric proton pump in gobiid and mudskipper fishes. Evidence of stomach loss?

Stomach loss has occurred independently multiple times during gnathostome evolution with notable frequency within the Teleostei. Significantly, this loss of acid-peptic digestion has been found to correlate with the secondary genomic loss of the gastric proton pump subunits (atp4a, atp4b) and pepsinogens/pepsins (pga, pgc). Gastric glands produce gastric juice containing the acid and pepsin and thus their presence is a hallmark feature of a digestive system capable of acid-peptic digestion. However, in gobiid fishes although oesogaster and gastric glands have been identified histologically, their functional significance has been questioned. In the present study we address whether the gastric proton pump is present and expressed in gastric glands of the goby Neogobius species (Gobiidae) and in members of the family Oxudercidae, a group of amphibious gobiid fishes commonly known as mudskippers (genera: Periophthalmus, Boleophthalmus, Periophthalmodon and Scartelaos). We confirmed the presence of gastric glands and have immunohistochemically localized gastric proton pump expression to these glands in Neogobius fluviatilis and Periophthalmus novemradiatus, Periophthalmus barbarus and Boleophthalmus boddarti. Genome analysis in Neogobius melanostomus, Periophthalmus magnuspinnatus, Scartelaos histophorus, Boleophthalmus pectinirostris, and Periophthalmodon schlosseri revealed the presence of both atp4a and atp4b subunit orthologues in all species in a conserved genomic loci organization. Moreover, it was possible to deduce that the complete open reading frame and the key functional amino acid residues are present. The conserved expression of the gastric proton pump provides clear evidence of the potential for gastric acid secretion indicating that acid digestion is retained in these gobiid fishes and not lost.

Contrasting population differentiation in two sympatric Triplophysa loaches on the Qinghai-Tibet Plateau

Genetic differentiation in aquatic organisms is usually shaped by drainage connectivity. Sympatric aquatic species are thus expected to show similar population differentiation patterns and similar genetic responses to their habitats. Water bodies on the Qinghai-Tibet Plateau (QTP) have recently experienced dramatic physicochemical changes, threatening the biodiversity of aquatic organisms on the "roof of the world." To uncover ecological genetics in Tibetan loaches (Triplophysa)-the largest component of the QTP ichthyofauna-we characterized population differentiation patterns and adaptive mechanisms to salinity change in two sympatric and phylogenetically closely related Tibetan loaches, T. stewarti and T. stenura, by integrating population genomic, transcriptomic, and electron probe microanalysis approaches. Based on millions of genome-wide SNPs, the two Tibetan loach species show contrasting population differentiation patterns, with highly geographically structured and clear genetic differentiation among T. stewarti populations, whereas there is no such observation in T. stenura, which is also supported by otolith microchemistry mapping. While limited genetic signals of parallel adaption to salinity changes between the two species are found from either genetic or gene expression variation perspective, a catalog of genes involved in ion transport, energy metabolism, structural reorganization, immune response, detoxification, and signal transduction is identified to be related to adaptation to salinity change in Triplophysa loaches. Together, our findings broaden our understanding of the population characteristics and adaptive mechanisms in sympatric Tibetan loach species and would contribute to biodiversity conservation and management of aquatic organisms on the QTP.

Genomic data is missing for many highly invasive species, restricting our preparedness for escalating incursion rates

Biological invasions drive environmental change, potentially threatening native biodiversity, human health, and global economies. Population genomics is an increasingly popular tool in invasion biology, improving accuracy and providing new insights into the genetic factors that underpin invasion success compared to research based on a small number of genetic loci. We examine the extent to which population genomic resources, including reference genomes, have been used or are available for invasive species research. We find that 82% of species on the International Union for Conservation of Nature “100 Worst Invasive Alien Species” list have been studied using some form of population genetic data, but just 32% of these species have been studied using population genomic data. Further, 55% of the list’s species lack a reference genome. With incursion rates escalating globally, understanding how genome-driven processes facilitate invasion is critical, but despite a promising trend of increasing uptake, “invasion genomics” is still in its infancy. We discuss how population genomic data can enhance our understanding of biological invasion and inform proactive detection and management of invasive species, and we call for more research that specifically targets this area.

Recurrent expansions of B30.2-associated immune receptor families in fish

B30.2 domains, also known as PRY/SPRY, are key components of specific subsets of two large families of proteins involved in innate immunity: the tripartite motif proteins (TRIMs) and the Nod-like receptors (NLRs). TRIM proteins are important, often inducible factors of antiviral innate immunity, targeting multiple steps of viral cycles through a variety of mechanisms. NLRs prime and regulate systemic innate defenses, especially against bacteria, and control inflammation. Large TRIM and NLR subsets characterized by the presence of a B30.2 domain have been reported from a few fish species including zebrafish and seem to be strongly prone to gene duplication/expansion. Here, we performed a large-scale survey of these receptors across about 150 fish genomes, focusing on ray-finned fishes. We assessed the number and genomic distribution of domains and domain combinations associated with TRIMs, NLRs, and other genes containing B30.2 domains and looked for gene expansion patterns across fish groups. We then used a model to test the impact of taxonomy, genome size, and environmental variables on the copy numbers of these genes. Our findings reveal novel domain structures, clade-specific gains and losses. They also assist with the timing of the gene expansions, reveal patterns associated with the MHC, and lay the groundwork for further studies delving deeper into the forces that drive the copy number variation of immune genes on a species level.

Breathing space: deoxygenation of aquatic environments can drive differential ecological impacts across biological invasion stages

The influence of climate change on the ecological impacts of invasive alien species (IAS) remains understudied, with deoxygenation of aquatic environments often-overlooked as a consequence of climate change. Here, we therefore assessed how oxygen saturation affects the ecological impact of a predatory invasive fish, the Ponto-Caspian round goby (Neogobius melanostomus), relative to a co-occurring endangered European native analogue, the bullhead (Cottus gobio) experiencing decline in the presence of the IAS. In individual trials and mesocosms, we assessed the effect of high, medium and low (90%, 60% and 30%) oxygen saturation on: (1) functional responses (FRs) of the IAS and native, i.e. per capita feeding rates; (2) the impact on prey populations exerted; and (3) how combined impacts of both fishes change over invasion stages (Pre-invasion, Arrival, Replacement, Proliferation). Both species showed Type II potentially destabilising FRs, but at low oxygen saturation, the invader had a significantly higher feeding rate than the native. Relative Impact Potential, combining fish per capita effects and population abundances, revealed that low oxygen saturation exacerbates the high relative impact of the invader. The Relative Total Impact Potential (RTIP), modelling both consumer species’ impacts on prey populations in a system, was consistently higher at low oxygen saturation and especially high during invader Proliferation. In the mesocosm experiment, low oxygen lowered RTIP where both species were present, but again the IAS retained high relative impact during Replacement and Proliferation stages at low oxygen. We also found evidence of multiple predator effects, principally antagonism. We highlight the threat posed to native communities by IAS alongside climate-related stressors, but note that solutions may be available to remedy hypoxia and potentially mitigate impacts across invasion stages.

The Visual Opsin Gene Repertoires of Teleost Fishes: Evolution, Ecology, and Function

Visual opsin genes expressed in the rod and cone photoreceptor cells of the retina are core components of the visual sensory system of vertebrates. Here, we provide an overview of the dynamic evolution of visual opsin genes in the most species-rich group of vertebrates, teleost fishes. The examination of the rich genomic resources now available for this group reveals that fish genomes contain more copies of visual opsin genes than are present in the genomes of amphibians, reptiles, birds, and mammals. The expansion of opsin genes in fishes is due primarily to a combination of ancestral and lineage-specific gene duplications. Following their duplication, the visual opsin genes of fishes repeatedly diversified at the same key spectral-tuning sites, generating arrays of visual pigments sensitive from the ultraviolet to the red spectrum of the light. Species-specific opsin gene repertoires correlate strongly with underwater light habitats, ecology, and color-based sexual selection. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The genome of walking catfish Clarias magur (Hamilton, 1822) unveils the genetic basis that may have facilitated the development of environmental and terrestrial adaptation systems in air-breathing catfishes

The walking catfish Clarias magur (Hamilton, 1822) (magur) is an important catfish species inhabiting the Indian subcontinent. It is considered as a highly nutritious food fish and has the capability to walk to some distance, and survive a considerable period without water. Assembly, scaffolding and several rounds of iterations resulted in 3,484 scaffolds covering ∼94% of estimated genome with 9.88 Mb largest scaffold, and N50 1.31 Mb. The genome possessed 23,748 predicted protein encoding genes with annotation of 19,279 orthologous genes. A total of 166 orthologous groups represented by 222 genes were found to be unique for this species. The Computational Analysis of gene Family Evolution (CAFE) analysis revealed expansion of 207 gene families and 100 gene families have rapidly evolved. Genes specific to important environmental and terrestrial adaptation, viz. urea cycle, vision, locomotion, olfactory and vomeronasal receptors, immune system, anti-microbial properties, mucus, thermoregulation, osmoregulation, air-breathing, detoxification, etc. were identified and critically analysed. The analysis clearly indicated that C. magur genome possessed several unique and duplicate genes similar to that of terrestrial or amphibians’ counterparts in comparison to other teleostean species. The genome information will be useful in conservation genetics, not only for this species but will also be very helpful in such studies in other catfishes.

Whole-genome sequencing reveals sex determination and liver high-fat storage mechanisms of yellowstripe goby (Mugilogobius chulae)

As a promising novel marine fish model for future research on marine ecotoxicology as well as an animal model of human disease, the genome information of yellowstripe goby (Mugilogobius chulae) remains unknown. Here we report the first annotated chromosome-level reference genome assembly for yellowstripe goby. A 20.67-cM sex determination region was discovered on chromosome 5 and seven potential sex-determining genes were identified. Based on combined genome and transcriptome data, we identified three key lipid metabolic pathways for high-fat accumulation in the liver of yellowstripe goby. The changes in the expression patterns of MGLL and CPT1 at different development stage of the liver, and the expansion of the ABCA1 gene, innate immune gene TLR23, and TRIM family genes may help in balancing high-fat storage in hepatocytes and steatohepatitis. These results may provide insights into understanding the molecular mechanisms of sex determination and high-fat storage in the liver of marine fishes.

Immunoglobulin T genes in Actinopterygii

In teleost fishes, there are three immunoglobulin isotypes named immunoglobulin M (IgM), D (IgD), and T (IgT). IgT was the last to be described in teleost fishes, and it is specific to them. From recent fish genomes, we identified and studied the immunoglobulin heavy chain genes in Actinopterygii. For this analysis, a custom bioinformatics and machine learning pipeline, we call CHfinder, was developed that identifies the exons coding for the CH domains of fish immunoglobulins. Some IgT in teleost and holostean fish found from this systematic study have not been previously described. Phylogenetic analysis of the deduced amino acid sequences of the IgT CH1 exons reveals they are similar to the CH1 of IgM. This analysis also shows that the other three domains (CH2, CH3, and CH4) were not the result of recent IgM duplication processes in Actinopterygii, demonstrating that it is an immunoglobulin of earlier origin. The bioinformatics program, CHfinder, is publicly available at https://github.com/compimmuno/CHfinder.

De novo transcriptome assembly, functional annotation and characterization of the Atlantic bluefin tuna (Thunnus thynnus) larval stage

In the present work, we assembled and characterized a de novo larval transcriptome of the Atlantic bluefin tuna Thunnus thynnus by taking advantage of publicly available databases with the goal of better understanding its larval development. The assembled transcriptome comprised 37,117 protein-coding transcripts, of which 13,633 full-length (>80% coverage), with an Ex90N50 of 3061 bp and 76% of complete and single-copy core vertebrate genes orthologues. Of these transcripts, 34,980 had a hit against the EggNOG database and 14,983 with the KEGG database. Codon usage bias was identified in processes such as translation and muscle development. By comparing our data with a set of representative fish species, 87.1% of tuna transcripts were included in orthogroups with other species and 5.1% in assembly-specific orthogroups, which were enriched in terms related to muscle and bone development, visual system and ion transport. Following this comparative approach, protein families related to myosin, extracellular matrix and immune system resulted significantly expanded in the Atlantic bluefin tuna. Altogether, these results provide a glimpse of how the Atlantic bluefin tuna might have achieved early physical advantages over competing species in the pelagic environment. The information generated lays the foundation for future research on the more detailed exploration of physiological responses at the molecular level in different larval stages and paves the way to evolutionary studies on the Atlantic bluefin tuna.