Host population crashes disrupt the diversity of associated marine microbiomes.

Host-associated microbial communities are shaped by myriad factors ranging from host conditions, environmental conditions and other microbes. Disentangling the ecological impact of each of these factors can be particularly difficult as many variables are correlated. Here, we leveraged earthquake-induced changes in host population structure to assess the influence of population crashes on marine microbial ecosystems. A large (7.8 magnitude) earthquake in New Zealand in 2016 led to widespread coastal uplift of up to ~6 m, sufficient to locally extirpate some intertidal southern bull kelp populations. These uplifted populations are slowly recovering, but remain at much lower densities than at nearby, less-uplifted sites. By comparing the microbial communities of the hosts from disturbed and relatively undisturbed populations using 16S rRNA gene amplicon sequencing, we observed that disturbed host populations supported higher functional, taxonomic and phylogenetic microbial beta diversity than non-disturbed host populations. Our findings shed light on microbiome ecological assembly processes, particularly highlighting that large-scale disturbances that affect host populations can dramatically influence microbiome structure. We suggest that disturbance-induced changes in host density limit the dispersal opportunities of microbes, with host community connectivity declining with the density of host populations.

Macroalgal microbiome biogeography is shaped by environmental drivers rather than geographical distance.

BACKGROUND AND AIMS: Contrasting patterns of host and microbiome biogeography can provide insight into the drivers of microbial community assembly. Distance-decay relationships are a classic biogeographical pattern shaped by interactions between selective and non-selective processes. Joint biogeography of microbiomes and their hosts is of increasing interest owing to the potential for microbiome-facilitated adaptation. METHODS: In this study, we examine the coupled biogeography of the model macroalga Durvillaea and its microbiome using a combination of genotyping by sequencing (host) and 16S rRNA amplicon sequencing (microbiome). Alongside these approaches, we use environmental data to characterize the relationship between the microbiome, the host, and the environment. KEY RESULTS: We show that although the host and microbiome exhibit shared biogeographical structure, these arise from different processes, with host biogeography showing classic signs of geographical distance decay, but with the microbiome showing environmental distance decay. Examination of microbial subcommunities, defined by abundance, revealed that the abundance of microbes is linked to environmental selection. As microbes become less common, the dominant ecological processes shift away from selective processes and towards neutral processes. Contrary to expectations, we found that ecological drift does not promote structuring of the microbiome. CONCLUSIONS: Our results suggest that although host macroalgae exhibit a relatively 'typical' biogeographical pattern of declining similarity with increasing geographical distance, the microbiome is more variable and is shaped primarily by environmental conditions. Our findings suggest that the Baas Becking hypothesis of 'everything is everywhere, the environment selects' might be a useful hypothesis to understand the biogeography of macroalgal microbiomes. As environmental conditions change in response to anthropogenic influences, the processes structuring the microbiome of macroalgae might shift, whereas those governing the host biogeography are less likely to change. As a result, increasingly decoupled host-microbe biogeography might be observed in response to such human influences.

Southern Hemisphere coasts are biologically connected by frequent, long-distance rafting events.

Globally, species distributions are shifting in response to environmental change,1 and those that cannot disperse risk extinction.2 Many taxa, including marine species, are showing poleward range shifts as the climate warms.3 In the Southern Hemisphere, however, circumpolar oceanic fronts can present barriers to dispersal.4 Although passive, southward movement of species across this barrier has been considered unlikely,5,6 the recent discovery of buoyant kelp rafts on beaches in Antarctica7,8 demonstrates that such journeys are possible. Rafting is a key process by which diverse taxa-including terrestrial, e.g., Lindo,9 Godinot,10 and Censky et al.,11 and marine, e.g., Carlton et al.12 and Gillespie et al.13 species-can cross oceans.14 Kelp rafts can carry passengers7,15-17 and thus can act as vectors for long-distance dispersal of coastal organisms. The small numbers of kelp rafts previously found in Antarctica7,8 do not, however, shed much light on the frequency of such dispersal events.18 We use a combination of high-resolution phylogenomic analyses (>220,000 SNPs) and oceanographic modeling to show that long-distance biological dispersal events in Southern Ocean are not rare. We document tens of kelp (Durvillaea antarctica) rafting events of thousands of kilometers each, over several decades (1950-2019), with many kelp rafts apparently still reproductively viable. Modeling of dispersal trajectories from genomically inferred source locations shows that distant landmasses are well connected, for example South Georgia and New Zealand, and the Kerguelen Islands and Tasmania. Our findings illustrate the power of genomic approaches to track, and modeling to show frequencies of, long-distance dispersal events.

