High polymorphism in big defensin gene expression reveals presence-absence gene variability (PAV) in the oyster Crassostrea gigas

Cross-talk and mutual shaping between the immune system and the microbiota during an oyster's life

The Pacific oyster Crassostrea gigas lives in microbe-rich marine coastal systems subjected to rapid environmental changes. It harbours a diversified and fluctuating microbiota that cohabits with immune cells expressing a diversified immune gene repertoire. In the early stages of oyster development, just after fertilization, the microbiota plays a key role in educating the immune system. Exposure to a rich microbial environment at the larval stage leads to an increase in immune competence throughout the life of the oyster, conferring a better protection against pathogenic infections at later juvenile/adult stages. This beneficial effect, which is intergenerational, is associated with epigenetic remodelling. At juvenile stages, the educated immune system participates in the control of the homeostasis. In particular, the microbiota is fine-tuned by oyster antimicrobial peptides acting through specific and synergistic effects. However, this balance is fragile, as illustrated by the Pacific Oyster Mortality Syndrome, a disease causing mass mortalities in oysters worldwide. In this disease, the weakening of oyster immune defences by OsHV-1 µVar virus induces a dysbiosis leading to fatal sepsis. This review illustrates the continuous interaction between the highly diversified oyster immune system and its dynamic microbiota throughout its life, and the importance of this cross-talk for oyster health. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Gene presence/absence variation in Mytilus galloprovincialis and its implications in gene expression and adaptation

Presence/absence variation (PAV) is a well-known phenomenon in prokaryotes that was described for the first time in bivalves in 2020 in Mytilus galloprovincialis. The objective of the present study was to further our understanding of the PAV phenomenon in mussel biology. The distribution of PAV was studied in a mussel chromosome-level genome assembly, revealing a widespread distribution but with hotspots of dispensability. Special attention was given to the effect of PAV in gene expression, since dispensable genes were found to be inherently subject to distortions due to their sparse distribution among individuals. Furthermore, the high expression and strong tissue specificity of some dispensable genes, such as myticins, strongly supported their biological relevance. The significant differences in the repertoire of dispensable genes associated with two geographically distinct populations suggest that PAV is involved in local adaptation. Overall, the PAV phenomenon would provide a key selective advantage at the population level.

Intraspecific de novo gene birth revealed by presence-absence variant genes in Caenorhabditis elegans

Genes embed their evolutionary history in the form of various alleles. Presence-absence variants (PAVs) are extreme cases of such alleles, where a gene present in one haplotype does not exist in another. Because PAVs may result from either birth or death of a gene, PAV genes and their alternative alleles, if available, can represent a basis for rapid intraspecific gene evolution. Using long-read sequencing technologies, this study traced the possible evolution of PAV genes in the PD1074 and CB4856 C. elegans strains as well as their alternative alleles in 14 other wild strains. We updated the CB4856 genome by filling 18 gaps and identified 46 genes and 7,460 isoforms from both strains not annotated previously. We verified 328 PAV genes, out of which 46 were C. elegans-specific. Among these possible newly born genes, 12 had alternative alleles in other wild strains; in particular, the alternative alleles of three genes showed signatures of active transposons. Alternative alleles of three other genes showed another type of signature reflected in accumulation of small insertions or deletions. Research on gene evolution using both species-specific PAV genes and their alternative alleles may provide new insights into the process of gene evolution.

Identification, characterization, and antimicrobial activity of a novel big defensin discovered in a commercial bivalve mollusc, Tegillarca granosa

