# Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture

Isabelle Boutet, Homère Alves Monteiro, Lyam Baudry, Takeshi Takeuchi, Eric Bonnivard, Bernard Billoud, Sarah Farhat, Ricardo Gonzales-Araya, Benoit Salaun, Ann Andersen, et al.

► To cite this version:

Isabelle Boutet, Homère Alves Monteiro, Lyam Baudry, Takeshi Takeuchi, Eric Bonnivard, et al. Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture. Evolutionary Applications, 2022, 15 (11), pp.1730-1748. 10.1111/eva.13462. hal-03896508

HAL Id: hal-03896508

https://hal.science/hal-03896508

Submitted on 15 Dec 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
DOI: 10.1111/eva.13462

ORIGINAL ARTICLE

Evolutionary Applications
Open Access

WILEY

# Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture

Isabelle Boutet<fnref n="1" /> | Homère J. Alves Monteiro<fnref n="2" /> | Lyam Baudry<fnref n="3" /> | Takeshi Takeuchi<fnref n="4" /> |

Eric Bonnivard<fnref n="1" /> | Bernard Billoud<fnref n="5" /> | Sarah Farhat<fnref n="6" /> | Ricardo Gonzales-Araya<fnref n="7" /> ID |

Benoit Salaun<fnref n="7" /> | Ann C. Andersen<fnref n="1" /> | Jean-Yves Toullec<fnref n="1" /> | François H. Lallier<fnref n="1" /> ID |

Jean-François Flot<fnref n="8" /> ID | Nadège Guiglielmoni<fnref n="8" /> | Ximing Guo<fnref n="9" /> | Cui Li<fnref n="10" /> |

Bassem Allam<fnref n="6" /> | Emmanuelle Pales-Espinosa<fnref n="6" /> | Jakob Hemmer-Hansen<fnref n="2" /> ID |

Pierrick Moreau<fnref n="3" /> | Martial Marbouty<fnref n="3" /> | Romain Koszul<fnref n="3" /> ID | Arnaud Tanguy<fnref n="1" /> ID

Correspondence

Arnaud Tanguy, Sorbonne Université,
CNRS, UMR 7144, Station Biologique de
Roscoff, Place Georges Teissier, CS90074,
29688 Roscoff, France.

Email: atanguy@sb-rocoff.fr

Abstract

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites *Bonamia ostreae* and *Marteilia refringens*. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the

This is an open access article under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2022 The Authors. *Evolutionary Applications* published by John Wiley & Sons Ltd.

<footnote n="1">¹Sorbonne Université, CNRS, UMR 7144, Station Biologique de Roscoff, Roscoff, France</footnote>

<footnote n="2">²National Institute of Aquatic Resources, Technical University of Denmark, Silkeborg, Denmark</footnote>

<footnote n="3">³Institut Pasteur, Unité Régulation Spatiale des Génomes, CNRS, Paris, France</footnote>

<footnote n="4">⁴Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa, Japan</footnote>

<footnote n="5">⁵Sorbonne Université, CNRS, UMR 8227, Station Biologique de Roscoff, Roscoff, France</footnote>

<footnote n="6">⁶Marine Animal Disease Laboratory, School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, New York, USA</footnote>

<footnote n="7">⁷Centre Régional de la Conchyliculture Bretagne Nord, Morlaix, France</footnote>

<footnote n="8">⁸Evolutionary Biology and Ecology, Université Libre de Bruxelles, Brussels, Belgium</footnote>

<footnote n="9">⁹Haskin Shellfish Research Laboratory, Department of Marine and Coastal Sciences, Rutgers University, Port Norris, New Jersey, USA</footnote>

<footnote n="10">¹⁰Department of Marine Organism Taxonomy and Phylogeny, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China</footnote>
expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of *O. edulis* to *M. refringens*. This genome provides a reference for genomic studies on *O. edulis* to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

KEYWORDS

aquaculture, flat oyster, genome, Martelia, transposable elements

# 1 | BACKGROUND

The European flat oyster, *Ostrea edulis* (NCBI:txid37623), is the native oyster species in Europe and its distribution ranges from the Norwegian Sea in the north to Morocco in the south, and east through the Mediterranean Sea to the Black Sea. The flat oyster industry has been in decline since the 19th century mainly due to its habitat destruction, over-exploitation, irregular recruitment, and outbreaks of the parasites *Marteilia refringens* and *Bonamia ostreae* during the 1970s (Smaal et al., 2015). In France, the *O. edulis* production amounts to 1100 tons per year and is mainly located in Brittany. Because of its economic and ecological importance, *O. edulis* is subjected to extensive protection and restoration and has been identified as a priority species by the OSPAR Convention (OSPAR, 2009). Indeed, the flat oyster ecosystem is similar to calcareous biogenic reefs and contributes to substrate stabilization, water filtration, and habitat formation for other species (Beck et al., 2011). Unlike the *Crassostrea* oysters *Crassostrea gigas* or *Crassostrea virginica*, *O. edulis* is a larviparous species with asynchronous hermaphroditism and rhythmic consecutive sexuality. Molecular studies on the physiology of the flat oyster, including basic traits and adaptative response to abiotic parameters, are rare and mainly focused on the study of the response of *O. edulis* to its parasite *B. ostreae* in laboratory experiments (Cocci et al., 2020; Gervais et al., 2019; Ronza et al., 2018). Several programs have aimed to identify quantitative trait loci (QTLs) associated with *B. ostreae* resistance with the use of comparative transcriptomic approaches (Pardo et al., 2016; Ronza et al., 2018) and SNP genotyping (Vera et al., 2019).

Aquaculture is increasingly benefiting from the availability of genomic information, allowing for high-quality basic and applied genomic research to achieve various objectives such as improving aquaculture production and studying the effect of external factors on the physiology and genetics of the animal species concerned and their pathogens (Bernatchez et al., 2017; Hollenbeck & Johnston, 2018; Potts et al., 2021). While traditional selective breeding has improved disease resistance in some mollusk species when the heritability is high, resistance to some diseases may be also under the control of many small-effect genes where genomic selection is expected to be more effective (Guo, 2021). To partly resolve these limitations and provide genomic resources for aquaculture

breeding, genomes of many bivalves have been sequenced, including the oysters *C. gigas* (Peñaloza et al., 2021; Zhang et al., 2012), *Crassostrea hongkongensis* (Peng et al., 2020), *Crassostrea ariakensis* (Li et al., 2021), *C. virginica* (Modak et al., 2021), and *Saccostrea glomerata* (Powell et al., 2018), the pearl oyster *Pinctada fucata* (Du et al., 2017), the scallops *Patinopecten yessoensis* (Wang et al., 2017) and *Pecten maximus* (Kenny et al., 2020), the blood clam *Scapharca broughtonii* (Bai et al., 2019), the hard-shelled mussel *Mytilus coruscus* (Yang et al., 2021), and the hard clam *Mercenaria mercenaria* (Farhat et al., 2022; Song et al., 2021). Since 2018, we initiated a selective breeding program (European project PERLE 2) to characterize rusticity parameters, including survival, growth, and resistance to parasites in the flat oyster *O. edulis*. The main objectives of this project are to (1) produce hundreds of full-sib families tested in different environments (deep water and foreshore), (2) introduce the most performant families in natural oyster beds showing a decline in their population in order to enhance their dynamics and then contribute to their restoration, and (3) specifically select the most performing families on the foreshore in order to reintroduce this mode of culture for *O. edulis*. In the present study, an improved chromosome-level assembly of *O. edulis* genome was developed through a combination of high coverage Pacific Biosciences (PacBio) long-read sequencing, 10X chromium library sequencing, accurate Illumina short-read data and Hi-C library sequencing. We focused our analyses on the diversity of transposable elements through a comparison with other oyster genomes, highlighted some specificities in noncoding RNA diversity and provided the first RNAseq analysis on *M. refringens*-infected individuals revealing some original molecular responses.

# 2 | MATERIALS AND METHODS

## 2.1 | Tissue sampling, RNA isolation, and RNA-sequencing for genome annotation

Tissues from oyster organs (gill, mantle, digestive gland, adductor muscle, gonads, hemocytes and palp) were sampled from 12 individuals and frozen in liquid nitrogen. Total RNA was extracted using TRIzol Reagent (Invitrogen, Life Technologies) according to the manufacturer's instructions. The RNA quality was checked on agarose gel
and quantified using a NanoPhotometer spectrophotometer (Thermo Fisher Scientific). An equal amount of RNA of each tissue from four individuals was pool to generate three libraries per tissue. Libraries were constructed at the McGill Genome Center using NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina (New England Biolabs) following the manufacturer's recommendations and sequenced by Illumina HiSeq with 2x150 cycles (Illumina Inc.).

## 2.2 | RNA extraction in *M. refringens*-infected oyster

Fifty oysters (size: 5–8 cm) were collected at the "banc du Roz" (Bay of Brest, France). Hemolymph was retrieved from the adductor muscle using a syringe, and hemocytes were collected by centrifugation and immediately frozen in liquid nitrogen. On the same individual, a piece of the digestive gland was sampled for DNA extraction and parasite detection by PCR using two specific *M. refringens* primers (Le Roux et al., 1999), and the remaining tissues were frozen in liquid nitrogen. RNA was then extracted from hemocytes and the digestive gland of the six *M. refringens* positive individuals identified by PCR and on six negative individuals using TRIzol Reagent (Invitrogen). The RNA quality was checked on agarose gel and quantified using a NanoPhotometer spectrophotometer (Thermo Fisher Scientific). Two pools containing each an equal amount of RNA from three individuals showing a positive PCR amplification (six individuals in total) and two pools containing each an equal amount of RNA from three individuals showing no PCR amplification (six individuals in total) were generated for RNAseq library construction and sequencing using the same protocol as described above and sequenced by Illumina HiSeq with 2x150 cycles (Illumina Inc.).

## 2.3 | DNA library preparation and whole genome sequencing

### 2.3.1 | DNA extraction

Genomic DNA used for PacBio, Illumina shotgun, and 10× chromium libraries was extracted from fresh adductor muscle dissected off a single specimen of *O. edulis* (size = 18 cm) collected from the Bay of Morlaix (Brittany, France), using a standard phenol-chloroform-isoamyl alcohol (PCI 25:24:1) protocol and ethanol precipitation. DNA quality was checked by electrophoresis on 1% agarose gels and DNA concentration was measured using a Qubit® dsDNA HS Assay Kit in Qubit® 2.0 Fluorometer (Life Technologies).

### 2.3.2 | Shotgun library

Shotgun library was generated at the McGill Genome Center using the NxSeq® AmpFREE Low DNA Library Kit Library Preparation

Kit (Lucigen Corp.) according to the manufacturer's recommendations and sequenced by Illumina HiSeqX with 2x150 cycles (Illumina Inc.).

### 2.3.3 | Pacbio sequel libraries

The DNA library was prepared following the Pacific Biosciences 20kb Template Preparation using 7.5 µg of high molecular weight genomic DNA using the SMRTbell Template Prep Kit 1.0 reagents (Pacific Biosciences). The DNA library was size-selected on a BluePippin system (Sage Science Inc.) using a cutoff range of 12–50kb. The libraries were sequenced on a PacBio Sequel instrument at a loading concentration (on-plate) of 6 pM using the diffusion loading protocol, Sequel Sequencing Plate 2.1, SMRT cells 1M v2, and 10h movies.

### 2.3.4 | 10× chromium library

The gDNA was size-selected on a BluePippin system (Sage Science Inc.) using a cutoff range of 40–80kb. The 10× Chromium shotgun libraries were prepared following the Chromium Genome Reagent kits v2 User Guide RevB protocol, using the Chromium™ Genome Library & Gel Bead Kit v2, Chromium™ Genome Chip Kit, and Chromium™ i7 Multiplex Kit (10X Genomics Inc.) and sequenced by Illumina HiSeqX with 2x150 cycles (Illumina Inc.).

### 2.3.5 | Hi-C libraries

The Hi-C library construction protocol was adapted from Lieberman-Aiden et al. (2009) and Lazar-Stefanita et al. (2017). Hi-C library was made from a piece of adductor muscle from another individual from the same population and sequenced on a NextSeq 550 apparatus (2x75 bp, paired-end Illumina NextSeq). Contact maps were generated from reads using the hicstuff pipeline for processing generic 3C data, available at https://github.com/koszullab/hicstuff. The backend uses the bowtie2 (version 2.2.5) aligner run in paired-end mode (with the following options: --maxins 5 -very-sensitive-local). Alignments with a mapping quality lower than 30 were discarded. The output was in the form of a sparse matrix where each fragment of every chromosome was given a unique identifier and every pair of fragments was given a contact count if it was nonzero. The assembled genome generated by instaGRAAL was polished to remove misassemblies (Baudry et al., 2020). Initial and final assembly metrics (Nx, GC distribution) were obtained using QUAST-LG. Misassemblies were quantified using QUAST-LG with the minimap2 aligner in the backend. Ortholog and assembly completeness was computed with BUSCO (v3). The evolution of genome metrics between cycles was obtained using instaGRAAL's own implementation.
### 2.3.6 | Gene prediction

All RNAseq reads obtained from the different organ transcriptomes were quality-trimmed using Trimmomatic (version 0.36) (Bolger et al., 2014) and mapped to the genome assembly using HISAT2 (Kim et al., 2019). The alignment information was processed to generate genome-guided transcriptome assembly using Trinity (ver. 2.8.4) (Grabherr et al., 2011) and in parallel, de novo transcriptome assembly was also generated by the Trinity software. The genome-guided and de novo transcriptome assemblies were processed by the PASA pipeline (ver. 2.3.3) (Haas et al., 2003) to generate a training set for optimizing gene prediction parameters in AUGUSTUS (Stanke et al., 2006). RNAseq reads were also aligned to the genome using STAR (ver. 2.7.1a) (Dobin et al., 2013) to produce hint files of exon and intron information. Then, de novo gene prediction was performed using AUGUSTUS (ver. 2.5.5). The completeness of the predicted gene models was assessed by Benchmarking Universal Single-Copy Orthologs (BUSCO) using metazoa_odb9 database (Simão et al., 2015). Gene annotation was made using Trinotate (Grabherr et al., 2011).