Commonly used Hardy-Weinberg equilibrium filtering schemes impact population structure inferences using RADseq data.

Reduced representation sequencing (RRS) is a widely used method to assay the diversity of genetic loci across the genome of an organism. The dominant class of RRS approaches assay loci associated with restriction sites within the genome (restriction site associated DNA sequencing, or RADseq). RADseq is frequently applied to non-model organisms since it enables population genetic studies without relying on well-characterized reference genomes. However, RADseq requires the use of many bioinformatic filters to ensure the quality of genotyping calls. These filters can have direct impacts on population genetic inference, and therefore require careful consideration. One widely used filtering approach is the removal of loci that do not conform to expectations of Hardy-Weinberg equilibrium (HWE). Despite being widely used, we show that this filtering approach is rarely described in sufficient detail to enable replication. Furthermore, through analyses of in silico and empirical data sets we show that some of the most widely used HWE filtering approaches dramatically impact inference of population structure. In particular, the removal of loci exhibiting departures from HWE after pooling across samples significantly reduces the degree of inferred population structure within a data set (despite this approach being widely used). Based on these results, we provide recommendations for best practice regarding the implementation of HWE filtering for RADseq data sets.

Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod.

Most animal mitochondrial genomes are small, circular and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long- and short-read DNA sequencing suggest that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

MARES, a replicable pipeline and curated reference database for marine eukaryote metabarcoding.

The use of DNA metabarcoding to characterise the biodiversity of environmental and community samples has exploded in recent years. However, taxonomic inferences from these studies are contingent on the quality and completeness of the sequence reference database used to characterise sample species-composition. In response, studies often develop custom reference databases to improve species assignment. The disadvantage of this approach is that it limits the potential for database re-use, and the transferability of inferences across studies. Here, we present the MARine Eukaryote Species (MARES) reference database for use in marine metabarcoding studies, created using a transparent and reproducible pipeline. MARES includes all COI sequences available in GenBank and BOLD for marine taxa, unified into a single taxonomy. Our pipeline facilitates the curation of sequences, synonymization of taxonomic identifiers used by different repositories, and formatting these data for use in taxonomic assignment tools. Overall, MARES provides a benchmark COI reference database for marine eukaryotes, and a standardised pipeline for (re)producing reference databases enabling integration and fair comparison of marine DNA metabarcoding results.

Testing the advantages and disadvantages of short- and long- read eukaryotic metagenomics using simulated reads.

BACKGROUND: The first step in understanding ecological community diversity and dynamics is quantifying community membership. An increasingly common method for doing so is through metagenomics. Because of the rapidly increasing popularity of this approach, a large number of computational tools and pipelines are available for analysing metagenomic data. However, the majority of these tools have been designed and benchmarked using highly accurate short read data (i.e. Illumina), with few studies benchmarking classification accuracy for long error-prone reads (PacBio or Oxford Nanopore). In addition, few tools have been benchmarked for non-microbial communities. RESULTS: Here we compare simulated long reads from Oxford Nanopore and Pacific Biosciences (PacBio) with high accuracy Illumina read sets to systematically investigate the effects of sequence length and taxon type on classification accuracy for metagenomic data from both microbial and non-microbial communities. We show that very generally, classification accuracy is far lower for non-microbial communities, even at low taxonomic resolution (e.g. family rather than genus). We then show that for two popular taxonomic classifiers, long reads can significantly increase classification accuracy, and this is most pronounced for non-microbial communities. CONCLUSIONS: This work provides insight on the expected accuracy for metagenomic analyses for different taxonomic groups, and establishes the point at which read length becomes more important than error rate for assigning the correct taxon.

Concordant geographic and genetic structure revealed by genotyping-by-sequencing in a New Zealand marine isopod.

Population genetic structure in the marine environment can be influenced by life-history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct-developing species, the New Zealand endemic marine isopod Isocladus armatus. Because our sampling range is intersected by two well-known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation-by-distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large-scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north-south population genetic break occurring between Mahia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation-by-distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta-population detected in the Auckland region.