The transition of cetaceans (whales, dolphins, and porpoises) from terrestrial to aquatic lifestyles is a striking example of natural selection driving major phenotypic changes (Figure 1). For instance, cetaceans have evolved the ability to withstand high pressure and to store oxygen for long periods, among other adaptations (Das et al. 2023). Many phenotypic changes, such as shifts in organ structure, have been well-characterized through fossils (Thewissen et al. 2009). Although such phenotypic transitions are now well understood, we have only a partial understanding of the underlying genetic mechanisms. Scanning for signatures of adaptation in genes related to phenotypes of interest is one approach to better understand these mechanisms. This was the focus of Uribe and colleagues’ (2024) work, who tested for such signatures across cetacean protein-coding genes.

Figure 1: The skeletons of Ambulocetus (an early whale; top) and Pakicetus (the earliest known cetacean, which lived about 50 million years ago; bottom). Copyright: J. G. M. Thewissen. Displayed here with permission from the copyright holder.

The authors were specifically interested in investigating the evolution of ion channels, as these proteins play fundamental roles in physiological processes. An important aspect of their work was to develop a bioinformatic pipeline to identify orthologous ion channel genes across a set of genomes. After applying their bioinformatic workflow to 18 mammalian species (including nine cetaceans), they conducted tests to find out whether these genes showed signatures of positive selection in the cetacean lineage. For many ion channel genes, elevated ratios of non-synonymous to synonymous substitution rates were detected (for at least a subset of sites, and not necessarily the entire coding region of the genes). The genes concerned were enriched for several functions, including heart and nervous system-related phenotypes.

One top gene hit among the putatively selected genes was SCN5A, which encodes a sodium channel expressed in the heart. Interestingly, the authors noted a specific amino acid replacement, which is associated with sensitivity to the toxin tetrodotoxin in other lineages. This substitution appears to have occurred in the common ancestor of toothed whales, and then was reversed in the ancestor of bottlenose dolphins. The authors describe known bottlenose dolphin interactions with toxin-producing pufferfish that could result in high tetrodotoxin exposure, and thus perhaps higher selection for tetrodotoxin resistance. Although this observation is intriguing, the authors emphasize it requires experimental confirmation.

The authors also recapitulated the previously described observation (Yim et al. 2014; Huelsmann et al. 2019) that cetaceans have fewer protein-coding genes compared to terrestrial mammals, on average. This signal has previously been hypothesized to partially reflect adaptive gene loss. For example, specific gene loss events likely decreased the risk of developing blood clots while diving (Huelsmann et al. 2019). Uribe and colleagues also considered overall gene turnover rate, which encompasses gene copy number variation across lineages, and found the cetacean gene turnover rate to be three times higher than that of terrestrial mammals. Finally, they found that cetaceans have a higher proportion of ion channel genes (relative to all protein-coding genes in a genome) compared to terrestrial mammals. 

Similar investigations of the relative non-synonymous to synonymous substitution rates across cetacean and terrestrial mammal orthologs have been conducted previously, but these have primarily focused on dolphins as the sole cetacean representative (McGowen et al. 2012; Nery et al. 2013; Sun et al. 2013). These projects have also been conducted across a large proportion of orthologous genes, rather than a subset with a particular function. Performing proteome-wide investigations can be valuable in that they summarize the genome-wide signal, but can suffer from a high multiple testing burden. More generally, investigating a more targeted question, such as the extent of positive selection acting on ion channels in this case, or on genes potentially linked to cetaceans’ increased brain sizes (McGowen et al. 2011) or hypoxia tolerance (Tian et al. 2016), can be easier to interpret, as opposed to summarizing broader signals. However, these smaller-scale studies can also experience a high multiple testing burden, especially as similar tests are conducted across numerous studies, which often is not accounted for (Ioannidis 2005). In addition, integrating signals across the entire genome will ultimately be needed given that many genetic changes undoubtedly underlie cetaceans’ phenotypic diversification. As highlighted by the fact that past genome-wide analyses have produced some differing biological interpretations (McGowen et al. 2012; Nery et al. 2013; Sun et al. 2013), this is not a trivial undertaking. 

