Search

Blue mussels ( Mytilus edulis ) are ecologically and economically important marine invertebrates whose populations are at risk from climate change-associated variation in their environment, such as decreased coastal salinity. Blue mussels are osmoconfomers and use components of the metabolome (free amino acids) to help maintain osmotic balance and cellular function during low salinity exposure. However, little is known about the capacity of blue mussels during the planktonic larval stages to regulate metabolites during osmotic stress. Metabolite studies in species such as blue mussels can help improve our understanding of the species' physiology, as well as their capacity to respond to environmental stress. We used 1D ¹H nuclear magnetic resonance (NMR) and 2D total correlation spectroscopy (TOCSY) experiments to describe baseline metabolite pools in larval (veliger and pediveliger stages) and juvenile blue mussels (gill, mantle, and adductor tissues) under ambient conditions and to quantify changes in the abundance of common osmolytes in these stages during low salinity exposure. We found evidence for stage- and tissue-specific differences in the baseline metabolic profiles of blue mussels, which reflect variation in the function and morphology of each larval stage or tissue type of juveniles. These differences impacted the utilization of osmolytes during low salinity exposure, likely stemming from innate physiological variation. This study highlights the importance of foundational metabolomic studies that include multiple tissue types and developmental stages to adequately evaluate organismal responses to stress and better place these findings in a broader physiological context., Competing Interests: The authors declare no conflict of interest.

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.

The North Atlantic intertidal community provides a rich set of organismal and environmental material for the study of ecological genetics. Clearly defined environmental gradients exist at multiple spatial scales: there are broad latitudinal trends in temperature, meso-scale changes in salinity along estuaries, and smaller scale gradients in desiccation and temperature spanning the intertidal range. The geology and geography of the American and European coasts provide natural replication of these gradients, allowing for population genetic analyses of parallel adaptation to environmental stress and heterogeneity. Statistical methods have been developed that provide genomic neutrality tests of population differentiation and aid in the process of candidate gene identification. In this paper, we review studies of marine organisms that illustrate associations between an environmental gradient and specific genetic markers. Such highly differentiated markers become candidate genes for adaptation to the environmental factors in question, but the functional significance of genetic variants must be comprehensively evaluated. We present a set of predictions about locus-specific selection across latitudinal, estuarine, and intertidal gradients that are likely to exist in the North Atlantic. We further present new data and analyses that support and contradict these simple selection models. Some taxa show pronounced clinal variation at certain loci against a background of mild clinal variation at many loci. These cases illustrate the procedures necessary for distinguishing selection driven by internal genomic vs. external environmental factors. We suggest that the North Atlantic intertidal community provides a model system for identifying genes that matter in ecology due to the clarity of the environmental stresses and an extensive experimental literature on ecological function. While these organisms are typically poor genetic and genomic models, advances in comparative genomics have provided access to molecular tools that can now be applied to taxa with well-defined ecologies. As many of the organisms we discuss have tight physiological limits driven by climatic factors, this synthesis of molecular population genetics with marine ecology could provide a sensitive means of assessing evolutionary responses to climate change.

Chemoreception in marine invertebrates mediates a variety of ecologically important behaviors including defense, reproduction, larval settlement and recruitment, and feeding. The sensory pathways that regulate deposit-feeding activity by polychaetes living in sedimentary habitats are of particular interest because such feeding has profound effects on the physical and chemical properties of the habitat. Nevertheless, little is known concerning the molecular mechanisms of chemical signal transduction associated with deposit feeding and other behaviors in polychaetes. Chemosensory-based feeding behaviors are typically regulated by G-protein-coupled signal transduction pathways. However, the presence and role of such pathways have not been demonstrated in marine polychaetes. Methodologies involving degenerate primer-based reverse transcription with the polymerase chain reaction and rapid amplification of cDNA ends were used to identify and characterize a Galphaq subunit expressed in the feeding palps of the spionid polychaete Dipolydora quadrilobata. The D. quadrilobata Galphaq protein had high sequence similarity with previously reported Galphaq subunits from both invertebrate and vertebrate taxa. Immunhistochemistry and immunocytochemistry were used with confocal laser scanning microscopy and transmission electron microscopy to visualize the distribution of a Galphaq antibody in whole worms and in cilia of the feeding palps. Galphaq immunoreactivity was concentrated in the nuchal organs, food-groove cilia, and lateral/abfrontal cilia of the feeding palps. Because these structures are known to be involved in chemoreception, we propose that Galphaq isolated from D. quadrilobata is a key component of chemosensory signal transduction pathways in this species.

