Epigenetic variations are more substantial than genetic variations in rapid adaptation of oyster to Pacific oyster mortality syndrome.

Disease emergence is accelerating with global changes. Understanding by which mechanisms host populations can rapidly adapt will be crucial for management practices. Pacific oyster mortality syndrome (POMS) imposes a substantial and recurrent selective pressure on oyster populations, and rapid adaptation may arise through genetics and epigenetics. In this study, we used (epi)genome-wide association mapping to show that oysters differentially exposed to POMS displayed genetic and epigenetic signatures of selection. Consistent with higher resistance to POMS, the genes targeted included many genes in several pathways related to immunity. By combining correlation, DNA methylation quantitative trait loci, and variance partitioning, we revealed that a third of phenotypic variation was explained by interactions between the genetic and epigenetic information, ~14% by the genome, and up to 25% by the epigenome alone. Similar to genetically based adaptation, epigenetic mechanisms notably governing immune responses can contribute substantially to the rapid adaptation of hosts to emerging infectious diseases.

Cooperation and cheating orchestrate Vibrio assemblages and polymicrobial synergy in oysters infected with OsHV-1 virus.

Polymicrobial infections threaten the health of humans and animals but remain understudied in natural systems. We recently described the Pacific Oyster Mortality Syndrome (POMS), a polymicrobial disease affecting oyster production worldwide. In the French Atlantic coast, the disease involves coinfection with ostreid herpesvirus 1 (OsHV-1) and virulent Vibrio. However, it is unknown whether consistent Vibrio populations are associated with POMS in different regions, how Vibrio contribute to POMS, and how they interact with OsHV-1 during pathogenesis. By connecting field-based approaches in a Mediterranean ecosystem, laboratory infection assays and functional genomics, we uncovered a web of interdependencies that shape the structure and function of the POMS pathobiota. We show that Vibrio harveyi and Vibrio rotiferianus are predominant in OsHV-1-diseased oysters and that OsHV-1 drives the partition of the Vibrio community observed in the field. However only V. harveyi synergizes with OsHV-1 by promoting mutual growth and accelerating oyster death. V. harveyi shows high-virulence potential and dampens oyster cellular defenses through a type 3 secretion system, making oysters a more favorable niche for microbe colonization. In addition, V. harveyi produces a key siderophore called vibrioferrin. This important resource promotes the growth of V. rotiferianus, which cooccurs with V. harveyi in diseased oysters, and behaves as a cheater by benefiting from V. harveyi metabolite sharing. Our data show that cooperative behaviors contribute to synergy between bacterial and viral coinfecting partners. Additional cheating behaviors further shape the polymicrobial consortium. Controlling cooperative behaviors or countering their effects opens avenues for mitigating polymicrobial diseases.

A core of functional complementary bacteria infects oysters in Pacific Oyster Mortality Syndrome.

BACKGROUND: The Pacific oyster Crassostrea gigas is one of the main cultivated invertebrate species worldwide. Since 2008, oyster juveniles have been confronted with a lethal syndrome known as the Pacific Oyster Mortality Syndrome (POMS). POMS is a polymicrobial disease initiated by a primary infection with the herpesvirus OsHV-1 µVar that creates an oyster immunocompromised state and evolves towards a secondary fatal bacteremia. RESULTS: In the present article, we describe the implementation of an unprecedented combination of metabarcoding and metatranscriptomic approaches to show that the sequence of events in POMS pathogenesis is conserved across infectious environments. We also identified a core bacterial consortium which, together with OsHV-1 µVar, forms the POMS pathobiota. This bacterial consortium is characterized by high transcriptional activities and complementary metabolic functions to exploit host's resources. A significant metabolic specificity was highlighted at the bacterial genus level, suggesting low competition for nutrients between members of the core bacteria. CONCLUSIONS: Lack of metabolic competition between the core bacteria might favor complementary colonization of host tissues and contribute to the conservation of the POMS pathobiota across distinct infectious environments.