Nonetheless, the work performed in this preprint, and in related research, is valuable for (at least) three reasons. First, although it is a challenging task, a better understanding of the genetic basis of cetacean phenotypes could have benefits for many aspects of cetacean biology, including conservation efforts. In addition, the remarkable phenotypic shifts in cetaceans make the question of what genetic mechanisms underlie these changes intrinsically interesting to a wide audience. Last, since the cetacean fossil record is especially well-documented (Thewissen et al. 2009), cetaceans represent an appealing system to validate and further develop statistical methods for inferring adaptation from genetic data. Uribe and colleagues’ (2024) analyses provide useful insights relevant to each of these points, and have generated intriguing hypotheses for further investigation.

References

Das K, Sköld H, Lorenz A, Parmentier E (2023) Who are the marine mammals? In: “Marine Mammals: A Deep Dive into the World of Science”. Brennecke D, Knickmeier K, Pawliczka I, Siebert U, Wahlberg, M (editors). Springer, Cham. p. 1–14. https://doi.org/10.1007/978-3-031-06836-2_1

Huelsmann M, Hecker N, Springer MS, Gatesy J, Sharma V, Hiller M (2019) Genes lost during the transition from land to water in cetaceans highlight genomic changes associated with aquatic adaptations. Science Advances, 5, eaaw6671. https://doi.org/10.1126/sciadv.aaw6671

Ioannidis JPA (2005) Why most published research findings are false. PLOS Medicine, 2, e124. https://doi.org/10.1371/journal.pmed.0020124

McGowen MR, Montgomery SH, Clark C, Gatesy J (2011) Phylogeny and adaptive evolution of the brain-development gene microcephalin (MCPH1) in cetaceans. BMC Evolutionary Biology, 11, 98. https://doi.org/10.1186/1471-2148-11-98

McGowen MR, Grossman LI, Wildman DE (2012) Dolphin genome provides evidence for adaptive evolution of nervous system genes and a molecular rate slowdown. Proceedings of the Royal Society B: Biological Sciences, 279, 3643–3651. https://doi.org/10.1098/rspb.2012.0869

Nery MF, González DJ, Opazo JC (2013) How to make a dolphin: molecular signature of positive selection in cetacean genome. PLOS ONE, 8, e65491. https://doi.org/10.1371/journal.pone.0065491

Sun YB, Zhou WP, Liu HQ, Irwin DM, Shen YY, Zhang YP (2013) Genome-wide scans for candidate genes involved in the aquatic adaptation of dolphins. Genome Biology and Evolution, 5, 130–139. https://doi.org/10.1093/gbe/evs123

Tian R, Wang Z, Niu X, Zhou K, Xu S, Yang G (2016) Evolutionary genetics of hypoxia tolerance in cetaceans during diving. Genome Biology and Evolution, 8, 827–839. https://doi.org/10.1093/gbe/evw037

Thewissen JGM, Cooper LN, George JC, Bajpai  (2009) From land to water: the origin of whales, dolphins, and porpoises. Evolution: Education and Outreach, 2, 272–288. https://doi.org/10.1007/s12052-009-0135-2

Uribe C, Nery M, Zavala K, Mardones G, Riadi G, Opazo J (2024) Evolution of ion channels in cetaceans: A natural experiment in the tree of life. bioRxiv, ver. 8 peer-reviewed and recommended by Peer Community in Genomics. https://doi.org/10.1101/2023.06.15.545160

Yim HS, Cho YS, Guang X, Kang SG, Jeong JY, Cha SS, Oh HM, Lee JH, Yang EC, Kwon KK, et al. (2014) Minke whale genome and aquatic adaptation in cetaceans. Nature Genetics, 46, 88–92. https://doi.org/10.1038/ng.2835

