Using a metabolomics approach to investigate the sensitivity of a potential Arctic-invader and its Arctic sister-species to marine heatwaves and traditional harvesting disturbances.

Coastal species are threatened by fishing practices and changing environmental conditions, such as marine heatwaves (MHW). The mechanisms that confer tolerance to such stressors in marine invertebrates are poorly understood. However, differences in tolerance among different species may be attributed to their geographical distribution. To test the tolerance of species occupying different thermal ranges, we used two closely related bivalves the softshell clam Mya arenaria (Linnaeus, 1758), a cold-temperate invader with demonstrated potential for establishment in the Arctic, and the blunt gaper Mya truncata (Linnaeus, 1758), a native polar species. Clams were subjected to a thermal stress, mimicking a MHW, and harvesting stress in a controlled environment. Seven acute temperature changes (2, 7, 12, 17, 22, 27, and 32 °C) were tested at two harvesting disturbance intensities (with, without). Survival was measured after 12 days and three tissues (gills, mantle, and posterior adductor muscle) collected from surviving individuals for targeted metabolomic profiling. MHW tolerance differed significantly between species: 26.9 °C for M. arenaria and 17.8 °C for M. truncata, with a negligeable effect of harvesting. At the upper thermal limit, M. arenaria displayed a more profound metabolomic remodelling when compared to M. truncata, and this varied greatly between tissue types. Network analysis revealed differences in pathway utilization at the upper MHW limit, with M. arenaria displaying a greater reliance on multiple DNA repair and expression and cell signalling pathways, while M. truncata was limited to fewer pathways. This suggests that M. truncata is ill equipped to cope with warming environments. MHW patterning in the Northwest Atlantic may be a strong predictor of population survival and future range shifts in these two clam species. As polar environments undergo faster rates of warming compared to the global average, M. truncata may be outcompeted by M. arenaria expanding into its native range.

Centuries of genome instability and evolution in soft-shell clam, Mya arenaria, bivalve transmissible neoplasia.

Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. We investigated the evolutionary history of a cancer lineage that has spread though the soft-shell clam (Mya arenaria) population by assembling a chromosome-scale soft-shell clam reference genome and characterizing somatic mutations in transmissible cancer. We observe high mutation density, widespread copy-number gain, structural rearrangement, loss of heterozygosity, variable telomere lengths, mitochondrial genome expansion and transposable element activity, all indicative of an unstable cancer genome. We also discover a previously unreported mutational signature associated with overexpression of an error-prone polymerase and use this to estimate the lineage to be >200 years old. Our study reveals the ability for an invertebrate cancer lineage to survive for centuries while its genome continues to structurally mutate, likely contributing to the evolution of this lineage as a parasitic cancer.

The Genome of the Softshell Clam Mya arenaria and the Evolution of Apoptosis.

Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

Population Genetic Structure is Unrelated to Shell Shape, Thickness and Organic Content in European Populations of the Soft-Shell Clam Mya Arenaria.

The soft-shell clam Mya arenaria is one of the most ancient invaders of European coasts and is present in many coastal ecosystems, yet little is known about its genetic structure in Europe. We collected 266 samples spanning a latitudinal cline from the Mediterranean to the North Sea and genotyped them at 12 microsatellite loci. In parallel, geometric morphometric analysis of shell outlines was used to test for associations between shell shape, latitude and genotype, and for a selection of shells we measured the thickness and organic content of the granular prismatic (PR), the crossed-lamellar (CL) and the complex crossed-lamellar (CCL) layers. Strong population structure was detected, with Bayesian cluster analysis identifying four groups located in the Mediterranean, Celtic Sea, along the continental coast of the North Sea and in Scotland. Multivariate analysis of shell shape uncovered a significant effect of collection site but no associations with any other variables. Shell thickness did not vary significantly with either latitude or genotype, although PR thickness and calcification were positively associated with latitude, while CCL thickness showed a negative association. Our study provides new insights into the population structure of this species and sheds light on factors influencing shell shape, thickness and microstructure.

Biogeography of resistance to paralytic shellfish toxins in softshell clam, Mya arenaria (L.), populations along the Atlantic coast of North America.