### 2.3.7 | Noncoding RNA prediction

tRNA genes were predicted using tRNAscan-SE v1.3.1 (Lowe & Eddy, 1997) with eukaryote parameters, and rRNA with high conservation were predicted by aligning reads to the *Arabidopsis thaliana* template rRNA sequences using BLASTN (Altschul et al., 1997), with an e-value of $10^{-5}$. Additionally, INFERNAL (Nawrocki & Eddy, 2013) was used to predict miRNA and snRNA genes on the basis of the Rfam database (Griffiths-Jones, 2005).

## 2.4 | Repeat annotation

### 2.4.1 | LTR-retrotransposons identification

We specifically investigated LTR-retrotransposons in four genomes of Ostreidae (O. edulis, C. gigas (GCA_902806645.1), C. virginica (GCA_002022765.4), and S. glomerata (GCA_003671525.1)) and one genome of a Pteriidae species, the pearl oyster *Pinctada martensii* (GCA_002216045.1) using a detailed and precise pipeline previously customized for M. mercenaria (Farhat et al., 2022). To annotate these LTR-retrotransposons, we first extracted and translated the RT/RNaseH domain from the sequences obtained with LTRHarvest using BLASTx (Camacho et al., 2009) ($e$-value $<10^{-5}$) against an in-house database of RT/RNaseH (Thomas-Bulle et al., 2018) made from Gypsy Database (Llorens et al., 2011). Here, we kept the single sequences and/or a consensus sequence per previously defined cluster and used them in phylogenetic approaches to determine the position of the elements in each clade. Cladistic analyses were performed on amino acid sequences corresponding to the RT/RNaseH domains of the newly characterized sequences and reference

elements. Multiple alignments of these protein sequences were performed using MAFFT (Katoh et al., 2018). After manual curation of the alignments, phylogenetic analyses were conducted using neighbor-joining and the pairwise deletion option of the MEGA 5.2 software (Tamura et al., 2011). Using Topali2.3 (Milne et al., 2009), the best-fitted substitution model retained was the JTT model with a gamma distribution. Support for individual groups was evaluated with nonparametric bootstrapping using 100 replicates.

### 2.4.2 | Repeated sequences annotation

