Parallel recolonizations generate distinct genomic sectors in kelp following high-magnitude earthquake disturbance

Large-scale disturbance events have the potential to drastically reshape biodiversity patterns. Notably, newly vacant habitat space cleared by disturbance can be colonized by multiple lineages, which can lead to the evolution of distinct spatial "sectors" of genetic diversity within a species. We test for disturbance-driven sectoring of genetic diversity in intertidal southern bull kelp, Durvillaea antarctica (Chamisso) Hariot, following the high-magnitude 1855 Wairarapa earthquake in New Zealand. Specifically, we use genotyping-by-sequencing (GBS) to analyse fine-scale population structure across the uplift zone and apply machine learning to assess the fit of alternative recolonizaton models. Our analysis reveals that specimens from the uplift zone carry distinctive genomic signatures potentially linked to post-earthquake recolonization processes. Specifically, our analysis identifies two parapatric spatial-genomic sectors of D. antarctica at Turakirae Head, which experienced the most dramatic uplift. Based on phylogeographical modelling, we infer that bull kelp in the Wellington region was probably a source for recolonization of the heavily uplifted Turakirae Head coastline, via two parallel, eastward recolonization events. By identifying multiple parapatric genotypic sectors within a recently recolonized coastal region, the current study provides support for the hypothesis that competing lineage expansions can generate striking spatial structuring of genetic diversity, even in highly dispersive taxa.

Application of palaeogenetic techniques to historic mollusc shells reveals phylogeographic structure in a New Zealand abalone

Natural history collections worldwide contain a plethora of mollusc shells. Recent studies have detailed the sequencing of DNA extracted from shells up to thousands of years old and from various taphonomic and preservational contexts. However, previous approaches have largely addressed methodological rather than evolutionary research questions. Here we report the generation of DNA sequence data from mollusc shells using such techniques, applied to Haliotis virginea Gmelin, 1791, a New Zealand abalone, in which morphological variation has led to the recognition of several forms and subspecies. We successfully recovered near-complete mitogenomes from 22 specimens including 12 dry-preserved shells up to 60 years old. We used a combination of palaeogenetic techniques that have not previously been applied to shell, including DNA extraction optimized for ultra-short fragments and hybridization-capture of single-stranded DNA libraries. Phylogenetic analyses revealed three major, well-supported clades comprising samples from: 1) the Three Kings Islands; 2) the Auckland, Chatham and Antipodes Islands; and 3) mainland New Zealand and Campbell Island. This phylogeographic structure does not correspond to the currently recognized forms. Critically, our non-reliance on freshly collected or ethanol-preserved samples enabled inclusion of topotypes of all recognized subspecies as well as additional difficult-to-sample populations. Broader application of these comparatively cost-effective and reliable methods to modern, historical, archaeological and palaeontological shell samples has the potential to revolutionize invertebrate genetic research.

Population structure of the New Zealand whelk, Cominella glandiformis (Gastropoda: Buccinidae), suggests sporadic dispersal of a direct developer

We examined phylogeographic structure in the direct-developing New Zealand endemic intertidal mud whelk, Cominella glandiformis. Two hundred and ninety-six whelks from 12 sites were collected from sheltered shores around New Zealand’s four largest islands (North Island, South Island, Stewart Island and Chatham Island), encompassing the geographical range of this species. Despite being direct developers, gene flow among C. glandiformis populations may occur over short distances by adult floating, and over larger distances by rafting of egg masses. Primers were developed to amplify variable microsatellite regions at six loci. All loci were variable, with 8–34 alleles/loci. Observed and expected heterozygosities were high across all alleles, with minimal evidence of null alleles. The average number of alleles varied from 3.5 (Chatham Island) to 7.5 (Waitemata Harbour). Strong genetic structure was evident, with distinct ‘eastern’ and ‘western’ groups. Each group extended over a large geographic area, including regions of unsuitable habitat, but were linked by oceanic currents. We suggest that the intraspecific geographic genetic structure in C. glandiformis has arisen due a combination of ocean currents (promoting gene flow between geographically distant regions) and upwelling areas (limiting gene flow between certain regions).