Phylogenetic analysis of cryptic speciation in the polychaete Pygospio elegans

Development in marine invertebrate species can take place through a variety of modes and larval forms, but within a species, developmental mode is typically uniform. Poecilogony refers to the presence of more than one mode of development within a single species. True poecilogony is rare, however, and in some cases, apparent poecilogony is actually the result of variation in development mode among recently diverged cryptic species. We used a phylogenetic approach to examine whether poecilogony in the marine polychaete worm, Pygospio elegans, is the result of cryptic speciation. Populations of worms identified as P. elegansooded, and intermediate larvae; these modes are found both within and among populations. We examined sequence variation among partial mitochondrial cytochrome c oxidase subunit I sequences obtained for 279 individual worms sampled across broad geographic and environmental scales. Despite a large number of unique haplotypes (121 haplotypes from 279 individuals), sequence divergence among European samples was low (1.7%) with most of the sequence variation observed within populations, relative to the variation among regions. More importantly, we observed common haplotypes that were widespread among the populations we sampled, and the two most common haplotypes were shared between populations differing in developmental mode. Thus, our results support an earlier conclusion of poecilogony in P elegans. In addition, predominantly planktonic populations had a larger number of population-specific low-frequency haplotypes. This finding is largely consistent with interspecies comparisons showing high diversity for species with planktonic developmental modes in contrast to low diversity in species with brooded developmental modes.     

Evolution of poecilogony from planktotrophy: cryptic speciation, phylogeography, and larval development in the gastropod genus Alderia

Poecilogony, a rare phenomenon in marine invertebrates, occurs when alternative larval morphs differing in dispersal potential or trophic mode are produced from a single genome. Because both poecilogony and cryptic species are prevalent among sea slugs in the suborder Sacoglossa (Gastropoda: Opisthobranchia), molecular data are needed to confirm cases of variable development and to place them in a phylogenetic context. The nominal species Alderia modesta produces long-lived, feeding larvae throughout the North Atlantic and Pacific, but in California can also produce short-lived larvae that metamorphose without feeding. We collected morphological, developmental, and molecular data for Alderia from 17 sites spanning the eastern and western Pacific and North Atlantic. Estuaries south of Bodega Harbor, California, contained a cryptic species (hereafter Alderia sp.) with variable development, sister to the strictly planktotrophic A. modesta. The smaller Alderia sp. seasonally toggled between planktotrophy and lecithotrophy, with some individuals differing in development but sharing mitochondrial DNA haplotypes. The sibling species overlapped in Tomales Bay, California, but showed no evidence of hybridization; laboratory mating trials suggest postzygotic isolation has arisen. Intra- and interspecific divergence times were estimated using a molecular clock calibrated with geminate sacoglossans. Speciation occurred about 4.1 million years ago during a major marine radiation in the eastern Pacific, when large inland embayments in California may have isolated ancestral populations. Atlantic and Pacific A. modesta diverged about 1.7 million years ago, suggesting trans-Arctic gene flow was interrupted by Pleistocene glaciation. Both Alderia species showed evidence of late Pleistocene population expansion, but the southern Alderia sp. likely experienced a more pronounced bottleneck. Reduced body size may have incurred selection against obligate planktotrophy in Alderia sp. by limiting fecundity in the face of high larval mortality rates in warm months. Alternatively, poecilogony may be an adaptive response to seasonal opening of estuaries, facilitating dispersal by long-lived larvae. An improved understanding of the forces controlling seasonal shifts in development in Alderia sp. may yield insight into the evolutionary forces promoting transitions to nonfeeding larvae.     

Seasonal genetic variation associated with population dynamics of a poecilogonous polychaete worm