Parasitoid wasps lay their eggs inside another arthropod, whose body is physically consumed by the parasitoid larvae. Phylogenetic inference suggests that Parasitoida are monophyletic, and that this clade underwent a strong radiation shortly after branching off from the Apocrita stem, some 236 million years ago (Peters et al. 2017). The increase in taxonomic diversity during evolutionary radiations is usually concurrent with an increase in genetic/genomic diversity, and is often associated with an increase in phenotypic diversity. Gene (or genome) duplication provides the evolutionary potential for such increase of genomic diversity by neo/subfunctionalisation of one of the gene paralogs, and is often proposed to be related to evolutionary radiations (Ohno 1970; Francino 2005).

In their recent preprint, Dominique Colinet and coworkers have explored the genetic and functional diversity of a Rho GTPase activating protein (RhoGAP) multigene family in two very divergent wasp clades within Parasitoida, namely Leptopilina (Figitidae) and Venturia (Ichneumonidae) (Colinet et al. 2024). Some members of the RhoGAP family are present in the venom of the parasitoid wasp Leptopilina boulardi as well as in other Leptopilina species, and are probably involved in the parasitic lifestyle by binding and inactivating host’s Rho GTPases, thereby interfering with the host’s immune response (Colinet et al. 2007).

Venom protein composition is highly variable, even between very closely related species, and is subject to rapid evolutionary changes. Although gene duplication and subsequent neo/subfunctionalisation have been frequently proposed as the main mechanism underlying this evolutionary diversification, observations are often compatible with alternative explanations, such as horizontal gene transfer, gene co-option or multifunctionalisation (Martinson et al. 2017; Alvarado et al. 2020; Huang et al. 2021; Undheim and Jenner 2021). Furthermore, high mutation rates in venom protein-encoding genes hinder phylogenetic hypothesis testing, and venom proteomics can be needed to verify transcriptomic predictions (Smith and Undheim 2018; von Reumont et al. 2022).

Colinet and coworkers (2024) have applied a combined transcriptomic, proteomic and functional approach to i) identify potential transcripts of the RhoGAP family in Leptopilina species using experimental and bioinformatic approaches; ii) experimentally identify proteins of the RhoGAP family in the venom of three Leptopilina species; iii) identify transcripts and proteins of the RhoGAP family in the ovarian calyx of Venturia canescens; and iv) perform phylogenetic and selection analyses on the extant sequences of these RhoGAP family genes to propose an evolutionary scenario for their origin and diversification. The most striking results are first the large diversity of RhoGAP sequences retrieved in the transcriptomes and proteomes of Leptopilina and of V. canescens, and second the high number of branches and positions identified to have evolved under positive selection. All the retrieved hits share a RhoGAP domain, either alone or in tandem, preceded in the case of Leptopilina RhoGAPs by a signal peptide that may be responsible for protein vehiculation for venom secretion. Further, for some of the protein positions identified to have evolved under positive selection, the authors have experimentally verified the functional impact of the changes by reverse genetic engineering.

The authors propose an evolutionary scenario to interpret the phylogenetic relationships among extant RhoGAP diversity in the clades under study. They posit that two independent, incomplete duplication events from the respectively ancestral RacGAP gene, followed by subsequent, lineage- and paralog-specific duplication events, lie at the origin of the wealth of diversity of in the Leptopilina venom RhoGAPs and of V. canescens ovarian calyx RhoGAPs. Notwithstanding, the global relationships presented in the work are not systematically consistent with this interpretation, e.g. regarding the absence of monophyly for Leptopilina RhoGAPs and Leptopilina RacGAP, and the same holds true for the respective V. canescens sequences. It may very well be that the high evolutionary rate of these genes has eroded the phylogenetic signal and prevented proper reconstruction, as the large differences between codon-based and amino acid-based phylogenies and the low support suggest. Explicit hypothesis testing, together with additional data from other taxa, may shed light onto the evolution of this gene family.