Early life microbial exposures shape the Crassostrea gigas immune system for lifelong and intergenerational disease protection.

BACKGROUND: The interaction of organisms with their surrounding microbial communities influences many biological processes, a notable example of which is the shaping of the immune system in early life. In the Pacific oyster, Crassostrea gigas, the role of the environmental microbial community on immune system maturation - and, importantly, protection from infectious disease - is still an open question. RESULTS: Here, we demonstrate that early life microbial exposure durably improves oyster survival when challenged with the pathogen causing Pacific oyster mortality syndrome (POMS), both in the exposed generation and in the subsequent one. Combining microbiota, transcriptomic, genetic, and epigenetic analyses, we show that the microbial exposure induced changes in epigenetic marks and a reprogramming of immune gene expression leading to long-term and intergenerational immune protection against POMS. CONCLUSIONS: We anticipate that this protection likely extends to additional pathogens and may prove to be an important new strategy for safeguarding oyster aquaculture efforts from infectious disease. tag the videobyte/videoabstract in this section Video Abstract.

Seaweeds influence oyster microbiota and disease susceptibility.

A growing awareness of role that microbiota can play in mediating the effects of pathogens on hosts has given rise to the concept of the pathobiome. Recently, we demonstrated that the Pacific oyster mortality syndrome affecting Crassostrea gigas oysters is caused by infection with the Ostreid herpesvirus type 1 (OsHV-1) followed by infection with multiple bacterial taxa. Here we extend the concept of this pathobiome beyond the host species and its bacterial microbiota by investigating how seaweed living in association with oysters influences their response to the disease. We hypothesized that by their mere presence in the environment, different species of seaweeds can positively or negatively influence the risk of disease in oysters by shaping their bacterial microbiota and their immune response. Although seaweed and oysters do not have direct ecological interactions, they are connected by seawater and likely share microbes. To test our hypothesis, oysters were acclimated with green, brown or red algae for 2 weeks and then challenged with OsHV-1. We monitored host survival and pathogen proliferation and performed bacterial microbiota and transcriptome analyses. We found that seaweeds can alter the bacterial microbiota of the host and its response to the disease. More particularly, green algae belonging to the genus Ulva spp. induced bacterial microbiota dysbiosis in oyster and modification of its transcriptional immune response leading to increased susceptibility to the disease. This work provides a better understanding of a marine disease and highlights the importance of considering both macrobiotic and microbiotic interactions for conservation, management and exploitation of marine ecosystems and resources.

Recent advances in bivalve-microbiota interactions for disease prevention in aquaculture.

In bivalves, no clear-cut functional role of microbiota has yet been identified, although many publications suggest that they could be involved in nutrition or immunity of their host. In the context of climate change, integrative approaches at the crossroads of disciplines have been developed to explore the environment-host-pathogen-microbiota system. Here, we attempt to synthesize work on (1) the current methodologies to analyse bivalve microbiota, (2) the comparison of microbiota between species, between host compartments and their surrounding habitat, (3) how the bivalve microbiota are governed by environmental factors and host genetics and (4) how host-associated microorganisms act as a buffer against pathogens and/or promote recovery, and could thereby play a role in the prevention of disease or mortalities.

In situ characterisation of pathogen dynamics during a Pacific oyster mortality syndrome episode.

Significant mortality of Crassostrea gigas juveniles is observed systematically every year worldwide. Pacific Oyster Mortality Syndrome (POMS) is caused by Ostreid Herpesvirus 1 (OsHV-1) infection leading to immune suppression, followed by bacteraemia caused by a consortium of opportunistic bacteria. Using an in-situ approach and pelagic chambers, our aim in this study was to identify pathogen dynamics in oyster flesh and in the water column during the course of a mortality episode in the Mediterranean Thau lagoon (France). OsHV-1 concentrations in oyster flesh increased before the first clinical symptoms of the disease appeared, reached maximum concentrations during the moribund phase and the mortality peak. The structure of the bacterial community associated with oyster flesh changed in favour of bacterial genera previously associated with oyster mortality including Vibrio, Arcobacter, Psychrobium, and Psychrilyobacter. During the oyster mortality episode, releases of OsHV-1 and opportunistic bacteria were observed, in succession, in the water surrounding the oyster lanterns. These releases may favour the spread of disease within oyster farms and potentially impact other marine species, thereby reducing marine biodiversity in shellfish farming areas.

