Comprehensive phylogeny of ray-finned fishes (Actinopterygii) based on transcriptomic and genomic data

The Fish Family Poeciliidae as a Model to Study the Evolution and Diversification of Regenerative Capacity in Vertebrates

The capacity of regenerating a new structure after losing an old one is a major challenge in the animal kingdom. Fish have emerged as an interesting model to study regeneration due to their high and diverse regenerative capacity. To date, most efforts have focused on revealing the mechanisms underlying fin regeneration, but information on why and how this capacity evolves remains incomplete. Here, we propose the livebearing fish family Poeciliidae as a promising new model system to study the evolution of fin regeneration. First, we review the current state of knowledge on the evolution of regeneration in the animal kingdom, with a special emphasis on fish fins. Second, we summarize recent advances in our understanding of the mechanisms behind fin regeneration in fish. Third, we discuss potential evolutionary pressures that may modulate the regenerative capacity of fish fins and propose three new theories for how natural and sexual selection can lead to the evolution of fin regeneration: (1) signaling-driven fin regeneration, (2) predation-driven fin regeneration, and (3) matrotrophy-suppressed fin regeneration. Finally, we argue that fish from the family Poeciliidae are an excellent model system to test these theories, because they comprise of a large variety of species in a well-defined phylogenetic framework that inhabit very different environments and display remarkable variation in reproductive traits, allowing for comparative studies of fin regeneration among closely related species, among populations within species or among individuals within populations. This new model system has the potential to shed new light on the underlying genetic and molecular mechanisms driving the evolution and diversification of regeneration in vertebrates.

The SARS-CoV-2 receptor and other key components of the Renin-Angiotensin-Aldosterone System related to COVID-19 are expressed in enterocytes in larval zebrafish

People with underlying conditions, including hypertension, obesity, and diabetes, are especially susceptible to negative outcomes after infection with coronavirus SARS-CoV-2, which causes COVID-19. Hypertension and respiratory inflammation are exacerbated by the Renin-Angiotensin-Aldosterone System (RAAS), which normally protects from rapidly dropping blood pressure via Angiotensin II (Ang II) produced by the enzyme Ace. The Ace paralog Ace2 degrades Ang II, counteracting its chronic effects, and serves as the SARS-CoV-2 receptor. Ace, the coronavirus, and COVID-19 comorbidities all regulate Ace2, but we do not yet understand how. To exploit zebrafish (Danio rerio) to help understand the relationship of the RAAS to COVID-19, we must identify zebrafish orthologs and co-orthologs of human RAAS genes and understand their expression patterns. To achieve these goals, we conducted genomic and phylogenetic analyses and investigated single cell transcriptomes. Results showed that most human RAAS genes have one or more zebrafish orthologs or co-orthologs. Results identified a specific type of enterocyte as the specific site of expression of zebrafish orthologs of key RAAS components, including Ace, Ace2, Slc6a19 (SARS-CoV-2 co-receptor), and the Angiotensin-related peptide cleaving enzymes Anpep (receptor for the common cold coronavirus HCoV-229E), and Dpp4 (receptor for the Middle East Respiratory Syndrome virus, MERS-CoV). Results identified specific vascular cell subtypes expressing Ang II receptors, apelin, and apelin receptor genes. These results identify genes and cell types to exploit zebrafish as a disease model for understanding mechanisms of COVID-19.

MRAP2 Interaction with Melanocortin-4 Receptor in SnakeHead (Channa argus)

The melanocortin-4 receptor (MC4R) plays an important role in the regulation of food intake and energy expenditure. Melanocortin-2 receptor accessory protein 2 (MRAP2) modulates trafficking, ligand binding, and signaling of MC4R. The Northern snakehead (Channa argus) is an economically important freshwater fish native to East Asia. To explore potential interaction between snakehead MC4R and MRAP2, herein we cloned snakehead mc4r and mrap2. The snakehead mc4r consisted of a 984 bp open reading frame encoding a protein of 327 amino acids, while snakehead mrap2 contained a 693 bp open reading frame encoding a protein of 230 amino acids. Synteny analysis indicated that mc4r was highly conserved with similar gene arrangement, while mrap2 contained two isoforms in teleost with different gene orders. Snakehead mc4r was primarily expressed in the brain, whereas mrap2 was expressed in the brain and intestine. Snakehead mc4r and mrap2 expression was modulated by fasting and refeeding. Further pharmacological experiments showed that the cloned snakehead MC4R was functional, capable of binding to peptide agonists and increasing intracellular cAMP production in a dose-dependent manner. Snakehead MC4R exhibited high constitutive activity. MRAP2 significantly decreased basal and agonist-stimulated cAMP signaling. These findings suggest that snakehead MC4R might be involved in energy balance regulation by interacting with MRAP2. Further studies are needed to elucidate MC4R in regulating diverse physiological processes in snakehead.

From scales to armor: Scale losses and trunk bony plate gains in ray-finned fishes

Actinopterygians (ray‐finned fishes) are the most diversified group of vertebrates and are characterized by a variety of protective structures covering their integument, the evolution of which has intrigued biologists for decades. Paleontological records showed that the first mineralized vertebrate skeleton was composed of dermal bony plates covering the body, including odontogenic and skeletogenic components. Later in evolution, the exoskeleton of actinopterygian's trunk was composed of scale structures. Although scales are nowadays a widespread integument cover, some contemporary lineages do not have scales but bony plates covering their trunk, whereas other lineages are devoid of any such structures. To understand the evolution of the integument coverage and particularly the transition between different structures, we investigated the pattern of scale loss events along with actinopterygian evolution and addressed the functional relationship between the scaleless phenotype and the ecology of fishes. Furthermore, we examined whether the emergence of trunk bony plates was dependent over the presence or absence of scales. To this aim, we used two recently published actinopterygian phylogenies, one including >11,600 species, and by using stochastic mapping and Bayesian methods, we inferred scale loss events and trunk bony plate acquisitions. Our results reveal that a scaled integument is the most frequent state in actinopterygians, but multiple independent scale loss events occurred along their phylogeny with essentially no scale re‐acquisition. Based on linear mixed models, we found evidence supporting that after a scale loss event, fishes tend to change their ecology and adopt a benthic lifestyle. Furthermore, we show that trunk bony plates appeared independently multiple times along the phylogeny. By using fitted likelihood models for character evolution, we show that trunk bony plate acquisitions were dependent on a previous scale loss event. Overall, our findings support the hypothesis that integument cover is a key evolutionary trait underlying actinopterygian radiation.

Resolving the Early Divergence Pattern of Teleost Fish Using Genome-Scale Data

Regarding the phylogenetic relationship of the three primary groups of teleost fishes, Osteoglossomorpha (bonytongues and others), Elopomorpha (eels and relatives), Clupeocephala (the remaining teleost fish), early morphological studies hypothesized the first divergence of Osteoglossomorpha, whereas the recent prevailing view is the first divergence of Elopomorpha. Molecular studies supported all the possible relationships of the three primary groups. This study analyzed genome-scale data from four previous studies: 1) 412 genes from 12 species, 2) 772 genes from 15 species, 3) 1,062 genes from 30 species, and 4) 491 UCE loci from 27 species. The effects of the species, loci, and models used on the constructed tree topologies were investigated. In the analyses of the data sets (1)–(3), although the first divergence of Clupeocephala that left the other two groups in a sister relationship was supported by concatenated sequences and gene trees of all the species and genes, the first divergence of Elopomorpha among the three groups was supported using species and/or genes with low divergence of sequence and amino-acid frequencies. This result corresponded to that of the UCE data set (4), whose sequence divergence was low, which supported the first divergence of Elopomorpha with high statistical significance. The increase in accuracy of the phylogenetic construction by using species and genes with low sequence divergence was predicted by a phylogenetic informativeness approach and confirmed by computer simulation. These results supported that Elopomorpha was the first basal group of teleost fish to have diverged, consistent with the prevailing view of recent morphological studies.