Blooms of Alexandrium spp., the causative agent of paralytic shellfish poisoning (PSP), recur with varying frequency and intensity on the Northwest Atlantic coast of North America, from New York, USA, to northern Canadian waters. Along this latitudinal range blooms co-occur with abundant, intertidal populations of softshell clams, Mya arenaria. Prior work identified a naturally-occurring genetic mutation in Domain II α-subunit of the clams' voltage-gated sodium channels (NaV), which significantly reduces the binding affinity of the paralytic shellfish toxin, saxitoxin (STX). This mutation provides clams with resistance to the deleterious effects of STX, allowing them to continue feeding during Alexandrium spp. blooms and attain very high tissue toxicities. This study used genetic sequencing of the NaV mutation locus in clams from four coastal regions of the Bay of Fundy-Gulf of Maine and the mid-Atlantic to determine the percentage of clams in each region that possess the resistant NaV mutation. The genotype composition was related to the occurrence and magnitude of PSP outbreaks based on shellfish toxicity, primarily that of mussels, Mytilus edulis, used as a proxy for the prevalence and severity of Alexandrium blooms in each region. As hypothesized, the proportion of clams bearing the resistant mutation generally matched up well with the historical incidence and intensity of Alexandrium spp. blooms. The highest percentage of homozygote resistant clams (RR = 70.0%), and the lowest percentage of sensitive clams (SS = 4.5%) were found in eastern Gulf of Maine populations. Exceptions at a few sites where anomalously high numbers of M. arenaria with the resistant mutation were found despite the absence of blooms, may be attributable to larval gene flow. There was no evidence that Alexandrium blooms occurring in Northport Harbor, Long Island, have resulted in a shift in genotypic composition of the local clam population, presumably due to their low cell toxicity. Seasonal mismatch of highly vulnerable M. arenaria postset with toxic blooms at this latitude may also partly explain this result. This study provides strong supporting evidence that Alexandrium blooms can select for resistance to PSP-toxins in M. arenaria populations and proposes a mechanism for the persistence of the sensitive allele throughout the region. Implications for clam aquaculture (seeding) efforts, as well as for shellfish toxicity monitoring are discussed.

Cellular stress responses to chronic heat shock and shell damage in temperate Mya truncata.

Acclimation, via phenotypic flexibility, is a potential means for a fast response to climate change. Understanding the molecular mechanisms underpinning phenotypic flexibility can provide a fine-scale cellular understanding of how organisms acclimate. In the last 30 years, Mya truncata populations around the UK have faced an average increase in sea surface temperature of 0.7 °C and further warming of between 1.5 and 4 °C, in all marine regions adjacent to the UK, is predicted by the end of the century. Hence, data are required on the ability of M. truncata to acclimate to physiological stresses, and most notably, chronic increases in temperature. Animals in the present study were exposed to chronic heat-stress for 2 months prior to shell damage and subsequently, only 3, out of 20 damaged individuals, were able to repair their shells within 2 weeks. Differentially expressed genes (between control and damaged animals) were functionally enriched with processes relating to cellular stress, the immune response and biomineralisation. Comparative transcriptomics highlighted genes, and more broadly molecular mechanisms, that are likely to be pivotal in this lack of acclimation. This study demonstrates that discovery-led transcriptomic profiling of animals during stress-response experiments can shed light on the complexity of biological processes and changes within organisms that can be more difficult to detect at higher levels of biological organisation.

Shell matrix proteins of the clam, Mya truncata: Roles beyond shell formation through proteomic study.

Mya truncata, a soft shell clam, is presented as a new model to study biomineralization through a proteomics approach. In this study, the shell and mantle tissue were analysed in order to retrieve knowledge about the secretion of shell matrix proteins (SMPs). Out of 67 and 127 shell and mantle proteins respectively, 16 were found in both shell and mantle. Bioinformatic analysis of SMP sequences for domain prediction revealed the presence of several new domains such as fucolectin tachylectin-4 pentraxin-1 (FTP), scavenger receptor, alpha-2-macroglobulin (α2 M), lipocalin and myosin tail along with previously reported SMP domains such as chitinase, carbonic anhydrase, tyrosinase, sushi, and chitin binding. Interestingly, these newly predicted domains are attributed with molecular functions other than biomineralization. These findings suggest that shells may not only act as protective armour from predatory action, but could also actively be related to other functions such as immunity. In this context, the roles of SMPs in biomineralization need to be looked in a new perspective.

