High-quality genome assembly of channel catfish, Ictalurus punctatus
GigaScience volume 5, Article number: 39 (2016)
The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world, especially in the United States’ aquaculture industry. Since its introduction into China in 1984, both cultivation area and yield of this species have been dramatically increased such that China is now the leading producer of channel catfish. To aid genomic research in this species, data sets such as genetic linkage groups, long-insert libraries, physical maps, bacterial artificial clones (BAC) end sequences (BES), transcriptome assemblies, and reference genome sequences have been generated. Here, using diverse assembly methods, we provide a comparable high-quality genome assembly for a channel catfish from a breeding stock inbred in China for more than three generations, which was originally imported to China from North America.
Approximately 201.6 gigabases (Gb) of genome reads were sequenced by the Illumina HiSeq 2000 platform. Subsequently, we generated high quality, cost-effective and easily assembled sequences of the channel catfish genome with a scaffold N50 of 7.2 Mb and 95.6 % completeness. We also predicted that the channel catfish genome contains 21,556 protein-coding genes and 275.3 Mb (megabase pairs) of repetitive sequences.
We report a high-quality genome assembly of the channel catfish, which is comparable to a recent report of the “Coco” channel catfish. These generated genome data could be used as an initial platform for molecular breeding to obtain novel catfish varieties using genomic approaches.
Library construction, read sequencing and filtering
To generate genome sequence data, genomic DNA from mixed tissues (including muscle and skin) of channel catfish was extracted from a chosen individual cultured at a local base of the Freshwater Fisheries Research Institute (Jiangsu Province, Nanjing, China) using Qiagen GenomicTip100 (Qiagen, Hilden, DE) as per standard protocols. Isolated genomic DNA was subsequently used to construct short-insert libraries (250, 500 and 800 bp) and long-insert libraries (2, 5, 10 and 20 kb) with the standard protocol provided by Illumina (San Diego, USA). Paired-end sequencing was performed using the Illumina HiSeq 2000 platform to generate 125-bp reads using a whole genome shotgun sequencing (WGS) strategy .
To improve the quality of sequenced reads, we trimmed 4 bases with edges from the reads of short-insert libraries and long-insert libraries, discarded duplicated reads from the long-insert libraries, and removed reads containing 10 or more Ns and low-quality bases. Finally, a total of 201.6-Gb clean reads were generated for further genome assembly.
Genome assembly and quality assessments
At first, we estimated the channel catfish genome size using k-mer analysis  with the formula: G = N*(L − 17 + 1)/K_depth, where N is the total number of reads, and K_depth indicates the frequency of reads occurring more frequently than the others. The calculated genome size is 0.839 Gb, which is shorter than that (1 Gb) from a 2016 report of an American-native channel catfish .
Simultaneously, we employed SOAPdenovo2 (version 2.04.4) software  with optimized parameters (pregraph −K 27 −d 1; contig −M 1; scaff −F −b 1.5 −p 16) to link sequenced reads to contigs and original scaffolds. All reads were then aligned onto the contigs for scaffold construction by utilizing long-insert paired-end information, which was subsequently supplied to link contigs to scaffolds in a step-wise manner. Gaps were closed using approximately 480 million of Illumina paired-end reads generated from the three libraries with insert sizes of 250, 500 and 800 bp as the input for GapCloser (v1.12-r6, default parameters and −p set to 25) . A final genome assembly of 0.845 Gb in length was obtained (Table 1), which is slightly shorter than that (0.942 Gb) of a recently reported a American-native channel catfish genome . The calculated contig N50 was 48.5 kilobases (kb), and the scaffold N50 was 7.2 Mb (Table 1). These values are also comparable to those in  (see details in Table 2).
Two typical methods were then used to assess the quality and completeness of the generated assembly. First, transcriptome evaluation was used to assess the completeness of gene regions in the genome assembly. We carried out de novo assembly of the RNA sequences of skin and muscle tissues using Trinity software . The assembled fragments were then aligned to the genome assembly with BLAT  (E-value = 10e-6, identity = 90 % and coverage >90 %). Our results indicate that the catfish genome assembly covered more than 90 % of gene-coding regions. Subsequently, Core Eukaryotic Genes Mapping Approach (CEGMA) software (version 2.3)  was employed with 248 conserved core eukaryotic genes (CEGs) to assess the gene space completeness within the generated genome assembly. These results demonstrate that the genome assembly covered more than 95 % of the CEG sequences, suggesting a high level of completeness.
Total RNA was extracted from muscle and skin tissues of a channel catfish (the same individual used for the above-mentioned genome sequencing) using TRIzol reagent (Invitrogen, USA). After purification using RNeasy Animal Mini Kit (Qiagen, USA), equal amounts of total RNA from each tissue were subjected to transcriptome sequencing (RNA-seq) on the HiSeq 2000 platform.
Firstly, RepeatModeller (version 1.04) and LTR_FINDER  were used to build a de novo repeat library with default parameters. Subsequently, RepeatMasker  (version 3.2.9) was utilized to map our sequences against the Repbase  transposable element (TE) library (version 14.04) and the de novo repeat library, so as to search for known and novel TEs. Next, we annotated tandem repeats using Tandem Repeat Finder  (version 4.04) with core parameters set as “Match = 2, Mismatch = 7, Delta = 7, PM = 80, PI = 10, Minscore = 50, and MaxPerid = 2000”. Furthermore, TE-relevant proteins were identified in our assembly using RepeatProteinMask software  (version 3.2.2). These identified repeat sequences accounted for 32.56 % of the channel catfish genome, of which the single largest class of TEs (representing 9.35 % of the whole genome) was the Tc1-mariner family.
Annotations of gene structure and function
The channel catfish genome assembly was annotated using three independent approaches: homology, de novo and RNA-seq annotations. For homology annotation, the protein sequences from zebrafish, Japanese fugu, spotted green pufferfish, Japanese medaka (Ensembl release 75), blue spotted mudskipper  and golden arowana  were mapped on the channel catfish genome using TblastN with e-value ≤ 1E-5. Genewise 2.2.0 software  was then employed to predict the potential gene structures of all alignments. Short genes (with fewer than 150 bp) and prematurely terminated or frame-shifted genes were discarded. Next, de novo annotation was used to annotate the gene structure from the genome assembly. We randomly selected 1000 complete genes from the homology annotation set to train the parameters for AUGUSTUS 2.5 . Simultaneously, all repetitive regions were replaced in the channel catfish genome with ‘N’ to decline the ratio of pseudogene annotations. Subsequently, we utilized AUGUSTUS 2.5 and GENSCAN 1.0  for de novo prediction of repeat-masked genome sequences. The filtered processes performed on the de novo annotation were the same as those used for homology prediction. Simultaneously, the RNA-seq annotation pipeline was also used to detect gene regions. We employed Tophat 1.2 software  to map the RNA reads extracted from the skin and muscle transcriptomes onto the channel catfish genome sequences. We then sorted and integrated Tophat alignments, and used Cufflink software  to analyze potential gene structures. Results from all three of the above-mentioned annotation pipelines were merged to produce a comprehensive and non-redundant gene set using GLEAN . This gene set contained 21,556 genes with an average of 8.7 exons per gene (Table 1). Because different annotation pipelines were applied, the total gene number predicted here is lower than the 26,661 reported in the American-native channel catfish genome . The Cuffdiff package  of Cufflink software (version 2.0.2.Linux_×86_64) with core parameters (−FDR 0.05 –geometric-norm TRUE –compatible-hits-norm TRUE) was utilized to calculate expression levels according to the GLEAN gene set and Tophat alignments. About 93.4 % of genes were predicted from at least two types of evidence, and approximate 78 % of the genes showed expression activity (fragments per kilobase of exon model per million mapped reads >0) in the skin and muscle tissues.
Simultaneously, all protein sequences from GLEAN results were mapped to SwissProt and TrEMBL  (UniProt release 2011.06) databases using BlastP  with an E-value ≤ 1e-5 to find the best hit for each protein. We also used InterProScan 4.7 software  to align the protein sequences against public databases, including Pfam , PRINTS , ProDom  and SMART , to examine the known motifs and domains in our sequences. Over 94.5 % of these predicted genes possessed at least one related functional assignment from other public databases (SwissProt , Interpro , TrEMBL and KEGG ). In addition, the gene structures (including exon length, intron regions and mRNAs) and exon number distributions (Table 1) were predicted to be similar to other representative teleost species such as zebrafish and medaka.
We generated a channel catfish genome assembly with high quality and comparable structures to other published fish genomes, especially the Coco catfish genome . This new assembly is a valuable resource and reference for further construction of high-density genetic linkage maps and identification of quantitative trait loci for molecular breeding of catfishes.
Availability of supporting data
Supporting data are available in the GigaDB database . Raw whole genome sequencing and transcriptome data are deposited in the SRA under bioproject number PRJNA319455.
BAC, bacterial artificial clone; BES, BAC end sequences/sequencing; CEG, core eukaryotic genes; CEGMA, core eukaryotic genes mapping approach; Gb, gigabases; kb, kilobases; Mb, megabases; TE, transposable element; WGS, whole genome shotgun
You X, Bian C, Zan Q, Xu X, Liu X, Chen J, et al. Mudskipper genomes provide insights into the terrestrial adaptation of amphibious fishes. Nat Commun. 2014;5:5594.
Li R, Fan W, Tian G, Zhu H, He L, Cai J, et al. The sequence and de novo assembly of the giant panda genome. Nature. 2010;463(7279):311–7.
Liu Z, Liu S, Yao J, Bao L, Zhang J, Li Y, et al. The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts. Nat Commun. 2016;7:11757.
Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012;1:12.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29(7):644–52.
Kent WJ. BLAT–the BLAST-like alignment tool. Genome Res. 2002;12(4):656–64.
Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23(9):1061–7.
Xu Z, Wang H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007;35(Web Server issue):W265–8.
Tarailo-Graovac M, Chen N. Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009;Chapter 4:Unit 4. 10.
Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005;110(1–4):462–7.
Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.
Bian C, Hu Y, Ravi V, Kuznetsova IS, Shen X, Mu X, et al. The Asian arowana (Scleropages formosus) genome provides new insights into the evolution of an early lineage of teleosts. Sci Rep. 2016;6:24501.
Birney E, Clamp M, Durbin R. GeneWise and Genomewise. Genome Res. 2004;14(5):988–95.
Stanke M, Keller O, Gunduz I, Hayes A, Waack S, Morgenstern B. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006;34(Web Server issue):W435–9.
Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997;268(1):78–94.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31(1):46–53.
Elsik CG, Mackey AJ, Reese JT, Milshina NV, Roos DS, Weinstock GM. Creating a honey bee consensus gene set. Genome Biol. 2007;8(1):R13.
Bairoch A, Apweiler R. The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 2000;28(1):45–8.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Zdobnov EM, Apweiler R. InterProScan–an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17(9):847–8.
Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, Sconnhammer EL. The Pfam protein families database. Nucleic Acids Res. 2000;28(1):263–6.
Attwood TK, Cronig MD, Flower DR, Lewis AP, Mabey JE, Scordis P, et al. PRINTS-S: the database formerly known as PRINTS. Nucleic Acids Res. 2000;28(1):225–7.
Corpet F, Gouzy J, Kahn D. Recent improvements of the ProDom database of protein domain families. Nucleic Acids Res. 1999;27(1):263–7.
Schult J, Copley RR, Doerks T, Ponting CP, Bork P. SMART: a web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 2000;28(1):231–4.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27(1):29–34.
Chen X, Zhong L, Bian C, Xu P, Qiu Y, You X, Zhang S, Yu H, Li J, Wang M, Qin Q, Zhu X, Peng C, Wong A, Zhu Z, Wang M, Ruobo G, Xu J, Shi Q, Bian W. Supporting data for “High-quality genome assembly of channel catfish, Ictalurus punctatus”. 2016. GigaScience Database, http://dx.doi.org/10.5524/100212.
This study was supported by the National Key Technology R&D Program of China (No. 2012BAD26B03), Fund for Independent Innovation of Agricultural Science and Technology of Jiangsu Province (No. CX(15)1013), Human Resources and Social Security of Jiangsu Province (No. 2014-NY-008), Three-Side Innovation Projects for Aquaculture in Jiangsu Province (No. Y2014-25 &Y2015-12), Shenzhen Special Program for Future Industrial Development (No. JSGG20141020113728803), and Zhenjiang Leading Talent Program for Innovation and Entrepreneurship.
XC, QS, CB, JX and WB conceived the project. LZ, YH, SZ, MHW, QQ, XY, CP, AW, ZZ, MW and RG collected the samples and extracted the genomic DNA. CB, YQ, JL and YH performed the genome assembly and data analysis. CB, QS, XC, LZ, XP, XZ and WB wrote the paper and all authors read and approved the final manuscript.
The authors declare that they have no competing interests.
About this article
Cite this article
Chen, X., Zhong, L., Bian, C. et al. High-quality genome assembly of channel catfish, Ictalurus punctatus. GigaSci 5, 39 (2016). https://doi.org/10.1186/s13742-016-0142-5
- Channel catfish
- Whole genome sequencing
- Gene prediction
- Repetitive sequence