The Laboratory Domestication of Zebrafish: From Diverse Populations to Inbred Substrains

We know from human genetic studies that practically all aspects of biology are strongly influenced by the genetic background, as reflected in the advent of “personalized medicine.” Yet, with few exceptions, this is not taken into account when using laboratory populations as animal model systems for research in these fields. Laboratory strains of zebrafish (Danio rerio) are widely used for research in vertebrate developmental biology, behavior, and physiology, for modeling diseases, and for testing pharmaceutic compounds in vivo. However, all of these strains are derived from artificial bottleneck events and therefore are likely to represent only a fraction of the genetic diversity present within the species. Here, we use restriction site-associated DNA sequencing to genetically characterize wild populations of zebrafish from India, Nepal, and Bangladesh, and to compare them to previously published data on four common laboratory strains. We measured nucleotide diversity, heterozygosity, and allele frequency spectra, and find that wild zebrafish are much more diverse than laboratory strains. Further, in wild zebrafish, there is a clear signal of GC-biased gene conversion that is missing in laboratory strains. We also find that zebrafish populations in Nepal and Bangladesh are most distinct from all other strains studied, making them an attractive subject for future studies of zebrafish population genetics and molecular ecology. Finally, isolates of the same strains kept in different laboratories show a pattern of ongoing differentiation into genetically distinct substrains. Together, our findings broaden the basis for future genetic, physiological, pharmaceutic, and evolutionary studies in Danio rerio.

Inferring the number and position of changes in selective regime in a non-equilibrium mutation-selection framework

BACKGROUND

Recovering the historical patterns of selection acting on a protein coding sequence is a major goal of evolutionary biology. Mutation-selection models address this problem by explicitly modelling fixation rates as a function of site-specific amino acid fitness values.However, they are restricted in their utility for investigating directional evolution because they require prior knowledge of the locations of fitness changes in the lineages of a phylogeny.

RESULTS

We apply a modified mutation-selection methodology that relaxes assumptions of equlibrium and time-reversibility. Our implementation allows us to identify branches where adaptive or compensatory shifts in the fitness landscape have taken place, signalled by a change in amino acid fitness profiles. Through simulation and analysis of an empirical data set of [Formula: see text]-lactamase genes, we test our ability to recover the position of adaptive events within the tree and successfully reconstruct initial codon frequencies and fitness profile parameters generated under the non-stationary model.

CONCLUSION

We demonstrate successful detection of selective shifts and identification of the affected branch on partitions of 300 codons or more. We successfully reconstruct fitness parameters and initial codon frequencies in simulated data and demonstrate that failing to account for non-equilibrium evolution can increase the error in fitness profile estimation. We also demonstrate reconstruction of plausible shifts in amino acid fitnesses in the bacterial [Formula: see text]-lactamase family and discuss some caveats for interpretation.

Identification and Characterization of Zebrafish Tlr4 Coreceptor Md-2

The zebrafish (Danio rerio) is a powerful model organism for studies of the innate immune system. One apparent difference between human and zebrafish innate immunity is the cellular machinery for LPS sensing. In amniotes, the protein complex formed by TLR4 and myeloid differentiation factor 2 (Tlr4/Md-2) recognizes the bacterial molecule LPS and triggers an inflammatory response. It is believed that zebrafish have neither Md-2 nor Tlr4; Md-2 has not been identified outside of amniotes, whereas the zebrafish tlr4 genes appear to be paralogs, not orthologs, of amniote TLR4s We revisited these conclusions. We identified a zebrafish gene encoding Md-2, ly96 Using single-cell RNA sequencing, we found that ly96 is transcribed in cells that also transcribe genes diagnostic for innate immune cells, including the zebrafish tlr4-like genes. In larval zebrafish, ly96 is expressed in a small number of macrophage-like cells. In a functional assay, zebrafish Md-2 and Tlr4ba form a complex that activates NF-κB signaling in response to LPS. In larval zebrafish ly96 loss-of-function mutations perturbed LPS-induced cytokine production but gave little protection against LPS toxicity. Finally, by analyzing the genomic context of tlr4 genes in 11 jawed vertebrates, we found that tlr4 arose prior to the divergence of teleosts and tetrapods. Thus, an LPS-sensitive Tlr4/Md-2 complex is likely an ancestral feature shared by mammals and zebrafish, rather than a de novo invention on the tetrapod lineage. We hypothesize that zebrafish retain an ancestral, low-sensitivity Tlr4/Md-2 complex that confers LPS responsiveness to a specific subset of innate immune cells.

Paternal exposure to a common pharmaceutical (Ritalin) has transgenerational effects on the behaviour of Trinidadian guppies

Evidence is emerging that paternal effects, the nongenetic influence of fathers on their offspring, can be transgenerational, spanning several generations. Methylphenidate hydrochloride (MPH; e.g. Ritalin) is a dopaminergic drug that is highly prescribed to adolescent males for the treatment of Attention-deficit/hyperactivity disorder. It has been suggested that MPH could cause transgenerational effects because MPH can affect the male germline in rodents and because paternal effects have been observed in individuals taking similar drugs (e.g. cocaine). Despite these concerns, the transgenerational effects of paternal MPH exposure are unknown. Therefore, we exposed male and female Trinidadian guppies (Poecilia reticulata) to a low, chronic dose of MPH and observed that MPH affected the anxiety/exploratory behaviour of males, but not females. Because of this male-specific effect, we investigated the transgenerational effects of MPH through the paternal line. We observed behavioural effects of paternal MPH exposure on offspring and great-grandoffspring that were not directly administered the drug, making this the first study to demonstrate that paternal MPH exposure can affect descendants. These effects were not due to differential mortality or fecundity between control and MPH lines. These results highlight the transgenerational potential of MPH.

A Long-Term Conserved Satellite DNA That Remains Unexpanded in Several Genomes of Characiformes Fish Is Actively Transcribed

Eukaryotic genomes contain large amounts of repetitive DNA sequences, such as tandemly repeated satellite DNAs (satDNAs). These sequences are highly dynamic and tend to be genus- or species-specific due to their particular evolutionary pathways, although there are few unusual cases of conserved satDNAs over long periods of time. Here, we used multiple approaches to reveal that an satDNA named CharSat01-52 originated in the last common ancestor of Characoidei fish, a superfamily within the Characiformes order, ∼140-78 Ma, whereas its nucleotide composition has remained considerably conserved in several taxa. We show that 14 distantly related species within Characoidei share the presence of this satDNA, which is highly amplified and clustered in subtelomeric regions in a single species (Characidium gomesi), while remained organized as small clusters in all the other species. Defying predictions of the molecular drive of satellite evolution, CharSat01-52 shows similar values of intra- and interspecific divergence. Although we did not provide evidence for a specific functional role of CharSat01-52, its transcriptional activity was demonstrated in different species. In addition, we identified short tandem arrays of CharSat01-52 embedded within single-molecule real-time long reads of Astyanax paranae (536 bp-3.1 kb) and A. mexicanus (501 bp-3.9 kb). Such arrays consisted of head-to-tail repeats and could be found interspersed with other sequences, inverted sequences, or neighbored by other satellites. Our results provide a detailed characterization of an old and conserved satDNA, challenging general predictions of satDNA evolution.

