Unlocking the secret life of blue mussels: Exploring connectivity in the Skagerrak through biophysical modeling and population genomics.

Knowledge of functional dispersal barriers in the marine environment can be used to inform a wide variety of management actions, such as marine spatial planning, restoration efforts, fisheries regulations, and invasive species management. Locations and causes of dispersal barriers can be studied through various methods, including movement tracking, biophysical modeling, demographic models, and genetics. Combining methods illustrating potential dispersal, such as biophysical modeling, with realized dispersal through, e.g., genetic connectivity estimates, provides particularly useful information for teasing apart potential causes of observed barriers. In this study, we focus on blue mussels (Mytilus edulis) in the Skagerrak-a marginal sea connected to the North Sea in Northern Europe-and combine biophysical models of larval dispersal with genomic data to infer locations and causes of dispersal barriers in the area. Results from both methods agree; patterns of ocean currents are a major structuring factor in the area. We find a complex pattern of source-sink dynamics with several dispersal barriers and show that some areas can be isolated despite an overall high dispersal capability. Finally, we translate our finding into management advice that can be used to sustainably manage this ecologically and economically important species in the future.

Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s).

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

A metabarcoding analysis of the wrackbed microbiome indicates a phylogeographic break along the North Sea-Baltic Sea transition zone.

Sandy beaches are biogeochemical hotspots that bridge marine and terrestrial ecosystems via the transfer of organic matter, such as seaweed (termed wrack). A keystone of this unique ecosystem is the microbial community, which helps to degrade wrack and re-mineralize nutrients. However, little is known about this community. Here, we characterize the wrackbed microbiome as well as the microbiome of a primary consumer, the seaweed fly Coelopa frigida, and examine how they change along one of the most studied ecological gradients in the world, the transition from the marine North Sea to the brackish Baltic Sea. We found that polysaccharide degraders dominated both microbiomes, but there were still consistent differences between wrackbed and fly samples. Furthermore, we observed a shift in both microbial communities and functionality between the North and Baltic Sea driven by changes in the frequency of different groups of known polysaccharide degraders. We hypothesize that microbes were selected for their abilities to degrade different polysaccharides corresponding to a shift in polysaccharide content in the different seaweed communities. Our results reveal the complexities of both the wrackbed microbial community, with different groups specialized to different roles, and the cascading trophic consequences of shifts in the near shore algal community.

Global, regional, and cryptic population structure in a high gene-flow transatlantic fish.

Lumpfish (Cyclopterus lumpus) is a transatlantic marine fish displaying large population sizes and a high potential for dispersal and gene-flow. These features are expected to result in weak population structure. Here, we investigated population genetic structure of lumpfish throughout its natural distribution in the North Atlantic using two approaches: I) 4,393 genome wide SNPs and 95 individuals from 10 locations, and II) 139 discriminatory SNPs and 1,669 individuals from 40 locations. Both approaches identified extensive population genetic structuring with a major split between the East and West Atlantic and a distinct Baltic Sea population, as well as further differentiation of lumpfish from the English Channel, Iceland, and Greenland. The discriminatory loci displayed ~2-5 times higher divergence than the genome wide approach, revealing further evidence of local population substructures. Lumpfish from Isfjorden in Svalbard were highly distinct but resembled most fish from Greenland. The Kattegat area in the Baltic transition zone, formed a previously undescribed distinct genetic group. Also, further subdivision was detected within North America, Iceland, West Greenland, Barents Sea, and Norway. Although lumpfish have considerable potential for dispersal and gene-flow, the observed high levels of population structuring throughout the Atlantic suggests that this species may have a natal homing behavior and local populations with adaptive differences. This fine-scale population structure calls for consideration when defining management units for exploitation of lumpfish stocks and in decisions related to sourcing and moving lumpfish for cleaner fish use in salmonid aquaculture.

Ten years of marine evolutionary biology-Challenges and achievements of a multidisciplinary research initiative.

