Directed sequencing and annotation of three Dicentrarchus labrax L. chromosomes by applying Sanger- and pyrosequencing technologies on pooled DNA of comparatively mapped BAC clones

Molecular and cellular characterization of European sea bass CD3ε+ T lymphocytes and their modulation by microalgal feed supplementation

The CD3 coreceptor is a master T cell surface marker, and genes encoding CD3ζ, γδ, and ε chains have been reported in several teleost fish. Here, a complete cDNA sequence of CD3ɛ chain was identified from a sea bass (Dicentrarchus labrax L.) gill transcriptome. Its basal expression was quantified in both lymphoid and non-lymphoid organs of sea bass juveniles with real-time qPCR analysis. After either in vitro stimulation of head kidney leukocytes with the T-cell mitogen phytohaemagglutinin or in vivo stimulation with an orally administered Vibrio anguillarum vaccine, CD3ε expression levels increased in head kidney leukocytes, confirming that CD3ε T cells may play important roles in fish systemic protection against pathogens. Further, three peptides were designed on the CD3ɛ cytoplasmic tail region and employed as immunogens for antibody production in rabbit. One antiserum so obtained, named RACD3/1, immunostained a band of the expected size in a western blot of a sea bass thymocyte lysate. The distribution of CD3ε⁺ lymphocyte population in the lymphoid organs and mucosal tissues was addressed in healthy fish by IHC. In decreasing percentage order, CD3ε⁺ lymphocytes were detected by flow cytometry in thymus, peripheral blood leukocytes, gills, head kidney, gut, and spleen. Finally, a significant in vivo enhancement of CD3ε⁺ T intestinal lymphocytes was found in fish fed on diets in which 100% fish meal was replaced by the microalgae Nannochloropsis sp. biomass. These results indicate that CD3ε⁺ T cells are involved in nutritional immune responses.

Identification, molecular characterization and functional analysis of interleukin (IL)-2 and IL-2like (IL-2L) cytokines in sea bass (Dicentrarchus labrax L.)

In mammals, interleukin (IL)-2, initially known as a T-cell grow factor, is an immunomodulatory cytokine involved in the proliferation of T cells upon antigen activation. In bony fish, some IL-2 orthologs have been identified, but, recently, an additional IL-2like (IL-2L) gene has been found. In this paper, we report the presence of these two divergent IL-2 isoforms in sea bass (Dicentrarchus labrax L.). Genomic analyses revealed that they originated from a gene duplication event, as happened in most percomorphs. These two IL-2 paralogs show differences in the amino acid sequence and in the exon 4 size, and these features could be an indication that they bind preferentially to different specific IL-2 receptors. Sea bass IL-2 paralogs are highly expressed in gut and spleen, which are tissues and organs involved in fish T cell immune functions, and the two cytokines could be up-regulated by both PHA stimulation and vaccination with a bacterial vaccine, with IL-2L being more inducible. To investigate the functional activities of sea bass IL-2 and IL-2L we produced the corresponding recombinant molecules in E. coli and used them to in vitro stimulate HK and spleen leukocytes. IL-2L is able to up-regulate the expression of markers related to different T cell subsets (Th1, Th2 and Th17) and to Treg cells in HK, whereas it has little effect in spleen. IL-2 is not active on these markers in HK, but shows an effect on Th1 markers in spleen. Finally, the stimulation with recombinant IL-2 and IL-2L is also able to induce in vitro proliferation of HK- and spleen-derived leukocytes. In conclusion, we have demonstrated that sea bass possess two IL-2 paralogs that likely have an important role in regulating T cell development in this species and that show distinct bioactivities.

BAC-pool sequencing and analysis confirms growth-associated QTLs in the Asian seabass genome

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation

The European sea bass (Dicentrarchus labrax) is a temperate-zone euryhaline teleost of prime importance for aquaculture and fisheries. This species is subdivided into two naturally hybridizing lineages, one inhabiting the north-eastern Atlantic Ocean and the other the Mediterranean and Black seas. Here we provide a high-quality chromosome-scale assembly of its genome that shows a high degree of synteny with the more highly derived teleosts. We find expansions of gene families specifically associated with ion and water regulation, highlighting adaptation to variation in salinity. We further generate a genome-wide variation map through RAD-sequencing of Atlantic and Mediterranean populations. We show that variation in local recombination rates strongly influences the genomic landscape of diversity within and differentiation between lineages. Comparing predictions of alternative demographic models to the joint allele-frequency spectrum indicates that genomic islands of differentiation between sea bass lineages were generated by varying rates of introgression across the genome following a period of geographical isolation.

