RNA-Seq reveals expression signatures of genes involved in oxygen transport, protein synthesis, folding, and degradation in response to heat stress in catfish.

Temperature is one of the most prominent abiotic factors affecting ectotherms. Most fish species, as ectotherms, have extraordinary ability to deal with a wide range of temperature changes. While the molecular mechanism underlying temperature adaptation has long been of interest, it is still largely unexplored with fish. Understanding of the fundamental mechanisms conferring tolerance to temperature fluctuations is a topic of increasing interest as temperature may continue to rise as a result of global climate change. Catfish have a wide natural habitat and possess great plasticity in dealing with environmental variations in temperature. However, no studies have been conducted at the transcriptomic level to determine heat stress-induced gene expression. In the present study, we conducted an RNA-Seq analysis to identify heat stress-induced genes in catfish at the transcriptome level. Expression analysis identified a total of 2,260 differentially expressed genes with a cutoff of twofold change. qRT-PCR validation suggested the high reliability of the RNA-Seq results. Gene ontology, enrichment, and pathway analyses were conducted to gain insight into physiological and gene pathways. Specifically, genes involved in oxygen transport, protein folding and degradation, and metabolic process were highly induced, while general protein synthesis was dramatically repressed in response to the lethal temperature stress. This is the first RNA-Seq-based expression study in catfish in response to heat stress. The candidate genes identified should be valuable for further targeted studies on heat tolerance, thereby assisting the development of heat-tolerant catfish lines for aquaculture.

Alternative complement pathway of channel catfish (Ictalurus punctatus): molecular characterization, mapping and expression analysis of factors Bf/C2 and Df.

The complement system is important in both innate and adaptive host defense against microbial infection in vertebrates. It contains three pathways: the classical, alternative, and lectin pathways. Complement component factors B and D are two crucial proteases in the alternative pathway. In this study, the genes of complement factors Bf/C2 and Df from channel catfish, Ictalurus punctatus were identified and characterized. Two complement factor B-related genes, Bf/C2A and Bf/C2B, and factor D gene Df were identified. Phylogenetic analysis suggested that Bf/C2A and Bf/C2B is likely orthologous to factor B and factor C2, respectively. Southern blot results suggested that these three genes are all single-copy genes in the catfish genome. The catfish Bf/C2A, Bf/C2B and Df genes were genetically mapped on linkage group 3, 20 and 29, respectively. Bf/C2A and Bf/C2B are highly expressed in liver and kidney, while Df is highly expressed in gill and spleen. After infection with Edwardsiella ictaluri, the expression of Bf/C2A, Bf/C2B and Df genes were found to be remarkably induced in the gill, liver, spleen and kidney at some sampling times, indicating that these three complement factors play a pivotal role in immune responses after the bacterial infection in catfish.

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

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts.

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

Discovery and validation of gene-linked diagnostic SNP markers for assessing hybridization between Largemouth bass (Micropterus salmoides) and Florida bass (M. floridanus).

Efforts to improve recreational fisheries have included widespread stocking of Micropterus floridanus outside its native range of peninsular Florida. Hybridization of Florida bass (M. floridanus) with largemouth bass (Micropterus salmoides) has now dramatically expanded beyond a naturally occurring intergrade zone in the southeast U.S. In recent years, there has been growing interest in protecting the genetic integrity of native basses and assessing the impact and nature of M. salmoides/M. floridanus introgression from the standpoint of hatchery and sport-fishery managers, fish biologists, ecologists and evolutionary biologists. Here, we conducted RNA-seq-based sequencing of the transcriptomes of M. salmoides, M. floridanus and their F1 hybrid and identified a set of 3674 SNP markers with fixed-allelic differences from 2112 unique genes. We then developed a subset of 25 of these markers into a single diagnostic multiplex assay and validated its capacity for assessing integrity and hybridization in hatchery and wild populations of largemouth and Florida bass. The availability of this resource, high-quality transcriptomes and a large set of gene-linked SNPs, should greatly facilitate functional and population genomics studies in these key species and allow the identification of traits and processes under selection during introgressive hybridization.

Identification and analysis of genome-wide SNPs provide insight into signatures of selection and domestication in channel catfish (Ictalurus punctatus).

