The Most Developmentally Truncated Fishes Show Extensive Hox Gene Loss and Miniaturized Genomes

The reference genome of Macropodus opercularis (the paradise fish)

Amongst fishes, zebrafish (Danio rerio) has gained popularity as a model system over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species has a highly complex behavioral repertoire and has been the subject of many ethological investigations but lacks genomic resources. Here we report the reference genome assembly of M. opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of 483,077,705 base pairs (~483 Mb) on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ~90% of them to orthogroups.

Comparative genomics reveals insights into anuran genome size evolution

BACKGROUND

Amphibians, particularly anurans, display an enormous variation in genome size. Due to the unavailability of whole genome datasets in the past, the genomic elements and evolutionary causes of anuran genome size variation are poorly understood. To address this, we analyzed whole-genome sequences of 14 anuran species ranging in size from 1.1 to 6.8 Gb. By annotating multiple genomic elements, we investigated the genomic correlates of anuran genome size variation and further examined whether the genome size relates to habitat types.

RESULTS

Our results showed that intron expansions or contraction and Transposable Elements (TEs) diversity do not contribute significantly to genome size variation. However, the recent accumulation of transposable elements (TEs) and the lack of deletion of ancient TEs primarily accounted for the evolution of anuran genome sizes. Our study showed that the abundance and density of simple repeat sequences positively correlate with genome size. Ancestral state reconstruction revealed that genome size exhibits a taxon-specific pattern of evolution, with families Bufonidae and Pipidae experiencing extreme genome expansion and contraction events, respectively. Our result showed no relationship between genome size and habitat types, although large genome-sized species are predominantly found in humid habitats.

CONCLUSIONS

Overall, our study identified the genomic element and their evolutionary dynamics accounting for anuran genome size variation, thus paving a path to a greater understanding of the size evolution of the genome in amphibians.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information available, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprising 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system (immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes) exhibit expansion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum when compared to mammals. Our functional studies suggest that lack of the gc gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi, as well as the genome evolution in cartilaginous fishes.

Small Body, Large Chromosomes: Centric Fusions Shaped the Karyotype of the Amazonian Miniature Fish Nannostomus anduzei (Characiformes, Lebiasinidae)

Miniature refers to species with extraordinarily small adult body size when adult and can be found within all major metazoan groups. It is considered that miniature species have experienced severe alteration of numerous morphological traits during evolution. For a variety of reasons, including severe labor concerns during collecting, chromosomal acquisition, and taxonomic issues, miniature fishes are neglected and understudied. Since some available studies indicate possible relationship between diploid chromosome number (2n) and body size in fishes, we aimed to study one of the smallest Neotropical fish Nannostomus anduzei (Teleostei, Characiformes, Lebiasinidae), using both conventional (Giemsa staining, C-banding) and molecular cytogenetic methods (FISH mapping of rDNAs, microsatellites, and telomeric sequences). Our research revealed that N. anduzei possesses one of the lowest diploid chromosome numbers (2n = 22) among teleost fishes, and its karyotype is entirely composed of large metacentric chromosomes. All chromosomes, except for pair number 11, showed an 18S rDNA signal in the pericentromeric region. 5S rDNA signals were detected in the pericentromeric regions of chromosome pair number 1 and 6, displaying synteny to 18S rDNA signals. Interstitial telomeric sites (ITS) were identified in the centromeric region of pairs 6 and 8, indicating that centric fusions played a significant role in karyotype evolution of studied species. Our study provides further evidence supporting the trend of diploid chromosome number reduction along with miniaturization of adult body size in fishes.

Transposable element and host silencing activity in gigantic genomes

