Comparative genomics of tarakihi (Nemadactylus macropterus) and five New Zealand fish species: assembly contiguity affects the identification of genic features but not transposable elements

Comparative analysis of whole-genome sequences can provide valuable insights into the evolutionary patterns of diversification and adaptation of species, including the genome contents and the regions under selection. However, such studies are lacking for fishes in New Zealand. To supplement the recently sequenced genome of tarakihi (Nemadactylus macropterus), the genomes of five additional percomorph species native to New Zealand (king tarakihi (Nemadactylus n.sp.), blue moki (Latridopsis ciliaris), butterfish (Odax pullus), barracouta (Thyrsites atun), and kahawai (Arripis trutta)) were determined and assembled using Illumina sequencing. While the proportion of repeat elements was highly correlated with the genome size (R2 = 0.97, P < 0.01), most of the metrics for the genic features (e.g. number of exons or intron length) were significantly correlated with assembly contiguity (| R2| = 0.79-0.97). A phylogenomic tree including eight additional high-quality fish genomes was reconstructed from sequences of shared gene families. The radiation of Percomorpha was estimated to have occurred c. 112 mya (mid-Cretaceous), while the Latridae have diverged from true Perciformes c. 83 mya (late Cretaceous). Evidence of positive selection was found in 65 genes in tarakihi and 209 genes in Latridae: the largest portion of these are involved in the ATP binding pathway and the integral structure of membranes. These results and the de novo genome sequences can be used to (1) inform future studies on both the strength and shortcomings of scaffold-level assemblies for comparative genomics and (2) provide insights into the evolutionary patterns and processes of genome evolution in bony fishes.