Molluscs, the second largest animal phylum on earth, primarily rely on cellular and humoral immune responses to fight against pathogen infection. Although antimicrobial peptides (AMPs) such as big defensin play crucial roles in the humoral immune response, it remains largely unknown in the ecological and economic important blood clam (Tegillarca granosa). In this study, a novel big defensin gene (TgBD) was identified in T. granosa through transcripts and whole genome searching. Bioinformatic analyses were conducted to explore the molecular characteristics of TgBD, and comparisons of TgBD with those reported in other molluscs were performed by multiple alignments and phylogenetic analysis. In addition, the expression patterns of TgBD in various tissues and upon bacterial challenge were investigated while the antimicrobial activity of synthetic N-terminal domain of TgBD was confirmed in vitro by radial diffusion experiment. Results obtained showed TgBD had an open reading frame (ORF) of 369 bp, encoding a prepropeptide containing a signal peptide and a propeptide. Similar to big defensins reported in other species, TgBD consists of a hydrophobic N-terminal domain containing β1-α1-α2-β2 folds and a cysteine-rich cationic C-terminal domain with three disulfide bonds between C1-C5, C2-C4, and C3-C6. Phylogenetic analysis showed that TgBD shared 76.80% similarity to its close relative ark shell (Scapharca broughtoni). In addition, TgBD expression was observed in all tissues investigated under normal conditions and was significantly induced by injection of Vibrio parahaemolyticus. Furthermore, synthetic N-terminal peptide of TgBD exhibited strong antimicrobial activity against Gram-positive bacteria tested. Our results indicated that TgBD is a constitutive and inducible acute phase AMP, which provides a universal and prompt protection for T. granosa.

Transcriptome analysis reveals potential key immune genes of Hong Kong oyster (Crassostrea hongkongensis) against Vibrio parahaemolyticus infection

Hong Kong oyster (Crassostrea hongkongensis) is one of the main species of economic shellfish cultivated in the coastal areas of southern China. The cultivation of this shellfish may be adversely impacted by Vibrio parahaemolyticus, a harmful pathogenic bacterium for many mariculture species, as it usually exists on the surface of Hong Kong oysters. Although previous studies have discovered that oysters rely on non-specific immune system to fight pathogen invasion, the genes corresponding to the complex immune system against Vibrio is still not fully elucidated. Therefore, we conducted a transcriptome analysis on the gill from Hong Kong oysters at two time points (i.e., 12 h and 24 h after V. parahaemolyticus or PBS challenge) to identify potential immune genes against V. parahaemolyticus infection. A total of 61779 unigenes with the average length of 1221 bp were obtained, and the annotation information of 39917 unigenes were obtained from Nr, SwissProt, KEGG and COG/KOG. After a pairwise comparison between V. parahaemolyticus or PBS challenge at the two time points, three groups of differentially expressed genes induced by V. parahaemolyticus were captured and analyzed. GO and KEGG analyses showed that multiple immune-related genes played an important role in pathogen infection, including HSP70, PCDP3 and TLR4. Furthermore, genes annotation indicated that LITAF, TNFSF10, Duox2 and big defensin family are also involved in immune regulation. Our study provides a reference for further exploration the molecular mechanism that defenses the pathogen infection regarding the identified immune-related genes in Hong Kong oysters.

Detecting Structural Variants and Associated Gene Presence-Absence Variation Phenomena in the Genomes of Marine Organisms

As complete genomes become easier to attain, even from previously difficult-to-sequence species, and as genomic resequencing becomes more routine, it is becoming obvious that genomic structural variation is more widespread than originally thought and plays an important role in maintaining genetic variation in populations. Structural variants (SVs) and associated gene presence-absence variation (PAV) can be important players in local adaptation, allowing the maintenance of genetic variation and taking part in other evolutionarily relevant phenomena. While recent studies have highlighted the importance of structural variation in Mollusca, the prevalence of this phenomenon in the broader context of marine organisms remains to be fully investigated.Here, we describe a straightforward and broadly applicable method for the identification of SVs in fully assembled diploid genomes, leveraging the same reads used for assembly. We also explain a gene PAV analysis protocol, which could be broadly applied to any species with a fully sequenced reference genome available. Although the strength of these approaches have been tested and proven in marine invertebrates, which tend to have high levels of heterozygosity, possibly due to their lifestyle traits, they are also applicable to other species across the tree of life, providing a ready means to begin investigations into this potentially widespread phenomena.

Investigating the Transcriptomic and Expression Presence-Absence Variation Exist in Japanese Eel (Anguilla japonica), a Primitive Teleost