Transposable elements (TEs) and the silencing machinery of their hosts are engaged in a germline arms-race dynamic that shapes TE accumulation and, therefore, genome size. In animal species with extremely large genomes (>10 Gb), TE accumulation has been pushed to the extreme, prompting the question of whether TE silencing also deviates from typical conditions. To address this question, we characterize TE silencing via two pathways-the piRNA pathway and KRAB-ZFP transcriptional repression-in the male and female gonads of Ranodon sibiricus, a salamander species with a ∼21 Gb genome. We quantify 1) genomic TE diversity, 2) TE expression, and 3) small RNA expression and find a significant relationship between the expression of piRNAs and TEs they target for silencing in both ovaries and testes. We also quantified TE silencing pathway gene expression in R. sibiricus and 14 other vertebrates with genome sizes ranging from 1 to 130 Gb and find no association between pathway expression and genome size. Taken together, our results reveal that the gigantic R. sibiricus genome includes at least 19 putatively active TE superfamilies, all of which are targeted by the piRNA pathway in proportion to their expression levels, suggesting comprehensive piRNA-mediated silencing. Testes have higher TE expression than ovaries, suggesting that they may contribute more to the species' high genomic TE load. We posit that apparently conflicting interpretations of TE silencing and genomic gigantism in the literature, as well as the absence of a correlation between TE silencing pathway gene expression and genome size, can be reconciled by considering whether the TE community or the host is currently "on the attack" in the arms race dynamic.

Towards the phylogenetic placement of the enigmatic African genus Prolabeops Schultz, 1941

The small cyprinid genus Prolabeops Schultz, 1941 is restricted to the Nyong and Sanaga River systems in Cameroon. In the past, the genus had been suggested to be either a member of the Labeoninae, Torinae or the Smiliogastrinae mainly on the basis of morphological similarities, and it is nowadays considered as incertae sedis within the Cypriniformes. This study provides the first attempt to reveal the phylogenetic position of Prolabeops using molecular data. For this purpose, the authors sequenced a large fraction of the mitochondrial genome (c. 13,600 bp), including all mitochondrial protein coding genes, of two Prolabeops melanhypopterus specimens and an additional four Enteromius specimens. The large-scale phylogenetic analysis was based on an alignment including all mitochondrial protein coding genes of 902 specimens representing c. 899 cypriniform species. Prolabeops was clearly recovered within the African Smiliogastrinae, forming a weakly supported clade together with Enteromius jae, Enteromius hulstaerti and Barboides gracilis. The study data underline the urgent need of a thorough taxonomic revision of the small African barbs collectively placed in the genus Enteromius.

A complete inventory of long-chain polyunsaturated fatty acid biosynthesis pathway enzymes in the miniaturized cyprinid Paedocypris micromegethes

The capacity for long-chain polyunsaturated fatty acid (LC-PUFA) biosynthesis activity in a species depends on the enzymatic activities of fatty acyl desaturase (Fads) and elongation of very long-chain fatty acid (Elovl). The miniaturized fish Paedocypris micromegethes is a developmentally truncated cyprinid living in highly acidic water conditions in tropical peat swamps. The capacity for LC-PUFA biosynthesis in this species, which has a reduced genome size, is unknown. A high-quality de novo transcriptome assembly enabled the identification of a putative Fads2 and four Elovl. The Fads2 was verified as a P. micromegethes Fads2 ortholog with in vitro Δ5 and Δ6 activities. The Elovl sequences were established as an Elovl5, Elovl2, and two Elovl4 paralogs, namely Elovl4a and Elovl4b. These Elovl enzymes, mainly Elovl5 and Elovl2, fulfill the necessary C18, C20, and C22 PUFA elongation steps for LC-PUFA biosynthesis. Collectively, these results validate the presence of a complete repertoire of LC-PUFA biosynthesis enzymes in a peat swamp miniatured freshwater fish.

The comparative energetics of the ray-finned fish in an evolutionary context