Poecilogonous species show variation in developmental mode, with larvae that differ both morphologically and ecologically. The spionid polychaete Pygospio elegans shows variation in developmental mode not only between populations, but also seasonally within populations. We investigated the consequences of this developmental polymorphism on the spatial and seasonal genetic structure of P. elegans at four sites in the Danish Isefjord‐Roskilde‐Fjord estuary at six time points, from March 2014 until February 2015. We found genetic differentiation between our sampling sites as well as seasonal differentiation at two of the sites. The seasonal genetic shift correlated with the appearance of new size cohorts in the populations. Additionally, we found that the genetic composition of reproductive individuals did not always reflect the genetic composition of the entire sample, indicating that variance in reproductive success among individuals is a likely explanation for the patterns of chaotic genetic patchiness observed during this and previous studies. The heterogeneous, unpredictable character of the estuary might maintain poecilogony in P. elegans as a bet‐hedging strategy in the Isefjord‐Roskilde‐Fjord complex in comparison with other sites where P. elegans are expected to be fixed to a certain mode of development.     

Poecilogony and population genetic structure in Elysia pusilla (Heterobranchia: Sacoglossa), and reproductive data for five sacoglossans that express dimorphisms in larval development

Credible cases of poecilogony, the production of two distinct larval morphs within a species, are extremely rare in marine invertebrates, yet peculiarly common in a clade of herbivorous sea slugs, the Sacoglossa. Only five animal species have been reported to express dimorphic egg sizes that result in planktotrophic and lecithotrophic larvae: the spionid polychaete Streblospio benedicti and four sacoglossans distributed in temperate estuaries or the Caribbean. Here, we present developmental and genetic evidence for a fifth case of poecilogony via egg-size dimorphism in the Sacoglossa and the first example from the tropical Indo-Pacific. The sea slug Elysia pusilla produced both planktotrophic and lecithotrophic larvae in Guam and Japan. Levels of genetic divergence within populations were markedly low and rule out cryptic species. However, divergence among populations was exceptionally high (10-12% at the mitochondrial cytochrome c oxidase I locus), illustrating that extensive phylogeographic structure can persist in spite of the dispersal potential of planktotrophic larvae. We review reproductive, developmental, and ecological data for the five known cases of poecilogony in the Sacoglossa, including new data for Costasiella ocellifera from the Caribbean. We hypothesize that sacoglossans achieve lecithotrophy at smaller egg sizes than do related clades of marine heterobranchs, which may facilitate developmental plasticity that is otherwise vanishingly rare among animals. Insight into the environmental drivers and evolutionary results of shifts in larval type will continue to be gleaned from population-level studies of poecilogonous taxa like E. pusilla, and should inform life-history theory about the causes and consequences of alternative development modes in marine animals.     

Polymorphism in developmental mode and its effect on population genetic structure of a spionid polychaete, Pygospio elegans

Population genetic structure of sedentary marine species is expected to be shaped mainly by the dispersal ability of their larvae. Long-lived planktonic larvae can connect populations through migration and gene flow, whereas species with nondispersive benthic or direct-developing larvae are expected to have genetically differentiated populations. Poecilogonous species producing different larval types are ideal when studying the effect of developmental mode on population genetic structure and connectivity. In the spionid polychaete Pygospio elegans, different larval types have been observed between, and sometimes also within, populations. We used microsatellite markers to study population structure of European P. elegans from the Baltic Sea (BS) and North Sea (NS). We found that populations with planktonic larvae had higher genetic diversity than did populations with benthic larvae. However, this pattern may not be related to developmental mode, since in P. elegans, developmental mode may be associated with geography. Benthic larvae were more commonly seen in the brackish BS and planktonic larvae were predominant in the NS, although both larval types also are found from both areas. Significant isolation-by-distance (IBD) was found overall and within regions. Most of the pair-wise F(ST) comparisons among populations were significant, although some geographically close populations with planktonic larvae were found to be genetically similar. However, these results, together with the pattern of IBD, autocorrelation within populations, as well as high estimated local recruitment, suggest that dispersal is limited in populations with planktonic larvae as well as in those with benthic larvae. The decrease in salinity between the NS and BS causes a barrier to gene flow in many marine species. In P. elegans, low, but significant, differentiation was detected between the NS and BS (3.34% in AMOVA), but no clear transition zone was observed, indicating that larvae are not hampered by the change in salinity.     

Poecilogony in the caenogastropod Calyptraea lichen (Mollusca: Gastropoda)

