Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management attains its sustainability goals. Here we generated 55 Gb of short Illumina reads (92x coverage) and 73 Gb of long Nanopore reads (122x) to produce the first genome assembly of the marine teleost tarakihi (Nemadactylus macropterus), a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq RNA reads were obtained from four tissue types of another specimen to assist in gene annotation. The final genome assembly was 568 Mb long and consisted of 1,214 scaffolds with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii BUSCOs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared to the same values reported for other fishes. Repetitive elements covered 30.45% of the genome and 20,169 protein-coding genes were annotated. Iso-Seq analysis recovered 91,313 unique transcripts (isoforms) from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly along with the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics, as well as comparative eco-evolutionary studies in other teleost and related organisms.