The interactome of KRAB zinc finger proteins reveals the evolutionary history of their functional diversification.

Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) are encoded in the hundreds by the genomes of higher vertebrates, and many act with the heterochromatin-inducing KAP1 as repressors of transposable elements (TEs) during early embryogenesis. Yet, their widespread expression in adult tissues and enrichment at other genetic loci indicate additional roles. Here, we characterized the protein interactome of 101 of the ~350 human KZFPs. Consistent with their targeting of TEs, most KZFPs conserved up to placental mammals essentially recruit KAP1 and associated effectors. In contrast, a subset of more ancient KZFPs rather interacts with factors related to functions such as genome architecture or RNA processing. Nevertheless, KZFPs from coelacanth, our most distant KZFP-encoding relative, bind the cognate KAP1. These results support a hypothetical model whereby KZFPs first emerged as TE-controlling repressors, were continuously renewed by turnover of their hosts' TE loads, and occasionally produced derivatives that escaped this evolutionary flushing by development and exaptation of novel functions.

Diverse Eukaryotic CGG-Binding Proteins Produced by Independent Domestications of hAT Transposons.

The human transcription factor (TF) CGGBP1 (CGG-binding protein) is conserved only in amniotes and is believed to derive from the zf-BED and Hermes transposase DNA-binding domains (DBDs) of a hAT DNA transposon. Here, we show that sequence-specific DNA-binding proteins with this bipartite domain structure have resulted from dozens of independent hAT domestications in different eukaryotic lineages. CGGBPs display a wide range of sequence specificity, usually including preferences for CGG or CGC trinucleotides, whereas some bind AT-rich motifs. The CGGBPs are almost entirely nonsyntenic, and their protein sequences, DNA-binding motifs, and patterns of presence or absence in genomes are uncharacteristic of ancestry via speciation. At least eight CGGBPs in the coelacanth Latimeria chalumnae bind distinct motifs, and the expression of the corresponding genes varies considerably across tissues, suggesting tissue-restricted function.

The KRAB Domain of ZNF10 Guides the Identification of Specific Amino Acids That Transform the Ancestral KRAB-A-Related Domain Present in Human PRDM9 into a Canonical Modern KRAB-A Domain.

Krüppel-associated box (KRAB) zinc finger proteins are a large class of tetrapod transcription factors that usually exert transcriptional repression through recruitment of TRIM28/KAP1. The evolutionary root of modern KRAB domains (mKRAB) can be traced back to an ancestral motif (aKRAB) that occurs even in invertebrates. Here, we first stratified three subgroups of aKRAB sequences from the animal kingdom (PRDM9, SSX and coelacanth KZNF families) and defined ancestral subdomains for KRAB-A and KRAB-B. Using human ZNF10 mKRAB-AB as blueprints for function, we then identified the necessary amino acid changes that transform the inactive aKRAB-A of human PRDM9 into an mKRAB domain capable of mediating silencing and complexing TRIM28/KAP1 in human cells when employed as a hybrid with ZNF10-B. Full gain of function required replacement of residues KR by the conserved motif MLE (positionsA32-A34), which inserted an additional residue, and exchange of A9/S for F, A20/M for L, and A27/R for V. AlphaFold2 modelling documented an evolutionary conserved L-shaped body of two α-helices in all KRAB domains. It is transformed into a characteristic spatial arrangement typical for mKRAB-AB upon the amino acid replacements and in conjunction with a third helix supplied by mKRAB-B. Side-chains pointing outward from the core KRAB 3D structure may reveal a protein-protein interaction code enabling graded binding of TRIM28 to different KRAB domains. Our data provide basic insights into structure-function relationships and emulate transitions of KRAB during evolution.

A microanatomical and histological study of the scales of the Devonian sarcopterygian Miguashaia bureaui and the evolution of the squamation in coelacanths.