To address challenges in management and conservation of fishes and fisheries it is essential to understand their life histories and energetics. The Add-my-Pet (AmP) collection of data on energetics and Dynamic Energy Budget (DEB) parameters currently contains 1150 of the 40000 extant species of fish. It gives 250-280 traits per species, depending on the model type that was applied, such as maximum reserve capacity, lifespan, specific respiration and precociality index, based on which the ray-finned fish (Actinopterygii) was compared with the four other fish classes (Cyclostomata, Chondrichthyes, Actinistia, Dipnoi) and the Tetrapoda. The Actinopterygii are the only vertebrate class that shows metabolic acceleration, and clearly so in only three sub-clades. Different from chondrichthyans, quite a few species follow the waste-to-hurry strategy, especially small bodied freshwater fish such as tropical annual killifish, but also in small minnows and darters in continental climates. We briefly discuss links between waste-to-hurry, which is associated with a large specific somatic maintenance, and sensitivity for pesticides. We discuss why this interferes with the physical co-variation between maximum reserve capacity and ultimate structural length or weight and explains why maximum reserve capacity increases with body length in chondrichthyans, but not in actinopterygians. Reserve capacity has relevance, e.g. mass-specific maintenance, starvation and the kinetics of lipophyllic compounds (such as pesticides), since reserve is relatively rich in lipids in fish. Also, unlike chondrichthyans, the size at birth is very small and not linked to ultimate size; we discuss the implications. Actinopterygians allocate more to soma, compared with chondrichthyans; the latter allocate more to maturity or reproduction. Actinopterygians, Actinistia and Dipnoi are near the supply-end of the supply-demand spectrum, while chondrichthyans clearly show demand properties.

The little fishes that could: smaller fishes demonstrate slow body size evolution but faster speciation in the family Percidae

Body size affects numerous aspects of organismal biology and many factors have been invoked to explain body size distributions in a macroecological and macroevolutionary context. Body size in the freshwater fish family Percidae is strongly right-skewed (i.e. dominated by small sizes), with small body size potentially being associated with fast water habitats. We constructed a new species-level, multi-locus, time-calibrated phylogeny of Percidae, and used it to test for changes in the rate and pattern of maximum body size evolution. We also tested whether speciation rates varied as a function of body size. We found that Etheostomatinae evolved towards a smaller adaptive optimum in body size compared to the other subfamilies of Percidae, and that this shift was associated with a reduction in the rate of body size evolution. Speciation rates were associated with body size across percids, showing a peak around small to medium body size. Small body size appears to partially, but not fully, explain the diversity of small percids, as many darters fall well below the “optimum” body size. Reinforcement of selection for small body size via selection for novel morphologies or via sexual selection may help to fully explain the remarkable diversity of darter radiation.

The emerging vertebrate model species for neurophysiological studies is Danionella cerebrum, new species (Teleostei: Cyprinidae)

The four described species of Danionella are tiny, transparent fishes that mature at sizes between 10–15 mm, and represent some of the most extreme cases of vertebrate progenesis known to date. The miniature adult size and larval appearance of Danionella, combined with a diverse behavioral repertoire linked to sound production by males, have established Danionella as an important model for neurophysiological studies. The external similarity between the different species of Danionella has offered an important challenge to taxonomic identification using traditional external characters, leading to confusion over the identity of the model species. Using combined morphological and molecular taxonomic approaches, we show here that the most extensively studied species of Danionella is not D. translucida, but represents an undescribed species, D. cerebrum n. sp. that is externally almost identical to D. translucida, but differs trenchantly in several internal characters. Molecular analyses confirm the distinctiveness of D. cerebrum and D. translucida and suggest that the two species are not even sister taxa. Analysis of the evolution of sexual dimorphisms associated with the Weberian apparatus reveals significant increases in complexity from the simpler condition found in D. dracula, to most complex conditions in D. cerebrum, D. mirifica and D. translucida.

Of Traits and Trees: Probabilistic Distances under Continuous Trait Models for Dissecting the Interplay among Phylogeny, Model, and Data