The Centre for Marine Evolutionary Biology (CeMEB) at the University of Gothenburg, Sweden, was established in 2008 through a 10-year research grant of 8.7 m€ to a team of senior researchers. Today, CeMEB members have contributed >500 scientific publications, 30 PhD theses and have organised 75 meetings and courses, including 18 three-day meetings and four conferences. What are the footprints of CeMEB, and how will the centre continue to play a national and international role as an important node of marine evolutionary research? In this perspective article, we first look back over the 10 years of CeMEB activities and briefly survey some of the many achievements of CeMEB. We furthermore compare the initial goals, as formulated in the grant application, with what has been achieved, and discuss challenges and milestones along the way. Finally, we bring forward some general lessons that can be learnt from a research funding of this type, and we also look ahead, discussing how CeMEB's achievements and lessons can be used as a springboard to the future of marine evolutionary biology.

Population dynamics of Baltic herring since the Viking Age revealed by ancient DNA and genomics.

The world's oceans are currently facing major stressors in the form of overexploitation and anthropogenic climate change. The Baltic Sea was home to the first "industrial" fishery ∼800 y ago targeting the Baltic herring, a species that is still economically and culturally important today. Yet, the early origins of marine industries and the long-term ecological consequences of historical and contemporary fisheries remain debated. Here, we study long-term population dynamics of Baltic herring to evaluate the past impacts of humans on the marine environment. We combine modern whole-genome data with ancient DNA (aDNA) to identify the earliest-known long-distance herring trade in the region, illustrating that extensive fish trade began during the Viking Age. We further resolve population structure within the Baltic and observe demographic independence for four local herring stocks over at least 200 generations. It has been suggested that overfishing at Øresund in the 16th century resulted in a demographic shift from autumn-spawning to spring-spawning herring dominance in the Baltic. We show that while the Øresund fishery had a negative impact on the western Baltic herring stock, the demographic shift to spring-spawning dominance did not occur until the 20th century. Instead, demographic reconstructions reveal population trajectories consistent with expected impacts of environmental change and historical reports on shifting fishing targets over time. This study illustrates the joint impact of climate change and human exploitation on marine species as well as the role historical ecology can play in conservation and management policies.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species.

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double-digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long-term (~30-year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation-by-distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

Stabilizing selection on Atlantic cod supergenes through a millennium of extensive exploitation.

Life on Earth has been characterized by recurring cycles of ecological stasis and disruption, relating biological eras to geological and climatic transitions through the history of our planet. Due to the increasing degree of ecological abruption caused by human influences many advocate that we now have entered the geological era of the Anthropocene, or "the age of man." Considering the ongoing mass extinction and ecosystem reshuffling observed worldwide, a better understanding of the drivers of ecological stasis will be a requisite for identifying routes of intervention and mitigation. Ecosystem stability may rely on one or a few keystone species, and the loss of such species could potentially have detrimental effects. The Atlantic cod (Gadus morhua) has historically been highly abundant and is considered a keystone species in ecosystems of the northern Atlantic Ocean. Collapses of cod stocks have been observed on both sides of the Atlantic and reported to have detrimental effects that include vast ecosystem reshuffling. By whole-genome resequencing we demonstrate that stabilizing selection maintains three extensive "supergenes" in Atlantic cod, linking these genes to species persistence and ecological stasis. Genomic inference of historic effective population sizes shows continued declines for cod in the North Sea-Skagerrak-Kattegat system through the past millennia, consistent with an early onset of the marine Anthropocene through industrialization and commercialization of fisheries throughout the medieval period.

Not that clean: Aquaculture-mediated translocation of cleaner fish has led to hybridization on the northern edge of the species' range.

Translocation and introduction of non-native organisms can have major impacts on local populations and ecosystems. Nevertheless, translocations are common practices in agri- and aquaculture. Each year, millions of wild-caught wrasses are transported large distances to be used as cleaner fish for parasite control in marine salmon farms. Recently, it was documented that translocated cleaner fish are able to escape and reproduce with local wild populations. This is especially a challenge in Norway, which is the world's largest salmon producer. Here, a panel of 84 informative SNPs was developed to identify the presence of nonlocal corkwing wrasse (Symphodus melops) escapees and admixed individuals in wild populations in western Norway. Applying this panel to ~2000 individuals, escapees and hybrids were found to constitute up to 20% of the local population at the northern edge of the species' distribution. The introduction of southern genetic material at the northern edge of the species distribution range has altered the local genetic composition and could obstruct local adaptation and further range expansion. Surprisingly, in other parts of the species distribution where salmon farming is also common, few escapees and hybrids were found. Why hybridization seems to be common only in the far north is discussed in the context of demographic and transport history. However, the current lack of reporting of escapes makes it difficult to evaluate possible causes for why some aquaculture-dense areas have more escapees and hybrids than others. The results obtained in this study, and the observed high genomic divergence between the main export and import regions, puts the sustainability of mass translocation of nonlocal wild wrasse into question and suggests that the current management regime needs re-evaluation.

Post-glacial establishment of locally adapted fish populations over a steep salinity gradient.

Studies of colonization of new habitats that appear from rapidly changing environments are interesting and highly relevant to our understanding of divergence and speciation. Here, we analyse phenotypic and genetic variation involved in the successful establishment of a marine fish (sand goby, Pomatoschistus minutus) over a steep salinity drop from 35 PSU in the North Sea (NE Atlantic) to two PSU in the inner parts of the post-glacial Baltic Sea. We first show that populations are adapted to local salinity in a key reproductive trait, the proportion of motile sperm. Thereafter, we show that genome variation at 22,190 single nucleotide polymorphisms (SNPs) shows strong differentiation among populations along the gradient. Sequences containing outlier SNPs and transcriptome sequences, mapped to a draft genome, reveal associations with genes with relevant functions for adaptation in this environment but without overall evidence of functional enrichment. The many contigs involved suggest polygenic differentiation. We trace the origin of this differentiation using demographic modelling and find the most likely scenario is that at least part of the genetic differentiation is older than the Baltic Sea and is a result of isolation of two lineages prior to the current contact over the North Sea-Baltic Sea transition zone.

Genome wide analysis reveals genetic divergence between Goldsinny wrasse populations.

BACKGROUND: Marine fish populations are often characterized by high levels of gene flow and correspondingly low genetic divergence. This presents a challenge to define management units. Goldsinny wrasse (Ctenolabrus rupestris) is a heavily exploited species due to its importance as a cleaner-fish in commercial salmonid aquaculture. However, at the present, the population genetic structure of this species is still largely unresolved. Here, full-genome sequencing was used to produce the first genomic reference for this species, to study population-genomic divergence among four geographically distinct populations, and, to identify informative SNP markers for future studies. RESULTS: After construction of a de novo assembly, the genome was estimated to be highly polymorphic and of ~600Mbp in size. 33,235 SNPs were thereafter selected to assess genomic diversity and differentiation among four populations collected from Scandinavia, Scotland, and Spain. Global FST among these populations was 0.015-0.092. Approximately 4% of the investigated loci were identified as putative global outliers, and ~ 1% within Scandinavia. SNPs showing large divergence (FST > 0.15) were picked as candidate diagnostic markers for population assignment. One hundred seventy-three of the most diagnostic SNPs between the two Scandinavian populations were validated by genotyping 47 individuals from each end of the species' Scandinavian distribution range. Sixty-nine of these SNPs were significantly (p < 0.05) differentiated (mean FST_173_loci = 0.065, FST_69_loci = 0.140). Using these validated SNPs, individuals were assigned with high probability (≥ 94%) to their populations of origin. CONCLUSIONS: Goldsinny wrasse displays a highly polymorphic genome, and substantial population genomic structure. Diversifying selection likely affects population structuring globally and within Scandinavia. The diagnostic loci identified now provide a promising and cost-efficient tool to investigate goldsinny wrasse populations further.

A Darwinian Laboratory of Multiple Contact Zones.

Barriers to gene flow between divergent populations result in contact (hybrid) zones. Locations where multiple contact zones overlap can be used in comparative studies asking: what mechanisms maintain barriers; what is the origin of the genetic variation involved; and do differences in life history affect the nature of barriers? In a review of 23 marine species' genetic divergence over a postglacial salinity gradient, many showed steep genetic clines supported by divergent selection and/or temporal or spatial segregation. Contacts were primary or secondary and shaped by ancestral variation sometimes involving inversions. The dispersal potential of species seemed less important in shaping clines. Studies of multispecies contact zones will increase our understanding of speciation, but we need to address the taxonomic bias and focus more on postzygotic isolation.

"A cleaner break": Genetic divergence between geographic groups and sympatric phenotypes revealed in ballan wrasse (Labrus bergylta).

Capture and long-distance translocation of cleaner fish to control lice infestations on marine salmonid farms has the potential to influence wild populations via overexploitation in source regions, and introgression in recipient regions. Knowledge of population genetic structure is therefore required. We studied the genetic structure of ballan wrasse, a phenotypically diverse and extensively used cleaner fish, from 18 locations in Norway and Sweden, and from Galicia, Spain, using 82 SNP markers. We detected two very distinct genetic groups in Scandinavia, northwestern and southeastern. These groups were split by a stretch of sandy beaches in southwest Norway, representing a habitat discontinuity for this rocky shore associated benthic egg-laying species. Wrasse from Galicia were highly differentiated from all Scandinavian locations, but more similar to northwestern than southeastern locations. Distinct genetic differences were observed between sympatric spotty and plain phenotypes in Galicia, but not in Scandinavia. The mechanisms underlying the geographic patterns between phenotypes are discussed, but not identified. We conclude that extensive aquaculture-mediated translocation of ballan wrasse from Sweden and southern Norway to western and middle Norway has the potential to mix genetically distinct populations. These results question the sustainability of the current cleaner fish practice.

Spatial genetic structure in a crustacean herbivore highlights the need for local considerations in Baltic Sea biodiversity management.

Incorporating species' eco-evolutionary responses to human-caused disturbances remains a challenge in marine management efforts. A prerequisite is knowledge of geographic structure and scale of genetic diversity and connectivity-the so-called seascape genetic patterns. The Baltic Sea is an excellent model system for studies linking seascape genetics with effects of anthropogenic stress. However, seascape genetic patterns in this area are only described for a few species and are completely unknown for invertebrate herbivores, which constitute a critical part of the ecosystem. This information is crucial for sustainable management, particularly under future scenarios of rapid environmental change. Here, we investigate the population genetic structure among 31 locations throughout the Baltic Sea, of which 45% were located in marine protected areas, in one of the most important herbivores of this region, the isopod crustacean Idotea balthica, using an array of 33,774 genome-wide SNP markers derived from 2b-RAD sequencing. In addition, we generate a biophysical connectivity matrix for I. balthica from a combination of oceanographic current models and estimated life history traits. We find population structure on scales of hundreds of kilometers across the Baltic Sea, where genomic patterns in most cases closely match biophysical connectivity, indicating passive transport with oceanographic currents as an important mean of dispersal in this species. We also find a reduced genetic diversity in terms of heterozygosity along the main salinity gradient of the Baltic Sea, suggesting periods of low population size. Our results provide crucial information for the management of a key ecosystem species under expected changes in temperature and salinity following global climate change in a marine coastal area.

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe.

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (FST ~0.1), over a short coastline (<60 km) in the North Sea-Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post-glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.

Assessing SNP-markers to study population mixing and ecological adaptation in Baltic cod.

Atlantic cod (Gadus morhua) is a species of great ecological and economical importance in the Baltic Sea. Here, two genetically differentiated stocks, the western and the eastern Baltic cod, display substantial mechanical mixing, hampering our understanding of cod ecology and impeding stock assessments and management. Based on whole-genome re-sequencing data from reference samples obtained from the study area, we designed two different panels of Single Nucleotide Polymorphisms markers (SNPs), which take into account the exceptional genome architecture of cod. A minimum panel of 20 diagnostic SNPs and an extended panel (20 diagnostic and 18 biologically informative SNPs, 38 in total) were developed and validated to distinguish unambiguously between the western and the eastern Baltic cod stocks and to enable studies of local adaptation to the specific environment in the Baltic Sea, respectively. We tested both panels on cod sampled from the southern Baltic Sea (n = 603) caught in 2015 and 2016. Genotyping results showed that catches from the mixing zone in the Arkona Sea, were composed of similar proportions of individuals of the western and the eastern stock. Catches from adjacent areas to the east, the Bornholm Basin and Gdansk Deep, were exclusively composed of eastern Baltic cod, whereas catches from adjacent western areas (Belt Sea and Öresund) were composed of western Baltic cod. Interestingly, the two Baltic cod stocks showed strong genetic differences at loci associated with life-history trait candidate genes, highlighting the species' potential for ecological adaptation even at small geographical scales. The minimum and the extended panel of SNP markers presented in this study provide powerful tools for future applications in research and fisheries management to further illuminate the mixing dynamics of cod in the Baltic Sea and to better understand Baltic cod ecology.

Disentangling structural genomic and behavioural barriers in a sea of connectivity.

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole-genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord-type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord-type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.

Safety assessment of EPA+DHA canola oil by fatty acid profile comparison to various edible oils and fat-containing foods and a 28-day repeated dose toxicity study in rats.

The omega-3 long-chain polyunsaturated fatty acids (LC-PUFAs) eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) are recognized for their health-promoting qualities. Marine fish and fish oil currently provide the main sources of EPA and DHA for human consumption. An alternative plant-based source of EPA and DHA is provided by EPA + DHA canola event LBFLFK (LBFLFK). A comparative analysis and a 28-day toxicity study assessed the safety of LBFLFK refined, bleached, and deodorized (RBD) oil. Thirty-one different commercially-obtained fat and oil samples were tested, and principal component analysis showed that the overall fatty acid profile of LBFLFK RBD oil was most similar to Mortierella alpina oil and salmon flesh. Samples with the fewest differences in the presence or absence of individual fatty acids compared to LBFLFK RBD oil were menhaden oil and some other fish oils. In a 28-day toxicity study, LBFLFK RBD oil was administered by oral gavage to male and female Wistar rats. No signs of toxicity were evident and no adverse effects were noted in clinical observations, clinical pathology, or histopathology. Overall, these studies support the safety of LBFLFK RBD oil as a source of EPA and DHA for human consumption.

