Patterns of dreissenid mussel invasions in western US lakes within an integrated gravity model framework.

Freshwater invasive species, such as the quagga mussel (Dreissena rostriformis bugensis), are causing over $1 billion USD annually in damages to water infrastructure, recreation, and the environment. Once established, quagga and other dreissenid mussels are extremely difficult to eradicate. Preventing the spread of these invasives is critical and of high management concern. Invasive dreissenid establishment is predicated upon both successful dispersal from a source and suitable habitat in the uninfested waterbody to which they are transported. Recreational boaters have become predominant dispersal vectors making it possible to forecast the risk of invasion of waterbodies for more targeted management and prevention. We developed an integrated mussel dispersal model that couples a constrained gravity model and habitat suitability model to forecast future invasions. The model simulates boater movement between lakes, the likelihood of boats transporting mussels, and the likelihood that those mussels survive in the environmental conditions of the new lake. Model output was most sensitive to changes in boater threshold, then buffer zones, while not as sensitive to changes in habitat suitability. From an initial infested source pool of 11 among 402 Western inland US lakes, we forecast additional lakes infested in several possible simulation scenarios. Constraining movement reduced connectivity between waterbodies with amplifying effects at different distance levels. This model can be used to determine waterbodies most at risk for dreissenid mussel invasion and to highlight the importance of multifactor integrated models in environmental management.

Positioning a multifunctional basic helix-loop-helix transcription factor within the Ciona notochord gene regulatory network.

In a multitude of organisms, transcription factors of the basic helix-loop-helix (bHLH) family control the expression of genes required for organ development and tissue differentiation. The functions of different bHLH transcription factors in the specification of nervous system and paraxial mesoderm have been widely investigated in various model systems. Conversely, the knowledge of the role of these regulators in the development of the axial mesoderm, the embryonic territory that gives rise to the notochord, and the identities of their target genes, remain still fragmentary. Here we investigated the transcriptional regulation and target genes of Bhlh-tun1, a bHLH transcription factor expressed in the developing Ciona notochord as well as in additional embryonic territories that contribute to the formation of both larval and adult structures. We describe its possible role in notochord formation, its relationship with the key notochord transcription factor Brachyury, and suggest molecular mechanisms through which Bhlh-tun1 controls the spatial and temporal expression of its effectors.

Clustered brachiopod Hox genes are not expressed collinearly and are associated with lophotrochozoan novelties.

Temporal collinearity is often considered the main force preserving Hox gene clusters in animal genomes. Studies that combine genomic and gene expression data are scarce, however, particularly in invertebrates like the Lophotrochozoa. As a result, the temporal collinearity hypothesis is currently built on poorly supported foundations. Here we characterize the complement, cluster, and expression of Hox genes in two brachiopod species, Terebratalia transversa and Novocrania anomalaT. transversa has a split cluster with 10 genes (lab, pb, Hox3, Dfd, Scr, Lox5, Antp, Lox4, Post2, and Post1), whereas N. anomala has 9 genes (apparently missing Post1). Our in situ hybridization, real-time quantitative PCR, and stage-specific transcriptomic analyses show that brachiopod Hox genes are neither strictly temporally nor spatially collinear; only pb (in T. transversa), Hox3 (in both brachiopods), and Dfd (in both brachiopods) show staggered mesodermal expression. Thus, our findings support the idea that temporal collinearity might contribute to keeping Hox genes clustered. Remarkably, expression of the Hox genes in both brachiopod species demonstrates cooption of Hox genes in the chaetae and shell fields, two major lophotrochozoan morphological novelties. The shared and specific expression of Hox genes, together with Arx, Zic, and Notch pathway components in chaetae and shell fields in brachiopods, mollusks, and annelids provide molecular evidence supporting the conservation of the molecular basis for these lophotrochozoan hallmarks.

Brachyury, Foxa2 and the cis-Regulatory Origins of the Notochord.

