Generation of genome-scale gene-associated SNPs in catfish for the construction of a high-density SNP array.

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genome-wide association studies. In order to provide the best genome coverage for the analysis of performance and production traits, a large number of relatively evenly distributed SNPs are needed. Gene-associated SNPs may fulfill these requirements of large numbers and genome wide distribution. In addition, gene-associated SNPs could themselves be causative SNPs for traits. The objective of this project was to identify large numbers of gene-associated SNPs using high-throughput next generation sequencing. Transcriptome sequencing was conducted for channel catfish and blue catfish using Illumina next generation sequencing technology. Approximately 220 million reads (15.6 Gb) for channel catfish and 280 million reads (19.6 Gb) for blue catfish were obtained by sequencing gene transcripts derived from various tissues of multiple individuals from a diverse genetic background. A total of over 35 billion base pairs of expressed short read sequences were generated. Over two million putative SNPs were identified from channel catfish and almost 2.5 million putative SNPs were identified from blue catfish. Of these putative SNPs, a set of filtered SNPs were identified including 342,104 intra-specific SNPs for channel catfish, 366,269 intra-specific SNPs for blue catfish, and 420,727 inter-specific SNPs between channel catfish and blue catfish. These filtered SNPs are distributed within 16,562 unique genes in channel catfish and 17,423 unique genes in blue catfish. For aquaculture species, transcriptome analysis of pooled RNA samples from multiple individuals using Illumina sequencing technology is both technically efficient and cost-effective for generating expressed sequences. Such an approach is most effective when coupled to existing EST resources generated using traditional sequencing approaches because the reference ESTs facilitate effective assembly of the expressed short reads. When multiple individuals with different genetic backgrounds are used, RNA-Seq is very effective for the identification of SNPs. The SNPs identified in this report will provide a much needed resource for genetic studies in catfish and will contribute to the development of a high-density SNP array. Validation and testing of these SNPs using SNP arrays will form the material basis for genome association studies and whole genome-based selection in catfish.

Liu S ，Zhou Z ，Lu J ，Sun F ，Wang S ，Liu H ，Jiang Y ，Kucuktas H ，Kaltenboeck L ，Peatman E ，Liu Z ... -  《BMC GENOMICS》

Assembly of 500,000 inter-specific catfish expressed sequence tags and large scale gene-associated marker development for whole genome association studies.

Through the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energy's Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification. A total of 438,321 quality ESTs were generated from 8 channel catfish (Ictalurus punctatus) and 4 blue catfish (Ictalurus furcatus) libraries, bringing the number of catfish ESTs to nearly 500,000. Assembly of all catfish ESTs resulted in 45,306 contigs and 66,272 singletons. Over 35% of the unique sequences had significant similarities to known genes, allowing the identification of 14,776 unique genes in catfish. Over 300,000 putative SNPs have been identified, of which approximately 48,000 are high-quality SNPs identified from contigs with at least four sequences and the minor allele presence of at least two sequences in the contig. The EST resource should be valuable for identification of microsatellites, genome annotation, large-scale expression analysis, and comparative genome analysis. This project generated a large EST resource for catfish that captured the majority of the catfish transcriptome. The parallel analysis of ESTs from two closely related Ictalurid catfishes should also provide powerful means for the evaluation of ancient and recent gene duplications, and for the development of high-density microarrays in catfish. The inter- and intra-specific SNPs identified from all catfish EST dataset assembly will greatly benefit the catfish introgression breeding program and whole genome association studies.

Wang S ，Peatman E ，Abernathy J ，Waldbieser G ，Lindquist E ，Richardson P ，Lucas S ，Wang M ，Li P ，Thimmapuram J ，Liu L ，Vullaganti D ，Kucuktas H ，Murdock C ，Small BC ，Wilson M ，Liu H ，Jiang Y ，Lee Y ，Chen F ，Lu J ，Wang W ，Xu P ，Somridhivej B ，Baoprasertkul P ，Quilang J ，Sha Z ，Bao B ，Wang Y ，Wang Q ，Takano T ，Nandi S ，Liu S ，Wong L ，Kaltenboeck L ，Quiniou S ，Bengten E ，Miller N ，Trant J ，Rokhsar D ，Liu Z ，Catfish Genome Consortium ... -  《GENOME BIOLOGY》

Efficient assembly and annotation of the transcriptome of catfish by RNA-Seq analysis of a doubled haploid homozygote.