Population genomics of parallel evolution in gene expression and gene sequence during ecological adaptation.

Natural selection often produces parallel phenotypic changes in response to a similar adaptive challenge. However, the extent to which parallel gene expression differences and genomic divergence underlie parallel phenotypic traits and whether they are decoupled or not remains largely unexplored. We performed a population genomic study of parallel ecological adaptation among replicate ecotype pairs of the rough periwinkle (Littorina saxatilis) at a regional geographical scale (NW Spain). We show that genomic changes underlying parallel phenotypic divergence followed a complex pattern of both repeatable differences and of differences unique to specific ecotype pairs, in which parallel changes in expression or sequence are restricted to a limited set of genes. Yet, the majority of divergent genes were divergent either for gene expression or coding sequence, but not for both simultaneously. Overall, our findings suggest that divergent selection significantly contributed to the process of parallel molecular differentiation among ecotype pairs, and that changes in expression and gene sequence underlying phenotypic divergence could, at least to a certain extent, be considered decoupled processes.

Diet-dependent gene expression highlights the importance of Cytochrome P450 in detoxification of algal secondary metabolites in a marine isopod.

Isopods of the genus Idotea have an unusual ability to feed on algae containing high amounts of chemical defense molecules, such as species of the genera Fucus and Ulva. In this study, we compared gene expression patterns of Idotea balthica individuals fed with Fucus vesiculosus to individuals fed with Ulva lactuca. We generated the first-ever transcriptome assembly for this species, and found 3,233 differentially expressed genes across feeding regimes. However, only a handful of biological functions were enriched with regard to differentially expressed genes, the most notable being "alkaloid metabolic process". Within this category, we found eight differentially expressed cytochrome P450 (CYP) unigenes, all of which had a higher expression in the U. lactuca diet treatment. A phylogenetic analysis showed that the differentially expressed CYP genes are closely related to a CYP gene described from the hepatopancreas of the spiny lobster Panulirus argus, and we hypothesize that these transcripts are involved in metabolite detoxification. This is a first step in the understanding of this algae-grazer interaction, and will form a basis for future work to characterize cytochrome P450 functioning in marine crustaceans.