Characterisation of the mantle transcriptome and biomineralisation genes in the blunt-gaper clam, Mya truncata.

Members of the Myidae family are ecologically and economically important, but there is currently very little molecular data on these species. The present study sequenced and assembled the mantle transcriptome of Mya truncata from the North West coast of Scotland and identified candidate biomineralisation genes. RNA-Seq reads were assembled to create 20,106 contigs in a de novo transciptome, 18.81% of which were assigned putative functions using BLAST sequence similarity searching (cuttoff E-value 1E-10). The most highly expressed genes were compared to the Antarctic clam (Laternula elliptica) and showed that many of the dominant biological functions (muscle contraction, energy production, biomineralisation) in the mantle were conserved. There were however, differences in the constitutive expression of heat shock proteins, which were possibly due to the M. truncata sampling location being at a relatively low latitude, and hence relatively warm, in terms of the global distribution of the species. Phylogenetic analyses of the Tyrosinase proteins from M. truncata showed a gene expansion which was absent in L. elliptica. The tissue distribution expression patterns of putative biomineralisation genes were investigated using quantitative PCR, all genes showed a mantle specific expression pattern supporting their hypothesised role in shell secretion. The present study provides some preliminary insights into how clams from different environments - temperate versus polar - build their shells. In addition, the transcriptome data provides a valuable resource for future comparative studies investigating biomineralisation.

Mining the transcriptomes of four commercially important shellfish species for single nucleotide polymorphisms within biomineralization genes.

Transcriptional profiling not only provides insights into patterns of gene expression, but also generates sequences that can be mined for molecular markers, which in turn can be used for population genetic studies. As part of a large-scale effort to better understand how commercially important European shellfish species may respond to ocean acidification, we therefore mined the transcriptomes of four species (the Pacific oyster Crassostrea gigas, the blue mussel Mytilus edulis, the great scallop Pecten maximus and the blunt gaper Mya truncata) for single nucleotide polymorphisms (SNPs). Illumina data for C. gigas, M. edulis and P. maximus and 454 data for M. truncata were interrogated using GATK and SWAP454 respectively to identify between 8267 and 47,159 high quality SNPs per species (total=121,053 SNPs residing within 34,716 different contigs). We then annotated the transcripts containing SNPs to reveal homology to diverse genes. Finally, as oceanic pH affects the ability of organisms to incorporate calcium carbonate, we honed in on genes implicated in the biomineralization process to identify a total of 1899 SNPs in 157 genes. These provide good candidates for biomarkers with which to study patterns of selection in natural or experimental populations.

Complete mitochondrial genome of the soft-shell clam Mya arenaria.

We have sequenced and characterized the complete mitochondrial genome of the soft-shell clam, Mya arenaria, an important organism for environmental toxicology and aquaculture. Mya arenaria is located in the taxonomic order Myoida, which lacks any member with a completely annotated mitogenome. The M. arenaria mitochondrial genome is 17 947 bp in length. Like most marine bivalves, the circular mitogenome codes entirely on the heavy strand, with no introns. As with other bivalves, the gene order of the mitochondrion is highly rearranged. The mitogenome contains 12 protein-coding genes but ATP8 is missing, consistent with about half of all bivalve genera. Twenty-three tRNAs were identified. Phylogenetic analysis shows that M. arenaria is related most closely with the bivalves Sinonovacula constricta, and Moerella iridescens, of the infraclass Euheterodonta (unassigned). This, along with the close grouping of the phylogenetic trees, confirms a close tie between Myoida and Euheterodonta (unassigned).

Horizontal transmission of clonal cancer cells causes leukemia in soft-shell clams.