Domestication and selection for important performance traits can impact the genome, which is most often reflected by reduced heterozygosity in and surrounding genes related to traits affected by selection. In this study, analysis of the genomic impact caused by domestication and artificial selection was conducted by investigating the signatures of selection using single nucleotide polymorphisms (SNPs) in channel catfish (Ictalurus punctatus). A total of 8.4 million candidate SNPs were identified by using next generation sequencing. On average, the channel catfish genome harbors one SNP per 116 bp. Approximately 6.6 million, 5.3 million, 4.9 million, 7.1 million and 6.7 million SNPs were detected in the Marion, Thompson, USDA103, Hatchery strain, and wild population, respectively. The allele frequencies of 407,861 SNPs differed significantly between the domestic and wild populations. With these SNPs, 23 genomic regions with putative selective sweeps were identified that included 11 genes. Although the function for the majority of the genes remain unknown in catfish, several genes with known function related to aquaculture performance traits were included in the regions with selective sweeps. These included hypoxia-inducible factor 1ß. HIFιß.. and the transporter gene ATP-binding cassette sub-family B member 5 (ABCB5). HIF1ß. is important for response to hypoxia and tolerance to low oxygen levels is a critical aquaculture trait. The large numbers of SNPs identified from this study are valuable for the development of high-density SNP arrays for genetic and genomic studies of performance traits in catfish.

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

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Channel catfish hemoglobin genes: identification, phylogenetic and syntenic analysis, and specific induction in response to heat stress.

Hemoglobins transport oxygen from gill to inner organs in fish, and this process is affected by temperature, one of the major environmental factors for fish. The hemoglobin gene clusters have been well studied in humans and several model fish species, but remain largely unknown in catfish. Here, eight α- and six ß-hemoglobin genes were identified and characterized in channel catfish. Genomic synteny analysis showed that these hemoglobin genes were separated into two unlinked clusters, the MN cluster containing six α- and six ß-hemoglobin genes, and the LA cluster consisting of two α-hemoglobin genes. Channel catfish hemoglobin genes were ubiquitously expressed in all the 10 tested tissues from healthy fish, but exhibited higher expression level in spleen, head kidney, and trunk kidney. In response to heat stress, hemoglobin genes, especially MN Hbα4, MN Hbα5, MN Hbα6, MN Hbß4, MN Hbß5, MN Hbß6, LA Hbα1, and LA Hbα2, presumably the embryonic hemoglobin genes, were drastically up-regulated in the gill and head kidney of heat-tolerant fishes, but not in these tissues of the heat-intolerant fish, suggesting the importance of the embryonic hemoglobin genes in coping with the low oxygen conditions under heat stress.

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.

Male-biased genes in catfish as revealed by RNA-Seq analysis of the testis transcriptome.

BACKGROUND: Catfish has a male-heterogametic (XY) sex determination system, but genes involved in gonadogenesis, spermatogenesis, testicular determination, and sex determination are poorly understood. As a first step of understanding the transcriptome of the testis, here, we conducted RNA-Seq analysis using high throughput Illumina sequencing. METHODOLOGY/PRINCIPAL FINDINGS: A total of 269.6 million high quality reads were assembled into 193,462 contigs with a N50 length of 806 bp. Of these contigs, 67,923 contigs had hits to a set of 25,307 unigenes, including 167 unique genes that had not been previously identified in catfish. A meta-analysis of expressed genes in the testis and in the gynogen (double haploid female) allowed the identification of 5,450 genes that are preferentially expressed in the testis, providing a pool of putative male-biased genes. Gene ontology and annotation analysis suggested that many of these male-biased genes were involved in gonadogenesis, spermatogenesis, testicular determination, gametogenesis, gonad differentiation, and possibly sex determination. CONCLUSION/SIGNIFICANCE: We provide the first transcriptome-level analysis of the catfish testis. Our analysis would lay the basis for sequential follow-up studies of genes involved in sex determination and differentiation in catfish.

Second-generation genetic linkage map of catfish and its integration with the BAC-based physical map.

Construction of high-density genetic linkage maps is crucially important for quantitative trait loci (QTL) studies, and they are more useful when integrated with physical maps. Such integrated maps are valuable genome resources for fine mapping of QTL, comparative genomics, and accurate and efficient whole-genome assembly. Previously, we established both linkage maps and a physical map for channel catfish, Ictalurus punctatus, the dominant aquaculture species in the United States. Here we added 2030 BAC end sequence (BES)-derived microsatellites from 1481 physical map contigs, as well as markers from singleton BES, ESTs, anonymous microsatellites, and SNPs, to construct a second-generation linkage map. Average marker density across the 29 linkage groups reached 1.4 cM/marker. The increased marker density highlighted variations in recombination rates within and among catfish chromosomes. This work effectively anchored 44.8% of the catfish BAC physical map contigs, covering ~52.8% of the genome. The genome size was estimated to be 2546 cM on the linkage map, and the calculated physical distance per centimorgan was 393 Kb. This integrated map should enable comparative studies with teleost model species as well as provide a framework for ordering and assembling whole-genome scaffolds.

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

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.