Tracing genome size dynamics in sharks and rays with inclusive sequence analysis by the Squalomix Consortium

Genomes have maintained stable sets of protein-coding genes during evolution, while chromosome organization and genome size vary drastically. Changes in genome size are often attributed to variable amounts of repetitive sequences, including transposable elements. However, it remains poorly understood how such changes were accommodated while maintaining other genomic components. Elasmobranchs, including sharks, rays, and skates, exhibit high among-species variation of genome size and high within-species variation of chromosome length, offering a unique study system to address the question. In this study, we present the first whole genome sequences of the whitebelly skate with remarkably small genome size among elasmobranchs (2.2 Gb), and the red stingray. These chromosome-scale assemblies enabled the assessment of genomic compositions including centromeres and non-coding elements, which revealed notable profiles of tRNA loci and unbiased intragenomic distribution of transposons in elasmobranch genomes. Comparative analyses across these species revealed a shared genomic architecture characterized by correlations of intergenic and intronic sequence lengths with chromosome sizes, with repetitive element accumulation in elongated regions. We analyzed tandemly duplicated genes with high copy number variability. This genome-wide survey revealed the tendency for more frequent tandem gene duplications along with genome size expansion. We document the batoid HoxC cluster in the red stingray genome, which has undergone extensive repetitive element invasion and co-localizes with the HoxB cluster on a sex chromosome. Our study demonstrates an inclusive analysis encompassing both coding and non-coding regions, adaptable to more species in the taxon and a basis for molecular-level understanding on phenotypic diversity of elasmobranchs.