Coelacanths have traditionally been described as morphologically conservative throughout their long evolutionary history, which spans more than 400 million years. After an initial burst during the Devonian, a morphological stasis was long thought to have prevailed since the Carboniferous, as shown by the extant Latimeria. New fossil discoveries have challenged this view, with punctual and sometimes unusual departures from the general coelacanth Bauplan. The dermal skeleton is considered to represent one, if not the main, example of morphological stasis in coelacanth evolution and as a consequence, has remained poorly surveyed. The lack of palaeohistological data on the dermoskeleton has resulted in a poor understanding of the early establishment and evolution of the coelacanth squamation. Here we describe the scales of Miguashaia bureaui from the Upper Devonian of Miguasha, Québec (Canada), revealing histological data for a Palaeozoic coelacanth in great detail and adding to our knowledge on the dermal skeleton of sarcopterygians. Miguashaia displays rounded scales ornamented by tubercules and narrow ridges made of dentine and capped with enamel. At least two generations of superimposed odontodes occur, which is reminiscent of the primitive condition of stem osteichthyans like Andreolepis or Lophosteus, and onychodonts like Selenodus. The middle vascular layer is well developed and shows traces of osteonal remodelling. The basal plate consists of a fully mineralised lamellar bone with a repetitive rotation pattern every five layers indicating a twisted plywood-like arrangement of the collagen plies. Comparisons with the extant Latimeria and other extinct taxa show that these features are consistently conserved across coelacanth evolution with only minute changes in certain taxa. The morphological and histological features displayed in the scales of Miguashaia enable us to draw a comprehensive picture of the onset of the coelacanth squamation and to propose and discuss evolutionary scenarios for the coelacanth dermoskeleton.

LEAP2 has antagonized the ghrelin receptor GHSR1a since its emergence in ancient fish.

Recent studies have demonstrated that liver-expressed antimicrobial peptide 2 (LEAP2) antagonizes the ghrelin receptor GHSR1a in mammals. However, its antagonistic function in lower vertebrates has not yet been tested. LEAP2 orthologs have been identified from a variety of fish species; however, previous studies all focused on their antimicrobial activity. To test whether LEAP2 functions as a GHSR1a antagonist in the lowest vertebrates, we studied the antagonism of a fish LEAP2 from Latimeria chalumnae, an extant coelacanth that is one of the closest living fish relatives of tetrapods. Using binding assays, we demonstrated that the coelacanth LEAP2 and ghrelin bound to the coelacanth GHSR1a with IC50 values in the nanomolar range. Using activation assays, we demonstrated that the coelacanth ghrelin activated the coelacanth GHSR1a with an EC50 value in the nanomolar range, and this activation effect was efficiently antagonized by a nanomolar range of the coelacanth LEAP2. In addition, we also showed that the human LEAP2 and ghrelin were as effective as their coelacanth orthologs towards the coelacanth GHSR1a; however, the coelacanth peptides had moderately lower activity towards the human GHSR1a. Thus, LEAP2 serves as an endogenous antagonist of the ghrelin receptor GHSR1a in coelacanth and the ghrelin-LEAP2-GHSR1a system has evolved slowly since its emergence in ancient fish.

Production and Characterization of Recombinant Wild Type Uricase from Indonesian Coelacanth (L. menadoensis) and Improvement of Its Thermostability by In Silico Rational Design and Disulphide Bridges Engineering.

The ideal therapeutic uricase (UOX) is expected to have the following properties; high expression level, high activity, high thermostability, high solubility and low immunogenicity. The latter property is believed to depend largely on sequence identity to the deduced human UOX (dH-UOX). Herein, we explored L. menadoensis uricase (LM-UOX) and found that it has 65% sequence identity to dH-UOX, 68% to the therapeutic chimeric porcine-baboon UOX (PBC) and 70% to the resurrected ancient mammal UOX. To study its biochemical properties, recombinant LM-UOX was produced in E. coli and purified to more than 95% homogeneity. The enzyme had specific activity up to 10.45 unit/mg, which was about 2-fold higher than that of the PBC. One-litre culture yielded purified protein up to 132 mg. Based on homology modelling, we successfully engineered I27C/N289C mutant, which was proven to contain inter-subunit disulphide bridges. The mutant had similar specific activity and production yield to that of wild type (WT) but its thermostability was dramatically improved. Up on storage at -20 °C and 4 °C, the mutant retained ~100% activity for at least 60 days. By keeping at 37 °C, the mutant retained ~100% activity for 15 days, which was 120-fold longer than that of the wild type. Thus, the I27C/N289C mutant has potential to be developed for treatment of hyperuricemia.