A main challenge of modern biology is to understand how specific constellations of genes are activated to differentiate cells and give rise to distinct tissues. This study focuses on elucidating how gene expression is initiated in the notochord, an axial structure that provides support and patterning signals to embryos of humans and all other chordates. Although numerous notochord genes have been identified, the regulatory DNAs that orchestrate development and propel evolution of this structure by eliciting notochord gene expression remain mostly uncharted, and the information on their configuration and recurrence is still quite fragmentary. Here we used the simple chordate Ciona for a systematic analysis of notochord cis-regulatory modules (CRMs), and investigated their composition, architectural constraints, predictive ability and evolutionary conservation. We found that most Ciona notochord CRMs relied upon variable combinations of binding sites for the transcription factors Brachyury and/or Foxa2, which can act either synergistically or independently from one another. Notably, one of these CRMs contains a Brachyury binding site juxtaposed to an (AC) microsatellite, an unusual arrangement also found in Brachyury-bound regulatory regions in mouse. In contrast, different subsets of CRMs relied upon binding sites for transcription factors of widely diverse families. Surprisingly, we found that neither intra-genomic nor interspecific conservation of binding sites were reliably predictive hallmarks of notochord CRMs. We propose that rather than obeying a rigid sequence-based cis-regulatory code, most notochord CRMs are rather unique. Yet, this study uncovered essential elements recurrently used by divergent chordates as basic building blocks for notochord CRMs.

Functional Brachyury binding sites establish a temporal read-out of gene expression in the Ciona notochord.

The appearance of the notochord represented a milestone in Deuterostome evolution. The notochord is necessary for the development of the chordate body plan and for the formation of the vertebral column and numerous organs. It is known that the transcription factor Brachyury is required for notochord formation in all chordates, and that it controls transcription of a large number of target genes. However, studies of the structure of the cis-regulatory modules (CRMs) through which this control is exerted are complicated in vertebrates by the genomic complexity and the pan-mesodermal expression territory of Brachyury. We used the ascidian Ciona, in which the single-copy Brachyury is notochord-specific and CRMs are easily identifiable, to carry out a systematic characterization of Brachyury-downstream notochord CRMs. We found that Ciona Brachyury (Ci-Bra) controls most of its targets directly, through non-palindromic binding sites that function either synergistically or individually to activate early- and middle-onset genes, respectively, while late-onset target CRMs are controlled indirectly, via transcriptional intermediaries. These results illustrate how a transcriptional regulator can efficiently shape a shallow gene regulatory network into a multi-tiered transcriptional output, and provide insights into the mechanisms that establish temporal read-outs of gene expression in a fast-developing chordate embryo.

Aryl hydrocarbon receptor (AHR) in the cnidarian Nematostella vectensis: comparative expression, protein interactions, and ligand binding.

The aryl hydrocarbon receptor (AHR) is a member of the basic helix-loop-helix/Per-ARNT-Sim (bHLH-PAS) family of transcription factors and has diverse roles in development, physiology, and environmental sensing in bilaterian animals. Studying the expression of conserved genes and function of proteins in outgroups to protostomes and deuterostomes assists in understanding the antiquity of gene function and deciphering lineage-specific differences in these bilaterian clades. We describe the developmental expression of AHR from the sea anemone Nematostella vectensis and compare its expression with three other members of the bHLH-PAS family (AHR nuclear translocator (ARNT), Cycle, and a proto-Single-Minded/Trachealess). NvAHR expression was highest early in the larval stage with spatial expression in the basal portion of the ectoderm that became increasingly restricted to the oral pole with concentrated expression in tentacles of the juvenile polyp. The other bHLH-PAS genes showed a divergent expression pattern in later larval stages and polyps, in which gene expression was concentrated in the aboral end, with broader expression in the endoderm later in development. In co-immunoprecipitation assays, we found no evidence for heterodimerization of AHR with ARNT, contrary to the conservation of this specific interaction in all bilaterians studied to date. Similar to results with other invertebrate AHRs but in contrast to vertebrate AHRs, NvAHR failed to bind two prototypical xenobiotic AHR ligands (2,3,7,8-tetrachlorodibenzo-p-dioxin, ß-naphthoflavone). Together, our data suggest that AHR's original function in Eumetazoa likely involved developmental patterning, potentially of neural tissue. The role of heterodimerization in the function of AHR may have arisen after the cnidarian-bilaterian ancestor. The absence of xenobiotic binding to NvAHR further supports a hypothesis for a derived role of this protein in chemical sensing within the chordates.

