Natal Origin Identification of Green Turtles in the North Pacific by Genome-Wide Population Analysis With Limited DNA Samples

Although studies using genome-wide single nucleotide polymorphisms (SNPs) are growing in non-model species, it is still difficult to prepare sufficient high-quality genomic DNA (gDNA) for population genomic analyses in many wild species. In this study, we first analyzed the population structure of green turtles in the North Pacific; four regions in Japan (the Yaeyama, Okinawa, Amami, and Ogasawara Islands), the Southeast Asia, Taiwan, the Federated States of Micronesia, the Republic Marshall Islands, and the Central or Eastern Pacific (C/E Pacific) using ≃1 ng of gDNA from green turtle field samples – including ones from dead carcasses – and multiplexed inter-simple sequence repeat (ISSR) genotyping by sequencing (MIG-seq). In addition, we explored whether the genome-wide SNP data narrowed down the natal origins of foraging turtles that were difficult to identify by mitochondrial DNA (mtDNA) haplotype alone. The overall nesting population structure observed from genome-wide SNPs was consistent with results from mtDNA, and the population isolation of the C/E Pacific and Ogasawara Islands from the other North Pacific populations was highlighted. However, the population boundaries among the Northwestern Pacific, other than those of the Ogasawaras, were not clear. The uniqueness of the Ogasawara population in genome-wide SNP data enabled the identification of turtles that were more likely to have been born on the Ogasawara Islands. Our results show that genome-wide SNP data are more practical in identifying the natal origins of turtles that were difficult determine by the conventional mtDNA-basis mixed stock analysis. This study also showed that MIG-seq can be expected to meet the potential demand to utilize many preserved or fragmented gDNA samples for population genomics in many marine megafaunas for which sample collections are difficult.