Superoxide Dismutase Multigene Family from a Primitive Chondrostean Sturgeon, Acipenser baerii: Molecular Characterization, Evolution, and Antioxidant Defense during Development and Pathogen Infection

Three distinct superoxide dismutases (SODs)-copper/zinc-SOD (SOD1), manganese-SOD (SOD2), and extracellular copper/zinc-SOD (SOD3)-were identified from a primitive chondrostean fish, Acipenser baerii, enabling the comparison of their transcriptional regulation patterns during development, prelarval ontogeny, and immune stimulation. Each A. baerii SOD isoform (AbSOD) shared conserved structural features with its vertebrate orthologs; however, phylogenetic analyses hypothesized a different evolutionary history for AbSOD3 relative to AbSOD1 and AbSOD2 in the vertebrate lineage. The AbSOD isoforms showed different tissue distribution patterns; AbSOD1 was predominantly expressed in most tissues. The expression of the AbSOD isoforms showed isoform-dependent dynamic modulation according to embryonic development and prelarval ontogenic behaviors. Prelarval microinjections revealed that lipopolysaccharide only induced AbSOD3 expression, while Aeromonas hydrophila induced the expression of AbSOD2 and AbSOD3. In fingerlings, the transcriptional response of each AbSOD isoform to bacterial infection was highly tissue-specific, and the three isoforms exhibited different response patterns within a given tissue type; AbSOD3 was induced the most sensitively, and its induction was the most pronounced in the kidneys and skin. Collectively, these findings suggest isoform-dependent roles for the multigene SOD family in antioxidant defenses against the oxidative stress associated with development and immune responses in these endangered sturgeon fish.

A high-quality de novo genome assembly of one swamp eel (Monopterus albus) strain with PacBio and Hi-C sequencing data

The swamp eel (Monopterus albus) is one economically important fish in China and South-Eastern Asia and a good model species to study sex inversion. There are different genetic lineages and multiple local strains of swamp eel in China, and one local strain of M. albus with deep yellow and big spots has been selected for consecutive selective breeding due to superiority in growth rate and fecundity. A high-quality reference genome of the swamp eel would be a very useful resource for future selective breeding program. In the present study, we applied PacBio single-molecule sequencing technique (SMRT) and the high-throughput chromosome conformation capture (Hi-C) technologies to assemble the M. albus genome. A 799 Mb genome was obtained with the contig N50 length of 2.4 Mb and scaffold N50 length of 67.24 Mb, indicating 110-fold and ∼31.87-fold improvement compared to the earlier released assembly (∼22.24 Kb and 2.11 Mb, respectively). Aided with Hi-C data, a total of 750 contigs were reliably assembled into 12 chromosomes. Using 22,373 protein-coding genes annotated here, the phylogenetic relationships of the swamp eel with other teleosts showed that swamp eel separated from the common ancestor of Zig-zag eel ∼49.9 million years ago, and 769 gene families were found expanded, which are mainly enriched in the immune system, sensory system, and transport and catabolism. This highly accurate, chromosome-level reference genome of M. albus obtained in this work will be used for the development of genome-scale selective breeding.

Evolution of regulatory networks associated with traits under selection in cichlids

BACKGROUND

Seminal studies of vertebrate protein evolution speculated that gene regulatory changes can drive anatomical innovations. However, very little is known about gene regulatory network (GRN) evolution associated with phenotypic effect across ecologically diverse species. Here we use a novel approach for comparative GRN analysis in vertebrate species to study GRN evolution in representative species of the most striking examples of adaptive radiations, the East African cichlids. We previously demonstrated how the explosive phenotypic diversification of East African cichlids can be attributed to diverse molecular mechanisms, including accelerated regulatory sequence evolution and gene expression divergence.

RESULTS

To investigate these mechanisms across species at a genome-wide scale, we develop a novel computational pipeline that predicts regulators for co-extant and ancestral co-expression modules along a phylogeny, and candidate regulatory regions associated with traits under selection in cichlids. As a case study, we apply our approach to a well-studied adaptive trait-the visual system-for which we report striking cases of network rewiring for visual opsin genes, identify discrete regulatory variants, and investigate their association with cichlid visual system evolution. In regulatory regions of visual opsin genes, in vitro assays confirm that transcription factor binding site mutations disrupt regulatory edges across species and segregate according to lake species phylogeny and ecology, suggesting GRN rewiring in radiating cichlids.

CONCLUSIONS

Our approach reveals numerous novel potential candidate regulators and regulatory regions across cichlid genomes, including some novel and some previously reported associations to known adaptive evolutionary traits.

Metagenomic sequencing reveals a lack of virus exchange between native and invasive freshwater fish across the Murray-Darling Basin, Australia

Biological invasions are among the biggest threats to freshwater biodiversity. This is increasingly relevant in the Murray-Darling Basin, Australia, particularly since the introduction of the common carp (Cyprinus carpio). This invasive species now occupies up to ninety per cent of fish biomass, with hugely detrimental impacts on native fauna and flora. To address the ongoing impacts of carp, cyprinid herpesvirus 3 (CyHV-3) has been proposed as a potentially effective biological control agent. Crucially, however, it is unknown whether CyHV-3 and other cyprinid herpesviruses already exist in the Murray-Darling. Further, little is known about those viruses that naturally occur in wild freshwater fauna, and the frequency with which these viruses jump species boundaries. To document the evolution and diversity of freshwater fish viromes and better understand the ecological context to the proposed introduction of CyHV-3, we performed a meta-transcriptomic viral survey of invasive and native fish across the Murray-Darling Basin, covering over 2,200 km of the river system. Across a total of thirty-six RNA libraries representing ten species, we failed to detect CyHV-3 nor any closely related viruses. Rather, meta-transcriptomic analysis identified eighteen vertebrate-associated viruses that could be assigned to the Arenaviridae, Astroviridae, Bornaviridae, Caliciviridae, Coronaviridae, Chuviridae, Flaviviridae, Hantaviridae, Hepeviridae, Paramyxoviridae, Picornaviridae, Poxviridae, Reoviridae and Rhabdoviridae families, and a further twenty-seven that were deemed to be associated with non-vertebrate hosts. Notably, we revealed a marked lack of viruses that are shared among invasive and native fish sampled here, suggesting that there is little virus transmission from common carp to native fish species, despite co-existing for over fifty years. Overall, this study provides the first data on the viruses naturally circulating in a major river system and supports the notion that fish harbour a large diversity of viruses with often deep evolutionary histories.

Diverse phylogenomic datasets uncover a concordant scenario of laurasiatherian interordinal relationships