In marine invertebrates, polymorphism and polyphenism in mode of development are known as “poecilogony.” Understanding the environmental correlates of poecilogony and the developmental mechanisms that produce it could contribute to a better understanding of evolutionary transitions in mode of development. However, poecilogony is rare in marine invertebrates, with only ten recognized, well‐documented cases. Five examples occur in sacoglossan gastropods, and five occur in spionid polychaetes. Here, we document the eleventh case, and the first in a caenogastropod mollusc. Females of Calyptraea lichen collected in the field or reared in the laboratory often produce broods of planktotrophic larvae. They can also be collected with mixed broods, in which each capsule contains planktotrophic larvae, nurse embryos, and adelphophagic embryos. Adelphophages eat the nurse embryos and hatch as short‐lived lecithotrophic larvae, or even as juveniles. Mitochondrial COI and 16S DNA sequences for females with different types of broods differ by less than 0.5%, supporting conspecific status. Some females collected in the field with mixed broods subsequently produced planktotrophic broods, demonstrating that females can produce two different kinds of broods. Calyptraea lichen is therefore polyphenic in two ways: mode of development can vary among embryos within a capsule, and females can change the types of broods they produce.     

Relationships among development, ecology, and morphology in the evolution of Echinoderm larvae and life cycles

The evolution of life cycles involves transitions between discrete states in one or more of the characters that comprise a developmental pattern. In this paper, we examine three of the major life cycle characters and the states for these characters. Using examples from echinoderms, we discuss the evolutionary transitions that have occurred in the type of morphogenesis, developmental habitat, and mode of nutrition during development. We evaluate the functional requirements associated with these transitions to infer the likelihood (frequency or rapidity) of change in a given character and of biases in the polarity of character state transitions. Using comparisons of closely related species, we evaluate the change between states in one character for dependence on the state of, or correlated changes in, other characters. Based on our analysis of congeneric species that differ in developmental habitat, we conclude that the transition between pelagic and benthic development is an ecological change that is independent of changes in morphogenesis and should be reversible. In contrast, the transition from feeding to nonfeeding development has been considered to be irreversible because it involves marked changes in larval morphology. We re-examine the transition between different modes of larval nutrition in light of recent studies that show that there exists a continuum of nutritional strategies between planktotrophy and lecithotrophy. This continuum is largely determined by variation in maternal investment and does not involve alterations in larval morphology. We suggest that the boundary between planktotrophy and lecithotrophy is frequently crossed and that this transition is reversible. Ecological changes represent the crossing of a functional threshold. Only after crossing the threshold, do larvae experience qualitatively different selective pressures that can lead to subsequent changes in morphology and development. Two different changes have occurred in the type of morphogenesis: the simplification of larval morphology that is associated with obligate (nonfeeding) lecithotrophy and the loss of the larval body plan in the evolution from indirect to direct development. It is the modification of morphology independent of the ecological changes that requires alterations in developmental processes, constrains evolutionary options, imposes irreversibility, and establishes the discrete nature of larval patterns in marine invertebrates.     

Evolutionary changes in the timing of gut morphogenesis in larvae of the marine annelid Streblospio benedicti

The planktonic larvae of marine invertebrates are diverse in their nutritional modes, suggesting that evolutionary transitions in larval nutritional mode have been frequent. One approach to identifying the developmental changes that play important roles in such transitions is to compare "intermediate" larval forms to closely related larvae representative of their common ancestor. Here we make such a comparison between obligately planktotrophic and facultatively feeding larvae of the poecilogonous polychaete annelid Streblospio benedicti. We used feeding experiments to show that the derived, facultatively feeding larvae of this species develop the ability to feed at a later developmental stage (five muscle bands) than planktotrophic larvae (two to three muscle bands). This delay in the onset of feeding ability does not appear to be caused by delay in the formation of particle capture structures, but instead by delay in the development of a continuous, functional gut. These observations are consistent with the hypothesis that evolutionary increases in egg size in annelids lead predictably to heterochronic delays in gut development, and hence to transitions in larval nutritional mode.     

Larval development of Phoronis pallida (Phoronida): implications for morphological convergence and divergence among larval body plans