Stochastic models of character trait evolution have become a cornerstone of evolutionary biology in an array of contexts. While probabilistic models have been used extensively for statistical inference, they have largely been ignored for the purpose of measuring distances between phylogeny-aware models. Recent contributions to the problem of phylogenetic distance computation have highlighted the importance of explicitly considering evolutionary model parameters and their impacts on molecular sequence data when quantifying dissimilarity between trees. By comparing two phylogenies in terms of their induced probability distributions that are functions of many model parameters, these distances can be more informative than traditional approaches that rely strictly on differences in topology or branch lengths alone. Currently, however, these approaches are designed for comparing models of nucleotide substitution and gene tree distributions, and thus, are unable to address other classes of traits and associated models that may be of interest to evolutionary biologists. Here, we expand the principles of probabilistic phylogenetic distances to compute tree distances under models of continuous trait evolution along a phylogeny. By explicitly considering both the degree of relatedness among species and the evolutionary processes that collectively give rise to character traits, these distances provide a foundation for comparing models and their predictions, and for quantifying the impacts of assuming one phylogenetic background over another while studying the evolution of a particular trait. We demonstrate the properties of these approaches using theory, simulations, and several empirical data sets that highlight potential uses of probabilistic distances in many scenarios. We also introduce an open-source R package named PRDATR for easy application by the scientific community for computing phylogenetic distances under models of character trait evolution.[Brownian motion; comparative methods; phylogeny; quantitative traits.].

A bird-like genome from a frog: Mechanisms of genome size reduction in the ornate burrowing frog, Platyplectrum ornatum

The diversity of genome sizes across the tree of life is of key interest in evolutionary biology. Various correlates of variation in genome size, such as accumulation of transposable elements (TEs) or rate of DNA gain and loss, are well known, but the underlying molecular mechanisms driving or constraining genome size are poorly understood. Here, we study one of the smallest genomes among frogs characterized thus far, that of the ornate burrowing frog (Platyplectrum ornatum) from Australia, and compare it to other published frog and vertebrate genomes to examine the forces driving reduction in genome size. At ∼1.06 gigabases (Gb), the P. ornatum genome is like that of birds, revealing four major mechanisms underlying TE dynamics: reduced abundance of all major classes of TEs; increased net deletion bias in TEs; drastic reduction in intron lengths; and expansion via gene duplication of the repertoire of TE-suppressing Piwi genes, accompanied by increased expression of Piwi-interacting RNA (piRNA)-based TE-silencing pathway genes in germline cells. Transcriptomes from multiple tissues in both sexes corroborate these results and provide insight into sex-differentiation pathways in Platyplectrum Genome skimming of two closely related frog species (Lechriodus fletcheri and Limnodynastes fletcheri) confirms a reduction in TEs as a major driver of genome reduction in Platyplectrum and supports a macroevolutionary scenario of small genome size in frogs driven by convergence in life history, especially rapid tadpole development and tadpole diet. The P. ornatum genome offers a model for future comparative studies on mechanisms of genome size reduction in amphibians and vertebrates generally.

Insights into skull evolution in fossorial snakes, as revealed by the cranial morphology of Atractaspis irregularis (Serpentes: Colubroidea)