Resolving the interordinal relationships in the mammalian superorder Laurasiatheria has been among the most intractable problems in higher-level mammalian systematics, with many conflicting hypotheses having been proposed. The present study collected three different sources of genome-scale data with comprehensive taxon sampling of laurasiatherian species, including two protein-coding datasets (4,186 protein-coding genes for an amino acid dataset comprising 2,761,247 amino acid residues and a nucleotide dataset comprising 5,516,340 nucleotides from 1st and 2nd codon positions), an intronic dataset (1,210 introns comprising 1,162,723 nucleotides) and an ultraconserved elements (UCEs) dataset (1,246 UCEs comprising 1,946,472 nucleotides) from 40 species representing all six laurasiatherian orders and 7 non-laurasiatherian outgroups. Remarkably, phylogenetic trees reconstructed with the four datasets using different tree-building methods (RAxML, FastTree, ASTRAL and MP-EST) all supported the relationship (Eulipotyphla, (Chiroptera, ((Carnivora, Pholidota), (Cetartiodactyla, Perissodactyla)))). We find a resolution of interordinal relationships of Laurasiatheria among all types of markers used in the present study, and the likelihood ratio tests for tree comparisons confirmed that the present tree topology is the optimal hypothesis compared to other examined hypotheses. Jackknifing subsampling analyses demonstrate that the results of laurasiatherian tree reconstruction varied with the number of loci and ordinal representatives used, which are likely the two main contributors to phylogenetic disagreements of Laurasiatheria seen in previous studies. Our study provides significant insight into laurasiatherian evolution, and moreover, an important methodological strategy and reference for resolving phylogenies of adaptive radiation, which have been a long-standing challenge in the field of phylogenetics.

Evolutionary determinism and convergence associated with water-column transitions in marine fishes

Repeatable, convergent outcomes are prima facie evidence for determinism in evolutionary processes. Among fishes, well-known examples include microevolutionary habitat transitions into the water column, where freshwater populations (e.g., sticklebacks, cichlids, and whitefishes) recurrently diverge toward slender-bodied pelagic forms and deep-bodied benthic forms. However, the consequences of such processes at deeper macroevolutionary scales in the marine environment are less clear. We applied a phylogenomics-based integrative, comparative approach to test hypotheses about the scope and strength of convergence in a marine fish clade with a worldwide distribution (snappers and fusiliers, family Lutjanidae) featuring multiple water-column transitions over the past 45 million years. We collected genome-wide exon data for 110 (∼80%) species in the group and aggregated data layers for body shape, habitat occupancy, geographic distribution, and paleontological and geological information. We also implemented approaches using genomic subsets to account for phylogenetic uncertainty in comparative analyses. Our results show independent incursions into the water column by ancestral benthic lineages in all major oceanic basins. These evolutionary transitions are persistently associated with convergent phenotypes, where deep-bodied benthic forms with truncate caudal fins repeatedly evolve into slender midwater species with furcate caudal fins. Lineage diversification and transition dynamics vary asymmetrically between habitats, with benthic lineages diversifying faster and colonizing midwater habitats more often than the reverse. Convergent ecological and functional phenotypes along the benthic-pelagic axis are pervasive among different lineages and across vastly different evolutionary scales, achieving predictable high-fitness solutions for similar environmental challenges, ultimately demonstrating strong determinism in fish body-shape evolution.

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.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Phylogenetic reconstruction of early diverging tapeworms (Cestoda: Caryophyllidea) reveals ancient radiations in vertebrate hosts and biogeographic regions

Tapeworms of the order Caryophyllidea are the earliest diverging 'true' tapeworms (Eucestoda) and parasitise cypriniform and siluriform fishes almost exclusively. They are typified by a monozoic (non-proglottised) body plan, which is a characteristic shared with early diverging 'cestodarians' Gyrocotylidea and Amphilinidea. Here we present the most comprehensive multi-gene molecular phylogeny of this group, to date. Specimens of 63 species from 32 genera (~50% and ~75% of known species and genus diversity, respectively) were gathered during an intense and targeted 15-year collecting effort. Phylogenetic reconstructions provide high nodal support for three major lineages, which only partly correspond to currently recognised families. The three well-supported clades were as follows: Clade A was in an unsupported position at the base of the tree and was almost exclusively comprised of parasites of catfishes (Siluriformes) from the Afrotropical and Indomalayan regions, including the type genus of the Lytocestidae (Lytocestus). Clade B formed the sister group to the remaining taxa (Clade C) and was composed of species that parasitise cyprinids and loaches (Cypriniformes: Cyprinoidei and Cobitoidei) from the Palaearctic Region. This clade included the type genus of the Caryophyllaeidae (Caryophyllaeus). Clade C comprised Nearctic species from suckers and minnows (Cypriniformes: Catostomidae and Cyprinoidei), which were previously accommodated in two families, i.e. Capingentidae and Caryophyllaeidae. This clade included the type genus of the Capingentidae (Capingens). In addition to Clades A-C, Balanotaenia bancrofti from the monotypic Balanotaeniidae, which parasitises plotosid catfishes in Australia, and Lytocestoides tanganyikae, which parasitises African cichlids, formed a poorly supported clade at the base of the tree. Whereas morphological characteristics traditionally used to differentiate caryophyllidean families do not characterise molecular lineages, host association and biogeographical distribution play a key role in the circumscription of the three well-supported clades revealed by molecular data. Thus, the taxonomic rearrangement proposed herein was guided by the molecular clades. The names of all four extant families were preserved and family affinity was determined by topological clustering with the type genera of the families. The family diagnoses of the Lytocestidae, Caryophyllaeidae and Capingentidae are amended. Biogeographic patterns are indicative of separate Gondwanan and Laurasian radiations having taken place. Regarding the Gondwanan radiation in the Siluriformes, the topology in Clade A indicates an Asian origin with a subsequent African colonisation. Concerning Laurasia, separate radiations appear to have taken place in the Cypriniformes in the temperate zones of North America and Eurasia. Complete absence of caryophyllideans in the Neotropical Region, where numerous catfishes occur, may be due to the Gondwanan radiation having taken place after the continental separation of Africa and South America.

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

The genomic timeline of cichlid fish diversification across continents

Cichlid fishes are celebrated for their vast taxonomic, phenotypic, and ecological diversity; however, a central aspect of their evolution - the timeline of their diversification - remains contentious. Here, we generate draft genome assemblies of 14 species representing the global cichlid diversity and integrate these into a new phylogenomic hypothesis of cichlid and teleost evolution that we time-calibrate with 58 re-evaluated fossil constraints and a new Bayesian model accounting for fossil-assignment uncertainty. Our results support cichlid diversification long after the breakup of the supercontinent Gondwana and lay the foundation for precise temporal reconstructions of the exceptional continental cichlid adaptive radiations.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.

Postcranial skeletal development of Mugil cephalus (Teleostei: Mugiliformes): morphological and life-history implications for Mugiliformes