The first direct evidence of a Late Devonian coelacanth fish feeding on conodont animals.

We describe the first known occurrence of a Devonian coelacanth specimen from the lower Famennian of the Holy Cross Mountains, Poland, with a conodont element preserved in its digestive tract. A small spiral and phosphatic coprolite (fossil excrement) containing numerous conodont elements and other unrecognized remains was also found in the same deposits. The coprolite is tentatively attributed to the coelacanth. Although it is unclear whether the Late Devonian coelacanth from Poland was an active predator or a scavenger, these finds provide the first direct evidence of feeding on conodont animals by early coelacanth fish, and one of the few evidences of feeding on these animals known to date. It also expands our knowledge about the diet and trophic relations between the Paleozoic marine animals in general.

Conservation of globin genes in the "living fossil" Latimeria chalumnae and reconstruction of the evolution of the vertebrate globin family.

The (hemo-)globins are among the best-investigated proteins in biomedical sciences. These small heme-proteins play an important role in oxygen supply, but may also have other functions. In addition to well known hemoglobin and myoglobin, six other vertebrate globin types have been identified in recent years: neuroglobin, cytoglobin, globin E, globin X, globin Y, and androglobin. Analyses of the genome of the "living fossil" Latimeria chalumnae show that the coelacanth is the only known vertebrate that includes all eight globin types. Thus, Latimeria can also be considered as a "globin fossil". Analyses of gene synteny and phylogenetic reconstructions allow us to trace the evolution and the functional changes of the vertebrate globin family. Neuroglobin and globin X diverged from the other globin types before the separation of Protostomia and Deuterostomia. The cytoglobins, which are unlikely to be involved in O2 supply, form the earliest globin branch within the jawed vertebrates (Gnathostomata), but do not group with the agnathan hemoglobins, as it has been proposed before. There is strong evidence from phylogenetic reconstructions and gene synteny that the eye-specific globin E and muscle-specific myoglobin constitute a common clade, suggesting a similar role in intracellular O2 supply. Latimeria possesses two α- and two ß-hemoglobin chains, of which one α-chain emerged prior to the divergence of Actinopterygii and Sarcopterygii, but has been retained only in the coelacanth. Notably, the embryonic hemoglobin α-chains of Gnathostomata derive from a common ancestor, while the embryonic ß-chains - with the exception of a more complex pattern in the coelacanth and amphibians - display a clade-specific evolution. Globin Y is associated with the hemoglobin gene cluster, but its phylogenetic position is not resolved. Our data show an early divergence of distinct globin types in the vertebrate evolution before the emergence of tetrapods. The subsequent loss of globins in certain taxa may be associated with changes in the oxygen-dependent metabolism. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

One thousand two hundred ninety nuclear genes from a genome-wide survey support lungfishes as the sister group of tetrapods.

The only currently unresolved portion of the backbone phylogeny of the vertebrates involves the relationships among coelacanths, lungfishes, and tetrapods. Despite active research on this question over the past three decades, it is still difficult to determine statistically whether lungfishes alone or both lungfishes and coelacanths together are closely related to tetrapods. To resolve this controversy, we assembled a data set comprising 1,290 nuclear genes encoding 690,838 amino acid residues by analyzing available genome and transcriptome data. Phylogenetic analyses of this data set provided overwhelming evidence that the lungfishes are the closest living relatives of the land vertebrates. This result is strongly supported by high bootstrap values from maximum likelihood and maximum parsimony analyses, Bayesian posterior probabilities of CAT model analysis, and topological tests. Additionally, a species tree analysis without data concatenation also strongly supported this result.

Evolutionary active transposable elements in the genome of the coelacanth.