Toward invasive mussel genetic biocontrol: Approaches, challenges, and perspectives.

Invasive freshwater mussels, such as the zebra (Dreissena polymorpha), quagga (Dreissena rostriformis bugensis), and golden (Limnoperna fortunei) mussel have spread outside their native ranges throughout many regions of the North American, South American, and European continents in recent decades, damaging infrastructure and the environment. This review describes ongoing efforts by multiple groups to develop genetic biocontrol methods for invasive mussels. First, we provide an overview of genetic biocontrol strategies that have been applied in other invasive or pest species. Next, we summarize physical and chemical methods that are currently in use for invasive mussel control. We then describe the multidisciplinary approaches our groups are employing to develop genetic biocontrol tools for invasive mussels. Finally, we discuss the challenges and limitations of applying genetic biocontrol tools to invasive mussels. Collectively, we aim to openly share information and combine expertise to develop practical tools to enable the management of invasive freshwater mussels.

Cell proliferation is necessary for the regeneration of oral structures in the anthozoan cnidarian Nematostella vectensis.

BACKGROUND: The contribution of cell proliferation to regeneration varies greatly between different metazoan models. Planarians rely on pluripotent neoblasts and amphibian limb regeneration depends upon formation of a proliferative blastema, while regeneration in Hydra can occur in the absence of cell proliferation. Recently, the cnidarian Nematostella vectensis has shown potential as a model for studies of regeneration because of the ability to conduct comparative studies of patterning during embryonic development, asexual reproduction, and regeneration. The present study investigates the pattern of cell proliferation during the regeneration of oral structures and the role of cell proliferation in this process. RESULTS: In intact polyps, cell proliferation is observed in both ectodermal and endodermal tissues throughout the entire oral-aboral axis, including in the tentacles and physa. Following bisection, there is initially little change in proliferation at the wound site of the aboral fragment, however, beginning 18 to 24 hours after amputation there is a dramatic increase in cell proliferation at the wound site in the aboral fragment. This elevated level of proliferation is maintained throughout the course or regeneration of oral structures, including the tentacles, the mouth, and the pharynx. Treatments with the cell proliferation inhibitors hydroxyurea and nocodazole demonstrate that cell proliferation is indispensable for the regeneration of oral structures. Although inhibition of regeneration by nocodazole was generally irreversible, secondary amputation reinitiates cell proliferation and regeneration. CONCLUSIONS: The study has found that high levels of cell proliferation characterize the regeneration of oral structures in Nematostella, and that this cell proliferation is necessary for the proper progression of regeneration. Thus, while cell proliferation contributes to regeneration of oral structures in both Nematostella and Hydra, Nematostella lacks the ability to undergo the compensatory morphallactic mode of regeneration that characterizes Hydra. Our results are consistent with amputation activating a quiescent population of mitotically competent stem cells in spatial proximity to the wound site, which form the regenerated structures.

Temporal regulation of the muscle gene cascade by Macho1 and Tbx6 transcription factors in Ciona intestinalis.