The pan-genome was defined as the complete gene set across strains, and it is built upon genes displaying presence-absence variations (PAVs); the pan-transcriptome is defined by recalling the pan-genome. Indeed, a PAV is reflected from the expression presence-absence variation (ePAV). In this study, treated with androgen, eels, which are a primitive fish from the basal lineage of Teleost, with different ovarian developments were chosen and submitted to RAN-sequencing. Transcriptomes were the assembly against eel genome scaffolds; a pair was the unit (the same eel before and after treatment) to analyze DEGs (differentially expressed genes); the core, unique, or accessory genes were identified, and the list of DEGs was analyzed to investigate ePAV. The results suggest that there was ePAV in Japanese eel, and the ePAV of eel was analyzed by pathway enrichment. These results signify the importance of genetic differential expression on the variations of phenotypes by androgen, and a transcriptomic approach appears to enable extracting multiple layers of genomic data.

Single individual structural variant detection uncovers widespread hemizygosity in molluscs

The advent of complete genomic sequencing has opened a window into genomic phenomena obscured by fragmented assemblies. A good example of these is the existence of hemizygous regions of autosomal chromosomes, which can result in marked differences in gene content between individuals within species. While these hemizygous regions, and presence/absence variation of genes that can result, are well known in plants, firm evidence has only recently emerged for their existence in metazoans. Here, we use recently published, complete genomes from wild-caught molluscs to investigate the prevalence of hemizygosity across a well-known and ecologically important clade. We show that hemizygous regions are widespread in mollusc genomes, not clustered in individual chromosomes, and often contain genes linked to transposition, DNA repair and stress response. With targeted investigations of HSP70-12 and C1qDC, we also show how individual gene families are distributed within pan-genomes. This work suggests that extensive pan-genomes are widespread across the conchiferan Mollusca, and represent useful tools for genomic evolution, allowing the maintenance of additional genetic diversity within the population. As genomic sequencing and re-sequencing becomes more routine, the prevalence of hemizygosity, and its impact on selection and adaptation, are key targets for research across the tree of life. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Extensive genome-wide duplications in the eastern oyster (Crassostrea virginica)

Genomic structural variation is an important source of genetic and phenotypic diversity, playing a critical role in evolution. The recent availability of a high-quality reference genome for the eastern oyster, Crassostrea virginica, and whole-genome sequence data of samples from across the species range in the USA, provides an opportunity to explore structural variation across the genome of this species. Our analysis shows significantly greater individual-level duplications of regions across the genome than that of most model vertebrate species. Duplications are widespread across all ten chromosomes with variation in frequency per chromosome. The eastern oyster shows a large interindividual variation in duplications as well as particular chromosomal regions with a higher density of duplications. A high percentage of duplications seen in C. virginica lie completely within genes and exons, suggesting the potential for impacts on gene function. These results support the hypothesis that structural changes may play a significant role in standing genetic variation in C. virginica, and potentially have a role in their adaptive and evolutionary success. Altogether, these results suggest that copy number variation plays an important role in the genomic variation of C. virginica. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Massive gene presence-absence variation shapes an open pan-genome in the Mediterranean mussel

BACKGROUND

The Mediterranean mussel Mytilus galloprovincialis is an ecologically and economically relevant edible marine bivalve, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other bivalves. Although these traits have been recently linked with the maintenance of a high genetic variation within natural populations, the factors underlying the evolutionary success of this species remain unclear.

RESULTS

Here, after the assembly of a 1.28-Gb reference genome and the resequencing of 14 individuals from two independent populations, we reveal a complex pan-genomic architecture in M. galloprovincialis, with a core set of 45,000 genes plus a strikingly high number of dispensable genes (20,000) subject to presence-absence variation, which may be entirely missing in several individuals. We show that dispensable genes are associated with hemizygous genomic regions affected by structural variants, which overall account for nearly 580 Mb of DNA sequence not included in the reference genome assembly. As such, this is the first study to report the widespread occurrence of gene presence-absence variation at a whole-genome scale in the animal kingdom.

CONCLUSIONS

Dispensable genes usually belong to young and recently expanded gene families enriched in survival functions, which might be the key to explain the resilience and invasiveness of this species. This unique pan-genome architecture is characterized by dispensable genes in accessory genomic regions that exceed by orders of magnitude those observed in other metazoans, including humans, and closely mirror the open pan-genomes found in prokaryotes and in a few non-metazoan eukaryotes.

Functional Insights From the Evolutionary Diversification of Big Defensins