Contribution of Viral Genomic Diversity to Oyster Susceptibility in the Pacific Oyster Mortality Syndrome.

Juvenile Pacific oysters (Crassostrea gigas) are subjected to recurrent episodes of mass mortalities that constitute a threat for the oyster industry. This mortality syndrome named "Pacific Oyster Mortality Syndrome" (POMS) is a polymicrobial disease whose pathogenesis is initiated by a primary infection by a variant of an Ostreid herpes virus named OsHV-1 µVar. The characterization of the OsHV-1 genome during different disease outbreaks occurring in different geographic areas has revealed the existence of a genomic diversity for OsHV-1 µVar. However, the biological significance of this diversity is still poorly understood. To go further in understanding the consequences of OsHV-1 diversity on POMS, we challenged five biparental families of oysters to two different infectious environments on the French coasts (Atlantic and Mediterranean). We observed that the susceptibility to POMS can be different among families within the same environment but also for the same family between the two environments. Viral diversity analysis revealed that Atlantic and Mediterranean POMS are caused by two distinct viral populations. Moreover, we observed that different oyster families are infected by distinct viral populations within a same infectious environment. Altogether these results suggest that the co-evolutionary processes at play between OsHV-1 µVar and oyster populations have selected a viral diversity that could facilitate the infection process and the transmission in oyster populations. These new data must be taken into account in the development of novel selective breeding programs better adapted to the oyster culture environment.

Microbiota Composition and Evenness Predict Survival Rate of Oysters Confronted to Pacific Oyster Mortality Syndrome.

Pacific Oyster Mortality Syndrome (POMS) affects Crassostrea gigas oysters worldwide and causes important economic losses. Disease dynamic was recently deciphered and revealed a multiple and progressive infection caused by the Ostreid herpesvirus OsHV-1 µVar, triggering an immunosuppression followed by microbiota destabilization and bacteraemia by opportunistic bacterial pathogens. However, it remains unknown if microbiota might participate to protect oysters against POMS, and if microbiota characteristics might be predictive of oyster mortalities. To tackle this issue, we transferred full-sib progenies of resistant and susceptible oyster families from hatchery to the field during a period in favor of POMS. After 5 days of transplantation, oysters from each family were either sampled for individual microbiota analyses using 16S rRNA gene-metabarcoding or transferred into facilities to record their survival using controlled condition. As expected, all oysters from susceptible families died, and all oysters from the resistant family survived. Quantification of OsHV-1 and bacteria showed that 5 days of transplantation were long enough to contaminate oysters by POMS, but not for entering the pathogenesis process. Thus, it was possible to compare microbiota characteristics between resistant and susceptible oysters families at the early steps of infection. Strikingly, we found that microbiota evenness and abundances of Cyanobacteria (Subsection III, family I), Mycoplasmataceae, Rhodobacteraceae, and Rhodospirillaceae were significantly different between resistant and susceptible oyster families. We concluded that these microbiota characteristics might predict oyster mortalities.

A Sustained Immune Response Supports Long-Term Antiviral Immune Priming in the Pacific Oyster, Crassostrea gigas.