For over a century, muscle formation in the ascidian embryo has been representative of 'mosaic' development. The molecular basis of muscle-fate predetermination has been partly elucidated with the discovery of Macho1, a maternal zinc-finger transcription factor necessary and sufficient for primary muscle development, and of its transcriptional intermediaries Tbx6b and Tbx6c. However, the molecular mechanisms by which the maternal information is decoded by cis-regulatory modules (CRMs) associated with muscle transcription factor and structural genes, and the ways by which a seamless transition from maternal to zygotic transcription is ensured, are still mostly unclear. By combining misexpression assays with CRM analyses, we have identified the mechanisms through which Ciona Macho1 (Ci-Macho1) initiates expression of Ci-Tbx6b and Ci-Tbx6c, and we have unveiled the cross-regulatory interactions between the latter transcription factors. Knowledge acquired from the analysis of the Ci-Tbx6b CRM facilitated both the identification of a related CRM in the Ci-Tbx6c locus and the characterization of two CRMs associated with the structural muscle gene fibrillar collagen 1 (CiFCol1). We use these representative examples to reconstruct how compact CRMs orchestrate the muscle developmental program from pre-localized ooplasmic determinants to differentiated larval muscle in ascidian embryos.

Direct activation of a notochord cis-regulatory module by Brachyury and FoxA in the ascidian Ciona intestinalis.

The notochord is a defining feature of the chordate body plan. Experiments in ascidian, frog and mouse embryos have shown that co-expression of Brachyury and FoxA class transcription factors is required for notochord development. However, studies on the cis-regulatory sequences mediating the synergistic effects of these transcription factors are complicated by the limited knowledge of notochord genes and cis-regulatory modules (CRMs) that are directly targeted by both. We have identified an easily testable model for such investigations in a 155-bp notochord-specific CRM from the ascidian Ciona intestinalis. This CRM contains functional binding sites for both Ciona Brachyury (Ci-Bra) and FoxA (Ci-FoxA-a). By combining point mutation analysis and misexpression experiments, we demonstrate that binding of both transcription factors to this CRM is necessary and sufficient to activate transcription. To gain insights into the cis-regulatory criteria controlling its activity, we investigated the organization of the transcription factor binding sites within the 155-bp CRM. The 155-bp sequence contains two Ci-Bra binding sites with identical core sequences but opposite orientations, only one of which is required for enhancer activity. Changes in both orientation and spacing of these sites substantially affect the activity of the CRM, as clusters of identical sites found in the Ciona genome with different arrangements are unable to activate transcription in notochord cells. This work presents the first evidence of a synergistic interaction between Brachyury and FoxA in the activation of an individual notochord CRM, and highlights the importance of transcription factor binding site arrangement for its function.

Strategic considerations for invasive species managers in the utilization of environmental DNA (eDNA): steps for incorporating this powerful surveillance tool.

Invasive species surveillance programs can utilize environmental DNA sampling and analysis to provide information on the presence of invasive species. Wider utilization of eDNA techniques for invasive species surveillance may be warranted. This paper covers topics directed towards invasive species managers and eDNA practitioners working at the intersection of eDNA techniques and invasive species surveillance. It provides background information on the utility of eDNA for invasive species management and points to various examples of its use across federal and international programs. It provides information on 1) why an invasive species manager should consider using eDNA, 2) deciding if eDNA can help with the manager's surveillance needs, 3) important components to operational implementation, and 4) a high-level overview of the technical steps necessary for eDNA analysis. The goal of this paper is to assist invasive species managers in deciding if, when, and how to use eDNA for surveillance. If eDNA use is elected, the paper provides guidance on steps to ensure a clear understanding of the strengths and limitation of the methods and how results can be best utilized in the context of invasive species surveillance.

Genes with spiralian-specific protein motifs are expressed in spiralian ciliary bands.