Big defensins are antimicrobial polypeptides believed to be the ancestors of β-defensins, the most evolutionary conserved family of host defense peptides (HDPs) in vertebrates. Nevertheless, big defensins underwent several independent gene loss events during animal evolution, being only retained in a limited number of phylogenetically distant invertebrates. Here, we explore the evolutionary history of this fascinating HDP family and investigate its patchy distribution in extant metazoans. We highlight the presence of big defensins in various classes of lophotrochozoans, as well as in a few arthropods and basal chordates (amphioxus), mostly adapted to life in marine environments. Bivalve mollusks often display an expanded repertoire of big defensin sequences, which appear to be the product of independent lineage-specific gene tandem duplications, followed by a rapid molecular diversification of newly acquired gene copies. This ongoing evolutionary process could underpin the simultaneous presence of canonical big defensins and non-canonical (β-defensin-like) sequences in some species. The big defensin genes of mussels and oysters, two species target of in-depth studies, are subjected to gene presence/absence variation (PAV), i.e., they can be present or absent in the genomes of different individuals. Moreover, big defensins follow different patterns of gene expression within a given species and respond differently to microbial challenges, suggesting functional divergence. Consistently, current structural data show that big defensin sequence diversity affects the 3D structure and biophysical properties of these polypeptides. We discuss here the role of the N-terminal hydrophobic domain, lost during evolution toward β-defensins, in the big defensin stability to high salt concentrations and its mechanism of action. Finally, we discuss the potential of big defensins as markers for animal health and for the nature-based design of novel therapeutics active at high salt concentrations.

Presence-absence polymorphisms of single-copy genes in the stony coral Acropora digitifera

Background

Despite the importance of characterizing genetic variation among coral individuals for understanding phenotypic variation, the correlation between coral genomic diversity and phenotypic expression is still poorly understood.

Results

In this study, we detected a high frequency of genes showing presence–absence polymorphisms (PAPs) for single-copy genes in Acropora digitifera. Among 10,455 single-copy genes, 516 (5%) exhibited PAPs, including 32 transposable element (TE)-related genes. Five hundred sixteen genes exhibited a homozygous absence in one (102) or more than one (414) individuals (n = 33), indicating that most of the absent alleles were not rare variants. Among genes showing PAPs (PAP genes), roughly half were expressed in adults and/or larvae, and the PAP status was associated with differential expression among individuals. Although 85% of PAP genes were uncharacterized or had ambiguous annotations, 70% of these genes were specifically distributed in cnidarian lineages in eumetazoa, suggesting that these genes have functional roles related to traits related to cnidarians or the family Acroporidae or the genus Acropora. Indeed, four of these genes encoded toxins that are usually components of venom in cnidarian-specific cnidocytes. At least 17% of A. digitifera PAP genes were also PAPs in A. tenuis, the basal lineage in the genus Acropora, indicating that PAPs were shared among species in Acropora.

Conclusions

Expression differences caused by a high frequency of PAP genes may be a novel genomic feature in the genus Acropora; these findings will contribute to improve our understanding of correlation between genetic and phenotypic variation in corals.

The Ancestral N-Terminal Domain of Big Defensins Drives Bacterially Triggered Assembly into Antimicrobial Nanonets

β-Defensins are host defense peptides controlling infections in species ranging from humans to invertebrates. However, the antimicrobial activity of most human β-defensins is impaired at physiological salt concentrations. We explored the properties of big defensins, the β-defensin ancestors, which have been conserved in a number of marine organisms, mainly mollusks. By focusing on a big defensin from oyster (Cg-BigDef1), we showed that the N-terminal domain lost during evolution toward β-defensins confers bactericidal activity to Cg-BigDef1, even at high salt concentrations. Cg-BigDef1 killed multidrug-resistant human clinical isolates of Staphylococcus aureus. Moreover, the ancestral N-terminal domain drove the assembly of the big defensin into nanonets in which bacteria are entrapped and killed. This discovery may explain why the ancestral N-terminal domain has been maintained in diverse marine phyla and creates a new path of discovery to design β-defensin derivatives active at physiological and high salt concentrations.

