Shallow gene pools in the high intertidal: extreme loss of genetic diversity in viviparous sea stars (Parvulastra).

We document an extreme example of reproductive trait evolution that affects population genetic structure in sister species of Parvulastra cushion stars from Australia. Self-fertilization by hermaphroditic adults and brood protection of benthic larvae causes strong inbreeding and range-wide genetic poverty. Most samples were fixed for a single allele at nearly all nuclear loci; heterozygotes were extremely rare (0.18%); mitochondrial DNA sequences were more variable, but few populations shared haplotypes in common. Isolation-with-migration models suggest that these patterns are caused by population bottlenecks (relative to ancestral population size) and low gene flow. Loss of genetic diversity and low potential for dispersal between high-intertidal habitats may have dire consequences for extinction risk and potential for future adaptive evolution in response to climate and other selective agents.

Extraordinarily rapid life-history divergence between Cryptasterina sea star species.

Life history plays a critical role in governing microevolutionary processes such as gene flow and adaptation, as well as macroevolutionary processes such speciation. Here, we use multilocus phylogeographic analyses to examine a speciation event involving spectacular life-history differences between sister species of sea stars. Cryptasterina hystera has evolved a suite of derived life-history traits (including internal self-fertilization and brood protection) that differ from its sister species Cryptasterina pentagona, a gonochoric broadcast spawner. We show that these species have only been reproductively isolated for approximately 6000 years (95% highest posterior density of 905-22 628), and that this life-history change may be responsible for dramatic genetic consequences, including low nucleotide diversity, zero heterozygosity and no gene flow. The rapid divergence of these species rules out some mechanisms of isolation such as adaptation to microhabitats in sympatry, or slow divergence by genetic drift during prolonged isolation. We hypothesize that the large phenotypic differences between species relative to the short divergence time suggests that the life-history differences observed may be direct responses to disruptive selection between populations. We speculate that local environmental or demographic differences at the southern range margin are possible mechanisms of selection driving one of the fastest known marine speciation events.

Trans-Arctic vicariance in Strongylocentrotus sea urchins.

The sea urchins Strongylocentotus pallidus and S. droebachiensis first invaded the Atlantic Ocean from the Pacific following the opening of the Bering seaway in the late Miocene. While trans-Arctic dispersal during the Pleistocene is thought to have maintained species' integrity, a recent genomic analysis identified a reproductively isolated cryptic species within S. droebachiensis. Based on previous studies, the distribution of one of these lineages (S. droebachiensis W) includes the shallow water habitats of the northwest Atlantic and Pacific, while the other (S. droebachiensis E) is found throughout the shallow habitat in the northeast but is mostly restricted to deep habitats (>65 m) in the northwest Atlantic. However, since genetic variation within S. droebachiensis has been largely unstudied in the north Pacific and Arctic oceans, the biogeography of the cryptic species is not well known, and it is difficult to identify the mechanisms driving population subdivision and speciation. Here we use population genetic analyses to characterize the distribution of each species, and to test hypotheses about the role of vicariance in the evolution of systematic and genomic divergence within the genus. We collected individuals of all three Strongylocentrotus species (n = 365) from 10 previously unsampled locations in the northeast Pacific and north Atlantic (Labrador Sea and Norway), and generated mtDNA sequence data for a 418 bp fragment of cytochrome c oxidase subunit I (COI). To assess the biogeography of all three species, we combined our alignment with five previously published data sets (total n = 789) and used statistical parsimony and maximum likelihood to identify species and characterize their distribution within and among oceans. Patterns of haplotype sharing, pairwise F ST , and hierarchical analyses of molecular variance (AMOVA) identified trans-Arctic dispersal in S. pallidus and S. droebachiensis W, but other than 5 previously reported singletons we failed to detect additional mtDNA haplotypes of S. droebachiensis E in the north Pacific. Within the Atlantic, patterns of habitat segregation suggests that temperature may play a role in limiting the distribution of S. droebachiensis E, particularly throughout the warmer coastal waters along the coast of Nova Scotia. Our results are consistent with the cycles of trans-Arctic dispersal and vicariance in S. pallidus and S. droebachiensis W, but we suggest that the evolution of Atlantic populations of S. droebachiensis E has been driven by persistent trans-Arctic vicariance that may date to the initial invasion in the late Pliocene.