Spiralia is a large, ancient and diverse clade of animals, with a conserved early developmental program but diverse larval and adult morphologies. One trait shared by many spiralians is the presence of ciliary bands used for locomotion and feeding. To learn more about spiralian-specific traits we have examined the expression of 20 genes with protein motifs that are strongly conserved within the Spiralia, but not detectable outside of it. Here, we show that two of these are specifically expressed in the main ciliary band of the mollusc Tritia (also known as Ilyanassa). Their expression patterns in representative species from five more spiralian phyla-the annelids, nemerteans, phoronids, brachiopods and rotifers-show that at least one of these, lophotrochin, has a conserved and specific role in particular ciliated structures, most consistently in ciliary bands. These results highlight the potential importance of lineage-specific genes or protein motifs for understanding traits shared across ancient lineages.

Simple chordates exhibit a proton-independent function of acid-sensing ion channels.

Acid-sensing ion channels (ASICs) constitute a family of neuron-specific voltage-insensitive sodium channels gated by extracellular protons. Functions of ASICs in mammals include nociception, mechanosensation, and modulation of synaptic transmission. However, the role protons play in mediating the effects of ASICs remains elusive. We have examined ASICs from the ascidian Ciona intestinalis, a simple chordate organism whose nervous system in the larval stage exhibits high similarity to that of higher vertebrates. We found that the ascidian genome contains a single ASIC gene that gives rise to two splice forms analogous to the mammalian ASIC1 and ASIC2. CiASIC is expressed in most neurons of the larva but is absent in the adult. Despite high sequence similarity with mammalian counterparts, CiASIC is proton-insensitive when examined in heterologous systems or in larval neurons; the latter rules out the possibility that proton sensitivity is conferred by accessory proteins or particular factors present only in Ciona neurons. Down-regulation of the CiASIC transcript by double-stranded RNA disrupted the regular pattern of larval swimming, implying that proton-independent mechanisms mediate the effects of ASIC in vivo. Together the data identify ASIC as a highly conserved channel distinctive of chordate nervous systems and show that protons are not essential for ASIC function.

Molecular patterning during the development of Phoronopsis harmeri reveals similarities to rhynchonelliform brachiopods.

BACKGROUND: Phoronids, rhynchonelliform and linguliform brachiopods show striking similarities in their embryonic fate maps, in particular in their axis specification and regionalization. However, although brachiopod development has been studied in detail and demonstrated embryonic patterning as a causal factor of the gastrulation mode (protostomy vs deuterostomy), molecular descriptions are still missing in phoronids. To understand whether phoronids display underlying embryonic molecular mechanisms similar to those of brachiopods, here we report the expression patterns of anterior (otx, gsc, six3/6, nk2.1), posterior (cdx, bra) and endomesodermal (foxA, gata4/5/6, twist) markers during the development of the protostomic phoronid Phoronopsis harmeri. RESULTS: The transcription factors foxA, gata4/5/6 and cdx show conserved expression in patterning the development and regionalization of the phoronid embryonic gut, with foxA expressed in the presumptive foregut, gata4/5/6 demarcating the midgut and cdx confined to the hindgut. Furthermore, six3/6, usually a well-conserved anterior marker, shows a remarkably dynamic expression, demarcating not only the apical organ and the oral ectoderm, but also clusters of cells of the developing midgut and the anterior mesoderm, similar to what has been reported for brachiopods, bryozoans and some deuterostome Bilateria. Surprisingly, brachyury, a transcription factor often associated with gastrulation movements and mouth and hindgut development, seems not to be involved with these patterning events in phoronids. CONCLUSIONS: Our description and comparison of gene expression patterns with other studied Bilateria reveals that the timing of axis determination and cell fate distribution of the phoronid shows highest similarity to that of rhynchonelliform brachiopods, which is likely related to their shared protostomic mode of development. Despite these similarities, the phoronid Ph. harmeri also shows particularities in its development, which hint to divergences in the arrangement of gene regulatory networks responsible for germ layer formation and axis specification.

Conservation of notochord gene expression across chordates: insights from the Leprecan gene family.