Within the fish taxon Mugiliformes, the larval development of Mugil cephalus has been studied most intensively, because it has the widest range of distribution among all mugilids and is of interest to aquaculture all over the world. Although numerous studies have dealt with larval rearing, growth and development, the osteological development of M. cephalus and mugiliforms in general has largely been neglected. Herein, we describe the skeletal development of mullets for the first time. Cleared and double-stained specimens of aquaculture-reared M. cephalus and wild-caught mugilid larvae were examined to describe the early development of the pectoral and pelvic girdle, the vertebral column and the caudal and median fins. The description of four embryonic and six larval developmental steps within the embryonic and larval period enables us to compare larval sizes of reared and wild-caught larvae. Ontogenetic fusions of ural centra 1 and 2 into a compound centrum, in addition to the fusion of two pterygiophores in the anal fin, have implications for the perception of the adult morphology. Moreover, comparison of mugilid development with that of other ovalentarian taxa shows that recent phylogenetic hypotheses need further morphological investigation.

Tumor Necrosis Factor Superfamily: Ancestral Functions and Remodeling in Early Vertebrate Evolution

The evolution of the tumor necrosis factor superfamily (TNFSF) in early vertebrates is inferred by comparing the TNFSF genes found in humans and nine fishes: three agnathans, two chondrichthyans, three actinopterygians, and the sarcopterygian Latimeria chalumnae. By combining phylogenetic and synteny analyses, the TNFSF sequences detected are classified into five clusters of genes and 24 orthology groups. A model for their evolution since the origin of vertebrates is proposed. Fifteen TNFSF genes emerged from just three progenitors due to the whole-genome duplications (WGDs) that occurred before the agnathan/gnathostome split. Later, gnathostomes not only kept most of the genes emerged in the WGDs but soon added several tandem duplicates. More recently, complex, lineage-specific patterns of duplications and losses occurred in different gnathostome lineages. In agnathan species only seven to eight TNFSF genes are detected, because this lineage soon lost six of the genes emerged in the ancestral WGDs and additional losses in both hagfishes and lampreys later occurred. The orthologs of many of these lost genes are, in mammals, ligands of death-domain-containing TNFSF receptors, indicating that the extrinsic apoptotic pathway became simplified in the agnathan lineage. From the patterns of emergence of these genes, it is deduced that both the regulation of apoptosis and the control of the NF-κB pathway that depends in modern mammals on TNFSF members emerged before the ancestral vertebrate WGDs.

Anatomy and evolution of the pectoral filaments of threadfins (Polynemidae)

The most remarkable anatomical specialization of threadfins (Percomorphacea: Polynemidae) is the division of their pectoral fin into an upper, unmodified fin and a lower portion with rays highly modified into specialized filaments. Such filaments are usually elongate, free from interradial membrane, and move independently from the unmodified fin to explore the environment. The evolution of the pectoral filaments involved several morphological modifications herein detailed for the first time. The posterior articular facet of the coracoid greatly expands anteroventrally during development. Similar expansions occur in pectoral radials 3 and 4, with the former usually acquiring indentations with the surrounding bones and losing association with both rays and filaments. Whereas most percomorphs typically have four or five muscles serving the pectoral fin, adult polynemids have up to 11 independent divisions in the intrinsic pectoral musculature. The main adductor and abductor muscles masses of the pectoral system are completely divided into two muscle segments, each independently serving the pectoral-fin rays (dorsally) and the pectoral filaments (ventrally). Based on the innervation pattern and the discovery of terminal buds in the external surface of the filaments, we demonstrate for the first time that the pectoral filaments of threadfins have both tactile and gustatory functions.

Chromosome-level genome assembly for the largemouth bass Micropterus salmoides provides insights into adaptation to fresh and brackish water

Largemouth bass (LMB; Micropterus salmoides) has been an economically important fish in North America, Europe, and China. This study obtained a chromosome-level genome assembly of LMB using PacBio and Hi-C sequencing. The final assembled genome is 964 Mb, with contig N50 and scaffold N50 values of 1.23 Mb and 36.48 Mb, respectively. Combining with RNA sequencing data, we annotated a total of 23,701 genes. Chromosomal assembly and syntenic analysis proved that, unlike most Perciformes with the popular haploid chromosome number of 24, LMB has only 23 chromosomes (Chr), among which the Chr1 seems to be resulted from a chromosomal fusion event. LMB is phylogenetically closely related to European seabass and spotted seabass, diverging 64.1 million years ago (mya) from the two seabass species. Eight gene families comprising 294 genes associated with ionic regulation were identified through positive selection, transcriptome and genome comparisons. These genes involved in iron facilitated diffusion (such as claudin, aquaporins, sodium channel protein and so on) and others related to ion active transport (such as sodium/potassium-transporting ATPase and sodium/calcium exchanger). The claudin gene family, which is critical for regulating cell tight junctions and osmotic homeostasis, showed a significant expansion in LMB with 27 family members and 68 copies for salinity adaptation. In summary, we reported the first high-quality LMB genome, and provided insights into the molecular mechanisms of LMB adaptation to fresh and brackish water. The chromosome-level LMB genome will also be a valuable genomic resource for in-depth biological and evolutionary studies, germplasm conservation and genetic breeding of LMB.

New Environment, New Invaders-Repeated Horizontal Transfer of LINEs to Sea Snakes

Although numerous studies have found horizontal transposon transfer (HTT) to be widespread across metazoans, few have focused on HTT in marine ecosystems. To investigate potential recent HTTs into marine species, we searched for novel repetitive elements in sea snakes, a group of elapids which transitioned to a marine habitat at most 18 Ma. Our analysis uncovered repeated HTTs into sea snakes following their marine transition. The seven subfamilies of horizontally transferred LINE retrotransposons we identified in the olive sea snake (Aipysurus laevis) are transcribed, and hence are likely still active and expanding across the genome. A search of 600 metazoan genomes found all seven were absent from other amniotes, including terrestrial elapids, with the most similar LINEs present in fish and marine invertebrates. The one exception was a similar LINE found in sea kraits, a lineage of amphibious elapids which independently transitioned to a marine environment 25 Ma. Our finding of repeated horizontal transfer events into marine snakes greatly expands past findings that the marine environment promotes the transfer of transposons. Transposons are drivers of evolution as sources of genomic sequence and hence genomic novelty. We identified 13 candidate genes for HTT-induced adaptive change based on internal or neighboring HTT LINE insertions. One of these, ADCY4, is of particular interest as a part of the KEGG adaptation pathway “Circadian Entrainment.” This provides evidence of the ecological interactions between species influencing evolution of metazoans not only through specific selection pressures, but also by contributing novel genomic material.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Contingency between Historical Substitutions in the Acetylcholine Receptor Pore

Human adult muscle-type acetylcholine receptors incorporating a reconstructed ancestral β-subunit exhibit reduced single-channel conductance when compared to wild-type. The ancestral and wild-type β-subunits differ by 132 amino acids, including substitution of residues that line the lumen of the channel pore, near its narrowest constriction. Here we show that a single historical substitution in this region of the ancestral β-subunit accounts for the difference in conductance. Furthermore, the contribution of the substituted residue to conductance is dependent upon its ancestral or wild-type background, and it can be modulated by a neighboring residue that has also evolved throughout the β-subunit history. Using an expanded molecular phylogeny, we track the order in which these two mutations occurred and then show that the order in which they are installed upon the ancestral, but not the human, background determines their individual contribution to conductance. Our results show how the contribution of amino acids to acetylcholine receptor conductance is contingent upon their evolutionary history and that the order in which substitutions occurred was important for shaping conductance in the modern-day receptor.

Intragenic Conflict in Phylogenomic Data Sets

