- Research article
- Open Access
A preliminary integrated genetic map distinguishes every chromosome pair and locates essential genes related to abiotic adaptation of Crassostrea angulata/gigas
BMC Geneticsvolume 19, Article number: 104 (2018)
The re-sequencing of C. angulata has revealed many polymorphisms in candidate genes related to adaptation to abiotic stress that are not present in C. gigas; these genes, therefore, are probably related to the ability of this oyster to retain high concentrations of toxic heavy metals. There is, in addition, an unresolved controversy as to whether or not C. angulata and C. gigas are the same species or subspecies. Both oysters have 20 metacentric chromosomes of similar size that are morphologically indistinguishable. From a genomic perspective, as a result of the great variation and selection for heterozygotes in C. gigas, the assembly of its draft genome was difficult: it is fragmented in more than seven thousand scaffolds.
In this work sixty BAC sequences of C. gigas downloaded from NCBI were assembled in BAC-contigs and assigned to BACs that were used as probes for mFISH in C. angulata and C. gigas. In addition, probes of H3, H4 histone, 18S and 5S rDNA genes were also used. Hence we obtained markers identifying 8 out the 10 chromosomes constituting the karyotype. Chromosomes 1 and 9 can be distinguished morphologically. The bioinformatic analysis carried out with the BAC-contigs annotated 88 genes. As a result, genes associated with abiotic adaptation, such as metallothioneins, have been positioned in the genome. The gene ontology analysis has also shown many molecular functions related to metal ion binding, a phenomenon associated with detoxification processes that are characteristic in oysters. Hence the provisional integrated map obtained in this study is a useful complementary tool for the study of oyster genomes.
In this study 8 out of 10 chromosome pairs of Crassostrea angulata/gigas were identified using BAC clones as probes. As a result all chromosomes can now be distinguished. Moreover, FISH showed that H3 and H4 co-localized in two pairs of chromosomes different that those previously escribed. 88 genes were annotated in the BAC-contigs most of them related with Molecular Functions of protein binding, related to the resistance of the species to abiotic stress. An integrated genetic map anchored to the genome has been obtained in which the BAC-contigs structure were not concordant with the gene structure of the C. gigas scaffolds displayed in the Genomicus database.
The worldwide production of oysters is 5.15 million tons/year, second only by weight to carp and other cyprinids ; the genera Ostrea and Crassostrea (Ostreidae) make up most of the species of commercial interest. The Pacific oyster Crassostrea gigas (Thunberg, 1793), in particular, is now one of the most important bivalves produced commercially throughout the world. Global annual production volume is uncertain, owing to taxonomic confusion in reports from China and other countries , but is conservatively estimated at 626,000 metric tons in 2016 . Another Crassostrea species of interest is the Portuguese oyster, C. angulata (Lamarck, 1919); this species is an edible cupped oyster of major commercial importance and wide geographic distribution. It is assumed to be native to the Northwest Pacific region and has been introduced into many countries around the world . C. angulata was a species of considerable commercial importance in Europe up to the early 1970’s, but serious mortalities in those years almost wiped out the Portuguese oyster from Europe. As a result, the extent of remaining C. angulata populations is poorly known. Recent studies have shown that the Fujian oyster that has been harvested and cultivated for more than 2000 years in coastal river mouths and estuaries of southern China is actually C. angulata. The production of this Fujian oyster accounts for about 50% of total oyster production in China. At the present time, according to some authors, C. angulata should be considered as a subspecies of C. gigas  . Because of the commercial importance of Crassostrea species, a number of breeding programs have been initiated over the years, and substantial genetic and genomic knowledge has been accumulated [4,5,6,7,8,9], although some important characteristics of genome organization remain unclear.
Oysters have a diploid number of 20 chromosomes . Their karyotypes show most of their chromosomes to be submetacentric or metacentric, that are of similar size and morphologically indistinguishable; molecular markers are needed to differentiate them, mainly because of the many instances of aneuploidy occurring in natural populations of mollusks . To study the evolution and phylogenetic relationships among oyster species, it is necessary to characterize the chromosomes individually, and particular regions of them. The C-banding and restriction techniques, as well as the study of Ag-NOR and cytogenetic markers, have been applied to several different species of oysters describing several cytogenetic markers that allow identify chromosome pairs in some cases [12,13,14,15,16,17].
The introduction of the fluorescence in situ hybridization (FISH) technique opened the possibility of studying more deeply inside the chromosomes than had been possible before. FISH is also a powerful tool for defining the cytogenetic location and relative order of DNA sequences, thereby anchoring the genome sequence to the chromosomes. It has been successfully used to integrate genetic and cytogenetic maps in many plants and animals [18,19,20,21,22]. During the past two decades, there have been many FISH studies, based on species of interest in aquaculture, that have contributed to a better understanding of these genomes. The FISH technique has greatly facilitated comparative genomic analysis to study the evolution of genomes [23,24,25,26,27]. In oysters, analysis has been focused on the study of repetitive and multigene families [28,29,30,31].
Histones are basic structural proteins of the nucleosomes and, in eukaryotic organisms, are essential for DNA packaging. This group is subdivided into five protein classes, termed H1, located bound to the “linker DNA” region between nucleosomes, and H2A, H2B, H3, and H4, located in the core of the nucleosomes. The histones have functional consequences, as trans-acting regulatory factors, because the transcription machinery must gain access to their specific binding sites at the DNA level . Thus, knowing its location can help to understand its chromosome organization and evolution. In addition, the histone genes have been used as markers for taxonomic differentiation [33,34,35,36].
Going beyond the controversy of whether C. angulata and C. gigas are different species or subspecies, the whole genome re-sequencing of C. angulata showed many polymorphisms in candidate genes related to adaptation to abiotic stress that are not present in C. gigas, and hence are probably related to the ability of this oyster to retain high concentrations of toxic heavy metals .
Genetic maps are essential tools, useful for mapping quantitative-trait loci in marker-assisted selection, positional cloning, and genome assembly. In oyster species, linkage and physical maps have been described [38, 39]. Although the publication of the C. gigas genome draft is considered a great step forward in the knowledge of these organisms, at the present time the genome is assembled only to a level below that of the scaffold (N50 = 31,239) due to the complexity of its genome, its variation and selection for heterozygotes ; currently more than seven thousand scaffolds have been found. Mapping using BAC-FISH probes has allowed progress, in a straightforward and easy way, to final assembly and mapping of genes in their chromosomes. However, cytogenetic mapping of clones with large inserts is at present limited to a few species of interest for aquaculture [18, 23, 26, 27]. Knowing from gene libraries the position of clones that contain large inserts presents important advantages, since it is not limited to the location of repetitive DNA and allows the localization of multiple probes for each chromosome, thereby raising the number of markers per chromosome karyotype and in each species. From these clones, the inserts can be sequenced with next generation sequencing techniques, thereby enabling us to know the position and insert genes in the library hybridized in a straightforward manner [23, 26]. Moreover, the application of FISH based on the use of bacterial artificial chromosome clones (BAC) as probes provides an efficient approach for anchoring linkage and genomic data in the physical chromosomes . In mollusks, genetic and cytogenetic maps have been integrated using this BAC-FISH approach in the Zhikong Scallop, Chlamys farreri, selecting markers from microsatellite-based linkage maps from a BAC library . In oysters, to date, this approach has been applied only in the American oyster C. virginica, using inserts from a library based on phage P1 .
In order to gain insight into the genomes and chromosomes of C. angulata/gigas, an integrated genetic map, based on BACs and multigene families, has been designed with the aim of providing chromosomal markers and that could be anchored reliably to the sequenced genome of the species.
Material and methods
Selection and BAC clones getting from a C. gigas genetic library
Sixty BAC sequences downloaded from the National Center for Biotechnology Information (NCBI: http://www.ncbi.nlm.nih.gov/) (Acc. numbers from GU207404 to GU207462, and GU324325) were assembled with SeqMan software (Lasergene, DNASTAR). These sequences correspond to different BAC clones from a genetic library of C. gigas belonging to the Clemson University Genomics Institute (CUGI). From the contigs obtained (from now on named BAC-contigs), one BAC per contig and a single BAC (not included in any contig) were selected and ordered from the CUGI (Additional file 1: Table S1). The BAC clones were then mapped on C. angulata and C. gigas chromosomes by means of 2- and multi-color FISH.
Probes and labelling for FISH
DNA was extracted from BACs using the BACMAXtm DNA purification kit (Epicentre), following manufacturer’s instructions. For 2-color FISH, probes were labeled using the DIG/BIO Nick Translation Mix kit (Roche Molecular Biochemicals). In order to label 5S rDNA and 18S rDNA multigene families, PCR amplifications with biotin-16-dUTP or digoxigenin-11-dUTP (Roche Molecular Biochemicals) were made . For the H3 histone gene the primers H3F (5’-CGTAAATCCACTGGAGGCAAGG-3′) and H3R (5’-GGATGGCGCACAGGTTGGTGTC-3′) were used . The H4 histone genes were amplified by primers H4-F (5’-TGAGAGATAACATCCAGGGTATCAC-3’) and H4-R (5′- CTCTTGAGGGCGTAGACAACGTCCAT-3′).
Chromosome preparations and FISH mapping
Samples of C. angulata of 2–3 cm in length were obtained from naturally-occurring populations located in the mouth of the River Guadalquivir, in Sanlúcar de Barrameda (Cádiz, Spain). C. gigas samples came from the facilities of a local aquaculture company. After feeding the oysters with phytoplankton for several days, they were placed in beakers and cultured for 8 h in colchicine. The gills were then extracted and put under hypotonic shock conditions by incubating with KCl 0.4% for at least 1 h. After that, they were transferred to fresh Carnoy for 1 h .
Twelve BACs were used as probes to hybridize on C. angulata and C. gigas chromosomes using the FISH-BAC technique (Table 1; Additional file 1: Table S1). In addition, probes of H3 and H4 histones, 18S and 5S rDNA genes were also used to analyze their position in chromosomes relative to the hybridized BACs.
The 2-color FISH method used for BAC clones probes was performed following the protocol of García-Cegarra et al. . Multiple-FISH (mFISH) was carried out as described in Portela-Bens et al. . The 2-color FISH technique described by Cross et al.  was used to localize the H3 and H4 histone genes.
Bioinformatic analysis of BAC clones
The contigs containing hybridized BACs were annotated from the C. gigas genome draft, available at the NCBI . The BAC-contig sequences were subjected to Blast search against C. gigas Whole-Genome shotgun contigs from NCBI/BLAST/blastn suite (from now on named Cg-contigs) (AFTI: Crassostrea gigas strain 05 × 7-T-G4–1.051#20, whole genome shotgun sequencing project). In order to obtain genes from these Cg-contigs, a manual search of scaffolds (containing Cg-contigs) stored in the “EnsemblMetazoa C. gigas” database was conducted. To confirm the presence of this set of genes obtained after scaffolds data mining, the genes were then blasted (Blast2seqs) against every single BAC-contig; genes showing a P-value < 0.001 and a coverage > 90% were then annotated. Next, using the genomic editor Apollo , and Geneious basic 5.6.5 (http://www.geneious.com), the results were individually tested and adjusted in the final editing process.
The study of functions in genome regions covered by the BAC-contigs was carried out using data mining techniques based on the information stored in the Gene Ontology (GO) database. After getting GO terms from genes annotated in the BAC-contigs (C. gigas EnsemblMetazoa) we have used the QuickGo browser (https://www.ebi.ac.uk/QuickGO) to obtain GO terms found for each BAC-contig, and then analyzed the main Molecular Function carried out in those regions. Comparison charts from the Biological Process, Molecular Function and Cellular Component domains were also prepared.
Finally, the structure of the BAC-contigs, obtained after annotation, and the information provided by the Genomicus database v25.01 using C. gigas as reference species , were compared. The object of this analysis is to test the concordance between, on the one hand, the C. gigas genome assembly and structure, available in C. gigas databases (EnsemblMetazoa C. gigas and Genomicus), and on the other, the BAC-contigs sequences. To carry out this analysis, a search was made for every gene annotated finally in BAC-contigs in the Genomicus database (using C. gigas genome as reference), and the scaffolds containing the annotated genes were scored. Then, for each BAC-contig, a map of scaffolds containing the genes was created. This information was used to compare the BAC-contig genetic structure obtained after annotation vs the C. gigas scaffolds structure from C. gigas genomic databases.
Sixty BAC sequences of C. gigas were downloaded from NCBI and assembled in BAC-contigs by means of Seq-Man software (Lasergene, DNASTAR). Eleven BACs from these contigs and one BAC more, not belonging to contigs (BAC-unitig), were used to carry out this study.
No differences in BAC-FISH results between C. angulata and C. gigas were found during the study. Results show three BACs hybridized in chromosome 2: 177C21, 50F9 and 102 L1 (Fig. 1). The last two are located in one arm of the chromosome. BAC 177C21 is located in the other arm, in subcentromeric position. Chromosome 3 contains the H4 histone gene and, in the same arm, the clones Ba17L12 and Ba18I16 are found co-localized. BAC Ba106O6 is located in chromosome 4 next to the 5S rDNA genes. The 5S rDNA genes are located in two chromosome pairs, 4 and 5, in C. angulata and C. gigas (Cross et al. 2005; Benabdelmouna et al. 2008); this BAC is located in a pericentromeric position, in a more internal position than 5S rDNA. The BAC 3D14 is found in chromosome 5, next to another 5S rDNA locus. Following these results, these two BACs allow chromosome 4 to be distinguished from chromosome 5 in these Crassostrea oysters studied (Table 2; Figs. 1 and 2). BAC 127F5 is located in a pericentromeric position of chromosome pair 6. BACs 133H22 and 96P24 are located in chromosome pair 7; and the BACs 13P13 and 45B21 are in chromosome 8. Of the 12 BACs hybridized, none is located in chromosomes 1 or 9. The 18S rDNA gene family is located in chromosome 10 .
H3, H4 and ribosomal genes (5S and 18S) probes were also hybridized by standard FISH protocol. The results showed that H3 and H4 co-localized in two pairs of chromosomes of medium size. One of them should be chromosome pair 3, as described in the preceding paragraph, by means of FISH-BAC protocols. The additional signal of hybridization was observed in a second chromosome pair of medium-size. Analysis of the relative position of histone genes and ribosomal genes showed that they are not co-localized in either the Japanese or Portuguese oysters (Figs. 1 and 2).
Annotation and integrated genetic map
Following the bioinformatic analysis carried out to annotate BAC-contigs, 88 genes were annotated in the BAC-contigs and mapped in Crassostrea chromosomes. In Table 2, Figure 2 and Additional file 2: Table S2, all the information about gene annotation is shown. Figure 3 shows the integrated map obtained after BAC-FISH and annotation analysis.
Results show that contig-1 arises from the assembly of eight BACs and it is 271,505 bp long (Table 1; Additional file 1: Table S1). This contig-1 is located in chromosome 4 and presents 14 genes, including 3 metallothionein genes (MT) and 5 novel protein codings, which are predicted genes described in the C. gigas Ensemble Metazoa database as being of unknown function.
Contig-2 is located in chromosome 5; it is 254,895 bp long (Table 1; Additional file 1: Table S1) and it contains 11 genes, four of them with unknown function. Some of the genes are associated with the myeloid differentiation primary response protein. Genes for calcineurin-binding protein cabin-1 and acetylchitobiose are also present.
Contig-3 results from the assembly of two BACs: Ba133H22 and Ba183L3. The former was used as probe in the FISH analysis. Since this contig presented a wide region of non-homologies between the two BACs, the annotation was made mainly based on the BAC hybridized. After annotations, it can be seen that contig-3 includes 7 genes including a thyroid transcription factor and an orexin receptor. Three genes with unknown function are also found.
Contig-4 emerges from the assembly of 5 BAC sequences and has a length of 222,359 bp (Table 1; Additional file 1: Table S1). From the five BACs included in this contig, the BAC named Ba13P13 was used as a probe. Contig-4 includes the largest number of genes of all contigs analyzed. It contains 19 genes; four of them are novel protein genes. The 15 genes with known function include metalloproteinase inhibitor 1, tissue inhibitor of metalloproteinase (TIMP) and acetyltransferase (ats1), among others. Contig 4 is located in chromosome 8.
Contig-5 arises from the assembly of four BACs and it is 204,747 bp long (Table 1; Additional file 1: Table S1). To locate the contig in chromosomes of the two oyster species, BAC Ba17L12 was used as a probe. It localizes in chromosome 3 (Figs. 1 and 2). This contig has annotated 10 genes including ubiquitin-protein ligase MYLIP, pyruvate dehydrogenase phosphatase regulatory subunit, and sodium and chloride-dependent betaine transporter, among others. Two genes of unknown function were also annotated (Table 2; Additional file 1: Tables S1 and Additional file 2 Table S2).
Contig-6 includes 3 BAC sequences and it is 141,754 bp long (Table 1; Additional file 1 Table S1). BAC Ba96P24 was used as a probe to localize the contig and it hybridized in chromosome 7. The annotations of this contig showed that it contains 7 genes, five of them with known functions. Several of the genes, such as thyroid transcription factor, orexin receptor type 2 and cys-loop ligand-gated ionic channel genes, are also present in contig 3, but others like galanin receptor type 2 and alpha-2A adrenergic receptor, are not. Contig-3 is also present in chromosome 7 and co-localizes with this contig-6.
Contig-7 is 237,501 bp long (Table 1; Additional file 1: Table S1) and it emerges from 3 BACs sequences. It includes 12 genes, eight of them with functions including heavy metal-binding protein, thioredoxin-like protein 4A and actin-binding protein 1. Like contig-4, it localizes in chromosome 8.
Contig-8 comes from the assembly of two BAC sequences, Ba18I16 and Ba8F15, and it is 238,060 bp long (Table 1; Additional file 1: Table S1). The Ba18I16 BAC was used as a probe to localize it in chromosomes of the oysters C. angulata and C. gigas. Contig-8 co-localizes with contig-5 in chromosome 3. Genes annotated are the same in both contigs, except for ribonuclease P protein subunit p40, which is absent in contig-8.
Chromosome 2 showed a hybridization signal with 3 BACs, corresponding to three contigs: contig-9, 10 and 11. Contig-9 results from the assembly of 3 BAC sequences, and is 198,436 bp long (Table 1; Additional file 1: Table S1). The annotation of contig-9 showed 6 genes with known function and 1 novel protein coding. Among the six genes can be found DNA replication complex GINS protein PSF1, adenylosuccinate synthetase, and microtubule-associated proteins 1A/1B light chain 3C. Contig-10 is 195,233 bp long and arises from two BACs (Table 1; Additional file 1: Table S1). It has 5 genes, but only two of them (ADMAMTS-like protein and F-box protein 22) with known function. Contig-11 did not show a match with any genes annotated in the NCBI C. gigas database, so it was BLAST-searched against the nucleotide collection (nr/nt) and filtered for mollusks in the NCBI database. Two genes were found: metallothionein (C. virginica IIG) and filamin-A-like (C. gigas predicted).
The unitig BAC Ba127F5 is 132,710 bp long and it has 6 genes and one novel protein coding (Tables 1 and 2; Additional file 1: Table S1). Included in this BAC are the galanin receptor type 2, orexin receptor type 2, and cys-loop ligand-gated ionic channel genes.
The GO analysis showed that most genes in BAC-contigs showed Molecular Functions (GO term category) related to protein binding. In particular, we have observed an enrichment in the GO terms “metal ion binding” and “zinc ion binding” in most BAC-contigs: 1, 2, 5, 7–11 and unitig BAC 127F5 (Table 3; Additional file 3: Table S3; Additional file 4: Figure S1). Other functions include “hydrolase activity” (BAC-contigs 2 and 4), “GTP binding” (2, 7 and 9) or “G-protein coupled receptor activity” (BAC-contigs 3, 6 and BAC-unitig) were also enriched (Table 3).
Genomic structure of BAC-contigs
The structure of the BAC-contigs, obtained after gene annotations, and the information about C. gigas gene order obtained from the Genomicus database, were compared. After testing concordance between the BAC-contigs structure and the scaffold assembly available in the C. gigas genome database (Table 4; Additional file 5: Figure S2), it was observed that most of the BAC-contigs structure were not concordant with the gene structure of the C. gigas scaffolds displayed in the Genomicus database. BAC-contig 1 presented 15 genes belonging partially to 4 scaffolds displayed in the Genomicus database, and they were incompatible with each other (Table 3; Additional file 4: Figure S1). BAC-contig 2 contains 11 genes and their relative position and genetic structure were not comparable to what is observed in the scaffolds of the Genomicus databases, where genes were found in partial regions of scaffold 870 and scaffold 617. BAC-contigs 3, 4, 5, 7, 8 and BAC 127F5 also presented scaffold incompatibility, being contained in a range of 2–4 incompatible scaffolds. Only in scaffolds 6, 9 and 10 could be observed a genetic structure compatible with scaffolds from the C. gigas genome databases, because genes annotated were present in single or several compatible scaffolds.
FISH has been proved to be an efficient method for correlating genetic and cytogenetic maps using marker-anchored BAC clones .
In the cytogenetic map of C. angulata and C. gigas oysters, following our results, we were able to identify 8 out of 10 chromosome pairs using BAC clones as probes. These results are relevant because the karyotype of the Crassostrea species presents a small range of morphology and sizes [10, 28], and it has not been possible to identify every one. In the oysters C. angulata and C. gigas, 5S rDNA genes are located in chromosomes 4 and 5 [28, 45] and the 5.8–18-28S rDNA genes in chromosome 10 [15, 46]. Our results show that chromosomes 4 and 5 can be identified by the position of Ba106O6 and Ba3D14 BAC clones. In relation to the H4 histone gene, when using BAC-FISH conditions, this gene was observed in chromosome pair 3. However, after using standard FISH protocols that are usually applied to multigene families , a second signal was detected in a medium-size chromosome. In C. gigas the location of the H3 gene has been described , but the position described is not consistent with our findings, since those authors described two hybridization signals on chromosome pairs 4 and 10. It could be argued that the position 3 observed in our results is different from chromosome pair 4, owing to an absence of chromosomal markers in that work leading to a mistake being made between the two chromosome pairs (3 and 4) that are almost identical in size and morphology. Our results show that the position is not in chromosome 4 as described by Boully et al. , because that position is occupied by 5S rDNA genes [28, 45]. However an even more relevant result is that histone genes are not in chromosome 10 , because that is the chromosome that bears major ribosomal genes described in the Pacific and Portuguese oysters [15, 46], and the double FISH carried out in this study confirms that the two families are not co-localized. For that reason, the location of the histone genes should be revised because contradictory results have been obtained in those experiments.
The GO term analysis showed an over-representation in BAC-contigs of molecular function associated with metal ion binding and zinc ion binding. In previous studies, enrichment was found in genes coding for proteins implicated in the detoxification of metal ions in polluted areas, in C. angulata, in comparison with those of C. gigas . In this work we have localized a metallothionein gene in chromosome 2 and 4, in the BAC-contig 1. Metallothioneins (MTs) are a superfamily of ubiquitously-expressed metal-binding proteins that can be upregulated by exposure to metal, oxidative stress and immune system challenge [47, 48]. MTs are well-known small, cystein-rich proteins that have the ability to bind with high affinity to heavy metal ions, and whose synthesis is regulated by concentrations of metals . MTs play roles in the metabolism of non-toxic essential metals (zinc and copper), as well as toxic heavy metals (cadmium, mercury, etc) . Species of the genus Crassostrea have been used in ecotoxicological bioassays to assess the bioaccumulation of Cd and Zn in field conditions  and in transplant assays [52, 53]. When comparing the genome of C. angulata with that of C. gigas, the GO terms (which describe molecular function) associated with binding to ions of zinc and other metals present in C. angulata account for the largest number of polymorphic genes (Single nucleotide polymorphisms: SNP) observed . In that context the location of MT genes in Crassostrea gigas/angulata genome represents an important advance in the study of those genes in relation to chromosome mapping.
In this work we have assembled BAC sequences available in the GeneBank database of the C. gigas oyster and obtained BAC-contig sequences longer than BAC sequences. This approach has allowed us to annotate larger regions than can be done at the level of BACs, and provides knowledge of larger regions of the genome. After annotation we have elucidated the structure of BAC-contigs and, after comparison with the scaffolds of the C. gigas genome database, a large number of inconsistencies have been revealed. In many cases, several genes belonging to the same BAC-contigs were mapped in partial regions of different scaffolds in the C. gigas database. The main consequence is that a substantial fraction of the genome scaffolds reported by Zhang et al.  are incorrectly assembled. These results are in line with those previously published by Hedgecock et al. , who after making second-generation linkage maps, with SNP markers, for the Pacific oyster C. gigas, observed that more than 38% of the scaffolds with two or more SNPs contained SNPs that unexpectedly map to two or more linkage groups. In our study it is observed that the more genes per BAC-contig, the more scaffolds to which contigs map. Thus, the higher the level of assembly, the greater the number of mistakes (i.e. mis-assemblies) that are made. In C. gigas, it has been reported previously that classifying scaffolds by the number of SNPs they contain and by the number of linkage groups to which they map reveals a highly significant, positive association between the two factors. Our data reinforce previous results that claim that misassembled scaffolds comprise large blocks, with the implication that, if the ordering of contigs and genes within scaffolds cannot always be accurate, results about the clustering of gene families, like Hox genes , should be regarded as provisional until the assembly of genome scaffolds is either confirmed or improved . Currently, the C. gigas genome (oyster assembly v9) is assembled to scaffold level and more than 7000 scaffolds are presented . The use of this technique in other species of shellfish and oysters will open a huge potential for direct and fast analysis of genomes in species of high commercial value.
Molecular markers for 8 out 10 chromosomes have been developed in Crassostrea angulata/gigas that allow the differentiation of all chromosomes in the species. An integrated genetic map anchored to the genome has been obtained. Important genes related to resistance to abiotic stress have been mapped.
Food and Agriculture Organization of the United Nations, FAO. 2016. http://www.fao.org/fishery/statistics/en.
Wang Y, Guo X. ITS length polymorphism in oysters and ITS use in species identification. J Shellfish Res. 2008;27(3):489–93.
Wang H, Quian L, Liu X, Zhang G, Guo X. Classification of a common cupped oyster from southern China. J Shellfish Res. 2010;29(4):857e66.
Langdon C, Evans F, Jacobson D, Blouin M. Yields of cultured Pacific oysters Crassostrea gigas Thunberg improved after one generation of selection. Aquaculture. 2003;220:227–44.
Hedgecock D, Davis JP. Heterosis for yield and crossbreeding of the Pacific oyster Crassostrea gigas. Aquaculture. 2007;272(Supplementary 1):S17–29.
Dégremont L, Bédier E, Boudry P. Summer mortality of hatchery-produced Pacific oyster spat (Crassostrea gigas). II. Response to selection for survival and its influence on growth and yield. Aquaculture. 2010;299:21–9.
Hedgecock D, Gaffney PM, Goulletquer P, Guo X, Reece K, Warr GW. The case for sequencing the Pacific oyster genome. J Shellfish Res. 2005;24:429–41.
Curole JP, Hedgecock D. Bivalve genomics: complications, challenges, and future perspectives. In: Liu ZJ, editor. Aquaculture Genome Technologies. Ames, Iowa: Blackwell publishing; 2007. p. 525–43.
Wang J, Li L, Zhang G. A high-density SNP genetic linkage map and QTL analysis of growth-related traits in a hybrid family of oysters (Crassostrea gigas x Crassostrea angulata) using genotyping-by-sequencing. G3 (Bethesda). 2016;6(5):1417–26.
Thiriot-Quiévreux C. Review of the literature on bivalve cytogenetics in the last ten years. Cah Biol. 2002;43:17–26.
Michinina S, Rebordinos L. Genetic differentiation in marine and estuarine natural populations of Crassostrea angulata. Mar Ecol Prog Ser. 1997;154:167–74.
Insua A, Thiriot-Quiévreux C. Karyotype and nucleolus organizer regions in Ostrea puelchana (Bivalvia: Ostreidae). Veliger. 1993;36:215–9.
Li XX, Havenhand JN. Karyotype, nucleolus organiser regions and constitutive heterochromatin in Ostrea angasi (Molluscae: Bivalvia): evidence of taxonomic relationships within the Ostreidae. Mar Biol. 1997;127:443–8.
Thiriot-Quiévreux C, Insua A. Nucleolar organiser region variation in the chromosomes of three oyster species. J Exp Mar Biol Ecol. 1992;157:33–40.
Cross I, Vega L, Rebordinos L. Nucleolar organizing regions in Crassostrea angulata: Chomosomal location and polymorphism. Genetica. 2003;119:65–74.
Leitão A, Chaves R, Santos S, Guedes-Pinto H, Boudry P. Restriction enzyme digestion chromosome banding in Crassostrea gigas, Crassostrea angulata. Ostrea edulis and Ostrea conchaphila Comparative karyological analysis within Ostreidae Genome. 2004;47:781–8.
Leitão A, Chaves R, Santos R, Guedes-Pinto H, Boudry P. Interspecific hybridisation in oysters: restriction enzyme digestion chromosome banding confirms Crassostrea angulata × Crassostrea gigas F1 hybrids. J Exp Mar Biol Ecol. 2007;343:253–60.
Feng L, Hu L, Fu X, Liao H, Li X, Zhan A, et al. An integrated genetic and cytogenetic map for Zhikong scallop, Chlamys farreri, based on microsatellite markers. PLoS One. 2014;9(4):e92567.
Yoshido A, Bando H, Yasukochi Y, Sahara K. The Bombyx mori karyotype and the assignment of linkage groups. Genetics. 2005;170:675–85.
Huang Y, Zhao Y, Haley CS, Hu S, Hao J, Wu C, et al. A genetic and cytogenetic map for the duck (Anas platyrhynchos). Genetics. 2006;173(1):287–96.
Stratikopoulos EE, Augustinos AA, Petalas YG, Vrahatis MN, Mintzas A, Mathiopoulos KD, et al. An integrated genetic and cytogenetic map for the Mediterranean fruit fly, Ceratitis capitata, based on microsatellite and morphological markers. Genetica. 2008;133:147–57.
Ren Y, Zhang ZH, Liu JH, Staub JE, Han YH, Cheng Z, et al. An integrated genetic and cytogenetic map of the cucumber genome. PLoS One. 2009;4(6):e5795.
Portela-Bens S, Merlo MA, Rodríguez ME, Cross I, Manchado M, Kosyakova N, et al. Integrated gene mapping and synteny studies give insights into the evolution of a sex proto-chromosome in Solea senegalensis. Chromosoma. 2017;126:261–77.
Phillips RB, Amores A, Morasch MR, Wilson C, Postlethwait JH. Assignment of zebrafish genetic linkage groups to chromosomes. Cytogenet Genome Res. 2006;114:155–62.
Úbeda-Manzanaro M, Merlo MA, Palazón JL, Cross I, Sarasquete C, Rebordinos L. Chromosomal mapping of the major and minor ribosomal genes, (GATA)n and U2 snRNA gene by double-colour FISH in species of the Batrachoididae family. Genetica. 2010;138:787–94.
García-Cegarra A, Merlo MA, Ponce M, Portela-Bens S, Cross I, Manchado M, et al. Preliminary genetic map in Solea senegalensis (Pleuronectiformes, Soleidae) using BAC-FISH and next-generation sequencing. Cytogenet Genome Res. 2013;141:227–40.
Taboada X, Pansonato-Alves JC, Foresti F, Martínez P, Viñas A, Pardo BG, et al. Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones. Chromosoma. 2014;123:281–91.
Cross I, Díaz E, Sánchez I, Rebordinos L. Molecular and cytogenetic characterization of Crassostrea angulata chromosomes. Aquaculture. 2005;247:135–44.
López-Flores I, Ruiz-Rejón C, Cross I, Rebordinos L, Robles F, Navajas-Pérez R, et al. Molecular characterization and evolution of an interspersed repetitive DNA family of oysters. Genetica. 2010;138:1211–9.
Bouilly K, Chaves R, Leitão A, Benabdelmouna A, Guedes-Pinto H. Chromosomal organization of simple sequence repeats in the Pacific oyster (Crassostrea gigas): (GGAT)(4), (GT)(7) and (TA)(10) chromosome patterns. J Genet. 2008;87(2):119–25.
Wang Y, Xu Z, Guo X. Differences in the rDNA-bearing chromosome divide the Asian-Pacific and Atlantic species of Crassostrea (Bivalvia, Mollusca). Biol Bull. 2004;206:46–54.
Doenecke D, Albig W, Bode C, Drabent B, Franke K, Gavenis K, et al. Histones: genetic diversity and tissue-specific gene expression. Histochem Cell Biol. 1997;107(1):1–10.
Armbruster GFJ, Böhme M, Bernhard D, Schlegel M. The H3/H4 histone gene cluster of land snails (Gastropoda: Stylommatophora): TS/TV ratio, GC3 drive and signals in stylommatophoran phylogeny. J Mollusc Stud. 2005;71(4):339–48.
González –RR, Ausió J, Méndez J, Eirín-López JM. Early evolution of histone genes: prevalence of an “Orphon” H1 lineage in protostomes and birth-and-death process in the H2A family. J Mol Evol. 2008;66(5):505–18.
Oskars TR, Bouchet P, MAE M. A new phylogeny of the Cephalaspidea (Gastropoda: Heterobranchia) based on expanded taxon sampling and gene markers. Mol Phylogenet Evol. 2005;89:130–50.
Combosch DJ, Giribet G. Clarifying phylogenetic relationships and the evolutionary history of the bivalve order Arcida (Mollusca: Bivalvia: Pteriomorphia). Mol Phylogenet Evol. 2016;94:298–312.
Cross I, Merlo MA, Rodríguez ME, Portela-Bens S, Rebordinos L. Adaptation to abiotic stress in the oyster Crassostrea angulata relies on genetic polymorphisms. Fish Shellfish Inmunol. 2014;41:618–24.
Hedgecock D, Shin G, Gracey AY, Den Berg DV, Samanta MP. Second-generation linkage maps for the Pacific oyster Crassostrea gigas reveal errors in assembly of genome scaffolds. G3 (Bethesda). 2015;5(10):2007–19.
Wang Y, Xu Z, Pierce JC, Guo X. Characterization of eastern oyster (Crassostrea virginica Gmelin) chromosomes by fluorescence in situ hybridization with bacteriophage P1 clones. Mar Biotechnol. 2005;7:207–14.
Zhang G, Fang X, Guo X, Li L, Luo R, Xu F, et al. The oyster genome reveals stress adaptation and complexity of shell formation. Nature. 2012;490(7418):49–54.
Mazzuchelli J, Kocher TD, Yang F, Martins C. Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish. BMC Genomics. 2012;13:463.
Bouilly K, Chaves R, Fernandes M, Guedes-Pinto H. Histone H3 gene in the Pacific oyster, Crassostrea gigas Thunberg, 1793: molecular and cytogenetic characterisations. Comp Cytogenet. 2010;4(2):111–21.
Lewis SE, Searle SM, Harris N, Gibson M, Lyer V, Richter J, et al. Apollo: a sequence annotation. Genome Biol. 2002;3:0082.
Louis A, Nguyen NTT, Muffato M, Crollius HR. Genomicus update 2015: KaryoView and MatrixView provide a genome-wide perspective to multispecies comparative genomics. Nucleic Acids Res. 2015;43:682–9.
Benabdelmouna A, Nadolna K, Ledu C. Fluorescence In Situ hybridization (FISH) of rDNA genes (5S and 18–5.8-28S) as a tool for chromosomal tagging and to assess the occurrence of somatic and sperm aneuploidy in Pacific oyster Crassostrea gigas. XVI Plant and Animal Genome Conference (PAG). 2008; http://archimer.ifremer.fr/doc/00000/3502/.
Xu Z, Guo X, Gaffney PM, Pierce JC. Chromosomal location of the major ribosomal RNA genes in Crassostrea virginica and Crassostrea gigas. Veliger. 2001;44(1, 2):79–83.
Jenny MJ, Ringwood AH, Schey K, Warr GW, Chapman RW. Diversity of metallothioneins in the American oyster, Crassostrea virginica, revealed by transcriptomic and proteomic approaches. Eur J Biochem. 2014;271:1702–12.
Jenny MJ, Payton SL, Baltzegar DA, Lozier JD. Phylogenetic analysis of molluscan Metallothioneins: evolutionary insight from Crassostrea virginica. J Mol Evol. 2016;83(3):110–25.
Kojima Y, Metallothionein KJHR. Trends Biochem Sci. 1978;3:90–3.
Tanguy A, Moraga D. Cloning and characterization of a gene coding for a novel metallothionein in the Pacific oyster Crassostrea gigas (CgMT2): a case of adaptive response to metal-induced stress? Gene. 2001;273:123–30.
Rebelo MF, Amaral MCR, High Zn PWC. Cd contamination in oyster from a contaminated coastal lagoon and its role in contamination assessment. Mar Pollut Bull. 2003;46:1354–8.
Amaral MCR, Rebelo MF, Torres JPM, Pfeiffer WC. Bioaccumulation and depuration of Zn and cd in mangrove oysters (Crasssotrea rhizophorae, Guilding, 1828) transplanted to and from a contaminated tropical coastal lagoon. Mar Environ Res. 2005;59:277–85.
Macías-Mayorga D, Laiz I, Moreno-Garrido I, Blasco J. Is oxidative stress related to cadmium accumulation in the Mollusc Crassostrea angulata? Aquat Toxicol. 2015;161:231–41.
Paps J, Xu F, Zhang G, Holland PW. Reinforcing the egg-timer: recruitment of novel lophotrochozoa homeobox genes to early and late development in the Pacific oyster. Genome Biol Evol. 2015;7(3):677–88.
The authors acknowledge their gratitude to the Ensemble team and Genomicus Browser curators for providing a browser and database that greatly facilitated the use and analyses of the BACs used in this work. We also thank all the institutions, consortium leaders, coordinators, and supervisors of genome projects that have made genome information available.
This study has been supported by the AQUAGENET project (SOE 2/ P1/E287) from the INTERREG IVB SUDOE program and by grants from the Junta de Andalucía (Spain) to the PAI BIO-219 group. The role of the funding bodies is limited to direct funding of the project activities in the field and laboratory that results in this manuscript.
Availability of data and materials
The data that support the findings of this study are available in the supplementary materials.
Ethics approval and consent to participate
Not applicable. The animals studied in this work are bivalve mollusks, non-cephalopod invertebrates. No field permissions were required to obtain these samples for the research.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary of BAC-contigs statistics. (DOCX 13 kb)
Table S2. Data of BAC-contigs annotations. (XLSX 19 kb)
Table S3. Summary of Gene Ontology term enrichment statistics. (XLSX 2457 kb)
Figure S1. Comparison chart of significant Gene Ontology terms of genes located in BAC-contigs. The arrows indicate the relationship among the GO categories. Colored boxes indicate enriched Gene Ontology terms found in the study. (PDF 2475 kb)
Figure S2. Comparative analysis of genomic structure of BAC-contig and Genomicus browser. (PDF 356 kb)