Big defensins, ancestors of β-defensins, are composed of a β-defensin-like C-terminal domain and a globular hydrophobic ancestral N-terminal domain. This unique structure is found in a limited number of phylogenetically distant species, including mollusks, ancestral chelicerates, and early-branching cephalochordates, mostly living in marine environments. One puzzling evolutionary issue concerns the advantage for these species of having maintained a hydrophobic domain lost during evolution toward β-defensins. Using native ligation chemistry, we produced the oyster Crassostrea gigas BigDef1 (Cg-BigDef1) and its separate domains. Cg-BigDef1 showed salt-stable and broad-range bactericidal activity, including against multidrug-resistant human clinical isolates of Staphylococcus aureus. We found that the ancestral N-terminal domain confers salt-stable antimicrobial activity to the β-defensin-like domain, which is otherwise inactive. Moreover, upon contact with bacteria, the N-terminal domain drives Cg-BigDef1 assembly into nanonets that entrap and kill bacteria. We speculate that the hydrophobic N-terminal domain of big defensins has been retained in marine phyla to confer salt-stable interactions with bacterial membranes in environments where electrostatic interactions are impaired. Those remarkable properties open the way to future drug developments when physiological salt concentrations inhibit the antimicrobial activity of vertebrate β-defensins.

Simulated Marine Heat Wave Alters Abundance and Structure of Vibrio Populations Associated with the Pacific Oyster Resulting in a Mass Mortality Event

Marine heat waves are predicted to become more frequent and intense due to anthropogenically induced climate change, which will impact global production of seafood. Links between rising seawater temperature and disease have been documented for many aquaculture species, including the Pacific oyster Crassostrea gigas. The oyster harbours a diverse microbial community that may act as a source of opportunistic pathogens during temperature stress. We rapidly raised the seawater temperature from 20 °C to 25 °C resulting in an oyster mortality rate of 77.4%. Under the same temperature conditions and with the addition of antibiotics, the mortality rate was only 4.3%, strongly indicating a role for bacteria in temperature-induced mortality. 16S rRNA amplicon sequencing revealed a change in the oyster microbiome when the temperature was increased to 25 °C, with a notable increase in the proportion of Vibrio sequences. This pattern was confirmed by qPCR, which revealed heat stress increased the abundance of Vibrio harveyi and Vibrio fortis by 324-fold and 10-fold, respectively. Our findings indicate that heat stress-induced mortality of C. gigas coincides with an increase in the abundance of putative bacterial pathogens in the oyster microbiome and highlights the negative consequences of marine heat waves on food production from aquaculture.

High individual variability in the transcriptomic response of Mediterranean mussels to Vibrio reveals the involvement of myticins in tissue injury

Mediterranean mussels (Mytilus galloprovincialis) are sessile filter feeders that live in close contact with numerous marine microorganisms. As all invertebrates, they lack an adaptive immune response and how these animals are able to respond to a bacterial infection and discriminate it from their normal microbiome is difficult to understand. In this work, we conducted Illumina sequencing of the transcriptome of individual mussels before and after being infected with Vibrio splendidus. The control mussels were injected with filtered seawater. We demonstrate that a great variability exists among individual transcriptomes and that each animal showed an exclusive repertoire of genes not shared with other individuals. The regulated genes in both the control and infected mussels were also analyzed and, unexpectedly, the sampling before the injection was considered a stress stimulus strong enough to trigger and modulate the response in hemocytes, promoting cell migration and proliferation. We found a clear response against the injection of filtered seawater, suggesting a reaction against a tissue injury in which the myticins, the most expressed antimicrobial peptides in mussel, appeared significantly up regulated. Functional experiments with flow cytometry confirmed the transcriptomic results since a significant alteration of hemocyte structures and a decrease in the number of hemocytes positive for myticin C were found only after a Vibrio infection and not observed when mussels were bled before, generating a tissue injury. Therefore, we report the involvement of myticins in the response to a danger signal such as a simple injection in the adductor muscle.

Relationship Among Intron Length, Gene Expression, and Nucleotide Diversity in the Pacific Oyster Crassostrea gigas