Historical human activities reshape evolutionary trajectories across both native and introduced ranges.

The same vectors that introduce species to new ranges could move them among native populations, but how human-mediated dispersal impacts native ranges has been difficult to address because human-mediated dispersal and natural dispersal can simultaneously shape patterns of gene flow. Here, we disentangle human-mediated dispersal from natural dispersal by exploiting a system where the primary vector was once extensive but has since ceased. From 10th to 19th Centuries, ships in the North Atlantic exchanged sediments dredged from the intertidal for ballast, which ended when seawater ballast tanks were adopted. We investigate genetic patterns from RADseq-derived SNPs in the amphipod Corophium volutator (n = 121; 4,870 SNPs) and the annelid Hediste diversicolor (n = 78; 3,820 SNPs), which were introduced from Europe to North America, have limited natural dispersal capabilities, are abundant in intertidal sediments, but not commonly found in modern water ballast tanks. We detect similar levels of genetic subdivision among introduced North American populations and among native European populations. Phylogenetic networks and clustering analyses reveal population structure between sites, a high degree of phylogenetic reticulation within ranges, and phylogenetic splits between European and North American populations. These patterns are inconsistent with phylogeographic structure expected to arise from natural dispersal alone, suggesting human activity eroded ancestral phylogeographic structure between native populations, but was insufficient to overcome divergent processes between naturalized populations and their sources. Our results suggest human activity may alter species' evolutionary trajectories on a broad geographic scale via regional homogenization and global diversification, in some cases precluding historical inference from genetic data.

Multiple gene genealogies reveal asymmetrical hybridization and introgression among strongylocentrotid sea urchins.

The evolution of incompatibilities between eggs and sperm is thought to play important roles in establishing and maintaining reproductive isolation among species of broadcast-spawning marine invertebrates. However, the effectiveness of gametic isolation in initiating the speciation process and/or in limiting the introgression of genes among species at later stages of divergence remains largely unknown. In the present study, we collected DNA sequence data from five loci in four species of Strongylocentrotus sea urchins (S. droebachiensis, S. pallidus, S. purpuratus, and S. franciscanus) to test whether the susceptibility of S. droebachiensis eggs to fertilization by heterospecific sperm results in gene flow between species. Despite the potential for introgression, a small but statistically significant signal of introgression was observed only between the youngest pair of sister taxa (S. pallidus and S. droebachiensis) that was strongly asymmetrical (from the former into the latter). No significant gene flow was observed for either S. purpuratus or S. franciscanus despite the ability of their sperm to readily fertilize the eggs of S. droebachiensis. Our results demonstrate that asymmetrical gamete compatibilities in strongylocentrotids can give rise to asymmetrical patterns of introgression but suggest that gamete traits alone cannot be responsible for maintaining species integrities. The genetic boundaries between strongylocentrotid urchin species in the northeast Pacific appear to be related to postzygotic isolating mechanisms that scale with divergence times and not intrinsic gametic incompatibilities per se.

Introgression versus immigration in hybridizing high-dispersal echinoderms.

Phylogeographic studies designed to estimate rates and patterns of genetic differentiation within species often reveal unexpected and graphically striking cases of allele or haplotype sharing between species (introgression) via hybridization and backcrossing. Does introgression between species significantly influence population genetic structure relative to more conventional sources of differentiation (drift) and similarity (dispersal) among populations within species? Here we use mtDNA sequences from four species in two genera of sea urchins and sea stars to quantify the relative magnitude of gene flow across oceans and across species boundaries in the context of the trans-Arctic interchange of marine organisms between the Pacific and Atlantic oceans. In spite of the much smaller distances between sympatric congeners, rates of gene flow between sympatric species via heterospecific gamete interactions were small and significantly lower than gene flow across oceans via dispersal of planktonic larvae. We conclude that, in these cases at least, larvae are more effective than gametes as vectors of gene flow.

Genotype by sequencing identifies natural selection as a driver of intraspecific divergence in Atlantic populations of the high dispersal marine invertebrate, Macoma petalum.