The European sea bass is an economically important fish species, which is subject to intense selective breeding. Here, the authors sequence the genome of the European sea bass and highlight gene family expansions underlying adaptation to salinity change, as well as the genomic architecture of speciation between two divergent sea bass lineages.

Advances in European sea bass genomics and future perspectives

Only recently available sequenced and annotated teleost fish genomes were restricted to a few model species, none of which were for aquaculture. The application of marker assisted selection for improved production traits had been largely restricted to the salmon industry and genetic and Quantitative Trait Loci (QTL) maps were available for only a few species. With the advent of next generation sequencing the landscape is rapidly changing and today the genomes of several aquaculture species have been sequenced. The European sea bass, Dicentrarchus labrax, is a good example of a commercially important aquaculture species in Europe for which in the last decade a wealth of genomic resources, including a chromosomal scale genome assembly, physical and linkage maps as well as relevant QTL have been generated. The current challenge is to stimulate the uptake of the resources by the industry so that the full potential of this scientific endeavor can be exploited and produce benefits for producers and the public alike.

Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones

Bacterial artificial chromosomes (BAC) have been widely used for fluorescence in situ hybridization (FISH) mapping of chromosome landmarks in different organisms, including a few in teleosts. In this study, we used BAC-FISH to consolidate the previous genetic and cytogenetic maps of the turbot (Scophthalmus maximus), a commercially important pleuronectiform. The maps consisted of 24 linkage groups (LGs) but only 22 chromosomes. All turbot LGs were assigned to specific chromosomes using BAC probes obtained from a turbot 5× genomic BAC library. It consisted of 46,080 clones with inserts of at least 100 kb and <5 % empty vectors. These BAC probes contained gene-derived or anonymous markers, most of them linked to quantitative trait loci (QTL) related to productive traits. BAC clones were mapped by FISH to unique marker-specific chromosomal positions, which showed a notable concordance with previous genetic mapping data. The two metacentric pairs were cytogenetically assigned to LG2 and LG16, and the nucleolar organizer region (NOR)-bearing pair was assigned to LG15. Double-color FISH assays enabled the consolidation of the turbot genetic map into 22 linkage groups by merging LG8 with LG18 and LG21 with LG24. In this work, a first-generation probe panel of BAC clones anchored to the turbot linkage and cytogenetical map was developed. It is a useful tool for chromosome traceability in turbot, but also relevant in the context of pleuronectiform karyotypes, which often show small hardly identifiable chromosomes. This panel will also be valuable for further integrative genomics of turbot within Pleuronectiformes and teleosts, especially for fine QTL mapping for aquaculture traits, comparative genomics, and whole-genome assembly.

Generation of physical map contig-specific sequences useful for whole genome sequence scaffolding

Along with the rapid advances of the nextgen sequencing technologies, more and more species are added to the list of organisms whose whole genomes are sequenced. However, the assembled draft genome of many organisms consists of numerous small contigs, due to the short length of the reads generated by nextgen sequencing platforms. In order to improve the assembly and bring the genome contigs together, more genome resources are needed. In this study, we developed a strategy to generate a valuable genome resource, physical map contig-specific sequences, which are randomly distributed genome sequences in each physical contig. Two-dimensional tagging method was used to create specific tags for 1,824 physical contigs, in which the cost was dramatically reduced. A total of 94,111,841 100-bp reads and 315,277 assembled contigs are identified containing physical map contig-specific tags. The physical map contig-specific sequences along with the currently available BAC end sequences were then used to anchor the catfish draft genome contigs. A total of 156,457 genome contigs (~79% of whole genome sequencing assembly) were anchored and grouped into 1,824 pools, in which 16,680 unique genes were annotated. The physical map contig-specific sequences are valuable resources to link physical map, genetic linkage map and draft whole genome sequences, consequently have the capability to improve the whole genome sequences assembly and scaffolding, and improve the genome-wide comparative analysis as well. The strategy developed in this study could also be adopted in other species whose whole genome assembly is still facing a challenge.

Applications for next-generation sequencing in fish ecotoxicogenomics