Most phylogenetic analyses assume that a single evolutionary history underlies one gene. However, both biological processes and errors can cause intragenic conflict. The extent to which this conflict is present in empirical data sets is not well documented, but if common, could have far-reaching implications for phylogenetic analyses. We examined several large phylogenomic data sets from diverse taxa using a fast and simple method to identify well-supported intragenic conflict. We found conflict to be highly variable between data sets, from 1% to >92% of genes investigated. We analyzed four exemplar genes in detail and analyzed simulated data under several scenarios. Our results suggest that alignment error may be one major source of conflict, but other conflicts remain unexplained and may represent biological signal or other errors. Whether as part of data analysis pipelines or to explore biologically processes, analyses of within-gene phylogenetic signal should become common.

Development of the ventricular myocardial trabeculae in Scyliorhinus canicula (Chondrichthyes): evolutionary implications

The development of the ventricular myocardial trabeculae occurs in three steps: emergence, trabeculation and remodeling. The whole process has been described in vertebrates with two different myocardial structural types, spongy (zebrafish) and compact (chicken and mouse). In this context, two alternative mechanisms of myocardial trabeculae emergence have been identified: (1) in chicken and mouse, the endocardial cells invade the two-layered myocardium; (2) in zebrafish, cardiomyocytes from the monolayered myocardium invaginate towards the endocardium. Currently, the process has not been studied in detail in vertebrates having a mixed type of ventricular myocardium, with an inner trabecular and an outer compact layer, which is presumptively the most primitive morphology in gnathostomes. We studied the formation of the mixed ventricular myocardium in the lesser spotted dogfish (Scyliorhinus canicula, Elasmobranchii), using light, scanning and transmission electron microscopy. Our results show that early formation of the mixed ventricular myocardium, specifically the emergence and the trabeculation steps, is driven by an endocardial invasion of the myocardium. The mechanism of trabeculation of the mixed ventricular myocardium in chondrichthyans is the one that best reproduces how this developmental process has been established from the beginning of the gnathostome radiation. The process has been apparently preserved throughout the entire group of sarcopterygians, including birds and mammals. In contrast, teleosts, at least those possessing a mostly spongy ventricular myocardium, seem to have introduced notable changes in their myocardial trabeculae development.

An intestinal cell type in zebrafish is the nexus for the SARS-CoV-2 receptor and the Renin-Angiotensin-Aldosterone System that contributes to COVID-19 comorbidities

People with underlying conditions, including hypertension, obesity, and diabetes, are especially susceptible to negative outcomes after infection with the coronavirus SARS-CoV-2. These COVID-19 comorbidities are exacerbated by the Renin-Angiotensin-Aldosterone System (RAAS), which normally protects from rapidly dropping blood pressure or dehydration via the peptide Angiotensin II (Ang II) produced by the enzyme Ace. The Ace paralog Ace2 degrades Ang II, thus counteracting its chronic effects. Ace2 is also the SARS-CoV-2 receptor. Ace , the coronavirus, and COVID-19 comorbidities all regulate Ace2 , but we don't yet understand how. To exploit zebrafish ( Danio rerio ) as a disease model to understand mechanisms regulating the RAAS and its relationship to COVID-19 comorbidities, we must first identify zebrafish orthologs and co-orthologs of human RAAS genes, and second, understand where and when these genes are expressed in specific cells in zebrafish development. To achieve these goals, we conducted genomic analyses and investigated single cell transcriptomes. Results showed that most human RAAS genes have an ortholog in zebrafish and some have two or more co-orthologs. Results further identified a specific intestinal cell type in zebrafish larvae as the site of expression for key RAAS components, including Ace, Ace2, the coronavirus co-receptor Slc6a19, and the Angiotensin-related peptide cleaving enzymes Anpep and Enpep. Results also identified specific vascular cell subtypes as expressing Ang II receptors, apelin , and apelin receptor genes. These results identify specific genes and cell types to exploit zebrafish as a disease model for understanding the mechanisms leading to COVID-19 comorbidities.

SUMMARY STATEMENT

Genomic analyses identify zebrafish orthologs of the Renin-Angiotensin-Aldosterone System that contribute to COVID-19 comorbidities and single-cell transcriptomics show that they act in a specialized intestinal cell type.

Migratory patterns and evolutionary plasticity of cranial neural crest cells in ray-finned fishes

The cranial neural crest (CNC) arises within the developing central nervous system, but then migrates away from the neural tube in three consecutive streams termed mandibular, hyoid and branchial, respectively, according to the order along the anteroposterior axis. While the process of neural crest emigration generally follows a conserved anterior to posterior sequence across vertebrates, we find that ray-finned fishes (bichir, sterlet, gar, and pike) exhibit several heterochronies in the timing and order of CNC emergence that influences their subsequent migratory patterns. First, emigration of the cranial neural crest in these fishes occurs prematurely compared to other vertebrates, already initiating during early neurulation and well before neural tube closure. Second, delamination of the hyoid stream occurs prior to the more anterior mandibular stream; this is associated with early morphogenesis of key hyoid structures like external gills (bichir), a large opercular flap (gar) or first forming cartilage (pike). In sterlet, the hyoid and branchial CNC cells form a single hyobranchial sheet, which later segregates in concert with second pharyngeal pouch morphogenesis. Taken together, the results show that despite generally conserved migratory patterns, heterochronic alterations in the timing of emigration and pattern of migration of CNC cells accompanies morphological diversity of ray-finned fishes.

Genomic Analysis of the Only Blind Cichlid Reveals Extensive Inactivation in Eye and Pigment Formation Genes

Trait loss represents an intriguing evolutionary problem, particularly when it occurs across independent lineages. Fishes in light-poor environments often evolve “troglomorphic” traits, including reduction or loss of both pigment and eyes. Here, we investigate the genomic basis of trait loss in a blind and depigmented African cichlid, Lamprologus lethops, and explore evolutionary forces (selection and drift) that may have contributed to these losses. This species, the only known blind cichlid, is endemic to the lower Congo River. Available evidence suggests that it inhabits deep, low-light habitats. Using genome sequencing, we show that genes related to eye formation and pigmentation, as well as other traits associated with troglomorphism, accumulated inactivating mutations rapidly after speciation. A number of the genes affected in L. lethops are also implicated in troglomorphic phenotypes in Mexican cavefish (Astyanax mexicanus) and other species. Analysis of heterozygosity patterns across the genome indicates that L. lethops underwent a significant population bottleneck roughly 1 Ma, after which effective population sizes remained low. Branch-length tests on a subset of genes with inactivating mutations show little evidence of directional selection; however, low overall heterozygosity may reduce statistical power to detect such signals. Overall, genome-wide patterns suggest that accelerated genetic drift from a severe bottleneck, perhaps aided by directional selection for the loss of physiologically expensive traits, caused inactivating mutations to fix rapidly in this species.

The immunogenetics of sexual parasitism

Sexual parasitism has evolved as a distinctive mode of reproduction among deep-sea anglerfishes. The permanent attachment of males to host females observed in these species represents a form of anatomical joining, which is otherwise unknown in nature. Pronounced modifications to immune facilities are associated with this reproductive trait. The genomes of species with temporarily attaching males lack functional aicda genes that underpin affinity maturation of antibodies. Permanent attachment is associated with additional alterations, culminating in the loss of functional rag genes in some species, abolishing somatic diversification of antigen receptor genes, the hallmark of canonical adaptive immunity. In anglerfishes, coevolution of innate and adaptive immunity has been disentangled, implying that an alternative form of immunity supported the emergence of this evolutionarily successful group of vertebrates.