Upon the completion of whole genome sequencing, thorough genome annotation that associates genome sequences with biological meanings is essential. Genome annotation depends on the availability of transcript information as well as orthology information. In teleost fish, genome annotation is seriously hindered by genome duplication. Because of gene duplications, one cannot establish orthologies simply by homology comparisons. Rather intense phylogenetic analysis or structural analysis of orthologies is required for the identification of genes. To conduct phylogenetic analysis and orthology analysis, full-length transcripts are essential. Generation of large numbers of full-length transcripts using traditional transcript sequencing is very difficult and extremely costly. In this work, we took advantage of a doubled haploid catfish, which has two sets of identical chromosomes and in theory there should be no allelic variations. As such, transcript sequences generated from next-generation sequencing can be favorably assembled into full-length transcripts. Deep sequencing of the doubled haploid channel catfish transcriptome was performed using Illumina HiSeq 2000 platform, yielding over 300 million high-quality trimmed reads totaling 27 Gbp. Assembly of these reads generated 370,798 non-redundant transcript-derived contigs. Functional annotation of the assembly allowed identification of 25,144 unique protein-encoding genes. A total of 2,659 unique genes were identified as putative duplicated genes in the catfish genome because the assembly of the corresponding transcripts harbored PSVs or MSVs (in the form of pseudo-SNPs in the assembly). Of the 25,144 contigs with unique protein hits, around 20,000 contigs matched 50% length of reference proteins, and over 14,000 transcripts were identified as full-length with complete open reading frames. The characterization of consensus sequences surrounding start codon and the stop codon confirmed the correct assembly of the full-length transcripts. The large set of transcripts assembled in this study is the most comprehensive set of genome resources ever developed from catfish, which will provide the much needed resources for functional genome research in catfish, serving as a reference transcriptome for genome annotation, analysis of gene duplication, gene family structures, and digital gene expression analysis. The putative set of duplicated genes provide a starting point for genome scale analysis of gene duplication in the catfish genome, and should be a valuable resource for comparative genome analysis, genome evolution, and genome function studies.

Liu S ，Zhang Y ，Zhou Z ，Waldbieser G ，Sun F ，Lu J ，Zhang J ，Jiang Y ，Zhang H ，Wang X ，Rajendran KV ，Khoo L ，Kucuktas H ，Peatman E ，Liu Z ... -  《BMC GENOMICS》

Transcriptome analysis and SNP identification in SCC of horn in (Bos indicus) Indian cattle.

Single Nucleotide Polymorphisms (SNPs) have become the marker of choice for genome wide association studies. In order to provide the best genome coverage for the analysis of disease, production and performance traits, a large number of relatively evenly distributed SNPs are needed. The main objective of present work was to identify large numbers of gene-associated SNPs using high-throughput sequencing in squamous cell carcinoma of horn. RNA-seq analysis was conducted on 2 tissues viz. Horn Cancer (HC) and Horn Normal (HN) in Kankrej breed of cattle. A total of 909,362 reads with average read length of 405 bp for HC and 583,491 reads with average read length of 411 bp for HN were obtained. We found 9532 and 7065 SNPs as well as 1771 and 1172 Indels in HC and HN, respectively, from which, 7889 SNPs and 1736 Indels were uniquely present in HC, 5886 SNPs and 1146 Indels were uniquely present in HN and reported first time in Bos indicus, whereas the rest are already reported in Bos taurus dbSNP database. The gene-associated SNPs and Indels were high in upregulated genes of HC as compared to HN. Analysis of differentially expressed genes was identified, these genes are involved in regulation of cell proliferation, apoptosis, gene transcription, cell survival and metabolism through various metabolic pathways. The result of transcriptome expression profiling was validated using Real Time quantitative PCR in nine randomly selected genes. We identified numbers aberrant signaling pathways responsible for carcinogenesis in HC which are also commonly altered in squamous cell carcinoma (SCC) of lung in human being. We conclude that a large number of altered genes and dysfunction of multiple pathways are involved in the development of Horn Cancer. The present findings contribute to theoretical information for further screening of genes and identification of markers for early diagnosis of HC as well as SNPs identified in this report provide a much needed resource for genetic studies in B. indicus and shall contribute to the development of a high density SNP array. Validation and testing of these SNPs using SNP arrays will form the material basis for gene associated SNPs in HC.

Koringa PG ，Jakhesara SJ ，Bhatt VD ，Patel AB ，Dash D ，Joshi CG ... -  《-》

Development of the catfish 250K SNP array for genome-wide association studies.

Quantitative traits, such as disease resistance, are most often controlled by a set of genes involving a complex array of regulation. The dissection of genetic basis of quantitative traits requires large numbers of genetic markers with good genome coverage. The application of next-generation sequencing technologies has allowed discovery of over eight million SNPs in catfish, but the challenge remains as to how to efficiently and economically use such SNP resources for genetic analysis. In this work, we developed a catfish 250K SNP array using Affymetrix Axiom genotyping technology. The SNPs were obtained from multiple sources including gene-associated SNPs, anonymous genomic SNPs, and inter-specific SNPs. A set of 640K high-quality SNPs obtained following specific requirements of array design were submitted. A panel of 250,113 SNPs was finalized for inclusion on the array. The performance evaluated by genotyping individuals from wild populations and backcross families suggested the good utility of the catfish 250K SNP array. This is the first high-density SNP array for catfish. The array should be a valuable resource for genome-wide association studies (GWAS), fine QTL mapping, high-density linkage map construction, haplotype analysis, and whole genome-based selection.

Liu S ，Sun L ，Li Y ，Sun F ，Jiang Y ，Zhang Y ，Zhang J ，Feng J ，Kaltenboeck L ，Kucuktas H ，Liu Z ... -  《-》