Life history and population genetic structure of sea stars from the family Asterinidae

Life history can influence population genetic variation by altering patterns of gamete union and dispersal. Sea stars from the family Asterinidae have evolved similar life histories multiple times in parallel including planktonic feeding larvae, planktonic non-feeding larvae, development in benthic egg masses, and viviparity. In this thesis I first examine the population genetic structure of a widespread planktotrophic asterinid sea star from the East Pacific (Patiria miniata). I use mitochondrial sequence markers to determine whether extrinsic factors such as vicariance or intrinsic properties such as dispersal mode are driving patterns of population genetic variation in this species. I then examine patterns of population genetic variation among eight additional asterinid species from Australia using a mixed species pool of genomic microsatellite markers. I use these microsatellite markers to characterize the genetic variation within groups of brooded offspring associated with the unusual life histories of two live bearing asterinids from the genus Parvulastra. Lastly, I examine the evolution of life history among the asterinids and use the phylogenetic relationships among species to examine the correlation between life history and population genetic structure in this group. Ultimately, I find that the degree to which intrinsic life history properties of asterinids and extrinsic factors contribute to population genetic variation varies among species and among clades. In P. miniata patterns of population genetic variation are influenced by both intrinsic and extrinsic factors. Using microsatellite markers I find that in general between-population genetic variation is high in benthic species (benthic egg laying and live bearing) relative to species with planktonic larvae and that genetic variation within populations is lower in benthic species relative to planktonic species. Lastly, I find that the degree to which phylogeny constrains the coevolution of population genetic structure and life history varies among life history characters and among-population genetic parameters. This thesis suggests that in many cases variation in life histories among a closely related group of marine species can predict patterns of population genetic variation. However, extrinsic factors can in some cases, act with or override life history characteristics in driving patterns of population genetic variation

Morphological and genetic variation indicate cryptic species within Lamarck’s little sea star, Parvulastra (=Patiriella) exigua

The asterinid sea star Parvulastra exigua (Lamarck) is a common member of temperate intertidal marine communities from geographically widespread sites around the southern hemisphere. Individuals from Australian populations lay benthic egg masses (through orally directed gonopores) from which nonplanktonic offspring hatch and metamorphose without a dispersing planktonic larval phase. Scattered reports in the taxonomic literature refer to a similar form in southern Africa with aborally directed gonopores (and possibly broadcast spawning of planktonic eggs and larvae); such differences would be consistent with cryptic species variation. Surveys of morphology and mtDNA sequences have revealed cryptic species diversity in other asterinid genera. Here we summarize the taxonomic history of Lamarck’s "Astérie exiguë" and survey morphological variation (the location of the gonopores) for evidence that some P. exigua populations include cryptic species with a different mode of reproduction. We found strong evidence for multiple species in the form of two phenotypes and modes of reproduction (oral and aboral gonopore locations) in populations from southern Africa and islands in the Atlantic and Indian oceans. Both modes of reproduction have broad geographic ranges. These results are consistent with previously published genetic data that indicate multiple species in African and island (but not Australian) populations

Dunlin_HighSeq_2014_R1_001.fastq

Genomic data from Illumina sequencing for Dunli

Dunlin_HighSeq_2014_R2_001.fastq

Genomic data from Illumina sequencing for Dunli

Genotyping-by-sequencing reveals genomic homogeneity among overwintering Pacific Dunlin (Calidris alpina pacifica) aggregations along the Pacific coast of North America

Information on how migratory populations are genetically structured during the overwintering season of the annual cycle can improve our understanding of the strength of migratory connectivity and help identify populations as units for management. Here, we use a genotype-by-sequencing approach to investigate whether population genetic structure exists among overwintering aggregations of the Pacific Dunlin subspecies (Calidris alpina pacifica) sampled at 2 spatial scales (within and among overwintering sites) in the eastern Pacific Flyway. Genome-wide analyses of 874 single nucleotide polymorphisms across 80 sampled individuals revealed no evidence for genetic differentiation among aggregations overwintering at 3 locations within the Fraser River Estuary (FRE) of British Columbia. Similarly, comparisons of aggregations in the FRE and those overwintering in southern sites in California and Mexico indicated no genetic segregation between northern and southern overwintering areas. These results suggest that Pacific Dunlin within the FRE, Sacramento Valley (California), and Guerrero Negro (Mexico) are genetically homogeneous, with no evident genetic structure between sampled sites or regions across the overwintering range. Despite no evidence for differentiation among aggregations, we identified a significant effect of geographical distance between sites on the distribution of individual genotypes in a redundancy analysis. A small proportion of the total genotypic variance (R 0.036, P = 0.011) was explained by the combined effect of latitude and longitude, suggesting weak genomic patterns of isolation-by-distance that are consistent with chain-like migratory connectivity between breeding and overwintering areas. Our study represents the first genome-scale investigation of population structure for a Dunlin subspecies and provides essential baseline estimates of genomic diversity and differentiation within the Pacific Dunlin

Data from: Genotyping-by-sequencing reveals genomic homogeneity among overwintering Pacific Dunlin (Calidris alpina pacifica) aggregations along the Pacific coast of North America

Information on how migratory populations are genetically structured during the overwintering season of the annual cycle can improve our understanding of the strength of migratory connectivity and help identify populations as units for management. Here, we use a genotype-by-sequencing approach to investigate whether population genetic structure exists among overwintering aggregations of the Pacific Dunlin subspecies (Calidris alpina pacifica) sampled at two spatial scales (i.e. within and among overwintering sites) in the eastern Pacific Flyway. Genome-wide analyses of 874 single nucleotide polymorphisms across 80 sampled individuals revealed no evidence for genetic differentiation among aggregations overwintering at three locations within the Fraser River Estuary (FRE) of British Columbia. Similarly, comparisons of aggregations in the FRE and those overwintering in southern sites in California and Mexico indicated no genetic segregation between northern and southern overwintering areas. These results suggest that Pacific Dunlin residing within the FRE, Sacramento Valley (California) and Guerrero Negro (Mexico) are genetically homogeneous, with no evident genetic structure between sampled sites or regions across the overwintering range. Despite no evidence for differentiation among aggregations, we identified a significant effect of geographical distance between sites on the distribution of individual genotypes in a redundancy analysis; however, a small proportion of the total genotypic variance (R2 = 0.036, P = 0.011) was explained by the combined effect of latitude and longitude, suggesting weak genomic patterns of isolation-by-distance that are consistent with chain-like migratory connectivity between breeding and overwintering areas. Our study represents the first genome-scale investigation of population structure for a Dunlin subspecies and provides essential baseline estimates of genomic diversity and differentiation within the Pacific Dunlin

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

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

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

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