Gene duplication has been hypothesized to play a major role in the evolution of genes and genomes and in generating phenotypic diversity among the proteins those genes encode. We have identified duplicate genes for the glycolytic enzyme mannose-6-phosphate isomerase (MPI: EC 5.3.1.8) in marine mussels in the genus Mytilus. Overall, there was only 52% sequence identity (72% similarity) between the proteins encoded by these two genes, designated as Mpi-A and Mpi-B. Based on a comparison of the rate of non-synonymous substitution between orthologous and paralogous Mpi coding sequences obtained from Mytilus edulis and the congener M. trossulus we estimate that the duplication of Mpi in mussels occurred ~170 MYA. We detected paralog-specific differences in the ratio of non-synonymous to synonymous substitutions (k(a)/k(s)) and in the predicted net charge for MPI-A and MPI-B. Using a real-time quantitative RT-PCR assay we observed substantial changes in Mpi-A and Mpi-B transcript levels between tissue types; the strongest expression of Mpi-A was observed in mantle tissue while Mpi-B expression exceeded that of Mpi-A in gill and hepatopancreas tissues. Taken together, these observations suggest that different functional roles have evolved for these two Mytilus Mpi genes subsequent to gene duplication.

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.

Freeman and Byers (Reports, 11 August 2006, p. 831) presented evidence for the rapid evolution of antipredator defenses in the mussel Mytilus edulis. However, their analysis is confounded by three issues. Samples from some sites are likely to have included a second species, M. trossulus; their manipulation of chemical cues does not preclude other interpretations; and they failed to establish an adaptive significance to shell thickening.

The cytochrome c oxidase subunit 2 gene (COII) encodes a highly conserved protein that is directly responsible for the initial transfer of electrons from cytochrome c to cytochrome c oxidase (COX) crucial to the production of ATP during cellular respiration. Despite its integral role in electron transport, we have observed extensive intraspecific nucleotide and amino acid variation among 26 full-length COII sequences sampled from seven populations of the marine copepod, Tigriopus californicus. Although intrapopulation divergence was virtually nonexistent, interpopulation divergence at the COII locus was nearly 20% at the nucleotide level, including 38 nonsynonymous substitutions. Given the high degree of interaction between the cytochrome c oxidase subunit 2 protein (COX2) and the nuclear-encoded subunits of COX and cytochrome c (CYC), we hypothesized that some codons in the COII gene are likely to be under positive selection in order to compensate for amino acid substitutions in other subunits. Estimates of the ratio of nonsynonymous to synonymous substitution (omega), obtained using a series of maximum likelihood models of codon substitution, indicated that the majority of codons in T. californicus COII are under strong purifying selection (omega << 1), while approximately 4% of the sites in this gene appear to evolve under relaxed selective constraint (omega = 1). A branch-site maximum likelihood model identified three sites that may have experienced positive selection within the central California sequence clade in our COII phylogeny; these results are consistent with previous studies showing functional and fitness consequences among interpopulation hybrids between central and northern California populations.

Doubly uniparental inheritance of mtDNA (DUI) is commonly observed in several genera of bivalves. Under DUI, female offspring inherit mtDNA from their mothers, while male offspring inherit mtDNA from both parents but preferentially transmit the paternally inherited mtDNA to their sons. Several studies have shown that the female- and male-specific mtDNA lineages in blue mussels, Mytilus spp., vary by upward of 20% at the nucleotide level. In addition to high levels of nucleotide substitution, the present study observed substantial gender-based length polymorphism in the presumptive mitochondrial control region (=large unassigned region; LUR) of North American M. trossulus. In this species, female lineage LUR haplotypes are over 2 kb larger than male lineage LUR haplotypes. Analysis of sequence data for these length variants indicates that the F LUR haplotypes of North American M. trossulus contain sequences similar to the F lineage control region in the congeners M. edulis and M. galloprovincialis. Relative to the F LUR in the latter two species, however, the F lineage LUR haplotypes in M. trossulus contain two large sequence insertions, each nearly 1 kb in size. One of these insertions has high sequence similarity to the male lineage LUR of M. trossulus. The tandem arrangement of F and M control region sequences in the F lineage LUR of M. trossulus is most likely the result of nonhomologous recombination between the male and the female mitochondrial genomes in M. trossulus, a finding that has important implications regarding the transmission and evolution of blue mussel mitochondrial genomes.