Crassostrea gigas is a model mollusk, but its genetic features have not been studied comprehensively. In this study, we used whole-genome resequencing data to identify and characterize nucleotide diversity and population recombination rate in a diverse collection of 21 C. gigas samples. Our analyses revealed that C. gigas harbors both extremely high genetic diversity and recombination rates across the whole genome as compared with those of the other taxa. The noncoding regions, introns, intergenic spacers, and untranslated regions (UTRs) showed a lower level diversity than the synonymous sites. The larger introns tended to have lower diversity. Moreover, we found a negative association of the non-synonymous diversity with gene expression, which suggested that purifying selection played an important role in shaping genetic diversity. The nucleotide diversity at the 100- and 50-kb levels was positively correlated with population recombination rates, which was expected if the diversity was shaped by purifying selection or hitchhiking of advantageous mutants. Our work gives a general picture of the oyster's polymorphism pattern and its association with recombination rates.

Genetic and evolutionary patterns of innate immune genes in the Pacific oyster Crassostrea gigas

The invertebrate innate immune system functions in immune defence and the stress response. However, knowledge of the genetic and evolutionary patterns of innate immune genes in Mollusca is limited, especially for oysters. Such information would help clarify how oysters adapt to pathogen-rich environments. Here, we characterized the genetic and evolutionary patterns of the innate immune genes in Crassostrea gigas, using population diversity analysis and evolution rates comparison. Innate immune genes have higher median nucleotide diversity than non-immune genes. Nucleotide diversity varied with functional regions and different immune-related gene families. Evolutionary analysis of two Crassostrea species showed that the innate immune genes are less conserved and have higher rates of evolution in C. gigas. We also noted a positive association between nucleotide diversity and selective pressures for genes having orthologues. Our findings will help determine the evolutionary patterns of innate immune genes and the association of these genes with mollusc immunity.

The association between DNA methylation and exon expression in the Pacific oyster Crassostrea gigas

Background

DNA methylation is one of the most important epigenetic modifications of eukaryotic genomes and is believed to play integral roles in diverse biological processes. Although DNA methylation has been well studied in mammals, data are limited in invertebrates, particularly Mollusca. The Pacific oyster Crassostrea gigas is an emerging genetic model for functional analysis of DNA methylation in Mollusca. Recent studies have shown that there is a positive association between methylation status and gene expression in C. gigas; however, whether this association exists at the exon level remains to be determined.

Results

In this study, we characterized the genome-wide methylation pattern across two different tissues of C. gigas and found that methylated genes are expressed in more tissues and development stages than unmethylated genes. Furthermore, we found that different types of exons had different methylation levels, with the lowest methylation levels in the first exons, followed by the last exons, and the internal exons. We found that the exons included in the gene transcript contained significantly higher DNA methylation levels than skipped exons. We observed that the DNA methylation levels increased slowly after the start sites and end sites of exons seperately, and then decreased quickly towards the middle sites of exons. We also found that methylated exons were significantly longer than unmethylated exons.

Conclusion

This study constitutes the first genome-wide analysis to show an association between exon-level DNA methylation and mRNA expression in the oyster. Our findings suggest that exon-level DNA methylation may play a role in the construction of alternative splicing by positively influencing exon inclusion during transcription.

Molecular characterization and protein localization of the antimicrobial peptide big defensin from the scallop Argopecten purpuratus after Vibrio splendidus challenge

Big defensins are antimicrobial peptides (AMPs) that are proposed as important effectors of the immune response in mollusks, chelicerates and chordates. At present, only two members of the big defensin family have been identified in scallop. In the present work, a cDNA sequence encoding a new big defensin homologue was characterized from the scallop Argopecten purpuratus, namely ApBD1. ApBD1 cDNA sequence comprised 585 nucleotides, with an open reading frame of 375 bp and 5'- and 3'-UTRs of 41 and 167 bp, respectively. The deduced protein sequence contains 124 amino acids with a molecular weight of 13.5 kDa, showing characteristic motifs of the big defensin family and presenting 76% identity with the big defensin from the scallop A. irradians. Phylogenetic analysis revealed that ApBD1 is included into the cluster of big defensins from mollusks. Tissue-specific transcript expression analysis by RT-qPCR showed that ApBD1 was present in all tissues tested from non-immune challenged scallops but it was most strongly expressed in the mantle. The transcript levels of ApBD1 were significantly up-regulated in gills at 24 and 48 h post-injection with the heat-attenuated bacteria Vibrio splendidus. Additionally, immunofluorescence analysis using a polyclonal anti-ApBD1 antibody showed that this protein was abundantly located in epithelial linings of gills and mantle; and also in digestive gland showing ApBD1-infiltrating hemocytes from immune challenged scallops. This is the first time that a big defensin is detected and located at the protein level in a mollusk. These results suggest an important role of ApBD1 in the mucosal immune response of A. purpuratus.