The apparent morphological stasis in the lineage of the coelacanth, which has been called a "living fossil" by many, has been suggested to be causally related to a slow evolution of its genome, with strongly reduced activity of transposable elements (TEs). Analysis of the African coelacanth showed that at least 25% of its genome is constituted of transposable elements including retrotransposons, endogenous retroviruses and DNA transposons, with a strong predominance of non-Long Terminal Repeat (non-LTR) retrotransposons. The coelacanth genome has been shaped by four major general bursts of transposition during evolution, with major contributions of LINE1, LINE2, CR1, and Deu non-LTR retrotransposons. Many transposable elements are expressed in different tissues and might be active. The number of TE families in coelacanth, but also in lungfish, is lower than in teleost fish, but is higher than in chicken and human. This observation is in agreement with the hypothesis of a sequential elimination of many TE families in the sarcopterygian lineage during evolution. Taken together, our analysis indicates that the coelacanth contains more TE families than birds and mammals, and that these elements have been active during the evolution of the coelacanth lineage. Hence, at the level of transposable element activity, the coelacanth genome does not appear to evolve particularly slowly.

First insight into accumulation of characteristics and tissue distribution of PCBs, PBDEs, and other BFRs in the living Indonesian coelacanth (Latimeria menadoensis).

Persistent organic pollutants, including polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and other brominated flame retardants, were detected in the liver, muscle, and ovary tissues of the Indonesian coelacanth (Latimeria menadoensis) incidentally caught around Gangga Island, North Sulawesi Province, Indonesia, on November 5, 2014. Concentrations of total PCBs (209 congeners, 300-2600 ng g-1 lipid weight) in all tissues showed higher than those of PBDEs (41 congeners, 3.9-6.3 ng g-1 lw) and BTBPE (1.1-3.6 ng g-1 lw). The tissue-specific PCB and PBDE profiles were likely due to differences in the lipid composition. Toxic equivalent (TEQ) values of dioxin-like PCBs in coelacanth tissues were lower than the benchmark values for early-life fish. However, compared with the data reported for deep-sea fishes in the Pacific and Indian Oceans, the relatively high concentrations of PCBs detected in this study raise concerns regarding Indonesian coelacanth conservation and habitat conditions.

Revision of the Middle Triassic coelacanth Ticinepomis Rieppel 1980 (Actinistia, Latimeriidae) with paleobiological and paleoecological considerations.

Coelacanths form today an impoverished clade of sarcopterygian fishes, which were somewhat more diverse during their evolutionary history, especially in the Triassic. Since the first description of the coelacanth Ticinepomis peyeri from the Besano Formation of the UNESCO World Heritage Site of Monte San Giorgio (Canton Ticino, Switzerland), the diversity of coelacanths in the Middle Triassic of this area of the western Paleo-Tethys has been enriched with discoveries of other fossil materials. At Monte San Giorgio, two specimens of Heptanema paradoxum and several specimens of the unusual coelacanth Rieppelia heinzfurreri, have been reported from the Meride Limestone and the Besano Formation, respectively. Another unusual coelacanth, Foreyia maxkuhni, and two specimens referred to Ticinepomis cf. T. peyeri have been described from the isochronous and paleogeographical close Prosanto Formation at the Ducanfurgga and Strel sites (near Davos, Canton Graubünden). In the framework of the revision of the coelacanth material from the Besano Formation kept in the collection of the Paläontologisches Institut und Museum der Universität Zürich (Switzerland), we reviewed the genus Ticinepomis on the basis of the holotype and four new referred specimens. Several morphological traits that were little and/or not understood in T. peyeri are here clarified. We re-evaluate the taxonomic attribution of the material of Ticinepomis cf. T. peyeri from the Prosanto Formation. Morphological characters are different enough from the type species, T. peyeri, to erect a new species, Ticinepomis ducanensis sp. nov., which is shown to be also present in the Besano Formation of Monte San Giorgio, where it is represented by fragmentary bone elements. The recognition of a new coelacanth species indicates that the diversity of this slow-evolving lineage was particularly high in this part of the Western Tethys during the Middle Triassic, especially between 242 and 240 million years ago.

The impact of local genomic properties on the evolutionary fate of genes.

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Mandibular musculature constrains brain-endocast disparity between sarcopterygians.