The work by Colinet and coworkers communicates sound, novel transcriptomic, proteomic and functional data from complex gene targets, consolidated from an important amount of experimental and bioinformatic work, and related to evolutionarily intriguing and complex phenotypes. These results, and the evolutionary hypothesis proposed to account for them, will be instrumental for our understanding of the evolution and diversity of vesicle-associated RhoGAPs in divergent parasitoid wasps.

References

Alvarado, G., Holland, S., R., DePerez-Rasmussen, J., Jarvis, B., A., Telander, T., Wagner, N., Waring, A., L., Anast, A., Davis, B., Frank, A., et al. (2020). Bioinformatic analysis suggests potential mechanisms underlying parasitoid venom evolution and function. Genomics 112(2), 1096–1104. https://doi.org/10.1016/j.ygeno.2019.06.022

Colinet, D., Cavigliasso, F., Leobold, M., Pichon, A., Urbach, S., Cazes, D., Poullet, M., Belghazi, M., Volkoff, A-N., Drezen, J-M., Gatti, J-L., and Poirié, M. (2024). Convergent origin and accelerated evolution of vesicle-associated RhoGAP proteins in two unrelated parasitoid wasps. bioRxiv, ver. 3 peer-reviewed and recommended by Peer Community in Genomics. https://doi.org/10.1101/2023.06.05.543686

Colinet, D., Schmitz, A., Depoix, D., Crochard, D., and Poirié, M. (2007). Convergent Use of RhoGAP Toxins by eukaryotic parasites and bacterial pathogens. PLoS Pathogens 3(12), e203. https://doi.org/10.1371/journal.ppat.0030203

Francino, M.P. (2005). An adaptive radiation model for the origin of new gene functions. Nature Genetics 37, 573–577. https://doi.org/10.1038/ng1579

Huang, J., Chen, J., Fang, G., Pang, L., Zhou, S., Zhou, Y., Pan, Z., Zhang, Q., Sheng, Y., Lu, Y., et al. (2021). Two novel venom proteins underlie divergent parasitic strategies between a generalist and a specialist parasite. Nature Communications 12, 234. https://doi.org/10.1038/s41467-020-20332-8

Martinson, E., O., Mrinalini, Kelkar, Y. D., Chang, C-H., and Werren, J., H. 2017. The evolution of venom by co-option of single-copy genes. Current Biololgy 27(13), 2007-2013.e8. https://doi.org/10.1016/j.cub.2017.05.032

Ohno, S. (1970). Evolution by gene duplication. New-York: Springer-Verlag.

Peters, R., S., Krogmann, L., Mayer, C., Donath, A., Gunkel, S., Meusemann, K., Kozlov, A., Podsiadlowski, L., Petersen, M., Lanfear, R., et al. (2017). Evolutionary history of the Hymenoptera. Current Biology 27(7), 1013–1018. https://doi.org/10.1016/j.cub.2017.01.027

von Reumont, B., M., Anderluh, G., Antunes, A., Ayvazyan, N., Beis, D., Caliskan, F., Crnković, A., Damm, M., Dutertre, S., Ellgaard, L., et al. (2022). Modern venomics—Current insights, novel methods, and future perspectives in biological and applied animal venom research. GigaScience 11, giac048. https://doi.org/10.1093/gigascience/giac048

Smith, J., J., and Undheim, E., A., B. (2018). True lies: using proteomics to assess the accuracy of transcriptome-based venomics in centipedes uncovers false positives and reveals startling intraspecific variation in Scolopendra subspinipes. Toxins 10(3), 96. https://doi.org/10.3390/toxins10030096

Undheim, E., A., B., and Jenner, R., A. (2021). Phylogenetic analyses suggest centipede venom arsenals were repeatedly stocked by horizontal gene transfer. Nature Communications 12, 818. https://doi.org/10.1038/s41467-021-21093-8