Myticalins: A Novel Multigenic Family of Linear, Cationic Antimicrobial Peptides from Marine Mussels (Mytilus spp.)

The application of high-throughput sequencing technologies to non-model organisms has brought new opportunities for the identification of bioactive peptides from genomes and transcriptomes. From this point of view, marine invertebrates represent a potentially rich, yet largely unexplored resource for de novo discovery due to their adaptation to diverse challenging habitats. Bioinformatics analyses of available genomic and transcriptomic data allowed us to identify myticalins, a novel family of antimicrobial peptides (AMPs) from the mussel Mytilus galloprovincialis, and a similar family of AMPs from Modiolus spp., named modiocalins. Their coding sequence encompasses two conserved N-terminal (signal peptide) and C-terminal (propeptide) regions and a hypervariable central cationic region corresponding to the mature peptide. Myticalins are taxonomically restricted to Mytiloida and they can be classified into four subfamilies. These AMPs are subject to considerable interindividual sequence variability and possibly to presence/absence variation. Functional assays performed on selected members of this family indicate a remarkable tissue-specific expression (in gills) and broad spectrum of activity against both Gram-positive and Gram-negative bacteria. Overall, we present the first linear AMPs ever described in marine mussels and confirm the great potential of bioinformatics tools for the de novo discovery of bioactive peptides in non-model organisms.

Myticalins: A Novel Multigenic Family of Linear, Cationic Antimicrobial Peptides from Marine Mussels (Mytilus spp.)

The application of high-throughput sequencing technologies to non-model organisms has brought new opportunities for the identification of bioactive peptides from genomes and transcriptomes. From this point of view, marine invertebrates represent a potentially rich, yet largely unexplored resource for de novo discovery due to their adaptation to diverse challenging habitats. Bioinformatics analyses of available genomic and transcriptomic data allowed us to identify myticalins, a novel family of antimicrobial peptides (AMPs) from the mussel Mytilus galloprovincialis, and a similar family of AMPs from Modiolus spp., named modiocalins. Their coding sequence encompasses two conserved N-terminal (signal peptide) and C-terminal (propeptide) regions and a hypervariable central cationic region corresponding to the mature peptide. Myticalins are taxonomically restricted to Mytiloida and they can be classified into four subfamilies. These AMPs are subject to considerable interindividual sequence variability and possibly to presence/absence variation. Functional assays performed on selected members of this family indicate a remarkable tissue-specific expression (in gills) and broad spectrum of activity against both Gram-positive and Gram-negative bacteria. Overall, we present the first linear AMPs ever described in marine mussels and confirm the great potential of bioinformatics tools for the de novo discovery of bioactive peptides in non-model organisms.

Parental diuron-exposure alters offspring transcriptome and fitness in Pacific oyster Crassostrea gigas

One of the primary challenges in ecotoxicology is to contribute to the assessment of the ecological status of ecosystems. In this study, we used Pacific oyster Crassostrea gigas to explore the effects of a parental exposure to diuron, a herbicide frequently detected in marine coastal environments. The present toxicogenomic study provides evidence that exposure of oyster genitors to diuron during gametogenesis results in changes in offspring, namely, transcriptomic profile alterations, increased global DNA methylation levels and reduced growth and survival within the first year of life. Importantly, we highlighted the limitations to identify particular genes or gene expression signatures that could serve as biomarkers for parental herbicide-exposure and further for multigenerational and transgenerational effects of specific chemical stressors. By analyzing samples from two independent experiments, we demonstrated that, due to complex confounding effects with both tested solvent vehicles, diuron non-specifically affected the offspring transcriptome. These original results question the potential development of predictive genomic tools for detecting specific indirect impacts of contaminants in environmental risk assessments. However, our results indicate that chronic environmental exposure to diuron over several generations may have significant long term impacts on oyster populations with adverse health outcomes.

Effects of a parental exposure to diuron on Pacific oyster spat methylome