Outbreaks of fatal leukemia-like cancers of marine bivalves throughout the world have led to massive population loss. The cause of the disease is unknown. We recently identified a retrotransposon, Steamer, that is highly expressed and amplified to high copy number in neoplastic cells of soft-shell clams (Mya arenaria). Through analysis of Steamer integration sites, mitochondrial DNA single-nucleotide polymorphisms (SNPs), and polymorphic microsatellite alleles, we show that the genotypes of neoplastic cells do not match those of the host animal. Instead, neoplastic cells from dispersed locations in New York, Maine, and Prince Edward Island (PEI), Canada, all have nearly identical genotypes that differ from those of the host. These results indicate that the cancer is spreading between animals in the marine environment as a clonal transmissible cell derived from a single original clam. Our findings suggest that horizontal transmission of cancer cells is more widespread in nature than previously supposed.

Tracking larvae with molecular markers reveals high relatedness and early seasonal recruitment success in a partially spawning marine bivalve.

The partial synchronized spawning strategy adopted by some marine invertebrate broadcast-spawners can lead to the production of many distinct pools of larvae within a single reproductive cycle. Following the fate of these larval groups from birth to settlement with molecular markers might shed light on mechanisms regulating their population recruitment. Larvae and recruits of Mya arenaria, a partially spawning marine bivalve, were monitored and collected over 13 consecutive weeks during an entire reproductive cycle. Each sampled individual (n = 218) was sorted according to size (early veligers, late veligers, post-larval recruits) and genotyped at seven microsatellite loci for comparisons among samples and with adult reference samples (n = 270). While traditional differentiation statistics (e.g., pairwise Δ(ST), allelic richness) suggested the absence of sweepstakes reproductive success, the level of relatedness found within and among larvae and recruit samples suggested otherwise. Four samples out of ten were observed to have higher within-sample relatedness values than randomly expected, including the very first group of early veligers produced in the season (E1) and the last group of post-larvae who survived recruitment (P10). E1 individuals were also found to be more related than randomly expected to individuals of more than 80% of all other samples including the last surviving recruits (P8 and P10). These results suggest that the first larvae produced in the season were the most successful to survive recruitment. Results also show direct evidence for larval retention and demonstrate for the first time larval and post-larval kin aggregation in a marine bivalve.

Activation of transcription and retrotransposition of a novel retroelement, Steamer, in neoplastic hemocytes of the mollusk Mya arenaria.

Bivalve mollusks of the North Atlantic, most prominently the soft shell clam Mya arenaria, are afflicted with an epidemic transmissible disease of the circulatory system closely resembling leukemia. The disease is characterized by a dramatic expansion of blast-like cells in the hemolymph with high mitotic index. Examination of hemolymph of diseased clams revealed high levels of reverse transcriptase activity, the hallmark of retroviruses and retroelements. By deep sequencing of RNAs from hemolymph, we identified transcripts of a novel retroelement, here named Steamer. The DNA of the element is marked by long terminal repeats and encodes a single large protein with similarity to mammalian retroviral Gag-Pol proteins. Steamer mRNA levels were specifically elevated in diseased hemocytes, and high expression was correlated with disease status. DNA copy number per genome was present at enormously high levels in diseased hemocytes, indicative of extensive reverse transcription and retrotransposition. Steamer activation in M. arenaria is an example of a catastrophic induction of genetic instability that may initiate or advance the course of leukemia.

A mu class glutathione S-transferase from Manila clam Ruditapes philippinarum (RpGSTµ): cloning, mRNA expression, and conjugation assays.