High-quality genomes are currently being generated at an unprecedented speed powered by long-read sequencing technologies. However, sequencing effort is concentrated unequally across the tree of life and several key evolutionary and ecological groups remain largely unexplored. So is the case for fish species of the family Scorpaenidae (Perciformes). Kitsoulis et al. present the genome of the devil firefish, Pterois miles (1). Following current best practices, the assembly relies largely on Oxford Nanopore long reads, aided by Illumina short reads for polishing to increase the per-base accuracy. PacBio’s IsoSeq was used to sequence RNA from a variety of tissues as direct evidence for annotating genes. The reconstructed genome is 902 Mb in size and has high contiguity (N50=14.5 Mb; 660 scaffolds, 90% of the genome covered by the 83 longest scaffolds) and completeness (98% BUSCO completeness). The new genome is used to assess the phylogenetic position of P. miles, explore gene synteny against zebrafish, look at orthogroup expansion and contraction patterns in Perciformes, as well as to investigate the evolution of toxins in scorpaenid fish (2). In addition to its value for better understanding the evolution of scorpaenid and teleost fishes, this new genome is also an important resource for monitoring its invasiveness through the Mediterranean Sea (3) and the Atlantic Ocean, in the latter case forming the invasive lionfish complex with P. volitans (4).

REFERENCES

1. Kitsoulis CV, Papadogiannis V, Kristoffersen JB, Kaitetzidou E, Sterioti E, Tsigenopoulos CS, Manousaki T. (2023) Near-chromosome level genome assembly of devil firefish, Pterois miles. BioRxiv, ver. 6 peer-reviewed and recommended by Peer Community in Genomics. https://doi.org/10.1101/2023.01.10.523469

2. Kiriake A, Shiomi K. (2011) Some properties and cDNA cloning of proteinaceous toxins from two species of lionfish (Pterois antennata and Pterois volitans). Toxicon, 58(6-7):494–501.  https://doi.org/10.1016/j.toxicon.2011.08.010

3. Katsanevakis S, et al. (2020) Un- published Mediterranean records of marine alien and cryptogenic species. BioInvasions Records, 9:165–182. https://doi.org/10.3391/bir.2020.9.2.01

4. Lyons TJ, Tuckett QM, Hill JE. (2019) Data quality and quantity for invasive species: A case study of the lionfishes. Fish and Fisheries, 20:748–759. https://doi.org/10.1111/faf.12374

Since the opening of the Rhine-Main-Danube-Channel, four goby species are known to have invaded the River Rhine. Of these, the most recent and numerous is the Caucasian Dwarf Goby, which has been found in the Rhine since 2019. This study presents a new high-quality genome for this species (Knipowitschia cf. caucasica) (Schoenle et al. 2024). Currently, chromosome-scale genome assemblies represent a key first step in invasion biology, allowing the reconstruction of a species’ invasion history and monitoring its progress, as well as identifying and characterizing candidate genes that control invasiveness (McCartney et al. 2019). 

The authors sequenced the nuclear and mitochondrial genomes of this species using state-of-the art methods including long-read sequencing techniques, scaffolded based on chromatin conformation data, and annotated using both direct transcriptomic and protein homology evidence. Data analyses follow currently established pipelines for genome assembly, scaffolding, annotation, and downstream bioinformatic analyses. The quality of the final genome was thoroughly assessed and conforms to what is expected from other genomes of fishes in the family Gobiidae. This study follows other recent endeavors that generated high-quality genomes to improve our understanding of invasion biology (e.g. Shao et al. 2020 and Kitsoulis et al. 2023). These studies are successfully contributing to increasing the genomic resources for the world’s most damaging invasive species, which were not available for even a third of the top 100 invasive species just five years ago (McCarthy et al. 2019). Beyond invasion biology, the Dwarf Goby genome is also an important resource for many other applications, including evolutionary genomic analyses and phylogeography of this species and closely related ones in their native ranges.

                
References          