Environmental epigenetic is an emerging field that studies the cause-effect relationship between environmental factors and heritable trait via an alteration in epigenetic marks. This field has received much attentions since the impact of environmental factors on different epigenetic marks have been shown to be associated with a broad range of phenotypic disorders in natural ecosystems. Chemical pollutants have been shown to affect immediate epigenetic information carriers of several aquatic species but the heritability of the chromatin marks and the consequences for long term adaptation remain open questions. In this work, we investigated the impact of the diuron herbicide on the DNA methylation pattern of spat from exposed Crassotrea gigas genitors. This oyster is one of the most important mollusk species produced worldwide and a key coastal economic resource in France. The whole genome bisulfite sequencing (WGBS, BS-Seq) was applied to obtain a methylome at single nucleotide resolution on DNA extracted from spat issued from diuron exposed genitors comparatively to control spat. We showed that the parental diuron exposure has an impact on the DNA methylation pattern of its progeny. Most of the differentially methylated regions occurred within coding sequences and we showed that this change in methylation level correlates with RNA level only in a very small group of genes. Although the DNA methylation profile is variable between individuals, we showed conserved DNA methylation patterns in response to parental diuron exposure. This relevant result opens perspectives for the setting of new markers based on epimutations as early indicators of marine pollutions.

Effect of 4-nonylphenol on the immune response of the Pacific oyster Crassostrea gigas following bacterial infection with Vibrio campbellii

The xenoestrogen 4-nonylphenol (NP) is a ubiquitous aquatic pollutant and has been shown to impair reproduction, development, growth and, more recently, immune function in marine invertebrates. We investigated the effects of short-term (7 d) exposure to low (2 μg l^-1) and high (100 μg l^-1) levels of NP on cellular and humoral elements of the innate immune response of Crassostrea gigas to a bacterial challenge. To this end, we measured 1) total hemocyte counts (THC), 2) relative transcript abundance of ten immune-related genes (defh1, defh2, bigdef1, bigdef2, bpi, lysozyme-1, galectin, C-type lectin 2, timp, and transglutaminase) in the hemocytes, gill and mantle, and 3) hemolymph plasma lysozyme activity, following experimental Vibrio campbellii infection. Both low and high levels of NP were found to repress a bacteria-induced increase in THC observed in the control oysters. While several genes were differentially expressed following bacterial introduction (bigdef2, bpi, lysozyme-1, timp, transglutaminase), only two genes (bpi in the hemocytes, transglutaminase in the mantle) exhibited a different bacteria-induced expression profile following NP exposure, relative to the control oysters. Independently of infection-status, exposure to NP also altered mRNA transcript abundance of several genes (bpi, galectin, C-type lectin 2) in naïve, saline-injected oysters. Finally, plasma lysozyme activity levels were significantly higher in low dose NP-treated oysters (both naïve and bacteria challenged) relative to control oysters. Combined, these results suggest that exposure to ecologically-relevant (low) and extreme (high) levels of NP can alter both cellular and humoral elements of the innate immune response in C. gigas, an aquaculture species of global economic importance.

The new insights into the oyster antimicrobial defense: Cellular, molecular and genetic view

Oysters are sessile filter feeders that live in close association with abundant and diverse communities of microorganisms that form the oyster microbiota. In such an association, cellular and molecular mechanisms have evolved to maintain oyster homeostasis upon stressful conditions including infection and changing environments. We give here cellular and molecular insights into the Crassostrea gigas antimicrobial defense system with focus on antimicrobial peptides and proteins (AMPs). This review highlights the central role of the hemocytes in the modulation and control of oyster antimicrobial response. As vehicles for AMPs and other antimicrobial effectors, including reactive oxygen species (ROS), and together with epithelia, hemocytes provide the oyster with local defense reactions instead of systemic humoral ones. These reactions are largely based on phagocytosis but also, as recently described, on the extracellular release of antimicrobial histones (ETosis) which is triggered by ROS. Thus, ROS can signal danger and activate cellular responses in the oyster. From the current literature, AMP production/release could serve similar functions. We provide also new lights on the oyster genetic background that underlies a great diversity of AMP sequences but also an extraordinary individual polymorphism of AMP gene expression. We discuss here how this polymorphism could generate new immune functions, new pathogen resistances or support individual adaptation to environmental stresses.