Morphological variation among larval body plans must be placed into a phylogenetic and ecological context to assess whether similar morphologies are the result of phylogenetic constraints or convergent selective pressures. Investigations are needed of the diverse larval forms within the Lophotrochozoa, especially the larvae of phoronids and brachiopods. The actinotroch larva of Phoronis pallida (Phoronida) was reared in the laboratory to metamorphic competence. Larval development and growth were followed with video microscopy, SEM, and confocal microscopy. Early developmental features were similar to other phoronid species. Gastrulation was accomplished by embolic invagination of the vegetal hemisphere. Mesenchymal cells were found in the remaining blastocoelic space after invagination began. Mesenchymal cells formed the body wall musculature during the differentiation of larval features. Body wall musculature served as the framework from which all other larval muscles proliferated. Larval growth correlated best with developmental stage rather than age. Consistent with other phoronid species, differentiation of juvenile tissues occurred most rapidly at the latest stages of larval development. The minimum precompetency period of P. pallida was estimated to be approximately 4-6 weeks. Previously published studies have documented that the planktonic embryos of P. pallida develop faster than the brooded embryos of P. vancouverensis. However, these data showed that the difference in developmental rate between the two species decreased in succeeding larval stages. There may be convergent selective pressures that result in similar timing to metamorphic competence among phoronid and brachiopod planktotrophic larval types. Morphological differences between these larval types result from heterochronic developmental shifts in the differentiation of juvenile tissue. Similarities in the larval morphology of phoronids and basal deuterostomes are likely the result of functional and developmental constraints rather than a shared (recent) evolutionary origin. These constraints are imposed by the functional design of embryological stages, feeding structures, and swimming structures.     

Seasonal polyphenism in larval type: rearing environment influences the development mode expressed by adults in the sea slug Alderia willowi

Dimorphisms occur when alternative developmental pathways produce discrete phenotypes within a species, and may promote evolutionary novelty in morphology, life history, and behavior. Among marine invertebrates, intra-specific dimorphism in larval type (poecilogony) is notably rare, but should provide insight into the selective forces acting on larval strategies. Most established cases of poecilogony appear to be allelic polymorphisms, with local expression regulated by population-genetic processes. Here, we present evidence that dimorphic larval development in the sea slug Alderia willowi is a seasonal polyphenism; the type of larvae produced by an adult slug depends on the rearing environment in which that slug matured. In field surveys of 1996-1999 and 2007-2009, the population in Mission Bay, San Diego (California, USA) produced only short-lived lecithotrophic larvae in summer and early fall, but a varying proportion of slugs expressed planktotrophy in winter and spring. In laboratory experiments, slugs reared under summer conditions (high temperature, high salinity) produced the highest proportion of lecithotrophic offspring, whereas winter conditions (low temperature, low salinity) induced the lowest proportion of lecithotrophy. The shift to a nondispersive morph under summer conditions may be an adaptive response to historical closure of coastal wetlands during the dry season in southern California, which would inhibit dispersal by larvae of back-bay taxa. In most animal polyphenisms, a single larval type is produced and the rearing environment determines which adult phenotype develops. In contrast, alternative larval morphs are produced by A. willowi in response to seasonal cues experienced by the adult stage, varying the phenotype and dispersal potential of offspring. As the only known case of polyphenism in mode of larval development, A. willowi should become a model organism for mechanistic studies of dimorphism and the evolution of alternative life histories.     

Patterns of nuclear genetic variation in the poecilogonous polychaete Streblospio benedicti

The evolution of marine larvae is replete with transitions in trophic mode, but little is known about how, at the genetic level, these transitions are achieved. Basic parameters, including the number of underlying loci, their molecular characteristics, and the population-genetic processes that drive transitions remain unknown. Streblospio benedicti, an abundant benthic polychaete with heritable poecilogony, provides a unique genetically tractable system for addressing these issues. Individuals of S. benedicti vary in diverse aspects of development. Some females produce small, planktotrophic larvae, whereas others produce large, yolky larvae capable of settling without feeding. Here, I present estimates of basic features of nuclear genetic variation in S. benedicti to lay the foundations for subsequent efforts to understand the genetic basis of poecilogony. Sequence of ～20?kb of random nuclear DNA indicates that the nucleotide composition, at 62.1% A?+?T, is typical of lophotrochozoan genomes. Population-genetic data, acquired by sequencing two loci (～2500 bp) in multiple animals of each developmental morph, indicate that the morphs exhibit very little differentiation at random loci. Nucleotide heterozygosity (θπ) is ～0.5-1% per site, and linkage disequilibrium decays within a few kilobases (ρ ～?3?×?10(-3) per site). These data suggest that genetic mapping by association will require a high density of markers, but linkage mapping and identification of regions of elevated inter-morph differentiation hold great promise.     