Kitsoulis CV, Papadogiannis V, Kristoffersen JB, Kaitetzidou E, Sterioti A, Tsigenopoulos CS, Manousaki T (2023) Near-chromosome level genome assembly of devil firefish, Pterois miles. Peer Community Journal 3:e64. https://doi.org/10.24072/pcjournal.295  

McCartney MA, Mallez S, Gohl DM (2019) Genome projects in invasion biology. Conservation Genetics 20:1201–1222. https://doi.org/10.1007/s10592-019-01224-x

Schoenle A, Guiglielmoni N, Mainz T, Greve C, Hamadou AB, Heermann L, Borcherding J, Waldvogel A-M (2024) Chromosome level genome reference of the Caucasian dwarf goby Knipowitschia cf. caucasica, a new alien Gobiidae invading the River Rhine. bioRxiv, ver. 3 peer-reviewed and recommended by PCI Genomics. https://doi.org/10.1101/2024.04.22.590508

Shao F, Ludwig A, Mao Y, Liu N, Peng Z (2020). Chromosome-level genome assembly of the female western mosquitofish (Gambusia affinis). GigaScience 9:giaa092. https://doi.org/10.1093/gigascience/giaa092

The assessment by Brézellec (2024) of the quality of HHpred-based SARS-CoV-2 protein annotations against the traditional Pfam annotations is highly justified and valuable. HHpred’s ability to detect remote homologies offers an expanded view of viral protein similarities, potentially uncovering subtle functional mimicries that Pfam may miss due to its sensitivity limitations when dealing with divergent sequences. However, the accuracy and specificity of HHpred results can be compromised by false positives, especially when dealing with complex viral proteins that feature transmembrane or low-complexity regions prone to spurious matches.

To address this, the author made a thoughtful decision to implement a multi-step validation protocol. This approach included establishing progressively lower probability thresholds to capture weaker but biologically plausible hits, and organizing hits into “families” of similarly located alignments to validate the robustness of matches. They also cross-verified results by running SARS-CoV-2 protein queries against non-human proteomes (plants, fruit flies, bacteria, and archaea), allowing them to discern between biologically meaningful matches and potentially random alignments. By adding manual verification with InterPro domain annotations, the authors took additional steps to ensure that identified similarities were not only statistically significant but also biologically relevant.

This rigorous validation strategy adds a layer of reliability to HHpred results, demonstrating an effective maximization of sensitivity while maintaining specificity. This approach yielded biologically intriguing and previously undocumented similarities, such as between the Spike-prominin and ORF3a-GPCR, underscoring the quality and depth of the annotation process. These findings highlight a pathway for further experimental validation and illustrate the potential of HHpred to contribute high-quality insights when applied with careful quality control measures.

In summary, the decision to adopt HHpred (Gabler et al. 2020) and enhance its outputs with a robust quality validation process not only improved the depth of SARS-CoV-2 protein annotations but also established a high standard for future viral annotation projects, striking an effective balance between discovery potential and annotation quality​. The authors have conducted a study that is methodologically rigorous, well-detailed, and highly pertinent to the field. This work stands as a significant contribution to the scientific community, providing resources and insights that are likely to guide future research in this area. 

              
References

Brézellec, P (2024) Re-annotation of SARS-CoV-2 proteins using an HHpred-based approach opens new opportunities for a better understanding of this virus. bioRxiv, ver. 3 peer-reviewed and recommended by PCI Genomics. https://doi.org/10.1101/2023.06.06.543855

Gabler F, Nam S-Z, Till S, Mirdita M, Steinegger M, Söding J, Lupas AN, Alva V (2020) Protein Sequence Analysis Using the MPI Bioinformatics Toolkit. Current Protocols in Bioinformatics, 72, e108. https://doi.org/10.1002/cpbi.108