Mitochondrial DNA analyses indicate that the Bay of Fundy population of the intertidal tellinid bivalve Macoma petalum is genetically divergent from coastal populations in the Gulf of Maine and Nova Scotia. To further examine the evolutionary forces driving this genetic break, we performed double digest genotype by sequencing (GBS) to survey the nuclear genome for evidence of both neutral and selective processes shaping this pattern. The resulting reads were mapped to a partial transcriptome of its sister species, M. balthica, to identify single nucleotide polymorphisms (SNPs) in protein-coding genes. Population assignment tests, principle components analyses, analysis of molecular variance, and outlier tests all support differentiation between the Bay of Fundy genotype and the genotypes of the Gulf of Maine, Gulf of St. Lawrence, and Nova Scotia. Although both neutral and non-neutral patterns of genetic subdivision were significant, genetic structure among the regions was nearly 20 times higher for loci putatively under selection, suggesting a strong role for natural selection as a driver of genetic diversity in this species. Genetic differences were the greatest between the Bay of Fundy and all other population samples, and some outlier proteins were involved in immunity-related processes. Our results suggest that in combination with limited gene flow across the mouth of the Bay of Fundy, local adaptation is an important driver of intraspecific genetic variation in this marine species with high dispersal potential.

Life-history predicts past and present population connectivity in two sympatric sea stars.

Life-history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic-developing sea star, Parvulastra exigua, has lower levels of within- and among-population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic-developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic-developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post-Pleistocene introgression.

Colonization, dispersal, and hybridization influence phylogeography of North Atlantic sea urchins (Strongylocentrotus droebachiensis).

We used frequency-based and coalescent-based phylogeographic analysis of sea urchin (Strongylocentrotus droebachiensis) mitochondrial DNA (mtDNA) sequences and previously published microsatellite data to understand the relative influence of colonization and gene flow from older (north Pacific) and younger (northeast Atlantic) sea urchin populations on genetic variation in the northwest Atlantic. We found strong evidence of survival of northwestern Atlantic populations in local Pleistocene glacial refugia: most haplotypes were the same as or closely related to Pacific haplotypes, with deep gene genealogies that reflect divergence times within the northwestern Atlantic that are much older than the last glacial maximum. We detected gene flow across the North Atlantic in the form of haplotypes shared with or recently descended from European populations. We also found evidence of significant introgression of haplotypes from a closely related species (S. pallidus). The relative magnitude of gene flow estimated by coalescent methods (and the effective population size differences among oceanic regions) depended on the genetic marker used. In general, we found very small effective population size in the northeastern Atlantic and high trans-Arctic gene flow between the Pacific and northwestern Atlantic. Both analyses suggested significant back-migration to the Pacific. However, microsatellites more strongly reflected older Pacific migration (with similar effective population sizes across the Arctic), whereas mtDNA sequences appeared to be more sensitive to recent trans- Atlantic dispersal (with larger differences in effective population size). These differences across marker types might have several biological or methodological causes, and they suggest caution in interpretation of the results from a single locus or class of markers.

From puffins to plankton: a DNA-based analysis of a seabird food chain in the northern Gulf of Maine.

The predator-prey interactions within food chains are used to both characterize and understand ecosystems. Conventional methods of constructing food chains from visual identification of prey in predator diet can suffer from poor taxonomic resolution, misidentification, and bias against small or completely digestible prey. Next-generation sequencing (NGS) technology has become a powerful tool for diet reconstruction through barcoding of DNA in stomach content or fecal samples. Here we use multi-locus (16S and CO1) next-generation sequencing of DNA barcodes on the feces of Atlantic puffin (Fratercula arctica) chicks (n=65) and adults (n=64) and the stomach contents of their main prey, Atlantic herring (Clupea harengus, n=44) to investigate a previously studied food chain. We compared conventional and molecular-derived chick diet, tested the similarity between the diets of puffin adults and chicks, and determined whether herring prey can be detected in puffin diet samples. There was high variability in the coverage of prey groups between 16S and CO1 markers. We identified more unique prey with our 16S compared to CO1 barcoding markers (51 and 39 taxa respectively) with only 12 taxa identified by both genes. We found no significant difference between the 16S-identified diets of puffin adults (n=17) and chicks (n=41). Our molecular method is more taxonomically resolved and detected chick prey at higher frequencies than conventional field observations. Many likely planktonic prey of herring were detected in feces from puffin adults and chicks, highlighting the impact secondary consumption may have on the interpretation of molecular dietary analysis. This study represents the first simultaneous molecular investigation into the diet of multiple components of a food chain and highlights the utility of a multi-locus approach to diet reconstruction that is broadly applicable to food web analysis.