Morphogenesis and phenotypic divergence in two developmental morphs of Streblospio benedicti (Annelida,Spionidae)

Abstract. The morphology of marine invertebrate larvae is strongly correlated with egg size and larval feeding mode. Planktotrophic larvae typically have suites of morphological traits that support a planktonic, feeding life style, while lecithotrophic larvae often have larger, yolkier bodies, and in some cases, a reduced expression of larval traits. Poecilogonous species provide interesting cases for the analysis of early morphogenesis, as two morphs of larvae are produced by a single species. We compared morphogenesis in planktotrophic and lecithotrophic morphs of the poecilogonous annelid Streblospio benedicti from the trochophore stage through metamorphosis, using observations of individuals that were observed alive, with scanning electron microscopy, or in serial sections. Offspring of alternate developmental morphs of this species are well known to have divergent morphologies in terms of size, yolk content, and the presence of larval bristles. We found that some phenotypic differences between morphs occur as traits that are present in only one morph (e.g., larval bristles, bacillary cells on the prostomium and pygidium), but that much of the phenotypic divergence is based on heterochronic changes in the differentiation of shared traits (e.g., gut and coelom). Tissue and organ development are compared in both morphs in terms of their structure and ontogenetic change throughout early development and metamorphosis.     

It's about time: divergence, demography, and the evolution of developmental modes in marine invertebrates

Differences in larval developmental mode are predicted to affect ecological and evolutionary processes ranging from gene flow and population bottlenecks to rates of population recovery from anthropogenic disturbance and capacity for local adaptation. The most powerful tests of these predictions use comparisons among species to ask how phylogeographic patterns are correlated with the evolution and loss of prolonged planktonic larval development. An important and largely untested assumption of these studies is that interspecific differences in population genetic structure are mainly caused by differences in dispersal and gene flow (rather than by differences in divergence times among populations or changes in effective population sizes), and that species with similar patterns of spatial genetic variation have similar underlying temporal demographic histories. Teasing apart these temporal and spatial patterns is important for understanding the causes and consequences of evolutionary changes in larval developmental mode. New analytical methods that use the coalescent history of allelic diversity can reveal these temporal patterns, test the strength of traditional population-genetic explanations for variation in spatial structure based on differences in dispersal, and identify strongly supported alternative explanations for spatial structure based on demographic history rather than on gene flow alone. We briefly review some of these recent analytical developments, and show their potential for refining ideas about the correspondence between the evolution of larval developmental mode, population demographic history, and spatial genetic variation.     

The effects of nurse eggs and sibling interactions on the larval development of the poecilogonous annelid Boccardia proboscidea (Spionidae)

In poecilogony, different types of larvae are produced within the same species. Previous studies have suggested maternal control of the production of larval types; however, it is not clear which factors or mechanisms generate contrasting developmental patterns among siblings. The spionid polychaete Boccardia proboscidea produces within the same capsule adelphophagic larvae that eat nurse eggs and siblings and complete all or most of their development inside the capsule (Type A larvae), and larvae with little growth until they hatch as planktotrophic larvae (Type B larvae). In this study, we manipulated capsule content to explore the factors determining larval type in B. proboscidea and the role of extra‐embryonic maternal nutrition and sib–sib interaction in the developmental fate of offspring. When early larval stages were grown individually in vitro, with nurse eggs as the only food source, some of them remained small, while others continue developing into larger pre‐competent larvae by feeding on nurse eggs. This suggests that larval types in B. proboscidea are determined very early in development and are not solely the product of sib–sib interaction inside the capsule. However, our data also suggest that hatching size variability within larval types of a clutch depends on nurse egg availability. Type B larvae grew normally to metamorphosis when phytoplankton was available, but suffered high rates of cannibalism by Type A larvae. These results are consistent with the hypothesis that individual larval fates are determined very early in development and that once their fate is determined, hatching size and intracapsular survival are affected by maternal food provisioning and sibling interaction.     