Comparative osteological analyses of extant organisms provide key insight into major evolutionary transitions and phylogenetic hypotheses. This is especially true for snakes, given their unique morphology relative to other squamates and the persistent controversy regarding their evolutionary origins. However, the osteology of several major snake groups remains undescribed, thus hindering efforts to accurately reconstruct the phylogeny of snakes. One such group is the Atractaspididae, a family of fossorial colubroids. We herein present the first detailed description of the atractaspidid skull, based on fully segmented micro-computed tomography (micro-CT) scans of Atractaspis irregularis. The skull of Atractaspis presents a highly unique morphology influenced by both fossoriality and paedomorphosis. This paedomorphosis is especially evident in the jaws, palate, and suspensorium, the major elements associated with macrostomy (large-gaped feeding in snakes). Comparison to scolecophidians-a group of blind, fossorial, miniaturized snakes-in turn sheds light on current hypotheses of snake phylogeny. Features of both the naso-frontal joint and the morphofunctional system related to macrostomy refute the traditional notion that scolecophidians are fundamentally different from alethinophidians (all other extant snakes). Instead, these features support the controversial hypothesis of scolecophidians as "regressed alethinophidians," in contrast to their traditional placement as the earliest-diverging snake lineage. We propose that Atractaspis and scolecophidians fall along a morphological continuum, characterized by differing degrees of paedomorphosis. Altogether, a combination of heterochrony and miniaturization provides a mechanism for the derivation of the scolecophidian skull from an ancestral fossorial alethinophidian morphotype, exemplified by the nonminiaturized and less extreme paedomorph Atractaspis.

Antigen receptor repertoires of one of the smallest known vertebrates

The rules underlying the structure of antigen receptor repertoires are not yet fully defined, despite their enormous importance for the understanding of adaptive immunity. With current technology, the large antigen receptor repertoires of mice and humans cannot be comprehensively studied. To circumvent the problems associated with incomplete sampling, we have studied the immunogenetic features of one of the smallest known vertebrates, the cyprinid fish Paedocypris sp. "Singkep" ("minifish"). Despite its small size, minifish has the key genetic facilities characterizing the principal vertebrate lymphocyte lineages. As described for mammals, the frequency distributions of immunoglobulin and T cell receptor clonotypes exhibit the features of fractal systems, demonstrating that self-similarity is a fundamental property of antigen receptor repertoires of vertebrates, irrespective of body size. Hence, minifish achieve immunocompetence via a few thousand lymphocytes organized in robust scale-free networks, thereby ensuring immune reactivity even when cells are lost or clone sizes fluctuate during immune responses.

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.

Genome of the four-finger threadfin Eleutheronema tetradactylum (Perciforms: Polynemidae)

BACKGROUND

Teleost fish play important roles in aquatic ecosystems and aquaculture. Threadfins (Perciformes: Polynemidae) show a range of interesting biology, and are of considerable importance for both wild fisheries and aquaculture. Additionally, the four-finger threadfin Eleutheronema tetradactylum is of conservation relevance since its populations are considered to be in rapid decline and it is classified as endangered. However, no genomic resources are currently available for the threadfin family Polynemidae.

RESULTS

We sequenced and assembled the first threadfin fish genome, the four-finger threadfin E. tetradactylum. We provide a genome assembly for E. tetradactylum with high contiguity (scaffold N50 = 56.3 kb) and high BUSCO completeness at 96.5%. The assembled genome size of E. tetradactylum is just 610.5 Mb, making it the second smallest perciform genome assembled to date. Just 9.07-10.91% of the genome sequence was found to consist of repetitive elements (standard RepeatMasker analysis vs custom analysis), making this the lowest repeat content identified to date for any perciform fish. A total of 37,683 protein-coding genes were annotated, and we include analyses of developmental transcription factors, including the Hox, ParaHox, and Sox families. MicroRNA genes were also annotated and compared with other chordate lineages, elucidating the gains and losses of chordate microRNAs.

CONCLUSIONS

The four-finger threadfin E. tetradactylum genome presented here represents the first available genome sequence for the ecologically, biologically, and commercially important clade of threadfin fish. Our findings provide a useful genomic resource for future research into the interesting biology and evolution of this valuable group of food fish.

Insights into the Evolution of Neoteny from the Genome of the Asian Icefish Protosalanx chinensis