Glutathione S-transferases (GSTs) are enzymes that catalyze xenobiotic metabolism in the phase II detoxification process. GSTs have a potential for use as indicators or biomarkers to assess the presence of organic and inorganic contaminants in aquatic environments. In this study, a full-length cDNA of a mu (µ) class GST (RpGSTµ) was identified from Manila clam (Ruditapes philippinarum) and biochemically characterized. The 1356 bp of the cDNA included an open reading frame of 651 bp encoding a polypeptide of 217 amino acid residues with a molecular mass of 25.04 kDa and an estimated pI of 6.34. Sequence analysis revealed that the RpGSTµ possessed several characteristic features of µ class GSTs, such as a thioredoxin-like N-terminal domain containing binding sites for glutathione (GSH), a C-terminal domain containing substrate binding sites, and a µ loop. The recombinant RpGSTµ (rRpGSTµ) protein exhibited GSH-conjugating catalytic activity towards several substrates, and significantly strong activity was detected against 4-nitrophenethyl bromide (5.77 ± 0.55) and 1-chloro-2,4-dinitrobenzene (CDNB, 3.19 ± 0.05). Kinetic analysis as a function of GSH and CDNB concentrations revealed relatively low Km values of 1.03 ± 0.46 mM and 0.56 ± 0.20 mM, respectively, thereby indicating a GSH-conjugation attributed with high rates. The optimum pH and temperature for the catalytic activity of the rRpGSTµ protein were 7.7 and 37°C, respectively. The effect of two inhibitors, Cibacron blue and hematin, on the activity of rRpGSTµ was evaluated and the IC50 values of 0.65 µM and 9 µM, respectively, were obtained. While RpGSTµ transcripts were highly expressed in gills and hemocytes, a significant elevation in mRNA levels was detected in these tissues after lipopolysaccharide (LPS), polyinosinic-polycytidylic acid (poly I:C) and live bacterial (Vibrio tapetis) challenges. These findings collectively suggest that RpGSTµ functions as a potent detoxifier of xenobiotic toxicants present in the aquatic environment, and that its mRNA expression could be modulated by pathogenic stress signal(s).

DNA integrity assessment in hemocytes of soft-shell clams (Mya arenaria) in the Saguenay Fjord (Québec, Canada).

The purpose of this study was to examine the effects of pollution on DNA integrity in the feral soft-shell clam (Mya arenaria) in the Saguenay Fjord. Intertidal clams were collected downstream and upstream of the fjord at sites under anthropogenic pollution. DNA integrity was assessed by following changes in single- and double-stranded breaks, variation in DNA content and micro-nuclei (MN) incidence in hemocytes. The results revealed that clams collected at polluted sites had reduced DNA strand breaks (lower DNA repair activity), increased DNA content variation and MN frequency in hemocytes. The data revealed that DNA content variation was closely related to MN frequency and negatively with DNA strand breaks formation. Water conductivity was also related to reduced MN frequency and DNA content variation, indicating that, in addition to the effects of pollution, the gradual dilution of saltwater could compromise mussel health.

Isolation and characterization of new microsatellite markers for the invasive softshell clam, Mya arenaria (L.) (Bivalvia: Myidae).

The invasive softshell clam (Mya arenaria Linnaeus, 1758) is native to the northwestern region of the Atlantic Ocean. This species has been introduced in the northeast Pacific and along the European coasts, due to intense naval transports and aquaculture, and it is now present in all the European seas. In this paper we describe seven new microsatellite loci for Mya arenaria. The isolated loci are polymorphic with a number of alleles per locus between 6 and 14. The observed and expected heterozygosities ranged from 0.417 to 0.951, and from 0.643 to 0.895, with an average of 0.716 and 0.775, respectively. These microsatellite markers should be useful in analyzing this species' genetic diversity, which could explain various processes of its invasion history.

Effects of pesticide compounds (chlorothalonil and mancozeb) and benzo[a]pyrene mixture on aryl hydrocarbon receptor, p53 and ubiquitin gene expression levels in haemocytes of soft-shell clams (Mya arenaria).

The aim of this study is to investigate the effects of the pesticides/polycyclic aromatic hydrocarbon mixture on aryl hydrocarbon receptor (AhR), p53 and ubiquitin mRNA level in haemocytes of Mya arenaria exposed to a mixture of chlorothalonil, mancozeb and benzo[a]pyrene (BaP) for 48 and 72 h. AhR, p53 and ubiquitin gene expression levels were quantified using quantitative Real-time PCR. For robust and accurate quantification of transcripts, suitable housekeeping genes were selected from four sets of ribosomal and elongation factors transcripts previously sequenced from Mya arenaria using geNorm open source software. Quantitative Real-time PCR data exhibited a significantly high expression of AhR after 72 h of exposure (P ≤ 0.05). p53 gene expression seems to be up-regulated by the mixture after 48 h, however not significantly; but the level of p53 mRNA is down-regulated by the xenobiotics between 48 and 72 h after exposure. This study postulates that AhR mRNA levels could be used as an indicator of the exposure of clams' haemocytes to a mixture of xenobiotics such as chlorothalonil, mancozeb and BaP. However, further studies have to be pursued in order to unravel the molecular mechanisms involved in the p53 signaling pathway.