The new technologies for next-generation sequencing (NGS) and global gene expression analyses that are widely used in molecular medicine are increasingly applied to the field of fish biology. This has facilitated new directions to address research areas that could not be previously considered due to the lack of molecular information for ecologically relevant species. Over the past decade, the cost of NGS has decreased significantly, making it possible to use non-model fish species to investigate emerging environmental issues. NGS technologies have permitted researchers to obtain large amounts of raw data in short periods of time. There have also been significant improvements in bioinformatics to assemble the sequences and annotate the genes, thus facilitating the management of these large datasets.The combination of DNA sequencing and bioinformatics has improved our abilities to design custom microarrays and study the genome and transcriptome of a wide variety of organisms. Despite the promising results obtained using these techniques in fish studies, NGS technologies are currently underused in ecotoxicogenomics and few studies have employed these methods. These issues should be addressed in order to exploit the full potential of NGS in ecotoxicological studies and expand our understanding of the biology of non-model organisms.

DNA methylation of the gonadal aromatase (cyp19a) promoter is involved in temperature-dependent sex ratio shifts in the European sea bass

Sex ratio shifts in response to temperature are common in fish and reptiles. However, the mechanism linking temperature during early development and sex ratios has remained elusive. We show in the European sea bass (sb), a fish in which temperature effects on sex ratios are maximal before the gonads form, that juvenile males have double the DNA methylation levels of females in the promoter of gonadal aromatase (cyp19a), the enzyme that converts androgens into estrogens. Exposure to high temperature increased the cyp19a promoter methylation levels of females, indicating that induced-masculinization involves DNA methylation-mediated control of aromatase gene expression, with an observed inverse relationship between methylation levels and expression. Although different CpGs within the sb cyp19a promoter exhibited different sensitivity to temperature, we show that the increased methylation of the sb cyp19a promoter, which occurs in the gonads but not in the brain, is not a generalized effect of temperature. Importantly, these effects were also observed in sexually undifferentiated fish and were not altered by estrogen treatment. Thus, methylation of the sb cyp19a promoter is the cause of the lower expression of cyp19a in temperature-masculinized fish. In vitro, induced methylation of the sb cyp19a promoter suppressed the ability of SF-1 and Foxl2 to stimulate transcription. Finally, a CpG differentially methylated by temperature and adjacent to a Sox transcription factor binding site is conserved across species. Thus, DNA methylation of the aromatase promoter may be an essential component of the long-sought-after mechanism connecting environmental temperature and sex ratios in vertebrate species with temperature-dependent sex determination.

A pilot study for channel catfish whole genome sequencing and de novo assembly

Background

Recent advances in next-generation sequencing technologies have drastically increased throughput and significantly reduced sequencing costs. However, the average read lengths in next-generation sequencing technologies are short as compared with that of traditional Sanger sequencing. The short sequence reads pose great challenges for de novo sequence assembly. As a pilot project for whole genome sequencing of the catfish genome, here we attempt to determine the proper sequence coverage, the proper software for assembly, and various parameters used for the assembly of a BAC physical map contig spanning approximately a million of base pairs.

Results

A combination of low sequence coverage of 454 and Illumina sequencing appeared to provide effective assembly as reflected by a high N50 value. Using 454 sequencing alone, a sequencing depth of 18 X was sufficient to obtain the good quality assembly, whereas a 70 X Illumina appeared to be sufficient for a good quality assembly. Additional sequencing coverage after 18 X of 454 or after 70 X of Illumina sequencing does not provide significant improvement of the assembly. Considering the cost of sequencing, a 2 X 454 sequencing, when coupled to 70 X Illumina sequencing, provided an assembly of reasonably good quality. With several software tested, Newbler with a seed length of 16 and ABySS with a K-value of 60 appear to be appropriate for the assembly of 454 reads alone and Illumina paired-end reads alone, respectively. Using both 454 and Illumina paired-end reads, a hybrid assembly strategy using Newbler for initial 454 sequence assembly, Velvet for initial Illumina sequence assembly, followed by a second step assembly using MIRA provided the best assembly of the physical map contig, resulting in 193 contigs with a N50 value of 13,123 bp.

Conclusions

A hybrid sequencing strategy using low sequencing depth of 454 and high sequencing depth of Illumina provided the good quality assembly with high N50 value and relatively low cost. A combination of Newbler, Velvet, and MIRA can be used to assemble the 454 sequence reads and the Illumina reads effectively. The assembled sequence can serve as a resource for comparative genome analysis. Additional long reads using the third generation sequencing platforms are needed to sequence through repetitive genome regions that should further enhance the sequence assembly.