Salangids, known as Asian icefishes, represent a peculiar radiation within the bony fish order Protacanthopterygii where adult fish retain larval characteristics such as transparent and miniaturized bodies and a cartilaginous endoskeleton into adulthood. Here, we report a de novo genome of Protosalanx chinensis, the most widely distributed salangid lineage. The P. chinensis genome assembly is more contiguous and complete than a previous assembly. We estimate that P. chinensis, salmons, trouts, and pikes diverged from a common ancestor 185 million years ago. A juxtaposition with other fish genomes revealed loss of the genes encoding ectodysplasin-A receptor (EDAR), SCPP1, and four Hox proteins and likely lack of canonical fibroblast growth factor 5 (FGF5) function. We also report genomic variations of P. chinensis possibly reflecting the immune system repertoire of a species with a larval phenotype in sexually mature individuals. The new Asian icefish reference genome provides a solid foundation for future studies.

Extreme Genome and Nervous System Streamlining in the Invertebrate Parasite Intoshia variabili

Orthonectida is a small, rare, and in many aspects enigmatic group of organisms with a unique life cycle and a highly simplified adult free-living stage parasitizing various marine invertebrates [1, 2]. Phylogenetic relationships of Orthonectida have remained controversial for a long time. According to recent data, they are close to Annelida, specifically to Clitellata [3-5]. Several studies have shown that parasitism can not only lead to a dramatic reduction of the body plan and morphological structures but also affect organisms at the genomic level [6, 7]. Comparative studies of parasites and closely related non-parasitic species could clarify the genome reduction degree and evolution of parasitism. Here, we report on the morphology, genome structure, and content of the smallest known Orthonectida species Intoshia variabili, inhabiting the flatworm Graffiellus croceus. This orthonectid with an extremely simplified nervous system demonstrates the smallest known genome (15.3 Mbp) and one of the lowest reported so far gene numbers (5,120 protein-coding genes) among metazoans. The genome is extremely compact, due to a significant reduction of gene number, intergenic regions, intron length, and repetitive elements. The small genome size is probably a result of extreme genome reduction due to their parasitic lifestyle, as well as of simplification and miniaturization of the free-living stages. Our data could provide further insights into the evolution of parasitism and could help to define a minimal bilaterian gene set.

Hybrid genome assembly and annotation of Danionella translucida

Studying neuronal circuits at cellular resolution is very challenging in vertebrates due to the size and optical turbidity of their brains. Danionella translucida, a close relative of zebrafish, was recently introduced as a model organism for investigating neural network interactions in adult individuals. Danionella remains transparent throughout its life, has the smallest known vertebrate brain and possesses a rich repertoire of complex behaviours. Here we sequenced, assembled and annotated the Danionella translucida genome employing a hybrid Illumina/Nanopore read library as well as RNA-seq of embryonic, larval and adult mRNA. We achieved high assembly continuity using low-coverage long-read data and annotated a large fraction of the transcriptome. This dataset will pave the way for molecular research and targeted genetic manipulation of this novel model organism.

Lizards possess the most complete tetrapod Hox gene repertoire despite pervasive structural changes in Hox clusters

Hox genes are a remarkable example of conservation in animal development and their nested expression along the head-to-tail axis orchestrates embryonic patterning. Early in vertebrate history, two duplications led to the emergence of four Hox clusters (A-D) and redundancy within paralog groups has been partially accommodated with gene losses. Here we conduct an inventory of squamate Hox genes using the genomes of 10 lizard and 7 snake species. Although the HoxC1 gene has been hypothesized to be lost in the amniote ancestor, we reveal that it is retained in lizards. In contrast, all snakes lack functional HoxC1 and -D12 genes. Varying levels of degradation suggest differences in the process of gene loss between the two genes. The vertebrate HoxC1 gene is prone to gene loss and its functional domains are more variable than those of other Hox1 genes. We describe for the first time the HoxC1 expression patterns in tetrapods. HoxC1 is broadly expressed during development in the diencephalon, the neural tube, dorsal root ganglia, and limb buds in two lizard species. Our study emphasizes the value of revisiting Hox gene repertoires by densely sampling taxonomic groups and its feasibility owing to growing sequence resources in evaluating gene repertoires across taxa.