Over the last decade, innate immune priming has been evidenced in many invertebrate phyla. If mechanistic models have been proposed, molecular studies aiming to substantiate these models have remained scarce. We reveal here the transcriptional signature associated with immune priming in the oyster Crassostrea gigas Oysters were fully protected against Ostreid herpesvirus 1 (OsHV-1), a major oyster pathogen, after priming with poly(I·C), which mimics viral double-stranded RNA. Global analysis through RNA sequencing of oyster and viral genes after immune priming and viral infection revealed that poly(I·C) induces a strong antiviral response that impairs OsHV-1 replication. Protection is based on a sustained upregulation of immune genes, notably genes involved in the interferon pathway and apoptosis, which control subsequent viral infection. This persistent antiviral alert state remains active over 4 months and supports antiviral protection in the long term. This acquired resistance mechanism reinforces the molecular foundations of the sustained response model of immune priming. It further opens the way to applications (pseudovaccination) to cope with a recurrent disease that causes dramatic economic losses in the shellfish farming industry worldwide.IMPORTANCE In the last decade, important discoveries have shown that resistance to reinfection can be achieved without a functional adaptive immune system, introducing the concept of innate immune memory in invertebrates. However, this field has been constrained by the limited number of molecular mechanisms evidenced to support these phenomena. Taking advantage of an invertebrate species, the Pacific oyster (Crassostrea gigas), in which we evidenced one of the longest and most effective periods of protection against viral infection observed in an invertebrate, we provide the first comprehensive transcriptomic analysis of antiviral innate immune priming. We show that priming with poly(I·C) induced a massive upregulation of immune-related genes, which control subsequent viral infection, and it was maintained for over 4 months after priming. This acquired resistant mechanism reinforces the molecular foundations of the sustained response model of immune priming. It opens the way to pseudovaccination to prevent the recurrent diseases that currently afflict economically or ecologically important invertebrates.

Differential basal expression of immune genes confers Crassostrea gigas resistance to Pacific oyster mortality syndrome.

BACKGROUND: As a major threat to the oyster industry, Pacific Oyster Mortality Syndrome (POMS) is a polymicrobial disease affecting the main oyster species farmed across the world. POMS affects oyster juveniles and became panzootic this last decade, but POMS resistance in some oyster genotypes has emerged. While we know some genetic loci associated with resistance, the underlying mechanisms remained uncharacterized. So, we developed a comparative transcriptomic approach using basal gene expression profiles between different oyster biparental families with contrasted phenotypes when confronted to POMS (resistant or susceptible). RESULTS: We showed that POMS resistant oysters show differential expression of genes involved in stress responses, protein modifications, maintenance of DNA integrity and repair, and immune and antiviral pathways. We found similarities and clear differences among different molecular pathways in the different resistant families. These results suggest that the resistance process is polygenic and partially varies according to the oyster genotype. CONCLUSIONS: We found differences in basal expression levels of genes related to TLR-NFκB, JAK-STAT and STING-RLR pathways. These differences could explain the best antiviral response, as well as the robustness of resistant oysters when confronted to POMS. As some of these genes represent valuable candidates for selective breeding, we propose future studies should further examine their function.

Immune-suppression by OsHV-1 viral infection causes fatal bacteraemia in Pacific oysters.

Infectious diseases are mostly explored using reductionist approaches despite repeated evidence showing them to be strongly influenced by numerous interacting host and environmental factors. Many diseases with a complex aetiology therefore remain misunderstood. By developing a holistic approach to tackle the complexity of interactions, we decipher the complex intra-host interactions underlying Pacific oyster mortality syndrome affecting juveniles of Crassostrea gigas, the main oyster species exploited worldwide. Using experimental infections reproducing the natural route of infection and combining thorough molecular analyses of oyster families with contrasted susceptibilities, we demonstrate that the disease is caused by multiple infection with an initial and necessary step of infection of oyster haemocytes by the Ostreid herpesvirus OsHV-1 µVar. Viral replication leads to the host entering an immune-compromised state, evolving towards subsequent bacteraemia by opportunistic bacteria. We propose the application of our integrative approach to decipher other multifactorial diseases that affect non-model species worldwide.

Inefficient immune response is associated with microbial permissiveness in juvenile oysters affected by mass mortalities on field.