The barnacle, Chelonibia testudinaria, is a common inhabitant of the marine turtle epibiont community and plays a key role in the development of this community. Phylogeographic analysis of 79 cytochrome c oxidase I (COX1) sequences for barnacles collected from five populations found contrasting patterns of divergence for populations in the Atlantic vs. the Pacific Ocean. Our analysis indicates that the two Pacific populations, Senri Beach, Japan and Bahia Magdalena, Mexico, were not only highly divergent from the Atlantic populations but are highly divergent from one another. We suggest that barnacles from these populations may represent cryptic species. In contrast, sequence divergence was greatly reduced among barnacles collected from Wassaw Island, GA, USA, Keewaydin, FL, USA, and Kyparissia, Pèloponnésus Island, Greece. A reduction in sequence diversity at the latter site was attributed to a recent range expansion into the Mediterranean Sea. We examined historical patterns of migration among the three Atlantic and Mediterranean populations using the program migrate. This analysis indicates a high rate of migration from Keewaydin to Wassaw Island, contrasted with a much lower rate of migration in the opposite direction. The estimated migration rate from Kyparissia to Keewaydin was also non-negligible. We suggest that the association between C. testudinaria and loggerhead turtles and the patterns of turtle migration have played key roles in the expansion of the range of C. testudinaria into the Mediterranean Sea and the subsequent patterns of barnacle migration. We further propose that the difference between ocean basins, with respect to the impact of host migration on barnacle gene flow, probably stems from the fact that host-mediated dispersal in the Atlantic depends on advanced stage juveniles and adults while any host-mediated dispersal in the Pacific would have to involve early "pelagic" stage juvenile loggerheads.

Geographically isolated populations may accumulate alleles that function well on their own genetic backgrounds but poorly on the genetic backgrounds of other populations. Consequently, interpopulation hybridization may produce offspring of low fitness as a result of incompatibilities arising in allopatry. Genes participating in these epistatic incompatibility systems remain largely unknown. In fact, despite the widely recognized importance of epistatic interactions among gene products, few data directly address the functional consequences of such interactions among natural genetic variants. In the marine copepod, Tigriopus californicus, we found that the cytochrome c variants isolated from two different populations each had significantly higher activity with the cytochrome c oxidase derived from their respective source population. Three amino acid substitutions in the cytochrome c protein appear to be sufficient to confer population specificity. These results suggest that electron transport system (ETS) proteins form coadapted sets of alleles within populations and that disruption of the coadapted ETS gene complex leads to functional incompatibilities that may lower hybrid fitness.

Mitochondrial energy production requires complex interactions among proteins encoded in both the nuclear and mitochondrial genomes. The intergenomic coevolution of interacting gene products has been previously suggested based on interspecific comparisons of cytochrome c (encoded by the nuclear CYC gene) and cytochrome c oxidase (partly encoded in the mitochondrial DNA by the COX1, COX2 and COX3 genes). In the intertidal copepod, Tigriopus californicus, non-synonymous substitutions in the COX1 gene have previously been found in interpopulation comparisons. In order to determine if CYC also shows interpopulation variation, this gene was isolated from a cDNA library using a degenerate primer/polymerase chain reaction approach. Characterization of a cDNA sequence and 25 genomic DNA sequences derived from four T. californicus populations yielded the following results: (1) the T. californicus CYC gene is interrupted by an intron that occurs at the same position as the intron found in vertebrate CYC genes; (2) there is extensive sequence variation within both the coding region and intron of this gene and the vast majority of this variation occurs between sequences drawn from geographically distinct populations; (3) the coding sequence variation includes a minimum of five amino acid replacement substitutions; (4) segregation of length variants among offspring in an interpopulation cross revealed genotypic ratios consistent with the proposed allelic nature of the CYC variants. These results demonstrate that the requisite genetic variation required for intergenomic coevolution exists in the CYC-COX system in T. californicus.

-Mytilus edulis and M. galloprovincialis are two blue mussel species that coexist in western Europe. Previously, we reported that M. galloprovincialis populations contain female and male haplotypes that are fixed in M. edulis populations as well as unique haplotypes. This study assesses whether paraphyly for these species is due to introgression or incomplete lineage extinction. The lineage extinction hypothesis predicts that the shared mtDNA haplotypes in M. galloprovincialis will be significantly diverged from those in M. edulis and form distinct sequence clades. In contrast, the introgression hypothesis proposes that M. edulis haplotypes have only recently been introduced into M. galloprovincialis through hybridization with relatively little divergence accumulating between the shared RFLP haplotypes. We examined 80 mtl6S gene sequences for both the maternal and paternal mtDNA lineages from mussels sampled from various European populations and found strong support for the introgression hypothesis. In addition, we found that M. edulis mtDNA haplotypes appear to be introgressing into mussel populations in the Baltic Sea, which have predominantly M. trossulus nuclear genotypes, indicating that introgressive hybridization is prevalent among European mussel populations., (© 1998 The Society for the Study of Evolution.)