Testing adaptive hypotheses on the evolution of larval life history in acorn and stalked barnacles

Despite strong selective pressure to optimize larval life history in marine environments, there is a wide diversity with regard to developmental mode, size, and time larvae spend in the plankton. In the present study, we assessed if adaptive hypotheses explain the distribution of the larval life history of thoracican barnacles within a strict phylogenetic framework. We collected environmental and larval trait data for 170 species from the literature, and utilized a complete thoracican synthesis tree to account for phylogenetic nonindependence. In accordance with Thorson's rule, the fraction of species with planktonic‐feeding larvae declined with water depth and increased with water temperature, while the fraction of brooding species exhibited the reverse pattern. Species with planktonic‐nonfeeding larvae were overall rare, following no apparent trend. In agreement with the “size advantage” hypothesis proposed by Strathmann in 1977, egg and larval size were closely correlated. Settlement‐competent cypris larvae were larger in cold water, indicative of advantages for large juveniles when growth is slowed. Planktonic larval duration, on the other hand, was uncorrelated to environmental variables. We conclude that different selective pressures appear to shape the evolution of larval life history in barnacles.     

Perspectives on cardiac physiological ontogeny in crustaceans

Crustacean embryonic and larval systems offer a unique and valuable tool for furthering our understanding of both developmental processes and physiological regulatory mechanisms. The diverse array of developmental patterns exhibited by crustaceans allows species choice to be based on the specific questions being investigated, where defined larval forms are chosen based on their developmental pattern, degree of maturation or regulatory capabilities. However, this great diversity in developmental patterns, as well as crustacean diversity, can also confound ones ability to define or identify species for investigation. These issues are addressed and suggestions put forth to clarify some of the problems. The complexity and overlapping nature of adult cardio-regulatory systems makes teasing them apart difficult. Embryonic and larval systems exhibit varying degrees of regulatory complexity depending on developmental stage and ontogenetic pattern. This can allow complex adult regulatory systems to be teased apart temporally, as the developing animal builds regulatory pathways. Equally important is the nature of crustacean larvae; many undergo dramatic metamorphoses in cases where the larvae have adaptations to environments different to those of the adult. During environmental transitions physiological adaptations to immediate change should take precedence over long-term adult adaptations. It is therefore possible to look at physiological responses as a function of developmental/environmental adaptation, independent of adult functions.     

Phenotypic and genetic consequences of size selection at the larval stage in the Pacific oyster (Crassostrea gigas)

The life histories of oysters in the genus Crassostrea, like those of most marine bivalves, are typified by high fecundity and low survival in nature. Rearing conditions in hatcheries however ensure optimized density, diet, and temperature. Hatcheries are becoming increasingly important for the production of juveniles in aquaculture, and their culture practices often include culling of slow growing larvae to reduce and synchronize the time taken to reach settlement. Because previous studies have found substantial genetic variation for early life developmental traits in Crassostrea gigas, these culling practices are likely to cause highly different selective pressures in hatcheries from those in the natural environment. We studied the phenotypic and genetic impact of such culling practices in a factorial cross between 10 males and 3 females subjected to progressive culling of the smallest 50% of larvae, compared with a non-culled control. Measurements were made on larval growth, survival, time taken to attain pediveliger stage and settlement success. Culling had a larger effect on the variance of these larval traits than on their means. The larvae in culled cultures were approximately 10% larger than those in controls, whereas the coefficient of variation was reduced by 30-40%. Culling also reduced the mean time to settlement by 12% and its variance by 55%. Using a multiplexed set of microsatellite markers to trace parentage, we also estimated the variance in reproductive success in a controlled experiment to quantify the consequences of intensive hatchery rearing practices. We also focused on changes in effective population size and genetic structure over time (and developmental stages). Our results show a loss of genetic diversity following removal of the smallest larvae by culling, as well as temporally varying genetic structure of the larval population. This supports the existence of genetic variability in early life developmental traits in C. gigas. Culling in hatcheries, like size-related selective pressures in the wild, are likely to have a significant genetic impact, through their effects on the timing of settlement.