Next-generation phylogeography: a targeted approach for multilocus sequencing of non-model organisms.

The field of phylogeography has long since realized the need and utility of incorporating nuclear DNA (nDNA) sequences into analyses. However, the use of nDNA sequence data, at the population level, has been hindered by technical laboratory difficulty, sequencing costs, and problematic analytical methods dealing with genotypic sequence data, especially in non-model organisms. Here, we present a method utilizing the 454 GS-FLX Titanium pyrosequencing platform with the capacity to simultaneously sequence two species of sea star (Meridiastra calcar and Parvulastra exigua) at five different nDNA loci across 16 different populations of 20 individuals each per species. We compare results from 3 populations with traditional Sanger sequencing based methods, and demonstrate that this next-generation sequencing platform is more time and cost effective and more sensitive to rare variants than Sanger based sequencing. A crucial advantage is that the high coverage of clonally amplified sequences simplifies haplotype determination, even in highly polymorphic species. This targeted next-generation approach can greatly increase the use of nDNA sequence loci in phylogeographic and population genetic studies by mitigating many of the time, cost, and analytical issues associated with highly polymorphic, diploid sequence markers.

Discordant distribution of populations and genetic variation in a sea star with high dispersal potential.

Patiria miniata, a broadcast-spawning sea star species with high dispersal potential, has a geographic range in the intertidal zone of the northeast Pacific Ocean from Alaska to California that is characterized by a large range gap in Washington and Oregon. We analyzed spatial genetic variation across the P. miniata range using multilocus sequence data (mtDNA, nuclear introns) and multilocus genotype data (microsatellites). We found a strong phylogeographic break at Queen Charlotte Sound in British Columbia that was not in the location predicted by the geographical distribution of the populations. However, this population genetic discontinuity does correspond to previously described phylogeographic breaks in other species. Northern populations from Alaska and Haida Gwaii were strongly differentiated from all southern populations from Vancouver Island and California. Populations from Vancouver Island and California were undifferentiated with evidence of high gene flow or very recent separation across the range disjunction between them. The surprising and discordant spatial distribution of populations and alleles suggests that historical vicariance (possibly caused by glaciations) and contemporary dispersal barriers (possibly caused by oceanographic conditions) both shape population genetic structure in this species.

Range-wide genetic homogeneity in the California sea mussel (Mytilus californianus): a comparison of allozymes, nuclear DNA markers, and mitochondrial DNA sequences.

We tested for genetic differentiation among six populations of California sea mussels (Mytilus californianus) sampled across 4000 km of its geographical range by comparing patterns of variation at four independent types of genetic markers: allozymes, single-copy nuclear DNA markers, and DNA sequences from the male and female mitochondrial genomes. Despite our extensive sampling and genotyping efforts, we detected no significant differences among localities and no signal of isolation by distance suggesting that M. californianus is genetically homogeneous throughout its range. This concordance differs from similar studies on other mytilids, especially in the role of postsettlement selection generating differences between exposed coastal and estuarine habitats. To assess if this homogeneity was due to M. californianus not inhabiting estuarine environments, we reviewed studies comparing allozymes with other classes of nuclear DNA markers. Although both types of markers gave broadly consistent results, there was a bias favouring studies in which allozymes were more divergent than DNA markers (nine to three) and a disproportionate number of these cases involved marine taxa (seven). Furthermore, allozymes were significantly more heterogeneous than DNA markers in three of the four studies that sampled coastal and estuarine habitats. We conclude that the genetic uniformity exhibited by M. californianus may result from a combination of extensive gene flow and the lack of exposure to strong selective gradients across its range.