Repeated sequences were annotated in the five genomes by running RepeatMasker (http://repeatmasker.org) with the default parameter and the script "Concatenate_sequences.py" (Thomas-Bulle et al., 2018). Two different libraries were used: (i) one with only inserts-cleaned consensus from each LTR-retrotransposon cluster detected by LTRHarvest to retrieve putatively missed copies (with corrupted LTRs or deleted); (ii) the second with all other consensus transposable elements obtained with RepeatModeler v2.0.1 (Flynn et al., 2020) using REPBASE, version 2017-01-27 (Jurka et al., 2005).

## 2.5 | Phylogenetic analysis

To cluster families from protein-coding genes, proteins from the longest transcripts of each gene from O. edulis and other mollusc species, including *Argopecten irradians* (SRP174526), A. purpuratus (PRJNA418203), *Aplysia californica* (GCA_000002075.2), *Biomphalaria glabrata* (GCF_000457365.1), C. gigas (GCA_902806645.1), C. virginica (GCA_002022765.4), *Chlamys farreri* (in MolluskDB), *Elysia chlorotica* (GCA_003991915.1), *Modiolus philippinarum* (GCA_002080025.1), P. fucata (PRJNA283019), *Lottia gigantea* (PRJNA175706), *Octopus binaculoides* (GCF_001194135.1), *S. glomerata* (GCA_003671525.1), *Haliotis discus* (PRJNA317403) were collected in GENBANK and data for *Ostrea denselamellosa* and *Ostrea stentina* were specifically generated for this analysis. Orthofinder version 2.4.1 (Emms & Kelly, 2019) was run using all proteomes with default parameters. From Orthofinder results, we used the single-copy gene OGs (114OGs) and concatenated the proteins per species to generate the species tree using MAFFT online (Katoh et al., 2018) with default parameters. MEGA version 5.2 software was used to construct the phylogenetic tree with a bootstrap of 1000 on a neighbor-joining method.

## 2.6 | Differential expression analysis

RNAseq reads from *M. refringens*-infected and -noninfected libraries were checked for quality issues and adapter content with FastQC 0.11.7 and cleaning for sequencing adapters, trimming of low-quality bases (minimum mean quality score of 30), and filtering for length were performed with Trimmomatic 3.3 (Bolger et al., 2014).
Reads were aligned to the Roscoff_O.edulis-V1 genome assembly with Bowtie 2 (Langmead et al., 2019) to generate BAM files that were used to generate expression counts with the IidxStats software (Li et al., 2009). Raw read counts were first normalized using the FPKM method (Mortazavi et al., 2008), and FPKM values below 5 were eliminated when the gene was present in both conditions. Data were log2 transformed and then imported into R 3.5.0 (R Core Team, 2017) for analysis with DESeq2 package (Love et al., 2014). The genes with a threshold of false discovery rate (FDR) < 0.05 and |log2 (fold change)| > 1 were defined as differentially expressed genes (DEGs). Functional enrichment analysis for the Gene Ontology term was performed using DAVID software (Huang et al., 2009) using a Benjamini adjusted *p*-value of .05.

## 3 | RESULTS

### 3.1 | Genome assembly and annotation

PacBio reads (45× raw coverage) and Illumina PE shotgun reads (104× raw coverage) were first assembled into contigs using Masurca hybrid assembler v3.2.8 (Zimin et al., 2013) assuming a 1 Gb genome size estimated by flow cytometry using SYBGREEN labeling (data not shown). Then, the 10× Chromium reads (67× raw coverage) were used by the ARCS v1.0.4 (https://github.com/bcgsc/arcs, v1.0.4) together with LINKs pipeline v1.8.5 (https://github.com/bcgsc/LINKS, v1.8.5) to scaffold and make the Masurca assembly more contiguous. A total of 5417 contigs were obtained for an estimated size of 1028 Mb with an N50 length of 0.94 Mb. HaploMerger2 software (https://github.com/maplefores t/HaploMerger2/releases/; Huang et al., 2017) was used to remove potential duplicated sequences and significantly improved the as- sembly to 2846 contigs (N50 of 1.62 Mb). A chromosome-level assembly was then generated using Hi-C data and InstaGRAAL software. The final assembly of the Roscoff_O.edulis-V1 genome was 1018 Mb in size, with the chromosome-level scaffolds rep- resented in 10 super-scaffolds corresponding to 887.2 Mb of sequence length and 553 unplaced scaffolds with a total N50 of 97.1 Mb for scaffold lengths. The 10 expected chromosomes are in accordance with the chromosomes previously evidenced in a karyotype analysis (Thiriot-Quiévreux & Ayalaud, 1982) and link- age map (Lallias et al., 2007) (Figure S1). The Roscoff_O.edulis-V1 genome shows a GC content of 35.46% and the BUSCO analysis indicates that the gene models include 95.1% complete BUSCOs. A total of 35,962 protein-coding genes were predicted from which 24,302 genes were functionally annotated (e-value < 10⁻³). The parameters of completeness and main genome features of the as- sembly are presented in Table 1 and are very close to those ob- tained in the OE_Roslin_V1 assembly (Gundappa et al., 2022). A phylogenetic tree confirms that Roscoff_O.edulis-V1 clusters with other Ostreidae and has a closer relationship with the other two Chinese Ostrea species (Figure 1).

<table><caption>TABLE 1 Genome assembly statistics for Roscoff_O.edulis-V1 genome.</caption><thead><tr><th>Metric</th><th>Value</th><th>Number of genes</th></tr></thead><tbody><tr><td>10-chrom assembly size (bp)</td><td>887,207,661</td><td></td></tr><tr><td>Scaffold 1 (bp)</td><td>117,936,894</td><td>4971</td></tr><tr><td>Scaffold 2 (bp)</td><td>111,702,126</td><td>4882</td></tr><tr><td>Scaffold 3 (bp)</td><td>97,982,138</td><td>4282</td></tr><tr><td>Scaffold 4 (bp)</td><td>97,658,951</td><td>4001</td></tr><tr><td>Scaffold 5 (bp)</td><td>97,091,390</td><td>4170</td></tr><tr><td>Scaffold 6 (bp)</td><td>95,771,753</td><td>3908</td></tr><tr><td>Scaffold 7 (bp)</td><td>90,095,513</td><td>3357</td></tr><tr><td>Scaffold 8 (bp)</td><td>74,472,875</td><td>2356</td></tr><tr><td>Scaffold 9 (bp)</td><td>53,076,851</td><td>2190</td></tr><tr><td>Scaffold 10 (bp)</td><td>51,419,170</td><td>1846</td></tr><tr><td>GC content (%)</td><td>35.46</td><td></td></tr><tr><td>Complete BUSCOs* (C)</td><td>931 (95.1%)</td><td></td></tr><tr><td>Complete and single-copy BUSCOs (S)</td><td>911 (93.1%)</td><td></td></tr><tr><td>Complete and duplicated BUSCOs (D)</td><td>20 (2.0%)</td><td></td></tr><tr><td>Protein-coding genes number</td><td>35,963</td><td></td></tr><tr><td>mean transcript length</td><td>13,916</td><td></td></tr><tr><td>mean cds length</td><td>1626</td><td></td></tr><tr><td>mean exons per transcript</td><td>7.2</td><td></td></tr><tr><td>mean introns per transcript</td><td>5.6</td><td></td></tr></tbody></table>

### 3.2 | RNA prediction

Detection of miRNA using INFERNAL and Rfam database allowed the identification of 126 families of miRNAs (e-value < 10⁻⁵) in the Roscoff_O.edulis-V1 genome. A large variation in the copy number is also detected among miRNA families and a total of 99 miRNA families contains more than one paralog. The mir-148 family is the most represented family with 219 copies followed by the mir-821 family with 106 copies and 23 families exhibiting more than 10 copies (Table S1A). The distribution of the 25 most represented miRNA families shows a relatively homogeneous distribution along the chromosomes with a higher frequency observed in chromosomes 5, 6, 7, 8, and 9. Most of the conserved miRNA families identified in other mollusc species (let-7, lin-4, miR-2, miR-29, miR-87, miR-184, miR-8, miR-10, miR-15, miR-17) are present in Roscoff_O.edulis-V1 genome (Table S1B). The number of miRNA families identified in O. edulis (128) are similar to what is detected in C. gigas (120 families), C. virginica (121 families), or P. fucata (126 families) and many of the miRNA classically found in molluscs are present. We performed the same analysis on C. virginica and M. mercenaria genomes to identify enrichments in miRNA fami- lies in O. edulis. No miRNA family appears to be specific to O. edulis when comparing with C. virginica, but O. edulis presents a globally higher number of copies for the most represented families (Table 2). By contrast, the diversity of miRNAs in M. mercenaria appears to be
<figure><img src="image_2.png" /><figcaption>TABLE 2 Numbers of copies for the 14 most represented miRNA in the genomes of *Ostrea edulis*, *Mercenaria mercenaria*, and *Crassostrea virginica* (miRNA annotation according to RFam database).</figcaption></figure>

<table><thead><tr><th>miRNA code in RFam</th><th>O. edulis</th><th>C. virginica</th><th>M. mercenaria</th></tr></thead><tbody><tr><td>Mir-1253</td><td>0</td><td>0</td><td>29</td></tr><tr><td>Mir-148</td><td>219</td><td>157</td><td>0</td></tr><tr><td>Mir-1222</td><td>24</td><td>20</td><td>2</td></tr><tr><td>Mir-1027</td><td>76</td><td>62</td><td>0</td></tr><tr><td>Mir-1023</td><td>84</td><td>54</td><td>6</td></tr><tr><td>Let-7</td><td>17</td><td>17</td><td>2</td></tr><tr><td>Mir-1803</td><td>28</td><td>13</td><td>0</td></tr><tr><td>Mir-2118</td><td>43</td><td>27</td><td>3</td></tr><tr><td>Mir-598</td><td>6</td><td>3</td><td>26</td></tr><tr><td>Mir-1122</td><td>65</td><td>48</td><td>12</td></tr><tr><td>Mir807</td><td>23</td><td>19</td><td>4</td></tr><tr><td>Mir-944</td><td>50</td><td>47</td><td>4</td></tr><tr><td>Mir-355</td><td>31</td><td>22</td><td>0</td></tr><tr><td>Mir-2118</td><td>41</td><td>27</td><td>3</td></tr></tbody></table>

much lower and some families appear to be specific to oysters such as mir-148, mir-1027, mir-1803, and mir-355. Only mir-1253 is present in *M. mercenaria* and absent in oyster genomes (Table 2).

A total of 1985 annotated tRNAs representing the 20 amino acids, three copies of the 5.8S rRNA, six copies of the LSU-rRNA,

and six copies of SSU-rRNA, but 54 copies of the 5S-rRNA have been annotated (Table 3). Two tRNA families are over-represented with 1019 copies of the tRNA-Ser from which 907 correspond to anticodon AGA and 339 copies of the tRNA-Thr from which 282 copies correspond to AGU. tRNA-Ser gene distribution exhibits a homogenous repartition in the genome (Figure S2). The pairwise alignment scores confirm that the tRNA-Ser motif is not included in a repeated element that could explain the high number of copies observed. The comparison of tRNA composition in *M. mercenaria* and *C. virginica* genomes using the same parameters showed that enrichment in tRNA-Ser was only detected in *O. edulis*, but a higher content in tRNA-Thr (791 copies) and tRNA-Ala (465 copies) was detected in *M. mercenaria*. The density of tRNA is higher in *O. edulis* (2.275 tRNA/Mb) than in *M. mercenaria* (1.18 tRNA/Mb) or in *C. virginica* (0.575 tRNA/Mb). Comparison with available *C. gigas* tRNA composition showed a similar pattern with 641 tRNA genes and a density of 1.08 tRNA/Mb for the Pacific oyster.

### 3.3 | Transposable elements in *O. edulis* and other oyster genomes

Repeated elements annotated in the five assemblies using the same pipeline method for proper comparison (Table 4) showed a similar content of various repeated elements, except for *P. martensii* (which is not a true oyster of Ostreidae like the other four species) for which TEs represent <15% of the genome. Amon the five species, *O. edulis*
<table><caption>TABLE 3 Numbers of tRNA and rRNA in the genome of Ostrea edulis, *Mercenaria mercenaria*, and *Crassostrea virginica*.</caption><thead><tr><th></th><th>O. edulis</th><th>M. mercenaria</th><th>C. virginica</th></tr></thead><tbody><tr><td>tRNA_Ser</td><td>1019</td><td>96</td><td>53</td></tr><tr><td>tRNA_Phe</td><td>42</td><td>25</td><td>11</td></tr><tr><td>tRNA_As</td><td>43</td><td>26</td><td>37</td></tr><tr><td>tRNA_Glu</td><td>31</td><td>20</td><td>56</td></tr><tr><td>tRNA_iMet</td><td>42</td><td>30</td><td>36</td></tr><tr><td>tRNA_Tyr</td><td>30</td><td>24</td><td>17</td></tr><tr><td>tRNA_Leu</td><td>44</td><td>51</td><td>60</td></tr><tr><td>tRNA_Lys</td><td>37</td><td>45</td><td>62</td></tr><tr><td>tRNA_Gln</td><td>39</td><td>48</td><td>21</td></tr><tr><td>tRNA_Ile</td><td>34</td><td>43</td><td>26</td></tr><tr><td>tRNA_Cys</td><td>20</td><td>26</td><td>10</td></tr><tr><td>tRNA_Pro</td><td>30</td><td>42</td><td>28</td></tr><tr><td>tRNA_Arg</td><td>49</td><td>74</td><td>95</td></tr><tr><td>tRNA_Trp</td><td>8</td><td>14</td><td>9</td></tr><tr><td>tRNA-Sec</td><td>1</td><td>2</td><td>1</td></tr><tr><td>tRNAiks</td><td>28</td><td>56</td><td>34</td></tr><tr><td>tRNA_His</td><td>12</td><td>26</td><td>9</td></tr><tr><td>tRNA_Thr</td><td>339</td><td>791</td><td>153</td></tr><tr><td>tRNA_Gly</td><td>26</td><td>64</td><td>80</td></tr><tr><td>tRNA_Met</td><td>17</td><td>51</td><td>28</td></tr><tr><td>tRNA_As</td><td>21</td><td>72</td><td>21</td></tr><tr><td>tRNA_Ala</td><td>76</td><td>465</td><td>93</td></tr><tr><td>Total</td><td>1985</td><td>2091</td><td>940</td></tr><tr><td>5_8S_rRNA</td><td>3</td><td>9</td><td>5</td></tr><tr><td>5S_rRNA</td><td>54</td><td>18</td><td>99</td></tr><tr><td>LSU_rRNA</td><td>6</td><td>36</td><td>9</td></tr><tr><td>SSU_rRNA</td><td>6</td><td>18</td><td>6</td></tr><tr><td>Total</td><td>69</td><td>81</td><td>119</td></tr></tbody></table>

presents the highest coverage in TEs (38.05%). Although each element type is variably represented among species (e.g., Penelope elements are rare in both Crassostrea, as are SINEs in C. gigas and P. martensii), O. edulis almost always has the highest (LINEs, YR elements) or second-highest (Penelope, SINEs, DNA transposons, rolling circles) proportions of TEs for each element type. Only LTR-retrotransposons seem to be a little less frequent (1.61%) than in other species.

In order to characterize LTR-retrotransposons in oyster genomes, LTRharvest was used and its output was integrated into phylogenetic analyses. For this purpose, two types of sequences were used: (i) either a consensus devoid of possible insertions in the case of a cluster of elements previously defined by Uclust or (ii) or isolated (single) sequences that did not cluster with any other elements. A total of 797 subfamilies were detected by LTRharvest in the five genomes, including 14 Copia, 155 BEL/Pao, and 628 Gypsy (Table S2). The relative abundance of the three superfamilies in O. edulis (76% for Gypsy elements, 22% for BEL/Pao elements, and 2% for Copia elements) is quite similar among oyster species; except for Copia elements in P. martensii only represented by one single copy (Table 5). The phylogenetic tree of Copia elements revealed that all elements belong to the GalEa clades (Figure 2).

The O. edulis genome did not contain elements of the BEL, Pao, and Dan clades (Figure S3) but contained TEs in the six other clades in the BEL/Pao superfamily. Most of its elements belong to the Sailor lineage, mostly to the Sparrow, Sinbad, and Tas clades and only one subfamily in the Suzu clade and one element in the Flow clade are present. This distribution is globally observed in the other oysters, with a large Sparrow clade. Only O. edulis and C. gigas present all the six clades, in particular, the absence of the element of the Tas clade in both P. martensii and S. glomerata.

Gypsy superfamily tree reveals 12 clades in oysters (Figure 3, Figure S4). Ten correspond to the MolGy clades previously defined from mollusc elements (except for clades MolGy11 and MolGy13) and to the clades A-B and C, also known to have Gypsy elements of molluscs. Two new putative clades can be further identified with oyster elements: MolGy17, observed in all five species, and MolGy18 that is only absent in C. virginica. O. edulis does not present any element of the clade MolGy6, although it is quite well represented in the four other species; and it also has very few elements of the MolAn clade (Table 5). Three clades (Molgy1, Molgy2, and the C-clade) are well recovered, represent about 1.26% of the genome, and largely dominate with three-quarters of the subfamilies as observed in other oysters.

## 3.4 | Differential expression analysis in M. refringens-infected oysters

In the hemocytes, seventy-eight DEGs are upregulated (Figure 4a) where 4 terms in the Biological Process category (GO:0006259 DNA metabolic process, GO:0090305 nucleic acid phosphodiester bond hydrolysis, GO:0006302 double-strand break repair, GO:0007155 cell adhesion), and 2 terms in the Molecular Function category (GO:0003682 Chromatin Binding and GO:0003677 DNA binding) are enriched (Table S4). Among the most downregulated regulated genes, we identified IAPs, mannose receptor C-type, serine protease inhibitor, one SOD, and several genes encoding C1q complements, which are all involved in the immune response. Among upregulated genes, we identified several proteins involved in apoptosis such as ced-1 and other genes involved in cell adhesion such as a sushi, von Willebrand factor type A, an EGF, and the pentraxin domain-containing protein 1 and several TRIM proteins. A total of 143 DEGs are downregulated in M. refringens-infected oysters (Figure 4b). GOterm enrichment analysis shows that 4 terms in the Biological Process category (GO:0007155 cell adhesion, GO:0032715 negative regulation of interleukin-6 production, GO:0032088 negative regulation of NF-kappaB transcription factor, and GO:0070373 negative regulation of ERK1 and ERK2 cascade), 5 terms in the Cellular Component category (GO:0005581 collagen trimer, GO:0005576 extracellular region, GO:0005615
<table><caption>TABLE 4 Copy number and genomic content of the repetitive elements in *Ostrea edulis*, *Crassostrea gigas*, *Crassostrea virginica*, *Pinctada martensii*, and *Saccostrea glomerata* genomes generated following the annotation pipeline described in this paper.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Ostrea edulis</th><th colspan="3">Crassostrea gigas</th><th colspan="3">Crassostrea virginica</th><th colspan="3">Pinctada martensii</th><th colspan="3">Saccostrea glomerata</th></tr><tr><th>% genome</th><th>Copies Nb</th><th>Nb of bp</th><th>% genome</th><th>Copies Nb</th><th>Nb of bp</th><th>% genome</th><th>Copies Nb</th><th>Nb of bp</th><th>% genome</th><th>Copies Nb</th><th>Nb of bp</th><th>% genome</th><th>Copies Nb</th><th>Nb of bp</th></tr></thead><tbody><tr><td>Retroelements</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Penelope</td><td>4.46</td><td>57,601</td><td>38,234,354</td><td>0.49</td><td>4421</td><td>3,174,681</td><td>0.11</td><td>764</td><td>740,415</td><td>5.84</td><td>90,093</td><td>57,830,428</td><td>3.01</td><td>39,094</td><td>23,757,514</td></tr><tr><td>LINES</td><td>3.05</td><td>31,673</td><td>26,178,491</td><td>2.83</td><td>18,815</td><td>18,348,219</td><td>1.79</td><td>12,916</td><td>12,273,678</td><td>0.94</td><td>10,497</td><td>9,359,382</td><td>1.53</td><td>12,860</td><td>12,064,399</td></tr><tr><td>L2</td><td>0.43</td><td>4915</td><td>3,721,701</td><td>0.04</td><td>310</td><td>239,915</td><td>0.03</td><td>237</td><td>210,627</td><td>0.09</td><td>1076</td><td>925,824</td><td>0.28</td><td>1812</td><td>2,188,757</td></tr><tr><td>CR1</td><td>0.23</td><td>2999</td><td>2,005,708</td><td>0.37</td><td>2610</td><td>2,399,937</td><td>0.52</td><td>3921</td><td>3,564,568</td><td>0.20</td><td>2864</td><td>1,989,573</td><td>0.24</td><td>2135</td><td>1,900,378</td></tr><tr><td>RTE-X</td><td>0.82</td><td>8206</td><td>7,016,657</td><td>0.82</td><td>6001</td><td>5,331,699</td><td>0.40</td><td>3138</td><td>2,769,765</td><td>0.18</td><td>2023</td><td>1,799,210</td><td>0.39</td><td>3578</td><td>3,071,346</td></tr><tr><td>I</td><td>0.33</td><td>5595</td><td>2,804,522</td><td>0.00</td><td>0</td><td>0</td><td>0.13</td><td>1491</td><td>883,115</td><td>0.00</td><td>0</td><td>0</td><td>0.07</td><td>867</td><td>565,958</td></tr><tr><td>L1</td><td>0.84</td><td>6494</td><td>7,179,616</td><td>0.75</td><td>3694</td><td>4,841,718</td><td>0.59</td><td>3301</td><td>4,023,806</td><td>0.29</td><td>2821</td><td>2,901,900</td><td>0.2</td><td>2489</td><td>2,529,348</td></tr><tr><td>RTE</td><td>0.33</td><td>2846</td><td>2,829,694</td><td>0.82</td><td>6001</td><td>5,331,699</td><td>0.07</td><td>461</td><td>456,806</td><td>0.15</td><td>1487</td><td>1,494,237</td><td>0.17</td><td>1548</td><td>1,331,519</td></tr><tr><td>Others</td><td>0.07</td><td>618</td><td>620,593</td><td>0.03</td><td>199</td><td>203,251</td><td>0.05</td><td>367</td><td>364,852</td><td>0.03</td><td>226</td><td>248,638</td><td>0.06</td><td>431</td><td>477,093</td></tr><tr><td>SINES</td><td>1.13</td><td>17,644</td><td>9,710,552</td><td>0.06</td><td>396</td><td>370,098</td><td>1.27</td><td>15,408</td><td>8,696,250</td><td>0.28</td><td>3787</td><td>2,755,480</td><td>0.65</td><td>11,734</td><td>5,149,632</td></tr><tr><td>LTR elements</td><td>1.61</td><td>9362</td><td>14,030,935</td><td>2.51</td><td>10,438</td><td>16,258,283</td><td>1.63</td><td>8440</td><td>11,191,940</td><td>1.00</td><td>6938</td><td>9,867,404</td><td>1.94</td><td>11,497</td><td>15,278,733</td></tr><tr><td>Gypsy</td><td>1.33</td><td>8012</td><td>11,588,206</td><td>2.07</td><td>9159</td><td>13,397,133</td><td>1.41</td><td>7519</td><td>9,650,060</td><td>0.88</td><td>6273</td><td>8,757,301</td><td>1.70</td><td>10,292</td><td>13,403,278</td></tr><tr><td>BEL/PAO</td><td>0.26</td><td>1237</td><td>2,229,921</td><td>0.42</td><td>1231</td><td>2,734,597</td><td>0.21</td><td>795</td><td>1,410,415</td><td>0.07</td><td>503</td><td>706,175</td><td>0.22</td><td>1144</td><td>1,745,192</td></tr><tr><td>Copia</td><td>0.02</td><td>79</td><td>189,857</td><td>0.02</td><td>48</td><td>126,553</td><td>0.01</td><td>45</td><td>100,528</td><td>0.00</td><td>0</td><td>0</td><td>0.02</td><td>61</td><td>130,263</td></tr><tr><td>Others</td><td>0.00</td><td>34</td><td>22,951</td><td>0.00</td><td>0</td><td>0</td><td>0.00</td><td>81</td><td>30,987</td><td>0.04</td><td>62</td><td>403,928</td><td>0.00</td><td>0</td><td>0</td></tr><tr><td>YR elements</td><td>0.65</td><td>5111</td><td>5,562,901</td><td>0.51</td><td>2588</td><td>3,317,411</td><td>0.33</td><td>1613</td><td>2,226,523</td><td>0.12</td><td>1028</td><td>1,163,186</td><td>0.37</td><td>2657</td><td>2,887,324</td></tr><tr><td>Ngaro</td><td>0.44</td><td>3150</td><td>3,810,113</td><td>0.34</td><td>1768</td><td>2,230,023</td><td>0.19</td><td>982</td><td>1,287,955</td><td>0.02</td><td>196</td><td>212,583</td><td>0.27</td><td>1810</td><td>2,115,549</td></tr><tr><td>DIRS</td><td>0.20</td><td>1961</td><td>1,752,788</td><td>0.16</td><td>820</td><td>1,087,388</td><td>0.14</td><td>631</td><td>938,568</td><td>0.10</td><td>832</td><td>950,603</td><td>0.10</td><td>847</td><td>771,775</td></tr><tr><td>DNA transposons</td><td>7.08</td><td>83,982</td><td>60,748,905</td><td>7.32</td><td>74,071</td><td>47,420,818</td><td>5.34</td><td>46,434</td><td>36,532,449</td><td>3.56</td><td>57,625</td><td>35,264,480</td><td>5.88</td><td>69,471</td><td>46,338,311</td></tr><tr><td>Maverick</td><td>0.28</td><td>751</td><td>2,438,022</td><td>0.22</td><td>764</td><td>1,393,542</td><td>0.38</td><td>1133</td><td>2,580,104</td><td>0.19</td><td>991</td><td>1,897,408</td><td>0.08</td><td>701</td><td>600,613</td></tr><tr><td>TcMar/Pogo</td><td>1.38</td><td>14,369</td><td>11,848,554</td><td>1.19</td><td>12,709</td><td>7,710,105</td><td>0.49</td><td>4683</td><td>3,333,461</td><td>0.48</td><td>7725</td><td>4,471,920</td><td>1.38</td><td>13,911</td><td>10,870,057</td></tr><tr><td>Hobo/Ac/Tam</td><td>0.95</td><td>12,619</td><td>8,173,678</td><td>0.39</td><td>3002</td><td>2,520,162</td><td>0.46</td><td>4165</td><td>3,166,922</td><td>0.80</td><td>12,909</td><td>7,915,023</td><td>0.49</td><td>5663</td><td>3,866,272</td></tr><tr><td>Zator</td><td>0.38</td><td>4849</td><td>3,295,890</td><td>0.15</td><td>1481</td><td>960,144</td><td>0.04</td><td>425</td><td>293,547</td><td>0.15</td><td>3014</td><td>1487,162</td><td>0.14</td><td>1896</td><td>1,064,820</td></tr><tr><td>Crypton</td><td>1.37</td><td>20,899</td><td>11,730,201</td><td>1.78</td><td>22,194</td><td>11,503,648</td><td>1.09</td><td>12,912</td><td>7,489,797</td><td>0.75</td><td>15,571</td><td>7,448,289</td><td>1.32</td><td>19,279</td><td>10,442,231</td></tr><tr><td>Academ</td><td>0.25</td><td>2649</td><td>2,453,984</td><td>0.03</td><td>167</td><td>213,307</td><td>0.13</td><td>651</td><td>894,907</td><td>0.08</td><td>791</td><td>743,310</td><td>0.03</td><td>229</td><td>274,905</td></tr><tr><td>EnSpm</td><td>0.09</td><td>705</td><td>809,979</td><td>0.32</td><td>1798</td><td>2,080,891</td><td>0.47</td><td>3155</td><td>3,220,923</td><td>0.16</td><td>2902</td><td>1,579,436</td><td>0.31</td><td>2877</td><td>2,450,110</td></tr><tr><td>MULE/MuDR/<br/>I5905</td><td>0.13</td><td>1341</td><td>1,130,079</td><td>0.04</td><td>328</td><td>288,206</td><td>0.10</td><td>806</td><td>672,632</td><td>0.07</td><td>875</td><td>678,893</td><td>0.08</td><td>701</td><td>600,613</td></tr><tr><td>PHIS</td><td>1.22</td><td>12,327</td><td>10,474,972</td><td>1.56</td><td>14,069</td><td>10,077,041</td><td>1.14</td><td>9623</td><td>7,838,419</td><td>0.66</td><td>10,261</td><td>6,499,076</td><td>1.09</td><td>13,198</td><td>8,619,218</td></tr><tr><td>Polobok</td><td>0.58</td><td>8162</td><td>5,001,169</td><td>1.24</td><td>14,339</td><td>8,041,774</td><td>0.61</td><td>6107</td><td>4,177,778</td><td>0.03</td><td>385</td><td>334,047</td><td>0.65</td><td>7877</td><td>5,160,850</td></tr><tr><td>PiggyBac</td><td>0.11</td><td>1710</td><td>965,878</td><td>0.15</td><td>923</td><td>962,318</td><td>0.13</td><td>485</td><td>901,089</td><td>0.00</td><td>0</td><td>0</td><td>0.09</td><td>1064</td><td>748,137</td></tr><tr><td>Sola</td><td>0.16</td><td>1741</td><td>1,401,559</td><td>0.15</td><td>896</td><td>988,546</td><td>0.13</td><td>855</td><td>871,136</td><td>0.11</td><td>990</td><td>1,111,636</td><td>0.09</td><td>673</td><td>703,386</td></tr><tr><td>Others</td><td>0.15</td><td>1840</td><td>1,324,940</td><td>0.11</td><td>791</td><td>681,134</td><td>0.16</td><td>1434</td><td>1,091,734</td><td>0.08</td><td>1031</td><td>820,280</td><td>0.12</td><td>1402</td><td>937,099</td></tr><tr><td>Rolling circles</td><td>9.88</td><td>110,094</td><td>84,753,110</td><td>12.47</td><td>79,436</td><td>80,781,395</td><td>9.48</td><td>72,591</td><td>64,896,885</td><td>3.83</td><td>54,233</td><td>37,948,915</td><td>9.88</td><td>110,094</td><td>84,753,110</td></tr><tr><td>Unclassified</td><td>10.20</td><td>131,922</td><td>87,508,527</td><td>3.90</td><td>31,005</td><td>25,285,376</td><td>8.74</td><td>82,096</td><td>59,864,032</td><td>0.09</td><td>1296</td><td>867,014</td><td>11.24</td><td>131,498</td><td>88,620,649</td></tr></tbody></table>
<table><caption>TABLE 5 Genomic proportions of the different clades of LTR-retrotransposons detected in *Ostrea edulis*, *Crassostrea gigas*, *Crassostrea virginica*, *Pinctada martensii*, and *Saccostrea glomerata* genomes.</caption><thead><tr><th></th><th>O. edulis</th><th>C. gigas</th><th>C. virginica</th><th>P. martensii</th><th>S. glomerata</th></tr></thead><tbody><tr><td>Genome size</td><td>872,374,424 bp</td><td>647,887,097 bp</td><td>684,723,884 bp</td><td>99,0984,031 bp</td><td>788,100,799 bp</td></tr><tr><td>GalEa</td><td>0.022 (4)</td><td>0.020 (4)</td><td>0.015 (2)</td><td>s</td><td>0.017 (4)</td></tr><tr><td>Flow</td><td>s</td><td>0.015 (1)</td><td>-</td><td>0.003 (1)</td><td>s</td></tr><tr><td>Suzu</td><td>0.005 (1)</td><td>0.015 (2)</td><td>0.006 (1)</td><td>s</td><td>0.006 (1)</td></tr><tr><td>Tas</td><td>0.027 (4)</td><td>0.019 (2)</td><td>0.041 (3)</td><td>-</td><td>-</td></tr><tr><td>Sinbad</td><td>0.034 (5)</td><td>0.036 (3)</td><td>0.027 (2)</td><td>0.004 (2)</td><td>0.015 (2)</td></tr><tr><td>Sparrow</td><td>0.187 (26)</td><td>0.300 (33)</td><td>0.127 (16)</td><td>0.055 (9)</td><td>0.161 (25)</td></tr><tr><td></td><td>s</td><td>0.021 (2)</td><td>0.005 (1)</td><td>0.004 (1)</td><td>0.011 (2)</td></tr><tr><td>Unclassified</td><td>0.002 (1)</td><td>0.017 (2)</td><td>-</td><td>0.005 (1)</td><td>0.029 (6)</td></tr><tr><td>Pao/BEL</td><td>0.256</td><td>0.422</td><td>0.206</td><td>0.071</td><td>0.221</td></tr><tr><td>AB-clade</td><td>0.034 (6)</td><td>0.039 (7)</td><td>0.025 (3)</td><td>-</td><td>0.034 (8)</td></tr><tr><td>C-clade</td><td>0.142 (28)</td><td>0.117 (13)</td><td>0.138 (18)</td><td>0.053 (16)</td><td>0.133 (23)</td></tr><tr><td>MolGy1</td><td>0.347 (23)</td><td>0.376 (32)</td><td>0.290 (25)</td><td>0.075 (6)</td><td>0.319 (26)</td></tr><tr><td>MolGy2</td><td>0.584 (47)</td><td>0.631 (51)</td><td>0.450 (29)</td><td>0.208 (21)</td><td>0.672 (52)</td></tr><tr><td>MolGy3</td><td>0.062 (11)</td><td>0.039 (5)</td><td>0.039 (6)</td><td>0.024 (6)</td><td>0.053 (10)</td></tr><tr><td>MolGy4</td><td>0.014 (1)</td><td>0.250 (7)</td><td>s</td><td>0.213 (4)</td><td>0.044 (2)</td></tr><tr><td>MolGy5</td><td>0.050 (2)</td><td>0.005 (1)</td><td>0.099 (2)</td><td>0.220 (8)</td><td>0.121 (3)</td></tr><tr><td>MolGy6</td><td>-</td><td>0.216 (11)</td><td>0.079 (6)</td><td>0.018 (2)</td><td>0.070 (4)</td></tr><tr><td>MolGy10</td><td>0.006 (1)</td><td>0.018 (1)</td><td>0.011 (1)</td><td>-</td><td>0.013 (2)</td></tr><tr><td>MolGy17</td><td>0.0182 (2)</td><td>s</td><td>s</td><td>0.007 (1)</td><td>0.011 (2)</td></tr><tr><td>MolGy18</td><td>0.005 (1)</td><td>0.006 (1)</td><td>-</td><td>0.003 (1)</td><td>0.010 (1)</td></tr><tr><td>MolAn</td><td>0.026 (2)</td><td>0.262 (16)</td><td>0.223 (10)</td><td>s</td><td>0.132 (14)</td></tr><tr><td>Unclassified</td><td>0.041 (7)</td><td>0.109 (8)</td><td>0.056 (9)</td><td>0.064 (12)</td><td>0.089 (11)</td></tr><tr><td>Gypsy</td><td>1328</td><td>2068</td><td>1409</td><td>0.884</td><td>1701</td></tr><tr><td>LTR-retrotransposons</td><td>1606</td><td>2509</td><td>1630</td><td>0.955</td><td>1939</td></tr></tbody></table>

Note: The number of subfamilies of each clade is given in brackets. “Unclassified” = elements not linked to a clade; “-” = no element detected; “s” = only single sequence detected.

extracellular space, GO:0031012 extracellular matrix, and GO:0030054 cell junction) and 3 terms in the Molecular Function category (GO:0005102 receptor binding, GO:0030020 extracellular matrix structural constituent conferring tensile strength and GO:0005518 collagen binding) are significantly over-represented in the downregulated DEGs (Table S3).

In the digestive gland, a total of 926 DEGs have been identified as over-expressed in *M. refringens*-infected oysters (Figure 5a) and represent 29 terms in the Biological Process category, 28 terms in the Cellular Component category, and 12 terms in the Molecular Function category (Table S5). Among the genes showing the strongest upregulations in *M. refringens*-infected individuals (log2 fold >4), several genes are known to be involved in immune response processes or cancer-like pathologies (prominin, DMBT1, 2-protein convertase subtilisin/kexin), cell adhesion (contactin, laminin). Several genes involved in DNA, RNA, or histone methylation process are also upregulated (log2 fold >2.5) such as lysine methyltransferase. We also identified a set of proteins known to specifically

interact with microorganisms such as meprin (three genes), coiled-coil domain-containing proteins (18 genes), complement component 1q (11 genes), and GTPase, IMAP family member proteins (seven genes), and four genes encoding the interferon-induced protein. Specific genes involved in the late response to apoptosis such as caspase 3 and 6 are also identified. We also identified members of transposable elements (3 Pif transposases, 3 LINE-1 transposase, and 3 Tigger transposable element).

We identified 1640 DEGs downregulated in *M. refringens*-infected oysters (Figure 5b) representing 68 terms in the Biological Process category, 22 terms in the Cellular Component category, and 33 terms in the Molecular Function category (Table S6). Among those DEGs, we identified three genes encoding Cadherin 12 that are involved in cell adhesion homeostasis, one gene involved in the protection of germline (Mage protein), one gene encoding the hemocytin, which is involved in hemocyte aggregation, and a gene encoding ADAR that is involved in RNA editing. Several genes involved in the positive regulation of transcription from RNA
<figure><img src="image_3.png" /><figcaption>FIGURE 2 Phylogenetic relationships of Copia retrotransposons. The tree is based on neighbor-joining analysis of RT/RNaseH domain amino acid sequences. The Copia subfamilies from oysters are indicated in color: *Crassostrea gigas* in dark blue, *Crassostrea virginica* in light blue, *Ostrea edulis* in red, *Pinctada martensii* in purple, and *Saccostrea glomerata* in green. Node statistical support values (>70%) come from nonparametric bootstrapping using 100 replicates.</figcaption></figure>

polymerase II promoter and 14 genes encoding G-protein-coupled receptors that convert extracellular signals into intracellular responses have been identified. We also detected a total of 77 genes encoding several families of solute carrier proteins (SLCs) as being downregulated in the presence of *M. refringens*. A set of genes also known to be involved in relationships with microorganisms have been detected including ABC transporters (10 genes), ADAM proteins (six genes), several CYP proteins (five genes), several proteins containing a tyrosine phosphatase or kinase activity (35 genes), 22 members of the transmembrane proteins family (TMEM), nine genes containing a von Willebrand factor D and EGF domains (VWDE), several LDL receptor-related proteins (11 Lrp genes), and six genes encoding carboxypeptidases. Four genes belonging to the B-cell lymphoma two family that are involved in the apoptotic

response, including two proapoptotic members (Bad, Bak) and two antiapoptotic members (BCL2), have been evidenced.

## 3.5 | SLCs distribution in *O. edulis* genome

A total of 511 SLC genes classified in 48 families (e-value < 10⁻⁴) were identified, which is similar to what was observed in *C. gigas* (510 SLCs) but lower than the 673 genes identified in the scallop *P. yessoensis* (Xun et al., 2020). The distribution of SLC across the 10 chromosomes showed differences in SLC repartition with a lower number of SLC genes on chromosomes 6, 9, and 10 (Table S7). The number of SLC genes per family strongly varies and SLC5, SLC6, SLC16, SLC17, SLC25, and SLC39 families have more than 20 genes. An enrichment
<figure><img src="image_4.png" /></figure>

FIGURE 3 Phylogenetic relationships of Gypsy retrotransposons. The tree is based on neighbor-joining analysis of RT/ RNaseH domain amino acid sequences. The clades that have elements of oysters are indicated in color, and the clades not present in oysters are represented by the black lines. Node statistical support values (>70%) come from nonparametric bootstrapping using 100 replicates.

<figure><img src="image_5.png" /></figure>

FIGURE 4 Heat maps of the most upregulated DEGs (a) and the most downregulated DEGs (b) in hemocytes of *M. refringens*-infected *O. edulis* individuals. The heat map shows the matrix of fold changes that were calculated for each condition and each DEG by normalizing the expression of each condition and DEG to the expression of the DEG in all conditions. Positive fold changes are in red, and negative fold changes are in green.

in SLC16 (70 genes) is detected in *O. edulis* (versus 54 in *C. gigas* and 75 in *P. yessoensis*) and 28 of the copies are present in the chromosome 8 and 11 in the chromosome 1 including a cluster of 8 genes. The SLC6 family is mainly present in the chromosome 1 (24 copies) and 2 (16 copies) and the SLC5 is mainly present in the chromosome 5. We analyzed the expression of the 511 SLC in the different organs used for the genome annotation (Table S7) and showed that 91 SLCs are significantly upregulated in the digestive gland compared with all other organs (Figure 6). Among these 91 genes, 26 are among the downregulated SLCs identified in the *M. refringens*-infected digestive gland. The largest number of digestive gland over-expressed SLCs are from SLC16 (15), SLC6 (16), PySLC5 (9), and SLC22 (6).

4 | DISCUSSION

The use of the combined PacBio long reads, Illumina shotgun, and 10X sequencing approaches allowed us to obtain a high-quality reference genome assembled at a chromosomal level for the flat oyster *O. edulis*, which was used for subsequent genomic analysis. With a size of 887 Mb distributed along the 10 predicted chromosomes and good assembly statistics, the *Roscoff_O.edulis-V1* assembly ranks well among other sequenced genomes of bivalve molluscs (Farhat et al., 2022; Li et al., 2021; Peñaloza et al., 2021). The *Roscoff_O.edulis-V1* genome size is larger than that of the Crassostreidae (around 600 Mb), but the number of protein-coding genes predicted (35,962)
**FIGURE 5** Heat maps of the most upregulated DEGs (a) and the most downregulated DEGs (b) in the digestive gland of *M. refringens*-infected *O. edulis* individuals. The heat map shows the matrix of fold changes that were calculated for each condition and each DEG by normalizing the expression of each condition and DEG to the expression of the DEG in all conditions. Positive fold changes are in red, and negative fold changes are in green.

<figure><img src="image_6.png" /></figure>

is close to that obtained in other oyster species such as *C. gigas* (31,371), *C. virginica* (34,587), and *P. fucata* (31,477).

**4.1 | Analysis of tRNA and miRNA reveals specificity in Ostrea genome**

Transfer RNAs play an essential role in cellular life and in many metabolic processes (Francklyn & Minajigi, 2010). In many species, the number of genes encoding each tRNA generally correlates with the amino acid frequency with some exceptions such as in *Arabidopsis thaliana*, where a higher number of genes encoding tRNAs-Tyr, Ser, and Pro are detected (Michaud et al., 2011). Counting of serine residues in the *O. edulis* proteome does not reveal an excess of this residue over other amino acids, failing to explain the high tRNA-Ser copy number. The frequency of serine amino acid in the proteomes of three other molluscs (*C. gigas*, *C. virginica*, and *M. mercenaria*) does not show any significant difference (ranging from 7.9% to 8.28%). The same result was obtained when considering the proportion of threonine residues in the proteome despite a high number of tRNA-Thr genes in *O. edulis*. The number of tRNA genes identified in the three species we analyzed is consistent with what is found in other molluscan species such as *Scapharca broughtonii* with 1541 tRNA genes (Bai et al., 2019) or *C. gigas* with 641 tRNA genes (Peñaloza et al., 2021). In eukaryotes, no general trend in the number of tRNA genes is evident and the strong variations observed could reflect differences in the evolutionary history of lineages. While the study of tRNA genomic organization and

structural evolution remains partial (Goodenbour & Pan, 2006), tRNA genes have been shown to play a role in genome organization by acting as barriers to DNA replication fork progression (McFarlane & Whitehall, 2009) or by contributing to chromosomal instability (Admire et al., 2006). No clear relation between tDNA copy number and their organization in clusters in genomes is evident, and most tDNA clusters are small, containing only a few colocalized tRNA genes with the exception of Nematostella and some fish species (Bermudez-Santana et al., 2010). Despite a high number of the tRNA-Ser and tRNA-Thr that also present a high disequilibrium in the proportion of isoacceptors in the *Roscoff_O.edulis*-V1 genome, no specific clustering of those tRNA genes is observed. Variations in codon usage might influence the variation in tDNA copy number (Rocha, 2004), but no strong disequilibrium for Ser resides is evident in *O. edulis* (Gerdol et al., 2015). Differences in the number of rDNA loci between individuals of the same species (e.g., *M. galloprovincialis*), and in the location of the rDNA loci between different cells from the same individual have been found (Insua & Mendez, 1998). Additional analyses are needed to better understand whether and how the enrichment observed for the tRNA-Ser and tRNA-Thr plays a functional role in *O. edulis* physiology because this specific pattern is not observed in other oyster genomes and no difference in amino acid composition has been detected between proteomes.

Even though molluscs represent the second largest phylum in Metazoa, only a limited number of studies have been performed on the identification and the role of miRNA despite their essential role in cellular development, proliferation, apoptosis, oncogenesis,
<figure><img src="image_7.png" /><figcaption>FIGURE 6 Heat map of the 30 most upregulated SLCs in the digestive gland (Digl) of *Ostrea edulis* compared with palp (Palp), mantle (MT), adductor muscle (Mus), gills (Gill), gonads (Gon), and hemocytes (Hemo). The heat map shows the matrix of fold changes that were calculated for each tissue by normalizing the expression of each SLC within a tissue to the expression of the SLC in all tissues. Positive fold changes are in red, and negative fold changes are in green.</figcaption></figure>

differentiation, and disease (Bueno & Malumbres, 2011). While more
than 48,800 miRNA are available in all species (www.mirbase.org),
only 245 families have been identified in molluscs mainly due to the
low number of genomes available, a lack of specific miRNA annota-
tion, and a limited number of miRNA transcriptomic analysis. The
comparison of the miRNAs copy number shows differences within a
species between the different approaches used for their annotation
limiting the ability to make reliable comparisons. In C. gigas, eight
copies of miRNA-184 are identified by Huang et al. (2021) versus
44 by Rosani et al. (2021). Further, differences are observed when
comparing predicted miRNA and miRNA identified in small-RNAseq
sequencing, which do not confirm that the predicted miRNAs are
actually expressed. In the Roscoff_O.edulis-V1, we identified several
miRNA families that exhibit high copy numbers distributed along the
10 chromosomes without specific clustering. The most represented
family is miRNA-148 with 219 copies, which is known to act as a
negative regulator of MyD88-dependent NF-κB signaling in the te-
leost fish (Chu et al., 2017) and inhibited in the herbivorous carp
Ctenopharyngodon idella, in response to bacterial infection (Fang
et al., 2020). The Mir-841 gene with 106 copies in O. edulis genome is
involved in response to metabolic stress (Nischal et al., 2014) or virus
(Moyo et al., 2017) in the plants. The involvement of several miRNAs
in response to the herpes virus OsHV-1 has been demonstrated in
C. gigas (Rosani et al., 2020), Chlamys farreri (Chen et al., 2014). In
O. edulis, miRNA involved in immune system regulation (miRNA1,
miRNA10, miRNA8, miRNA33) has been shown to be regulated in
response to B. ostreae (Martín-Gómez et al., 2014). Further work

should focus on coupling miRNA prediction and miRNA transcrip-
tomic profiling in O. *edulis*.

4.2 | Transposable elements

A comparison of the TE composition in the four Ostreidae spe-
cies and the Pteriidae species analyzed shows a relative homoge-
neity of their content. The percentages of genome coverage are
close for most types of elements (LINEs, SINEs, DNA transposons,
YR elements) and the only remarkable difference in O. edulis is the
over-representation of Penelope elements (4.46% vs. 0.5% in other
species), and a global under-representation of LTR-retrotransposons
(1.61% in O. edulis vs. 4.37%). With the exception of the MolGy6
clade, O. edulis possesses all the clades detected in oysters, and 17
of the 23 clades described in molluscs are present in oysters. A tar-
geted study of Gypsy elements within 20 bivalve genomes revealed
20 branches within the Gypsy clade C (Farhat et al., 2022), and 18
are found in oyster genomes. Analysis of the five genomes has fur-
ther defined two new clades, MolGy17 and MolGy 18 that appear to
be specific to Ostreidae and Pteriidae. Specific elements such as the
Steamer LTR-retrotransposon family, which is associated with neo-
plasia in other bivalve species (Metzger et al., 2018) have not been
detected in O. edulis genome while it is present in other oysters like
C. gigas, C. virginica, and S. glomerata (Farhat et al., 2022). Contrary to
C. gigas, no specific enrichment in the rolling-circle transposable ele-
ments Helitrons is detected in O. edulis. However, the high percent-
age of unclassified TEs in both O. edulis and S. glomerate genomes
suggests that flat oysters could contain new families of ETs not yet
described. Transposable elements are known to have a large impact
on genome structure and stability, and are therefore considered
to play an important role in many evolutionary and ecological pro-
cesses such as biodiversity generation, stress response, adaptation,
speciation, and colonization mechanisms in eukaryotes (Biemont &
Vieira, 2006; Finnegan, 2012). Environmental variations can pro-
mote genome plasticity through transcriptional activation and TEs
mobilization, often in response to specific stimuli such as biotic
stresses (pathogen) or abiotic environmental changes (temperature,
heavy metals, UV) (Capy et al., 2000; Melayah et al., 2001). This vari-
ability transmitted to the offspring would also favor a better adapta-
tion of the organisms to environmental changes (long-term effect
detectable in the genome). Since several TEs appear to be regulated
in response to M. refringens, an extensive analysis of TE regulation in
response to biotic and abiotic environmental parameters should be
conducted to better understand the role of TE in flat oyster biology.

4.3 | Transcriptomic regulation in response to M. refringens

Studies of the responses of *Ostrea edulis* to the presence of *M. re-
fringens* have not yet been carried out, despite the strong impact of
*M. refringens* in terms of mortality on the oyster beds. The present
study carried out on individuals naturally infected by *M. refringens*, although preliminary, highlighted some metabolic pathways classically involved in parasite interaction but also new pathways that could be more specifically associated with the response to a eukaryote parasite. The low number of genes regulated in hemocytes in response to *M. refringens* indicates that despite its key role in immunity, this tissue appears to play only a minor role in the response to this parasite. On the contrary, the high number of differentially regulated genes and corresponding biological functions in the digestive gland emphasizes the systemic response of *O. edulis* to *M. refringens*. In hemocytes, a downregulation occurred for genes encoding IAPs, mannose receptor C-type, serine protease inhibitor, and C1q complement that all contribute to the molecular response to parasites (Chan et al., 2021; Li et al., 2019). In the digestive gland, 11 genes coding for C1q different from those identified in hemocytes are instead upregulated demonstrating differential regulation of this protein family. C1q domain-containing proteins are characterized by a large diversity of gene number across species (Farhat et al., 2022; Gerdol et al., 2019; Mun et al., 2017; Peng et al., 2020) and have been shown to be over-represented in *O. edulis* genome (Gundappa et al., 2022). Further studies will therefore be necessary to better understand the recognition specificities of this set of molecules. Serine protease inhibitor is involved in *Perkinsus marinus* resistance in oyster *C. virginica* (He et al., 2012; Yu et al., 2011) and their downregulation in *M. refringens*-infected oysters may reflect a lower ability to inhibit the serine proteases produced by the parasite. In both tissues, the increased expression of genes involved in apoptosis is coherent with previous results obtained on *O. edulis* in response to *B. ostreae* (Gervais et al., 2019) or *C. virginica* in response to *Perkinsus marinus* (Lau et al., 2018) In the digestive gland, several effector caspases (casp 3 and casp 6) are over-expressed in infected individuals while paradoxically, two pro-apoptotic genes involved in the early phases of regulation of this process (Bak and Bad genes) and two antiapoptotic genes (BCL2) are inhibited, reflecting the existence of control on the pro-anti apoptosis balance. The apoptotic response in molluscs that involves many molecular actors (Kiss, 2010; Witkop et al., 2022) is under the control of a wide diversity of molecular modulators and understanding their regulation by *M. refringens* remains an area of research to be explored further.

The identification of genes previously shown to respond to other models of parasite infection confirms that the transcriptomic response observed may result from the *M. refringens* effects, highlighting the complexity of the molecular targets affected/modulated by this parasite. We particularly found that genes involved in cilium motility or assembly or genes encoding G-protein-coupled receptors (14 genes) that are involved in the regulation of several cellular processes by sensing molecular cues outside the cells may play a role in the immune response (Bezares-Calderon et al., 2020; He et al., 2015; Song et al., 2015). The identification of increased expression in infected oysters of several genes involved in DNA, RNA, or histone methylation process suggests that epigenetic processes could also participate in the response to *M. refringens*. Indeed, a growing number of studies suggested that protozoan parasites such as *Leishmania*, *Toxoplasma*,

and *Theileria* manipulate host cells via epigenetic modification of host gene expression mainly resulting in permanent downregulation of host defense mechanisms to promote intracellular replication and survival of the pathogen (McMaster et al., 2016). As the number of studies on the regulation of methylation processes by parasites in molluscs is limited, the case of the *M. refringens*-O. edulis interaction could pave the way for a new research axis. In response to *M. refringens* infection, the downregulation of a large set of genes encoding SLC is detected. The immune system coordinates complex signals to support the proliferation, differentiation, and effector function to protect cells against tumors and infections and SLCs play an important role in these regulations and participate in the modulation of various metabolic pathways by import and export of a large variety of small molecules (nutrients, drugs...) across biological membranes (Rives et al., 2017). Specific members of SLC families (SLC16A4, SLC39) play a role in viral replication resistance factor, providing a potential protective effect against the disease (Fisel et al., 2018; Liao et al., 2016) and members of both families are regulated in response to *M. refringens*. Changes in SLC expression are consistent with a disruption of many processes that depend on it and would explain the cascade of associated regulations. Strong alterations of these processes could conduct to the decrease in the condition index, loss of glycogen, discoloration of the digestive gland, growth arrest, tissue necrosis, and finally death of *M. refringens*-infected oysters. SLC genes have also been shown to be under directional selection (Li et al., 2021), suggesting their role in adaptive response to environmental parameters (Kenkel et al., 2013; Martinez Barrio et al., 2016; Zhou et al., 2018). *M. refringens* seems to act by inhibiting certain processes normally associated with an activation of the immune response. This is the case of the repression of the ADAM proteins which belong to the matrix-metalloproteinases family. They allow the detachment of adhesion proteins which play a central role in attracting inflammatory cells to sites of tissue damage and regulate the activation and phagocytic activity of immune cells (Rahn & Becker-Pauly, 2022). Some ADAM genes might have a protective role against some parasites (Geurts et al., 2012). We observed an upregulation of several GTPase IMAP family members in response to *M. refringens*. Those genes have been shown to be upregulated in *Oncorhynchus mykiss* in response to the myxozoan *Ceratonova shasta* (Barrett & Bartholomew, 2021) or in the mice in response to *Toxoplasma gondii* infection (Kim et al., 2018). If the way these genes might mediate resistance is unclear, one possible mechanism may be through mediating the effects of interferons (IFN pathway), which orchestrates many cellular pathways and regulates the expression of hundreds of genes. In *O. edulis*, associated with the over-expression of IMAP genes, we observed the over-expression of several interferons induced proteins which are known to act as a primary way of defending against some parasite infections (Baerwald, 2013). In the same order, many proteins exhibiting a coiled-coil domain are upregulated in response to *M. refringens* suggesting the presence of a strong cytoskeletal remodeling related to the presence of the parasite in the tissue. The coiled-coil domain-containing (CCDC) proteins have been implicated in many physiological processes including interactions with molecular components of signaling pathways and cytoskeletal
polymerization, and they are implicated in the pathogenesis of a large number of cancers (Priyanka & Yenugu, 2021). Some parasites appear to use the low-density lipoprotein receptor (LDLR) pathway of their host to facilitate their development and invasion of host tissue (Nagajyothi et al., 2011). The low-density lipoprotein receptor (LDLR) family is composed of a class of transmembrane glycoproteins that bind and internalize extracellular ligands for degradation by lysosomes. Their downregulation in infected oysters may reflect an alteration in the functioning of lysosomes, limiting their intervention in the digestion of parasites that could be endocytosed before their degradation. Interestingly, we also observed a downregulation of seven genes encoding proteins associated with lysosomal membrane or enzymes such as acid phosphatases. A large set of transmembrane proteins (TMEM) are downregulated in response to *M. refringens*. TMEM acts as channels to allow the transport of specific substances across the biological membranes. Their specific functions remain poorly described, but some of them are involved in modulating cancer cell dissemination and metastasis formation (Marx et al., 2020; Qiao et al., 2016). No data on their role in response to parasites are available in molluscs but with 170 putative genes detected in *O. edulis* genome, the understanding of their role in response to parasites will require additional studies. Complementary analyses involving a larger number of individuals could be undertaken to better understand the entire process of response to this parasite and to identify putative genetic loci associated with *M. refringens* resistance/tolerance particularly with the objective of re-establishing the culture of *O. edulis* on the foreshore.

## 5 | CONCLUSION

In this study, we present a high-quality genome assembly for the flat oyster *O. edulis* as a tool for further genomic studies. The analysis of tRNA, miRNA, and TEs revealed some specificities in this genome. The preliminary transcriptomic analysis made on the response of *O. edulis* to *M. refringens*, opens new research opportunities to better understand the molecular mechanisms involved in the oyster-parasite interactions, which are not currently understood despite the highly lethal nature of this parasite. Our study provides also a valuable reference genome for further comparative genomics and population genetic analysis or to detect the signature of local adaptation and genomic re-arrangement such as chromosome inversion as observed in several marine taxa. This genome provides a useful resource for genome-wide studies of production traits and genomic selection which may be essential for the sustainable development of flat oyster aquaculture and restoration.

## ACKNOWLEDGMENTS

This study was funded by the Roscoff_O.edulis-V1 Fonds Européen pour les Affaires Maritimes et la Pêche (FEAMP) (grant code: PFEA470017FA1000016). We thank Génome Québec Centre d'Expertise et de Services for performing the genome sequencing and providing associated advice leading up to the work. We

thank Dominique Marie for his help in flux cytometry, the Genomer Platform of Roscoff for technical support, and Yvon Madec and Stéphane Pouvreau for providing oyster populations.

## CONFLICT OF INTEREST

The authors declare no conflicts of interest.

## DATA AVAILABILITY STATEMENT

Raw sequence reads can be downloaded at the National Center for Biotechnology Information (NCBI) under BioProject PRJNA772088. The Roscoff_O.edulis-V1 genome is available under the accession JAMSGK010000000.

## ORCID

Ricardo Gonzales-Araya https://orcid.org/0000-0002-6518-4096

François H. Lallier https://orcid.org/0000-0001-8979-8917

Jean-François Flot https://orcid.org/0000-0003-4091-7916

Jakob Hemmer-Hansen https://orcid.org/0000-0002-1033-2565

Romain Koszul https://orcid.org/0000-0002-3086-1173

Arnaud Tanguy https://orcid.org/0000-0002-9670-2693

## REFERENCES

Admire, A., Shanks, L., Danzl, N., Wang, M., Weier, U., Stevens, W., Hunt, E., & Weinert, T. (2006). Cycles of chromosome instability are associated with a fragile site and are increased by defects in DNA replication and checkpoint controls in yeast. *Genes and Development*, 20(2), 159-173. https://doi.org/10.1101/gad.1392506

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., & Lipman, D. J. (1997). Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. *Nucleic Acids Research*, 25(17), 3389-3402. https://doi.org/10.1093/nar/25.17.3389

Baerwald, M. R. (2013). Temporal expression patterns of rainbow trout immune-related genes in response to *Myxobolus cerebralis* exposure. *Fish Shellfish Immunology*, 35, 965-971. https://doi.org/10.1016/j.fsi.2013.07.008

Bai, C. M., Xin, L. S., Rosani, U., Wu, B., Wang, Q. C., Duan, X. K., Liu, Z. H., & Wang, C. M. (2019). Chromosomal-level assembly of the blood clam, Scapharca (Anadara) broughtonii, using long sequence reads and Hi-C. *Gigascience*, 8(7), giz067. https://doi.org/10.1093/gigascience/giz067

Barrett, D. E., & Bartholomew, J. L. (2021). A tale of two fish: Comparative transcriptomics of resistant and susceptible steelhead following exposure to *Ceratonova shasta* highlights differences in parasite recognition. *PLoS One*, 16(2), e0234837.

Baudry, L., Guiglielmoni, N., Marie-Nelly, H., Cormier, A., Marbouty, M., Avia, K., Mie, Y. L., Godfroy, O., Sterck, L., Cock, J. M., Zimmer, C., Coelho, S. M., & Koszul, R. (2020). instaGRAAL: Chromosome-level quality scaffolding of genomes using a proximity ligation-based scaffolder. *Genome Biology*, 21(1), 148. https://doi.org/10.1186/s13059-020-02041-z

Beck, M. W., Brumbaugh, R. D., Airoldi, L., Carranza, A., Coen, L. D., Crawford, C., Defeo, O., Edgar, G. J., Hancock, B., Kay, M. C., Lenihan, H. S., Luckenbach, M. W., Toropova, C. L., Zhang, G., & Guo, X. (2011). Oyster reefs at risk and recommendations for conservation, restoration, and management. *Bioscience*, 61, 107-116.

Bermudez-Santana, C., Stephan-Otto Attolini, C., Kirsten, T., Engelhardt, J., Prohaska, S. J., Steigele, S., & Stadler, P. F. (2010). RGenomic organization of eukaryotic tRNAs. *BMC Genomics*, 11, 270.
Bernatchez, L., Wellenreuther, M., Araneda, C., Ashton, D. T., Barth, J. M. I., Beacham, T. D., & Withler, R. E. (2017). Harnessing the power of genomics to secure the future of seafood. *Trends in Ecology & Evolution*, 32(9), 665-680. https://doi.org/10.1016/j.tree.2017.06.010

Bezares-Calderon, L. A., Berger, J., & Jekely, G. (2020). Diversity of cilia-based mechanosensory systems and their functions in marine animal behaviour. *Philosophical Transactions Royal Society B*, 375, 20190376. https://doi.org/10.1098/rstb.2019.0376

Biemont, C., & Vieira, C. (2006). Genetics: Junk DNA as an evolutionary force. *Nature*, 443(7111), 521-524. https://doi.org/10.1038/443521a

Bolger, A. M., Lohse, M., & Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics*, 30, 2114-2120. https://doi.org/10.1093/bioinformatics/btu170

Bueno, M. J., & Malumbres, M. (2011). MicroRNAs and the cell cycle. *Biochimica et Biophysica Acta*, 1812(5), 592-601. https://doi.org/10.1016/j.bbadis.2011.02.002

Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., & Madden, T. L. (2009). BLAST+: Architecture and applications. *BMC Bioinformatics*, 10(1), 1-9. https://doi.org/10.1186/1471-2105-10-421

Capy, P., Gasperi, G., Biemont, C., & Bazin, C. (2000). Stress and transposable elements: Co-evolution or useful parasites? *Heredity*, 85(Pt 2), 101-106. https://doi.org/10.1046/j.1365-2540.2000.00751.x

Chan, J., Wang, L., Li, L., Mu, K., Bushek, D., Xu, Y., Guo, X., Zhang, G., & Zhang, L. (2021). Transcriptomic response to *Perkinsus marinus* in two *Crassostrea* oysters reveals evolutionary dynamics of host-parasite interactions. *Frontiers in Genetics*, 12, 795706. https://doi.org/10.3389/fgene.2021.795706

Chen, G., Zhang, C., Jiang, F., Wang, Y., Xu, Z., & Wang, C. (2014). Bioinformatics analysis of hemocyte miRNAs of scallop *Chlamys farreri* against acute viral necrobiotic virus (AVNV). *Fish and Shellfish Immunology*, 37(1), 75-86. https://doi.org/10.1016/j.fsi.2014.01.002

Chu, Q., Sun, Y., Cui, J., & Xu, T. (2017). MicroRNA-3570 modulates the NF-κB Pathway in Teleost Fish by Targeting MyD88. *Journal Immunology*, 198(8), 3274-3282. https://doi.org/10.4049/jimmunol.1602064

Cocci, P., Roncarati, A., Capriotti, M., Mosconi, G., & Palermo, F. A. (2020). Transcriptional alteration of gene biomarkers in hemocytes of wild *Ostrea edulis* with molecular evidence of infections with *Bonamia spp.* and/or *Marteilia refringens* parasites. *Pathogens*, 9(5), 323. https://doi.org/10.3390/pathogens9050323

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., Batut, P., Chaisson, M., & Gingeras, T. R. (2013). STAR: Ultrafast universal RNA-seq aligner. *Bioinformatics*, 29(1), 15-21. https://doi.org/10.1093/bioinformatics/bts635

du, X., Fan, G., Jiao, Y., Zhang, H., Guo, X., Huang, R., Zheng, Z., Bian, C., Deng, Y., Wang, Q., Wang, Z., Liang, X., Liang, H., Shi, C., Zhao, X., Sun, F., Hao, R., Bai, J., Liu, J., ... Liu, X. (2017). The pearl oyster *Pinctada fucata martensii* genome and multi-omic analyses provide insights into biomineralization. *Gigascience*, 6, 1-12. https://doi.org/10.1093/gigascience/gix059

Emms, D. M., & Kelly, S. (2019). OrthoFinder: Phylogenetic orthology inference for comparative genomics. *Genome Biology*, 20(1), 1-14. https://doi.org/10.1186/s13059-019-1832-y

Fang, Y., Xu, X. Y., Shen, Y., & Li, J. (2020). miR-148 targets CiGadd45ba and CiGadd45bb to modulate the inflammatory response to bacterial infection in grass carp. *Developmental Comparative Immunology*, 106, 103611. https://doi.org/10.1016/j.dci.2020.103611

Farhat, S., Bonnivard, E., Pales Espinosa, E., Tanguy, A., Boutet, I., Guiglielmoni, N., Flot, J. F., & Allam, B. (2022). Comparative analysis of the *Mercenaria mercenaria* genome provides insights into the diversity of transposable elements and immune molecules in bivalve molluscs. *BMC Genomics*, 23(1), 192. https://doi.org/10.1186/s12864-021-08262-1

Finnegan, D. J. (2012). Retrotransposons. *Current Biology*, 22(11), R432-R437. https://doi.org/10.1016/j.cub.2012.04.025

Fisel, P., Schaeffeffe, E., & Schwab, M. (2018). Clinical and functional relevance of the monocarboxylate transporter family in disease pathology and drug therapy. *Clinical and Translational Science*, 11, 352-364. https://doi.org/10.1111/cts.12551

Flynn, J. M., Hubley, R., Goubert, C., Rosen, J., Clark, A. G., Feschotte, C., & Smit, A. F. (2020). RepeatModeler2 for automated genomic discovery of transposable element families. *Proceedings of the National Academy of Sciences of the United States of America*, 117(17), 9451-9457. https://doi.org/10.1073/pnas.1921046117

Francklyn, C. S., & Minajigi, A. (2010). tRNA as an active chemical scaffold for diverse chemical transformations. *FEBS Letters*, 584, 366-375. https://doi.org/10.1016/j.febsslet.2009.11.045

Gerdol, M., De Moro, G., Venier, P., & Pallavicini, A. (2015). Analysis of synonymous codon usage patterns in sixty-four different bivalve species. *PeerJ*, 3, e1520. https://doi.org/10.7717/peerj.1520

Gerdol, M., Greco, S., & Pallavicini, A. (2019). Extensive tandem duplication events drive the expansion of the C1q-domain-containing gene family in bivalves. *Marine Drugs*, 17(10), 583. https://doi.org/10.3390/md17100583

Gervais, O., Chollet, B., Dubreuil, C., Durante, S., Feng, C., Hénard, C., Lecadet, C., Serpin, D., Tristan, R., & Arzul, I. (2019). Involvement of apoptosis in the dialogue between the parasite *Bonamia ostreae* and the flat oyster *Ostrea edulis*. *Fish and Shellfish Immunology*, 93, 958-964. https://doi.org/10.1016/j.fsi.2019.08.035

Geurts, N., Opdenakker, G., & Van den Steen, P. E. (2012). Matrix metalloproteinases as therapeutic targets in protozoan parasitic infections. *Pharmacology and Therapeutics*, 125:279, 257-279. https://doi.org/10.1016/j.pharmthera.2011.11.008

Goodenbour, J. M., & Pan, T. (2006). Diversity of tRNA genes in eukaryotes. *Nucleic Acids Research*, 34(21), 6137-6146. https://doi.org/10.1093/nar/gkl725

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., Adiconis, X., Fan, L., Raychowdhury, R., Zeng, Q., Chen, Z., Mauceli, E., Hacohen, N., Gnrake, A., Rhind, N., di Palma, F., Birren, B. W., Nusbaum, C., Lindblad-Toh, K., ... Regev, A. (2011). Full-length transcriptome assembly from RNA-seq data without a reference genome. *Nature Biotechnology*, 29(7), 644-652. https://doi.org/10.1038/nbt.1883

Griffiths-Jones, S. (2005). Annotating non-coding RNAs with Rfam. *Current Protocols Bioinformatics*, 12, 12.5.1-12.5.12. https://doi.org/10.1002/0471250953.bi1205s9

Gundappa, M. K., Peñaloza, C., Regan, T., Boutet, I., Tanguy, A., Houston, R. D., Bean, T. P., & Macqueen, D. J. (2022). Chromosome level reference genome for European flat oyster (*Ostrea edulis L*). *bioRxiv*. https://doi.org/10.1101/2022.06.26.497633

Guo, X. (2021). Genetics in shellfish culture. In S. Shumway (Ed.), *Molluscan shellfish aquaculture: A practical guide* (pp. 393-413). 5m Books Ltd.

Haas, B. J., Delcher, A. L., Mount, S. M., Wortman, J. R., Smith, R. K. Jr., Hannick, L. I., Maiti, R., Ronning, C. M., Rusch, D. B., Town, C. D., Salzberg, S. L., & White, O. (2003). Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. *Nucleic Acids Research*, 31(19), 5654-5666. https://doi.org/10.1093/nar/gkg770

He, Y., Jouaux, A., Ford, S. E., Lelong, C., Sourdaine, P., Mathieu, M., & Guo, X. (2015). Transcriptome analysis reveals strong and complex antiviral response in a mollusc. *Fish and Shellfish Immunology*, 46, 131-144. https://doi.org/10.1016/j.fsi.2015.05.023

He, Y., Yu, H., Bao, Z., Zhang, Q., & Guo, X. (2012). Mutation in promoter region of a serine protease inhibitor confers *Perkinsus marinus* resistance in the eastern oyster (*Crassostrea virginica*). *Fish and Shellfish Immunology*, 33(2), 411-417.

Hollenbeck, C. M., & Johnston, I. A. (2018). Genomic tools and selective breeding in Molluscs. *Frontiers in Genetics*, 9, 253. https://doi.org/10.3389/fgene.2018.00253
Huang, D. W., Sherman, B. T., & Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. *Nature Protocols*, 4(1), 44–57.

Huang, S., Asaduzzaman, M., Yoshitake, K., & Kinoshita, S. (2021). Discovery and functional understanding of miRNAs in molluscs: A genome wide profiling approach. *RNA Biology*, 18, 1702–1715. https://doi.org/10.1080/15476286.2020.1867798

Huang, S., Kang, M., & Xu, A. (2017). HaploMerger2: Rebuilding both haploid sub-assemblies from high-heterozygosity diploid genome assembly. *Bioinformatics*, 33(16), 2577–2579. https://doi.org/10.1093/bioinformatics/btx220

Insua, A., & Mendez, J. (1998). Physical mapping and activity of ribosomal RNA genes in mussel *Mytilus Galloprovincialis*. *Hereditas*, 128(3), 189–194.

Jurka, J., Kapitonov, V. V., Pavlicek, A., Klonowski, P., Kohany, O., & Walichiewicz, J. (2005). Repbase update, a database of eukaryotic repetitive elements. *Cytogenetic and Genome Research*, 110(1–4), 462–467. https://doi.org/10.1159/000084979

Katoh, K., Rozewicki, J., & Yamada, K. D. (2018). MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. *Briefing in Bioinformatics*, 20(4), 1160–1166. https://doi.org/10.1093/bib/bbx108

Kenkel, C. D., Meyer, E., & Matz, M. V. (2013). Gene expression under chronic heat stress in populations of the mustard hill coral (*Porites astreoides*) from different thermal environments. *Molecular Ecology*, 22, 4322–4334. https://doi.org/10.1111/mec.12390

Kenny, N. J., McCarthy, S. A., Dudchenko, O., James, K., Betteridge, E., Corton, C., Dolucan, J., Mead, D., Oliver, K., Omer, A. D., Pelan, S., Ryan, Y., Sims, Y., Skelton, J., Smith, M., Torrance, J., Weisz, D., Wipat, A., Aiden, E. L., ... Williams, S. T. (2020). The gene-rich genome of the scallop *Pecten maximus*. *Gigascience*, 9(5), giaa037. https://doi.org/10.1093/gigascience/giaa037

Kim, C. Y., Zhang, X., & Witola, W. H. (2018). Small GTPase immunity-associated proteins mediate resistance to *Toxoplasma gondii* infection in Lewis rat. *Infection Immunity*, 86(4), e00582-17. https://doi.org/10.1128/IAI.00582-17

Kim, D., Paggi, J. M., Park, C., Bennett, C., & Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. *Nature Biotechnology*, 37(8), 907–915. https://doi.org/10.1038/s41587-019-0201-4

Kiss, T. (2010). Apoptosis and its functional significance in molluscs. *Apoptosis*, 15, 313–321. https://doi.org/10.1007/s10495-009-0446-3

Lallias, D., Beaumont, A. R., Haley, C. S., Boudry, P., Heurtebise, S., & Lapègue, S. (2007). A first-generation genetic linkage map of the European flat oyster *Ostrea edulis* (L.) based on AFLP and microsatellite markers. *Animal Genetics*, 38(6), 560–568. https://doi.org/10.1111/j.1365-2052.2007.01647.x

Langmead, B., Trapnell, C., Pop, M., & Salzberg, S. L. (2019). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. *Genome Biology*, 30(3), R25. https://doi.org/10.1186/gb-2009-10-3-r25

Lau, Y. T., Santos, B., Barbosa, M., Pales Espinosa, E., & Allam, B. (2018). Regulation of apoptosis-related genes during interactions between oyster hemocytes and the alveolate parasite *Perkinsus marinus*. *Fish and Shellfish Immunology*, 83, 180–189. https://doi.org/10.1016/j.fsi.2018.09.006

Lazar-Stefanita, L., Scolari, V. F., Mercy, G., Muller, H., Guérin, T. M., Thierry, A., Mozziconacci, J., & Koszul, R. (2017). Cohesins and condensins orchestrate the 4D dynamics of yeast chromosomes during the cell cycle. *EMBO Journal*, 36(18), 2684–2697. https://doi.org/10.15252/embj.201797342

Le Roux, F., Audemard, C., Barnaud, A., & Berthe, F. (1999). DNA probes as potential tools for the detection of *Marteilia refringens*. *Marine Biotechnology*, 1(6), 588–597. https://doi.org/10.1007/pl00011814

Li, A., Dai, H., Guo, X., Zhang, Z., Zhang, K., Wang, C., Wang, X., Wang, W., Chen, H., Li, X., Zheng, H., Li, L., & Zhang, G. (2021). Genome of the estuarine oyster provides insights into climate impact and adaptive plasticity. *Communications Biology*, 4, 1287. https://doi.org/10.1038/s42003-021-02823-6

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G., Durbin, R., & 1000 Genome Project Data Processing Subgroup. (2009). The sequence alignment/map format and SAMtools. *Bioinformatics*, 25(16), 2078–2079. https://doi.org/10.1093/bioinformatics/btp352

Li, Y., Niu, D., Bai, Y., Lan, T., Peng, M., Dong, Z., & Li, J. (2019). Characterization of the ScghC1q-1 gene in *Sinonovacula constricta* and its role in innate immune responses. *Developmental and Comparative Immunology*, 94, 16–21. https://doi.org/10.1016/j.dci.2019.01.004

Liao, Q. D., Wang, C. G., Zhu, Y. D., Chen, W. H., Shao, S. L., Jiang, F. N., & Xu, X. M. (2016). Decreased expression of SLC39A14 is associated with tumor aggressiveness and biochemical recurrence of human prostate cancer. *Oncotargets and Therapy*, 9, 4197–4205. https://doi.org/10.2147/OTT.S103640

Lieberman-Aiden, E., van Berkum, N. L., Williams, L., Imakaev, M., Ragoczy, T., Telling, A., Amit, I., Lajoie, B. R., Sabo, P. J., Dorschner, M. O., Sandstrom, R., Bernstein, B., Bender, M. A., Groudine, M., G nirke, A., Stamatoyannopoulos, J., Mirny, L. A., Lander, E. S., & Dekker, J. (2009). Comprehensive mapping of long-range inter-actions reveals folding principles of the human genome. *Science*, 326(5950), 289–293. https://doi.org/10.1126/science.1181369

Llorens, C., Futami, R., Covelli, L., Domínguez-Escribá, L., Viu, J. M., Tamarit, D., Aguilar-Rodríguez, J., Vicente-Ripolles, M., Fuster, G., Bernet, G. P., Maumus, F., Munoz-Pomer, A., Sempere, J. M., Latorre, A., & Moya, A. (2011). The gypsy database (GyDB) of Mobile genetic elements: Release 2.0. *Nucleic Acids Research*, 39(Suppl 1), D70–D74. https://doi.org/10.1093/nar/gkq1061

Love, M. I., Huber, W., & Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biology*, 15(12), 550. https://doi.org/10.1186/s13059-014-0550-8

Lowe, T. M., & Eddy, S. R. (1997). tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. *Nucleic Acids Research*, 25(5), 955–964. https://doi.org/10.1093/nar/25.5.955

Martinez Barrio, A., Lamichhaney, S., Fan, G., Rafati, N., Pettersson, M., Zhang, H., Dainat, J., Ekman, D., Höppner, M., Jern, P., Martin, M., Nystedt, B., Liu, X., Chen, W., Liang, X., Shi, C., Fu, Y., Ma, K., Zhan, X., ... Andersson, L. (2016). The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing. *eLife*, 5, e12081. https://doi.org/10.7554/eLife.12081.001

Martín-Gómez, L., Villalba, A., Kerkhoven, R. H., & Abollo, E. (2014). Role of microRNAs in the immunity process of the flat oyster *Ostrea edulis* against bonamiosis. *Infection Genetics Evolution*, 27, 40–50. https://doi.org/10.1016/j.meegid.2014.06.026

Marx, S., Dal Maso, T., Chen, J. W., Bury, M., Wouters, J., Michiels, C., & Le Calvé, B. (2020). Transmembrane (TMEM) protein family members: Poorly characterized even if essential for the metastatic process. *Seminars in Cancer Biology*, 60, 96–106. https://doi.org/10.1016/j.semcancer.2019.08.018

McFarlane, R. J., & Whitehall, S. K. (2009). tRNA genes in eukaryotic genome organization and reorganization. *Cell Cycle*, 8(19), 3102–3106. https://doi.org/10.4161/cc.8.19.9625

McMaster, W. R., Morrison, C. J., & Kobor, M. S. (2016). Epigenetics: A new model for intracellular parasite-host cell regulation. *Trends Parasitology*, 32(7), 515–521. https://doi.org/10.1016/j.pt.2016.04.002

Melayah, D., Bonnivard, E., Chalhoub, B., Audeon, C., & Grandbastien, M. A. (2001). The mobility of the tobacco Tnt1 retrotransposon correlates with its transcriptional activation by functional factors. *Plant Journal*, 28(2), 159–168. https://doi.org/10.1046/j.1365-313x.2001.01141.x
Metzger, M. J., Paynter, A. N., Siddall, M. E., & Goff, S. P. (2018). Horizontal transfer of retrotransposons between bivalves and other aquatic species of multiple phyla. *Proceedings of the National Academy of Sciences of the United States of America*, 115(18), E4227-E4235. https://doi.org/10.1073/pnas.17172

Michaud, M., Cognat, V., Duchêne, A. M., & Maréchal-Drouard, L. (2011). A global picture of tRNA genes in plant genomes. *Plant Journal*, 66(1), 80-93. https://doi.org/10.1111/j.1365-313X.2011.04490.x

Milne, I., Lindner, D., Bayer, M., Husmeier, D., Mcguire, G., Marshall, D. F., & Wright, F. (2009). TOPALi v2: A rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops. *Bioinformatics*, 25(1), 126-127. https://doi.org/10.1093/bioinformatics/btn575

Modak, T. H., Literman, R., Puritz, J. B., Johnson, K. M., Roberts, E. M., Proestou, D., Guo, X., Gomez-Chiarri, M., & Schwartz, R. (2021). Extensive genome-wide duplications in the eastern oyster (*Crassostrea virginica*). *Philosophical Transactions of the Royal Society B*, 376(1825), 20200164.

Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., & Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-seq. *Nature Methods*, 5, 621-628. https://doi.org/10.1038/nmeth.1226

Moyo, L., Ramesh, S. V., Kappagantam, M., Mitter, N., Sathuvalli, V., & Pappu, H. R. (2017). The effects of potato virus Y-derived virus small interfering RNAs of three biologically distinct strains on potato (*Solanum tuberosum*) transcriptome. *Virology Journal*, 14(1), 129. https://doi.org/10.1186/s12985-017-0803-8

Mun, S., Kim, Y. J., Markkandan, K., Shin, W., Oh, S., Woo, J., Yoo, J., An, H., & Han, K. (2017). The whole-genome and transcriptome of the Manila clam (*Ruditapes philippinarum*). *Genome Biology and Evolution*, 9(6), 1487-1498. https://doi.org/10.1093/gbe/evx096

Nagajyothi, F., Weiss, L. M., Silver, D. L., Desruisseaux, M. S., Scherer, P. E., Herz, J., & Tanowitz, H. B. (2011). *Trypanosoma cruzi* utilizes the host low density lipoprotein receptor in invasion. *PLoS Neglected Tropical Diseases*, 5(2), e953. https://doi.org/10.1371/journal.pntd.0000953

Nawrocki, E. P., & Eddy, S. R. (2013). Infernal 1.1: 100-fold faster RNA homology searches. *Bioinformatics*, 29(22), 2933-2935. https://doi.org/10.1093/bioinformatics/btt509

Nischal, L., Mohsin, M., Khan, I., Kardam, H., Wadhwa, A., Abrol, Y. P., Iqbal, M., & Ahmad, A. (2014). Identification and comparative analysis of microRNAs associated with low-N tolerance in rice genotypes. *PLoS One*, 7(12), 1-13. https://doi.org/10.1371/journal.pone.0050261

OSPAR. (2009). Background document for Ostrea edulis and Ostrea edulis beds. In Convention, O. (Ed.), https://www.ospar.org/documents?d=7183

Pardo, B. G., Álvarez-Dios, J. A., Cao, A., Ramilo, A., Gómez-Tato, A., Planas, J. V., Villalba, A., & Martínez, P. (2016). Construction of an *Ostrea edulis* database from genomic and expressed sequence tags (ESTs) obtained from *Bonamia ostreae* infected haemocytes: Development of an immune-enriched oligo-microarray. *Fish and Shellfish Immunology*, 59, 331-344. https://doi.org/10.1016/j.fsi.2016.10.047

Peñaloza, C., Gutierrez, A. P., Eöry, L., Wang, S., Guo, X., Archibald, A. L., Bean, T. P., & Houston, R. D. (2021). A chromosome-level genome assembly for the Pacific oyster *Crassostrea gigas*. *GigaScience*, 10(3), giab020. https://doi.org/10.1093/gigascience/giab020

Peng, J., Li, Q., Xu, L., Wei, P., He, P., Zhang, X., Zhang, L., Guan, J., Zhang, X., Lin, Y., Gui, J., & Chen, X. (2020). Chromosome-level analysis of the *Crassostrea hongkongensis* genome reveals extensive duplication of immune-related genes in bivalves. *Molecular Ecology Resources*, 20(4), 980-994. https://doi.org/10.1111/1755-0998.13157

Potts, R. W. A., Gutierrez, A. P., Penaloza, C. S., Regan, T., Bean, T. P., & Houston, R. D. (2021). Potential of genomic technologies to improve disease resistance in molluscan aquaculture. *Philosophical

Transactions of the Royal Society B, 376, B3762020016820200168
10.1098/rstb.2020.0168.

Powell, D., Subramanian, S., Suwansa-Ard, S., Zhao, M., O'Connor, W., Raftos, D., & Elizur, A. (2018). The genome of the oyster Saccostrea offers insight into the environmental resilience of bivalves. *DNA Research*, 25(6), 655-665.

Priyanka, P., & Yenugu, S. (2021). Coiled-coil domain-containing (CCDC) proteins: Functional roles in general and male reproductive physiology. *Reproductive Sciences*, 28(10), 2725-2734. https://doi.org/10.1007/s43032-021-00595-2

Qiao, W., Han, Y., Jin, W., Tian, M., Chen, P., Min, J., Hu, H., Xu, B., Zhu, W., Xiong, L., & Lin, Q. (2016). Overexpression and biological function of TMEM48 in non-small cell lung carcinoma. *Tumour Biology*, 37, 2575-2586. https://doi.org/10.1007/s13277-015-4014-2

R Core Team. (2017). R: A language and environment for statistical computing. R Foundation for Statistical Computing. https://www.R-project.org/

Rahn, S., & Becker-Pauly, C. (2022). Meprin and ADAM proteases as triggers of systemic inflammation in sepsis. *FEBS Letters*, 596(5), 534-556. https://doi.org/10.1002/1873-3468.14225

Rives, M. L., Javitch, J. A., & Wickenden, A. D. (2017). Potentiating SLC transporter activity: Emerging drug discovery opportunities. *Biochemical Pharmacology*, 135, 1-11. https://doi.org/10.1016/j.bcp.2017.02.010

Rocha, E. P. C. (2004). Codon usage bias from tRNA's point of view: Redundancy, specialization, and efficient decoding for translation optimization. *Genome Research*, 14, 2279-2286.

Ronza, P., Cao, A., Robledo, D., Gómez-Tato, A., Álvarez-Dios, J. A., Hasanuzzaman, A. F. M., Quiroga, M. I., Villalba, A., Pardo, B. G., & Martínez, P. (2018). Long-term affected flat oyster (*Ostrea edulis*) haemocytes show differential gene expression profiles from naïve oysters in response to *Bonamia ostreae*. *Genomics*, 110(6), 390-398. https://doi.org/10.1016/j.ygeno.2018.04.002

Rosani, U., Abbadi, M., Green, T. J., & Bai, C. M. (2020). Parallel analysis of miRNAs and mRNAs suggests distinct regulatory networks in *Crassostrea gigas* infected by Ostreid herpesvirus 1. *BMC Genomics*, 21(1), 620. https://doi.org/10.1186/s12864-020-07026-7

Rosani, U., Borteletto, E., Bai, C. M., Novoa, B., Figueras, A., Venier, P., & Fromm, B. (2021). Digging into bivalve miRNAomes: Between conservation and innovation. *Philosophical Transactions of the Royal Society B*, 376, 20200165.

Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., & Zdobnov, E. M. (2015). BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. *Bioinformatics*, 31(19), 3210-3212. https://doi.org/10.1093/bioinformatics/btv351

Smaal, A. C., Kamermans, P., van der Have, T. M., Engelsma, M., & Sas, H. J. W. (2015). Feasibility of flat oyster (*Ostrea edulis* L.). *Restoration in the Dutch part of the North Sea* (pp. 1-58). Institute for Marine Resources & Ecosystem Studies.

Song, H., Guo, X., Sun, L., Wang, Q., Han, F., Wang, H., Wray, G. A., Davidson, P., Wang, Q., Hu, Z., Zhou, C., Yu, Z., Yang, M., Feng, J., Shi, P., Zhou, Y., Zhang, L., & Zhang, T. (2021). The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia. *BMC Biology*, 19(1), 15. https://doi.org/10.1186/s12915-020-00943-9

Song, L., Wang, L., Zhang, H., & Wang, M. (2015). The immune system and its modulation mechanism in scallop. *Fish Shellfish Immunology*, 46, 65-78. https://doi.org/10.1016/j.fsi.2015.03.013

Stanke, M., Schöffmann, O., Morgenstern, B., & Waack, S. (2006). Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. *BMC Bioinformatics*, 7(1), 62. https://doi.org/10.1186/1471-2105-7-62

Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., & Kumar, S. (2011). MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum
parsimony methods. *Molecular Biology and Evolution*, 28(10), 2731-2739. https://doi.org/10.1093/molbev/msr121

Thiriot-Quiévreux, C., & Ayraud, N. (1982). Les caryotypes de quelques espèces de bivalves et de gastéropodes marins. *Marine Biology*, 70, 165-172. https://doi.org/10.1007/BF00397681

Thomas-Bulle, C., Piednoël, M., Donnart, T., Filée, J., Jollivet, D., & Bonnivard, E. (2018). Mollusc genomes reveal variability in patterns of LTR-retrotransposons dynamics. *BMC Genomics*, 19(1), 1-18. https://doi.org/10.1186/s12864-018-5200-1

Vera, M., Pardo, B. G., Cao, A., Vilas, R., Fernández, C., Blanco, A., Gutierrez, A. P., Bean, T. P., Houston, R. D., Villalba, A., & Martínez, P. (2019). Signatures of selection for bonamiosis resistance in European flat oyster (*Ostrea edulis*): New genomic tools for breeding programs and management of natural resources. *Evolutionary Applications*, 12(9), 1781-1796. https://doi.org/10.1111/eva.12832

Wang, S., Zhang, J., Jiao, W., Li, J., Xun, X., Sun, Y., Guo, X., Huan, P., Dong, B., Zhang, L., Hu, X., Sun, X., Wang, J., Zhao, C., Wang, Y., Wang, D., Huang, X., Wang, R., Lv, J., ... Bao, Z. (2017). Scallop genome provides insights into evolution of bilaterian karyotype and development. *Nature Ecology & Evolution*, 1(5), 1-12.

Witkop, E. M., Proestou, D. A., & Gomez-Chiarri, M. (2022). The expanded inhibitor of apoptosis gene family in oysters possesses novel domain architectures and may play diverse roles in apoptosis following immune challenge. *BMC Genomics*, 23, 201. https://doi.org/10.1186/s12864-021-08233-6

Xun, X., Cheng, J., Wang, J., Li, Y., Li, X., Li, M., Lou, J., Kong, Y., Bao, Z., & Hu, X. (2020). Solute carriers in scallop genome: Gene expansion and expression regulation after exposure to toxic dinoflagellate. *Chemosphere*, 241, 124968. https://doi.org/10.1016/j.chemosphere.2019.124968

Yang, J. L., Feng, D. D., Liu, J., Xu, J. K., Chen, K., Li, Y. F., Zhu, Y. T., Liang, X., & Lu, Y. (2021). Chromosome-level genome assembly of the hard-shelled mussel *Mytilus coruscus*, a widely distributed species from the temperate areas of East Asia. *GigaScience*, 10(4), giab024. https://doi.org/10.1093/gigascience/giab024

Yu, H., He, Y., Wang, X., Zhang, Q., Bao, Z., & Guo, X. (2011). Polymorphism in a serine protease inhibitor gene and its association with disease

resistance in the eastern oyster (*Crassostrea virginica Gmelin*). *Fish & Shellfish Immunology*, 30(3), 757-762.

Zhang, G., Fang, X., Guo, X., Li, L., Luo, R., Xu, F., Yang, P., Zhang, L., Wang, X., Qi, H., Xiong, Z., Que, H., Xie, Y., Holland, P. W., Paps, J., Zhu, Y., Wu, F., Chen, Y., Wang, J., ... Wang, J. (2012). The oyster genome reveals stress adaptation and complexity of shell formation. *Nature*, 490(7418), 49-54.

Zhou, X., Guang, X., Sun, D., Xu, S., Li, M., Seim, I., Jie, W., Yang, L., Zhu, Q., Xu, J., Gao, Q., Kaya, A., Dou, Q., Chen, B., Ren, W., Li, S., Zhou, K., Gladyshev, V. N., Nielsen, R., ... Yang, G. (2018). Population genomics of finless porpoises reveal an incipient cetacean species adapted to freshwater. *Nature Communications*, 9(1), 1276. https://doi.org/10.1038/s41467-018-03722-x

Zimin, A. V., Marçais, G., Puiu, D., Roberts, M., Salzberg, S. L., & Yorke, J. A. (2013). The MaSuRCA genome assembler. *Bioinformatics*, 29(21), 2669-2677. https://doi.org/10.1093/bioinformatics/btt476

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Supporting Information section at the end of this article.

**How to cite this article:** Boutet, I., Alves Monteiro, H. J., Baudry, L., Takeuchi, T., Bonnivard, E., Billoud, B., Farhat, S., Gonzales-Araya, R., Salaun, B., Andersen, A. C., Toullec, J.-Y., Lallier, F. H., Flot, J.-F., Guiglielmoni, N., Guo, X., Li, C., Allam, B., Pales-Espinosa, E., Hemmer-Hansen, J. ... Tanguy, A. (2022). Chromosomal assembly of the flat oyster (*Ostrea edulis* L.) genome as a new genetic resource for aquaculture. *Evolutionary Applications*, 15, 1730-1748. https://doi.org/10.1111/eva.13462