Since 2008, juvenile Crassostrea gigas oysters have suffered from massive mortalities in European farming areas. This disease of complex etiology is still incompletely understood. Triggered by an elevated seawater temperature, it has been associated to infections by a herpes virus named OsHV-1 as well as pathogenic vibrios of the Splendidus clade. Ruling out the complexity of the disease, most of our current knowledge has been acquired in controlled experiments. Among the many unsolved questions, it is still ignored what role immunity plays in the capacity oysters have to survive an infectious episode. Here we show that juvenile oysters susceptible to the disease mount an inefficient immune response associated with microbial permissiveness and death. We found that, in contrast to resistant adult oysters having survived an earlier episode of mortality, susceptible juvenile oysters never exposed to infectious episodes died by more than 90% in a field experiment. Susceptible oysters were heavily colonized by OsHV-1 herpes virus as well as bacteria including vibrios potentially pathogenic for oysters, which proliferated in oyster flesh and body fluids during the mortality event. Nonetheless, susceptible oysters were found to sense microbes as indicated by an overexpression of immune receptors and immune signaling pathways. However, they did not express important immune effectors involved in antimicrobial immunity and apoptosis and showed repressed expression of genes involved in ROS and metal homeostasis. This contrasted with resistant oysters, which expressed those important effectors, controlled bacterial and viral colonization and showed 100% survival to the mortality event. Altogether, our results demonstrate that the immune response mounted by susceptible oysters lacks some important immune functions and fails in controlling microbial proliferation. This study opens the way to more holistic studies on the "mass mortality syndrome", which are now required to decipher the sequence of events leading to oyster mortalities and determine the relative weight of pathogens, oyster genetics and oyster-associated microbiota in the disease.

Response of the pearl oyster Pinctada margaritifera to cadmium and chromium: Identification of molecular biomarkers.

This study was designed to identify in the pearl oyster Pinctada margaritifera, used as a bio-accumulator, molecular biomarkers for the presence of heavy metals in the lagoon environment. Pearl oysters were exposed to 2 concentrations (1 and 10µgL-1) of cadmium (Cd) and chromium (Cr) compared to a control. Twelve target genes encoding proteins potentially involved in the response to heavy metal contamination with antioxidant, detoxification or apoptosis activities were selected. P. margaritifera accumulated Cd but not Cr, and mortality was related to the amount of Cd accumulated in tissues. In response to Cd-Cr contamination, metallothionein (MT) was significantly up-regulated by Cd-Cr at both concentrations, while 7 others (SOD, CAT, GPX, GSTO, GSTM, CASP, MDR) were down-regulated. Based on the development of these molecular tools, we propose that the pearl oyster, P. margaritifera, could be used as a sentinel species for heavy metal contamination in the lagoons of tropical ecosystems.

An updated assessment of Symbiodinium spp. that associate with common scleractinian corals from Moorea (French Polynesia) reveals high diversity among background symbionts and a novel finding of clade B.

The adaptative bleaching hypothesis (ABH) states that, depending on the symbiotic flexibility of coral hosts (i.e., the ability of corals to "switch" or "shuffle" their algal symbionts), coral bleaching can lead to a change in the composition of their associated Symbiodinium community and, thus, contribute to the coral's overall survival. In order to determine the flexibility of corals, molecular tools are required to provide accurate species delineations and to detect low levels of coral-associated Symbiodinium. Here, we used highly sensitive quantitative (real-time) PCR (qPCR) technology to analyse five common coral species from Moorea (French Polynesia), previously screened using only traditional molecular methods, to assess the presence of low-abundance (background) Symbiodinium spp. Similar to other studies, each coral species exhibited a strong specificity to a particular clade, irrespective of the environment. In addition, however, each of the five species harboured at least one additional Symbiodinium clade, among clades A-D, at background levels. Unexpectedly, and for the first time in French Polynesia, clade B was detected as a coral symbiont. These results increase the number of known coral-Symbiodinium associations from corals found in French Polynesia, and likely indicate an underestimation of the ability of the corals in this region to associate with and/or "shuffle" different Symbiodinium clades. Altogether our data suggest that corals from French Polynesia may favor a trade-off between optimizing symbioses with a specific Symbiodinium clade(s), maintaining associations with particular background clades that may play a role in the ability of corals to respond to environmental change.