We have used the polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) techniques to design two genetic markers for blue mussels in the Mytilus edulis species complex. Both of these markers target the gene encoding the mussel polyphenolic adhesive protein. The first marker, Glu-5', is highly differentiated among and can be used to identify the three blue mussel species, M. edulis, M. galloprovincialis and M. trossulus. The second marker, Glu-3', can identify M. edulis and M. galloprovincialis. Using these markers we have demonstrated that hybrid mussels from Whitsand Bay, UK carry alleles for this gene that are the products of intragenic recombination. The high frequency (10 per cent) of these recombinant alleles within the hybrid population suggests that recombination is fairly frequent within this gene or that hybridization between M. edulis and M. galloprovincialis is substantial and has been occurring over considerable evolutionary time. The two novel genetic markers, Glu-5' and Glu-3' will be invaluable in additional studies regarding the importance of hybridization among blue mussels.

Blue mussels in the Mytilus edulis species complex have a doubly uniparental mode of mtDNA inheritance with separate maternal and paternal mtDNA lineages. Female mussels inherit their mtDNA solely from their mother, while males inherit mtDNA from both parents. In the male gonad the paternal mtDNA is preferentially replicated so that only paternal mtDNA is transmitted from fathers to sons. Hybridization is common among differentiated blue mussel taxa; whenever it involves M. trossulus, doubly uniparental mtDNA inheritance is disrupted. We have found high frequencies of males without and females with paternal mtDNA among hybrid mussels produced by interspecific matings between M. galloprovincialis and M. trossulus. In contrast, hybridization between M. galloprovincialis and M. edulis does not affect doubly uniparental inheritance, indicating a difference in the divergence of the mechanisms regulating mtDNA inheritance among the three blue mussel taxa. Our data indicate a high frequency of disrupted mtDNA transmission in F1 hybrids and suggest that two separate mechanisms, one regulating the transmission of paternal mtDNA to males and another inhibiting the establishment of paternal mtDNA in females, act to regulate doubly uniparental inheritance. We propose a model for the regulation of doubly uniparental inheritance that is consistent with these observations.

A novel form of mitochondrial DNA (mtDNA) inheritance has previously been documented for the blue mussel (Mytilus edulis). Female mussels inherit their mtDNA solely from their mother while males inherit mtDNA from both their mother and their father. In males, the paternal mtDNA is preferentially amplified so that the male gonad is highly enriched for the paternal mtDNA that is then transmitted from fathers to sons. We demonstrate that this mode of mtDNA inheritance also operates in the closely related species M. galloprovincialis and M. trossulus. The evolutionary relationship between the male and female mtDNA lineages is estimated by phylogenetic analysis of 455 nucleotides from the large subunit ribosomal RNA gene. We have found that the male and female lineages are highly divergent; the divergence of these lineages began prior to the speciation of the three species of blue mussels. Further, the separation between the male and female lineages is estimated to have occurred between 5.3 and 5.7 MYA.

Offspring from half-sib and full-sib families of the hard clam, Mercenaria mercenaria were reared in five locations along the Atlantic Coast to test for the presence of genotype-environment interaction for juvenile growth rate. Location effects upon growth rate variation were prevalent; of the genetic effects, the additive genetic by location variance was predominant with the nonadditive genetic by location component contributing to a lesser degree to the interaction variance. The additive and nonadditive variation over all environments was negligible. Genotype-environment interaction was found to be at least partially due to a change in the amount of genetic variation expressed at each location; with significant additive variation detected at Charleston and Georgetown, SC sites and significant nonadditive variation at Millsboro, DE. Genetic covariance/correlation analysis indicated that reversals in relative family performance across locations were prevalent, implying the possibility of habitat specialization among genotypes. In addition, graphical analysis produced no evidence of a ubiquitously superior genotype. These analyses suggest that genotype-environment interaction should act to constrain the evolution of juvenile growth rate in Mercenaria, preserve any heritable variation associated with this trait and may lead to the development of phenotypic plasticity for growth., (© 1991 The Society for the Study of Evolution.)