Contrasting gene flow at different spatial scales revealed by genotyping-by-sequencing in Isocladus armatus, a massively colour polymorphic New Zealand marine isopod

Understanding how genetic diversity is maintained within populations is central to evolutionary biology. Research on colour polymorphism (CP), which typically has a genetic basis, can shed light on this issue. However, because gene flow can homogenise genetic variation, understanding population connectivity is critical in examining the maintenance of polymorphisms. In this study we assess the utility of genotyping-by-sequencing to resolve gene flow, and provide a preliminary investigation into the genetic basis of CP in Isocladus armatus, an endemic New Zealand marine isopod. Analysis of the genetic variation in 4,000 single nucleotide polymorphisms (SNPs) within and among populations and colour morphs revealed large differences in gene flow across two spatial scales. Marine isopods, which lack a pelagic larval phase, are typically assumed to exhibit greater population structuring than marine invertebrates possessing a biphasic life cycle. However, we found high gene flow rates and no genetic subdivision between two North Island populations situated 8 km apart. This suggests that I. armatus is capable of substantial dispersal along coastlines. In contrast, we identified a strong genetic disjunction between North and South Island populations. This result is similar to those reported in other New Zealand marine species, and is congruent with the presence of a geophysical barrier to dispersal down the east coast of New Zealand. We also found some support for a genetic basis to colouration evidenced by positive F(ST) outlier tests, with two SNPs in particular showing strong association to the expression of a striped morph. Our study provides one of the first population genomic studies of a marine organism in New Zealand, and suggests that genotyping-by-sequencing can be a good alternative to more traditional investigations based on traditional markers such as microsatellites. Our study provides a foundation for further development of a highly tractable system for research on the evolutionary maintenance of CP.