Bacterial community characterization of water and intestine of the shrimp Litopenaeus stylirostris in a biofloc system.

BACKGROUND: Biofloc technology (BFT), a rearing method with little or no water exchange, is gaining popularity in aquaculture. In the water column, such systems develop conglomerates of microbes, algae and protozoa, together with detritus and dead organic particles. The intensive microbial community presents in these systems can be used as a pond water quality treatment system, and the microbial protein can serve as a feed additive. The current problem with BFT is the difficulty of controlling its bacterial community composition for both optimal water quality and optimal shrimp health. The main objective of the present study was to investigate microbial diversity of samples obtained from different culture environments (Biofloc technology and clear seawater) as well as from the intestines of shrimp reared in both environments through high-throughput sequencing technology. RESULTS: Analyses of the bacterial community identified in water from BFT and "clear seawater" (CW) systems (control) containing the shrimp Litopenaeus stylirostris revealed large differences in the frequency distribution of operational taxonomic units (OTUs). Four out of the five most dominant bacterial communities were different in both culture methods. Bacteria found in great abundance in BFT have two principal characteristics: the need for an organic substrate or nitrogen sources to grow and the capacity to attach to surfaces and co-aggregate. A correlation was found between bacteria groups and physicochemical and biological parameters measured in rearing tanks. Moreover, rearing-water bacterial communities influenced the microbiota of shrimp. Indeed, the biofloc environment modified the shrimp intestine microbiota, as the low level (27 %) of similarity between intestinal bacterial communities from the two treatments. CONCLUSION: This study provides the first information describing the complex biofloc microbial community, which can help to understand the environment-microbiota-host relationship in this rearing system.

Antimicrobial peptides in marine invertebrate health and disease.

Aquaculture contributes more than one-third of the animal protein from marine sources worldwide. A significant proportion of aquaculture products are derived from marine protostomes that are commonly referred to as 'marine invertebrates'. Among them, penaeid shrimp (Ecdysozosoa, Arthropoda) and bivalve molluscs (Lophotrochozoa, Mollusca) are economically important. Mass rearing of arthropods and molluscs causes problems with pathogens in aquatic ecosystems that are exploited by humans. Remarkably, species of corals (Cnidaria) living in non-exploited ecosystems also suffer from devastating infectious diseases that display intriguing similarities with those affecting farmed animals. Infectious diseases affecting wild and farmed animals that are present in marine environments are predicted to increase in the future. This paper summarizes the role of the main pathogens and their interaction with host immunity, with a specific focus on antimicrobial peptides (AMPs) and pathogen resistance against AMPs. We provide a detailed review of penaeid shrimp AMPs and their role at the interface between the host and its resident/pathogenic microbiota. We also briefly describe the relevance of marine invertebrate AMPs in an applied context.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Effect of temperature, food availability, and estradiol injection on gametogenesis and gender in the pearl oyster Pinctada margaritifera.