The notochord is a defining character of the chordates, and the T-box transcription factor Brachyury has been shown to be required for notochord development in all chordates examined. In the ascidian Ciona intestinalis, at least 44 notochord genes have been identified as bona fide transcriptional targets of Brachyury. We examined the embryonic expression of a subset of murine orthologs of Ciona Brachyury target genes in the notochord to assess its conservation throughout chordate evolution. We focused on analyzing the Leprecan gene family, which in mouse is composed of three genes, as opposed to the single-copy Ciona gene. We found that all three mouse Leprecan genes are expressed in the notochord. Additionally, while Leprecan expression in C. intestinalis is confined to the notochord, expression of its mouse orthologs includes dorsal root ganglia, limb buds, branchial arches, and developing kidneys. These results have interesting implications for the evolution and development of chordates.

Evolutionary conservation of vertebrate notochord genes in the ascidian Ciona intestinalis.

To reconstruct a minimum complement of notochord genes evolutionarily conserved across chordates, we scanned the Ciona intestinalis genome using the sequences of 182 genes reported to be expressed in the notochord of different vertebrates and identified 139 candidate notochord genes. For 66 of these Ciona genes expression data were already available, hence we analyzed the expression of the remaining 73 genes and found notochord expression for 20. The predicted products of the newly identified notochord genes range from the transcription factors Ci-XBPa and Ci-miER1 to extracellular matrix proteins. We examined the expression of the newly identified notochord genes in embryos ectopically expressing Ciona Brachyury (Ci-Bra) and in embryos expressing a repressor form of this transcription factor in the notochord, and we found that while a subset of the genes examined are clearly responsive to Ci-Bra, other genes are not affected by alterations in its levels. We provide a first description of notochord genes that are not evidently influenced by the ectopic expression of Ci-Bra and we propose alternative regulatory mechanisms that might control their transcription.

Live imaging of fluorescent proteins in chordate embryos: from ascidians to mice.

Although we have advanced in our understanding of the molecular mechanisms intrinsic to the morphogenesis of chordate embryos, the question of how individual developmental events are integrated to generate the final morphological form is still unresolved. Microscopic observation is a pivotal tool in developmental biology, both for determining the normal course of events and for contrasting this with the results of experimental and pathological perturbations. Since embryonic development takes place in three dimensions over time, to fully understand the events required to build an embryo we must investigate embryo morphogenesis in multiple dimensions in situ. Recent advances in the isolation of naturally fluorescent proteins, and the refinement of techniques for in vivo microscopy offer unprecedented opportunities to study the cellular and molecular events within living, intact embryos using optical imaging. These technologies allow direct visual access to complex events as they happen in their native environment, and thus provide greater insights into cell behaviors operating during embryonic development. Since most fluorescent protein probes and modes of data acquisition are common across species, we have chosen the mouse and the ascidian, two model organisms at opposite ends of the chordate clade, to review the use of some of the current genetically-encoded fluorescent proteins and their visualization in vivo in living embryos for the generation of high-resolution imaging data.

The developmental basis for the recurrent evolution of deuterostomy and protostomy.

The mouth opening of bilaterian animals develops either separate from (deuterostomy) or connected to (protostomy) the embryonic blastopore, the site of endomesoderm internalization. Although this distinction preluded the classification of bilaterian animals in Deuterostomia and Protostomia, and has influenced major scenarios of bilaterian evolution, the developmental basis for the appearance of these different embryonic patterns remains unclear. To identify the underlying mechanisms, we compared the development of two brachiopod species that show deuterostomy (Novocrania anomala) and protostomy (Terebratalia transversa), respectively. We show that the differential activity of Wnt signalling, together with the timing and location of mesoderm formation, correlate with the differential behaviour and fate of the blastopore. We further assess these principles in the spiral-cleaving group Annelida, and propose that the developmental relationships of mouth and blastoporal openings are secondary by-products of variations in axial and mesoderm development. This challenges the previous evolutionary emphasis on extant blastoporal behaviours to explain the origin and diversification of bilaterian animals.