Unravelling the Complex Duplication History of Deuterostome Glycerol Transporters

Transmembrane glycerol transport is an ancient biophysical property that evolved in selected subfamilies of water channel (aquaporin) proteins. Here, we conducted broad level genome (>550) and transcriptome (>300) analyses to unravel the duplication history of the glycerol-transporting channels (glps) in Deuterostomia. We found that tandem duplication (TD) was the major mechanism of gene expansion in echinoderms and hemichordates, which, together with whole genome duplications (WGD) in the chordate lineage, continued to shape the genomic repertoires in craniates. Molecular phylogenies indicated that aqp3-like and aqp13-like channels were the probable stem subfamilies in craniates, with WGD generating aqp9 and aqp10 in gnathostomes but aqp7 arising through TD in Osteichthyes. We uncovered separate examples of gene translocations, gene conversion, and concerted evolution in humans, teleosts, and starfishes, with DNA transposons the likely drivers of gene rearrangements in paleotetraploid salmonids. Currently, gene copy numbers and BLAST are poor predictors of orthologous relationships due to asymmetric glp gene evolution in the different lineages. Such asymmetries can impact estimations of divergence times by millions of years. Experimental investigations of the salmonid channels demonstrated that approximately half of the 20 ancestral paralogs are functional, with neofunctionalization occurring at the transcriptional level rather than the protein transport properties. The combined findings resolve the origins and diversification of glps over >800 million years old and thus form the novel basis for proposing a pandeuterostome glp gene nomenclature.

Liver transcriptome resources of four commercially exploited teleost species

The generation of omic resources is central to develop adequate management strategies for species with economic value. Here, we provide high-coverage RNA-seq datasets of liver tissue (containing between 80,2 and 88,4 million of paired-end reads) from four wildtype teleost species with high commercial value: Trachurus trachurus (TTR; Atlantic horse mackerel), Scomber scombrus (SSC; Atlantic mackerel), Trisopterus luscus (TLU; pout), and Micromesistius poutassou (MPO; blue whiting). A comprehensive assembly pipeline, using de novo single and multi-kmer assembly approaches, produced 64 single high-quality liver transcriptomes - 16 per species. The final assemblies, with N50 values ranging from 2,543-3,700 bp and BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness values between 81.8-86.5% of the Actinopterygii gene set, were subjected to open reading frame (ORF) prediction and functional annotation. Our study provides the first transcriptomic resources for these species and offers valuable tools to evaluate both neutral and selected genetic variation among populations, and to identify candidate genes for environmental adaptation assisting in the investigation of the effects of global changes in fisheries.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

A fish cytokine related to human IL-3, IL-5, and GM-CSF, induces development of eosinophil/basophil/mast-cell type (EBM) granulocytes

Interleukin-3 (IL-3), IL-5, and granulocyte-macrophage colony-stimulating factor (GM-CSF) are related cytokines that signal through receptors possessing the β common (βc) chain. As a family, these cytokines combine rather non-specific hematopoietic growth factor properties with a special importance for eosinophils, basophils, and mast cells. In fish the cytokines of this family are called IL-5fam, and the present study, using carp, constitutes their first functional analysis. Carp il-5fam expression was enhanced by stimulation with phytohemagglutinin and killed bacteria. Reminiscent of mammalian IL-3/IL-5/GM-CSF family members, recombinant carp IL-5fam (rcIL-5fam) induced activation of transcription factor STAT5 and efficiently promoted proliferation and colony-formation of eosinophil/basophil/mast-cell type (EBM) granulocytes. Upon addition of recombinant carp βc the growth effect of rcIL-5fam was reduced, suggesting βc participation in the signaling route. In summary, despite differences in individual cytokines and cell populations, fish and mammalian IL-3/IL-5/GM-CSF family members share growth factor functions for non-neutrophil granulocytes.

Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis

Medaka (Oryzias latipes) has become an important vertebrate model widely used in genetics, developmental biology, environmental sciences, and many other fields. A high-quality genome sequence and a variety of genetic tools are available for this model organism. However, existing genome annotation is still rudimentary, as it was mainly based on computational prediction and short-read RNA-seq data. Here we report a dynamic transcriptome landscape of medaka embryogenesis profiled by long-read RNA-seq, short-read RNA-seq, and ATAC-seq. By integrating these data sets, we constructed a much-improved gene model set including about 17,000 novel isoforms and identified 1600 transcription factors, 1100 long noncoding RNAs, and 150,000 potential cis-regulatory elements as well. Time-series data sets provided another dimension of information. With the expression dynamics of genes and accessibility dynamics of cis-regulatory elements, we investigated isoform switching, as well as regulatory logic between accessible elements and genes, during embryogenesis. We built a user-friendly medaka omics data portal to present these data sets. This resource provides the first comprehensive omics data sets of medaka embryogenesis. Ultimately, we term these three assays as the minimum ENCODE toolbox and propose the use of it as the initial and essential profiling genomic assays for model organisms that have limited data available. This work will be of great value for the research community using medaka as the model organism and many others as well.

Visual pigment evolution in Characiformes: The dynamic interplay of teleost whole-genome duplication, surviving opsins and spectral tuning

Vision represents an excellent model for studying adaptation, given the genotype-to-phenotype map that has been characterized in a number of taxa. Fish possess a diverse range of visual sensitivities and adaptations to underwater light, making them an excellent group to study visual system evolution. In particular, some speciose but understudied lineages can provide a unique opportunity to better understand aspects of visual system evolution such as opsin gene duplication and neofunctionalization. In this study, we showcase the visual system evolution of neotropical Characiformes and the spectral tuning mechanisms they exhibit to modulate their visual sensitivities. Such mechanisms include gene duplications and losses, gene conversion, opsin amino acid sequence and expression variation, and A₁ /A₂ -chromophore shifts. The Characiforms we studied utilize three cone opsin classes (SWS2, RH2, LWS) and a rod opsin (RH1). However, the characiform's entire opsin gene repertoire is a product of dynamic evolution by opsin gene loss (SWS1, RH2) and duplication (LWS, RH1). The LWS- and RH1-duplicates originated from a teleost specific whole-genome duplication as well as characiform-specific duplication events. Both LWS-opsins exhibit gene conversion and, through substitutions in key tuning sites, one of the LWS-paralogues has acquired spectral sensitivity to green light. These sequence changes suggest reversion and parallel evolution of key tuning sites. Furthermore, characiforms' colour vision is based on the expression of both LWS-paralogues and SWS2. Finally, we found interspecific and intraspecific variation in A₁ /A₂ -chromophores proportions, correlating with the light environment. These multiple mechanisms may be a result of the diverse visual environments where Characiformes have evolved.

The osteology of the golden grey mullet Liza aurata (Teleostei: Mugiliformes: Mugilidae) including interactive three-dimensional reconstructions