The black-lip pearl oyster Pinctada margaritifera is a protandrous hermaphrodite species. Its economic value has led to the development of controlled hatchery reproduction techniques, although many aspects remain to be optimized. In order to understand reproductive mechanisms and their controlling factors, two independent experiments were designed to test hypotheses of gametogenesis and sex ratio control by environmental and hormonal factors. In one, pearl oysters were exposed under controlled conditions at different combinations of temperature (24 and 28°C) and food level (10,000 and 40,000 cells mL(-1) ); whereas in the other, pearl oysters were conditioned under natural conditions into the lagoon and subjected to successive 17ß-estradiol injections (100 µg per injection). Gametogenesis and sex ratio were assessed by histology for each treatment. In parallel, mRNA expressions of nine marker genes of the sexual pathway (pmarg-foxl2, pmarg-c43476, pmarg-c45042, pmarg-c19309, pmarg-c54338, pmarg-vit6, pmarg-zglp1, pmarg-dmrt, and pmarg-fem1-like) were investigated. Maximum maturation was observed in the treatment combining the highest temperature (28°C) and the highest microalgae concentration (40,000 cells mL(-1) ), where the female sex tended to be maintained. Injection of 17ß-estradiol induced a significant increase of undetermined stage proportion 2 weeks after the final injection. These results suggest that gametogenesis and gender in adult pearl oysters can be controlled by environmental factors and estrogens. While there were no significant effects on relative gene expression, the 3-gene-pair expression ratio model of the sexual pathway of P. margaritifera, suggest a probable dominance of genetic sex determinism without excluding a mixed sex determination mode (genetic + environmental). J. Exp. Zool. 325A:13-24, 2016. © 2015 Wiley Periodicals, Inc.

Yes, it turns: experimental evidence of pearl rotation during its formation.

Cultured pearls are human creations formed by inserting a nucleus and a small piece of mantle tissue into a living shelled mollusc, usually a pearl oyster. Although many pearl observations intuitively suggest a possible rotation of the nucleated pearl inside the oyster, no experimental demonstration of such a movement has ever been done. This can be explained by the difficulty of observation of such a phenomenon in the tissues of a living animal. To investigate this question of pearl rotation, a magnetometer system was specifically engineered to register magnetic field variations with magnetic sensors from movements of a magnetic nucleus inserted in the pearl oyster. We demonstrated that a continuous movement of the nucleus inside the oyster starts after a minimum of 40 days post-grafting and continues until the pearl harvest. We measured a mean angular speed of 1.27° min(-1) calculated for four different oysters. Rotation variability was observed among oysters and may be correlated to pearl shape and defects. Nature's ability to generate so amazingly complex structures like a pearl has delivered one of its secrets.

Identification of genes associated with shell color in the black-lipped pearl oyster, Pinctada margaritifera.

BACKGROUND: Color polymorphism in the nacre of pteriomorphian bivalves is of great interest for the pearl culture industry. The nacreous layer of the Polynesian black-lipped pearl oyster Pinctada margaritifera exhibits a large array of color variation among individuals including reflections of blue, green, yellow and pink in all possible gradients. Although the heritability of nacre color variation patterns has been demonstrated by experimental crossing, little is known about the genes involved in these patterns. In this study, we identify a set of genes differentially expressed among extreme color phenotypes of P. margaritifera using a suppressive and subtractive hybridization (SSH) method comparing black phenotypes with full and half albino individuals. RESULTS: Out of the 358 and 346 expressed sequence tags (ESTs) obtained by conducting two SSH libraries respectively, the expression patterns of 37 genes were tested with a real-time quantitative PCR (RT-qPCR) approach by pooling five individuals of each phenotype. The expression of 11 genes was subsequently estimated for each individual in order to detect inter-individual variation. Our results suggest that the color of the nacre is partially under the influence of genes involved in the biomineralization of the calcitic layer. A few genes involved in the formation of the aragonite tablets of the nacre layer and in the biosynthesis chain of melanin also showed differential expression patterns. Finally, high variability in gene expression levels were observed within the black phenotypes. CONCLUSIONS: Our results revealed that three main genetic processes were involved in color polymorphisms: the biomineralization of the nacreous and calcitic layers and the synthesis of pigments such as melanin, suggesting that color polymorphism takes place at different levels in the shell structure. The high variability of gene expression found within black phenotypes suggests that the present work should serve as a basis for future studies exploring more thoroughly the expression patterns of candidate genes within black phenotypes with different dominant iridescent colors.