Mesodermal gene expression during the embryonic and larval development of the articulate brachiopod Terebratalia transversa.

BACKGROUND: Brachiopods undergo radial cleavage, which is distinct from the stereotyped development of closely related spiralian taxa. The mesoderm has been inferred to derive from the archenteron walls following gastrulation, and the primary mesoderm derivative in the larva is a complex musculature. To investigate the specification and differentiation of the mesoderm in the articulate brachiopod Terebratalia transversa, we have identified orthologs of genes involved in mesoderm development in other taxa and investigated their spatial and temporal expression during the embryonic and larval development of T. transversa. RESULTS: Orthologs of 17 developmental regulatory genes with roles in the development of the mesoderm in other bilaterian animals were found to be expressed in the developing mesoderm of T. transversa. Five genes, Tt.twist, Tt.GATA456, Tt.dachshund, Tt.mPrx, and Tt.NK1, were found to have expression throughout the archenteron wall at the radial gastrula stage, shortly after the initiation of gastrulation. Three additional genes, Tt.Pax1/9, Tt.MyoD, and Tt.Six1/2, showed expression at this stage in only a portion of the archenteron wall. Tt.eya, Tt.FoxC, Tt.FoxF, Tt.Mox, Tt.paraxis, Tt.Limpet, and Tt.Mef2 all showed initial mesodermal expression during later gastrula or early larval stages. At the late larval stage, Tt.dachshund, Tt.Limpet, and Tt.Mef2 showed expression in nearly all mesoderm cells, while all other genes were localized to specific regions of the mesoderm. Tt.FoxD and Tt.noggin both showed expression in the ventral mesoderm at the larval stages, with gastrula expression patterns in the archenteron roof and blastopore lip, respectively. CONCLUSIONS: Expression analyses support conserved roles for developmental regulators in the specification and differentiation of the mesoderm during the development of T. transversa. Expression of multiple mesodermal factors in the archenteron wall during gastrulation supports previous morphological observations that this region gives rise to larval mesoderm. Localized expression domains during gastrulation and larval development evidence early regionalization of the mesoderm and provide a basis for hypotheses regarding the molecular regulation underlying the complex system of musculature observed in the larva.

Glutathione S-Transferase (GST) Gene Diversity in the Crustacean Calanus finmarchicus--Contributors to Cellular Detoxification.

Detoxification is a fundamental cellular stress defense mechanism, which allows an organism to survive or even thrive in the presence of environmental toxins and/or pollutants. The glutathione S-transferase (GST) superfamily is a set of enzymes involved in the detoxification process. This highly diverse protein superfamily is characterized by multiple gene duplications, with over 40 GST genes reported in some insects. However, less is known about the GST superfamily in marine organisms, including crustaceans. The availability of two de novo transcriptomes for the copepod, Calanus finmarchicus, provided an opportunity for an in depth study of the GST superfamily in a marine crustacean. The transcriptomes were searched for putative GST-encoding transcripts using known GST proteins from three arthropods as queries. The identified transcripts were then translated into proteins, analyzed for structural domains, and annotated using reciprocal BLAST analysis. Mining the two transcriptomes yielded a total of 41 predicted GST proteins belonging to the cytosolic, mitochondrial or microsomal classes. Phylogenetic analysis of the cytosolic GSTs validated their annotation into six different subclasses. The predicted proteins are likely to represent the products of distinct genes, suggesting that the diversity of GSTs in C. finmarchicus exceeds or rivals that described for insects. Analysis of relative gene expression in different developmental stages indicated low levels of GST expression in embryos, and relatively high expression in late copepodites and adult females for several cytosolic GSTs. A diverse diet and complex life history are factors that might be driving the multiplicity of GSTs in C. finmarchicus, as this copepod is commonly exposed to a variety of natural toxins. Hence, diversity in detoxification pathway proteins may well be key to their survival.