The transition from water to land by the earliest tetrapods in the Devonian Period is seen as one of the greatest steps in evolution. However, little is understood concerning changes in brain morphology over this transition. Here, we determine the brain-braincase relationship in fishes and basal lissamphibians as a proxy to elucidate the changes that occurred over the fish-tetrapod transition. We investigate six basal extant sarcopterygians spanning coelacanths to salamanders (Latimeria chalumnae, Neoceratodus, Protopterus aethiopicus, P. dolloi, Cynops, Ambystoma mexicanum) using micro-CT and MRI and quantify the brain-braincase relationship in these extant taxa. Our results show that regions of lowest brain-endocast disparity are associated with regions of bony reinforcement directly adjacent to masticatory musculature for the mandible except in Neoceratodus and Latimeria. In Latimeria this deviation from the trend can be accounted for by the possession of an intracranial joint and basicranial muscles, whereas in Neoceratodus difference is attributed to dermal bones contributing to the overall neurocranial reinforcement. Besides Neoceratodus and Latimeria, regions of low brain-endocast disparity occur where there is less reinforcement away from high mandibular muscle mass, where the trigeminal nerve complex exits the braincase and where endolymphatic sacs occupy space between the brain and braincase wall. Despite basal tetrapods possessing reduced adductor muscle mass and a different biting mechanism to piscine sarcopterygians, regions of the neurocranium lacking osteological reinforcement in the basal tetrapods Lethiscus and Brachydectes broadly correspond to regions of high brain-endocast disparity seen in extant taxa.

The hepatic architecture of the coelacanth differs from that of the lungfish in portal triad formation.

The liver architecture of vertebrates can be classified into two types, the portal triad type (having periportal bile ducts) and the non-portal triad type (having non-periportal bile ducts). The former is detectable in the tetrapod liver whereas the lungfish liver has the latter. It remains to be revealed which type of hepatic architecture the coelacanth, which together with the lungfish belongs to the Sarcopterygii, possesses. The present study was undertaken to determine the histological characteristics of the coelacanth liver, and to compare with those of other vertebrates. The coelacanth liver had periportal bile ducts and ductules as detected in mammalian livers. The hepatic artery was found around large portal veins. Hagfish, shark, bichir, sturgeon, bowfin and frog livers had periportal bile ducts and bile ductules, whereas most intrahepatic bile ducts of the lungfish were independent of the distribution of the portal veins as seen in the Otocephala and Euteleostei. The lungfish liver developed duct and ductule structures in the parenchyma. These data indicate that the coelacanth liver had a mammalian-type hepatic architecture with a portal triad, and that the ancestors of tetrapods may have had a portal triad-type liver architecture.

Coelacanth-specific adaptive genes give insights into primitive evolution for water-to-land transition of tetrapods.

Coelacanth is a group of extant lobe-finned fishes in Sarcopterygii that provides evolutionary information for the missing link between ray-finned fish and tetrapod vertebrates. Its phenotypes, different from actinopterygian fishes, have been considered as primitive terrestrial traits such as cartilages in their fatty fins which are homologous with the humerus and femur. To investigate molecular evolution of coelacanth which led to its divergence into Sarcopterygii, we compared its protein coding sequences with 11 actinopterygian fishes. We identified 47 genes under positive selection specific to coelacanth, when compared to Holostei and Teleostei. Out of these, NCDN and 14 genes were associated with spatial learning and nitrogen metabolism, respectively. In homeobox gene superfamily, we identified coelacanth-specific amino acid substitutions, and also observed that one of replacements in SHOX was shared with extant tetrapods. Such molecular changes may cause primordial morphological change in the common ancestor of sarcopterygians. These results suggest that certain genes such as NCDN, MMS19, TRMT1, ALX1, DLX5 and SHOX might have played a role in the evolutionary transition between aquatic and terrestrial vertebrates.

Support for Lungfish as the Closest Relative of Tetrapods by Using Slowly Evolving Ray-Finned Fish as the Outgroup.