Grey mullets are remarkably characterized by their overall uniform external morphology. Identifying species as well as positioning the Mugiliformes in a phylogenetic context is rather difficult. Most recently they were placed in the newly erected Ovalentaria, but more detailed relationships to potential sister taxa were not resolved. Studying the internal morphology, especially the osteology, might provide new insights into the evolution of the Mugiliformes as well as help clarify their phylogenetic position within the Ovalentaria. A detailed osteology of the golden grey mullet Liza aurata is presented. The use of cleared and stained specimens allowed for a complete examination of bony and cartilaginous structures, and a 3D reconstruction from a μCT data set provided additional information on the positional relationships of the bones. Following this, the data obtained were compared with different mugilid species, particularly with the flathead grey mullet Mugil cephalus. Several differences between these species could be identified, such as the position of the basisphenoid, the shape of the hyomandibular and the composition of the branchial arches. These characters might help in understanding the evolutionary changes happening within the mugiliforms and will provide the basis to study this taxon in detail, finally allowing the reconstruction of the body plan of grey mullets.

Extreme genomic volatility characterizes the evolution of the immunoglobulin heavy chain locus in cyprinodontiform fishes

The evolution of the adaptive immune system has provided vertebrates with a uniquely sophisticated immune toolkit, enabling them to mount precise immune responses against a staggeringly diverse range of antigens. Like other vertebrates, teleost fishes possess a complex and functional adaptive immune system; however, our knowledge of the complex antigen-receptor genes underlying its functionality has been restricted to a small number of experimental and agricultural species, preventing systematic investigation into how these crucial gene loci evolve. Here, we analyse the genomic structure of the immunoglobulin heavy chain (IGH) gene loci in the cyprinodontiforms, a diverse and important group of teleosts present in many different habitats across the world. We reconstruct the complete IGH loci of the turquoise killifish (Nothobranchius furzeri) and the southern platyfish (Xiphophorus maculatus) and analyse their in vivo gene expression, revealing the presence of species-specific splice isoforms of transmembrane IGHM. We further characterize the IGH constant regions of 10 additional cyprinodontiform species, including guppy, Amazon molly, mummichog and mangrove killifish. Phylogenetic analysis of these constant regions suggests multiple independent rounds of duplication and deletion of the teleost-specific antibody class IGHZ in the cyprinodontiform lineage, demonstrating the extreme volatility of IGH evolution. Focusing on the cyprinodontiforms as a model taxon for comparative evolutionary immunology, this work provides novel genomic resources for studying adaptive immunity and sheds light on the evolutionary history of the adaptive immune system.

The genome of Mekong tiger perch (Datnioides undecimradiatus) provides insights into the phylogenetic position of Lobotiformes and biological conservation

Mekong tiger perch (Datnioides undecimradiatus) is an ornamental and vulnerable freshwater fish native to the Mekong basin in Indochina, belonging to the order Lobotiformes. Here, we generated 121X stLFR co-barcode clean reads and 18X Oxford Nanopore MinION reads and obtained a 595 Mb Mekong tiger perch genome, which is the first whole genome sequence in the order Lobotiformes. Based on this genome, the phylogenetic tree analysis suggested that Lobotiformes is more closely related to Sciaenidae than to Tetraodontiformes, resolving a long-time dispute. We depicted the genes involved in pigment development in Mekong tiger perch and results confirmed that the four rate-limiting genes of pigment synthesis had been retained after fish-specific genome duplication. We also estimated the demographic history of Mekong tiger perch, which showed that the effective population size suffered a continuous reduction possibly related to the contraction of immune-related genes. Our study provided a reference genome resource for the Lobotiformes, as well as insights into the phylogenetic position of Lobotiformes and biological conservation.

Basal teleosts provide new insights into the evolutionary history of teleost-duplicated aromatase

Duplicated cyp19a1 genes (cyp19a1a encoding aromatase a and cyp19a1b encoding aromatase b) have been identified in an increasing number of teleost species. Cyp19a1a is mainly expressed in the gonads, while cyp19a1b is mainly expressed in the brain, specifically in radial glial cells, as largely investigated by Kah and collaborators. The third round of whole-genome duplication that specifically occurred in the teleost lineage (TWGD or 3R) is likely at the origin of the duplicated cyp19a1 paralogs. In contrast to the situation in other teleosts, our previous studies identified a single cyp19a1 in eels (Anguilla), which are representative species of a basal group of teleosts, Elopomorpha. In the present study, using genome data mining and phylogenetic and synteny analyses, we confirmed that the whole aromatase genomic region was duplicated in eels, with most aromatase-neighboring genes being conserved in duplicate in eels, as in other teleosts. These findings suggest that specific gene loss of one of the 3R-duplicated cyp19a1 paralogs occurred in Elopomorpha after TWGD. Similarly, a single cyp19a1 gene was found in the arowana, which is a representative species of another basal group of teleosts, Osteoglossomorpha. In eels, the single cyp19a1 is expressed in both the brain and the gonads, as observed for the single CYP19A1 gene present in other vertebrates. The results of phylogenetic, synteny, closest neighboring gene, and promoter structure analyses showed that the single cyp19a1 of the basal teleosts shared conserved properties with both teleost cyp19a1a and cyp19a1b paralogs, which did not allow us to conclude which of the 3R-duplicated paralogs (cyp19a1a or cyp19a1b) was lost in Elopomorpha. Elopomorpha and Osteoglossomorpha cyp19a1 genes exhibited preserved ancestral functions, including expression in both the gonad and brain. We propose that the subfunctionalization of the 3R-duplicated cyp19a1 paralogs expressed specifically in the gonad or brain occurred in Clupeocephala, after the split of Clupeocephala from Elopomorpha and Osteoglossomorpha, which represented a driving force for the conservation of both 3R-duplicated paralogs in all extant Clupeocephala. In contrast, the functional redundancy of the undifferentiated 3R-duplicated cyp19a1 paralogs in elopomorphs and osteoglossomorphs would have favored the loss of one 3R paralog in basal teleosts.

Linking Virus Discovery to Immune Responses Visualized during Zebrafish Infections

The discovery of new viruses currently outpaces our capacity for experimental examination of infection biology. To better couple virus discovery with immunology, we genetically modified zebrafish to visually report on virus infections. After generating a strain that expresses green fluorescent protein (GFP) under an interferon-stimulated gene promoter, we repeatedly observed transgenic larvae spontaneously expressing GFP days after hatching. RNA sequencing comparisons of co-housed GFP-positive and GFP-negative zebrafish revealed a naturally occurring picornavirus that induced a canonical interferon-mediated response and hundreds of antiviral defense genes not observed following immunostimulatory treatments or experimental infections with other viruses. Among the many genes induced by picornavirus infection was a large set encoding guanosine triphosphatase (GTPase) of immunity-associated proteins (GIMAPs). The GIMAP gene family is massively expanded in fish genomes and may also play a crucial role in antiviral responses in mammals, including humans. We subsequently detected zebrafish picornavirus in publicly available sequencing data from seemingly asymptomatic zebrafish in many research institutes and found that it altered gene expression in a previous study of zebrafish development. Experiments revealed a horizontal mode of virus transmission, highlighting a system for studying the spread of picornavirus infections within and between individuals. Our study describes a naturally occurring picornavirus that elicits strong antiviral responses in zebrafish and provides new strategies for simultaneously discovering viruses and their impact on vertebrate hosts.