Characterization of Gonadal Transcriptomes from the Turbot (Scophthalmus Maximus)
Total Page:16
File Type:pdf, Size:1020Kb
Genome Characterization of Gonadal Transcriptomes from the Turbot (Scophthalmus maximus) Journal: Genome Manuscript ID gen-2014-0190.R3 Manuscript Type: Article Date Submitted by the Author: 02-Oct-2015 Complete List of Authors: Hu, Yulong; Yellow Sea Fisheries Research Institute (YSFRI), Huang, Meng; cornell university, Wang, Weiji; Yellow Sea Fisheries Research Institute (YSFRI), Guan, Jiantao;Draft Yellow Sea Fisheries Research Institute (YSFRI), Kong, Jie; Yellow Sea Fisheries Research Institute (YSFRI), Keyword: Turbot, Transcriptome, sex, SNP Note: The following files were submitted by the author for peer review, but cannot be converted to PDF. You must view these files (e.g. movies) online. gen-2014-0190.R3 Supplement1-Dataset S1.rar gen-2014-0190.R3 Supplement2-Dataset S2.zip https://mc06.manuscriptcentral.com/genome-pubs Page 1 of 701 Genome Characterization of Gonadal Transcriptomes from the Turbot ( Scophthalmus maximus ) Yulong Hu #, Meng Huang #, Weiji Wang, Jiantao Guan, Jie Kong* Yellow Sea Fisheries Research Institute (YSFRI), Chinese Academy of Fishery Sciences, Qingdao, China Draft Corresponding author: Jie Kong, PhD Yellow Sea Fisheries Research Institute (YSFRI) Chinese Academy of Fishery Sciences Qingdao, China Phone: 0086-532-85821650 Fax: 0086 053285811514 Email: [email protected] 1 https://mc06.manuscriptcentral.com/genome-pubs Genome Page 2 of 701 Abstract The mechanisms underlying sexual reproduction and sex ratio determination remains unclear in turbot, a flatfish of great commercial value. And there is limited information in the turbot database regarding genes related to the reproductive system. Here, we conducted high-throughput transcriptome profiling of turbot gonad tissues to better understand their reproductive functions and to supply essential gene sequence information for marker-assisted selection programs in the turbot industry. In this study, two gonad libraries representing sex differences in S.maximus yielded 453,818 high-quality reads that were assembled into 24,611 contigs and 33,713 singletons by using 454 pyrosequencing, 13,936Draft CS (contigs and singletons) of which were annotated using BLASTx. GO (Gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses, revealed that various biological functions and processes were associated with many of the annotated CS (contigs and singletons). Expression analyses showed that 510 genes were differentially expressed in males versus females; 80% of these genes were annotated. In addition, 6,484 and 6,036 single nucleotide polymorphisms (SNPs) were identified in males and females libraries, respectively. This transcriptome resource will serve as the foundation for cDNA or SNP microarray construction, gene expression characterization, and sex-specific linkage mapping in turbot. Key words: Turbot; Transcriptome; sex; SNP 2 https://mc06.manuscriptcentral.com/genome-pubs Page 3 of 701 Genome Introduction The turbot ( Scophthalmus maximus ), also known as Psetta maxima, is an important economic flatfish species belonging to the family Scophthalmidae (order Pleuronectiformes). It is a marine flatfish species naturally distributed all along the European coasts, from Norway to Morocco in the Atlantic area, and across the northern coasts of the Mediterranean up to the Black Sea (Blanquer et al. 1992). Following, it was introduced to China as well for farming purposes in 1992 (Lei 2003). In 2012, turbot aquaculture production was 12,676 tons in Europe, while 64,000 tons in China, which represents a world leader in adult turbot production (FAO 2014). During intensive farming, turbot females can exhibit over two-fold increase in their growth rate as compared to malesDraft (Husebye et al. 1994). This is among the highest sex-related differential growth rates in cultured marine fish (Piferrer et al. 1995). The difference in the growth rate could be observed as early as 8 months after hatching (Imsland et al. 1997), and it is maintained throughout the production cycle, even following sexual maturation (Piferrer et al. 2004). Under culture conditions, the males also mature earlier than females, which can reduce somatic growth and increase susceptibility to diseases and mortality (Imsland et al. 1997; Piferrer et al. 2003). Altogether, these factors result in loss of production during the grow-out phase and, therefore, all-female populations should be established to optimize turbot production. The sex ratio depends on processes related to both genetic and environmental factors, which are directly related to the reproductive potential of the population (Penman et al. 2008). Flatfish have been reported to possess genetic systems ranging from polygenic to dominant sex-determining factors that are influenced by autosomal controls in male (XX/XY) or female heterogametic systems (ZZ/ZW) (Purdom 1972; Lincoln 1981; Tabata 1991; Howell et al. 1995; Aida and Arai 1998; Goto et al. 1999; 3 https://mc06.manuscriptcentral.com/genome-pubs Genome Page 4 of 701 Tvedt et al. 2006 ). The most common mode of GSD is determined by the inherited combination of sex-determining genes or polygenic systems (Oldfield 2005 ). Recently, the number of identified expressed sequence tags (ESTs) has increased very rapidly in turbot, it remains insufficient when compared with some other teleost (e.g. 15,598 vs 301,380; Salmon, GenBank/NCBI). Although several sex-associated markers have been discovered (Casas et al. 2005 ) and some linkage maps constructed or enriched (Bouza et al. 2007; Bouza et al. 2008; Bouza et al. 2012; Martínez et al. 2008 ), the sex-determining gene remains not defined and precisely sex-associated linkage maps are not available. However, development of a database of genes sequences has led to the identification of the major sex-determining region (Martínez et al. 2009 ) and mapping of DNA sex-specific markers and genes related to sex differentiation in turbot (Viñas et al. 2012). Draft RNA sequencing (RNA-seq) was used for transcriptome profiling, allowing us to rapidly identify the vast majority of transcriptomes and new splicing variants in a cost-effective manner. The Genome Sequencer FLX (GS FLX) System, powered by 454 sequencing technology, was selected as the sequencing platform because this technology is more appropriate than other next-generation sequencing platforms (such as the Illumina Solexa and ABI SOLiD) for de novo transcriptome sequencing in non-model organisms (Vera et al. 2008; Meyer et al. 2009; O’Neil et al. 2010 ). In particular, many transcriptome profiles of some marine species have been successfully mapped by using 454 sequencing (Kawahara-Miki et al. 2011; Milan et al. 2011; Hou et al. 2011; Du et al. 2012; Pereiro et al. 2012 ). Recent studies have reported large-scale transcriptome profiling of S.maximus using 454 sequencing to identify candidate genes involved in immunity and sex differentiation (Pereiro et al. 2012; Rabis et al. 2013 ). Although the data (more than 4 https://mc06.manuscriptcentral.com/genome-pubs Page 5 of 701 Genome 120,876 contigs and 237,580 singletons) obtained from those tissues (brain, hypophysis, gonad, gill, liver, spleen and head kidney) had found some genes involved in sex differentiation and maturation, more sequencing information is still needed to identify novel, key genes involved in gonad-specific or sex-associated genes. To supplement the genomic resources of the turbot, especially with respect to the sex transcriptome, we used a 454 pyrosequencing method to identify differentially expressed genes and potential genetic markers in gonad tissues. The data generated in this study will facilitate further research examining sex- and reproductive-related genes in the turbot, which will be useful for linkage map construction. Materials and methods Draft Ethics Statement The Institutional Animal Care and Use Committee (IACUC) of the Yellow Sea Fisheries Research Institute (YSFRI) approved this animal study. All of the fishes were sacrificed by an overdose of anesthesia to minimize suffering, and the experimental procedures were performed according to the institutional guidelines of the State Council of the People’s Republic of China. Turbot sampling Gonad tissues from 12 male and 12 female turbot (aged 600 days) were collected for 454 sequencing. Artificial fertilization and maintenance of the larval cultures were performed by the Huanghai Aquatic Product Group Co., Ltd. (Haiyang, Shandong Province, China). Twenty juvenile fishes (aged 150 days) were provided by the Laboratory of Genetic Resources and Breeding for independent SNP validation at the Yellow Sea Fisheries Research Institute (Qingdao, Shandong Province, China). The analyzed turbot represented a mixture of unrelated families, and the number, age and 5 https://mc06.manuscriptcentral.com/genome-pubs Genome Page 6 of 701 mean value of morphometric data (standard length, width and weight) for each animal group were as follows: adult males(n=12; 600 days post-hatch; 29.5±1.6 cm, 24.0±1.0 cm, 1398.7±88.4 g); adult females ( n=12; 600 days post-hatch; 30.7±0.9 cm, 25.0±0.6 cm, 1519.8±165.7 g); juvenile males ( n=10; 150 days post-hatch; 8.8±1.1 cm, 6.7±1.0 cm, 28.2±11.0 g); juvenile females ( n=10; 150 days post-hatch; 9.3±0.6 cm, 7.3±0.4 cm, 36.2±7.8 g). Gonads were removed from all turbot individuals for 454 sequencing and sex identification, and pterygiophore tissues were removed from juvenile turbot for SNP validation. For sex identification, the tissues were stained with hematoxylinand eosin, embedded in paraffin and sliced into 6-µm sections. Phenotypic sex identification was conducted by microscopic observation. All of the samples subjected to 454 sequencing and SNP validation were flash frozen in liquid nitrogen and stored at -80 °C untilDraft RNA extraction. Total RNA extraction, cDNA library construction and 454 GS-FLX sequencing Total RNA was extracted and purified from each tissue using the RNA prep pure Tissue Kit (Tiangen Biotech (Beijing) Co., Ltd., China) according to the manufacturer’s protocol. The quantity and quality of the total RNA was assessed using an Ultrospec TM 2100 pro UV/Visible Spectrophotometer (Amersham Biosciences, Uppsala, Sweden) and by gel electrophoresis.