In a previous analysis of the phylogenetic relationships of coelacanths, lungfishes and tetrapods, using cartilaginous fish (CF) as the outgroup, the sister relationship of lungfishes and tetrapods was constructed with high statistical support. However, using as the outgroup ray-finned fish (RF), which are more taxonomically closely related to the three lineages than CF, the sister relationship of coelacanths and tetrapods was most often constructed depending on the methods and the data sets, but the statistical support was generally low except in the cases in which the data set including a small number of species was analyzed. In this study, instead of the fast evolving ray-finned fish, teleost fish (TF), in the previous data sets, by using two slowly evolving RF, gar and bowfin, as the outgroup, we showed that the sister relationship of lungfishes and tetrapods was reconstructed with high statistical support. In our analysis the evolutionary rates of gar and bowfin were similar to each other and one third to one half of TF. The difference of the amino acid frequencies of the two species with other lineages was larger than those of TF. This study provides a strong support for lungfishes as the closest relative of tetrapods and indicates the importance of using an appropriate outgroup with small divergence in phylogenetic construction.

A comparative study of piscine defense: The scales of Arapaima gigas, Latimeria chalumnae and Atractosteus spatula.

We compare the characteristics of the armored scales of three large fish, namely the Arapaima gigas (arapaima), Latimeria chalumnae (coelacanth), and Atractosteus spatula (alligator gar), with specific focus on their unique structure-mechanical property relationships and their specialized ability to provide protection from predatory pressures, with the ultimate goal of providing bio-inspiration for manmade materials. The arapaima has flexible and overlapping cycloid scales which consist of a tough Bouligand-type arrangement of collagen layers in the base and a hard external mineralized surface, protecting it from piranha, a predator with extremely sharp teeth. The coelacanth has overlapping elasmoid scales that consist of adjacent Bouligand-type pairs, forming a double-twisted Bouligand-type structure. The collagenous layers are connected by collagen fibril struts which significantly contribute to the energy dissipation, so that the scales have the capability to defend from predators such as sharks. The alligator gar has inflexible articulating ganoid scales made of a hard and highly mineralized enamel-like outer surface and a tough dentine-like bony base, which resist powerful bite forces of self-predation and attack by alligators. The structural differences between the three scales correspond with the attack of their predators, and show refined mechanisms which may be imitated and incorporated into superior bioinspired and biomimetic designs that are specialized to resist specific modes of predation.

Alternative adaptive immunity strategies: coelacanth, cod and shark immunity.

The advent of high throughput sequencing has permitted to investigate the genome and the transcriptome of novel non-model species with unprecedented depth. This technological advance provided a better understanding of the evolution of adaptive immune genes in gnathostomes, revealing several unexpected features in different fish species which are of particular interest. In the present paper, we review the current understanding of the adaptive immune system of the coelacanth, the elephant shark and the Atlantic cod. The study of coelacanth, the only living extant of the long thought to be extinct Sarcopterygian lineage, is fundamental to bring new insights on the evolution of the immune system in higher vertebrates. Surprisingly, coelacanths are the only known jawed vertebrates to lack IgM, whereas two IgD/W loci are present. Cartilaginous fish are of great interest due to their basal position in the vertebrate tree of life; the genome of the elephant shark revealed the lack of several important immune genes related to T cell functions, which suggest the existence of a primordial set of TH1-like cells. Finally, the Atlantic cod lacks a functional major histocompatibility II complex, but balances this evolutionary loss with the expansion of specific gene families, including MHC I, Toll-like receptors and antimicrobial peptides. Overall, these data point out that several fish species present an unconventional adaptive immune system, but the loss of important immune genes is balanced by adaptive evolutionary strategies which still guarantee the establishment of an efficient immune response against the pathogens they have to fight during their life.

Endogenous Retroviruses in Fish Genomes: From Relics of Past Infections to Evolutionary Innovations?

The increasing availability of fish genome sequences has allowed to gain new insights into the diversity and host distribution of retroviruses in fish and other vertebrates. This distribution can be assessed through the identification and analysis of endogenous retroviruses, which are proviral remnants of past infections integrated in genomes. Retroviral sequences are probably important for evolution through their ability to induce rearrangements and to contribute regulatory and coding sequences; they may also protect their host against new infections. We argue that the current mass of genome sequences will soon strongly improve our understanding of retrovirus diversity and evolution in aquatic animals, with the identification of new/re-emerging elements and host resistance genes that restrict their infectivity.