Induction of transposase and polyprotein RNA levels in disseminated neoplastic hemocytes of soft-shell clams: Mya arenaria.

In Prince Edward Island, a high mortality of soft-shell clams Mya arenaria was found to be related to the disease known as disseminated neoplasia (DN). However, the molecular mechanisms by which hemocytes of clams are transformed in the course of DN remain by far unknown. This study aims at identifying the transcripts involved in the development of the disease. Four subtractive cDNA sequence libraries were generated and more than 200,000 reads were obtained. Following similarity searches in genome databases, the transcripts were assigned to cellular functions including mitochondrial respiration, structural proteins, cytoskeleton, nucleic acid regulation, general metabolism, signal transduction, apoptosis, cell cycle regulation, as well as virus transcripts. The expression levels of transposase and polyprotein genes were evaluated in clams with various percentages of tetraploid hemocytes. Data have shown that expression levels were significantly higher in clams with a high percentage of tetraploid hemocytes. These results reinforce the hypothesis of endogenous retrotransposon involvement in the etiology of the disease. Further investigations are needed, however, to elucidate the role of transposase and polyprotein in the disease development.

Identification and expression of immune-related genes in hemocytes of soft-shell clams, Mya arenaria, challenged with Vibrio splendidus.

Although the mollusc immune system has been studied at the cellular level, the response to pathogens at gene expression level has not been thoroughly investigated. This study aimed to investigate the early molecular response of hemocytes of soft-shell clams, Mya arenaria, to Vibrio splendidus strain LGP32 by identification of transcripts involved in immune defense. The Suppression Subtractive Hybridization (SSH) was used to selectively identify differentially expressed genes in hemocytes exposed to V. splendidus at a ratio 1:1 for 2 h. Both forward and reverse subtracted cDNA were constructed and a total of 16,000 reads were obtained and analyzed. Identity searches in genome databases were performed using BlastX program and transcripts were clustered to cellular functions including structural proteins, immunity, stress proteins, apoptosis, cell process, metabolism and signal transduction. Among the differentially expressed immune associated genes were ficolin, killer cell lectin-like receptor, natural resistance-associated macrophage protein 1 (Nramp-1), mitogen-activated protein kinases (MAPK), ferritin, heat shock proteins 90 (HSP90) and cathepsin and their expressions were quantified using Reverse Transcription quantitative Polymerase Chain Reaction (RT-qPCR) at 1, 2 and 3 h post-Vibrio challenge. These genes showed similar expression patterns, up-regulation at 1 h, followed by a down-regulation at 2 and 3 h. These data corroborates our previous observations of cell rounding, reduced phagocytosis and respiratory burst activity. To our knowledge, this is the first study to demonstrate an effect of V. splendidus on expression of genes related to immune system in soft-shell clams M. arenaria. However, further investigations are needed to unravel the molecular mechanisms of hemocytes subjected to V. splendidus.

Differential gene expression of gamma-actin, Toll-like receptor 2 (TLR-2) and interleukin-1 receptor-associated kinase 4 (IRAK-4) in Mya arenaria haemocytes induced by in vivo infections with two Vibrio splendidus strains.

Immune function gene expression in Mya arenaria haemocytes was evaluated following in vivo infection with Vibrio splendidus LGP32-GFP and 7SHRW. Elongation factor 1alpha (EF-1alpha) with 2 (EF-2), after challenge with LGP32-GFP, and EF-1alpha with the ribosomal protein S-18, after challenge with 7SHRW, were found to be the most stable housekeeping genes. Using these internal controls and comparing the regulation induced by both strains, up-regulation of gamma-actin, down-regulation of TLR-2 and up-regulation of IRAK-4 was significantly higher after challenge with LGP32-GFP (p<0.001, p=0.001 and p<0.05, respectively). These results suggest specific responses at a molecular level modulated by the bacterial strains. LGP32-GFP induced marked responses which coincide with a similar trend previously found on phenotypic responses under our experimental model.